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: 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: 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: There is an increasing need to understand interfaces between two-dimensional materials to realize an energy efficient boundary with low contact resistance and small heat dissipation. In this respect, we investigated the impact of charge and substitutional atom doping on the electronic transport properties of the hybrid metallic-semiconducting lateral junctions, formed between metallic (1T and 1T(d)) and semiconducting (1H) phases of MoS2 by means of first-principles and non-equilibrium Green function formalism based calculations. Our results clearly revealed the strong influence of the type of interface and crystallographic orientation of the metallic phase on the transport properties of these systems. The Schottky barrier height, which is the dominant mechanism for contact resistance, was found to be as large as 0.63 eV and 1.19 eV for holes and electrons, respectively. We found that armchair interfaces are more conductive as compared to zigzag termination due to the presence of the metallic Mo zigzag chains that are directed along the transport direction. In order to manipulate these barrier heights we investigated the influence of electron doping of the metallic part (i.e. 1T(d) -MoS2). We observed that the Fermi level of the hybrid system moves towards the conduction band of semiconducting 1H-MoS2 due to filling of 4d-orbital of metallic MoS2, and thus the Schottky barrier for electrons decreases considerably. Besides electron doping, we also investigated the effect of substitutional doping of metallic MoS2 by replacing Mo atoms with either Re or Ta. Due to its valency, Re (Ta) behaves as a donor (acceptor) and reduces the Schottky barrier for electrons (holes). Since Re and Ta based transition metal dichalcogenides crystallize in either the 1T(d) or 1T phase, substitutional doping with these atom favors the stabilization of the 1T(d) phase of MoS2. Co-doping of hybrid structure results in an electronic structure, which facilities easy dissociation of excitons created in the 1H part.

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 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: 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: MXenes are an emerging class of two-dimensional materials, which in their thinnest limit consist of a monolayer of carbon or nitrogen (X) sandwiched between two transition metal (M) layers. We have systematically searched for superconductivity among MXenes for a range of transition metal elements, based on a full first-principles characterization in combination with the Eliashberg formalism. Thus, we identified six superconducting MXenes: three carbides (Mo2C, W2C and Sc2C) and three nitrides (Mo2N, W2N and Ta2N). The highest critical temperature of similar to 16 K is found in Mo2N, for which a successful synthesis method has been established [Urbankowskiet al.,Nanoscale, 2017,9, 17722-17730]. Moreover, W2N presents a novel case of competing superconducting and charge density wave phases.

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 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 the growth and optimization of ultralow loss Si-based Al2O3 planar waveguides, which have a high potential to boost the performance of rare-earth ion doped waveguide devices operating at visible and C-band wavelength ranges. The planar waveguide structures are grown using thermal atomic layer deposition. Systematic characterization of the obtained thin films is performed by spectroscopic ellipsometry, X-ray diffraction, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy analyses, and the optimum parameters are identified. The optical loss measurements for both transverse electric (TE) and transverse magnetic polarized light at 633, 829, and 1549 nm are performed. The lowest propagation loss value of 0.04 +/- 0.02 dB/cm for the Al2O3 waveguides for TE polarization at 1549 nm is demonstrated.

Abstract: X-ray diffraction and energy dispersive x-ray spectroscopic analyses showed a natural galena (PbS) crystal from Freiberg in Saxony (Germany) to be a single phase specimen [rock salt (NaCl) structure type, space group Fm3m, a = 5.932(1) angstrom] with stoichiometric composition and an enhanced dislocation density (8 approximate to 1011 cm-2). The latter parameter leads to an increase of the electrical resistivity in the high-temperature regime, as well as to the appearance of phonon resonance with a characteristic frequency coPR = 3.8(1) THz. Being in the same range (i.e., 3-5.5 THz) with the sulfur optical modes of highest group velocities, it results in a drastic reduction (by similar to 75%) of thermal conductivity (K) at lower temperatures (i.e., < 100 K), as well as in the appearance of a characteristic minimum in K at T approximate to 30 K. Furthermore, the studied galena is characterized by phonon-drag behavior and by temperature dependent switch of the charge carrier scattering mechanism regime (i.e., scattering on dislocations for T < 100 K, on acoustic phonons for 100 K < T < 170 K and on both acoustic and optical phonons for 170 K < T < 300 K). The combined theoretical calculation and optical spectroscopic study confirm this mineral to be a direct gap degenerate semiconductor. The possible origins of the second-order Raman spectrum are discussed.

Abstract: Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalco-genides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory-based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber-type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectralenergy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.