Abstract: We report a high-throughput first-principles characterization of the structural, mechanical, electronic, and vibrational properties of tetragonal single-layer transition metal Xenes (t-TMXs). Our calculations revealed 22 dynamically, mechanically and chemically stable structures among the 96 possible free-standing layers present in the t-TMX family. As a fingerprint for their structural identification, we identified four characteristic Raman active phonon modes, namely three in-plane and one out-of-plane optical branches, with various intensities and frequencies depending on the material in question. Spin-polarized electronic calculations demonstrated that anti-ferromagnetic (AFM) metals, ferromagnetic (FM) metals, AFM semiconductors, and non-magnetic semiconductor materials exist within this family, evidencing the potential of t-TMXs for further use in multifunctional heterostructures.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

Abstract: Herein, we investigate the DDT (1-dodecanethiol) functionalization of exfoliated MoS2 by using experimental and theoretical tools. For the functionalization of MoS2, DDT treatment was incorporated into the conventional NMP (N-methyl pyrrolidone) exfoliation procedure. Afterward, it has been demonstrated that the functionalization process is successful through optical, morphological and theoretical analysis. The D, G and 2LA peaks seen in the Raman spectrum of exfoliated NMP-MoS2 particles, indicate the formation of graphitic species on MoS2 sheets. In addition, as the DDT ratio increases, the vacant sites on MoS2 sheets diminish. Moreover, at an optimized ratio of DDT-NMP, the maximum number of graphitic quantum dots (GQDs) is observed on MoS2 nanosheets. Specifically, the STEM and AFM data confirm that GQDs reside on the MoS2 nano-sheets and also that the particle size of the DDT-MoS2 is mostly fixed, while the NMP-MoS2 show many smaller and distributed sizes. The comparison of PL intensities of the NMP-MoS2 and DDT-MoS2 samples states a 10-fold increment is visible, and a 60-fold increment in NIR region photoluminescent properties. Moreover, our results lay out understanding and perceptions on the surface and edge chemistry of exfoliated MoS2 and open up more opportunities for MoS2 and GQD particles with broader applications.

Abstract: Graphene oxide (GO) is one of the most promising candidates for next generation of atomically thin membranes. Nevertheless, one of the major issues for real world application of GO membranes is their undesirable swelling in an aqueous environment. Recently, we demonstrated that generation of H3O+ and OH- ions (e.g., with an external electric field) in the interlayer gallery could impart aqueous stability to the layered GO membranes (A. Gogoi, ACS Appl. Mater. Interfaces, 2022, 14, 34946). This, however, compromises the water flux through the membrane. In this study, we report on reducing the GO nanosheets as a solution to this issue. With the reduction of the GO nanosheets, the water flux through the layered GO membrane initially increases and then decreases again beyond a certain degree of reduction. Here, two key factors are at play. Firstly, the instability of the H-bond network between water molecules and the GO nanosheets, which increases the water flux. Secondly, the pore size reduction in the interlayer gallery of the membranes, which decreases the water flux. We also observe a significant improvement in the salt rejection of the membranes, due to the dissociation of water molecules in the interlayer gallery. In particular, for the case of 10% water dissociation, the water flux through the membranes can be enhanced without altering its selectivity. This is an encouraging observation as it breaks the traditional tradeoff between water flux and salt rejection of a membrane.

Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: The vast biomedical potential of cold atmospheric pressure plasmas (CAPs) is governed by the formation of reactive species. These biologically active species are formed upon the interaction of CAPs with the surroundings. In biological milieu, water plays an essential role. The development of biomedical CAPs thus requires understanding of the sources of the reactive species in aqueous media exposed to the plasma. This is especially important in case of the COST RF plasma jet, which is developed as a reference microplasma system. In this work, we investigated the formation of the OH radicals, H atoms and H2O2 in aqueous solutions exposed to the COST plasma jet. This was done by combining experimental and modelling approaches. The liquid phase species were analysed using UV-Vis spectroscopy and spin trapping with hydrogen isotopes and electron paramagnetic resonance (EPR) spectroscopy. The discrimination between the species formed from the liquid phase and the gas phase molecules was performed by EPR and 1H-NMR analyses of the liquid samples. The concentrations of the reactive species in the gas phase plasma were obtained using a zero-dimensional (0D) chemical kinetics computational model. A three-dimensional (3D) fluid dynamics model was developed to provide information on the induced humidity in the plasma effluent. The comparison of the experimentally obtained trends for the formation of the species as a function of the feed gas and effluent humidity with the modelling results suggest that all reactive species detected in our system are mostly formed in the gas phase plasma inside the COST jet, with minor amounts arising from the plasma effluent humidity.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: Cold atmospheric plasma in contact with solutions has many applications, but its chemistry contains many unknowns such as the undescribed reactions with solutes. By combining experiments and modelling, we report the first direct demonstration of the reaction of chloride with oxygen atoms in aqueous solutions exposed to cold plasma.

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.