Abstract: Dibismuth dioxychalcogenides, Bi2O2Ch (Ch = S, Se, Te), are a promising class of materials for next-generation electronics and thermoelectrics due to their ultrahigh carrier mobility and excellent air stability. An interesting member of this family is Bi2O2S, which has a stereochemically active 6s2 lone pair of Bi3+ cations, heterogeneous bonding, and a high mass contrast between its constituent elements. In the present study, we have used first-principles calculations in combination with Boltzmann transport theory to systematically investigate the effect of hydrostatic pressure on lattice dynamics and phonon transport properties of Bi2O2S. We found that the ambient Pnmn phase has a low average lattice thermal conductivity (kappa l) of 1.71 W/(m K) at 300 K. We also predicted that Bi2O2S undergoes a structural phase transition from a low-symmetry (Pnmn) to a high-symmetry (I4/mmm) structure at around 4 GPa due to centering of Bi3+ cations with pressure. Upon compression, the lone pair activity of Bi3+ cations is suppressed, which increases kappa l by almost 3 times to 4.92 W/ (m K) at 5 GPa for the I4/mmm phase. The computed phonon lifetimes and Gru''neisen parameters show that anharmonicity decreases with increasing pressure due to further suppression of the lone pair activity and strengthening of intra-and intermolecular interactions, leading to an average room-temperature kappa l of 12.82 W/(m K) at 20 GPa. Overall, this study provides a comprehensive understanding of the effect of hydrostatic pressure on the stereochemical activity of the lone pair of Bi3+ cations and its implications on the phonon transport properties of Bi2O2S.

Abstract: Owing to its inherent non-trivial geometry, the unique structural motif of the recently discovered kagome topological superconductor AV(3)Sb(5) (A = K, Rb, Cs) is an ideal host of diverse topologically non-trivial phenomena, including giant anomalous Hall conductivity, topological charge order, charge density wave (CDW), and unconventional superconductivity. Despite possessing a normal-state CDW order in the form of topological chiral charge order and diverse superconducting gaps structures, it remains unclear how fundamental atomic-level properties and many-body effects including Fermi surface nesting, electron-phonon coupling, and orbital hybridization contribute to these symmetry-breaking phenomena. Here, the direct participation of the V3d-Sb5p orbital hybridization in mediating the CDW phase transition in CsV3Sb5 is reported. The combination of temperature-dependent X-ray absorption and first-principles studies clearly indicates the inverse Star-of-David structure as the preferred reconstruction in the low-temperature CDW phase. The results highlight the critical role that Sb orbitals play and establish orbital hybridization as the direct mediator of the CDW states and structural transition dynamics in kagome unconventional superconductors. This is a significant step toward the fundamental understanding and control of the emerging correlated phases from the kagome lattice through the orbital interactions and provides promising approaches to novel regimes in unconventional orders and topology.

Abstract: Experimental realizations of two-dimensional materials are hardly free of structural defects such as e.g. vacancies, which, in turn, modify drastically its pristine physical defect-free properties. In this work, we explore effects due to point defect clustering on the electronic and transport properties of bilayer graphene nanoribbons, for AA and AB stacking and zigzag and armchair boundaries, by means of the tight-binding approach and scattering matrix formalism. Evident vacancy concentration signatures exhibiting a maximum amplitude and an universality regardless of the system size, stacking and boundary types, in the density of states around the zero-energy level are observed. Our results are explained via the coalescence analysis of the strong sizeable vacancy clustering effect in the system and the breaking of the inversion symmetry at high vacancy densities, demonstrating a similar density of states for two equivalent degrees of concentration disorder, below and above the maximum value.

Abstract: 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 = c(1) + c(2)/m(2), 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.

Abstract: Using first-principles calculations, we demonstrated that high lithium storage capacity and fast kinetics are achieved for Ti3C2O2 by preintercalating organic molecules. As a proof-of-concept, two different quinone molecules, namely 1,4-benzoquinone (C6H4O2) and tetrafluoro-1,4-benzoquinone (C6F4O2) were selected as the molecular linkers to demonstrate the feasibility of this interlayer engineering strategy for energy storage. As compared to Ti3C2O2 bilayer without linker molecules, our pillared structures facilitate a much faster ion transport, promising a higher charge/discharge rate for Li. For example, while the diffusion barrier of a single Li ion within pristine Ti3C2O2 bilayer is at least 1.0 eV, it becomes 0.3 eV in pillared structures, which is comparable and even lower than that of commercial materials. At high Li concentrations, the calculated diffusion barriers are as low as 0.4 eV. Out-of-plane migration of Li ions is hindered due to large barrier energy with a value of around 1-1.35 eV. Concerning storage capacity, we can only intercalate one monolayer of Li within pristine Ti3C2O2 bilayer. In contrast, pillared structures offer significantly higher storage capacity. Our calculations showed that at least two layers of Li can be intercalated between Ti3C2O2 layers without forming bulk Li and losing the pillared structure upon Li loading/unloading. A small change in the in-plane lattice parameters (<0.5%) and volume (<1.0%) and ab initio molecular dynamics simulations prove the stability of the pillared structures against Li intercalation and thermal effects. Intercalated molecules avoid the large contraction/expansion of the whole structure, which is one of the key problems in electrochemical energy storage. Pillared structures allow us to realize electrodes with high capacity and fast kinetics. Our results open new research paths for improving the performance of not only MXenes but also other layered materials for supercapacitor and battery applications.

Abstract: Van der Waals quantum wells, naturally formed in two-dimensional layered materials with nanoscale thickness, possess many inherent advantages over conventional molecular beam epitaxy grown counterparts, and could bring up intriguing physics and applications. However, optical transitions originated from the series of quantized states in these emerging quantum wells are still elusive. Here, we show that multilayer black phosphorus appears to be an excellent candidate for van der Waals quantum wells with well-defined subbands and high optical quality. Using infrared absorption spectroscopy, we probe subband structures of multilayer black phosphorus with tens of atomic layers, revealing clear signatures for optical transitions with subband index as high as 10, far from what was attainable previously. Surprisingly, in addition to allowed transitions, an unexpected series of “forbidden” transitions is also evidently observed, which enables us to determine energy spacings separately for conduction and valence subbands. Furthermore, the linear tunability of subband spacings by temperature and strain is demonstrated. Our results are expected to facilitate potential applications for infrared optoelectronics based on tunable van der Waals quantum wells.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: 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).

Abstract: High-performance thermoelectric materials are critical in recuperating the thermal losses in various machinery and promising in renewable energy applications. In this respect, the search for novel thermoelectric materials has attracted considerable attention. In particular, low dimensional materials have been proposed as potential candidates due to their unique and controllable thermal and electronic transport properties. The considerable potential of several two-dimensional materials as thermoelectric devices has already been uncovered and many new candidates that merit further research have been suggested. In this regard, we comprehensively investigate the thermoelectric coefficients and electronic fitness function (EFF) of a large family of structurally isotropic and anisotropic two-dimensional layered materials using density functional theory combined with semi-classical Boltzmann transport theory. With this high-throughput screening, we bring to light additional 2D crystals that haven't been previously classified as favorable TE materials. We predict that Pb2Se2, GeS2, As-2, NiS2, Hf2O6, Zr2O6, AsBrS, ISbTe, ISbSe, AsISe, and AsITe are promising isotropic thermoelectric materials due to their considerably high EFF values. In addition to these materials, Hf2Br4, Zr2Br4, Hf2Cl4, Zr2Cl4, Hf2O6, Zr(2)O(6)and Os(2)O(4)exhibit strong anisotropy and possess prominently high EFF values.

Abstract: Semiconductor monolayer transition metal dichalcogenides (TMDs) have brought a new paradigm by introducing optically addressable valley degree of freedom. Concomitantly, their high flexibility constitutes a unique platform that links optics to mechanics via valleytronics. With the intention to expedite the research in this direction, we investigated ten TMDs, namely MoS2, MoSe2, MoTe2, WS2, WSe2, WTe2, MoSSe, MoSeTe, WSSe, and WSeTe, which particularly includes their so-called janus types (JTMDs). First, we obtained their electronic band structures using regular and hybrid density functional theory (DFT) calculations in the presence of the spin-orbit coupling and biaxial or uniaxial strain. Our DFT results indicated that against the expectations based on their reported piezoelectric behavior, JTMDs typically interpolated between the standard band properties of the constituent TMDs without producing a novel feature. Next, by fitting to our DFT data we generated both spinless and spinful k center dot p parameter sets which are quite accurate over the K valley where the optical activity occurs. As an important application of this parametrization, we considered the circular and linear dichroism under strain. Among the studied (J)TMDs, WTe2 stood out with its largest linear dichroism under uniaxial strain because of its narrower band gap and large K valley uniaxial deformation potential. This led us to suggest WTe2 monolayer membranes for optical polarization-based strain measurements, or conversely, as strain tunable optical polarizers.

Abstract: The newest members of a two-dimensional material family, involving transition metal carbides and nitrides (called MXenes), have garnered increasing attention due to their tunable electronic and thermal properties depending on the chemical composition and functionalization. This flexibility can be exploited to fabricate efficient electrochemical energy storage (batteries) and energy conversion (thermoelectric) devices. In this study, we calculated the Seebeck coefficients and lattice thermal conductivity values of oxygen terminated M2CO2 (where M = Ti, Zr, Hf, Sc) monolayer MXene crystals in two different functionalization configurations (model-II (MD-II) and model-III (MD-III)), using density functional theory and Boltzmann transport theory. We estimated the thermoelectric figure-of-merit, zT, of these materials by two different approaches, as well. First of all, we found that the structural model (i.e. adsorption site of oxygen atom on the surface of MXene) has a paramount impact on the electronic and thermoelectric properties of MXene crystals, which can be exploited to engineer the thermoelectric properties of these materials. The lattice thermal conductivity kappa(l), Seebeck coefficient and zT values may vary by 40% depending on the structural model. The MD-III configuration always has the larger band gap, Seebeck coefficient and zT, and smaller kappa(l) as compared to the MD-II structure due to a larger band gap, highly flat valence band and reduced crystal symmetry in the former. The MD-III configuration of Ti2CO2 and Zr2CO2 has the lowest kappa(l) as compared to the same configuration of Hf2CO2 and Sc2CO2. Among all the considered structures, the MD-II configuration of Hf2CO2 has the highest kappa(l), and Ti2CO2 and Zr2CO2 in the MD-III configuration have the lowest kappa(l). For instance, while the band gap of the MD-II configuration of Ti2CO2 is 0.26 eV, it becomes 0.69 eV in MD-III. The zT(max) value may reach up to 1.1 depending on the structural model of MXene.

Abstract: Molybdenum carbide has various applications for which studying the material using classical molecular dynamics simulations would be valuable. Here, we develop an empirical potential within the Tersoff formalism using particle swarm optimization for the orthorhombic phase of Mo2C. The developed potential is shown to predict lattice constants, elastic properties, and equation of state results that are consistent with current and previously reported results from experiments and first principles calculations. We demonstrate the potential with simulations of indentation using multiple indenter sizes that load and unload in three different directions relative to the crystallographic lattice of orthorhombic Mo2C. Direction-dependent force-displacement trends are analyzed and explained in terms of the spatial distributions of stress and strain within the material during indentation. This study reveals the anisotropic elasticity of orthorhombic Mo2C and, more generally, provides researchers with a new empirical potential that can be used to explore the properties and behavior of the material going forward.

Abstract: We detail the effects of Cl and S functionalization on the superconducting properties of layered (bulk) and monolayer niobium carbide (Nb2C) MXene crystals, based on first-principles calculations combined with Eliashberg theory. For bulk layered Nb2CCl2, the calculated superconducting transition temperature (T-c) is in very good agreement with the recently measured value of 6 K. We show that T-c is enhanced to 10 K for monolayer Nb2CCl2, due to an increase in the density of states at the Fermi level, and the corresponding electron-phonon coupling. We further demonstrate feasible gate- and strain-induced enhancements of T-c for both bulk-layered and monolayer Nb2CCl2 crystals, resulting in T-c values of around 38 K. In the S-functionalized Nb2CCl2 crystals, our calculations reveal the importance of phonon softening in understanding their superconducting properties. Finally, we predict that Nb3C2S2 in bulk-layered and monolayer forms is also superconducting, with a T-c of around 28 K. Considering that Nb2C is not superconducting in pristine form, our findings promote functionalization as a pathway towards robust superconductivity in MXenes.

Abstract: Recent experiments demonstrated that proton transport through graphene electrodes can be accelerated by over an order of magnitude with low intensity illumination. Here we show that this photo-effect can be suppressed for a tuneable fraction of the infra-red spectrum by applying a voltage bias. Using photocurrent measurements and Raman spectroscopy, we show that such fraction can be selected by tuning the Fermi energy of electrons in graphene with a bias, a phenomenon controlled by Pauli blocking of photo-excited electrons. These findings demonstrate a dependence between graphene's electronic and proton transport properties and provide fundamental insights into molecularly thin electrode-electrolyte interfaces and their interaction with light. Recent experiments have shown that proton transport through graphene electrodes can be promoted by light, but the understanding of this phenomenon remains unclear. Here, the authors report the electrical tunability of this photo-effect, showing a connection between graphene electronic and proton transport properties.

Abstract: Defect-free graphene is impermeable to all atoms(1-5) and ions(6,7) under ambient conditions. Experiments that can resolve gas flows of a few atoms per hour through micrometre-sized membranes found that monocrystalline graphene is completely impermeable to helium, the smallest atom(2,5). Such membranes were also shown to be impermeable to all ions, including the smallest one, lithium(6,7). By contrast, graphene was reported to be highly permeable to protons, nuclei of hydrogen atoms(8,9). There is no consensus, however, either on the mechanism behind the unexpectedly high proton permeability(10-14) or even on whether it requires defects in graphene's crystal lattice(6,8,15-17). Here, using high-resolution scanning electrochemical cell microscopy, we show that, although proton permeation through mechanically exfoliated monolayers of graphene and hexagonal boron nitride cannot be attributed to any structural defects, nanoscale non-flatness of two-dimensional membranes greatly facilitates proton transport. The spatial distribution of proton currents visualized by scanning electrochemical cell microscopy reveals marked inhomogeneities that are strongly correlated with nanoscale wrinkles and other features where strain is accumulated. Our results highlight nanoscale morphology as an important parameter enabling proton transport through two-dimensional crystals, mostly considered and modelled as flat, and indicate that strain and curvature can be used as additional degrees of freedom to control the proton permeability of two-dimensional materials. A study using high-resolution scanning electrochemical cell microscopy attributes proton permeation through defect-free graphene and hexagonal boron nitride to transport across areas of the structure that are under strain.