Abstract: We developed a microkinetic model to reveal the effects of plasma-generated radicals, intermediates, and vibrationally excited species on the catalytic hydrogenation of CO2 to CH3OH on a Cu(111) surface. As a benchmark, we first present the mechanisms of thermal catalytic CH3OH formation. Our model predicts that the reverse water-gas shift reaction followed by CO hydrogenation, together with the formate path, mainly contribute to CH3OH formation in thermal catalysis. Adding plasma-generated radicals and intermediates results in a higher CH3OH turnover frequency (TOF) by six to seven orders of magnitude, showing the potential of plasma-catalytic CO2 hydrogenation into CH3OH, in accordance with the literature. In addition, CO2 vibrational excitation further increases the CH3OH TOF, but the effect is limited due to relatively low vibrational temperatures under typical plasma catalysis conditions. The predicted increase in CH3OH formation by plasma catalysis is mainly attributed to the increased importance of the formate path. In addition, the conversion of plasma-generated CO to HCO* and subsequent HCOO* or H2CO* formation contribute to CH3OH formation. Both pathways bypass the HCOO* formation from CO2, which is the main bottleneck in the process. Hence, our model points toward the important role of CO, but also O, OH, and H radicals, as they influence the reactions that consume CO2 and CO. In addition, our model reveals that the H pressure should not be smaller than ca. half of the O pressure in the plasma as this would cause O* poisoning, which would result in very small product TOFs. Thus, plasma conditions should be targeted with a high CO and H content as this is favorable for CH3OH formation, while the O content should be minimized.

Abstract: Assemblies of closely separated gold nanoparticles exhibit a strong collective plasmonic response due to coupling of the plasmon modes of the individual nanostructures. In the context of self-assembly of nanoparticles, close-packed two-dimensional (2D) clusters of spherical nanoparticles present an important composite system that promises numerous applications. The present study probes the collective plasmonic characteristics and resulting photothermal behavior of close-packed 2D Au nanoparticle clusters to delineate the effects of the cluster size, interparticle distance, and particle size. Smaller nanoparticles (20 and 40 nm in diameter) that exhibit low individual scattering and high absorption were considered for their relevance to photothermal applications. In contrast to typical literature studies, the present study compares the optical response of clusters of different sizes ranging from a single nanoparticle up to large assemblies of 61 nanoparticles. Increasing the cluster size induces significant changes to the spectral position and optophysical characteristics. Based on the model outcome, an optimal cluster size for maximum absorption per nanoparticle is also determined for enhanced photothermal effects. The effect of the particle size and interparticle distance is investigated to elucidate the nature of interaction in terms of near-field and far-field coupling. The photothermal effect resulting from absorption is compared for different cluster sizes and interparticle distances considering a homogeneous water medium. A strong dependence of the steady-state temperature of the nanoparticles on the cluster size, particle position in the cluster, incident light polarization, and interparticle distance provides new physical insight into the local temperature control of plasmonic nanostructures.

Abstract: Bimetallic plasmonic nanoparticles enable tuning of the optical response and chemical stability by variation of the composition. The present numerical simulation study compares Ag–Au alloy, Ag@Au core–shell, and Au@Ag core–shell bimetallic plasmonic nanoparticles of both spherical and anisotropic (nanotriangle and nanorods) shapes. By studying both spherical and anisotropic (with LSPR in the near-infrared region) shapes, cases with and without interband transitions of Au can be decoupled. Explicit comparisons are facilitated by numerical models supported by careful validation and examination of optical constants of Au–Ag alloys reported in the literature. Although both Au–Ag core–shell and alloy nanoparticles exhibit an intermediary optical response between that of pure Ag and Au nanoparticles, there are noticeable differences in the spectral characteristics. Also, the effect of the bimetallic constitution in anisotropic nanoparticles is starkly different from that in spherical nanoparticles due to the absence of Au interband transitions in the former case. In general, the improved chemical stability of Ag nanoparticles by incorporation of Au comes with a cost of reduction in plasmonic enhancement, also applicable to anisotropic nanoparticles with a weaker effect. A photothermal heat transfer study confirms that increased absorption by the incorporation of Au in spherical Ag nanoparticles also results in an increased steady-state temperature. On the other hand, anisotropic nanoparticles are inherently better absorbers and hence better photothermal sources, and their photothermal properties are apparently not strongly affected by the incorporation of one metal in the other. This study of the optical/spectral and photothermal characteristics of bimetallic Au–Ag alloy versus core–shell nanoparticles provides detailed physical insight for development of new taylor-made plasmonic nanostructures.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: The quartic force fields of FNO and CINO have been computed at the CCSD(T)/cc-pVTZ level. Using an ''augmented'' basis set dramatically improves results for FNO but has no significant effect for CINO. The best computed force field for FNO yields harmonic frequencies and fundamentals in excellent agreement with experiment. Overall, the force fields proposed in the present work are probably the most reliable ones ever published for these molecules. Total atomization energies have been computed using basis sets of spdfg quality: our best estimates are Sigma D-0 = 208.5 +/- 1 and 185.4 +/- 1 kcal/mol for FN0 and CINO, respectively. The computed value for FNO suggests a problem with the established experimental heat of formation. Thermodynamic tables in JANAF style at 100-2000 K are presented for both FNO and CINO.

Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: The gliding arc plasmatron (GAP) is a highly efficient atmospheric plasma source, which is very promising for CO2 conversion applications. To understand its operation principles and to improve its application, we present here comprehensive modeling results, obtained by means of computational fluid dynamics simulations and plasma modeling. Because of the complexity of the CO2 plasma, a full 3D plasma model would be computationally impractical. Therefore, we combine a 3D turbulent gas flow model with a 2D plasma and gas heating model in order to calculate the plasma parameters and CO2 conversion characteristics. In addition, a complete 3D gas flow and plasma model with simplified argon chemistry is used to evaluate the gliding arc evolution in space and time. The calculated values are compared with experimental data from literature as much as possible in order to validate the model. The insights obtained in this study are very helpful for improving the application of CO2 conversion, as they allow us to identify the limiting factors in the performance, based on which solutions can be provided on how to further improve the capabilities of CO2 conversion in the GAP.

Abstract: Doping two single transition-metal (TM) atoms on a substrate host opens numerous possibilities for catalyst design. However, what if the substrate contains more than one vacancy site? Then, the combination of two TMs along with their distribution patterns becomes a design parameter potentially complementary to the substrate itself and the bimetal composition. In this study, we investigate ammonia synthesis under mild electrocatalytic conditions on a transition-metal-doped porous C24N24 catalyst using density functional theory (DFT). The TMs studied include Ti, Mn, and Cu in a 2:4 dopant ratio (Ti2Mn4@C24N24 and Ti2Cu4@N-24(24)). Our computations show that a single Ti atom in both catalysts exhibits the highest selectivity for N-2 fixation at ambient conditions. This work is a good theoretical model to establish the structure-activity relationship, and the knowledge earned from the metal-N-4 moieties may help studies of related nanomaterials, especially those with curved structures.

Abstract: Supersonic flows provide a high thermodynamic
nonequilibrium, which is crucial for energy-efficient conversion of
CO 2 in microwave plasmas and are therefore of great interest.
However, the effect of the flow on the chemical reactions is poorly
understood. In this work, we present a combined flow and plasma
chemical kinetics model of a microwave CO 2 plasma in a Laval
nozzle setup. The effects of the flow field on the different dissociation
and recombination mechanisms, the vibrational distribution, and the
vibrational transfer mechanism are discussed. In addition, the effect
of experimental parameters, like position of power deposition, outlet
pressure, and specific energy input, on the CO 2 conversion and
energy efficiency is examined. The short residence time of the gas in
the plasma region, the shockwave, and the maximum critical heat,
and thus power, that can be added to the flow to avoid thermal
choking are the main obstacles to reaching high energy efficiencies.

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.