Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2025

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2024

Magnetism and finite-temperature effects in UZr 2 : A density functional theory analysis

The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries

The formation of a solid electrolyte interphase (SEI) at the electrode/electrolyte interface substantially affects the stability and lifetime of lithium-ion batteries (LIBs). One of the methods to improve the lifetime of LIBs is by the inclusion of additive molecules to stabilize the SEI. To understand the effect of additive molecules on the initial stage of SEI formation, in this work we compare the decomposition and oligomerization reactions of a fluoroethylene carbonate (FEC) additive on a range of oxygen-functionalized graphitic anodes to those of an ethylene carbonate (EC) organic electrolyte. A series of density functional theory (DFT) calculations augmented by ab initio molecular dynamics (AIMD) simulations reveal that EC decomposition on an oxygen-functionalized graphitic ($11\bar{2}0$) edge facet through a nucleophilic attack on an ethylene carbon site (C E ) of an EC molecule (S2 mechanism) is spontaneous during the initial charging process of LIBs. However, decomposition of EC through a nucleophilic attack on a carbonyl carbon (C C ) site (S1 mechanism) results in alkoxide species regeneration that is responsible for continual oligomerization along the graphitic surface. In contrast, FEC prefers to decompose through an S1 pathway, which does not promote alkoxide regeneration. Including FEC as an additive is thus able to suppress alkoxide regeneration and results in a smaller and thinner SEI layer that is more flexible toward lithium intercalation during the charging/discharging process. In addition, we find that the presence of different oxygen functional groups at the surface of graphite dictates the oligomerization products and the LiF formation mechanism in the SEI.

25 ENERGY STORAGE↗
DOE OSTI2023

First-Principles-Based Machine Learning Models for Phase Behavior and Transport Properties of CO 2

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN and SCAN-rvv10 approximations of density functional theory. Here we employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

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

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2025

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

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

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Cation Effects on CO 2 Delivery to Cu Electrode in Reactive Capture of CO 2

The direct electrochemical conversion of captured CO 2 , known as reactive capture of CO 2 (RCC), remains a formidable challenge in heterogeneous catalysis. Given that amines are one of the most widely used capture agents for CO 2 , it would be desirable to electrochemically reduce the resultant adducts, such as carbamate, directly in RCC. However, current understanding suggests that the primary species undergoing reduction in RCC with amines is the CO 2 dissociated from the sorbent. Herein, we employ ab initio molecular dynamics (AIMD) with DFT to analyze how the nature of alkali metal cations in the electrolyte affects carbamate at the Cu surface, thereby assessing the possibility of promoting RCC by cation effects. The simulations show that the carbamate’s orientation with respect to the electrode is governed by the optimal distance between the carbamate and the cation, specifically how this distance aligns with the cation’s hydration spheres. Moreover, the slow-growth AIMD results indicate that the CO 2 dissociation barrier correlates with the orientation of carbamate at the interface. When the carbamate resides beyond the cation’s first hydration sphere, it adopts a flat orientation with respect to the surface that promotes the release of CO 2 from the capture agent. In contrast, when the carbamate disrupts the first hydration sphere and exhibits a strong cation−π interaction, it adopts an upright orientation that is less conducive to CO 2 release. These findings reveal a nontrivial cation effect in RCC, suggesting that it should be possible to optimize RCC via the choice of the electrolyte.

ab initio molecular dynamics↗
DOE OSTI2022

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe–N–C Materials

The systematic improvement of Fe–N–C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe–N–C materials is the formation of the Fe III –OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of Fe III –OH to Fe II and desorbed H 2 O requires an axial H 2 O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. Importantly, these findings highlight the importance of solvent–substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

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

Ordering in liquid and its heredity impact on phase transformation of Mg-Al-Ca alloys

It is a long-sought goal to achieve desired mechanical properties through tailoring phase formation in alloys, especially for complicated multi-phase alloys. In fact, unveiling nucleation of competitive crystalline phases during solidification hinges on the nature of liquid. Here we employ ab initio molecular dynamics simulations (AIMD) to reveal liquid configuration of the Mg-Al-Ca alloys and explore its effect on the transformation of Ca-containing Laves phase from Al 2 Ca to Mg 2 Ca with increasing Ca/Al ratio (r Ca/Al ). There is structural similarity between liquid and crystalline phase in terms of the local arrangement environment, and the connection schemes of polyhedras. The forming signature of Mg 2 Ca, as hinted by the topological and chemical short-range order originating from liquid, ascends monotonically with increasing r Ca/Al . However, Al 2 Ca crystal-like order increase at first and then decrease at the crossover of r Ca/Al = 0.74, corresponding to experimental composition of phase transition from Al 2 Ca to Mg 2 Ca. The origin of phase transformation across different compositions lies in the dense packing of atomic configurations and preferential bonding of chemical species in both liquid and solid. The present finding provides a feasible scenario for manipulating phase formation to achieve high performance alloys by tailoring the crystal-like order in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of functional groups on the electrocatalytic hydrogenation of aromatic carbonyls to alcohols

We report that electrocatalytic hydrogenation (ECH) of biomass-derived feedstocks has a critical dependence on the molecular structure of the organic and its adsorption on the electrode surface. In this study, we investigated the role of functional groups in the adsorption of the organic molecule on the charged Pd (111) surface and its subsequent effect on organic reduction in electrochemical hydrogenation of organic molecules. With three aromatic carbonyls of benzaldehyde (BZD), acetophenone (ACE), and vanillin (VAN), we rationalize molecular-scale adsorption and interfacial charge transfer processes by employing density-functional-theory based ab initio molecular dynamics simulations. We observe that the functional group and electrode charge strongly affect the proximity of organic molecule to the Pd (111) surface, where distances of aromatic ring and carbonyl group of the organic on the electrode change distinctively with functional groups and charge state of electrode, which strongly impact reduction of organics on the surface. Calculations of differential electron density show the strongest reduction with benzaldehyde via interfacial electron transfer from the charged Pd surface. We also observe that the interaction between the functional groups and solvent (VAN > BZD > ACE) significantly influence the organic interaction with the charged electrode (BZD > VAN > ACE), resulting in the net interaction energy between the organic and the electrode in the order of BZD > ACE > VAN. Experimental measurement of ECH rate also show the same trend of the net interaction energy. These results demonstrate the significance of solvent effect on the reducibility of organic molecules on electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Superionic-like diffusion in yttrium dihydride

For the next-generation high temperature microreactors, yttrium dihydride (YH 2 ) is an attractive solid state neutron moderator. Despite a number of recent investigations, the mechanism of hydrogen transport remains poorly understood. Experimental evaluations of diffusivity are inconclusive with large variations in diffusivities and activation energies. In this work, we perform ab initio molecular dynamics (AIMD) simulations on YH 2 for temperatures spanning 300 K to 1200 K. Our main finding is that YH 2 shows a superionic-like behavior with hydrogen atoms hopping from one native site to another above a characteristic temperature of 800 K. This correlated motion results in quasi-one-dimensional string-like displacements that enable the hydrogen atoms to diffuse rapidly. We confirm that the octahedral sites are mostly unoccupied, although channeling through them is the most favored pathway between lattice hops above 800 K. At the highest temperature of 1200 K, the string relaxation time is merely of the order of a few picoseconds, which indicates a liquid-like diffusive behavior. Based on the formation of spontaneous thermal vacancies, an order-disorder crossover temperature T α ~ 800 K is established for YH 2 with an activation energy of 0.83 eV for hydrogen diffusion in the superionic-like state.

Superionic-like Diffusion↗
DOE OSTI2023

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile

Lanthanide–ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu 3+ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu 3+ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu 3+ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu 3+ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu 3+ coordinating molecules as the crystal structure. Finally, this study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

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

Cooperative Effects Associated with High Electrolyte Concentrations in Driving the Conversion of CO2 to C2H4 on Copper

Compared to a conventional electrolyte concentration of 1 M HCOOK, the use of a highly concentrated 7.1 M HCOOK electrolyte increases the Faradaic efficiency (FE) ratio of C2H4/CO from 2.2 +- 0.3 to 18.3 +- 4.8 at -1.08 V vs. reversible hydrogen electrode (RHE) on a Cu gas-diffusion electrode. Based on electrochemical analysis and ab initio molecular dynamics (AIMD) simulation, the identity and concentration of the cation and anion play more important roles in controlling the CO2R reaction pathway than the bulk CO2 solubility and the bulk pH of electrolytes. In situ attenuated reflectance surface enhanced infrared absorption spectroscopy (ATR-SEIRAS) suggests that, unlike 1 M HCOOK, the *CO-bridge-binding mode on Cu is dominant in 7.1 M HCOOK electrolyte, which potentially results in less CO release and higher yield of C2H4. This study demonstrates that although we can tailor the electrolyte composition to shift product selectivity, the factors that control the product selectivity are numerous and cannot be distilled down into one correlated property-reactivity relationship.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗