Abstract: Recently, an extremely low lattice thermal conductivity value has been reported for the alkali-based telluride material BaIn2Te4. The value is comparable with low-thermal conductivity metal chalcogenides, and the glass limit is highly intriguing. Therefore, to shed light on this issue, we performed first-principles phonon thermal transport calculations. We predicted highly anisotropic lattice thermal conductivity along different directions via the solution of the linearized phonon Boltzmann transport equation. More importantly, we determined several different factors as the main sources of the predicted ultralow lattice thermal conductivity of this crystal, such as the strong interactions between low-frequency optical phonons and acoustic phonons, small phonon group velocities, and lattice anharmonicity indicated by large negative mode Gruneisen parameters. Along with thermal transport calculations, we also investigated the electronic transport properties by accurately calculating the scattering mechanisms, namely the acoustic deformation potential, ionized impurity, and polar optical scatterings. The inclusion of spin-orbit coupling (SOC) for electronic structure is found to strongly affect the p-type Seebeck coefficients. Finally, we calculated the thermoelectric properties accurately, and the optimal ZT value of p-type doping, which originated from high Seebeck coefficients, was predicted to exceed unity after 700 K and have a direction averaged value of 1.63 (1.76 in the y-direction) at 1000 K around 2 x 1020 cm-3 hole concentration. For n-type doping, a ZT around 3.2 x 1019 cm-3 concentration was predicted to be a direction-averaged value of 1.40 (1.76 in the z-direction) at 1000 K, mostly originating from its high electron mobility. With the experimental evidence of high thermal stability, we showed that the BaIn2Te4 compound has the potential to be a promising mid- to high-temperature thermoelectric material for both p-type and n-type systems with appropriate doping.

Abstract: Different allotropes of phosphorene are possible of which black and blue phosphorus are the most stable. While blue phosphorus has isotropic properties, black phosphorus is strongly anisotropic in its electronic and optical properties due to its anisotropic crystal structure. In this work, we systematically investigated the lattice thermal properties of black and blue phosphorene by using first-principles calculations based on the quasiharmonic approximation approach. Similar to the optoelectronic and electronic properties, we predict that black phosphorene has highly anisotropic thermal properties, in contrast to the blue phase. The linear thermal expansion coefficients along the zigzag and armchair direction differ up to 20% in black phosphorene. The armchair direction of black phosphorene is more expandable as compared to the zigzag direction and the biaxial expansion of blue phosphorene under finite temperature. Our comparative analysis reveals that the inclusion of finite-temperature effects makes the blue phase thermodynamically more stable over the black phase above 135 K.

Abstract: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

Abstract: Thermal expansion behavior of two-dimensional (2D) nitrides and graphene were studied by ab initio molecular dynamics (MD) simulations as well as quasiharmonic approximation (QHA). Anharmonicity of the acoustic phonon modes are related to the unusual negative thermal expansion (NTE) behavior of the nitrides. Our results also hint that direct ab initio MD simulations are a more elaborate method to investigate thermal expansion behavior of 2D materials than the QHA. Nevertheless, giant NTE coefficients are found for h-GaN and h-AlN within the covered temperature range 100-600 K regardless of the chosen computational method. This unusual NTE of 2D nitrides is reasoned with the out-of-plane oscillations related to the rippling behavior of the monolayers.

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: 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 atomically laminated layered ternary transition-metal borides (the MAB phases) have demonstrated outstanding properties and have been applied in various fields. Understanding their thermal and mechanical properties is critical to determining their applicability in various fields such as high-temperature applications. To achieve this, we conducted first-principles calculations based on density-functional theory and the quasiharmonic approximation to determine the thermal expansion coefficients, Gruneisen parameters, bulk moduli, hardness, thermal conductivity, electron-phonon coupling parameters, and the structural and vibrational properties of MoAlB, WAlB, Cr2AlB2, and Ti2InB2. We found varying degrees of anisotropy in the thermal expansion and mechanical properties in spite of similarities in their crystal structures. MoAlB has a mild degree of anisotropy in its thermal expansion coefficient (TEC), while Cr2AlB2 and WAlB display the highest level of TEC anisotropy. We assessed various empirical models to calculate hardness and thermal conductivity, and correlated the calculated values with the material properties such as elastic moduli, Gruneisen parameter, Debye temperature, and type of bonding. Owing to their higher Gruneisen parameters, implying a greater degree of anharmonicity in lattice vibrations and lower phonon group velocities, MoAlB and WAlB have significantly lower lattice thermal conductivity values than those of Cr2AlB2 and Ti2InB2. The hardness and lattice thermal conductivity of MAB phases can be predicted with high accuracy if one utilizes an appropriate model.

Abstract: Seminal Goodenough-Kanamori-Anderson (GKA) rules provide an inceptive understanding of the superexchange interaction of two magnetic metal ions bridged with an anion, and suggest fostered ferromagnetic interaction for orthogonal bridging bonds. However, there are no examples of two-dimensional (2D) materials with structure that optimizes the GKA arguments towards enhanced ferromagnetism and its critical temperature. Here we reveal that an ideally planar GKA ferromagnetism is indeed stable in selected tetragonal transition-metal xenes (tTMXs), with Curie temperature above 300 K found in CrC and MnC. We provide the general orbitally-resolved analysis of magnetic interactions that supports the claims and sheds light at the mechanisms dominating the magnetic exchange process in these structures. Furthermore, we propose the set of three GKA-like rules that will guarantee room temperature ferromagetnism. With recent advent of epitaxially-grown tetragonal 2D materials, our findings earmark tTMXs for facilitated spintronic and magnonic applications, or as a desirable magnetic constituent of functional 2D heterostructures.

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: 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: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

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