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: 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 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: We present here a new interatomic potential parameter set to predict the thermal conductivity of zirconium trisulfide monolayers. The generated Tersoff-type force field is parameterized using data collected with first-principles calculations. We use non-equilibrium molecular dynamics simulations to predict the thermal conductivity. The generated parameters result in very good agreement in structural, mechanical, and dynamical parameters. The room temperature lattice thermal conductivity ( kappa) of the considered crystal is predicted to be kappa x x = 25.69Wm – 1K – 1 and kappa y y = 42.38Wm – 1K – 1, which both agree well with their corresponding first-principles values with a discrepancy of less than 5%. Moreover, the calculated kappa variation with temperature (200 and 400 K) are comparable within the framework of the accuracy of both first-principles and molecular dynamics simulations.

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: The structure and thermal boundary conductance of the wurtzite GaN/AlN (0001) interface are investigated using molecular dynamics simulation. Simulation results with three different empirical interatomic potentials have produced similar misfit dislocation networks and dislocation core structures. Specifically, the misfit dislocation network at the GaN/AlN interface is found to consist of pure edge dislocations with a Burgers vector of 1/3(1 (2) over bar 10) and the misfit dislocation core has an eight-atom ring structure. Although different interatomic potentials lead to different dislocation properties and thermal conductance values, all have demonstrated a significant effect of misfit dislocations on the thermal boundary conductance of the GaN/AlN (0001) interface. Published under an exclusive license by AIP Publishing.

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 report a high-throughput first-principles characterization of the structural, mechanical, electronic, and vibrational properties of tetragonal single-layer transition metal Xenes (t-TMXs). Our calculations revealed 22 dynamically, mechanically and chemically stable structures among the 96 possible free-standing layers present in the t-TMX family. As a fingerprint for their structural identification, we identified four characteristic Raman active phonon modes, namely three in-plane and one out-of-plane optical branches, with various intensities and frequencies depending on the material in question. Spin-polarized electronic calculations demonstrated that anti-ferromagnetic (AFM) metals, ferromagnetic (FM) metals, AFM semiconductors, and non-magnetic semiconductor materials exist within this family, evidencing the potential of t-TMXs for further use in multifunctional heterostructures.

Abstract: Two-dimensional materials have been attracting increasing interests because of their outstanding properties for Lithium-ion battery applications. In particular, a material family called MXenes (Mn+1Cn, where n = 1, 2, 3) have been recently attracted immense interest in this respect due to their incomparable fast-charging properties and high capacity promises. In this article, we review the state-of-the-art computational progress on Li-ion battery applications of MXene materials in accordance with our systematical DFT calculations. Structural, mechanical, dynamical, and electrical properties of 20 distinct MXene (M: Sc, Ti, V, Cr, Nb, Mo, Hf, Ta, W, and Zr) have been discussed. The battery performances of these MXene monolayers are further investigated by Li-ion binding energies, open circuit voltage values, and Li migration energy barriers. The experimental and theoretical progress up to date demonstrates particularly the potential of non-terminated or pristine MXene materials in Li ion-storage applications. Stability analyses show most of the pristine MXenes should be achievable, however susceptible to the development progress on the experimental growth procedures. Among pristine MXenes, Ti2C, V2C, Sc2C, and Zr2C compounds excel with their high charge/discharge rate prospect due to their extremely low Li diffusion energy barriers. Considering also their higher predicted gravimetric capacities, Sc, Ti, V, and Zr containing MXenes are more promising for their utilization in energy storage 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: The theoretical investigation on structural, vibrational, and electronic properties of zinc-blende (ZB) AgI were carried out employing first principles density functional theory calculations. Thermoelectric properties then were predicted through semi-classical Boltzmann transport equations within the constant relaxation time approximation. Equilibrium lattice parameter, bulk modulus, elastic constants, and vibrational properties were calculated by using generalized gradient approximation. Calculated properties are in good agreement with available experimental values. Electronic and thermoelectric properties were investigated both with and without considering spin-orbit coupling (SOC) effect which is found to have a strong influence on p-type Seebeck coefficient as well as the power factor of the ZB-AgI. By inclusion of SOC, a reduction of the band-gap and p-type Seebeck coefficients as well as the power factor was found which is the indication of that spin-orbit interaction cannot be ignored for p-type thermoelectric properties of the ZB-AgI. By using deformation potential theory for electronic relaxation time and experimentally predicted lattice thermal conductivity, we obtained aZTvalue 1.69 (0.89) at 400 K for n-type (p-type) carrier concentration of 1.5 x 10(18)(4.6 x10(19)) cm(-3)that makes ZB-AgI as a promising room temperature thermoelectric material.

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: Various theoretical and experimental methods are utilized to investigate the thermal conductivity of nanostructured materials; this is a critical parameter to increase performance of thermoelectric devices. Among these methods, equilibrium molecular dynamics (EMD) is an accurate technique to predict lattice thermal conductivity. In this study, by means of systematic EMD simulations, thermal conductivity of bulk Si-Ge structures (pristine, alloy and superlattice) and their nanostructured one dimensional forms with square and circular cross-section geometries (asymmetric and symmetric) are calculated for different crystallographic directions. A comprehensive temperature analysis is evaluated for selected structures as well. The results show that one-dimensional structures are superior candidates in terms of their low lattice thermal conductivity and thermal conductivity tunability by nanostructuring, such as by diameter modulation, interface roughness, periodicity and number of interfaces. We find that thermal conductivity decreases with smaller diameters or cross section areas. Furthermore, interface roughness decreases thermal conductivity with a profound impact. Moreover, we predicted that there is a specific periodicity that gives minimum thermal conductivity in symmetric superlattice structures. The decreasing thermal conductivity is due to the reducing phonon movement in the system due to the effect of the number of interfaces that determine regimes of ballistic and wave transport phenomena. In some nanostructures, such as nanowire superlattices, thermal conductivity of the Si/Ge system can be reduced to nearly twice that of an amorphous silicon thermal conductivity. Additionally, it is found that one crystal orientation, <100>, is better than the <111> crystal orientation in one-dimensional and bulk SiGe systems. Our results clearly point out the importance of lattice thermal conductivity engineering in bulk and nanostructures to produce high-performance thermoelectric materials.

Abstract: In recent years, transition metal dichalcogenides (TMDs) displaying astonishing properties are emerged as a new class of two-dimensional layered materials. The understanding and characterization of thermal transport in these materials are crucial for efficient engineering of 2D TMD materials for applications such as thermoelectric devices or overcoming general overheating issues. In this work, we obtain accurate Stillinger-Weber type empirical potential parameter sets for single-layer WS2 and WSe2 crystals by utilizing particle swarm optimization, a stochastic search algorithm. For both systems, our results are quite consistent with first-principles calculations in terms of bond distances, lattice parameters, elastic constants and vibrational properties. Using the generated potentials, we investigate the effect of temperature on phonon energies and phonon linewidth by employing spectral energy density analysis. We compare the calculated frequency shift with respect to temperature with corresponding experimental data, clearly demonstrating the accuracy of the generated inter-atomic potentials in this study. Also, we evaluate the lattice thermal conductivities of these materials by means of classical molecular dynamics simulations. The predicted thermal properties are in very good agreement with the ones calculated from first-principles. (C) 2017 Elsevier B.V. All rights reserved.