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: 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: Metallic contacts are critical components of electronic devices and the presence of a large Schottky barrier is detrimental for an optimal device operation. Here, we show by using first-principles calculations that a self-assembled monolayer (SAM) of polar molecules between the metal electrode and MoSe2 monolayer is able to convert the Schottky contact into an almost Ohmic contact. We choose -CH3 and -CF3 terminated short-chain alkylthiolate (i.e. SCH3 and fluorinated alkylthiolates (SCF3)) based SAMs to test our approach. We consider both high (Au) and low (Sc) work function metals in order to thoroughly elucidate the role of the metal work function. In the case of Sc, the Fermi level even moves into the conduction band of the MoSe2 monolayer upon SAM insertion between the metal surface and the MoSe2 monolayer, and hence possibly switches the contact type from Schottky to Ohmic. The usual Fermi level pinning at the metal-transition metal dichalcogenide (TMD) contact is shown to be completely removed upon the deposition of a SAM. Systematic analysis indicates that the work function of the metal surface and the energy level alignment between the metal electrode and the TMD monolayer can be tuned significantly by using SAMs as a buffer layer. These results clearly indicate the vast potential of the proposed interface engineering to modify the physical and chemical properties of MoSe2.

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: 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 adsorption and diffusion of Na, K, and Ca atoms on MXene/graphene heterostructures of MXene systems Sc2C(OH)(2), Ti2CO2, and V2CO2 are systematically investigated by using first-principles methods. We found that alkali metal intercalation is energetically favorable and thermally stable for Ti2CO2/graphene and V2CO2/graphene heterostructures but not for Sc2C(OH)(2). Diffusion kinetics calculations showed the advantage of MXene/graphene heterostructures over sole MXene systems as the energy barriers are halved for the considered alkali metals. Low energy barriers are found for Na and K ions, which are promising for fast charge/discharge rates. Calculated voltage profiles reveal that estimated high capacities can be fully achieved for Na ion in V2CO2/graphene and Ti2CO2/graphene heterostructures. Our results indicate that Ti2CO2/graphene and V2CO2/graphene electrode materials are very promising for Na ion battery applications. The former could be exploited for low voltage applications while the latter will be more appropriate for higher voltages.

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

Abstract: In this study, by means of classical molecular dynamics simulations, we investigate the thermal-transport properties of hexagonal single-layer, zinc-blend, and wurtzite phases of BN, AlN, and GaN crystals, which are very promising for the application and design of high-quality electronic devices. With this in mind, we generate fully transferable Tersoff-type empirical interatomic potential parameter sets by utilizing an optimization procedure based on particle-swarm optimization. The predicted thermal properties as well as the structural, mechanical, and vibrational properties of all materials are in very good agreement with existing experimental and first-principles data. The impact of isotopes on thermal transport is also investigated and between approximately 10 and 50% reduction in phonon thermal transport with random isotope distribution is observed in BN and GaN crystals. Our investigation distinctly shows that the generated parameter sets are fully transferable and very useful in exploring the thermal properties of systems containing these nitrides.

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