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: The reduction of dimensionality is a very effective way to achieve appealing properties in two-dimensional materials (2DMs). First-principles calculations can greatly facilitate the prediction of 2DM properties and find possible approaches to enhance their performance. We employed first-principles calculations to gain insight into the impact of different types of point defects (vacancies and substitutional dopants) on the electronic and magnetic properties of a ZnSnN2 (ZSN) monolayer. We show that Zn, Sn, and N + Zn vacancy-defected structures are p-type conducting, while the defected ZSN with a N vacancy is n-type conducting. For substitutional dopants, we found that all doped structures are thermally and energetically stable. The most stable structure is found to be B-doping at the Zn site. The highest work function value (5.0 eV) has been obtained for Be substitution at the Sn site. Li-doping (at the Zn site) and Be-doping (at the Sn site) are p-type conducting, while B-doping (at the Zn site) is n-type conducting. We found that the considered ZSN monolayer-based structures with point defects are magnetic, except those with the N vacancy defects and Be-doped structures. The ab initio molecular dynamics simulations confirm that all substitutionally doped and defected structures are thermally stable. Thus, our results highlight the possibility of tuning the magnetism in ZnSnN2 monolayers through defect engineering.

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: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).

Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

Abstract: We evaluate the importance of quantum effects on the sticking of H2 on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.

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: 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: Dislocations are always observed during crystal growth, and it is usually desirable to reduce the dislocation density in high-quality crystals. Here, the annihilation process of the 30 degrees Shockley partial dislocation pairs in CdTe is studied by first-principles calculations. We found that the dislocations can glide relatively easily due to the weak local bonding. Our systematic study of the slipping mechanism of the dislocations suggests that the energy barrier for the annihilation process is low. Band structure calculations reveal that the band bending caused by the charge transfer between the two dislocation cores depends on the core-core distance. A simple linear model is proposed to describe the mechanism of formation of the dislocation pair. More importantly, we demonstrate that hole injection can affect the core structure, increase the mobility, and eventually trigger a spontaneous mutual annihilation, which could be employed as a possible facile way to reduce the dislocation density.

Abstract: We present a first-principles investigation on the stability, electronic structure, and mechanical response of ultrathin heterostructures composed of single layers of InSe and SiGe. First, by performing total energy optimization and phonon calculations, we show that single layers of InSe and SiGe can form dynamically stable heterostructures in 12 different stacking types. Valence and conduction band edges of the heterobilayers form a type-I heterojunction having a tiny band gap ranging between 0.09 and 0.48 eV. Calculations on elastic-stiffness tensor reveal that two mechanically soft single layers form a heterostructure which is stiffer than the constituent layers because of relatively strong interlayer interaction. Moreover, phonon analysis shows that the bilayer heterostructure has highly Raman active modes at 205.3 and 43.7 cm(-1), stemming from the out-of-plane interlayer mode and layer breathing mode, respectively. Our results show that, as a stable type-I heterojunction, ultrathin heterobilayer of InSe/SiGe holds promise for nanoscale device applications.

Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.

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: Plasma technology is gaining increasing interest for CO2 conversion. The gas temperature in (and after) the plasma reactor largely affects the performance. Therefore, we examine the effect of cooling/quenching, during and after the plasma, on the CO2 conversion and energy efficiency, for typical “warm” plasmas, by means of chemical kinetics modeling. For plasmas at low specific energy input (SEI ∼ 0.5 eV/molecule), it is best to quench at the plasma end, while for high-SEI plasmas (SEI ∼ 4 eV/molecule), quenching at maximum conversion is better. For low-SEI plasmas, quenching can even increase the conversion beyond the dissociation in the plasma, known as superideal quenching. To better understand the effects of quenching at different plasma conditions, we study the dissociation and recombination rates, as well as the vibrational distribution functions (VDFs) of CO2, CO, and O2. When a high vibrational−translational (VT) nonequilibrium exists at the moment of quenching, the dissociation and recombination reaction rates both increase. Depending on the conversion degree at the moment of quenching, this can lead to a net increase or decrease of CO2 conversion. In general, however, and certainly for equilibrium plasmas at high temperature, quenching after the plasma helps prevent recombination reactions and clearly enhances the final CO2 conversion. We also investigate the effect of different quenching cooling rates on the CO2 conversion and energy efficiency. Finally, we compare plasma-based conversion to purely thermal conversion. For warm plasmas with typical temperatures of 3000−4000 K, the conversion is roughly thermal.

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.