Search

Altermagnet (AM) is a new proposed magnetic state with collinear antiferromagnetic ground state but presents some transport properties that were only believed to exist in ferromagnets or non-collinear antiferromagnets. To have a comprehensive understanding of the transport properties of AMs, especially from the experimental point of view, a promising altermagnetic metal is crucial. In all the proposed altermagnetic metals, RuO$_2$ has a special position, since it is the first proposed AM with the largest spin splitting and several important altermagnetism featured experiments were first performed based on it. However, a very recent report based on sensitive muon-spin measurements suggest a super small local magnetization from Ru, i.e. a nonmagnetic ground state in RuO$_2$. Therefore, a determination of the existence of the altermagnetic ground state is the basic starting point for all the previously altermagnetic transport properties in RuO$_2$. In this work, we propose to identify its magnetic ground state from the Fermi surface (FS) via the electronic transport property of quantum oscillation (QO). We systematically analyzed the FSs of RuO$_2$ in both nonmagnetic and altermagnetic states via first principles calculations. Our work should be helpful for future experiments on QO measurements to confirm its ground state by the interplay between transport measurements and computations., Comment: 11 pages, 6 figures

Nonlinear Hall effect (NHE) can be generated via Berry curvature dipole (BCD) on nonequilibrium Fermi surface in a non-magnetic system without inversion symmetry.To achieve a large BCD, strong local Berry curvatures and their variation with respect to momentum are necessary and hence topological materials with strong inter-band coupling emerge as promising candidates. In this study, we propose a switchable and robust BCD in the newlydiscovered dual quantum spin Hall insulator (QSHI) $\text{TaIrTe}_4$ by applying out-of-plane electric fields. Switchable BCD could be found along with topological phase transitions or insulator-metal transition in the primitive cell and CDW phases of $\text{TaIrTe}_4$ monolayer. This work presents an instructive strategy for achieving a switchable and robust BCD within dual QSHIs, which should be helpful for designing the NHE-based THz radiations detector.

Dual quantum spin Hall insulator (QSHI) is a newly discovered topological state in the 2D material TaIrTe$_4$, exhibiting both a traditional $Z_2$ band gap at charge neutrality point and a van Hove singularity (VHS) induced correlated $Z_2$ band gap with weak doping. Inspired by the recent progress in theoretical understanding and experimental measurements, we predicted a promising dual QSHI in the counterpart material of the NbIrTe4 monolayer by first-principles calculations. In addition to the well-known band inversion at the charge neutrality point, two new band inversions were found after CDW phase transition when the chemical potential is near the VHS, one direct and one indirect $Z_2$ band gap. The VHS-induced non-trivial band gap is around 10 meV, much larger than that from TaIrTe$_4$. Furthermore, since the new generated band gap is mainly dominated by the $4d$ orbitals of Nb, electronic correlation effects should be relatively stronger in NbIrTe$_4$ as compared to TaIrTe$_4$. Therefore, the dual QSHI state in the NbIrTe$_4$ monolayer is expected to be a good platform for investigating the interplay between topology and correlation effects.

Topological insulating states in two-dimensional (2D) materials are ideal systems to study different types of quantized response signals due to their in gap metallic states. Very recently, the quantum spin Hall (QSH) effect was discovered in monolayer $\text{TaIrTe}_4$ via the observation of quantized longitudinal conductance that rarely exists in other 2D topological insulators. The non-trivial $Z_2$ topological charges can exist at both charge neutrality point and the van Hove singularity point with correlation effect induced band gap. Based on this model 2D material, we studied the switch of quantized signals between longitudinal conductance and transversal Hall conductance via tuning external magnetic field. In $Z_2$ topological phase of monolayer $\text{TaIrTe}_4$, the zero Chern number can be understood as 1-1=0 from the double band inversion from spin-up and spin-down channels. After applying a magnetic field perpendicular to the plane, the Zeeman split changes the band order for one branch of the band inversion from spin-up and spin-down channels, along with a sign charge of the Berry phase. Then the net Chern number of 1-1=0 is tuned to 1+1=2 or -1-1=-2 depending on the orientation of the magnetic field. The quantized signal not only provides another effective method for the verification of topological state in monolayer $\text{TaIrTe}_4$, but also offers a strategy for the utilization of the new quantum topological states based on switchable quantized responses.

Reconstructive phase transitions involving breaking and reconstruction of primary chemical bonds are ubiquitous and important for many technological applications. In contrast to displacive phase transitions, the dynamics of reconstructive phase transitions are usually slow due to the large energy barrier. Nevertheless, the reconstructive phase transformation from $\beta$- to $\lambda$-Ti$_3$O$_5$ exhibits an ultrafast and reversible behavior. Despite extensive studies, the underlying microscopic mechanism remains unclear. Here, we discover a kinetically favorable in-plane nucleated layer-by-layer transformation mechanism through metadynamics and large-scale molecular dynamics simulations. This is enabled by developing an efficient machine learning potential with near first-principles accuracy through an on-the-fly active learning method and an advanced sampling technique. Our results reveal that the $\beta$-$\lambda$ phase transformation initiates with the formation of two-dimensional nuclei in the $ab$-plane and then proceeds layer-by-layer through a multistep barrier-lowering kinetic process via intermediate metastable phases. Our work not only provides important insight into the ultrafast and reversible nature of the $\beta$-$\lambda$ transition, but also presents useful strategies and methods for tackling other complex structural phase transitions., Comment: 26 pages,23 figures (including Supporting Information)

The recent experimental claim of room-temperature ambient-pressure superconductivity in a Cu-doped lead-apatite (LK-99) has ignited substantial research interest in both experimental and theoretical domains. Previous density functional theory (DFT) calculations with the inclusion of an on-site Hubbard interaction $U$ consistently predict the presence of flat bands crossing the Fermi level. This is in contrast to DFT plus dynamical mean field theory calculations, which reveal the Mott insulating behavior for the stoichiometric Pb$_{9}$Cu(PO$_4$)$_6$O compound. Nevertheless, the existing calculations are all based on the $P6_3/m$ structure, which is argued to be not the ground-state structure. Here, we revisit the electronic structure of Pb$_{9}$Cu(PO$_4$)$_6$O with the energetically more favorable $P\bar{3}$ structure, fully taking into account electronic symmetry breaking. We examine all possible configurations for Cu substituting the Pb sites. Our results show that the doped Cu atoms exhibit a preference for substituting the Pb2 sites than the Pb1 sites. In both cases, the calculated substitutional formation energies are large, indicating the difficulty in incorporating Cu at the Pb sites. We find that most of structures with Cu at the Pb2 site tend to be insulating, while the structures with both two Cu atoms at the Pb1 sites (except one configuration) are predicted to be metallic by DFT+$U$ calculations. However, when accounting for the electronic symmetry breaking, some Cu-doped configurations previously predicted to be metallic (including the structure studied in previous DFT+$U$ calculations) become insulating. Our work highlights the importance of symmetry breaking in obtaining correct electronic state for Pb$_{9}$Cu(PO$_4$)$_6$O, thereby reconciling previous DFT+$U$ and DFT+DMFT calculations., Comment: 19 pages, 9 figures (including Supplementary Material)

Recently, Lee et al. reported the experimental discovery of room-temperature ambient-pressure superconductivity in a Cu-doped lead-apatite (LK-99) (arXiv:2307.12008, arXiv:2307.12037). Remarkably, the superconductivity persists up to 400 K at ambient pressure. Despite strong experimental evidence, the electronic structure of LK-99 has not yet been studied. Here, we investigate the electronic structures of LK-99 and its parent compound using first-principles calculations, aiming to elucidate the doping effects of Cu. Our results reveal that the parent compound Pb$_{10}$(PO$_4$)$_6$O is an insulator, while Cu doping induces an insulator-metal transition and thus volume contraction. The band structures of LK-99 around the Fermi level are featured by a half-filled flat band and a fully-occupied flat band. These two flat bands arise from both the $2p$ orbitals of $1/4$-occupied O atoms and the hybridization of the $3d$ orbitals of Cu with the $2p$ orbitals of its nearest-neighboring O atoms. Interestingly, we observe four van Hove singularities on these two flat bands. Furthermore, we show that the flat band structures can be tuned by including electronic correlation effects or by doping different elements. We find that among the considered doping elements (Ni, Cu, Zn, Ag, and Au), both Ni and Zn doping result in the gap opening, whereas Au exhibits doping effects more similar to Cu than Ag. Our work provides a foundation for future studies on the role of unique electronic structures of LK-99 in superconductivity.

Superhard materials with good fracture toughness have found wide industrial applications, which necessitates the development of accurate hardness and fracture toughness models for efficient materials design. Although several macroscopic models have been proposed, they are mostly semiempirical based on prior knowledge or assumptions, and obtained by fitting limited experimental data. Here, through an unbiased and explanatory symbolic regression technique, we built a macroscopic hardness model and fracture toughness model, which only require shear and bulk moduli as inputs. The developed hardness model was trained on an extended dataset, which not only includes cubic systems, but also contains non-cubic systems with anisotropic elastic properties. The obtained models turned out to be simple, accurate, and transferable. Moreover, we assessed the performance of three popular deep learning models for predicting bulk and shear moduli, and found that the crystal graph convolutional neural network and crystal explainable property predictor perform almost equally well, both better than the atomistic line graph neural network. By combining the machine-learned bulk and shear moduli with the hardness and fracture toughness prediction models, potential superhard materials with good fracture toughness can be efficiently screened out through high-throughput calculations.

Excellent high-temperature mechanical properties of Ni-based single crystal superalloys (NSCSs) are attributed to the yield strength anomaly of Ni$_{3}$Al that is intimately related to generalized stacking fault energies (GSFEs). Therefore, clarifying the effects of alloying elements on the GSFEs is of great significance for alloys design. Here, by means of ab initio density functional theory calculations, we systematically calculated the GSFEs of different slip systems of Ni and Ni$_{3}$Al without and with alloying elements using the alias shear method. We obtained that for Ni, except for magnetic elements Mn, Fe, and Co, most of alloying elements decrease the unstable stacking fault energy ($\gamma_{usf}$) of the $[01\bar{1}](111)$ and $[11\bar{2}](111)$ slip systems and also decrease the stable stacking fault energy ($\gamma_{sf}$) of the $[11\bar{2}](111)$ slip system. For Ni$_{3}$Al, most of alloying elements in groups IIIB-VIIB show a strong Al site preference. Except for Mn and Fe, the elements in groups VB-VIIB and the first column of group VIII increase the values of $\gamma_{usf}$ of different slip systems of Ni$_{3}$Al. On the other hand, the elements in groups IIIB-VIIB also increase the value of $\gamma_{sf}$. We found that Re is an excellent strengthening alloying element that significantly increases the slip barrier of the tailing slip process for Ni, and also enhances the slip barrier of the leading slip process of three slip systems for Ni$_{3}$Al. W and Mo exhibit similar effects as Re. We predicted that Os, Ru, and Ir are good strengthening alloying elements as well, since they show the strengthening effects on both the leading and tailing slip process for Ni and Ni$_{3}$Al.

Recently, near-ambient superconductivity has been experimentally evidenced in a nitrogen-doped lutetium hydride by Dasenbrock-Gammon \emph{et al.} [Nature 615, 244 (2023)], which yields a remarkable maximum $T_c$ of 294 K at just 10 kbar. However, due to the difficulty of x-ray diffraction (XRD) in identifying light elements such as hydrogen and nitrogen, the crystal structure of the superconductor remains elusive, in particular for the actual stoichiometry of hydrogen and nitrogen and their atomistic positions. This holds even for its parent structure. Here, we set out to address this issue by performing a thorough density functional theory study on the structural, electronic, dynamical, and optical properties of lutetium hydrides. Through thermal and lattice dynamic analysis as well as XRD and superconductor color comparisons, we unambiguously clarified that the parent structures are a mixture of dominant LuH$_2$ phase of the CaF$_2$-type (instead of originally proposed LuH$_3$ structure of $Fm\bar{3}m$ space group) and minor LuH phase of the NaCl-type., Comment: 4 pages, 3 figures

Ferroelectric control of two-dimensional magnetism is promising in fabricating electronic devices with high speed and low energy consumption. The newly discovered layered MnBi$_2$Te$_4$(Bi$_2$Te$_3$)$_n$ and their Sb counterparts exhibit A-type antiferromagnetism with intriguing topological properties. Here, we propose to obtain tunable magnetic multistates in their thin films by ferroelectrically manipulating the interlayer magnetic couplings (IMCs) based on the Heisenberg model and first-principles calculations. Our strategy relies on that interfacing the thin films with appropriate ferroelectric materials can switch on/off an interlayer hopping channel between Mn-$e_g$ orbitals as the polarizations reversed, thus resulting in a switchable interlayer antiferromagnetism-to-ferromagnetism transition. On the other hand, the interface effect leads to asymmetric energy barrier heights for the two polarization states. These properties allow us to build ferroelectrically switchable triple and quadruple magnetic states with multiple Chern numbers in thin films. Our study reveals that ferroelectrically switchable magnetic and topological multistates in MnBi$_2$Te$_4$ family can be obtained by rational design for multifunctional electronic devices, which can also be applied to other two-dimensional magnetic materials., Comment: 8.2 pages and 3 figures

Over the past decade, utilizing spin currents in the linear response of electric field to manipulate magnetization states via spin-orbit torques (SOTs) is one of the core concepts for realizing a multitude of spintronic devices. Besides the linear regime, recently, nonlinear charge-spin conversion under the square of electric field has been recognized in a wide variety of materials with nontrivial spin textures, opening an emerging field of nonlinear spintronics. Here, we report the investigation of both linear and nonlinear charge-spin conversion in one single topological semimetal SrIrO3(110) thin film that hosts strong spin-orbit coupling and nontrivial spin textures in the momentum space. In the nonlinear regime, the observation of crystalline direction dependent response indicates the presence of anisotropic surface states induced spin-momentum locking near the Fermi level. Such anisotropic spin textures also give rise to spin currents in the linear response regime, which mainly contributes to the fieldlike SOT component. Our work demonstrates the power of combination of linear and nonlinear approaches in understanding and utilizing charge-spin conversion in topological materials., Comment: 18 pages, 5 figures

Adsorption of carbon monoxide (CO) on transition-metal surfaces is a prototypical process in surface sciences and catalysis. Despite its simplicity, it has posed great challenges to theoretical modeling. Pretty much all existing density functionals fail to accurately describe surface energies, CO adsorption site preference, as well as adsorption energies simultaneously. Although the random phase approximation (RPA) cures these density functional theory failures, its large computational cost makes it prohibitive to study the CO adsorption for any but the simplest ordered cases. Here, we address these challenges by developing a machine-learned force field (MLFF) with near RPA accuracy for the prediction of coverage-dependent adsorption of CO on the Rh(111) surface through an efficient on-the-fly active learning procedure and a $\Delta$-machine learning approach. We show that the RPA-derived MLFF is capable to accurately predict the Rh(111) surface energy, CO adsorption site preference as well as adsorption energies at different coverages that are all in good agreement with experiments. Moreover, the coverage-dependent ground-state adsorption patterns and adsorption saturation coverage are identified., Comment: 6 pages, 2 figures

Two-dimensional (2D) materials have gained lots of attention due to the potential applications. In this work, we propose that based on first-principles calculations, the (2$\times$2) patterned PtTe$_2$ monolayer with kagome lattice formed by the well-ordered Te vacancy (PtTe$_{1.75}$) hosts large spin Hall conductivity (SHC) and excellent hydrogen evolution reaction (HER) activity. The unconventional nature relies on the $A1@1b$ band representation (BR) of the highest valence band without SOC. The large SHC comes from the Rashba spin-orbit coupling (SOC) in the noncentrosymmetric structure induced by the Te vacancy. Even though it has a metallic SOC band structure, the $\mathbb Z_2$ invariant is well defined due to the existence of the direct band gap and is computed to be nontrivial. The calculated SHC is as large as 1.25$\times 10^3 \frac{\hbar}{e} (\Omega~cm)^{-1}$ at the Fermi level ($E_F$). By tuning the chemical potential from $E_F-0.3$ to $E_F+0.3$ eV, it varies rapidly and monotonically from $-1.2\times 10^3$ to 3.1$\times 10^3 \frac{\hbar}{e} (\Omega~cm)^{-1}$. In addition, we also find the Te vacancy in the patterned monolayer can induce excellent HER activity. Our results not only offer a new idea to search 2D materials with large SHC, i.e., by introducing inversion-symmetry breaking vacancies in large SOC systems, but also provide a feasible system with tunable SHC (by applying gate voltage) and excellent HER activity.

Antiferromagnets (AFMs) with anomalous quantum responses have lead to new progress for the understanding of their magnetic and electronic structures from symmetry and topology points of view. Two typical topological states are the collinear antiferromagnetic Weyl semimetal (WSM) and Weyl nodal loop semimetal (WNLSM). In comparison with the counterparts in ferromagnets and non-collinear AFMs, the WSMs and WNLSMs in collinear AFMs are still waiting for experimental verification. In this work, we theoretically predicted the coexistence of Weyl points (WPs) and Weyl nodal loops (WNLs) in transition metal oxide RuO2. Owing to the small magnetocrystalline anisotropy energy, the WPs and WNLs can transform to each other via tuning the Neel vector. Moreover, since the WPs are very close to Fermi level and the WNLs are even crossing Fermi level, the topological states in RuO2 can be easily probed by photoemission and STM methods. Our result provides a promising material platform for the study of WSM and WNLSM states in collinear AFMs., Comment: 6 pages, 4 figures

According to topological quantum chemistry, a class of electronic materials have been called obstructed atomic insulators (OAIs), in which a portion of valence electrons necessarily have their centers located on some empty $\textit{Wyckoff}$ positions without atoms occupation in the lattice. The obstruction of centering these electrons coinciding with their host atoms is nontrivial and results in metallic boundary states when the boundary is properly cut. Here, on basis of first-principles calculations in combination with topological quantum chemistry analysis, we propose two dimensional MA$_2$Z$_4$ (M = Cr, Mo and W; A = Si and Ge, Z = N, P and As) monolayer family are all OAIs. A typical case is the recently synthesized MoSi$_2$N$_4$. Although it is a topological trivial insulator with the occupied electronic states being integer combination of elementary band representations, it has valence electrons centering empty $\textit{Wyckoff}$ positions. It exhibits unique OAI-induced metallic edge states along the (1$\bar{1}$0) edge of MoSi$_2$N$_4$ monolayer and the in-gap corner states at three vertices of certain hexagonal nanodisk samples respecting C$_3$ rotation symmetry. The readily synthesized MoSi$_2$N$_4$ is quite stable and has a large bulk band gap of 1.94 eV, which makes the identification of these edge and corner states most possible for experimental clarification., Comment: 13 pages, 10 figures and 6 tables

With the miniaturization and integration of nanoelectronic devices, efficient heat removal becomes a key factor affecting the reliable operation of the nanoelectronic device. With the high intrinsic thermal conductivity, good mechanical flexibility, and precisely controlled growth, two-dimensional (2D) materials are widely accepted as ideal candidates for thermal management materials. In this work, by solving the phonon Boltzmann transport equation (BTE) based on first-principles calculations, we comprehensively investigated the thermal conductivity of novel 2D layered MSi$_2$N$_4$ (M = Mo, W). Our results point to competitive thermal conductivities (162 W/mK) of monolayer MoSi$_2$N$_4$, which is around two times larger than that of WSi$_2$N$_4$ and seven times larger than that of silicene despite their similar non-planar structures. It is revealed that the high thermal conductivity arises mainly from its large group velocity and low anharmonicity. Our result suggests that MoSi$_2$N$_4$ could be a potential candidate for 2D thermal management materials.

The realization of multifunctional two-dimensional (2D) materials is fundamentally intriguing, such as combination of piezoelectricity with topological insulating phase or ferromagnetism. In this work, a Janus monolayer $\mathrm{SrAlGaSe_4}$ is built from 2D $\mathrm{MA_2Z_4}$ family with dynamic, mechanical and thermal stabilities, which is piezoelectric due to lacking inversion symmetry. The unstrained $\mathrm{SrAlGaSe_4}$ monolayer is a narrow gap normal insulator (NI) with spin orbital coupling (SOC). However, the NI to topological insulator (TI) phase transition can be induced by the biaxial strain, and a piezoelectric quantum spin Hall insulator (PQSHI) can be achieved. More excitingly, the phase transformation point is only about 1.01 tensile strain, and nontrivial band topology can hold until considered 1.16 tensile strain. Moreover, a Rashba spin splitting in the conduction bands can exit in PQSHI due to the absence of a horizontal mirror symmetry and the presence of SOC. For monolayer $\mathrm{SrAlGaSe_4}$, both in-plane and much weak out-of-plane piezoelectric polarizations can be induced with a uniaxial strain applied. The calculated piezoelectric strain coefficients $d_{11}$ and $d_{31}$ of monolayer $\mathrm{SrAlGaSe_4}$ are -1.865 pm/V and -0.068 pm/V at 1.06 tensile strain as a representative TI. In fact, many PQSHIs can be realized from 2D $\mathrm{MA_2Z_4}$ family. To confirm that, similar to $\mathrm{SrAlGaSe_4}$, the coexistence of piezoelectricity and topological orders can be realized by strain (about 1.04 tensile strain) in the $\mathrm{CaAlGaSe_4}$ monolayer. Our works suggest that Janus monolayer $\mathrm{SrAlGaSe_4}$ is a pure 2D system for PQSHI, enabling future studies exploring the interplay between piezoelectricity and topological orders, which can lead to novel applications in electronics and spintronics., Comment: 9 pages,10 figures

The recently experimentally synthesized monolayer $\mathrm{MoSi_2N_4}$ and $\mathrm{WSi_2N_4}$ (\textcolor[rgb]{0.00,0.00,1.00}{Science 369, 670-674 (2020})) lack inversion symmetry, which allows them to become piezoelectric. In this work, based on ab initio calculations, we report structure effect on intrinsic piezoelectricity in septuple-atomic-layer $\mathrm{MSi_2N_4}$ (M=Mo and W), and six structures ($\alpha_i$ ($i$=1 to 6)) are considered with the same space group.It is found that $\mathrm{MSi_2N_4}$ (M=Mo and W) with $\alpha_i$ ($i$=1 to 6) all are indirect band gap semiconductors. Calculated results show that $\mathrm{MoSi_2N_4}$ and $\mathrm{WSi_2N_4}$ monolayers have the same structural dependence on piezoelectric strain and stress coefficients ($d_{11}$ and $e_{11}$), together with the ionic and electronic contributions to $e_{11}$.Finally, we investigate the intrinsic piezoelectricity of monolayer $\mathrm{MA_2Z_4}$ (M=Cr, Mo and W; A=Si and Ge; Z=N and P) with $\alpha_1$ and $\alpha_2$ phases expect $\mathrm{CrGe_2N_4}$, because they all are semiconductors and their enthalpies of formation between $\alpha_1$ and $\alpha_2$ phases are very close. The most important result is that monolayer $\mathrm{MA_2Z_4}$ containing P atom have more stronger piezoelectric polarization than one including N atom. The largest $d_{11}$ among $\mathrm{MA_2N_4}$ materials is 1.85 pm/V, which is close to the smallest $d_{11}$ of 1.65 pm/V in $\mathrm{MA_2P_4}$ monolayers. For $\mathrm{MA_2P_4}$, the largest $d_{11}$ is up to 6.12 pm/V. Among the 22 monolayers, $\alpha_1$-$\mathrm{CrSi_2P_4}$, $\alpha_1$-$\mathrm{MoSi_2P_4}$, $\alpha_1$-$\mathrm{CrGe_2P_4}$, $\alpha_1$-$\mathrm{MoGe_2P_4}$ and $\alpha_2$-$\mathrm{CrGe_2P_4}$ have large $d_{11}$, which are greater than or close to 5 pm/V, a typical value for bulk piezoelectric materials., Comment: 8 pages, 11 figures

The septuple-atomic-layer $\mathrm{VSi_2P_4}$ with the same structure of experimentally synthesized $\mathrm{MoSi_2N_4}$ is predicted to be a spin-gapless semiconductor (SGS). In this work, the biaxial strain is applied to tune electronic properties of $\mathrm{VSi_2P_4}$, and it spans a wide range of properties upon the increasing strain from ferromagnetic metal (FMM) to SGS to ferromagnetic semiconductor (FMS) to SGS to ferromagnetic half-metal (FMHM). Due to broken inversion symmetry, the coexistence of ferromagnetism and piezoelectricity can be achieved in FMS $\mathrm{VSi_2P_4}$ with strain range of 0\% to 4\%. The calculated piezoelectric strain coefficients $d_{11}$ for 1\%, 2\% and 3\% strains are 4.61 pm/V, 4.94 pm/V and 5.27 pm/V, respectively, which are greater than or close to a typical value of 5 pm/V for bulk piezoelectric materials. Finally, similar to $\mathrm{VSi_2P_4}$, the coexistence of piezoelectricity and ferromagnetism can be realized by strain in the $\mathrm{VSi_2N_4}$ monolayer. Our works show that $\mathrm{VSi_2P_4}$ in FMS phase with intrinsic piezoelectric properties can have potential applications in spin electronic devices., Comment: 6 pages,7 figures

Motivated by the fact that septuple-atomic-layer MnBi$_2$Te$_4$ can be structurally viewed as the combination of double-atomic-layer MnTe intercalating into quintuple-atomic-layer Bi$_2$Te$_3$, we present a general approach of constructing twelve septuple-atomic-layer $\alpha_i$- and $\beta_i$-$MA_2Z_4$ monolayer family (\emph{i} = 1 to 6) by intercalating MoS$_2$-type $MZ$$_2$ monolayer into InSe-type A$_2$Z$_2$ monolayer. Besides reproducing the experimentally synthesized $\alpha_1$-MoSi$_2$N$_4$, $\alpha_1$-WSi$_2$N$_4$ and $\beta_5$-MnBi$_2$Te$_4$ monolayer materials, another 66 thermodynamically and dynamically stable $MA_2Z_4$ were predicted, which span a wide range of properties upon the number of valence electrons (VEC). $MA_2Z_4$ with the rules of 32 or 34 VEC are mostly semiconductors with direct or indirect band gap and, however, with 33 VEC are generally metal, half-metal ferromagnetism, or spin-gapless semiconductor upon whether or not an unpaired electron is spin polarized. Moreover, we propose $\alpha_2$-WSi$_2$P$_4$ for the spin-valley polarization, $\alpha_1$-TaSi$_2$N$_4$ for Ising superconductor and $\beta_2$-SrGa$_2$Se$_4$ for topological insulator., Comment: Maintext 9 pages; 5 figures; Supplementary Materials 8 figures and 4 tables

The discovery of topological quantum states marks a new chapter in both condensed matter physics and materials sciences. By analogy to spin electronic system, topological concepts have been extended into phonons, boosting the birth of topological phononics (TPs). Here, we present a high-throughput screening and data-driven approach to compute and evaluate TPs among over 10,000 materials. We have clarified 5014 TP materials and classified them into single Weyl, high degenerate Weyl, and nodal-line (ring) TPs. Among them, three representative cases of TPs have been discussed in detail. Furthermore, we suggest 322 TP materials with potential clean nontrivial surface states, which are favorable for experimental characterizations. This work significantly increases the current library of TP materials, which enables an in-depth investigation of their structure-property relations and opens new avenues for future device design related to TPs., Comment: Maintext: 7 pages and 4 figures; Supplementary Materials: 14,246 pages, 5100 supplementary figures and 11,782 supplementary tables

Recent experimental realizations of the topological semimetal states in several monolayer systems are very attractive because of their exotic quantum phenomena and technological applications. Based on first-principles density-functional theory calculations including spin-orbit coupling, we here explore the drastically different two-dimensional (2D) topological semimetal states in three monolayers Cu$_2$Ge, Fe$_2$Ge, and Fe$_2$Sn, which are isostructural with a combination of the honeycomb Cu or Fe lattice and the triangular Ge or Sn lattice. We find that (i) the nonmagnetic (NM) Cu$_{2}$Ge monolayer having a planar geometry exhibits the massive Dirac nodal lines, (ii) the ferromagentic (FM) Fe$_2$Ge monolayer having a buckled geometry exhibits the massive Weyl points, and (iii) the FM Fe$_2$Sn monolayer having a planar geometry and an out-of-plane magnetic easy axis exhibits the massless Weyl nodal lines. It is therefore revealed that mirror symmetry cannot protect the four-fold degenerate Dirac nodal lines in the NM Cu$_{2}$Ge monolayer, but preserves the doubly degenerate Weyl nodal lines in the FM Fe$_{2}$Sn monolayer. Our findings demonstrate that the interplay of crystal symmetry, magnetic easy axis, and band topology is of importance for tailoring various 2D topological states in Cu$_2$Ge, Fe$_2$Ge, and Fe$_2$Sn monlayers.

Beryllium was recently discovered to harbor a Dirac nodal line (DNL) in its bulk phase and the DNL-induced non-trivial drumhead-like surface states (DNSSs) on its (0001) surface, rationalizing several already-existing historic puzzles [Phys. Rev. Lett., \textbf{117}, 096401 (2016)]. However, to date the underlying mechanism, as to why its (0001) surface exhibits an anomalously large electron-phonon coupling effect ($\lambda_{e-ph}^s$ $\approx$ 1.0), remains unresolved. Here, by means of first-principles calculations we have evidenced that the coupling of the DNSSs with the phononic states mainly contributes to its novel surface \emph{e-ph} enhancement. Besides that the experimentally observed $\lambda_{e-ph}^s$ and the main Eliashberg coupling function (ECF) peaks have been reproduced well, we have decomposed the ECF, $\alpha^{2}$$F$(\emph{k},\textbf{\emph{q}};\emph{v}), and the \emph{e-ph} coupling strength $\lambda(\emph{k},\textbf{\emph{q}};\emph{v})$ as a function of each electron momentum (\emph{k}), each phonon momentum (\textbf{\emph{q}}) and each phonon mode ($v$), evidencing the robust connection between the DNSSs and both $\alpha^{2}$$F$(\emph{k},\textbf{\emph{q}};\emph{v}) and $\lambda(\emph{k},\textbf{\emph{q}};\emph{v})$. The results reveal the strong \emph{e-ph} coupling between the DNSSs and the phonon modes, which contributes over 80$\%$ of the $\lambda_{e-ph}^s$ coefficient on the Be (0001) surface. It highlights that the anomalously large \emph{e-ph} coefficient on the Be (0001) surface can be attributed to the presence of its DNL-induced DNSSs, clarifying the long-term debated mechanism., Comment: 6 pages, 3 figures

By means of first-principles calculations and modeling analysis, we have predicted that the traditional 2D-graphene hosts the topological phononic Weyl-like points (PWs) and phononic nodal line (PNL) in its phonon spectrum. The phonon dispersion of graphene hosts three type-I PWs (both PW1 and PW2 at the BZ corners \emph{K} and \emph{K}', and PW3 locating along the $\Gamma$-\emph{K} line), one type-II PW4 locating along the $\Gamma$-\emph{M} line, and one PNL surrounding the centered $\Gamma$ point in the $q_{x,y}$ plane. The calculations further reveal that Berry curvatures are vanishingly zero throughout the whole BZ, except for the positions of these four pairs of Weyl-like phonons, at which the non-zero singular Berry curvatures appear with the Berry phase of $\pi$ or -$\pi$, confirming its topological non-trivial nature. The topologically protected non-trivial phononic edge states have been also evidenced along both the zigzag-edged and armchair-edged boundaries. These results would pave the ways for further studies of topological phononic properties of graphene, such as phononic destructive interference with a suppression of backscattering and intrinsic phononic quantum Hall-like effects., Comment: 5 pages, 4 figures

Electrides, with their excess electrons distributed in crystal cavities playing the role of anions, exhibit a variety of unique electronic and magnetic properties. In this work, we employ the first-principles crystal structure prediction to identify a new prototype of A$_3$B electride in which both interlayer spacings and intralayer vacancies provide channels to accommodate the excess electrons in the crystal. This A$_3$B type of structure is calculated to be thermodynamically stable for two alkaline metals oxides (Rb$_3$O and K$_3$O). Remarkably, the unique feature of multiple types of cavities makes the spatial arrangement of anionic electrons highly flexible via elastic strain engineering and chemical substitution, in contrast to the previously reported electrides characterized by a single topology of interstitial electrons. More importantly, our first-principles calculations reveal that Rb$_3$O is a topological Dirac nodal line semimetal, which is induced by the Rb-$s$ $\rightarrow$ O-$p$ band inversion at the general electronic k momentums in the Brillouin zone associated with the intersitial electric charges. The discovery of flexible electride in combining with topological electronic properties opens an avenue for electride design and shows great promises in electronic device applications., Comment: 15 pages, 12 figures

We study the optical properties of the Ruddlesden-Popper series of iridates Sr$_{n+1}$Ir$_n$O$_{3n+1}$ ($n$=1, 2 and $\infty$) by solving the Bethe-Salpeter equation (BSE), where the quasiparticle (QP) energies and screened interactions $W$ are obtained by the $GW$ approximation including spin-orbit coupling. The computed optical conductivity spectra show strong excitonic effects and reproduce very well the experimentally observed double-peak structure, in particular for the spin-orbital Mott insulators Sr$_2$IrO$_4$ and Sr$_3$Ir$_2$O$_7$. However, $GW$ does not account well for the correlated metallic state of SrIrO$_3$ owing to a much too small band renormalization, and this affects the overall quality of the optical conductivity. Our analysis describes well the progressive redshift of the main optical peaks as a function of dimensionality ($n$), which is correlated with the gradual decrease of the electronic correlation (quantified by the constrained random phase approximation) towards the metallic $n=\infty$ limit. We have also assessed the quality of a computationally cheaper BSE approach that is based on a model dielectric function and conducted on top of DFT+$U$ one-electron energies. Unfortunately, this model BSE approach does not accurately reproduce the outcome of the full $GW$+BSE method and leads to larger deviations to the measured spectra., Comment: 13 pages, 8 figures

Based on first-principles calculations, we predict that the superconducting MgB$_2$ with a AlB$_2$-type centrosymmetric lattice host the so-called phononic topological Weyl nodal lines (PTWNLs) on its bulk phonon spectrum. These PTWNLs can be viewed as countless Weyl points (WPs) closely aligned along the straight lines in the $-$H-K-H direction within the three-dimensional Brillouin zone (BZ). Their topological non-trivial natures are confirmed by the calculated Berry curvature distributions on the planes perpendicular to these lines. These lines are highly unique, because they exactly locate at the high-symmetry boundary of the BZ protected by the mirror symmetry and, simultaneously, straightly transverse the whole BZ, in different from known classifications including nodal rings, nodal chains or nets, and nodal loops. On the (10$\bar{1}$0) crystal surface, the PTWNLs-induced drumhead-like non-trivial surface states appear within the rectangular area confined by the projected lines of the PTWNLs with opposite chirality. Moreover, when the mirror symmetry is broken, the double-degenerate PTWNLs are further lifted to form a pair WPs with opposite chirality. Our results pave the ways for future experimental study on topological phonons on MgB$_2$ and highlights similar results in a series of isostructural AlB$_2$-type metallic diborides., Comment: 4 figures

In analogy to various fermions of electrons in topological semimetals, topological mechanical states with two type of bosons, Dirac and Weyl bosons, were reported in some macroscopic systems of kHz frequency and those with a type of doubly-Weyl phonons in atomic vibrational framework of THz frequency of solid crystals were recently predicted. However, to date no three-component bosons of phonon has been reported. Here, through first-principles calculations, we have reported that the phonon spectra of the WC-type TiS, ZrSe, and HfTe commonly host the unique triply degenerate nodal points (TDNPs) and single two-component Weyl points (WPs) in THz frequency. Quasiparticle excitations near TDNPs of phonons are three-component bosons, beyond the conventional and known classifications of Dirac, Weyl and doubly-Weyl phonons. {Moreover, we have found that both TiS and ZrSe have five pairs of type-I Weyl phonons and a pair of type-II Weyl phonons, whereas HfTe only has four pairs of type-I Weyl phonons. They carry non-zero topological charges. On the (10$\bar{1}$0) crystal surfaces, we observe topological protected surface arc states connecting two WPs with opposite charges, which host modes that propagate nearly in one direction on the surface., Comment: 10 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:1706.06876

Topological semimetals with several types of three-dimensional (3D) fermion of electrons, such as Dirac fermions, Weyl fermions, Dirac nodal lines and triply degenerate nodal points have been theoretically predicted and then experimentally discovered in the electronic structures of a series of solid crystals. In analogy of various typical fermions, topological mechanical states with two type of bosons, Dirac and Weyl bosons, were also experimentally reported in some macroscopic systems of kHz frequency and with a type of doubly-Weyl phonons in atomic vibrational framework of THz frequency of solid crystal was also recently predicted. However, to date no triply degenerate nodal point of phonon beyond the conventional Dirac, Weyl and doubly-Weyl phonons has been reported. Here, through first-principles calculations, we have reported on the prediction that the WC-type TiS, ZrSe, and HfTe commonly host the unique triply degenerate nodal point of phonon in THz frequency due to the occurrence of the phonon band inversion between the doubly degenerate planar vibrational mode and the singlet vertical vibrational mode at the boundary A point of the bulk Brillouin zone. Quasiparticle excitations near this triply degenerate nodal point of phonons are three-component bosons, different from the known classifications. The underlying mechanism can be attributed to the leading role of the comparable atomic masses of constituent elements in compounds in competition with the interatomic interaction. Additionally, the electronic structures in their bulk crystals exhibit the coexisted triply degenerate nodal point and Weyl fermions. The novel coexistence of three-component bosons, three-component fermions and Weyl fermions in these materials thus suggest an enriched platform for studying the interplay between different types of fermions and bosons., Comment: 9 figures, 1 table, 10 pages

Topological materials in crystal solids, including topological insulators (TIs), topological crystalline insulators (TCIs), topological Dirac semimetals (DSMs), topological Weyl semimetals (WSMs), topological Dirac or Weyl nodal line semimetals (NLSMs) and beyond, are mainly featured with topological, protected non-trivial surface states, and their bulk phases are insulators or semimetals with the proper presence of Dirac cones, Weyl nodes or Dirac nodal lines around the Fermi level. The author suggests a two-step computational screening strategy of 3D topological materials by mixing chemistry with physics with the considerations of fully filled bands and band inversion., Comment: 3 pages, 1 figure, accepted by National Science Review

Topological nodal line (DNL) semimetals, formed by a closed loop of the inverted bands in the bulk, result in the nearly flat drumhead-like surface states with a high electronic density near the Fermi level. The high catalytic active sites associated with the high electronic densities, the good carrier mobility, and the proper thermodynamic stabilities with $\Delta G_{H^*}$$\approx$0 are currently the prerequisites to seek the alternative candidates to precious platinum for catalyzing electrochemical hydrogen (HER) production from water. Within this context, it is natural to consider whether or not the DNLs are a good candidate for the HER because its non-trivial surface states provide a robust platform to activate possibly chemical reactions. Here, through first-principles calculations we reported on a new DNL TiSi-type family with a closed Dirac nodal line consisting of the linear band crossings in the $k_y$ = 0 plane. The hydrogen adsorption on the (010) and (110) surfaces yields the $\Delta G_{H^*}$ to be almost zero. The topological charge carries have been revealed to participate in this HER. The results are highlighting that TiSi not only is a promising catalyst for the HER but also paves a new routine to design topological quantum catalyst utilizing the topological DNL-induced surface bands as active sites, rather than edge sites-, vacancy-, dopant-, strain-, or heterostructure-created active sites., Comment: 8 pages, 5 figures

The recent discoveries of ferroelectric metal and Weyl semimetal (WSM) have stimulated a natural question: whether these two exotic states of matter can coexist in a single material or not. These two discoveries ensure us that physically it is possible since both of them share the same necessary condition, the broken inversion symmetry. Here, by using first-principles calculations, we demonstrate that the experimentally synthesized nonmagnetic HgPbO$_3$ represents a unique example of such hybrid "\emph{Weyl ferroelectric semimetal}". Its centrosymmetric $R\bar{3}c$ phase will undergo a ferroelectric phase transition to the ferroelectric $R3c$ structure. Both phases are metallic and the ferroelectric phase owns a spontaneous polarization of 33 $\mu$C/cm$^2$. Most importantly, it also harbors six pairs of chiral Weyl nodes around the Fermi level to be an oxide WSM. The structural symmetry broken phase transition induces a topological phase transition. The coexistence of ferroelectricity and Weyl nodes in HgPbO$_3$ is an ideal platform for exploring multiphase interaction and mutual control. The Weyl nodes can be tuned by external pulse electric field, which is promising for potential applications of integrated topotronic and ferroelectric devices., Comment: 6 pages, 6 figures

We study the effects of dilute La and Rh doping on the electronic structure of the relativistic Mott insulator Sr$_2$IrO$_4$ using fully relativistic and magnetically non-collinear density functional theory with the inclusion of an on-site Hubbard $U$. To model doping effects, we have adopted the supercell approach that allows for a realistic treatment of structural relaxations and electronic effects beyond a purely rigid band approach. By means of the band unfolding technique we have computed the spectral function and constructed the effective band structure and Fermi surface (FS), which are readily comparable with available experimental data. Our calculations clearly indicate that La and Rh doping can be interpreted as effective electron and (fractional) hole doping, respectively. In Sr$_{2-x}$La$_x$IrO$_4$ the IMT is accompanied by a moderate renormalization of the electronic correlation substantiated by a reduction of the effective on-site Coulomb repulsion $U-J$ from 1.6 eV ($x=0$) to 1.4 eV (metallic regime $x=12.5\%$). The progressive closing of the relativistic Mott gap leads to the emergence of connected elliptical electron pockets at ($\pi$/2,$\pi$/2) and less intense features at $X$ in the FS. The substitution of Ir with the nominally isovalent Rh is accompanied by a substantial hole transfer from the Rh site to the nearest neighbor Ir sites. This shifts down the chemical potential, creates almost circular disconnected hole pockets in the FS and establishes the emergence of a two-dimensional metallic state formed by conducting Rh-planes intercalated by insulating Ir-planes. Finally, our data indicate that hole doping causes a flipping of the in-plane net ferromagnetic moment in the Rh plane and induces a magnetic transition from the AF-I to the AF-II ordering.

Beryllium is a simple alkali earth metal, but has been the target of intensive studies for decades because of its unusual electron behaviors at surfaces. Puzzling aspects include (i) severe deviations from the description of the nearly free electron picture, (ii) anomalously large electron-phonon coupling effect, and (iii) giant Friedal oscillations. The underlying origins for such anomalous surface electron behaviors have been under active debate, but with no consensus. Here, by means of first-principle calculations, we discover that this pure metal system, surprisingly, harbors the Dirac node line (DNL) that in turn helps to rationalize many of the existing puzzles. The DNL is featured by a closed line consisting of linear band crossings and its induced topological surface band agrees well with previous photoemission spectroscopy observation on Be (0001) surface. We further reveal that each of the elemental alakali earth metals of Mg, Ca, and Sr also harbors the DNL, and speculate that the fascinating topological property of DNL might naturally exist in other elemental metals as well., Comment: 5 pages, 4 figures

The Mg-Si-O system is the major Earth and rocky planet-forming system. Here, through quantum variable-composition evolutionary structure explorations, we have discovered several unexpected stable binary and ternary compounds in the Mg-Si-O system. Besides the well-known SiO2 phases, we have found two extraordinary silicon oxides, SiO3 and SiO, which become stable at pressures above 0.51 TPa and 1.89 TPa, respectively. In the Mg-O system, we have found one new compound, MgO3, which becomes stable at 0.89 TPa. We find that not only the (MgO)x(SiO2)y compounds, but also two (MgO3)x(SiO3)y compounds, MgSi3O12 and MgSiO6, have stability fields above 2.41 TPa and 2.95 TPa, respectively. The highly oxidized MgSi3O12 can form in deep mantles of mega-Earths with masses above 20 M+ (M+:Earth's mass). Furthermore, the dissociation pathways of pPv-MgSiO3 are also clarified, and found to be different at low and high temperatures. The low-temperature pathway is MgSiO3 -> Mg2SiO4 + MgSi2O5 -> SiO2 + Mg2SiO4 -> MgO + SiO2, while the high-temperature pathway is MgSiO3 -> Mg2SiO4 + MgSi2O5 -> MgO + MgSi2O5 -> MgO + SiO2. Present results are relevant for models of the internal structure of giant exoplanets, and for understanding the high-pressure behavior of materials., Comment: 24 pages, 6 figures

By means of first-principles calculations, we identified the structural phase transition of Na$_3$Bi from hexagonal ground state to cubic $cF$16 phase above 0.8 GPa, in agreement with the experimental findings. Upon the releasing of pressure, \emph{cF}16 phase of Na$_3$Bi is mechanically stable at ambient condition. The calculations revealed that the $cF$16 phase is topological semimetal, in similarity to well-known HgTe and it even exhibits an unusually low $C^\prime$ modulus (only about 1.9 GPa) and a huge anisotropy, $A^u$ of as high as 11, the third highest value among all known cubic crystals in their elastic behaviors. These facts render \emph{cF}16-type Na$_3$Bi very soft with a liquid-like elastic deformation in the (110)$<$1$\overline{1}$0$>$ slip system. Importantly, as accompanied with this deformation, Na$_3$Bi shows a topological phase transition from a topological semimetal state at its strain-free cubic phase to a topological insulating state at its distorted phase. Because the $C^\prime$ elastic deformation almost costs no energy in a reversible and liquid-like soft manner, \emph{cF}16-type Na$_3$Bi would potentially provide a fast on/off switching way between topological insulator and topological semimetal, which would be beneficial to the quantum electronic devices for practical applications., Comment: 5 pages, 4 figures

We put forward a scheme to study the anisotropic magnetic couplings in Sr2IrO4 by mapping fully relativistic constrained noncollinear density functional theory including an on-site Hubbard U correction onto a general spin model Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the interplay between structural distortions and the competition between isotropic exchange and DM interactions. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation.

By means of first-principles and \emph{ab initio} tight-binding calculations, we found that the compound of NaBi is a three-dimensional non-trivial topological metal. Its topological feature can be confirmed by the presence of band inversion, the derived effective Z$_2$ invariant and the non-trivial surface states with the presence of Dirac cones. Interestingly, our calculations further demonstrated that NaBi exhibits the uniquely combined properties between the electron-phonon coupling superconductivity in nice agreement with recent experimental measurements and the obviously anisotropic but extremely low thermal conductivity. The spin-orbit coupling effects greatly affect those properties. NaBi may provide a rich platform to study the relationship among metal, topology, superconductivity and thermal conductivity., Comment: 5 pages, 4 figures

By means of first-principles calculations coupled with the kinetic Monte Carlo simulations, we have systematically investigated the effects of dilute substitutional solutes on the behaviors of carbon in $\alpha$-Fe. Our results uncover that: ($i$) Without the Fe vacancy the interactions between most solutes and carbon are repulsive due to the strain relief, whereas Mn has a weak attractive interaction with its nearest-neighbor carbon due to the local ferromagnetic coupling effect. ($ii$) The presence of the Fe vacancy results in attractive interactions of all the solutes with carbon. In particular, the Mn-vacancy pair shows an exceptionally large binding energy of -0.81 eV with carbon. ($iii$) The alloying addition significantly impacts the atomic-scale concentration distributions and chemical potential of carbon in the Fe matrix. Among them, Mn and Cr increase the carbon chemical potential whereas Al and Si reduce it. ($iv$) Within the dilute scale of the alloying solution, the solute concentration and temperature dependent carbon diffusivities demonstrate that Mn has a little impact on the carbon diffusion whereas Cr (Al or Si) remarkably retards the carbon diffusion. Our results provide certain implication for better understanding the experimental observations related with the carbon solubility limit, carbon micro-segregation and carbide precipitations in the ferritic steels., Comment: 13 pages, 14 figure, Phys. Rev. B, accepted

By means of first-principles calculations, we found that the early characterized three-dimensional topological bulk Dirac semimetal of the $P$6$_3$/$mmc$ Na$_3$Bi is dynamically unstable at the ground state due to the presence of the large imaginary phonon frequencies around the $K$ point. Alternatively, our calculations suggest a new ground state phase which crystalizes in the $P$$\overline3$$c$1 structure with a buckled graphite-like Na/Bi sheets, which is both energetically and dynamically more stable. Moreover, the calculations also uncovered that the $P$$\overline3$$c$1 phase at the ground state is also a topological 3D Dirac semimetal, being exactly the same as the electronic states of the metastable $P$6$_3$/$mmc$ phase., Comment: 5 pages, 4 figures

Interaction between hydrogen (H) and metals is central to many materials problems of scientific and technological importance. Chief among them is the development of H storage and H-resistant materials. H segregation or trapping at lattice defects, including vacancies, dislocations, grain boundaries, etc, plays a crucial role in determining the properties of these materials. Here, through first-principles simulations, we propose a unified mechanism involving charge transfer induced strain destabilization to understand H segregation behavior at vacancies. We discover that H prefers to occupy interstitials with high pre-existing charge densities and the availability of such interstitials sets the limit on H trapping capacity at a vacancy. Once the maximum H capacity is reached, the dominant charge donors switch from the nearest-neighbor (NN) to the next-nearest-neighbor (NNN) metal atoms. Accompanying with this long-range charge transfer, a significant reorganization energy would occur, leading to instability of the H-vacancy complex. The physical picture unveiled here appears universal across the BCC series and is believed to be relevant to other metals/defects as well. The insight gained from this study is expected to have important implications for the design of H storage and H-resistant materials., Comment: 5 pages, 3 figures

By means of variable-composition evolutionary algorithm coupled with density functional theory and in combination with aberration-corrected high-resolution transmission electron microscopy experiments, we have studied and characterized the composition, structure and hardness properties of WB$_{3+x}$ ($x$ $<$ 0.5). We provide robust evidence for the occurrence of stoichiometric WB$_3$ and non-stoichiometric WB$_{3+x}$ both crystallizing in the metastable $hP$16 ($P6_3/mmc$) structure. No signs for the formation of the highly debated WB$_4$ (both $hP$20 and $hP$10) phases were found. Our results rationalize the seemingly contradictory high-pressure experimental findings and suggest that the interstitial boron atom is located in the tungsten layer and vertically interconnect with four boron atoms, thus forming a typical three-center boron net with the upper and lower boron layers in a three-dimensional covalent network, which thereby strengthen the hardness., Comment: 5 pages, 5 figures

Unlike time-reversal topological insulators, surface metallic states with Dirac cone dispersion in the recently discovered topological crystalline insulators (TCIs) are protected by crystal symmetry. To date, TCI behaviors have been observed in SnTe and the related alloys Pb$_{1-x}$Sn$_{x}$Se/Te, which incorporate heavy elements with large spin-orbit coupling (SOC). Here, by combining first-principles and {\it ab initio} tight-binding calculations, we report the formation of a TCI in the relatively lighter rock-salt SnS and SnSe. This TCI is characterized by an even number of Dirac cones at the high-symmetry (001), (110) and (111) surfaces, which are protected by the reflection symmetry with respect to the ($\bar{1}$10) mirror plane. We find that both SnS and SnSe have an intrinsically inverted band structure and the SOC is necessary only to open the bulk band gap. The bulk band gap evolution upon volume expansion reveals a topological transition from an ambient pressure TCI to a topologically trivial insulator. Our results indicate that the SOC alone is not sufficient to drive the topological transition., Comment: 5 pages, 5 figures

The phenomenon, channel segregates (CS) as a result of gravity-driven flow due to density contrast occurred in the solid-liquid mushy zones1during solidification, often causes the severe destruction of homogeneity and even some fatal damages. Investigation on its mechanism sheds light on the understanding and controlling of the formation of solidifying metals,earth's core, igneous rock and sea ice. Until now, it still remains controversial what composes the density contrasts and, to what extent, how it affects channel segregates. Here, we show that in experimental 500kg and 100 ton commercial cast steel ingots CS originates from oxide Al2O3/MnS impurity clusters (OICs) initially nucleated from the oxide (Al2O3) particles, which induce an extra flow due to sharp density contrast between clusters and melt. The results uncover that, as OICs enrich and grow, their driven flow becomes stronger than the traditionally recognized inter-dendritic thermo-solutal convection, dominating the subsequent opening of the channels. This study extends the classical macrosegregation theory, highlights a significant technological breakthrough to control CS, and could quickly yield practical benefits to the worldwide manufacture of over 50 million tons of ingots, super-thick slab and heavy castings annually, as well as has general implications for the elaboration of other related natural phenomena., Comment: 18 pages, 5 figures

In integrating topological insulators (TIs) with conventional materials, one crucial issue is how the topological surface states (TSS) will behave in such heterostructures. We use first-principles approaches to establish accurate tunability of the vertical location of the TSS via intriguing dual-proximity effects. By depositing a conventional insulator (CI) overlayer onto a TI substrate (Bi2Se3 or Bi2Te3), we demonstrate that, the TSS can float to the top of the CI film, or stay put at the CI/TI interface, or be pushed down deeper into the otherwise structurally homogeneous TI substrate. These contrasting behaviors imply a rich variety of possible quantum phase transitions in the hybrid systems, dictated by key material-specific properties of the CI. These discoveries lay the foundation for accurate manipulation of the real space properties of TSS in TI heterostructures of diverse technological significance.

By electron and X-ray diffraction we establish that the CrB$_4$ compound discovered over 40 years ago crystallizes in the $oP10$ (\emph{Pnnm}) structure, in disagreement with previous experiments but in agreement with a recent first-principles prediction. The 3D boron network in the new structure is a distorted version of the rigid carbon $sp^3$ network proposed recently for the high-pressure C$_4$ allotrope. According to our density functional theory calculations and the analysis of the bonding, CrB$_4$ is a potentially superhard material. In fact, the calculated weakest shear and tensile stresses exceed 50 GPa and its Vickers hardness is estimated to be 48 GPa., Comment: 5 pages, 4 figures

We report a general scheme to systematically construct two classes of structural families of superhard sp3 carbon allotropes of cold compressed graphite through the topological analysis of odd 5+7 or even 4+8 membered carbon rings stemmed from the stacking of zigzag and armchair chains. Our results show that the previously proposed M, bct-C4, W and Z allotropes belong to our currently proposed families and that depending on the topological arrangement of the native carbon rings numerous other members are found that can help us understand the structural phase transformation of cold-compressed graphite and carbon nanotubes (CNTs). In particular, we predict the existence of two simple allotropes, R- and P-carbon, which match well the experimental X-ray diffraction patterns of cold-compressed graphite and CNTs, respectively, display a transparent wide-gap insulator ground state and possess a large Vickers hardness comparable to diamond., Comment: 5 pages, 4 figures, accepted by Phys. Rev. Lett

We revisit and interpret the mechanical properties of the recently proposed allotrope of carbon, T-carbon [Sheng \emph{et al.}, Phys. Rev. Lett., \textbf{106}, 155703 (2011)], using density functional theory in combination with different empirical hardness models. In contrast with the early estimation based on the Gao's model, which attributes to T-carbon an high Vickers hardness of 61 GPa comparable to that of superhard cubic boron nitride (\emph{c}-BN), we find that T-carbon is not a superhard material, since its Vickers hardenss does not exceed 10 GPa. Besides providing clear evidence for the absence of superhardenss in T-carbon, we discuss the physical reasons behind the failure of Gao's and \v{S}im$\rm\mathring{u}$nek and Vack\'a\v{r}'s (SV) models in predicting the hardness of T-carbon, residing on their improper treatment of the highly anisotropic distribution of quasi-\emph{sp}$^3$-like C-C hybrids. A possible remedy to the Gao and SV models based on the concept of superatom is suggest, which indeed yields a Vickers hardness of about 8 GPa., Comment: 5 pages and 3 figures, accepted by Phys. Rev. B as a Rapid Communication

We explore the topological behavior of the binary Zintl phase of the alkaline earth metals based compounds Sr$_2$Pb and Sr$_2$Sn using both standard and hybrid density functional theory. It is found that Sr$_2$Pb lies on the verge of a topological instability which can be suitably tuned through the application of a small uniaxial expansion strain ($>$ 3%). The resulting non-trivial topologically insulating state display well-defined metallic states in the Sr$_{2}$Pb(010) surface, whose evolution is studied as a function of the film thickness., Comment: 5 pages, 4 figures