No Arabic abstract
In this work, we report on the mechanical responses and fracture behavior of pristine and defected monolayer 1T-Titanium Disulfide using classical molecular dynamics simulation. We investigated the effect of temperature, strain rate and defect ratio on the uniaxial tensile properties in both armchair and zigzag direction. We found that monolayer TiS2 shows isotropic uniaxial tensile properties except for failure strain which is greater in zigzag direction than armchair direction. We also observed a negative correlation of ultimate tensile strength, failure strain and youngs modulus with temperature and defect ratio. Results depicts that strain rate has no effect on the youngs modulus of monolayer TiS2 but higher strain rate results in higher ultimate tensile strength and failure strain.
In this study, we report the mechanical properties and fracture mechanism of pre-cracked and defected InSe nanosheet samples using molecular dynamics (MD) simulations. We noticed that the failure of pre-cracked and defected InSe nanosheet is governed by brittle type fracture. Armchair directional bonds exhibit a greater resistance for crack propagation relative to the zigzag directional ones. Thus, fracture strength of the pre-cracked sheet is slightly higher for zigzag directional loading than that for armchair. We evaluated the limitation of the applicability of Griffiths criterion for single layer (SL) InSe sheet for nano-cracks as the brittle failure of Griffith prediction demonstrates significant differences with the MD fracture strength. We inspected the effect of temperature on the mechanical properties of the pre-cracked samples of SLInSe. We also discussed the fracture mechanism of both defected and pre-cracked structure at length.
Halide perovskites make efficient solar cells due to their exceptional optoelectronic properties, but suffer from several stability issues. The characterization of the degradation processes is challenging because of the limitations in the spatio-temporal resolution in experiments and the absence of efficient computational methods to study the reactive processes. Here, we present the first effort in developing reactive force fields for large scale molecular dynamics simulations of the phase instability and the defect-induced degradation reactions in inorganic CsPbI$_{3}$. We find that the phase transitions are driven by a combination of the anharmonicity of the perovskite lattice with the thermal entropy. At relatively low temperatures, the Cs cations tend to move away from the preferential positions with good contacts with the surrounding metal halide framework, potentially causing its conversion to a non-perovskite phase. Our simulations of defective structures reveal that, although both iodine vacancies and interstitials are very mobile in the perovskite lattice, the vacancies have a detrimental effect on the stability, initiating the decomposition reactions of perovskites to PbI$_{2}$. Our work puts ReaxFF forward as an effective computational framework to study reactive processes in halide perovskites.
Cerium oxide (ceria, CeO2) is one of the most promising mixed ionic and electronic conducting materials. Previous atomistic analysis has covered widely the effects of substitution on oxygen vacancy migration. However, an in-depth analysis of the role of cation substitution beyond trivalent cations has rarely been explored. Here, we investigate soluble monovalent, divalent, trivalent and tetravalent cation substituents. By combining classical simulations and quantum mechanical calculations, we provide an insight into defect association energies between substituent cations and oxygen vacancies as well as their effects on the diffusion mechanisms. Our simulations indicate that oxygen ionic diffusivity of subvalent cation-substituted systems follows the order Gd>Ca>Na. With the same charge, a larger size mismatch with Ce cation yields a lower oxygen ionic diffusivity, i.e., Na>K, Ca>Ni, Gd>Al. Based on these trends, we identify species that could tune the oxygen ionic diffusivity: we estimate that the optimum oxygen vacancy concentration for achieving fast oxygen ionic transport is 2.5% for GdxCe1-xO2-x/2, CaxCe1-xO2-x and NaxCe1-xO2-3x/2 at 800 K. Remarkably, such a concentration is not constant and shifts gradually to higher values as the temperature is increased. We find that co-substitutions can enhance the impact of the single substitutions beyond that expected by their simple addition. Furthermore, we identify preferential oxygen ion migration pathways, which illustrate the electro-steric effects of substituent cations in determining the energy barrier of oxygen ion migration. Such fundamental insights into the factors that govern the oxygen diffusion coefficient and migration energy would enable design criteria to be defined for tuning the ionic properties of the material, e.g., by co-doping.
Pentadiamond is a recently proposed new carbon allotrope consisting of a network of pentagonal rings where both sp$^2$ and sp$^3$ hybridization are present. In this work we investigated the mechanical and electronic properties, as well as, the thermal stability of pentadiamond using DFT and fully atomistic reactive molecular dynamics (MD) simulations. We also investigated its properties beyond the elastic regime for three different deformation modes: compression, tensile and shear. The behavior of pentadiamond under compressive deformation showed strong fluctuations in the atomic positions which are responsible for the strain softening at strains beyond the linear regime, which characterizes the plastic flow. As we increase temperature, as expected, Youngs modulus values decrease, but this variation (up to 300 K) is smaller than 10% (from 347.5 to 313.6 GPa), but the fracture strain is very sensitive, varying from $sim$44% at 1K to $sim$5% at 300K.
The electronic and thermal transport properties have been systematically investigated in monolayer C$_4$N$_3$H with first-principles calculations. The intrinsic thermal conductivity of monolayer C$_4$N$_3$H was calculated coupling with phonons Boltzmann transport equation. For monolayer C$_4$N$_3$H, the thermal conductivity (k{appa}) (175.74 and 157.90 W m-1K-1 with a and b-plane, respectively) is significantly lower than that of graphene (3500 Wm$^{-1}$K$^{-1}$) and C3N(380 Wm$^{-1}$K$^{-1}$). Moreover, it is more than the second time higher than C$_2$N (82.88 Wm$^{-1}$K$^{-1}$) at 300 K. Furthermore, the group velocities, relax time, anharmonicity, as well as the contribution from different phonon branches, were thoroughly discussed in detail. A comparison of the thermal transport characters among 2D structure for monolayer C$_4$N$_3$H, graphene, C$_2$N and C$_3$N has been discussed. This work highlights the essence of phonon transport in new monolayer material.