Source author record

Peitao Liu

Peitao Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Two-dimensional Obstructed Atomic Insulators with Fractional Corner Charge in MA$_2$Z$_4$ Family

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.

preprint2021arXiv

Phase Transitions of Zirconia: Machine-Learned Force Fields Beyond Density Functional Theory

We present an approach to generate machine-learned force fields (MLFF) with beyond density functional theory (DFT) accuracy. Our approach combines on-the-fly active learning and $Δ$-machine learning in order to generate an MLFF for zirconia based on the random phase approximation (RPA). Specifically, an MLFF trained on-the-fly during DFT based molecular dynamics simulations is corrected by another MLFF that is trained on the differences between RPA and DFT calculated energies, forces and stress tensors. Thanks to the relatively smooth nature of the differences, the expensive RPA calculations are performed only on a small number of representative structures of small unit cells. These structures are determined by a singular value decomposition rank compression of the kernel matrix with low spatial resolution. This dramatically reduces the computational cost and allows us to generate an MLFF fully capable of reproducing high-level quantum-mechanical calculations beyond DFT. We carefully validate our approach and demonstrate its success in studying the phase transitions of zirconia.

preprint2020arXiv

Accurate optical spectra through time-dependent density functional theory based on screening-dependent hybrid functionals

We investigate optical absorption spectra obtained through time-dependent density functional theory (TD-DFT) based on nonempirical hybrid functionals that are designed to correctly reproduce the dielectric function. The comparison with state-of-the-art $GW$ calculations followed by the solution of the Bethe-Sapeter equation (BSE-$GW$) shows close agreement for both the transition energies and the main features of the spectra. We confront TD-DFT with BSE-$GW$ by focusing on the model dielectric function and the local exchange-correlation kernel. The present TD-DFT approach achieves the accuracy of BSE-$GW$ at a fraction of the computational cost.

preprint2020arXiv

Comparative ab initio study of the structural, electronic, magnetic, and dynamical properties of LiOsO$_3$ and NaOsO$_3$

Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means of comparative first-principles DFT+$U$ calculations with the inclusion of the spin-orbit coupling, we ($i$) identify the origin of the different structural ($R3c$ vs. $Pnma$) properties using a symmetry-adapted soft mode analysis, ($ii$) provide evidence that all considered exchange-correlation functionals (LDA, PBE, PBEsol, SCAN, and HSE06) and the spin disordered polymorphous descriptions are unsatisfactory to accurately describe the electronic and magnetic properties of both systems simultaneously, and ($iii$) clarify that the distinct electronic (metallic vs. insulating) properties originates mainly from a cooperative steric and magnetic effect. Finally, we find that although at ambient pressure LiOsO$_3$ with a $Pnma$ symmetry and NaOsO$_3$ with a $R\bar{3}c$ symmetry are energetically unfavorable, they do not show soft phonons and therefore are dynamically stable. A pressure-induced structural phase transition from $R3c$ to $Pnma$ for LiOsO$_3$ is predicted, whereas for NaOsO$_3$ no symmetry change is discerned in the considered pressure range.

preprint2019arXiv

Assessing model-dielectric-dependent hybrid functionals on the antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO

Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common model dielectric function approach, but differ in the way how to non-empirically obtain the range-separation parameter. By retaining the full short-range Fock exchange and a fraction of the long-range Fock exchange that equals the inverse of the dielectric constant, both DD-RSH-CAM and DSH turn out to perform very well in predicting the band gaps for a large variety of semiconductors and insulators. Here, we assess how these two hybrid functionals perform on challenging antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO by comparing them to other conventional hybrid functionals and the $GW$ method. We find that single-shot DD0-RSH-CAM and DSH0 improve the band gaps towards experiments as compared to conventional hybrid functionals. The magnetic moments are slightly increased, but the predicted dielectric constants are decreased. The valence band density of states (DOS) predicted by DD0-RSH-CAM and DSH0 are as satisfactory as HSE03 in comparison to experimental spectra, however, the conduction band DOS are shifted to higher energies by about 2 eV compared to HSE03. Self-consistent DD-RSH-CAM and DSH deteriorate the results with a significant overestimation of band gaps.

preprint2016arXiv

$GW$100: a plane wave perspective for small molecules

In a recent work, van Setten and coworkers have presented a carefully converged $G_0W_0$ study of 100 closed shell molecules [J. Chem. Theory Comput. 11, 5665 (2015)]. For two different codes they found excellent agreement to within few 10 meV if identical Gaussian basis sets were used. We inspect the same set of molecules using the projector augmented wave method and the Vienna ab initio simulation package (VASP). For the ionization potential, the basis set extrapolated plane wave results agree very well with the Gaussian basis sets, often reaching better than 50 meV agreement. In order to achieve this agreement, we correct for finite basis set errors as well as errors introduced by periodically repeated images. For electron affinities below the vacuum level differences between Gaussian basis sets and VASP are slightly larger. We attribute this to larger basis set extrapolation errors for the Gaussian basis sets. For quasi particle (QP) resonances above the vacuum level, differences between VASP and Gaussian basis sets are, however, found to be substantial. This is tentatively explained by insufficient basis set convergence of the Gaussian type orbital calculations as exemplified for selected test cases.

preprint2016arXiv

Cubic scaling $GW$: towards fast quasiparticle calculations

Within the framework of the full potential projector-augmented wave methodology, we present a promising low-scaling $GW$ implementation. It allows for quasiparticle calculations with a scaling that is cubic in the system size and linear in the number of $k$ points used to sample the Brillouin zone. This is achieved by calculating the polarizability and self-energy in the real space and imaginary time domain. The transformation from the imaginary time to the frequency domain is done by an efficient discrete Fourier transformation with only a few nonuniform grid points. Fast Fourier transformations are used to go from real space to reciprocal space and vice versa. The analytic continuation from the imaginary to the real frequency axis is performed by exploiting Thiele's reciprocal difference approach. Finally, the method is applied successfully to predict the quasiparticle energies and spectral functions of typical semiconductors (Si, GaAs, SiC, and ZnO), insulators (C, BN, MgO, and LiF), and metals (Cu and SrVO$_3$). The results are compared with conventional $GW$ calculations. Good agreement is achieved, highlighting the strength of the present method.

preprint2016arXiv

Electron and hole doping in the relativistic Mott insulator Sr$_2$IrO$_4$: a first-principles study using band unfolding technique

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 ($π$/2,$π$/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.

preprint2016arXiv

Lifshitz transition driven by spin fluctuations and spin-orbit renormalization in NaOsO$_3$

In systems where electrons form both dispersive bands and small local spins, we show that changes of the spin configuration can tune the bands through a Lifshitz transition, resulting in a continuous metal-insulator transition associated with a progressive change of the Fermi surface topology. In contrast to a Mott-Hubbard and Slater pictures, this spin-driven Lifshitz transition appears in systems with small electron-electron correlation and large hybridization. We show that this situation is realized in 5$d$ distorted perovskites with an half-filled $t_{2g}$ bands such as NaOsO$_3$, where the strong $p-d$ hybridization reduces the local moment, and spin-orbit coupling causes a large renormalization of the electronic mobility. This weakens the role of electronic correlations and drives the system towards an itinerant magnetic regime which enables spin-fluctuations.

preprint2015arXiv

Anisotropic Magnetic Couplings and Structure-Driven Canted to Collinear Transitions in Spin-orbit Coupled Sr2IrO4

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.

preprint2014arXiv

Effects of dilute substitutional solutes on carbon in $α$-Fe: interactions and associated carbon diffusion from first-principles calculations

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 $α$-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.

preprint2013arXiv

Interstitial-Boron Solution Strengthened WB$_{3+x}$

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.