Interactions of Water with Pristine and Defective MoS2

Molybdenum disulfide (MoS₂) is a lamellar solid lubricant often used in aerospace applications because of its extremely low friction coefficient (∼0.01) in inert environments. The lubrication performance of MoS₂ is significantly impaired by exposure to even small amounts of water and oxygen, and the mechanisms behind this remain poorly understood. We use density functional theory calculations to study the binding of water on MoS₂ sheets with and without defects. In general, we find that pristine MoS₂ is slightly hydrophilic but that defects greatly increase the binding affinity for water. Intercalated water disrupts the crystal structure of bulk MoS₂ due to the limited space between lamellae (∼3.4 Å), and this leads to generally unfavorable adsorption, except in the cases where water molecules are located on the sites of sulfur vacancies. We also find that water adsorption is more favorable directly below a surface layer of MoS₂ compared to in the bulk.

Stick or Spill? Scaling Relationships for the Binding Energies of Adsorbates on Single-Atom Alloy Catalysts

Single-atom alloy catalysts combine catalytically active metal atoms, present as dopants, with the selectivity of coinage metal hosts. Determining whether adsorbates stick at the dopant or spill over onto the host is key to understanding catalytic mechanisms on these materials. Despite a growing body of work, simple descriptors for the prediction of spillover energies (SOEs), i.e., the relative stability of an adsorbate on the dopant versus the host site, are not yet available. Using Density Functional Theory (DFT) calculations on a large set of adsorbates, we identify the dopant charge and the SOE of carbon as suitable descriptors. Combining them into a linear surrogate model, we can reproduce DFT-computed SOEs within 0.06 eV mean absolute error. More importantly, our work provides an intuitive theoretical framework, based on the concepts of electrostatic interactions and covalency, that explains SOE trends and can guide the rational design of future single-atom alloy catalysts.

Exploration and validation of force field design protocols through QM-to-MM mapping

The scale of the parameter optimisation problem in traditional molecular mechanics force field construction means that design of a new force field is a long process, and sub-optimal choices made in the early stages can persist for many generations. We hypothesise that careful use of quantum mechanics to inform molecular mechanics parameter derivation (QM-to-MM mapping) should be used to significantly reduce the number of parameters that require fitting to experiment and increase the pace of force field development. Here, we design and train a collection of 15 new protocols for small, organic molecule force field derivation, and test their accuracy against experimental liquid properties. Our best performing model has only seven fitting parameters, yet achieves mean unsigned errors of just 0.031 g cm^-3 and 0.69 kcal mol^-1 in liquid densities and heats of vaporisation, compared to experiment. The software required to derive the designed force fields is freely available at https://github.com/qubekit/QUBEKit.

Two-Dimensional Metal-Organic Framework with Ultrahigh Water Stability for Separation of Acetylene from Carbon Dioxide and Ethylene

Highly selective separation and purification of acetylene (C2H2) from ethylene (C2H4) and carbon dioxide (CO2) are daunting challenges in light of their similar molecule sizes and physical properties. Herein, we report a two-dimensional (2D) stable metal-organic framework (MOF), NUM-11 ([Cu(Hmpba)2]·1.5DMF) (H2mpba = 4-(3,5-dimethyl-1H-pyrazol-4-yl)benzoic acid), with sql topology, stacked together through π-π interactions for efficient separation of C2H2 from C2H4 and CO2. The 2D-MOF material offers high hydrolytic stability and good purification capacity; especially, it could survive in water for 7 months, even longer. This stable MOF selectively captures C2H2 from mixtures containing C2H4 and CO2, as determined by adsorption isotherms. The ideal adsorbed solution theory selectivity calculations and transient breakthrough experiments were performed to verify the separation capacity. The low isosteric heat of NUM-11a (desolvated NUM-11) (18.24 kJ mol-1 for C2H2) validates the feasibility of adsorbent regeneration with low energy footprint consumption. Furthermore, Grand Canonical Monte Carlo simulations confirmed that the pore surface of the NUM-11 framework enabled preferential binding of C2H2 over C2H4 and CO2 via multiple C-H···O, C-H···π, and C-H···C interactions. This work provides some insights to prepare stable MOF materials toward the purification of C2H2, and the water-stable structure, low isosteric heat, and good cycling stability of NUM-11 make it very promising for practical industrial application.

Benchmarking various types of partial atomic charges for classical all-atom simulations of metal-organic frameworks

The density derived electrostatic and chemical (DDEC) approach for calculating the charges of atoms in a metal-organic framework (MOF) is considered to be the most accurate (yet computationally costly) one among many charge-assignment methods. Here, we conducted a comparative study on five different types of atomic partial charges (namely CM5, Mulliken, Qeq, EQeq and PACMOF) prepared for a subset of MOFs with affordable computational costs and benchmarked them with respect to the DDEC charges, which is particularly relevant because currently most databases lack MOFs with pre-calculated DDEC charges. To find a suitable charge type alternative to the DDEC approach, we statistically ranked the five charge types based on two metrics, the relative standard deviation of charges and relative dipole moment difference, based on which we provide general guidance as well as suggestions for specific MOFs according to bond polarity analyses. Finally, we recommend a possible and more accurate parametrization scheme for future studies.

Multi-Level Computational Screening of in Silico Designed MOFs for Efficient SO2 Capture

SO₂ presence in the atmosphere can cause significant harm to the human and environment through acid rain and/or smog formation. Combining the operational advantages of adsorption-based separation and diverse nature of metal-organic frameworks (MOFs), cost-effective separation processes for SO₂ emissions can be developed. Herein, a large database of hypothetical MOFs composed of >300,000 materials is screened for SO₂/CH₄, SO₂/CO₂, and SO₂/N₂ separations using a multi-level computational approach. Based on a combination of separation performance metrics (adsorption selectivity, working capacity, and regenerability), the best materials and the most common functional groups in those most promising materials are identified for each separation. The top bare MOFs and their functionalized variants are determined to attain SO₂/CH₄ selectivities of 62.4-16899.7, SO₂ working capacities of 0.3-20.1 mol/kg, and SO₂ regenerabilities of 5.8-98.5%. Regarding SO₂/CO₂ separation, they possess SO₂/CO₂ selectivities of 13.3-367.2, SO₂ working capacities of 0.1-17.7 mol/kg, and SO₂ regenerabilities of 1.9-98.2%. For the SO₂/N₂ separation, their SO₂/N₂ selectivities, SO₂ working capacities, and SO₂ regenerabilities span the ranges of 137.9-67,338.9, 0.4-20.6 mol/kg, and 7.0-98.6%, respectively. Besides, using breakdowns of gas separation performances of MOFs into functional groups, separation performance limits of MOFs based on functional groups are identified where bare MOFs (MOFs with multiple functional groups) tend to show the smallest (largest) spreads.

DFT Investigations of KTiOPO4Mx (M = K, Na, Li) Anodes for Alkali-Ion Battery

Properties of KTiOPO4M x (M = K, Na, and Li, x = 0.000 to 1.000) as anodes for potassium-ion battery (PIB), sodium-ion battery (SIB), and lithium-ion battery (LIB) are investigated by density functional theory (DFT) calculations. Our work reveals that electrochemical performance of KTiOPO4 as an anode for PIB is superior to that for SIB and LIB, in terms of average voltage and ion diffusion kinetics. The ab initio molecular dynamics (AIMD) simulations indicate that KTiOPO4M x anode exhibits high structural stability, and alkali ion intercalation contributes to accelerate ion diffusion during charging process. Particularly, the low activation energy 0.406 eV of K migration on surface KTP(210) obtained by the climbing-image nudged elastic band (Cl-NEB) method, which suggests a high-rate capability. The systematical comparison of the performance of KTiOPO4 as an anode for PIB, SIB, and LIB on the theoretical perspective clarifies that a large channel is not always promising for small radius ion intercalation and diffusion.

Atomistic Insight into the Hydration States of Layered Double Hydroxides

Effective protective coatings are an essential component of lightweight engineering materials in a large variety of applications as they ensure structural integrity of the base material throughout its whole service life. Layered double hydroxides (LDHs) loaded with corrosion inhibitors depict a promising approach to realize an active corrosion protection for aluminum and magnesium. In this work, we employed a combination of density functional theory and molecular dynamics simulations to gain a deeper understanding of the influence of intercalated water content on the structure, the stability, and the anion-exchange capacity of four different LDH systems containing either nitrate, carbonate, or oxalate as potential corrosion inhibiting agents or chloride as a corrosion initiator. To quantify the structural change, we studied the atom density distribution, radial distribution function, and orientation of the intercalated anions. Additionally, we determined the stability of the LDH systems by calculating their respective hydration energies, hydrogen-bonded network connected to the intercalated water molecules, as well as the self-diffusion coefficients of the intercalated anions to provide an estimate for the probability of their release after intercalation. The obtained computational results suggest that the hydration state of LDHs has a significant effect on their key properties like interlayer spacing and self-diffusion coefficients of the intercalated anions. Furthermore, we conclude from our simulation results that a high self-diffusion coefficient which is linked to the mobility of the intercalated anions is vital for its release via an anion-exchange mechanism and to subsequently mitigate corrosion reactions. Furthermore, the presented theoretical study provides a robust force field for the computer-assisted design of further LDH-based active anticorrosion coatings.

Role of Oxygen and Fluorine in Passivation of the GaSb(111) Surface Depending on Its Termination

The mechanism of the chemical bonding of oxygen and fluorine on the GaSb(111) surface depending on its termination is studied by the projector augmented-waves method within density functional theory. It is shown that on an unreconstructed (111) surface with a cation termination, the adsorption of fluorine leads to the removal of surface states from the band gap. The binding energy of fluorine on the cation-terminated surface in the most preferable Ga-T position is lower by ~0.4 eV than that of oxygen, but it is significantly lower (by ~0.8 eV) on the anion-terminated surface. We demonstrate that the mechanism of chemical bonding of electronegative adsorbates with the surface has an ionic–covalent character. The covalence of the O–Sb bond is higher than the F–Sb one, and it is higher than both O–Ga and F–Ga bonds. Trends in the change in the electronic structure of the GaSb(111) surface upon adsorption of fluorine and oxygen are discussed. It is found that an increase in the oxygen concentration on the Sb-terminated GaSb(111) surface promotes a decrease in the density of surface states in the band gap.

Property-activity relations of multifunctional reactive ensembles in cation-exchanged zeolites: a case study of methane activation on Zn2+-modified zeolite BEA

The reactivity theories and characterization studies for metal-containing zeolites are often focused on probing the metal sites. We present a detailed computational study of the reactivity of Zn-modified BEA zeolite towards C-H bond activation of the methane molecule as a model system that highlights the importance of representing the active site as the whole reactive ensemble integrating the extra-framework Zn_EF²⁺ cations, framework oxygens (O_F^2-), and the confined space of the zeolite pores. We demonstrate that for our model system the relationship between the Lewis acidity, defined by the probe molecule adsorption energy, and the activation energy for methane C-H bond cleavage performs with a determination coefficient R² = 0.55. This suggests that the acid properties of the localized extra-framework cations can be used only for a rough assessment of the reactivity of the cations in the metal-containing zeolites. In turn, studying the relationship between the activation energy and pyrrole adsorption energy revealed a correlation, with R² = 0.80. This observation was accounted for by the similarity between the local geometries of the pyrrole adsorption complexes and the transition states for methane C-H bond cleavage. The inclusion of a simple descriptor for zeolite local confinement allows transferability of the obtained property-activity relations to other zeolite topologies. Our results demonstrate that the representation of the metal cationic species as a synergistically cooperating active site ensembles allows reliable detection of the relationship between the acid properties and reactivity of the metal cation in zeolite materials.

Atomistic Studies on Water-Induced Lithium Corrosion

It is well known that lithium reacts violently with water under the release of molecular hydrogen and the formation of lithium hydroxide. In this work, the initial mechanisms for the surface reactions of metallic lithium with water from the gas phase were investigated by means of periodic density functional theory calculations. For this purpose, adsorption/absorption structures and diffusion and dissociation processes of hydrogen, OH, and H₂ O on low-index metallic lithium surfaces were investigated. Through thermodynamic and kinetic considerations, negatively charged centers on the surface were identified as the origin of hydrogen formation. The strikingly low reaction barriers for the reaction at these centers implied a self-supporting effect of hydrogen evolution and the associated lithium degradation.

Studying the adsorption mechanisms of nanoplastics on covalent organic frameworks via molecular dynamics simulations

Covalent organic frameworks (COFs) with well-defined supramolecular structures and high surface-area-to-volume ratio have received extensive attention on their adsorption of contaminants from micro- to nano-size. Here, we studied the adsorption mechanisms of three typical nanoplastics (NP), including polyethylene (PE), nylon-6 (PA 6), and polyethylene terephthalate (PET) on chemically stable COFs (TpPa-X, X = H, CH₃, OH, NO₂ and F) by molecular dynamics simulations. Depending on molecular structure and surface composition, two distinct interactions-electrostatic interaction and van der Waals (vdW) interaction-are identified to be responsible for the adsorption of different NP pollutants on TpPa-X. The vdW interaction is dominant during the adsorption process, while polar groups in polymers and COFs can enhance the adsorption because of the electrostatic interaction. Compared with other functional COFs, we found that TpPa-OH shows the strongest adsorption with the NP pollutants employed in this study. This work reveals the COF-polymer adsorption behavior and properties at atomic scale, which is crucial to the development of promising COF materials to deal with NP pollution.

Assessing doping strategies for monolayer MoS2 towards non-enzymatic detection of cortisol: a first-principles study

In this work, we investigate by means of atomistic density functional theory simulations the interaction between cortisol (the target molecule) and monolayer MoS₂ (the substrate). The aim is to assess viable strategies for the non-enzymatic chemical sensing of cortisol. Metal doping of the sensing material could offer a way to improve the device response upon analyte adsorption, and could also enable novel and alternative detection mechanisms. For such reasons, we explore metal doping of MoS₂ with Ni, Pd, and Pt, as these are metal elements commonly used in experiments. Then, we study the material response from the structural, electronic, and charge-transfer points of view. Based on our results, we propose two possible sensing mechanisms and device architectures: (i) a field-effect transistor, and (ii) an electrochemical sensor. In the former, Ni-doped MoS₂ would act as the FET channel, and the sensing mechanism involves the variation of the surface electrostatic charge upon the adsorption of cortisol. In the latter, MoS₂ decorated with Pt nanoparticles could act as the working electrode, and the sensing mechanism would involve the reduction of cortisol. In addition, our findings may suggest the suitability of both doped and metal-doped MoS₂ as sensing layers in an optical sensor.

Apples to apples comparison of standardized to unstandardized principal component analysis of methods that assign partial atomic charges in molecules

Articles by Cho et al. (ChemPhysChem, 2020, 21, 688–696) and Manz (RSC Adv., 2020, 10, 44121–44148) performed unstandardized and standardized, respectively, principal component analysis (PCA) to study atomic charge assignment methods for molecular systems. Both articles used subsets of atomic charges computed by Cho et al.; however, the data subsets employed were not strictly identical. Herein, an element by element analysis of this dataset is first performed to compare the spread of charge values across individual chemical elements and charge assignment methods. This reveals an underlying problem with the reported Becke partial atomic charges in this dataset. Due to their unphysical values, these Becke charges were not included in the subsequent PCA. Standardized and unstandardized PCA are performed across two datasets: (i) 19 charge assignment methods having a complete basis set limit and (ii) all 25 charge assignment methods (excluding Becke) for which Cho et al. computed atomic charges. The dataset contained ∼2000 molecules having a total of 29 907 atoms in materials. The following five methods (listed here in alphabetical order) showed the greatest correlation to the first principal component in standardized and unstandardized PCA: DDEC6, Hirshfeld-I, ISA, MBIS, and MBSBickelhaupt (note: MBSBickelhaupt does not appear in the 19 methods dataset). For standardized PCA, the DDEC6 method ranked first followed closely by MBIS. For unstandardized PCA, Hirshfeld-I (19 methods) or MBSBickelhaupt (25 methods) ranked first followed by DDEC6 in second place (both 19 and 25 methods).

Standardized and unstandardized principal component analyses are directly compared for a dataset of ∼2000 molecules across various charge assignment methods.

Experimental Evidence for the Molecular Molybdenum Fluorides MoF to MoF6: a matrix isolation and DFT investigation

All of the molecular molybdenum fluorides, MoF to MoF6, have been synthesised from the reaction of molybdenum atoms with fluorine molecules and atoms, trapped in argon matrices, and characterised by...

Neural-network-backed evolutionary search for SrTiO3(110) surface reconstructions

The determination of atomic structures in surface reconstructions has typically relied on structural models derived from intuition and domain knowledge. Evolutionary algorithms have emerged as powerful tools for such structure searches. However, when density functional theory is used to evaluate the energy the computational cost of a thorough exploration of the potential energy landscape is prohibitive. Here, we drive the exploration of the rich phase diagram of TiO_x overlayer structures on SrTiO₃(110) by combining the covariance matrix adaptation evolution strategy (CMA-ES) and a neural-network force field (NNFF) as a surrogate energy model. By training solely on SrTiO₃(110) 4×1 overlayer structures and performing CMA-ES runs on 3×1, 4×1 and 5×1 overlayers, we verify the transferability of the NNFF. The speedup due to the surrogate model allows taking advantage of the stochastic nature of the CMA-ES to perform exhaustive sets of explorations and identify both known and new low-energy reconstructions.

The covariance matrix adaptation evolution strategy (CMA-ES) and a fully automatically differentiable, transferable neural-network force field are combined to explore TiO_x overlayer structures on SrTiO₃(110) 3×1, 4×1 and 5×1 surfaces.

Incorporating Flexibility Effects into Metal-Organic Framework Adsorption Simulations Using Different Models

High-throughput calculations based on molecular simulations to predict the adsorption of molecules inside metal-organic frameworks (MOFs) have become a useful complement to experimental efforts to identify promising adsorbents for chemical separations and storage. For computational convenience, all existing efforts of this kind have relied on simulations in which the MOF is approximated as rigid. In this paper, we use extensive adsorption-relaxation simulations that fully include MOF flexibility effects to explore the validity of the rigid framework approximation. We also examine the accuracy of several approximate methods to incorporate framework flexibility that are more computationally efficient than adsorption-relaxation calculations. We first benchmark various models of MOF flexibility for four MOFs with well-established CO₂ experimental consensus isotherms. We then consider a range of adsorption properties, including Henry's constants, nondilute loadings, and adsorption selectivity, for seven adsorbates in 15 MOFs randomly selected from the CoRE MOF database. Our results indicate that in many MOFs adsorption-relaxation simulations are necessary to make quantitative predictions of adsorption, particularly for adsorption at dilute concentrations, although more standard calculations based on rigid structures can provide useful information. Finally, we investigate whether a correlation exists between the elastic properties of empty MOFs and the importance of including framework flexibility in making accurate predictions of molecular adsorption. Our results did not identify a simple correlation of this type.

Parametrization of Nonbonded Force Field Terms for Metal-Organic Frameworks Using Machine Learning Approach

The enormous structural and chemical diversity of metal-organic frameworks (MOFs) forces researchers to actively use simulation techniques as often as experiments. MOFs are widely known for their outstanding adsorption properties, so a precise description of the host-guest interactions is essential for high-throughput screening aimed at ranking the most promising candidates. However, highly accurate ab initio calculations cannot be routinely applied to model thousands of structures due to the demanding computational costs. Furthermore, methods based on force field (FF) parametrization suffer from low transferability. To resolve this accuracy-efficiency dilemma, we applied a machine learning (ML) approach: extreme gradient boosting. The trained models reproduced the atom-in-material quantities, including partial charges, polarizabilities, dispersion coefficients, quantum Drude oscillator, and electron cloud parameters, with accuracy similar to the reference data set. The aforementioned FF precursors make it possible to thoroughly describe noncovalent interactions typical for MOF-adsorbate systems: electrostatic, dispersion, polarization, and short-range repulsion. The presented approach can also readily facilitate hybrid atomistic simulation/ML workflows.

Infrared Activities of Adsorbed Species on Metal Surfaces: The Puzzle of Adsorbed Methyl (CH3)

Reflection-absorption infrared spectroscopy (RAIRS) is widely used to identify molecular adsorbates on metals during surface chemical reactions, but the interpretation of RAIRS data can be difficult with experiment alone. Here, we reveal from first-principles calculations the origin of the contrasting RAIRS spectra of methyl adsorbed on Pt(111) and Ni(111). We find that the dynamic dipole associated with the symmetric C-H stretch vibration of CH₃ along surface normal is significant on Pt(111) but negligibly small on Ni(111), explaining the strong IR activity in the former and the absence of any RAIRS peaks in the latter. This difference is correlated to different charge transfer patterns between metals and the adsorbate, which are determined by the different preferred adsorption sites of methyl on the two surfaces. This work highlights the need of electronic structure calculations in interpreting RAIRS spectra of adsorbates on metal surfaces.

The Catalytic Mechanism of Pdx2 Glutaminase Driven by a Cys-His-Glu Triad: A Computational Study

The catalytic mechanism of Pdx2 was studied with atomic detail employing the computational ONIOM hybrid QM/MM methodology. Pdx2 employs a Cys-His-Glu catalytic triad to deaminate glutamine to glutamate and ammonia - the source of the nitrogen of pyridoxal 5'-phosphate (PLP). This enzyme is, therefore, a rate-limiting step in the PLP biosynthetic pathway of Malaria and Tuberculosis pathogens that rely on this mechanism to obtain PLP. For this reason, Pdx2 is considered a novel and promising drug target to treat these diseases. The results obtained show that the catalytic mechanism of Pdx2 occurs in six steps that can be divided into four stages: (i) activation of Cys₈₇ , (ii) deamination of glutamine with the formation of the glutamyl-thioester intermediate, (iii) hydrolysis of the formed intermediate, and (iv) enzymatic turnover. The kinetic data available in the literature (19.1-19.5 kcal mol^-1 ) agree very well with the calculated free energy barrier of the hydrolytic step (18.2 kcal.mol^-11 ), which is the rate-limiting step of the catalytic process when substrate is readily available in the active site. This catalytic mechanism differs from other known amidases in three main points: i) it requires the activation of the nucleophile Cys₈₇ to a thiolate; ii) the hydrolysis occurs in a single step and therefore does not require the formation of a second tetrahedral reaction intermediate, as it is proposed, and iii) Glu₁₉₈ does not have a direct role in the catalytic process. Together, these results can be used for the synthesis of new transition state analogue inhibitors capable of inhibiting Pdx2 and impair diseases like Malaria and Tuberculosis.

Electromechanical Characteristics by a Vertical Flip of C70 Fullerene Prolate Spheroid in a Single-Electron Transistor: Hybrid Density Functional Methods

In this study, the B3LYP hybrid density functional theory was used to investigate the electromechanical characteristics of C70 fullerene with and without point charges to model the effect of the surface of the gate electrode in a C70 single-electron transistor (SET). To understand electron tunneling through C70 fullerene species in a single-C70 transistor, descriptors of geometrical atomic structures and frontier molecular orbitals were analyzed. The findings regarding the node planes of the lowest unoccupied molecular orbitals (LUMOs) of C70 and both the highest occupied molecular orbitals (HOMOs) and the LUMO of the C70 anion suggest that electron tunneling of pristine C70 prolate spheroidal fullerene could be better in the major axis orientation when facing the gate electrode than in the major (longer) axis orientation when facing the Au source and drain electrodes. In addition, we explored the effect on the geometrical atomic structure of C70 by a single-electron addition, in which the maximum change for the distance between two carbon sites of C70 is 0.02 Å.

Identifying the Critical Surface Descriptors for the Negative Slopes in the Adsorption Energy Scaling Relationships via Density Functional Theory and Compressed Sensing

Adsorption energy scaling relationships have progressed beyond their original form, which was primarily focused on optimizing catalytic sites and lowering computational costs in simulations. The recent rise in interest in adsorption energy scaling relations is to investigate surfaces other than transition metals (TMs) as well as interactions involving complex compounds. In this work, we report our extensive study on the scaling relation (SR) between oxygen (O), with elements of neighboring groups such as boron (B), aluminum (Al), carbon (C), silicon (Si), nitrogen (N), phosphorus (P), and fluorine (F) on magnetic bimetallic surfaces. We observed that only O versus N and F seems to have a positive slope; the other slopes are negative. We present new theoretical model in terms of multiple surface descriptors using density functional theory and compressed sensing, whereas the original scaling theory was based on a single adsorbate descriptor: adsorbate valency.

Effects of external electric field on the sensing property of volatile organic compounds over Janus MoSSe monolayer: a first-principles investigation

Janus 2D transition metal dichalcogenide (TMD) is a new generation 2D material with a unique asymmetric structure. This asymmetric structure (or out-off plane symmetric geometry) of Janus 2D TMD has been reported to yield tunable electronic properties through strain and electric field, which can also be applied in gas sensing. In this work, we performed DFT calculations to investigate the gas sensing property of cyclohexane and acetone on MoS2 and Janus MoSSe monolayers under external electric fields. Our results show that cyclohexane possesses slightly larger adsorption energy on pristine MoS2 and Janus MoSSe monolayers than acetone without external electric fields. After applying the external electric fields, the adsorption energy for cyclohexane on MoS2 shows no enhancement. However, the adsorption energy of acetone shows the most substantial enhancement on the Janus MoSSe monolayer. We found that the dipole moment orientations of adsorbates and the monolayer can strongly interact with the external electric fields. Hence, the combination of polar adsorbate and polar material, i.e., acetone and Janus MoSSe, demonstrates the most vital sensitivity under the applied bias. On the other hand, the non-polar adsorbate and non-polar material combination show a negligible effect on external bias. These findings can be applied to the design of gas sensors in the future through polar materials.

First-Principles Density Functional Theory Calculations of Bilayer Membranes Heterostructures of Ti3C2T2 (MXene)/Graphene and AgNPs

The properties of two-dimensional (2D) layered membrane systems can be medullated by the stacking arrangement and the heterostructure composition of the membrane. This largely affects the performance and stability of such membranes. Here, we have used first-principle density functional theory calculations to conduct a comparative study of two heterostructural bilayer systems of the 2D-MXene (Ti3C2T2, T = F, O, and OH) sheets with graphene and silver nanoparticles (AgNPs). For all considered surface terminations, the binding energy of the MXene/graphene and MXene/AgNPs bilayers increases as compared with graphene/graphene and MXene/MXene bilayer structures. Such strong interlayer interactions are due to profound variations of electrostatic potential across the layers. Larger interlayer binding energies in MXene/graphene systems were obtained even in the presence of water molecules, indicating enhanced stability of such a hybrid system against delamination. We also studied the structural properties of Ti3C2X2 MXene (X = F, O and OH) decorated with silver nanoclusters Agn (n ≤ 6). We found that regardless of surface functionalization, Ag nanoclusters were strongly adsorbed on the surface of MXene. In addition, Ag nanoparticles enhanced the binding energy between MXene layers. These findings can be useful in enhancing the structural properties of MXene membranes for water purification applications.

Modeling Molecular Emitters in Organic Light-Emitting Diodes with the Quantum Mechanical Bespoke Force Field

Combined molecular dynamics (MD) and quantum mechanics (QM) simulation procedures have gained popularity in modeling the spectral properties of functional organic molecules. However, the potential energy surfaces used to propagate long-time scale dynamics in these simulations are typically described using general, transferable force fields designed for organic molecules in their electronic ground states. These force fields do not typically include spectroscopic data in their training, and importantly, there is no general protocol for including changes in geometry or intermolecular interactions with the environment that may occur upon electronic excitation. In this work, we show that parameters tailored for thermally activated delayed fluorescence (TADF) emitters used in organic light-emitting diodes (OLEDs), in both their ground and electronically excited states, can be readily derived from a small number of QM calculations using the QUBEKit (QUantum mechanical BEspoke toolKit) software and improve the overall accuracy of these simulations.

Molecular Sieving of Acetylene from Ethylene in a Rigid Ultra-microporous Metal Organic Framework

Rigid molecular sieving materials are the ideal candidates for gas separation (e. g., C₂ H₂ /C₂ H₄ ) due to their ultrahigh adsorption selectivity and the absence of gas co-adsorption. However, the absolute molecular sieving effect for C₂ H₂ /C₂ H₄ separation has rarely been realized because of their similar physicochemical properties. Herein, we demonstrate the absolute molecular sieving of C₂ H₂ from C₂ H₄ by a rigid ultra-microporous metal-organic framework (F-PYMO-Cu) with 1D regular channels (pore size of ca. 3.4 Å). F-PYMO-Cu exhibited moderate acetylene uptake (35.5 cm³ /cm³ ), but very low ethylene uptake (0.55 cm³ /cm³ ) at 298 K and 1 bar, yielding the second highest C₂ H₂ /C₂ H₄ uptake ratio of 63.6 up to now. One-step C₂ H₄ production from a binary mixture of C₂ H₂ /C₂ H₄ and a ternary mixture of C₂ H₂ /CO₂ /C₂ H₄ at 298 K was achieved and verified by dynamic breakthrough experiments. Coupled with excellent thermal and water stability, F-PYMO-Cu could be a promising candidate for industrial C₂ separation tasks.

Impact of interstitial impurities on the trapping of dislocation loops in tungsten

Ab initio simulations are employed to assess the interaction of typical interstitial impurities with self-interstitial atoms, dislocation loops and edge dislocation lines in tungsten. These impurities are present in commercial tungsten grades and are also created as a result of neutron transmutation or the plasma in-take process. The relevance of the study is determined by the application of tungsten as first wall material in fusion reactors. For the defects with dislocation character, the following ordering of the interaction strength was established: H < N < C < O < He. The magnitude of the interaction energy was rationalized by decomposing it into elastic (related to the lattice strain) and chemical (related to local electron density) contributions. To account for the combined effect of impurity concentration and pinning strength, the impact of the presence of these impurities on the mobility of isolated dislocation loops was studied for DEMO relevant conditions in the non-elastic and dilute limit.

Electronic and vibrational properties of the highTcsuperconductor Bi2Sr2CaCu2O8: anab initiostudy

In this work,ab initiocalculations were performed in order to study the vibrational spectra of the Bi₂Sr₂CaCu₂O₈(Bi2212) compound. A structural modulation correction on some atomic positions, producing a distorted structure with lower symmetry, is used for the calculation. We argue that this correction allows to account for an average effect of the incommensurate superstructure, generating a more accurate representation of the real unit cell observed in this compound. A complete and conclusive vibrational assignment is performed, discussing the correspondences with previous experimental and theoretical reports. A brief analysis of the electronic density of states and band structure comparing the tetragonal and distorted unit cell is also included.

Phase Transformation Dynamics in Sulfate-Loaded Lanthanide Triphosphonates. Proton Conductivity and Application as Fillers in PEMFCs

Phase transformation dynamics and proton conduction properties are reported for cationic layer-featured coordination polymers derived from the combination of lanthanide ions (Ln3+) with nitrilo-tris(methylenephosphonic acid) (H6NMP) in the presence of sulfate ions. Two families of materials are isolated and structurally characterized, i.e., [Ln2(H4NMP)2(H2O)4](HSO4)2·nH2O (Ln = Pr, Nd, Sm, Eu, Gd, Tb, Er, Yb; n = 4-5, Series I) and [Ln(H5NMP)]SO4·2H2O (Ln = Pr, Nd, Eu, Gd, Tb; Series II). Eu/Tb bimetallic solid solutions are also prepared for photoluminescence studies. Members of families I and II display high proton conductivity (10-3 and 10-2 S·cm-1 at 80 °C and 95% relative humidity) and are studied as fillers for Nafion-based composite membranes in PEMFCs, under operating conditions. Composite membranes exhibit higher power and current densities than the pristine Nafion membrane working in the range of 70-90 °C and 100% relative humidity and with similar proton conductivity.

Partial charges for molecular-mechanical models of heterocyclic compounds: pyridine nitrogen

This paper deals with the effect of introducing an additional interaction site onto molecular-mechanical models of nitrogen-containing heterocyclic compounds. The introduction of only one additional site next to nitrogen atoms is shown to result in significant improvement of the quality of the models along with negligible slowdown of calculation speed. Concretely, it was proposed to introduce the site inside the aromatic ring at a distance of 0.4 Å from the nitrogen atom center. All the parametrization can be completely automated. The proposed force field allows predicting heats of evaporation of liquids of the compound under investigation with an accuracy of 1 kcal/mol.

Reaction pathways in the solid state and the Hubbard U correction

We investigate how the Hubbard U correction influences vacancy defect migration barriers in transition metal oxide semiconductors. We show that, depending on the occupation of the transition metal d orbitals, the Hubbard U correction can cause severe instabilities in the migration barrier energies predicted using generalized gradient approximation density functional theory (GGA DFT). For the d⁰ oxide SrTiO₃, applying a Hubbard correction to the Ti⁴⁺ 3d orbitals below 4-5 eV yields a migration barrier of ∼0.4 eV. However, above this threshold, the barrier increases suddenly to ∼2 eV. This sudden increase in the transition state barrier arises from the Hubbard U correction changing the Ti⁴⁺ t_2g/e_g orbital occupation, and hence electron density localization, along the migration pathway. Similar results are observed in the d¹⁰ oxide ZnO; however, significantly larger Hubbard U corrections must be applied to the Zn²⁺ 3d orbitals for the same instability to be observed. These results highlight important limitations to the application of the Hubbard U correction when modeling reactive pathways in solid state materials using GGA DFT.

Interface Interaction of Benzohydroxamic Acid with Lead Ions on Oxide Mineral Surfaces: A Coordination Mechanism Study

Surface coordination chemistry is important in areas such as adsorption, separation, and catalysts. In this work, surface coordination interactions of benzohydroxamic acid (BHA) with the lead ion [Pb(II)] adsorbed on the cassiterite surface have been investigated by first-principles calculations due to its great significance in froth flotation. Cluster calculations show that BHA possesses the weakest chelation with Pb(II) due to the electron withdrawal ability of the benzyl ring in comparison with other hydroxamic acids. Pb(II) thermodynamically prefers to react with the cassiterite surface rather than BHA. On the other hand, the partial density of states and the atomic overlap populations have consistently verified that the adsorption of BHA results in a better symmetry in electron densities than the hydrated Pb(II). The electron density maps and the electronic localization functions have further visualized the rearrangement of the 6s² lone pair around the lead atom. It can be concluded that the surface coordination mechanisms of Pb(II) on oxide minerals can be attributed to the coordination ability of BHA and the unique electronic structure of Pb(II), which accounts for the reported better flotation performance of the pre-assemble strategy than the pre-activating approach. This work sheds some new light on the unique coordination activation mechanism of metal ions on oxide mineral surfaces. It should be instructive to design and screen new environment-friendly flotation reagents and flotation flowsheets.

The role of sodium in stabilizing tin-lead (Sn-Pb) alloyed perovskite quantum dots

Narrow-bandgap CsSn x Pb1-x I3 perovskite quantum dots (QDs) show great promise for optoelectronic applications owing to their reduced use of toxic Pb, improved phase stability, and tunable band gaps in the visible and near-infrared range. The use of small ions has been proven beneficial in enhancing the stability and photoluminescence quantum yield (PLQY) of perovskite QDs. The introduction of sodium (Na) has succeeded in boosting the PLQY of CsSn0.6Pb0.4I3 QDs. Unfortunately, the initial PLQY of the Na-doped QDs undergoes a fast degradation after one-day storage in solution, hindering their practical applications. Using density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations, we study the effect of Na ions on the strength of surface bonds, defect formation energies, and the interactions between surface ligands and perovskite QDs. Our results suggest that Na ions enhance the covalent bonding of surface tin-iodine bonds and form strong ionic bonding with the neighboring iodine anions, thus suppressing the formation of I and Sn vacancies. Furthermore, Na ions also enhance the binding strength of the surface ligands with the perovskite QD surface. However, according to our AIMD simulations, the enhanced surface ligand binding is only effective on a selected surface configuration. While the position of Na ions remains intact on a CsI-terminated surface, they diffuse vigorously on an MI2-terminated surface. As a result, the positive effect of Na vanishes with time, explaining the relatively short lifetime of the experimentally obtained high PLQYs. Our results indicate that engineering the surface termination of the QDs could be the next step in maintaining the favorable effect of Na doping for a high and stable PLQY of Sn-Pb QDs.

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal-Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal-organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO₂. Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal-Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. A strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

Phenylalkylammonium passivation enables perovskite light emitting diodes with record high-radiance operational lifetime: the chain length matters

Perovskite light emitting diodes suffer from poor operational stability, exhibiting a rapid decay of external quantum efficiency within minutes to hours after turn-on. To address this issue, we explore surface treatment of perovskite films with phenylalkylammonium iodide molecules of varying alkyl chain lengths. Combining experimental characterization and theoretical modelling, we show that these molecules stabilize the perovskite through suppression of iodide ion migration. The stabilization effect is enhanced with increasing chain length due to the stronger binding of the molecules with the perovskite surface, as well as the increased steric hindrance to reconfiguration for accommodating ion migration. The passivation also reduces the surface defects, resulting in a high radiance and delayed roll-off of external quantum efficiency. Using the optimized passivation molecule, phenylpropylammonium iodide, we achieve devices with an efficiency of 17.5%, a radiance of 1282.8 W sr-1 m-2 and a record T50 half-lifetime of 130 h under 100 mA cm-2.

Lignin Intermediates on Palladium: Insights into Keto-Enol Tautomerization from Theoretical Modelling

It has been suggested in the literature that keto-to-enol tautomerization plays a vital role for lignin fragmentation under mild conditions. On the other hand, previous modelling has shown that the adsorbed keto form is more stable than enol on the Pd(111) catalyst. The current density functional theory study of lignin model molecules shows that, in the gas-phase, keto is more stable than enol, but on the Pd surface, we find enol conformers that are at least as stable as keto. This supports the experimental result that the favourable reaction pathway for lignin depolymerization involves keto-enol tautomerization. An energy decomposition analysis gives insights concerning the origin of the fine energy balance between the keto and enol forms, where the molecule-surface interaction (-7 eV) and the molecular strain energy (+3 eV) are the main contributors to the adsorption energy.

Computational Study of Janus Transition Metal Dichalcogenide Monolayers for Acetone Gas Sensing

Recently, Janus two-dimensional (2D) transition metal dichalcogenides (TMDs) have been widely investigated and have provided exciting prospects in many fields such as photoelectric materials, photocatalysis, and gas sensors. In this study, we performed density functional theory (DFT) calculations to study the sensitivity of four volatile organic compounds (VOCs), including acetone, methanol, ethanol, and formyl aldehyde, over pristine 2D TMDs and 2D Janus TMD monolayers. We found that MoS2, Janus MoSSe, and Janus MoSTe demonstrated greater sensitivity toward acetone than other VOCs. Furthermore, the band gap values of the Janus MoSSe and Janus MoSTe monolayers dramatically changed after acetone adsorption on their sulfur layers, which was quite larger than the band gap change after acetone adsorption on the MoS2 monolayer. This result also leads to the extremely large conductivity change of Janus MoSSe and Janus MoSTe after sensing acetone. Hence, Janus MoSSe and Janus MoSTe monolayers show much higher sensitivity toward acetone in comparison with the pristine MoS2 monolayer. Finally, our finding indicates that Janus MoSSe and Janus MoSTe monolayers can be proposed as ultrahigh-sensitivity 2D TMD materials for acetone sensors.

Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations

Inspired by the analysis of Kohn-Sham energy densities by Nakai and coworkers, we extended the energy density analysis to linear-response time-dependent density functional theory (LR-TDDFT) calculations. Using ethylene-tetrafluoroethylene and oxyluciferin-water complexes as examples, distinctive distribution patterns were demonstrated for the excitation energy densities of local excitations (within a molecular fragment) and charge-transfer excitations (between molecular fragments). It also provided a simple way to compute the effective energy of both hot carriers (particle and hole) from charge-transfer excitations via an integration of the excitation energy density over the donor and acceptor grid points.