Gaussian approximation potentials for accurate thermal properties of two-dimensional materials.

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.

First-principles discovery of stable two-dimensional materials with high-level piezoelectric response.

The rational design of two-dimensional (2D) piezoelectric materials has recently garnered great interest due to their increasing use in technological applications, including sensor technology, actuating devices, energy harvesting, and medical applications. Several materials possessing high piezoelectric response have been reported so far, but a high-throughput first-principles approach to estimate the piezoelectric potential of layered materials has not been performed yet. In this study, we systematically investigated the piezoelectric (e 11, d 11) and elastic (C11 and C12) properties of 128 thermodynamically stable 2D semiconductor materials by employing first-principle methods. Our high-throughput approach demonstrates that the materials containing Group-V elements produce significantly high piezoelectric strain constants, d 11 > 40 pm V-1, and 49 of the materials considered have the e 11 coefficient higher than MoS2 insomuch as BrSSb has one of the largest d 11 with a value of 373.0 pm V-1. Moreover, we established a simple empirical model in order to estimate the d 11 coefficients by utilizing the relative ionic motion in the unit cell and the polarizability of the individual elements in the compounds.

Determination of Dynamically Stable Electrenes toward Ultrafast Charging Battery Applications.

Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first-principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).

A distinct correlation between the vibrational and thermal transport properties of group VA monolayer crystals.

The investigation of thermal transport properties of novel two-dimensional materials is crucially important in order to assess their potential to be used in future technological applications, such as thermoelectric power generation. In this respect, the lattice thermal transport properties of the monolayer structures of group VA elements (P, As, Sb, Bi, PAs, PSb, PBi, AsSb, AsBi, SbBi, P3As1, P3Sb1, P1As3, and As3Sb1) with a black phosphorus like puckered structure were systematically investigated by first-principles calculations and an iterative solution of the phonon Boltzmann transport equation. Phosphorene was found to have the highest lattice thermal conductivity, κ, due to its low average atomic mass and strong interatomic bonding character. As a matter of course, anisotropic κ was obtained for all the considered materials, owing to anisotropy in frequency values and phonon group velocities calculated for these structures. However, the determined linear correlation between the anisotropy in the κ values of P, As, and Sb is significant. The results corresponding to the studied compound structures clearly point out that thermal (electronic) conductivity of pristine monolayers might be suppressed (improved) by alloying them with the same group elements. For instance, the room temperature κ of PBi along the armchair direction was predicted to be as low as 1.5 W m-1 K-1, whereas that of P was predicted to be 21 W m-1 K-1. In spite of the apparent differences in structural and vibrational properties, we peculiarly revealed an intriguing correlation between the κ values of all the considered materials as κ = c1 + c2/m2, in particular along the zigzag direction. Furthermore, our calculations on compound structures clearly showed that the thermoelectric potential of these materials can be improved by suppressing their thermal properties. The presence of ultra-low κ values and high electrical conductivity (especially along the armchair direction) makes this class of monolayers promising candidates for thermoelectric applications.