Abstract: CsPbI3 is a promising perovskite material for photovoltaic applications in its photoactive perovskite or black phase. However, the material degrades to a photovoltaically inactive or yellow phase at room temperature. Various mitigation strategies are currently being developed to increase the lifetime of the black phase, many of which rely on inducing strains in the material that hinder the black-to-yellow phase transition. Physical insight into how these strategies exactly induce strain as well as knowledge of the spatial extent over which these strains impact the material is crucial to optimize these approaches but is still lacking. Herein, we combine machine learning potential-based molecular dynamics simulations with our in silico strain engineering approach to accurately quantify strained large-scale atomic structures on a nanosecond time scale. To this end, we first model the strain fields introduced by atomic substitutions as they form the most elementary strain sources. We demonstrate that the magnitude of the induced strain fields decays exponentially with the distance from the strain source, following a decay rate that is largely independent of the specific substitution. Second, we show that the total strain field induced by multiple strain sources can be predicted to an excellent approximation by summing the strain fields of each individual source. Finally, through a case study, we illustrate how this additive character allows us to explain how complex strain fields, induced by spatially extended strain sources, can be predicted by adequately combining the strain fields caused by local strain sources. Hence, the strain additivity proposed here can be adopted to further our insight into the complex strain behavior in perovskites and to design strain from the atomic level onward to enhance their sought-after phase stability.

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: The synthesis of gold nanorods requires the presence of symmetry-breaking and shape-directing additives, among which bromide ions and quaternary ammonium surfactants have been reported as essential. As a result, hexadecyltrimethylammonium bromide (CTAB) has been selected as the most efficient surfactant to direct anisotropic growth. One of the difficulties arising from this selection is the low solubility of CTAB in water at room temperature, and therefore the seeded growth of gold nanorods is usually performed at 25 degrees C or above, which has restricted so far the analysis of kinetic effects derived from lower temperatures. We report a systematic study of the synthesis of gold nanorods from pentatwinned seeds using hexadecyltrimethylammonium chloride (CTAC) as the principal surfactant and a low concentration of bromide as shape-directing agent. Under these conditions, the synthesis can be performed at temperatures as low as 8 degrees C, and the corresponding kinetic effects can be studied, resulting in temperature-controlled aspect ratio tunability.

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: 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: A deep understanding of the fine structure at the atomic scale of halide perovskite materials has been limited by their sensitivity to the electron beam that is widely used for structural characterization. The sensitivity of a gamma-CsPbIBr2 perovskite thin film under electron beam irradiation is revealed by scanning transmission electron microscopy (STEM) through a universal large-range electron dose measurement, which is based on discrete single-electron events in the STEM mode. Our research indicates that the gamma-CsPbIBr2 thin film undergoes structural changes with increasing electron overall dose (e(-).A(-2)) rather than dose rate (e(-).A(-2).s(-1)), which suggests that overall dose is the key operative parameter. The electron beam-induced structural evolution of gamma-CsPbIBr2 is monitored by fine control of the electron beam dose, together with the analysis of high-resolution (S)TEM, diffraction, and energy-dispersive X-ray spectroscopy. Our results show that the gamma-CsPbIBr2 phase first forms an intermediate phase [e.g., CsPb(1-x)(IBr)((3-y))] with a superstructure of ordered vacancies in the pristine unit cell, while a fraction of Pb2+ is reduced to Pb-0. As the electron dose increases, Pb nanoparticles precipitate, while the remaining framework forms the Cs2IBr phase, accompanied by some amorphization. This work provides guidelines to minimize electron beam irradiation artifacts for atomic-resolution imaging on CsPbIBr2 thin films.

Abstract: Semiconductor monolayer transition metal dichalcogenides (TMDs) have brought a new paradigm by introducing optically addressable valley degree of freedom. Concomitantly, their high flexibility constitutes a unique platform that links optics to mechanics via valleytronics. With the intention to expedite the research in this direction, we investigated ten TMDs, namely MoS2, MoSe2, MoTe2, WS2, WSe2, WTe2, MoSSe, MoSeTe, WSSe, and WSeTe, which particularly includes their so-called janus types (JTMDs). First, we obtained their electronic band structures using regular and hybrid density functional theory (DFT) calculations in the presence of the spin-orbit coupling and biaxial or uniaxial strain. Our DFT results indicated that against the expectations based on their reported piezoelectric behavior, JTMDs typically interpolated between the standard band properties of the constituent TMDs without producing a novel feature. Next, by fitting to our DFT data we generated both spinless and spinful k center dot p parameter sets which are quite accurate over the K valley where the optical activity occurs. As an important application of this parametrization, we considered the circular and linear dichroism under strain. Among the studied (J)TMDs, WTe2 stood out with its largest linear dichroism under uniaxial strain because of its narrower band gap and large K valley uniaxial deformation potential. This led us to suggest WTe2 monolayer membranes for optical polarization-based strain measurements, or conversely, as strain tunable optical polarizers.

Abstract: We present an extension and revision of the spectroscopic and structural data of the mixed-stack charge-transfer (CT) crystal 3,3 ',5,5 '-tetramethylbenzidine-tetrafluorotetracyano-quinodimethane (TMB-TCNQF4), associated with new electric and dielectric measurements. Refinement of synchrotron structural data at low temperature has led to revise the previously reported C2/m structure. The revised structure is P2(1)/m, with two dimerized stacks per unit cell, and is consistent with the low temperature vibrational data. However, polarized Raman data in the low-frequency region also indicate that by increasing temperature above 200 K, the structure presents an increasing degree of disorder, mainly along the stack axis. X-ray diffraction data at room temperature have confirmed that the correct structure is P2(1)/ m -no phase transitions -but did not allow substantiating the presence of disorder. On the other hand, dielectric measurements have evidenced a typical relaxor ferroelectric behavior already at room temperature, with a peak in the real part of dielectric constant epsilon'(T,v) around 200 K and 0.1 Hz. The relaxor behavior is explained in terms of the presence of spin solitons separating domains of opposite polarity that yield to ferroelectric nanodomains. TMB-TCNQF(4) is confirmed to be a narrow-gap band semiconductor (Ea similar to 0.3 eV) with a room-temperature conductivity of similar to 10(-4) Omega(-1) cm(-1).

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: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

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 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: 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: Plasma is gaining increasing interest for N2 fixation, being a flexible, electricity-driven alternative for the current conventional fossil fuel-based N2 fixation processes. As the vibrational-induced dissociation of N2 is found to be an energy-efficient pathway to acquire atomic N for the fixation processes, plasmas that are in vibrational nonequilibrium seem promising for this application. However, an important challenge in using nonequilibrium plasmas lies in preventing vibrational−translational (VT) relaxation processes, in which vibrational energy crucial for N2 dissociation is lost to gas heating. We present here both experimental and modeling results for the vibrational and gas temperature in a microsecond-pulsed microwave (MW) N2 plasma, showing how power pulsing can suppress this unfavorable VT relaxation and achieve a maximal vibrational nonequilibrium. By means of our kinetic model, we demonstrate that pulsed plasmas take advantage of the long time scale on which VT processes occur, yielding a very pronounced nonequilibrium over the whole N2 vibrational ladder. Additionally, the effect of pulse parameters like the pulse frequency and pulse width are investigated, demonstrating that the advantage of pulsing to inhibit VT relaxation diminishes for high pulse frequencies (around 7000 kHz) and long power pulses (above 400 μs). Nevertheless, all regimes studied here demonstrate a clear vibrational nonequilibrium while only requiring a limited power-on time, and thus, we may conclude that a pulsed plasma seems very interesting for energyefficient vibrational excitation.

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.