Materials Science (cond-mat.mtrl-sci)

Abstract: We theoretically study nearly uniform electron models with weak crystalline potentials. In particular, we theorize the modulation of the plane-wave branches at linear regions where multiple Bragg planes intersect. Any such linear intersections involve three or more plane-wave branches diffracted by the periodic potential. Small inter-branch interactions can yield various crossing and anticrossing singularities with promised breakdown of the quadratic approximation, extending alongside the intersection lines. Most of the intersections run in low-symmetric paths in the Brillouin zone and therefore we cannot completely characterize their electronic states with standard band structure plotting methods. The present theory reveals a general mechanism in nearly uniform systems to induce approximately degenerate and linearly continuous Dirac and van-Hove singularities in three dimensions, which may host a variety of anomalous low-energy electronic properties. We apply the theory to a recently discovered high temperature superconductor H$_{3}$S to interpret the enigmatic density-of-state peaking therein.

Abstract: (2D) lateral heterostructures (LH) combining Ti$_2$C and Ta$_2$C MXenes were investigated by means of first-principles calculations. Our structural and elastic properties calculations show that the lateral Ti$_2$C/Ta$_2$C heterostructure results in a 2D material that is stronger than the original isolated MXenes and other 2D monolayers such as germanene or MoS$_2$. The analysis of the evolution of the charge distribution with the size of the LH shows that, for small systems, it tends to distribute homogeneously between the two monolayers, whereas for larger systems electrons tend to accumulate in a region of $\sim$~6 {\AA} around the interface. The work function of the heterostructure, one crucial parameter in the design of electronic nanodevices, is found to be lower than that of some conventional 2D LH. Remarkably, all the heterostructures studied show a very high Curie temperature (between 696 K and 1082 K), high magnetic moments %present in the ferromagnetic ground state, and high magnetic anisotropy energies. These features make (Ti$_2$C)/(Ta$_2$C) lateral heterostructures very suitable candidates for spintronic, photocatalysis, and data storage applications based upon 2D magnetic materials.

Abstract: Analytical gradients of the total energy are provided for local density and generalized-gradient hybrid approximations to generalized Kohn-Sham spin-current density functional theory (SCDFT). An implementation is presented in the public \textsc{crystal} program. It is shown that gradients may be determined analytically, in a two-component framework, including spin-orbit coupling (SOC), with high accuracy. We demonstrate that renormalization of the electron-electron potential by SOC-induced spin-currents can account for considerable modification of crystal structures. In the case of iodine-based molecular crystals, the effect amounts to about half or more of the total modification of the structure by SOC. Such effects necessitate an SCDFT, rather than DFT, formulation, in which exchange-correlation functionals are endowed with an explicit dependence on spin-current densities.

Abstract: Phase-field models are widely employed to simulate microstructure evolution during processes such as solidification or heat treatment. The resulting partial differential equations, often strongly coupled together, may be solved by a broad range of numerical methods, but this often results in a high computational cost, which calls for advanced numerical methods to accelerate their resolution. Here, we quantitatively test the efficiency and accuracy of semi-implicit Fourier spectral-based methods, implemented in Python programming language and parallelized on a graphics processing unit (GPU), for solving a phase-field model coupling Cahn-Hilliard and Allen-Cahn equations. We compare computational performance and accuracy with a standard explicit finite difference (FD) implementation with similar GPU parallelization on the same hardware. For a similar spatial discretization, the semi-implicit Fourier spectral (FS) solvers outperform the FD resolution as soon as the time step can be taken 5 to 6 times higher than afforded for the stability of the FD scheme. The accuracy of the FS methods also remains excellent even for coarse grids, while that of FD deteriorates significantly. Therefore, for an equivalent level of accuracy, semi-implicit FS methods severely outperform explicit FD, by up to 4 orders of magnitude, as they allow much coarser spatial and temporal discretization.

Abstract: PtSe$_2$ is a promising 2D material for nanoelectromechanical sensing and photodetection in the infrared regime. One of its most compelling features is the facile synthesis at temperatures below 500 {\deg}C, which is compatible with current back-end-of-line semiconductor processing. However, this process generates polycrystalline thin films with nanoflake-like domains of 5 to 100 nm size. To investigate the lateral quantum confinement effect in this size regime, we train a deep neural network to obtain an interatomic potential at DFT accuracy and use that to model ribbons, surfaces, nanoflakes, and nanoplatelets of PtSe$_2$ with lateral widths between 5 to 15 nm. We determine which edge terminations are the most stable and find evidence that the electrical conductivity is localized on the edges for lateral sizes below 10 nm. This suggests that the transport channels in thin films of PtSe$_2$ might be dominated by networks of edges, instead of transport through the layers themselves.

Abstract: Machine learning is used to generate empirical pseudopotentials that characterize the local screened interactions in the Kohn-Sham Hamiltonian. Our approach incorporates momentum-range-separated rotation-covariant descriptors to capture crystal symmetries as well as crucial directional information of bonds, thus realizing accurate descriptions of anisotropic solids. Trained empirical potentials are shown to be versatile and transferable such that the calculated energy bands and wave functions without cumbersome self-consistency reproduce conventional ab initio results even for semiconductors with defects, thus fostering faster and faithful data-driven materials researches.

Abstract: Actinide-bearing intermetallics display unusual electronic, magnetic, and physical properties which arise from the complex behavior of their 5$f$ electron orbitals. Temperature ($T$) effects on actinide intermetallics are well studied, but high pressure ($P$) properties and phase stabilities are known for only a handful of compositions. Furthermore, almost no data exist for simultaneous high $P$ and high $T$. We performed ambient-$T$ diamond anvil cell X-ray diffraction experiments to study the behavior of the intermetallic U$_6$Fe upon compression up to 82 GPa. U$_6$Fe remains stable in the tetragonal $I4/mcm$ structure over this pressure range. We also performed ambient $P$, low-$T$ diffraction and heat capacity measurements to constrain U$_6$Fe's thermal behavior. These data were combined with calculations and fitted to a Mie-Gruneisen/Birch-Murnaghan thermal equation of state with the following parameter values at ambient $P$: bulk modulus $B_0$ = 124.0 GPa, pressure derivative $B'_0$ = 5.6, Gruneisen parameter $\Gamma_0$ = 2.028, volume exponent $q$ = 0.934, Debye temperature $\theta_0$ = 175 K, and unit cell volume $V_0$ = 554.4 angstrom$^3$. We report $T$-dependent thermal expansion coefficients and bond lengths of U$_6$Fe, which demonstrate the anisotropic compressibility and negative thermal expansion of the crystallographic $c$ axis. Additionally, density-functional theory calculations indicate increased delocalization of U$_6$Fe bonds at high $P$.

Abstract: We present an \textit{ab-initio} study of the quasi-2D layered perovskite Sr$_3$Hf$_2$O$_7$ com\-pound, performed within the framework of the Density Functional Theory and lattice dynamics analysis. At high temperatures, this compound takes a \textit{I4/mmm} centrosym\-met\-ric structure (S.G. n. 139); as the temperature is lowered, the symmetry is broken into other intermediate polymorphs before reaching the ground state structure, which is the \textit{Cmc2$_1$} ferroelectric phase (S.G. n. 36). One of these intermediate polymorphs is the \textit{Ccce} structural phase (S.G. n. 68). Additionally, we have probed the \textit{C2/c} system (S.G n. 15), which was obtained by following the atomic displacements corresponding to the eigenvectors of the imaginary frequency mode localized at the $\mathbf{\Gamma}$-point of the \textit{Ccce} phase. By observing the enthalpies at low pressures, we found that the \textit{Cmc2$_1$} phase is thermodynamically the most stable. Our results show that the \textit{I4/mmm} and \textit{C2/c} phases never stabilize in the 0-20 GPa range of pressure values. On the other hand, the \textit{Ccce} phase becomes energetically more stable at around 17 GPa, surpassing the \textit{Cmc2$_1$} structure. By considering the effect of entropy and the constant-volume free energies, we observe that the \textit{Cmc2$_1$} polymorph is energetically the most stable phase at low temperature; however, at 350 K the \textit{Ccce} system becomes the most stable. By probing the volume-dependent free energies at 19 GPa, we see that \textit{Ccce} is always the most stable phase between the two structures and also throughout the studied temperature range. When analyzing the phonon dispersion frequencies, we conclude that the \textit{Ccce} system becomes dynamically stable only around 19-20 GPa, and that the \textit{Cmc2$_1$} phase, is metastable up to 30 GPa.