In this study, non-equilibrium molecular dynamics (NEMD) simulations with a reactive force field were used to investigate the tribochemical properties of glycerol, with and without water, confined between two ferrous surfaces. The results demonstrated that glycerol significantly reduced friction on α-Fe slabs more effectively than on functionalized amorphous magnetite. A numerical method was introduced to identify the interface region and evaluate the dissociated surface atoms. It was found that the dissociation rate of glycerol molecules increased with applied normal pressure and shear stress. Additionally, the production rate of water molecules from glycerol dissociation was consistently positive for all solutions above 80 % wt. The assumption of linear velocity distribution across the film thickness was validated for all systems studied.

In this work, we investigate the ion pair tetramethylphosphonium cation, [P1,1,1,1]+, and bis(oxalato)borate anion, [BOB]−, as a model system for the study of ionic liquids interacting with both hydroxylated and oxygen terminated α-SiO2 (001) surfaces, using first-principles electronic structure theory. We use a single ionic pair and clusters of ion pairs, in order to have exclusively neutral supercell slab models. We use dispersion-corrected density functional theory (DFT) to ascertain that both the strong physical binding between the ions, dominated by ionic binding, and the weaker physical binding of ions to the different surfaces are correctly described. We have found that the binding of ion pairs is stronger to the hydroxylated α-SiO2 (001) surface compared to the oxygen terminated surface, which is attributed to the formation of H-binding with the oxygen atom(s) of the [BOB]− anion. Through rotation of ionic pair(s), we estimate the surface-ions energy barrier for translational movement and, thus, the strength of H-binding of the ions. At the surface of hydroxylated α-SiO2 (001), we have studied how water molecules form a network of H-binding with the OH groups of the surface and the [BOB]− anion, which offers an explanation for the reduction in the friction of ionic liquids on the inclusion of water. We suggest modeling protocols for simulation of ion pairs on surfaces, which can open up the possibility to use DFT to aid in designing and understanding the physicochemical mechanism of interactions of ionic materials (including ionic liquids) in various technological applications.

We present first-principles results on monolayer (ML) and bulk titanium carbide (TiC) exhibiting dual topological characteristics: a quantum spin Hall insulator state in two-dimensional ML and a type-II Dirac semimetal state in its three-dimensional bulk form. The nontrivial nature of ML TiC is confirmed through the calculation of a Z2 invariant, spin Hall conductivity, and edge states. The edge band structure of ML TiC displays a single pair of gapless edge states at the 𝑀 point. The insulating topological phase in ML TiC is driven by a band inversion around the 𝑀 point involving Ti-𝑑 orbitals, with a nontrivial band gap of 0.47 eV. Our findings also indicate that ML TiC possesses excellent dynamical and thermal stability. Moreover, the bulk, ML hollow sphere arrays, and a 4-nm-thick film of TiC have already been synthesized, suggesting the high feasibility of the experimental synthesizing of ML TiC. On the other hand, we demonstrate that the bulk TiC hosts both a type-II Dirac cone and a topological nodal surface without spin-orbit coupling, protected by a combined symmetry of inversion and time reversal. The presence of spin-orbit coupling removes the nodal surface while preserving the type-II Dirac cone. Surface states connecting bulk Dirac nodes in bulk TiC can be readily observed, making the surface characteristics easily detectable in experiments.

Carbon nanotubes (CNTs), hollow cylinders of carbon, hold great promise for advanced technologies, provided their structure remains uniform throughout their length. Their growth takes place at high temperatures across a tube-catalyst interface. Structural defects formed during growth alter CNT properties. These defects are believed to form and heal at the tube-catalyst interface but an understanding of these mechanisms at the atomic-level is lacking. Here we present DeepCNT-22, a machine learning force field (MLFF) to drive molecular dynamics simulations through which we unveil the mechanisms of CNT formation, from nucleation to growth including defect formation and healing. We find the tube-catalyst interface to be highly dynamic, with large fluctuations in the chiral structure of the CNT-edge. This does not support continuous spiral growth as a general mechanism, instead, at these growth conditions, the growing tube edge exhibits significant configurational entropy. We demonstrate that defects form stochastically at the tube-catalyst interface, but under low growth rates and high temperatures, these heal before becoming incorporated in the tube wall, allowing CNTs to grow defect-free to seemingly unlimited lengths. These insights, not readily available through experiments, demonstrate the remarkable power of MLFF-driven simulations and fill long-standing gaps in our understanding of CNT growth mechanisms.

We present first-principles insights into the electrical and electrochemical properties of Cu2N, a newly synthesized two-dimensional material that features a planar, checkerboard lattice structure [Hu et al., Nano Lett. 2023, 23 (12), 5610–5616]. We evaluate the suitability of monolayer Cu2N as an anode material for Li and Na-ion batteries by examining its storage capacity, diffusion barrier, open-circuit voltage (OCV), volume expansion, and the impact of defects on its electrochemical performance. The monolayer Cu2N demonstrates a storage capacity of 379.88 mAh.g−1 for both Li and Na, comparable to that of commercial graphite for Li (372 mAh.g−1) and significantly higher for Na (less than 35 mAh.g−1). The migration barriers for Li and Na are found to be 0.1 eV and 0.01 eV, respectively, substantially lower than those theoretically reported for commercial anodes TiO2 (0.4–1.0 eV) and graphite (∼0.4 eV), which imply that monolayer Cu2N demonstrates excellent charge/discharge capabilities. Moreover, the volume growth of monolayer Cu2N is 4.14 % with maximal Li adsorption, which is 2.4 times less than graphite. The analysis of vacancy defects reveals a significant enhancement in the binding energies of Li and Na atoms, accompanied by minimal changes in diffusion barriers. Since monolayer Cu2N has already been successfully synthesized, these findings would pave the way for large-scale experimental fabrication of monolayer Cu2N as a battery anode.

Recently, triboelectric nanogenerators (TENGs) have been widely used for energy harvesting and self-powered sensing due to their significant and unique advantages. However, the intrinsic mechanisms that contribute to tribo-electricification (TE) between two materials remain as a subject of rigorous debate. In addition to predicting the qualitative charge transfer in solid-solid contacts based on the difference in the work functions of the two moieties constituting the interface, we argue that it is essential to obtain atomic-level, first principles, insights into the bonding properties, quantitative charge transfer, and the possible presence of a electrostatic potential barrier at the interface to fully understand the TE mechanism of a system. We have utilized dispersion-corrected density functional theory (DFT) calculations in this study to systematically investigate the TE potential of bare surface Ti3C2 and Ti3N2 2D MXene monolayers and their surface functionalized modifications Ti3C2R2 and Ti3N2R2 (where R = -O, -OH, or -F) in contact with Al(111). For these heterostructures, we have analyzed the adhesive energy of the interfaces, the nature of interaction through the electron localization function (ELF), and the charge distribution, which have revealed distinct characteristics of MXene/Al contacts for these monolayer/metal interfaces at their equilibrium distance and the changes in their properties under uniaxial pressure. Among all the metallic 2D MXene variants investigated in this study, we have determined that Ti3C2F2/Al and Ti3N2F2/Al interfaces show exceptional potential for TE.

Two-dimensional Dirac materials have stimulated substantial research interest as binder-free anodes in metal-ion batteries, owing to their ultrahigh electronic conductivity, large specific area, and higher energy density. Here, using first-principles density functional theory calculations, we have investigated the feasibility of monolayer TiC as a potential anode material for Li/Na-ion batteries. The results indicate that monolayer TiC exhibits excellent dynamical and thermal stability. The electronic structure of monolayer TiC shows semimetallic characteristics with a Dirac cone at the M high symmetry point and the formation of Ti or C vacancies transforms the Dirac cone into a nodal loop or a nodal surface, respectively. Thus, monolayer TiC possesses superior electrical conductivity, which can be further enhanced by the formation of Ti or C vacancies in the material. Furthermore, the calculated adsorption energy values of -0.85 and -0.46 eV for Li-ion and Na-ion, respectively, indicate that Li/Na atom adsorption over monolayer TiC is a favorable process. The density of states plots show that after the adsorption of a single Li/Na atom, monolayer TiC maintains its metallic state, which is advantageous for the diffusion of stored electrons. Most remarkably, monolayer TiC exhibits energy densities of 2684 and 2015 mWh/g for Li and Na, respectively, which are significantly higher than commercial graphite and most other 2D anode materials. The fully loaded TiC anode exhibits excellent cycle stability with volume expansions as low as 0.13 and 0.11%, for Li and Na, respectively. Furthermore, an ultrafast diffusivity with low energy barriers of 0.02 and 0.10 eV is found in monolayer TiC for Li-ion and Na-ion, respectively, which suggests that it has an excellent charge/discharge capability. These exceptional properties make monolayer TiC an excellent candidate as an anode material for Li-ion and Na-ion batteries. Finally, SiC(111) has been proposed as a candidate substrate for monolayer TiC due to its minimal lattice mismatch.

Crystalline thermoelectric materials, especially SnSe crystals, have emerged as promising candidates for power generation and electronic cooling. In this study, significant enhancement in ZT is achieved through the combined effects of lattice distortions and band convergence in multiple electronic valence bands. Density functional theory (DFT) calculations demonstrate that cation vacancies together with Pb substitutional doping promote the band convergence and increase the density of states (DOS) near the Fermi surface of SnSe, leading to a notable increase in the Seebeck coefficient (S). The complex defects formed by Sn vacancies and Pb doping not only boost the Seebeck coefficient but also optimize carrier concentration (nH) and enhance electrical conductivity (σ), resulting in a higher power factor (PF). Furthermore, the localized lattice distortions induced by these defects increase phonon scattering, significantly reducing lattice thermal conductivity (κlat) to as low as 0.29 W m−1 K−1at 773 K in Sn0.92Pb0.03Se. Consequently, these synergistic effects on phonon and electron transport contribute to a high ZT of 1.8. This study provides a framework for rational design of high-performance thermoelectric materials based on first-principles insights and experimental validation.

The adsorption of organic molecules to surfaces is a central issue to achieve fully-functional molecular devices, for which porphyrins are well-studied due to their chemical stability and functional diversity. Herein, we investigate both the physical and the chemical adsorption of the free-base tetraphenylporphyrin 2H-TPP on the Cu(111) surface within the framework of density functional theory and find that the most stable physisorbed configuration is more weakly bound by -0.36 eV than the chemisorbed configuration. We use the electron localization function to investigate the difference in binding mechanisms between strong physisorption and weak chemisorption. We have computed a reaction barrier of 0.12 eV in going from physical binding to chemical bonding to the surface, and a barrier of 50 meV in going between neighboring physical binding sites. Our results support the possibility of realizing free-base porphyrins either physisorbed or chemisorbed on Cu(111) depending on the deposition procedure and experimental conditions.