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: 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: 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: 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: 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: 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: Single-layer, large-scale two-dimensional material growth is still a challenge for their wide-range usage. Therefore, we carried out a comprehensive study of monolayer MoS2 growth by CVD investigating the influence of growth zone configuration and precursors ratio. We first compared the two commonly used approaches regarding the relative substrate and precursor positions, namely, horizontal and face-down configurations where facedown approach is found to be more favorable to obtain larger flakes under identical growth conditions. Secondly, we used different types of substrate holders to investigate the influence of the Mo and S vapor confinement on the resulting diffusion environment. We suggest that local changes of the S to Mo vapor ratio in the growth zone is a key factor for the change of shape, size and uniformity of the resulting MoS2 formations, which is also confirmed by performing depositions under different precursor ratios. Therefore, to obtain continuous monolayer films, the S to Mo vapor ratio is needed to be kept within a certain range throughout the substrate. As a conclusion, we obtained monolayer triangles with a side length of 90 mu m and circles with a diameter of 500 mu m and continuous films with an area of 85 0 mu m x 1 cm when the S-to-Mo vapor ratio is optimized. (C) 2017 The Japan Society of Applied Physics

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: 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: 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: The structural stability of 2D transition metal dichalcogenide (TMD) formations is of particular importance for their reliable device performance in nano-electronics and opto-electronics. Recent observations show that the CVD-grown TMD monolayers are likely to encounter stability problems such as cracking or fracturing when they are kept under ambient conditions. Here, two different growth configurations are investigated and a favorable growth geometry is proposed, which also sheds light onto the growth mechanism and provides a solution for the stability and fracture formation issues for TMDs specifically for MoS2 monolayers. It is shown that 18 months naturally and thermally aged MoS2 monolayer flakes grown using specifically developed conditions, retain their stability. To understand the mechanism of the structural deterioration, two possible effective mechanisms, S vacancy defects and growth-induced tensile stress, are assessed by the first principle calculations where the role of S vacancy defects in obtaining oxidation resistant MoS2 monolayer flakes is revealed to be rather more critical. Hence, these simulations, time-dependent observations and thermal aging experiments show that durability and stability of 2D MoS2 flakes can be controlled by CVD growth configuration.

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: Using first-principles calculations, we evaluate the electrochemical performance of heterostructures made up of Ti2CO2 and chemically modified graphene for Li batteries. We find that heteroatom doping and molecule intercalation have a significant impact on the storage capacity and Li migration barrier energies. While N and S doping do not improve the storage capacity, B doping together with molecule interaction make it possible to intercalate two layers of Li, which stick separately to the surface of Ti2CO2 and B-doped graphene. The calculated diffusion-barrier energies (E-diff), which are between 0.3 and 0.4 eV depending on Li concentration, are quite promising for fast charge and discharge rates. Besides, the predicted E-diff as much as 2 eV for the diffusion of the Li atom from the Ti2CO2 surface to the B-doped graphene surface significantly suppresses the interlayer Li migration, which diminishes the charge and discharge rates. The calculated volume and lattice parameter changes indicate that Ti2CO2/graphene hybrid structures exhibit cyclic stability against Li loading and unloading. Consequently, first-principles calculations we perform evidently highlight the favorable effect of molecular intercalation on the capacity improvement of ion batteries.

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.