Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

Abstract: X-ray diffraction and energy dispersive x-ray spectroscopic analyses showed a natural galena (PbS) crystal from Freiberg in Saxony (Germany) to be a single phase specimen [rock salt (NaCl) structure type, space group Fm3m, a = 5.932(1) angstrom] with stoichiometric composition and an enhanced dislocation density (8 approximate to 1011 cm-2). The latter parameter leads to an increase of the electrical resistivity in the high-temperature regime, as well as to the appearance of phonon resonance with a characteristic frequency coPR = 3.8(1) THz. Being in the same range (i.e., 3-5.5 THz) with the sulfur optical modes of highest group velocities, it results in a drastic reduction (by similar to 75%) of thermal conductivity (K) at lower temperatures (i.e., < 100 K), as well as in the appearance of a characteristic minimum in K at T approximate to 30 K. Furthermore, the studied galena is characterized by phonon-drag behavior and by temperature dependent switch of the charge carrier scattering mechanism regime (i.e., scattering on dislocations for T < 100 K, on acoustic phonons for 100 K < T < 170 K and on both acoustic and optical phonons for 170 K < T < 300 K). The combined theoretical calculation and optical spectroscopic study confirm this mineral to be a direct gap degenerate semiconductor. The possible origins of the second-order Raman spectrum are discussed.

Abstract: This study focuses on optimizing the composition and firing temperature of porcelain tiles using statistical analysis techniques. A full factorial design, including model adequacy checking, analysis of variance, Pareto charts, interaction plots, regression model, and response optimizer is employed. The key factors were the Seger ratios of SiO2/Al2O3, Na2O/K2O, MgO/CaO, and firing temperature. The response variables investigated were bulk density, water absorption, linear shrinkage, coefficient of thermal expansion (at 500 degrees C), and strength. The statistical analysis revealed highly significant results, which were further validated, confirming their reliability for practical use in the production of porcelain tiles. The study demonstrated the effectiveness of utilizing Seger formulas and properties of typical raw materials to accurately predict the final properties of ceramic tiles. By employing SiO2/Al2O3 = 5.2, Na2O/K2O = 1.50, MgO/CaO = 3.0, and firing temperature of 1180 degrees C, optimized properties, such as maximum strength, maximum bulk density, and minimum water absorption, was achieved with a composite desirability of 0.9821.

Abstract: We demonstrate the possibility of engineering the optical properties of transition metal dichalcogenide heterobilayers when one of the constitutive layers has a Janus structure. We investigate different MoS2@Janus layer combinations using first-principles methods including excitons and exciton-phonon coupling. The direction of the intrinsic electric field from the Janus layer modifies the electronic band alignments and, consequently, the energy separation between dark interlayer exciton states and bright in-plane excitons. We find that in-plane lattice vibrations strongly couple the two states, so that exciton-phonon scattering may be a viable generation mechanism for interlayer excitons upon light absorption. In particular, in the case of MoS2@WSSe, the energy separation of the low-lying interlayer exciton from the in-plane exciton is resonant with the transverse optical phonon modes (40 meV). We thus identify this heterobilayer as a prime candidate for efficient generation of charge-separated electron-hole pairs.

Abstract: Two-dimensional transition metal carbides and nitrides (MXenes) are an emerging class of atomically-thin superconductors, whose characteristics are highly prone to tailoring by surface functionalization. Here we explore the use of hydrogen adatoms to enhance phonon-mediated superconductivity in MXenes, based on first-principles calculations combined with Eliashberg theory. We first demonstrate the stability of three different structural models of hydrogenated Mo- and W-based MXenes. Particularly high critical temperatures of over 30 K are obtained for hydrogenated Mo2N and W2N. Several mechanisms responsible for the enhanced electron-phonon coupling are uncovered, namely (i) hydrogen-induced changes in the phonon spectrum of the host MXene, (ii) emerging hydrogen-based phonon modes, and (iii) charge transfer from hydrogen to the MXene layer, boosting the density of states at the Fermi level. Finally, we demonstrate that hydrogen adatoms are moreover able to induce superconductivity in MXenes that are not superconducting in pristine form, such as Nb2C.

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: Metallic contacts are critical components of electronic devices and the presence of a large Schottky barrier is detrimental for an optimal device operation. Here, we show by using first-principles calculations that a self-assembled monolayer (SAM) of polar molecules between the metal electrode and MoSe2 monolayer is able to convert the Schottky contact into an almost Ohmic contact. We choose -CH3 and -CF3 terminated short-chain alkylthiolate (i.e. SCH3 and fluorinated alkylthiolates (SCF3)) based SAMs to test our approach. We consider both high (Au) and low (Sc) work function metals in order to thoroughly elucidate the role of the metal work function. In the case of Sc, the Fermi level even moves into the conduction band of the MoSe2 monolayer upon SAM insertion between the metal surface and the MoSe2 monolayer, and hence possibly switches the contact type from Schottky to Ohmic. The usual Fermi level pinning at the metal-transition metal dichalcogenide (TMD) contact is shown to be completely removed upon the deposition of a SAM. Systematic analysis indicates that the work function of the metal surface and the energy level alignment between the metal electrode and the TMD monolayer can be tuned significantly by using SAMs as a buffer layer. These results clearly indicate the vast potential of the proposed interface engineering to modify the physical and chemical properties of MoSe2.

Abstract: Group II-VI semiconductor honeycomb monolayers have a noncentrosymmetric crystal structure and therefore are expected to be important for nano piezoelectric device applications. This motivated us to perform first principles calculations based on density functional theory to unveil the piezoelectric properties (i.e., piezoelectric stress (e(11)) and piezoelectric strain (d(11)) coefficients) of these monolayer materials with chemical formula MX (where M = Be, Mg, Ca, Sr, Ba, Zr, Cd and X = S, Se, Te). We found that these two-dimensional materials have peculiar piezoelectric properties with d(11) coefficients 1 order of magnitude larger than those of commercially utilized bulk materials. A clear trend in their piezoelectric properties emerges, which

Abstract: By performing density functional theory based first-principles calculations, the electronic, vibrational, elastic, and piezoelectric properties of two dynamically stable crystal phases of monolayer Janus MoSTe, namely 1H-MoSTe and 1T'-MoSTe, are investigated. Vibrational frequency analysis reveals that the other possible crystal structure, 1T-MoSTe, of this Janus monolayer does not exhibit dynamical stability. The 1H-MoSTe phase is found to be an indirect band-gap semiconductor while 1T'-MoSTe is predicted as small-gap semiconductor. Notably, in contrast to the direct band-gap nature of monolayers 1H-MoS2 and 1H-MoTe2, 1H-MoSTe is found to be an indirect gap semiconductor driven by the induced surface strains on each side of the structure. The calculated Raman spectrum of each structure shows unique character enabling us to clearly distinguish the stable crystal phases via Raman measurements. The systematic piezoelectric stress and strain coefficient analysis reveals that out-of-plane piezoelectricity appears in 1H-MoSTe and the noncentral symmetric 1T'-MoSTe has large piezoelectric coefficients. Static total-energy calculations show clearly that the formation of 1T'-MoSTe is feasible by using 1T'-MoTe2 as a basis monolayer. Therefore, we propose that the Janus MoSTe structure can be fabricated in two dynamically stable phases which possess unique electronic, dynamical, and piezoelectric properties.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: MXenes are the newest class of two-dimensional (2D) materials, and they offer great potential in a wide range of applications including electronic devices, sensors, and thermoelectric and energy storage materials. In this work, we combined the outstanding electrical conductivity, that is essential for battery applications, of graphene with MXene monolayers (M2CX2 where M = Sc, Ti, V and X = OH, O) to explore its potential in Li battery applications. Through first principles calculations, we determined the stable stacking configurations of M2CX2/graphene bilayer heterostructures and their Li atom intercalation by calculating the Li binding energy, diffusion barrier and voltage. We found that: (1) for the ground state stacking, the interlayer binding is strong, yet the interlayer friction is small; (2) Li binds more strongly to the O-terminated monolayer, bilayer and heterostructure MXene systems when compared with the OHterminated MXenes due to the H+ induced repulsion to the Li atoms. The binding energy of Li decreases as the Li concentration increases due to enhanced repulsive interaction between the positively charged Li ions; (3) Ti2CO2/graphene and V2CO2/graphene heterostructures exhibit large Li atom binding energies making them the most promising candidates for battery applications. When fully loaded with Li atoms, the binding energy is -1.43 eV per Li atom and -1.78 eV per Li atom for Ti2CO2/graphene and V2CO2/graphene, respectively. These two heterostructures exhibit a nice compromise between storage capacity and kinetics. For example, the diffusion barrier of Li in Ti2CO2/graphene is around 0.3 eV which is comparable to that of graphite. Additionally, the calculated average voltages are 1.49 V and 1.93 V for Ti2CO2/graphene and V2CO2/graphene structures, respectively; (4) a small change in the in-plane lattice parameters (<1%), interatomic bond lengths and interlayer distances (<0.5 angstrom) proves the stability of the heterostructures against Li intercalation, and the impending phase separation into constituent layers and capacity fading during charge-discharge cycles in real battery applications; (5) as compared to bare M2CX2 bilayers, M2CX2/graphene heterostructures have lower molecular mass, offering high storage capacity; (6) the presence of graphene ensures good electrical conductivity that is essential for battery applications. Given these advantages, Ti2CO2/graphene and V2CO2/graphene heterostructures are predicted to be promising for lithium-ion battery applications.

Abstract: The adsorption and diffusion of Li, Na, K and Ca atoms on a Mo2C monolayer are systematically investigated by using first principles methods. We found that the considered metal atoms are strongly bound to the Mo2C monolayer. However, the adsorption energies of these alkali and earth alkali elements decrease as the coverage increases due to the enhanced repulsion between the metal ions. We predict a significant charge transfer from the ad-atoms to the Mo2C monolayer, which indicates clearly the cationic state of the metal atoms. The metallic character of both pristine and doped Mo2C ensures a good electronic conduction that is essential for an optimal anode material. Low migration energy barriers are predicted as small as 43 meV for Li, 19 meV for Na and 15 meV for K, which result in the very fast diffusion of these atoms on Mo2C. For Mo2C, we found a storage capacity larger than 400 mA h g(-1) by the inclusion of multilayer adsorption. Mo2C expands slightly upon deposition of Li and Na even at high concentrations, which ensures the good cyclic stability of the atomic layer. The calculated average voltage of 0.68 V for Li and 0.30 V for Na ions makes Mo2C attractive for low charging voltage applications.