Abstract: Research progress on single layer group III monochalcogenides has been increasing rapidly owing to their interesting physics. Herein, we investigate the dynamically stable single layer forms of XBi (X = Ge, Si or Sn) using density functional theory calculations. Phonon band dispersion calculations and ab initio molecular dynamics simulations reveal the dynamical and thermal stability of the considered monolayers. Raman spectra calculations indicate the existence of 5 Raman active phonon modes, 3 of which are prominent and can be observed in possible Raman measurements. The electronic band structures of the XBi single layers were investigated with and without the effects of spin-orbit coupling (SOC). Our results show that XBi single layers show semiconducting properties with narrow band gap values without SOC. However, only single layer SiBi is an indirect band gap semiconductor, while GeBi and SnBi exhibit metallic behaviors when adding spin-orbit coupling effects. In addition, the calculated linear elastic parameters indicate the soft nature of the predicted monolayers. Moreover, our predictions for the thermoelectric properties of single layer XBi reveal that SiBi is a good thermoelectric material with increasing temperature. Overall, it is proposed that single layer XBi structures can be alternative, stable 2D single layers with varying electronic and thermoelectric properties.

Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: Inspired by the recent successful growth of Ti2C and Ti3C2 monolayers, here, we investigate the structural, electronic, and mechanical properties of functionalized Ti2C and Ti3C2 monolayers by means of density functional theory calculations. The results reveal that monolayers of Ti2C and Ti3C2 are dynamically stable metals. Phonon band dispersion calculations demonstrate that two-surface functionalization of Ti2C and Ti(3)C(2)via chalcogenides (S, Se, and Te), halides (F, Cl, Br, and I), and oxygen atoms results in dynamically stable novel functionalized monolayer materials. Electronic band dispersions and density of states calculations reveal that all functionalized monolayer structures preserve the metallic nature of both Ti2C and Ti3C2 except Ti2C-O-2, which possesses the behavior of an indirect semiconductor via full-surface oxygen passivation. In addition, it is shown that although halide passivated Ti3C2 structures are still metallic, there exist multiple Dirac-like cones around the Fermi energy level, which indicates that semi-metallic behavior can be obtained upon external effects by tuning the energy of the Dirac cones. In addition, the computed linear-elastic parameters prove that functionalization is a powerful tool in tuning the mechanical properties of stiff monolayers of bare Ti2C and Ti3C2. Our study discloses that the electronic and structural properties of Ti2C and Ti3C2 MXene monolayers are suitable for surface modification, which is highly desirable for material property engineering and device integration.

Abstract: Motivated by the recent successful synthesis of highly crystalline ultrathin BiTeCl and BiTeBr layered sheets [Debarati Hajra et al., ACS Nano, 2020, 14, 15626], herein for the first time, we carry out a comprehensive study on the structural and electronic properties of BiTeCl and BiTeBr Janus monolayers using density functional theory (DFT) calculations. Different structural and electronic parameters including the lattice constant, bond lengths, layer thickness in the z-direction, different interatomic angles, work function, charge density difference, cohesive energy and Rashba coefficients are determined to acquire a deep understanding of these monolayers. The calculations show good stability of the studied single layers. BiTeCl and BiTeBr monolayers are semiconductors with electronic bandgaps of 0.83 and 0.80 eV, respectively. The results also show that the semiconductor-metal transformation can be induced by increasing the number of layers. In addition, the engineering of the electronic structure is also studied by applying an electric field, and mechanical uniaxial and biaxial strain. The results show a significant change of the bandgaps and that an indirect-direct band-gap transition can be induced. This study highlights the positive prospect for the application of BiTeCl and BiTeBr layered sheets in novel electronic and energy conversion systems.

Abstract: In the present work, the identification of two hexagonal phases of germanium oxides (namely GeO2 and GeO) through the vibrational and optical properties is reported using density functional theory calculations. While structural optimizations show that single-layer GeO2 and GeO crystallize in 1T and buckled phases, phonon band dispersions reveal the dynamical stability of each structure. First-order off-resonant Raman spectral predictions demonstrate that each free-standing single-layer possesses characteristic peaks that are representative for the identification of the germanium oxide phase. On the other hand, electronic band dispersion analysis shows the insulating and large-gap semiconducting nature of single-layer GeO2 and GeO, respectively. Moreover, optical absorption, reflectance, and transmittance spectra obtained by means of G(0)W(0)-BSE calculations reveal the existence of tightly bound excitons in each phase, displaying strong optical absorption. Furthermore, the excitonic gaps are found to be at deep UV and visible portions of the spectrum, for GeO2 and GeO crystals, with energies of 6.24 and 3.10 eV, respectively. In addition, at the prominent excitonic resonances, single-layers display high reflectivity with a zero transmittance, which is another indication of the strong light-matter interaction inside the crystal medium.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.