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: Optimization of catalytic materials for a given application is greatly constrained by linear scaling relations. Recently, however, it has been demonstrated that it is possible to reversibly modulate the chemisorption of molecules on nanomaterials by charging (i.e., injection or removal of electrons) and hence reversibly and selectively modify catalytic activity beyond structure−activity correlations. The fundamental physical relation between the properties of the material, the charging process, and the chemisorption energy, however, remains unclear, and a systematic exploration and optimization of charge-switchable sorbent materials is not yet possible. Using hybrid DFT calculations of CO2 chemisorption on hexagonal boron nitride nanosheets with several types of defects and dopants, we here reveal the existence of fundamental correlations between the electron affinity of a material and charge-induced chemisorption, show how defect engineering can be used to modulate the strength and efficiency of the adsorption process, and demonstrate that excess electrons stabilize many topological defects. We then show how these insights could be exploited in the development of new electrocatalytic materials and the synthesis of doped nanomaterials. Moreover, we demonstrate that calculated chemical properties of charged materials are highly sensitive to the employed computational methodology because of the self-interaction error, which underlines the theoretical challenge posed by such systems.

Abstract: In this paper we employ first-principles calculations to investigate the effect of substitutional Si doping in the amorphous calcium-phosphate (a-HAP) structure on the work of adhesion, integral charge transfer, charge density difference and theoretical tensile strengths between an a-HAP coating and amorphous titanium dioxide (a-TiO2) substrate systemically. Our calculations demonstrate that substitution of a P atom by a Si atom in a-HAP (a-Si-HAP) with the creation of OH-vacancies as charge compensation results in a significant increase of the bonding strength of the coating to the substrate. The work of adhesion of the optimized Si-doped interfaces reaches a value of up to -2.52 J m(-2), which is significantly higher than for the stoichiometric a-HAP/a-TiO2. Charge density difference analysis indicates that the dominant interactions at the interface have significant covalent character, and in particular two Ti-O and three Ca-O bonds are formed for a-Si-HAP/a-TiO2 and one Ti-O and three Ca-O bonds for a-HAP/a-TiO2. From the stress-strain curve, the Young's modulus of a-Si-HAP/a-TiO2 is calculated to be about 25% higher than that of the a-HAP/a-TiO2, and the yielding stress is about 2 times greater than that of the undoped model. Our calculations therefore demonstrate that the presence of Si in the a-HAP structure strongly alters not only the bioactivity and resorption rates, but also the mechanical properties of the a-HAP/a-TiO2 interface. The results presented here provide an important theoretical insight into the nature of the chemical bonding at the a-HAP/a-TiO2 interface, and are particularly significant for the practical medical applications of HAP-based biomaterials.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.

Abstract: The ever increasing global production and dispersion of methane requires novel chemistry to transform it into easily condensable energy carriers that can be integrated into the chemical infrastructure. In this context, single atom catalysts have attracted considerable interest due to their outstanding catalytic activity. We here use density functional theory (DFT) computations to compare the reaction and activation energies of M and MN4 embedded graphene (M = Ni and Si) on the methane-to-methanol conversion near room temperature. Thermodynamically, conversion of methane to methanol is energetically favorable at ambient conditions. Both singlet and triplet spin state of the studied systems are considered in all of the calculations. The DFT results show that the barriers are significantly lower when the complexes are in the triplet state than in the singlet state. In particular, Si-G with the preferred spin multiplicity of triplet seems to be viable catalysts for methane oxidation thanks to the corresponding lower energy barriers and higher stability of the obtained configurations. Our results provide insights into the nature of methane conversion and may serve as guidance for fabricating cost-effective graphene-based single atom catalysts.

Abstract: In this work, the effects of N-doping into the Co-doped single vacancy (Co-SV-G) and di-vacancy graphene flake (Co-dV-G) are investigated and compared toward direct oxidation of methane to methanol (DOMM) employing two different oxidants (N2O and O-2) using density functional theory (DFT) calculation. We found that DOMM on CoN3-G utilizing the N2O molecule as oxygen-donor proceeds via a two-step reaction with low activation energies. In addition, we found that although CoN3-G might be a good catalyst for methane conversion, it can also catalyze the oxidation of methanol to CO2 and H2O due to the required low activation barriers. Moreover, the adsorption behaviors of CHx (x = 0-4) species and dehydrogenation of CHx (x = 1-4) species on CoN3-G are investigated. We concluded that CoN3-G can be used as an efficient catalyst for DOMM and N-2O reduction at ambient conditions which may serve as a guide for fabricating effective C/N catalysts in energy-related devices.

Abstract: To gain full control over the growth of carbon nanotubes (CNTs) using plasma-enhanced chemical vapor deposition (PECVD), a thorough understanding of the underlying plasma-catalyst mechanisms is required. Oxygen-containing species are often used as or added to the growth precursor gas, but these species also yield various radicals and ions, which may simultaneously etch the CNT during the growth. At present, the effect of these reactive species on the growth onset has not yet been thoroughly investigated. We here report on the etching mechanism of incipient CNT structures from OH and O radicals as derived from combined (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations. Our results indicate that the oxygen-containing radicals initiate a dissociation process. In particular, we show how the oxygen species weaken the interaction between the CNT and the nanocluster. As a result of this weakened interaction, the CNT closes off and dissociates from the cluster in the form of a fullerene. Beyond the specific systems studied in this work, these results are generically important in the context of PECVD-based growth of CNTs using oxygen-containing precursors.

Abstract: Immobilization of single metal atoms on a solid host opens numerous possibilities for catalyst designs. If that host is a two-dimensional sheet, sheet curvature becomes a design parameter potentially complementary to host and metal composition. Here, we use a combination of density functional theory calculations and microkinetic modeling to compare the mechanisms and kinetics of formic acid decomposition and formation, chosen for their relevance as a potential hydrogen storage medium, over single Ru atoms anchored to pyridinic nitrogen in a planar graphene flake (RuN4-G) and curved carbon nanotube (RuN4-CNT). Activation barriers are lowered and the predicted turnover frequencies are increased over RuN4-CNT relative to RuN4-CNT. The results highlight the potential of curvature control as a means to achieve high performance and robust catalysts.

Abstract: Plasma-catalytic CO, hydrogenation for methanol production is gaining increasing interest, but our understanding of its reaction mechanism remains primitive. We present a combined experimental/computational study on plasma-catalytic CO, hydrogenation to CH3OH over a size-selected Cu/gamma-Al2O3 catalyst. Our experiments demonstrate a synergistic effect between the Cu/gamma-Al2O3 catalyst and the CO2/H-2 plasma, achieving a CO2 conversion of 10% at 4 wt % Cu loading and a CH3OH selectivity near 50% further rising to 65% with H2O addition (for a H2O/CO2 ratio of 1). Furthermore, the energy consumption for CH3OH production was more than 20 times lower than with plasma only. We carried out density functional theory calculations over a Cu-13/gamma-Al2O3 model, which reveal that the interfacial sites of the Cu-13 cluster and gamma-Al2O3 support show a bifunctional effect: they not only activate the CO2 molecules but also strongly adsorb key intermediates to promote their hydrogenation further. Reactive plasma species can regulate the catalyst surface reactions via the Eley-Rideal (E-R) mechanism, which accelerates the hydrogenation process and promotes the generation of the key intermediates. H2O can promote the CH3OH desorption by competitive adsorption over the Cu-13/gamma-Al2O3 surface. This study provides new insights into CO2 hydrogenation through plasma catalysis, and it provides inspiration for the conversion of some other small molecules (CH4, N-2, CO, etc.) by plasma catalysis using supported-metal clusters.

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: Although the heterogeneity of pyrolyzed Fe???N???C materials is known and has been reported previously, the atomic structure of the active sites and their detailed reaction mechanisms are still unknown. Here, we identified two pyridinic Fe???N4-like centers with different local C coordinates, i.e., FeN4C8 and FeN4C10, and studied their electrocatalytic activity for the nitrogen reduction reaction (NRR) based on density functional theory (DFT) calculations. We also discovered the influence of the adsorption of NH2 as a functional ligand on catalyst performance on the NRR. We confirmed that the NRR selectivity of the studied catalysts is essentially governed either by the local C coordination or by the dynamic structure associated with the FeII/FeIII. Our investigations indicate that the proposed traditional pyridinic FeN4C10 has higher catalytic activity and selectivity for the NRR than the robust FeN4C8 catalyst, while it may have outstanding activity for promoting other (electro)catalytic reactions. <comment>Superscript/Subscript Available</comment

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: New & beta;-stabilizedTi-based alloys are highly promising forbone implants, thanks in part to their low elasticity. The natureof this elasticity, however, is as yet unknown. We here present combinedfirst-principles DFT calculations and experiments on the microstructure,structural stability, mechanical characteristics, and electronic structureto elucidate this origin. Our results suggest that the studied & beta;Ti-35Nb-7Zr-5Ta wt % (TNZT) alloy manufacturedby the electron-beam powder bed fusion (E-PBF) method has homogeneousmechanical properties (H = 2.01 & PLUSMN; 0.22 GPa and E = 69.48 & PLUSMN; 0.03 GPa) along the building direction,which is dictated by the crystallographic texture and microstructuremorphologies. The analysis of the structural and electronic properties,as the main factors dominating the chemical bonding mechanism, indicatesthat TNZT has a mixture of strong metallic and weak covalent bonding.Our calculations demonstrate that the softening in the Cauchy pressure(C & PRIME; = 98.00 GPa) and elastic constant C ̅ ( 44 ) = 23.84 GPa is the originof the low elasticity of TNZT. Moreover, the nature of this softeningphenomenon can be related to the weakness of the second and thirdneighbor bonds in comparison with the first neighbor bonds in theTNZT. Thus, the obtained results indicate that a carefully designedTNZT alloy can be an excellent candidate for the manufacturing oforthopedic internal fixation devices. In addition, the current findingscan be used as guidance not only for predicting the mechanical propertiesbut also the nature of elastic characteristics of the newly developedalloys with yet unknown properties.

Abstract: Purpose Natural deep eutectic solvents (NADES) represent a green alternative to conventional organic solvents as reaction medium, offering more benign properties. To efficiently design NADES for biocatalysis, a better understanding of their effect on these reactions is needed. We hypothesize that this effect can be described by separately considering (1) the solvent interactions with the substrates, (2) the solvent viscosities and (3) the enzyme stability in NADES. Methods We investigated the effect of substrate solvation and viscosity on the reaction rate; and the stability of the enzyme in NADES. To this end, we monitored the conversion over time of the transesterification of vinyl laurate with 1- butanol by the lipase enzyme Candida antarctica B in NADES of different compounds and molar ratios. Results The initial reaction rate is higher in most NADES ( varying between 1.14 and 15.07 mu mol min(-1) mg(-1)) than in the reference n-hexane (4.0 mu mol min(-1) mg(-1))), but no clear relationship between viscosity and initial reaction rate was found. The increased reaction rate is most likely related to the solvation of the substrate due to a change in the activation energy of the reaction or a change in the conformation of the substrate. The enzyme retained part of its activity after the first 2 h of reaction (on average 20 % of the substrate reacted in the 2-24 h period). Enzyme incubation in ethylene glycol-based NADES resulted in a reduced reaction rate ( 15.07 vs. 3.34 mu mol min(-1) mg(-1)), but this may also be due to slow dissolution of the substrate. Conclusions The effect of viscosity seems to be marginal next to the effect of solvation and possible enzyme-NADES interaction. The enzyme retains some of its activity during the 24-hour measurements, but the enzyme incubation experiments did not yield accurate, comparable values. [GRAPHICS] .