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: 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: 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: 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: Using density functional theory, we obtain the mechanical and thermal properties of MX2 monolayers (where M = Cr, Mo, W and X = O, S, Se, Te). The C-centered phonon frequencies (i.e., A(1), A(2)'', E ', and E ''), relative frequency values of A(1), and E ' modes, and mechanical properties (i.e., elastic constants, Young modulus, and Poisson's ratio) display a strong dependence on the type of metal and chalcogenide atoms. In each chalcogenide (metal) group, transition-metal dichalcogenides (TMDCs) with W (O) atom are found to be much stiffer. Consistent with their stability, the thermal expansion of lattice constants for TMDCs with O (Te) is much slower (faster). Furthermore, in a heterostructure of these materials, the difference of the thermal expansion of lattice constants between the individual components becomes quite tiny over the whole temperature range. The calculated mechanical and thermal properties show that TMDCs are promising materials for heterostructures. (C) 2014 AIP Publishing LLC.

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