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: 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 use microkinetic modeling to examine the potential of plasma-catalytic partial oxidation (POX) of CH4 as a promising new approach to produce oxygenates. We study how different plasma species affect POX of CH4 on the Pt(111) surface, and we discuss the associated kinetic and mechanistic changes. We discuss the effect of vibrationally excited CH4 and O2, as well as plasma-generated radicals and stable intermediates. Our results show that vibrational excitation enhances the turnover frequency (TOF) of catalytic CH4 dissociation and has good potential for improving the selectivities toward CH3OH, HCOOH, and C2 hydrocarbons. Nevertheless, when also considering plasma-generated radicals, we find that these species mainly govern the surface chemistry. Additionally, we find that plasma-generated radicals and stable intermediates enhance the TOFs of COx and oxygenates, increase the selectivity toward oxygenates, and make the formation of HCOOH more significant on Pt(111). We also briefly illustrate the potential impact of Eley−Rideal reactions that involve plasma-generated radicals. Finally, we reveal how various radicals affect the catalyst surface chemistry and we link this to the formation of different products. This allows us to make suggestions on how the plasma composition should be altered to improve the formation of desired products.

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: 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: 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: 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: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities.spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.

Abstract: This work contributes to new and complementary experimental viscosity data for blended amine mixtures of aqueous N-methyldiethanolamine + 2-amino-2-methyl-1-propanol (MDEA + AMP) solutions with and without CO2 at different temperatures and mass fractions. For the unloaded MDEA + AMP solutions, measurements were conducted with total amine mass fractions ranging from 0.30 to 0.60. In the case of CO2-loaded aqueous MDEA + AMP solutions, experiments were performed at CO2 loadings ranging from 0.11 to 0.80. Proposed correlations were used to represent viscosity at the unloaded and CO2-loaded solutions within experimental uncertainty.

Abstract: The chemical accessibility of the Celv oxidation state enables redox chemistry to be performed on the naturally coinagemetal -deficient phases CeM1-xSO (M = Cu, Ag). A metastable black compound with the PbFC1 structure type (space group P4/nmm: a = 3.8396(1) angstrom, c = 6.607(4) angstrom, V = 97.40(6) angstrom(3)) and a composition approaching CeSO is obtained by deintercalation of Ag from CeAg0.8SO. High-resolution transmission electron microscopy reveals the presence of large defect-free regions in CeSO, but stacking faults are also evident which can be incorporated into a quantitative model to account for the severe peak anisotropy evident in all the highresolution X-ray and neutron diffractograms of bulk CeSO samples; these suggest that a few percent of residual Ag remains. A strawcolored compound with the filled PbFCI (i.e., ZrSiCuAs- or HfCuSi2type) structure (space group P4/nmm: a = 3.98171(1) angstrom, c = 8.70913(5) angstrom, V = 138.075(1) angstrom 3) and a composition close to LiCeSO, but with small amounts of residual Ag, is obtained by direct reductive lithiation of CeAga8S0 or by insertion of Li into CeSO using chemical or electrochemical means. Computation of the band structure of pure, stoichiometric CeSO predicts it to be a Ce' compound with the 4f-states lying approximately 1 eV above the sulfide-dominated valence band maximum. Accordingly, the effective magnetic moment per Ce ion measured in the CeSO samples is much reduced from the value found for the Ce3+-containing LiCeSO, and the residual paramagnetism corresponds to the Ce3+ ions remaining due to the presence of residual Ag, which presumably reflects the difficulty of stabilizing Ce' in the presence of sulfide (S2-). Comparison of the behavior of CeCu0.8SO with that of CeCu0.8SO reveals much slower reaction kinetics associated with the Cu,_xS layers, and this enables intermediate CeCui LixSO phases to be isolated.

Abstract: Layered Li(M,Li)O2 (where M is a transition metal) ordered rock-salt-type structures are used in advanced metal-ion batteries as one of the best hosts for the reversible intercalation of Li ions. Besides the conventional redox reaction involving oxidation/reduction of the M cation upon Li extraction/insertion, creating oxygen-located holes because of the partial oxygen oxidation increases capacity while maintaining the oxidized oxygen species in the lattice through high covalency of the M–O bonding. Typical degradation mechanism of the Li(M,Li)O2 electrodes involves partially irreversible M cation migration toward the Li positions, resulting in gradual capacity/voltage fade. Here, using LiRhO2 as a model system (isostructural and isoelectronic to LiCoO2), for the first time, we demonstrate an intimate coupling between the oxygen redox and M cation migration. A formation of the oxidized oxygen species upon electrochemical Li extraction coincides with transformation of the layered Li1–xRhO2 structure into the γ-MnO2-type rutile–ramsdellite intergrowth LiyRh3O6 structure with rutile-like [1 × 1] channels along with bigger ramsdellite-like [2 × 1] tunnels through massive and concerted Rh migration toward the empty positions in the Li layers. The oxidized oxygen dimers with the O–O distances as short as 2.26 Å are stabilized in this structure via the local Rh–O configuration reminiscent to that in the μ-peroxo-μ-hydroxo Rh complexes. The LiyRh3O6 structure is remarkably stable upon electrochemical cycling illustrating that proper structural implementation of the oxidized oxygen species can open a pathway toward deliberate employment of the anion redox chemistry in high-capacity/high-voltage positive electrodes for metal-ion batteries. Upon chemical or electrochemical oxidation, layered LiRhO2 shows a unique structural transformation that involves both cation migration and oxidation of oxygen resulting in a stable tunnel-like rutile−ramsdellite intergrowth LiyRh3O6 structure. This structure demonstrates excellent performance with the steady and reversible capacity of ∼200 mAh/g. The stability of LiyRh3O6 is rooted in the accommodation of partially oxidized oxygen species through the formation of short O−O distances that are compatible with the connectivity of RhO6 octahedra.

Abstract: Atomic arrangements in the mixed-valence oxide U3O7 are refined from high-resolution neutron scattering data. The crystallographic model describes a long-range structural order in a U60O140 primitive cell (space group P42/n) containing distorted cuboctahedral oxygen clusters. By combining experimental data and electronic structure calculations accounting for spin–orbit interactions, we provide robust evidence of an interplay between charge localization and the magnetic moments carried by the uranium atoms. The calculations predict U3O7 to be a semiconducting solid with a band gap of close to 0.32 eV, and a more pronounced charge-transfer insulator behavior as compared to the well-known Mott insulator UO2. Most uranium ions (56 out of 60) occur in 9-fold and 10-fold coordinated environments, surrounding the oxygen clusters, and have a tetravalent (24 out of 60) or pentavalent (32 out of 60) state. The remaining uranium ions (4 out of 60) are not contiguous to the oxygen cuboctahedra and have a very compact, 8-fold coordinated environment with two short (2 × 1.93(3) Å) “oxo-type” bonds. The higher Hirshfeld charge and the diamagnetic character point to a hexavalent state for these four uranium ions. Hence, the valence state distribution corresponds to 24/60 × U(IV) + 32/60 U(V) + 4/60 U(VI). The tetravalent and pentavalent uranium ions are predicted to carry noncollinear magnetic moments (with amplitudes of 1.6 and 0.8 μB, respectively), resulting in canted ferromagnetic order in characteristic layers within the overall fluorite-related structure.