Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

Abstract: The presence of iodide ions during the growth of gold nanorods strongly affects the shape of the final products, which is proposed to be due to selective iodide adsorption on certain crystallographic facets. Therefore, a detailed structural and morphological characterization of the starting rods is crucial toward understanding this effect. Electron tomography is used to determine the crystallographic indices of the lateral facets of gold nanorods, as well as those present at the tips. Based on this information, density functional theory calculations are used to determine the surface and interface energies of the observed facets and provide insight into the relationship between the amount of iodide ions in the growth solution and the final morphology of anisotropic gold nanoparticles.

Abstract: Atomic scale simulations have been shown to be a powerful tool for elucidating the growth mechanisms of carbon nanotubes. The growth picture is however not entirely clear yet due to the gap between current simulations and real experiments. We here simulate for the first time the nucleation and subsequent growth of single-wall carbon nanotubes (SWNTs) from oxygen-containing hydrocarbon feedstocks using the hybrid Molecular Dynamics/Monte Carlo technique. The underlying nucleation mechanisms of Ni-catalysed SWNT growth are discussed in detail. Specifically, we find that as a function of the feedstock, different carbon fractions may emerge as the main growth species, due to a competition between the feedstock decomposition, its rehydroxylation and its contribution to etching of the growing SWNT. This study provides a further understanding of the feedstock effects in SWNT growth in comparison with available experimental evidence as well as with<italic>ab initio</italic>and other simulation data, thereby reducing the simulation–experiment gap.

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.

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.

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.

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.

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.

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(³P_J) + C₆H₆and C⁺(²P) + C₆H₆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(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆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(³P_J) + C₆H₆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.

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.

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.

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.