“Ensemble-Based Molecular Simulation of Chemical Reactions under Vibrational Nonequilibrium”. Bal KM, Bogaerts A, Neyts EC, Journal Of Physical Chemistry Letters 11, 401 (2020). http://doi.org/10.1021/acs.jpclett.9b03356
Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 5.7
DOI: 10.1021/acs.jpclett.9b03356
|
“Accelerated molecular dynamics simulation of large systems with parallel collective variable-driven hyperdynamics”. Fukuhara S, Bal KM, Neyts EC, Shibuta Y, Computational Materials Science 177, 109581 (2020). http://doi.org/10.1016/j.commatsci.2020.109581
Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.3
DOI: 10.1016/j.commatsci.2020.109581
|
“Activation of CO2on Copper Surfaces: The Synergy between Electric Field, Surface Morphology, and Excess Electrons”. Jafarzadeh A, Bal KM, Bogaerts A, Neyts EC, Journal Of Physical Chemistry C 124, 6747 (2020). http://doi.org/10.1021/acs.jpcc.0c00778
Abstract: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.7
DOI: 10.1021/acs.jpcc.0c00778
|
“Predicted Influence of Plasma Activation on Nonoxidative Coupling of Methane on Transition Metal Catalysts”. Engelmann Y, Mehta P, Neyts EC, Schneider WF, Bogaerts A, Acs Sustainable Chemistry &, Engineering 8, 6043 (2020). http://doi.org/10.1021/acssuschemeng.0c00906
Abstract: The combination of catalysis and nonthermal plasma holds promise for enabling difficult chemical conversions. The possible synergy between both depends strongly on the nature of the reactive plasma species and the catalyst material. In this paper, we show how vibrationally excited species and plasma-generated radicals interact with transition metal catalysts and how changing the catalyst material can improve the conversion rates and product selectivity. We developed a microkinetic model to investigate the impact of vibrational excitations and plasma-generated radicals on the nonoxidative coupling of methane over transition metal surfaces. We predict a significant increase in ethylene formation for vibrationally excited methane. Plasma-generated radicals have a stronger impact on the turnover frequencies with high selectivity toward ethylene on noble catalysts and mixed selectivity on non-noble catalysts. In general, we show how the optimal catalyst material depends on the desired products as well as the plasma conditions.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT); Movement Antwerp (MOVANT)
Impact Factor: 8.4
DOI: 10.1021/acssuschemeng.0c00906
|
“Mechanisms of selective nanocarbon synthesis inside carbon nanotubes”. Khalilov U, Neyts EC, Carbon 171, 72 (2021). http://doi.org/10.1016/j.carbon.2020.08.060
Abstract: The possibility of confinement effects inside a carbon nanotube provides new application opportunities, e.g., growth of novel carbon nanostructures. However, the understanding the precise role of catalystfeedstock in the nanostructure synthesis is still elusive. In our simulation-based study, we investigate the Ni-catalyzed growth mechanism of encapsulated carbon nanostructures, viz. double-wall carbon nanotube and graphene nanoribbon, from carbon and hydrocarbon growth precursors, respectively. Specifically, we find that the tube and ribbon growth is determined by a catalyst-vs-feedstock competition effect. We compare our results, i.e., growth mechanism and structure morphology with all available theoretical and experimental data. Our calculations show that all encapsulated nanostructures contain metal (catalyst) atoms and such structures are less stable than their pure counterparts. Therefore, we study the purification mechanism of these structures. In general, this study opens a possible route to the controllable synthesis of tubular and planar carbon nanostructures for today’s nanotechnology.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 6.337
DOI: 10.1016/j.carbon.2020.08.060
|
“Free energy barriers from biased molecular dynamics simulations”. Bal KM, Fukuhara S, Shibuta Y, Neyts EC, Journal Of Chemical Physics 153, 114118 (2020). http://doi.org/10.1063/5.0020240
Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.4
DOI: 10.1063/5.0020240
|
“Entropic and enthalpic factors determining the thermodynamics and kinetics of carbon segregation from transition metal nanoparticles”. Fukuhara S, Bal KM, Neyts EC, Shibuta Y, Carbon 171, 806 (2021). http://doi.org/10.1016/j.carbon.2020.09.059
Abstract: The free energy surface (FES) for carbon segregation from nickel nanoparticles is obtained from advanced molecular dynamics simulations. A suitable reaction coordinate is developed that can distinguish dissolved carbon atoms from segregated dimers, chains and junctions on the nanoparticle surface. Because of the typically long segregation time scale (up to ms), metadynamics simulations along the developed reaction coordinate are used to construct FES over a wide range of temperatures and carbon concentrations. The FES revealed the relative stability of different stages in the segregation process, and free energy barriers and rates of the individual steps could then be calculated and decomposed into enthalpic and entropic contributions. As the carbon concentration in the nickel nanoparticle increases, segregated carbon becomes more stable in terms of both enthalpy and entropy. The activation free energy of the reaction also decreases with the increase of carbon concentration, which can be mainly attributed to entropic effects. These insights and the methodology developed to obtain them improve our understanding of carbon segregation process across materials science in general, and the nucleation and growth of carbon nanotube in particular.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 6.337
DOI: 10.1016/j.carbon.2020.09.059
|
“Reactive plasma cleaning and restoration of transition metal dichalcogenide monolayers”. Marinov D, de Marneffe J-F, Smets Q, Arutchelvan G, Bal KM, Voronina E, Rakhimova T, Mankelevich Y, El Kazzi S, Nalin Mehta A, Wyndaele P-J, Heyne MH, Zhang J, With PC, Banerjee S, Neyts EC, Asselberghs I, Lin D, De Gendt S, npj 2D Materials and Applications 5, 17 (2021). http://doi.org/10.1038/s41699-020-00197-7
Abstract: The cleaning of two-dimensional (2D) materials is an essential step in the fabrication of future devices, leveraging their unique physical, optical, and chemical properties. Part of these emerging 2D materials are transition metal dichalcogenides (TMDs). So far there is limited understanding of the cleaning of “monolayer” TMD materials. In this study, we report on the use of downstream H<sub>2</sub>plasma to clean the surface of monolayer WS<sub>2</sub>grown by MOCVD. We demonstrate that high-temperature processing is essential, allowing to maximize the removal rate of polymers and to mitigate damage caused to the WS<sub>2</sub>in the form of sulfur vacancies. We show that low temperature in situ carbonyl sulfide (OCS) soak is an efficient way to resulfurize the material, besides high-temperature H<sub>2</sub>S annealing. The cleaning processes and mechanisms elucidated in this work are tested on back-gated field-effect transistors, confirming that transport properties of WS<sub>2</sub>devices can be maintained by the combination of H<sub>2</sub>plasma cleaning and OCS restoration. The low-damage plasma cleaning based on H<sub>2</sub>and OCS is very reproducible, fast (completed in a few minutes) and uses a 300 mm industrial plasma etch system qualified for standard semiconductor pilot production. This process is, therefore, expected to enable the industrial scale-up of 2D-based devices, co-integrated with silicon technology.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
DOI: 10.1038/s41699-020-00197-7
|
“Plasma-Catalytic Ammonia Reforming of Methane over Cu-Based Catalysts for the Production of HCN and H2at Reduced Temperature”. Yi Y, Wang X, Jafarzadeh A, Wang L, Liu P, He B, Yan J, Zhang R, Zhang H, Liu X, Guo H, Neyts EC, Bogaerts A, Acs Catalysis 11, 1765 (2021). http://doi.org/10.1021/acscatal.0c04940
Abstract: Industrial production of HCN from NH3 and CH4 not only uses precious Pt or Pt−Rh catalysts but also requires extremely high temperatures (∼1600 K). From an energetic, operational, and safety perspective, a drastic decrease in temperature is highly desirable. Here, we report ammonia reforming of methane for the production of HCN and H2 at 673 K by the combination of CH4/NH3 plasma and a supported Cu/silicalite-1 catalyst. 30% CH4 conversion has been achieved with 79% HCN selectivity. Catalyst characterization and plasma diagnostics reveal that the excellent reaction performance is attributed to metallic Cu active sites. In addition, we propose a possible reaction pathway, viz. E-R reactions with N, NH, NH2, and CH radicals produced in the plasma, for the production of HCN, based on density functional theory calculations. Importantly, the Cu/silicalite-1 catalyst costs less than 5% of the commercial Pt mesh catalyst.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 10.614
DOI: 10.1021/acscatal.0c04940
|
“Effect of chemical modification on electronic transport properties of carbyne”. Berdiyorov GR, Khalilov U, Hamoudi H, Neyts EC, Journal Of Computational Electronics (2021). http://doi.org/10.1007/s10825-020-01639-7
Abstract: Using density functional theory in combination with the Green’s functional formalism, we study the effect of surface functionalization on the electronic transport properties of 1D carbon allotrope—carbyne. We found that both hydrogenation and fluorination result in structural changes and semiconducting to metallic transition. Consequently, the current in the functionalization systems increases significantly due to strong delocalization of electronic states along the carbon chain. We also study the electronic transport in partially hydrogenated carbyne and interface structures consisting of pristine and functionalized carbyne. In the latter case, current rectification is obtained in the system with rectification ratio up to 50%. These findings can be useful for developing carbyne-based structures with tunable electronic transport properties.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 1.526
DOI: 10.1007/s10825-020-01639-7
|
“Reaction mechanisms of C(3PJ) and C+(2PJ) with benzene in the interstellar medium from quantum mechanical molecular dynamics simulations”. Izadi ME, Bal KM, Maghari A, Neyts EC, Physical Chemistry Chemical Physics 23, 4205 (2021). http://doi.org/10.1039/D0CP04542J
Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>and C<sup>+</sup>(<sup>2</sup>P) + C<sub>6</sub>H<sub>6</sub>systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C<sub>7</sub>H<sub>5</sub>) and benzocyclopropenyl radical (C<sub>6</sub>H<sub>5</sub>–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C<sup>+</sup>(<sup>2</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(<sup>3</sup>P<sub>J</sub>) + C<sub>6</sub>H<sub>6</sub>system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.123
DOI: 10.1039/D0CP04542J
|
Bal KM, Neyts EC (2021) Quantifying the impact of vibrational nonequilibrium in plasma catalysis: insights from a molecular dynamics model of dissociative chemisorption. 394004
Abstract: The rate, selectivity and efficiency of plasma-based conversion processes is strongly affected by nonequilibrium phenomena. High concentrations of vibrationally excited molecules are such a plasma-induced effect. It is frequently assumed that vibrationally excited molecules are important in plasma catalysis because their presence lowers the apparent activation energy of dissociative chemisorption reactions and thus increases the conversion rate. A detailed atomic-level understanding of vibrationally stimulated catalytic reactions in the context of plasma catalysis is however lacking. Here, we couple a recently developed statistical model of a plasma-induced vibrational nonequilibrium to molecular dynamics simulations, enhanced sampling methods, and machine learning techniques. We quantify the impact of a vibrational nonequilibrium on the dissociative chemisorption barrier of H2 and CH4 on nickel catalysts over a wide range of vibrational temperatures. We investigate the effect of surface structure and compare the role of different vibrational modes of methane in the dissociation process. For low vibrational temperatures, very high vibrational efficacies are found, and energy in bend vibrations appears to dominate the dissociation of methane. The relative impact of vibrational nonequilibrium is much higher on terrace sites than on surface steps. We then show how our simulations can help to interpret recent experimental results, and suggest new paths to a better understanding of plasma catalysis.
Keywords: A1 Journal Article;plasma catalysis; vibrational nonequilibrium; dissociative chemisorption; free energy barriers; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 2.588
DOI: 10.1088/1361-6463/ac113a
|
Engelmann Y, van &rsquo,t Veer K, Gorbanev Y, Neyts EC, Schneider WF, Bogaerts A (2021) Plasma Catalysis for Ammonia Synthesis: A Microkinetic Modeling Study on the Contributions of Eley–Rideal Reactions. 13151–13163
Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.
Keywords: A1 Journal Article;Plasma catalysis; Eley−Rideal reactions; Volcano plots; Vibrational excitation; Radical reactions; Dielectric barrier discharge; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 5.951
DOI: 10.1021/acssuschemeng.1c02713
|
“Foundations of plasma catalysis for environmental applications”. Bogaerts A, Neyts EC, Guaitella O, Murphy AB, Plasma Sources Science &, Technology (2022). http://doi.org/10.1088/1361-6595/ac5f8e
Abstract: Plasma catalysis is gaining increasing interest for various applications, but the underlying mechanisms are still far from understood. Hence, more fundamental research is needed to understand these mechanisms. This can be obtained by both modelling and experiments. This foundations paper describes the fundamental insights in plasma catalysis, as well as efforts to gain more insights by modelling and experiments. Furthermore, it discusses the state-of-the-art of the major plasma catalysis applications, as well as successes and challenges of technology transfer of these applications.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.8
DOI: 10.1088/1361-6595/ac5f8e
|
“Identification of a Robust and Durable FeN4CxCatalyst for ORR in PEM Fuel Cells and the Role of the Fifth Ligand”. Nematollahi P, Barbiellini B, Bansil A, Lamoen D, Qingying J, Mukerjee S, Neyts EC, ACS catalysis , 7541 (2022). http://doi.org/10.1021/acscatal.2c01294
Abstract: Although recent studies have advanced the understanding of pyrolyzed
Fe−N−C materials as oxygen reduction reaction (ORR) catalysts, the atomic and
electronic structures of the active sites and their detailed reaction mechanisms still remain unknown. Here, based on first-principles density functional theory (DFT) computations, we discuss the electronic structures of three FeN4 catalytic centers with different local topologies of the surrounding C atoms with a focus on unraveling the mechanism of their ORR activity in acidic electrolytes. Our study brings back a forgotten, synthesized pyridinic Fe−N coordinate to the community’s attention, demonstrating that this catalyst can exhibit excellent activity for promoting direct four-electron ORR through the addition of a fifth ligand such as −NH2, −OH, and −SO4. We also identify sites with good stability properties through the combined use of our DFT calculations and Mössbauer spectroscopy data.
Keywords: A1 Journal article; Electron microscopy for materials research (EMAT); Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 12.9
DOI: 10.1021/acscatal.2c01294
|
“Comment on “Misinterpretation of the Shuttleworth equation””. Faraji F, Neyts EC, Milošević, MV, Peeters FM, Scripta Materialia 250, 116186 (2024). http://doi.org/10.1016/j.scriptamat.2024.116186
Keywords: A1 Journal Article; CMT
Impact Factor: 6
DOI: 10.1016/j.scriptamat.2024.116186
|
“Capillary Condensation of Water in Graphene Nanocapillaries”. Faraji F, Neyts EC, Milošević, MV, Peeters FM, Nano Letters 24, 5625 (2024). http://doi.org/10.1021/acs.nanolett.4c01088
Abstract: Recent experiments have revealed that the macroscopic Kelvin equation remains surprisingly accurate even for nanoscale capillaries. This phenomenon was so far explained by the oscillatory behavior of the solid−liquid interfacial free energy. We here demonstrate thermodynamic and capillarity inconsistencies with this explanation. After revising the Kelvin equation, we ascribe its validity at nanoscale confinement to the effect of disjoining pressure.
To substantiate our hypothesis, we employed molecular dynamics simulations to evaluate interfacial heat transfer and wetting properties. Our assessments unveil a breakdown in a previously established proportionality between the work of adhesion and the Kapitza conductance at capillary heights below 1.3 nm, where the dominance of the work of adhesion shifts primarily from energy to entropy. Alternatively, the peak density of the initial water layer can effectively probe the work of adhesion. Unlike under bulk conditions, high confinement renders the work of adhesion entropically unfavorable.
Keywords: A1 Journal Article; CMT
Impact Factor: 10.8
DOI: 10.1021/acs.nanolett.4c01088
|
“Improving Molecule–Metal Surface Reaction Networks Using the Meta-Generalized Gradient Approximation: CO2Hydrogenation”. Cai Y, Michiels R, De Luca F, Neyts E, Tu X, Bogaerts A, Gerrits N, The Journal of Physical Chemistry C 128, 8611 (2024). http://doi.org/10.1021/acs.jpcc.4c01110
Abstract: Density functional theory is widely used to gain insights into molecule−metal surface reaction networks, which is important for a better understanding of catalysis. However, it is well-known that generalized gradient approximation (GGA)
density functionals (DFs), most often used for the study of reaction networks, struggle to correctly describe both gas-phase molecules and metal surfaces. Also, GGA DFs typically underestimate reaction barriers due to an underestimation of the selfinteraction energy. Screened hybrid GGA DFs have been shown to reduce this problem but are currently intractable for wide usage. In this work, we use a more affordable meta-GGA (mGGA) DF in combination with a nonlocal correlation DF for the first time to study and gain new insights into a catalytically important surface
reaction network, namely, CO2 hydrogenation on Cu. We show that the mGGA DF used, namely, rMS-RPBEl-rVV10, outperforms typical GGA DFs by providing similar or better predictions for metals and molecules, as well as molecule−metal surface adsorption
and activation energies. Hence, it is a better choice for constructing molecule−metal surface reaction networks.
Keywords: A1 Journal Article; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 3.7
DOI: 10.1021/acs.jpcc.4c01110
|