Abstract: In this study, we predicted two new stable metallic Re-C based monolayer structures with a rectangular (r-ReC2) and a hexagonal (h-Re2C) crystal symmetry using first-principle calculations based on density functional theory. Our results obtained from mechanical and phonon calculations and high-temperature molecular dynamic simulations clearly proved the stability of these two-dimensional (2D) crystals. Interestingly, Re-C monolayers in common transition metal carbide structures (i.e. MXenes) were found to be unstable, contrary to expectations. We found that the stable structures, i.e. r-ReC2 and h-Re2C, display superior mechanical properties over the well-known 2D materials. The Young's modulus for r-ReC2 and h-Re2C are extremely high and were calculated as 351 (1310) and 617 (804) N/m (GPa), respectively. Both materials have larger Young's modulus values than the most of the well-known 2D materials. We showed that the combination of the short strong directional p-d bonds, the high coordination number of atoms in the unit-cell and high valence electron density result in strong mechanical properties. Due to its crystal structure, the r-ReC2 monolayer has anisotropic mechanical properties and the crystallographic direction parallel to the C-2 dimers is stiffer compared to perpendicular direction due to strong covalent bonding within C-2 dimers. h-Re2C was derived from the corresponding bulk structure for which we determined the critical thickness for the dynamically stable bulk-derived monolayer structures. In addition, we also investigated the electronic of these two stable structures. Both exhibit metallic behavior and Re-5d orbitals dominate the states around the Fermi level. Due to their ultra high mechanical stability and stiffness, these novel Re-C monolayers can be exploited in various engineering applications.

Abstract: The peculiar potential transition metal dichalcogenides in regard to sensor and device applications have been exhibited by both experimental and theoretical studies. The use of these materials, thermodynamically stable even at elevated temperatures, particularly in nano- and optoelectronic technology, is about to come true. On the other hand, the distinct electronic and thermal transport properties possessing unique coherency, which may result in higher thermoelectric efficiency, have also been reported. However, exploiting this potential in terms of power generation and cooling applications requires a deeper understanding of these materials in this regard. This perspective study, concentrated with this intention, summarizes thermoelectric research based on transition metal dichalcogenides from a broad perspective and also provides a general evaluation of future theoretical investigations inevitable to shed more light on the physics of electronic and thermal transport in these materials and to lead future experimental research.

Abstract: Efficient thermal transport control is a fundamental issue for electronic device applications such as information, communication, and energy storage technologies in modern electronics in order to achieve desired thermal conditions. Structural defects in materials provide a mechanism to adjust the thermal transport properties of these materials on demand. In this context, the effect of structural defects on lattice thermal conductivities of two-dimensional hexagonal binary group-III nitride (XN, X = B, Al, and Ga) semiconductors is systematically investigated by means of classical molecular dynamics simulations performed with recently developed transferable inter-atomic potentials accurately describing defect energies. Here, two different Green-Kubo based approaches and another approach based on non-equilibrium molecular dynamics are compared in order to get an overall understanding. Our investigation clearly shows that defect concentrations of 3% decrease the thermal conductivity of systems containing these nitrites up to 95%. Results hint that structural defects can be used as effective adjustment parameters in controlling thermal transport properties in device applications associated with these materials. Published under an exclusive license by AIP Publishing.

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 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: Piezoelectricity is a unique material property that allows one to convert mechanical energy into electrical one or vice versa. Transition metal dichalcogenides (TMDC) and transition metal dioxides (TMDO) are expected to have great potential for piezoelectric device applications due to their noncentrosymmetric and two-dimensional crystal structure. A detailed theoretical investigation of the piezoelectric stress (e 11 ) and piezoelectric strain (d 11 ) coefficients of single layer TMDCs and TMDOs with chemical formula MX 2 (where M= Cr, Mo, W, Ti, Zr, Hf, Sn and X = O, S, Se, Te) is presented by using first-principles calculations based on density func- tional theory. We predict that not only the Mo- and W-based members of this family but also the other materials with M= Cr, Ti, Zr and Sn exhibit highly promising piezoelectric properties. CrTe 2 has the largest e 11 and d 11 coefficients among the group VI elements (i.e., Cr, Mo, and W). In addition, the relaxed-ion e 11 and d 11 coefficients of SnS 2 are almost the same as those of CrTe 2 . Furthermore, TiO 2 and ZrO 2 pose comparable or even larger e 11 coefficients as compared to Mo- and W-based TMDCs and TMDOs. Our calculations reveal that TMDC and TMDO structures are strong candidates for future atomically thin piezoelectric applications such as transducers, sensors, and energy harvesting devices due to their piezoelectric coefficients that are comparable (even larger) to currently used bulk piezoelectric materials.

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: 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: Combining first-principles calculations with Landauer-Mittiker formalism, ballistic thermoelectric transport properties of semiconducting two-dimensional transition metal dichalcogenides (TMDs) and oxides (TMOs) (namely MX2 with M = Cr, Mo, W, Ti, Zr, Hf; X = O, S, Se, Te) are investigated in their 2H and 1T phases. Having computed structural, as well as ballistic electronic and phononic transport properties for all structures, we report the thermoelectric properties of the semiconducting ones. We find that 2H phases of four of the studied structures have very promising thermoelectric properties, unlike their 1T phases. The maximum room temperature p-type thermoelectric figure of merit (ZT) of 1.57 is obtained for 2H-HfSe2, which can be as high as 3.30 at T = 800 K. Additionally, 2H-ZrSe2, 2H-ZrTe2, and 2H-HfS2 have considerable ZT values (both nand p-type), that are above 1 at room temperature. The 1T phases of Zr and Hf-based oxides possess relatively high power factors, however their high lattice thermal conductance values limit their ZT values to below 1 at room temperature.

Abstract: This study focuses on optimizing the composition and firing temperature of porcelain tiles using statistical analysis techniques. A full factorial design, including model adequacy checking, analysis of variance, Pareto charts, interaction plots, regression model, and response optimizer is employed. The key factors were the Seger ratios of SiO2/Al2O3, Na2O/K2O, MgO/CaO, and firing temperature. The response variables investigated were bulk density, water absorption, linear shrinkage, coefficient of thermal expansion (at 500 degrees C), and strength. The statistical analysis revealed highly significant results, which were further validated, confirming their reliability for practical use in the production of porcelain tiles. The study demonstrated the effectiveness of utilizing Seger formulas and properties of typical raw materials to accurately predict the final properties of ceramic tiles. By employing SiO2/Al2O3 = 5.2, Na2O/K2O = 1.50, MgO/CaO = 3.0, and firing temperature of 1180 degrees C, optimized properties, such as maximum strength, maximum bulk density, and minimum water absorption, was achieved with a composite desirability of 0.9821.

Abstract: Large area, single layer WS2 has a high potential for use in optoelectrical devices with its high photo-luminescence intensity and low response time. In this work, we demonstrate a systematic study of controlled tungsten disulfide (WS2) monolayer growth using chemical vapor deposition (CVD) technique. With a detailed investigation of process parameters such as H-2 gas inclusion into the main carrier gas, growth temperature and duration, we have gained insight into two-dimensional (2D) WS2 synthesis through controlling the seed formations and the radical vapor density associated with WO3. We confirm that H-2 gas, when included to the carrier gas, is directly involved in WO3 reduction due to its reductive reagent nature, which provides a more effective sulfurization and monolayer formation process. Additionally, by changing the CVD growth configuration, hence, increasing the tungsten related vapor density and confining the reactant radicals, we succeed in realizing larger WS(2 )monolayers, which is still a technological challenge in order to utilize these structures for practical applications. Further optimization of the growth procedure is demonstrated by tuning the growth duration to prevent the excess seed formations and additional layers which will possibly limit the device performance of the monolayer flakes or films when applied.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: High-performance rechargeable batteries are becoming very important for high-end technologies with their ever increasing application areas. Hence, improving the performance of such batteries has become the main bottleneck to transferring high-end technologies to end users. In this study, we propose an argon intercalation strategy to enhance battery performance via engineering the interlayer spacing of honeycomb structures such as graphite, a common electrode material in lithium-ion batteries (LIBs). Herein, we systematically investigated the LIB performance of graphite and hexagonal boron nitride (h-BN) when argon atoms were sent into between their layers by using first-principles density-functional-theory calculations. Our results showed enhanced lithium binding for graphite and h-BN structures when argon atoms were intercalated. The increased interlayer space doubles the gravimetric lithium capacity for graphite, while the volumetric capacity also increased by around 20% even though the volume was also increased. The ab initio molecular dynamics simulations indicate the thermal stability of such graphite structures against any structural transformation and Li release. The nudged-elastic-band calculations showed that the migration energy barriers were drastically lowered, which promises fast charging capability for batteries containing graphite electrodes. Although a similar level of battery promise was not achieved for h-BN material, its enhanced battery capabilities by argon intercalation also support that the argon intercalation strategy can be a viable route to enhance such honeycomb battery electrodes.

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: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.