BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hongjun Xiang

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallcond-mat.str-elcond-mat.supr-conphysics.comp-phMachine Learningcond-mat.other

Catalog footprint

What is connected

41works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

41 published item(s)

preprint2024arXiv

Origin of zigzag antiferromagnetic orders in XPS3 (X= Fe, Ni) monolayers

Recently, two monolayer magnetic materials, i.e., FePS3 and NiPS3, have been successfully fabricated. Despite that they have the same atomic structure, the two monolayers exhibit distinct magnetic properties. FePS3 holds an out-of-plane zigzag antiferromagnetic (AFM-ZZ) structure, while NiPS3 exhibits an in-plane AFM-ZZ structure. However, there is no theoretical model which can properly describe its magnetic ground state due to the lack of a full understanding of its magnetic interactions. Here, by combining the first-principles calculations and the newly developed machine learning method, we construct an exact spin Hamiltonian of the two magnetic materials. Different from the previous studies which failed to fully consider the spin-orbit coupling effect, we find that the AFM-ZZ ground state in FePS3 is stabilized by competing ferromagnetic nearest-neighbor and antiferromagnetic third nearest-neighbor exchange interactions, and combining single-ion anisotropy. Whereas, the often ignored nearest-neighbor biquadratic exchange is responsible for the in-plane AFM-ZZ ground state in NiPS3. We additionally calculate spin-wave spectrum of AFM-ZZ structure in the two monolayers based on the exact spin Hamiltonian, which can be directly verified by the experimental investigation. Our work provides a theoretical framework for the origin of AFM-ZZ ground state in two-dimensional materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Convert widespread paraelectric perovskite to ferroelectrics

While nature provides a plethora of perovskite materials, only a few exhibits large ferroelectricity and possibly multiferroicity. The majority of perovskite materials have the non-polar CaTiO$_3$(CTO)structure, limiting the scope of their applications. Based on effective Hamiltonian model as well as first-principles calculations, we propose a general thin-film design method to stabilize the functional BiFeO$_3$(BFO)-type structure, which is a common metastable structure in widespread CaTiO$_3$-type perovskite oxides. It is found that the improper antiferroelectricity in CTO-type perovskite and ferroelectricity in BFO-type perovskite have distinct dependences on mechanical and electric boundary conditions, both of which involve oxygen octahedral rotation and tilt. The above difference can be used to stabilize the highly polar BFO-type structure in many CTO-type perovskite materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Edge-based Tensor prediction via graph neural networks

Message-passing neural networks (MPNN) have shown extremely high efficiency and accuracy in predicting the physical properties of molecules and crystals, and are expected to become the next-generation material simulation tool after the density functional theory (DFT). However, there is currently a lack of a general MPNN framework for directly predicting the tensor properties of the crystals. In this work, a general framework for the prediction of tensor properties was proposed: the tensor property of a crystal can be decomposed into the average of the tensor contributions of all the atoms in the crystal, and the tensor contribution of each atom can be expanded as the sum of the tensor projections in the directions of the edges connecting the atoms. On this basis, the edge-based expansions of force vectors, Born effective charges (BECs), dielectric (DL) and piezoelectric (PZ) tensors were proposed. These expansions are rotationally equivariant, while the coefficients in these tensor expansions are rotationally invariant scalars which are similar to physical quantities such as formation energy and band gap. The advantage of this tensor prediction framework is that it does not require the network itself to be equivariant. Therefore, in this work, we directly designed the edge-based tensor prediction graph neural network (ETGNN) model on the basis of the invariant graph neural network to predict tensors. The validity and high precision of this tensor prediction framework were shown by the tests of ETGNN on the extended systems, random perturbed structures and JARVIS-DFT datasets. This tensor prediction framework is general for nearly all the GNNs and can achieve higher accuracy with more advanced GNNs in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flat-band based ferromagnetic semiconducting state in the graphitic C$_4$N$_3$ monolayer

A new set of lattice-models based on the hexagonal $\sqrt{N}\times\sqrt{N}$ super-cells of the well-known honeycomb lattice with single-hole defect (HL-D-1/2N) are proposed to realize the nontrivial isolated flat-bands. Through performing both tight-binding and density functional theory calculations, we demonstrate that the experimentally realized graphitic carbon nitride (Adv. Mater., 22, 1004, 2010; Nat. Commun., 9, 3366, 2018), the HL-D-1/8 based C$_4$N$_3$, is a perfect system to host such flat bands. For the flat high-energy P-6m2 C$_4$N$_3$ structure, it displays the ferromagnetic half-metallicity which is not related to the isolated flat bands. However, the P-6m2 C$_4$N$_3$ structure is dynamically unstable. Using a structure searching method based on group and graph theory, we find that a new corrugated Pca21 C4N3 structure has the lowest energy among all known C$_4$N$_3$ structures. This Pca21 C$_4$N$_3$ structure is an intrinsic ferromagnetic half-semiconductor (Tc$\approx$241 K) with one semiconducting spin-channel (1.75 eV) and one insulating spin-channel (3.64 eV), which is quite rare in the two-dimensional (2D) systems. Its ferromagnetic semiconducting property originates from the isolated p$_z$-state flat-band as the corrugation shift the flat band upward to the Fermi level. Interestingly, this Pca21 C$_4$N$_3$ structure is found to be piezoelectric and ferroelectric, which makes C$_4$N$_3$ an unusual transition-metal-free 2D multiferroic.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Light-induced Magnetic Phase Transition in van der Waals Antiferromagnets

Based on a simple tight-binding model, we propose a general theory of light-induced magnetic phase transition (MPT) in antiferromagnets based on the general conclusion that the bandgap of antiferromagnetic (AFM) phase is usually larger than that of ferromagnetic (FM) one in a given system. Light-induced electronic excitation prefers to stabilize the FM state over the AFM one, and once the critical photocarrier concentration (α_c) is reached, an MPT from AFM phase to FM phase takes place. This theory has been confirmed by performing first-principles calculations on a series of two-dimensional (2D) van der Waals (vdW) antiferromagnets and a linear relationship between α_c and the intrinsic material parameters is obtained. Importantly, our conclusion is still valid even considering the strong exciton effects during photoexcitation. Our general theory provides new ideas to realize reversible read-write operations for future memory devices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Two-Dimensional Organic-Inorganic Room-Temperature Multiferroic

Organic-inorganic multiferroics are promising for the next generation of electronic devices. To date, dozens of organic-inorganic multiferroics have been reported; however, most of them show magnetic Curie temperature much lower than room temperature, which drastically hampers their application. Here, by performing first-principle calculations and building effective model Hamiltonians, we reveal a molecular orbital-mediated magnetic coupling mechanism in two-dimensional Cr(pyz)2 (pyz=pyrazine), and the role that the valence state of the molecule plays in determining the magnetic coupling type between metal ions. Based on these, we demonstrate that a two-dimensional organic-inorganic room-temperature multiferroic, Cr(h-fpyz)2 (h-fpyz= half-fluoropyrazine), can be rationally designed by introducing ferroelectricity in Cr(pyz)2 while keeping the valence state of the molecule unchanged. Our work not only reveals the origin of magnetic coupling in 2D organic-inorganic systems, but also provides a way to design room temperature multiferroic materials rationally.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Complex Spin Hamiltonian Represented by Artificial Neural Network

The effective spin Hamiltonian method is widely adopted to simulate and understand the behavior of magnetism. However, the magnetic interactions of some systems, such as itinerant magnets, are too complex to be described by any explicit function, which prevents an accurate description of magnetism in such systems. Here, we put forward a machine learning (ML) approach, applying an artificial neural network (ANN) and a local spin descriptor to develop effective spin potentials for any form of interaction. The constructed Hamiltonians include an explicit Heisenberg part and an implicit non-linear ANN part. Such a method successfully reproduces artificially constructed models and also sufficiently describe the itinerant magnetism of bulk Fe3GeTe2. Our work paves a new way for investigating complex magnetic phenomena (e.g., skyrmions) of magnetic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Berezinskii-Kosterlitz-Thouless Phase in Two-dimensional Ferroelectrics

The celebrated Berezinskii-Kosterlitz-Thouless (BKT) phase transition refers to a topological transition characterized, e.g., by the dissociation of vortex-antivortex pairs in two-dimensional (2D) systems. Such unusual phase has been reported in various types of materials, but never in the new class of systems made by one-unit-cell-thick (1UC) ferroelectrics (also coined as 2D ferroelectrics). Here, the use of a first-principles-based effective Hamiltonian method leads to the discovery of many fingerprints of a BKT phase existing in-between the ferroelectric and paraelectric states of 1UC tin tellurium being fully relaxed. Moreover, epitaxial strain is found to have dramatic consequences on the temperature range of such BKT phase for the 1UC SnTe. Consequently, our predictions extend the playground of BKT theory to a novel class of functional materials, and demonstrate that strain is an effective tool to alter BKT characteristics there.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linear versus nonlinear electro-optic effects in materials

Two schemes are proposed to compute the nonlinear electro-optic (EO) tensor for the first time. In the first scheme, we compute the linear EO tensor of the structure under a finite electric field, while we compute the refractive index of the structure under a finite electric field in the second scheme. Such schemes are applied to Pb(Zr,Ti)O$_{3}$ and BaTiO$_{3}$ ferroelectric oxides. It is found to reproduce a recently observed feature, namely why Pb(Zr$_{0.52}$Ti$_{0.48}$)O$_{3}$ adopts a mostly linear EO response while BaTiO$_{3}$ exhibits a strongly nonlinear conversion between electric and optical properties. Furthermore, the atomistic insight provided by the proposed ab-initio scheme reveals the origin of such qualitatively different responses, in terms of the field-induced behavior of the frequencies of some phonon modes and of some force constants.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Possible Kitaev Quantum Spin Liquid State in 2D Materials with S=3/2

Quantum spin liquids (QSLs) form an extremely unusual magnetic state in which the spins are highly correlated and fluctuate coherently down to the lowest temperatures, but without symmetry breaking and without the formation of any static long-range-ordered magnetism. Such intriguing phenomena are not only of great fundamental relevance in themselves, but also hold the promise for quantum computing and quantum information. Among different types of QSLs, the exactly solvable Kitaev model is attracting much attention, with most proposed candidate materials, e.g., RuCl$_3$ and Na$_2$IrO$_3$, having an effective $S$=1/2 spin value. Here, via extensive first-principle-based simulations, we report the investigation of the Kitaev physics and possible Kitaev QSL state in epitaxially strained Cr-based monolayers, such as CrSiTe$_3$, that rather possess a $S$=3/2 spin value. Our study thus extends the playground of Kitaev physics and QSLs to 3$d$ transition metal compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tunable spin textures in polar antiferromagnetic hybrid organic inorganic perovskites by electric and magnetic fields

The hybrid organic inorganic perovskites (HOIPs) have attracted much attention for their potential applications as novel optoelectronic devices. Remarkably, the Rashba band splitting, together with specific spin orientations in k space (i.e., spin texture), has been found to be relevant for the optoelectronic performances. In this work, by using first principles calculations and symmetry analyses, we study the electric polarization, magnetism, and spin texture properties of the antiferromagnetic (AFM) HOIP ferroelectric TMCM_MnCl3 (TMCM = (CH3)3NCH2Cl, trimethylchloromethyl ammonium). This recently synthesized compound is a prototype of order disorder and displacement-type ferroelectric with a large piezoelectric response, high ferroelectric transition temperature, and excellent photoluminescence properties [You et al., Science 357, 306 (2017)]. The most interesting result is that the inversion symmetry breaking coupled to the spin orbit coupling gives rise to a Rashba-like band splitting and a related robust persistent spin texture (PST) and/or typical spiral spin texture, which can be manipulated by tuning the ferroelectric or, surprisingly, also by the AFM magnetic order parameter. The tunability of spin texture upon switching of AFM order parameter is largely unexplored and our findings not only provide a platform to understand the physics of AFM spin texture but also support the AFM HOIP ferroelectrics as a promising class of optoelectronic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-dimensional Topological Semimetals Protected by Symmorphic Symmetries

Two-dimensional (2D) band crossing semimetals (BCSMs) could be used to build a range of novel nanoscale devices such as superlenses and transistors. We find that symmorphic symmetry can protect a new type of robust 2D BCSMs, unlike the previously proposed 2D essential BCSMs protected by non-symmorphic symmetry [Young et al., Phys. Rev. Lett. 115, 126803 (2015)]. This type of symmorphic symmetry protected (SSP) 2D essential BCSMs cannot be pair annihilated without destroying the crystalline symmetries, as opposed to the 2D BCSMs caused by the accidental band crossing. Through group theory analysis, we find that 2D SSP BCSMs can only exist at the K (K') point of Brillouin zone (BZ) of four layer groups and identify nonmagnetic 2D FeB2 as a candidate. Interestingly, nonmagnetic 2D SSP BCSMs can host a single pair of band crossing points (BCPs), whereas nonmagnetic three-dimensional (3D) Weyl semimetals (WSMs) have at least two pairs of band crossing Weyl points. It is found that the single pair of BCPs are robust against any kinds of strain. Furthermore, our calculation suggests that essential 2D SSP BCSMs can be used to realize electric field control of spin-texture, thus are promising candidates for spintronic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Band splitting with vanishing spin polarizations in noncentrosymmetric crystals

The Dresselhaus and Rashba effects are well-known phenomena in solid-state physics, in which spin-orbit coupling (SOC) splits spin-up and spin-down energy bands of nonmagnetic non-centrosymmetric crystals. Here, we discover a new phenomenon, dubbed as band splitting with vanishing spin polarizations (BSVSP), in which, as usual, SOC splits the energy bands in nonmagnetic non-centrosymmetric systems; surprisingly, however, both split bands show no net spin polarization along certain high-symmetry lines in the Brillouin zone. In order to rationalize this phenomenon, we propose a new classification of point groups into pseudo-polar and non-pseudo-polar groups. By means of first-principles simulations, we demonstrate that BSVSP can take place in both symmorphic (e.g., bulk GaAs) and non-symmorphic systems (e.g., two dimensional ferroelectric SnTe). Furthermore, we propose a novel linear magnetoelectric coupling in reciprocal space, which could be employed to tune the spin polarization with an external electric field. The BSVSP effect and its manipulation could therefore pave a new way to novel spintronic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magnetic Skyrmion State in Janus Monolayers of Chromium Trihalides Cr(I,X)3

Magnetic skyrmions are nano-scale spin structures that are promising for ultra-dense memory and logic devices. Recent progresses in two-dimensional magnets encourage the idea to realize skyrmionic states in freestanding monolayers. However, monolayers such as CrI3 lack Dzyaloshinskii-Moriya interactions (DMI) and thus do not naturally exhibit skyrmions but rather a ferromagnetic state. Here we propose the fabrication of Cr(I,X)3 Janus monolayers, in which the Cr atoms are covalently bonded to the underlying I ions and top-layer Br or Cl atoms. By performing first-principles calculations and Monte-Carlo simulations, we identify strong enough DMI, which leads to not only helical cycloid phases, but also to intrinsic skyrmionic states in Cr(I,Br)3 and magnetic-field-induced skyrmions in Cr(I,Cl)3.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Magnetic Properties from the Viewpoints of Electronic Hamiltonian: Spin Exchange Parameters, Spin Orientation and Spin-Half Misconception

In this chapter we review the quantitative and qualitative aspects of describing the properties of magnetic solids on the basis of electronic Hamiltonian. We show that a spin Hamiltonian approach becomes consistent with an electronic Hamiltonian approach if the spin lattice and its associated spin exchange parameters, to be used for the spin Hamiltonian, are determined by the energy-mapping analysis based on DFT calculations. The preferred spin orientation (i.e., the magnetic anisotropy) of a magnetic ion is not predicted by a spin Hamiltonian because it does not include the orbital degree of freedom explicitly. In contrast, the magnetic anisotropy is readily predicted by electronic structure theories employing both orbital and spin degrees of freedom, if one takes into consideration the spin-orbit coupling (SOC). It was shown that the preferred spin orientation of a magnetic ion can be predicted and understood in terms of the HOMO-LUMO interactions of the magnetic ion by taking SOC as perturbation. A spin Hamiltonian gives rise to the spin-half misconception, namely, the blind belief that spin-half magnetic ions do not possess magnetic anisotropy that arise from SOC. This misconception is a direct consequence from the limitedness of a spin Hamiltonian that it lacks the orbital degree of freedom. We show that the magnetic properties of 5d ion oxides are better explained by the LS-coupling than by the jj-coupling scheme of SOC, that the spin-orbital entanglement of 5d ions is not as strong as has been assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin orientations of the spin-half Ir4+ ions in Sr3NiIrO6, Sr2IrO4 and Na2IrO3: Density functional, perturbation theory and Madelung potential analyses

The spins of the low-spin Ir4+ (S = 1/2, d5) ions at the octahedral sites of the oxides Sr3NiIrO6, Sr2IrO4 and Na2IrO3 exhibit preferred orientations with respect to their IrO6 octahedra. We evaluated the magnetic anisotropies of these S = 1/2 ions on the basis of DFT calculations including spin-orbit coupling (SOC), and probed their origin by performing perturbation theory analyses with SOC as perturbation within the LS coupling scheme. The observed spin orientations of Sr3NiIrO6 and Sr2IrO4 are correctly predicted by DFT calculations, and are accounted for by the perturbation theory analysis. As for the spin orientation of Na2IrO3, both experimental studies and DFT calculations have not been unequivocal. Our analysis reveals that the Ir4+ spin orientation of Na2IrO3 should have nonzero components along the c- and a-axes directions. The spin orientations determined by DFT calculations are sensitive to the accuracy of the crystal structures employed, which is explained by perturbation theory analyses when interactions between adjacent Ir4+ ions are taken into consideration. There are indications implying that the 5d electrons of Na2IrO3 are less strongly localized compared with those of Sr3NiIrO6 and Sr2IrO4. This implication was confirmed by showing that the Madelung potentials of the Ir4+ ions are less negative in Na2IrO3 than in Sr3NiIrO6, Sr2IrO4. Most transition-metal S = 1/2 ions do have magnetic anisotropies because the SOC induces interactions among their crystal-field split d-states, and the associated mixing of the states modifies only the orbital parts of the states. This finding cannot be mimicked by a spin Hamiltonian because this model Hamiltonian lacks the orbital degree of freedom, thereby leading to the spin-half syndrome. The spin-orbital entanglement for the 5d spin-half ions Ir4+ is not as strong as has been assumed lately.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structure and Composition of the 200 K-Superconducting Phase of H2S under Ultrahigh Pressure: The Perovskite (SH-)(H3S+)

H2S is converted under ultrahigh pressure (> 110 GPa) to a metallic phase that becomes superconducting with a record Tc of 200 K. It has been proposed that the superconducting phase is body-centered cubic H3S ( Im3m , a = 3.089 Å) resulting from a decomposition reaction 3H2S --> 2H3S + S. The analogy of H2S and H2O leads us to a very different conclusion. The well-known dissociation of water into H3O+ and OH- increases by orders of magnitude under pressure. An equivalent behavior of H2S is anticipated under pressure with the dissociation, 2H2S --> H3S+ + SH- forming a perovskite structure (SH-)(H3S+), which consists of corner-sharing SH6 octahedra with SH- at each A-site (i.e., the center of each S8 cube). Our DFT calculations show that the perovskite (SH-)(H3S+) is thermodynamically more stable than the Im3m structure of H3S, and suggest that the A-site H atoms are most likely fluxional even at Tc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two Dimensional Phosphorus Oxides as Energy and Information Materials

Two-dimensional (2D) black phosphorus (i.e., phosphorene) has become a rising star in electronics. Recently, 2D phosphorus oxides with higher stability have been synthesized. In this work, we systematically explore the structures and properties of 2D phosphorus oxides on the basis of global optimization approach and first-principles calculations. We find that the structural features of 2D phosphorus oxides PxOy vary with the oxygen concentration. When the oxygen content is poor, the most stable 2D PxOy can be obtained by adsorbing O atoms on phosphorene. However, when the oxygen concentration becomes rich, stable structures are no longer based on phosphorene and will contain P-O-P motifs. For the 2D P4O4, we find that it has a direct band gap (about 2.24 eV), good optical absorption, and high stability in water, suggesting that it may be good candidate for photochemical water splitting application. Interestingly, 2D P2O3 adopt two possible stable ferroelectric structures (P2O3-I and P2O3-II) as the lowest energy configurations within a given layer thickness. The electric polarizations of P2O3-I and P2O3-II are perpendicular and parallel to the lateral plane, respectively. We propose that 2D P2O3 could be used in a novel nanoscale multiple-state memory device. Our work suggests that 2D phosphorus oxides may be excellent functional materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Giant biquadratic interaction induced magnetic anisotropy in the iron-based superconductor AxFe2-ySe2

The emergence of the electron-pocket only iron-based superconductor AxFe2-ySe2 (A = alkali metal) challenges the Fermi-surface nesting picture established in iron-pnictides. It was widely believed that magnetism is correlated with the superconductivity in AxFe2-ySe2. Unfortunately, the highly anisotropic exchange parameters and the disagreement between theoretical calculations and experimental results triggered a fierce debate on the nature of magnetism in AxFe2-ySe2. Here we find that the strong magnetic anisotropy is from the anisotropic biquadratic interaction. In order to accurately obtain the magnetic interaction parameters, we propose a universal method, which does not need including other high energy configurations as did in conventional energy mapping method. We show that our model successfully captures the magnetic interactions in AxFe2-ySe2 and correctly predicts the spin wave spectrum, in quantitative agreement with the experimental observation. These results suggest that the local moment picture, including the biquadratic term, can describe accurately the magnetic properties and spin excitations in AxFe2-ySe2, which sheds new light on the future study of the high-Tc iron-based superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Oxygen Vacancy Induced Flat Phonon Mode at FeSe /SrTiO3 interface

A high-frequency optical phonon mode of SrTiO3 (STO) was found to assist the high-temperature superconductivity observed recently at the interface between monolayer FeSe and STO substrate. However, the origin of this mode is not clear. Through first-principles calculations, we find that there is a novel polar phonon mode on the surface layers of the STO substrate, which does not exist in the STO crystals. The oxygen vacancies near the FeSe/STO interface drives the dispersion of this phonon mode to be flat and lowers its energy, whereas the charge transfer between STO substrate and FeSe monolayer further reduces its energy to 81 meV. This energy is in good agreement with the experimental value fitted by Lee et al. for the phonon mode responsible for the observed replica band separations and the increased superconducting gap. The oxygen-vacancy-induced flat and polar phonon mode provides clues for understanding the origin of high Tc superconductivity at the FeSe/STO interface.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of Silicon-Based Layered Structures for Optoelectronic Applications

A method based on the particle swarm optimization (PSO) algorithm is presented to design quasi-two-dimensional (Q2D) materials. With this development, various single-layer and bi-layer materials in C, Si, Ge, Sn, and Pb were predicted. A new Si bi-layer structure is found to have a much-favored energy than the previously widely accepted configuration. Both single-layer and bi-layer Si materials have small band gaps, limiting their usages in optoelectronic applications. Hydrogenation has therefore been used to tune the electronic and optical properties of Si layers. We discover two hydrogenated materials of layered Si8H2 and Si6H2 possessing quasi-direct band gaps of 0.75 eV and 1.59 eV, respectively. Their potential applications for light emitting diode and photovoltaics are proposed and discussed. Our study opened up the possibility of hydrogenated Si layered materials as next-generation optoelectronic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of spin orientations in terms of HOMO-LUMO interactions using spin-orbit coupling as perturbation

The preferred spin orientation of a magnetic ion can be predicted on the basis of density functional theory (DFT) calculations including electron correlation and spin-orbit coupling (SOC). However, most chemists and physicists are unaware of how the observed and/or calculated spin orientations are related to the local electronic structures of the magnetic ions. The objective of this article is to provide a conceptual framework of thinking about and predicting the preferred spin orientation of a magnetic ion by examining the relationship between the spin orientation and the local electronic structure of the ion. In general, a magnetic ion (i.e., an ion possessing unpaired spins) in a solid or a molecule is surrounded with main-group ligand atoms to form a polyhedron, and the d-states of the polyhedron are split because the antibonding interactions of the metal d-orbitals with the p orbitals of the surrounding ligands depend on the symmetries of the orbitals involved. The magnetic ion of the polyhedron has a certain preferred spin direction because its split d-states interact among themselves under SOC and because the energy lowering associated with the SOC-induced interactions depends on spin orientation. The preferred spin direction can be readily predicted on the basis of perturbation theory, in which the SOC is taken as perturbation and the split d-states as unperturbed states, by inspecting the magnetic quantum numbers of its d-orbitals present in the HOMO and LUMO of the polyhedron. Experimentally, the determination of the preferred spin orientations of magnetic ions requires a sophisticated level of experiments. Theoretically, it requires an elaborate level of DFT electronic structure calculations. We show that the outcomes of such experimental measurements and theoretical calculations can be predicted by a simple perturbation theory analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structural Evolution and Optoelectronic Applications of Multilayer Silicene

Despite the recent progress on two-dimensional multilayer materials (2DMM) with weak interlayer interactions, the investigation on 2DMM with strong interlayer interactions is far from its sufficiency. Here we report on first-principles calculations that clarify the structural evolution and optoelectronic properties of such a 2DMM, multilayer silicene. With our newly developed global optimization algorithm, we discover the existence of rich dynamically stable multilayer silicene phases, the stability of which is closely related to the extent of sp3 hybridization that can be evaluated by the average bonds and effective bond angles. The stable Si(111) surface structures are obtained when the silicene thickness gets up to four, showing the critical thickness for the structural evolution. We also find that the multilayer silicene with pi-bonded surfaces present outstanding optoelectronic properties for the solar cells and optical fiber communications due to the incorporation of sp2-type bonds in the sp3-type bonds dominated system. This study is helpful to complete the picture of structure and related property evolution of 2DMM with strong interlayer interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Antiferromagnetic ground state with pair-checkboard order in FeSe

Monolayer FeSe thin film grown on SrTiO_{3}(001) (STO) shows the sign of T_{c}> 77 K, which is higher than the T_{c}-record of 56 K for the bulk FeAs-based superconductors. However, little is known about the magnetic ground state of FeSe, which should be closely related to its unusual superconductivity. Previous studies presume the collinear stripe antiferromagnetic (AFM) state as the ground state of FeSe, same to that in FeAs superconductors. Here we find a novel magnetic order named "pair-checkboard AFM" as the magnetic ground state of tetragonal FeSe. The novel pair-checkboard order results from the interplay between the nearest, the next-nearest and the unnegligible next-next-nearest neighbor magnetic exchange couplings of Fe atoms. The monolayer FeSe in pair-checkbord order shows an unexpected insulating behavior with a Dirac-cone-like band structure related to the specific orbital order of d_{xz} and d_{yz} characters of Fe atoms, which could explain recently observed insulator-superconductor transition. The present results cast new insights on the magnetic ordering in FeSe monolayer and its derived superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Most spin-1/2 transition-metal ions do have single ion anisotropy

The cause for the preferred spin orientation in magnetic systems containing spin-1/2 transition-metal ions was explored by studying the origin of the easy-plane anisotropy of the spin-1/2 Cu2+ ions in CuCl2.2H2O, LiCuVO4, CuCl2 and CuBr2 on the basis of density functional theory and magnetic dipole-dipole energy calculations as well as a perturbation theory treatment of the spin-orbit coupling. We find that the spin orientation observed for these spin-1/2 ions is not caused by their anisotropic spin exchange interactions, nor by their magnetic dipole-dipole interactions, but by the spin-orbit coupling associated with their crystal-field split d-states. Our study also predicts in-plane anisotropy for the Cu2+ ions of Bi2CuO4 and Li2CuO2. The results of our investigations dispel the mistaken belief that magnetic systems with spin-1/2 ions have no magnetic anisotropy induced by spin-orbit coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Origin of Polar Distortion in LiNbO3-type "Ferroelectric" Metals: Role of A-site Instability and Short-Range Interactions

Since conduction electrons of a metal screen effectively the local electric dipole moments, it was widely believed that the ferroelectric-like distortion cannot occur in metals. Recently, metallic LiOsO3, was discovered to be the first clear-cut example of an Anderson-Blount "ferroelectric" metal, which at 140 K undergoes a ferroelectric-like structural transition similar to insulating LiNbO3. This is very surprising because the mechanisms for structural phase transitions are usually quite distinct in metals and insulators. Through performing first principles calculations, here we reveal that the local polar distortion in LiOsO3 is solely due to the instability of the A-site Li atom, in contrast to the LiNbO3 case where the second order Jahn-Teller effect of the B-site Nb ion also plays an additional role. More importantly, the "ferroelectric"-like long range order of the local polar distortion is found to be due to the predominantly ferroelectric short-range pair interactions between the local polar modes which are not screened by conduction electrons. Furthermore, we predict that LiNbO3-type MgReO3 is also a "ferroelectric" metal, but with a much higher structural transition temperature by 391 K than LiOsO3. Our work not only unravels the origin of FE-like distortion in LiNbO3-type "ferroelectric" metals, but also provides clue for designing other multi-functional "ferroelectric" metals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Understanding the Magnetic Puzzles of Double Perovskites A2FeOsO6 (A=Ca, Sr)

Double perovskites Sr2FeOsO6 and Ca2FeOsO6 show puzzling magnetic properties, the former a low-temperature antiferromagnet while the later a high-temperature insulating ferrimagnet. Here, in order to understand the underlying mechanism, we have investigated the frustrated magnetism of A2FeOsO6 by employing density functional theory and maximally-localized Wannier functions. We find that lattice distortion enhances the antiferromagnetic nearest-neighboring Fe-O-Os interaction but weakens the antiferromagnetic interactions through the Os-O-O-Os and Fe-O-Os-O-Fe paths, which is responsible for the magnetic transition from the low-temperature antiferromagnetism to the high-temperature ferrimagnetism with the decrease of the radius of the A2+ ions. We also discuss the 5d3-3d5 superexchange and propose such superexchange is intrinsically antiferromagnetic instead of the expected ferromagnetic. Our work illustrate that the magnetic frustration can be effectively relieved by lattice distortion, which provides another dimension to tune the complex magnetism in other 3d-5d (4d) double perovskites.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Prediction of (TiO2)x(Cu2O)y Alloys for Photoelectrochemical Water Splitting

The formation of (TiO2)x(Cu2O)y solid-solutions are investigated using a global optimization evolutionary algorithm. First-principles calculations based on density functional theory are then used to gain insight into the electronic properties of these alloys. We find that: (i) Ti and Cu in (TiO2)x(Cu2O)y alloys have similar local environments as in bulk TiO2 and Cu2O except for (TiO2)(Cu2O) which has some trigonal-planar Cu ions. (ii) The predicted optical band gaps are around 2.1 eV (590 nm), thus having much better performance for the absorption of visible light compared with both binary oxides. (iii) (TiO2)2(Cu2O) has the lowest formation energy amongst all studied alloys and the positions of its band edges are found to be suitable for solar-driven water splitting applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The interfacial effects on the spin density wave in FeSe/SrTiO3 thin film

Recently, the signs of both superconducting transition temperature (Tc) beyond 60 K and spin density wave (SDW) have been observed in FeSe thin film on SrTiO3 (STO) substrate, which suggests a strong interplay between superconductivity and magnetism. With the first-principles calculations, we find that the substrate-induced tensile strain tends to stabilize the SDW state in FeSe thin film by enhancing of the next-nearest-neighbor superexchange antiferromagnetic interaction bridged through Se atoms. On the other hand, we find that when there are oxygen vacancies in the substrate, the significant charge transfer from the substrate to the first FeSe layer would suppress the magnetic order there, and thus the high-temperature superconductivity could occur. In addition, the stability of the SDW is lowered when FeSe is on a defect-free STO substrate due to the redistribution of charges among the Fe 3d-orbitals. Our results provide a comprehensive microscopic explanation for the recent experimental findings, and build a foundation for the further exploration of the superconductivity and magnetism in this novel superconducting interface.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unveiling the Origin of the Insulating Ferromagnetism in LaMnO3 Thin Film

By combining genetic algorithm optimizations, first-principles calculations and the double-exchange model studies, we have unveiled that the exotic insulating ferromagnetism in LaMnO3 thin film originates from the previously unreported G-type d_{3z^2-r^2}/d_{x^2-y^2} orbital ordering. An insulating gap opens as a result of both the orbital ordering and the strong electron-phonon coupling. Therefore, there exist two strain induced phase transitions in the LaMnO3 thin film, from the insulating A-type antiferromagnetic phase to the insulating ferromagnetic phase and then to the metallic ferromagnetic phase. These phase transitions may be exploited in tunneling magnetoresistance and tunneling electroresistance related devices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

What Is Armchair-Zigzag Grain Boundary Structure in Graphene?

We have developed a new global optimization method for the determination of interface structure based on the differential evolution algorithm. Here, we applied this method to search for the ground state atomic structures of the grain boundary between the armchair and zigzag oriented graphene. We find two new grain boundary structures with considerably lower formation energy of about 1 eV/nm than that of the previously widely used structural models. These newly predicted structures show better mechanical property under external uniaxial strain, and distinguishable scanning tunneling microscope features, compared with the previous structural models. Our results provide important new information for the determination of grain boundary structures and henceforth the electronic properties of defected graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ordered Semiconducting Nitrogen-Graphene Alloys

The interaction between substitutional nitrogen atoms in graphene is studied by performing first principles calculations. The nearest neighbor interaction between nitrogen dopants is highly repulsive because of the strong electrostatic repulsion between nitrogen atoms, which prevents the full phase separation in nitrogen doped graphene. Interestingly, there are two relatively stable nitrogen-nitrogen pairs due to the anisotropy charge redistribution induced by nitrogen doping. We reveal two stable semiconducting ordered N doped graphene structures C3N and C12N through the cluster expansion technique and particle swarm optimization method. In particular, C12N has a direct band gap of 0.98 eV. The heterojunctions between C12N and graphene nanoribbons might be promising organic solar cells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Si3AlP: A new promising material for solar cell absorber

First-principles calculations are performed to study the structural and optoelectronic properties of the newly synthesized nonisovalent and lattice-matched (Si2)0.6(AlP)0.4 alloy [T. Watkins et al., J. Am. Chem. Soc. 2011, 133, 16212.] We find that the ordered CC-Si3AlP with a basic unit of one P atom surrounded by three Si atoms and one Al atom is the most stable one within the experimentally observed unit cell.1 Si3AlP has a larger fundamental band gap and a smaller direct band gap than Si, thus it has much higher absorption in the visible light region. The calculated properties of Si3AlP suggest that it is a promising candidate for improving the performance of the existing Si-based solar cells. The understanding on the stability and band structure engineering obtained in this study is general and can be applied for future study of other nonisovalent and lattice-matched semiconductor alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin reorientation in the square-lattice antiferromagnets RMnAsO (R = Ce, Nd): Density functional analysis of the spin exchange interactions between the rare-earth and transition-metal ions

The spin reorientation (SR) phenomenon of the square-lattice antiferromagnets RMnAsO (R = Ce, Nd) was investigated by analyzing the spin exchange interactions between the rare-earth and transition-metal ions (R3+ and Mn2+, respectively) on the basis of density functional calculations. It is found that the symmetry and strength of the Dzyaloshinskii-Moriya (DM) interaction are determined primarily by the partially filled 4f states of the R3+ ions, and that the DM and biquadratic (BQ) exchanges between the R3+ and Mn2+ ions are unusually strong and control the observed spin reorientation phenomenon. Below their SR temperature, the Mn2+ and Ce3+ moments are orthogonal in CeMnAsO but are collinear in NdMnAsO, because the DM interaction dominates over the BQ interaction for CeMnAsO while the opposite is the case for NdMnAsO. Experiments designed to test the implications of our findings are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unexpected large thermal rectification in asymmetric grain boundary of graphene

We have investigated the lattice thermal transport across the asymmetric tilt grain boundary between armchair and zigzag graphene by nonequilibrium molecular dynamics (NEMD). We have observed significant temperature drop and ultra-low temperature-dependent thermal boundary resistance. More importantly, we find an unexpected thermal rectification phenomenon. The thermal conductivity and Kapitza conductance is direction-dependent. The effect of thermal rectification could be amplified by increasing the difference of temperature imposed on two sides. Our results propose a promising kind of thermal rectifier and phonon diodes based on polycrystalline graphene without delicate manipulation of the atomic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Controlling doping in graphene through a SiC substrate: A first-principles study

Controlling the type and density of charge carriers by doping is the key step for developing graphene electronics. However, direct doping of graphene is rather a challenge. Based on first-principles calculations, a concept of overcoming doping difficulty in graphene via substrate is reported.We find that doping could be strongly enhanced in epitaxial graphene grown on silicon carbide substrate. Compared to free-standing graphene, the formation energies of the dopants can decrease by as much as 8 eV. The type and density of the charge carriers of epitaxial graphene layer can be effectively manipulated by suitable dopants and surface passivation. More importantly, contrasting to the direct doping of graphene, the charge carriers in epitaxial graphene layer are weakly scattered by dopants due to the spatial separation between dopants and the conducting channel. Finally, we show that a similar idea can also be used to control magnetic properties, for example, induce a half-metallic state in the epitaxial graphene without magnetic impurity doping.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Density functional theory analysis of the interplay between Jahn-Teller instability, uniaxial magnetism, spin arrangement, metal-metal interaction and spin-orbit coupling in Ca3CoMO6 (M = Co, Rh, Ir)

In the isostructural oxides Ca3CoMO6 (M = Co, Rh, Ir), the CoMO6 chains made up of face-sharing CoO6 trigonal prisms and MO6 octahedra are separated by Ca atoms. We analyzed the magnetic and electronic properties of these oxides on the basis of density functional theory calculations including on-site repulsion and spin-orbit coupling, and examined the essential one-electron pictures hidden behind results of these calculations. Our analysis reveals an intimate interplay between Jahn-Teller instability, uniaxial magnetism, spin arrangement, metal-metal interaction, and spin-orbit coupling in governing the magnetic and electronic properties of these oxides. These oxides undergo a Jahn-Teller distortion but their distortions are weak, so that their trigonal-prism Con+ (n = 2, 3) ions still give rise to strong easy-axis anisotropy along the chain direction. As for the d-state split pattern of these ions, the electronic and magnetic properties of Ca3CoMO6 (M = Co, Rh, Ir) are consistent with d0 < (d2, d-2) < (d1, d-1), but not with (d2, d-2) < d0 < (d1, d-1). The trigonal-prism Co3+ ion in Ca3Co2O6 has the L = 2 configuration (d0)1(d2, d-2)3(d1, d-1)2 because of the metal-metal interaction between adjacent Co3+ ions in each Co2O6 chain, which is mediated by their z2 orbitals, and the spin-orbit coupling of the trigonal-prism Co3+ ion. The spins in each CoMO6 chain of Ca3CoMO6 prefer the ferromagnetic arrangement for M = Co and Rh, but the antiferromagnetic arrangement for M = Ir. The octahedral M4+ ion of Ca3CoMO6 has the (1a)1(1e)4 configuration for M = Rh but the (1a)2(1e)3 configuration for M = Ir, which arises from the difference in the spin-orbit coupling of the M4+ ions and the Co...M metal-metal interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First Principles Study of Adsorption of $O_{2}$ on Al Surface with Hybrid Functionals

Adsorption of $O_{2}$ molecule on Al surface has been a long standing puzzle for the first principles calculation. We have studied the adsorption of $O_{2}$ molecule on the Al(111) surface using hybrid functionals. In contrast to the previous LDA/GGA, the present calculations with hybrid functionals successfully predict that $O_{2}$ molecule can be absorbed on the Al(111) surface with a barrier around 0.2$\thicksim$0.4 eV, which is in good agreement with experiments. Our calculations predict that the LUMO of $O_{2}$ molecule is higher than the Fermi level of the Al(111) surface, which is responsible for the barrier of the $O_{2}$ adsorption.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Layer and size dependence of thermal conductivity in multilayer graphene nanoribbons

Using non-equilibrium molecular dynamics method(NEMD), we have found that the thermal conductivity of multilayer graphene nanoribbons monotonously decreases with the increase of the number of layers, such behavior can be attributed to the phonon resonance effect of out-of-plane phonon modes. The reduction of thermal conductivity is found to be proportional to the layer size, which is caused by the increase of phonon resonance. Our results are in agreement with recent experiment on dimensional evolution of thermal conductivity in few layer graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Predicting Two-Dimensional Boron-Carbon Compounds by the Global Optimization Method

We adopt a global optimization method to predict two-dimensional (2D) nanostructures through the particle-swarm optimization (PSO) algorithm. By performing PSO simulations, we predict new stable structures of 2D boron-carbon (B-C) compounds for a wide range of boron concentrations. Our calculations show that: (1) All 2D B-C compounds are metallic except for BC3 which is a magic case where the isolation of carbon six-membered ring by boron atoms results in a semiconducting behavior. (2) For C-rich B-C compounds, the most stable 2D structures can be viewed as boron doped graphene structures, where boron atoms typically form 1D zigzag chains except for BC3 in which boron atoms are uniformly distributed. (3) The most stable 2D structure of BC has alternative carbon and boron ribbons with strong in-between B-C bonds, which possesses a high thermal stability above 2000K. (4) For B-rich 2D B-C compounds, there is a novel planar-tetracoordinate carbon motif with an approximate C2v symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetic structure and ferroelectric polarization of MnWO4 investigated by density functional calculations and classical spin analysis

The ordered magnetic states of MnWO4 at low temperatures were examined by evaluating the spin exchange interactions between the Mn2+ ions of MnWO4 on the basis of first principles density functional calculations and by performing classical spin analysis with the resulting spin exchange parameters. Our work shows that the spin exchange interactions are frustrated within each zigzag chain of Mn2+ ions along the c-direction and between such chains of Mn2+ ions along the a-direction. This explains the occurrence of a spiral-spin order along the c- and a-directions in the incommensurate magnetic state AF2, and that of a uudd spin order along the c- and a-directions in the commensurate magnetic state AF1. The ferroelectric polarization of MnWO4 in the spiral-spin state AF2 was examined by performing Berry phase calculations for a model superstructure to find that the ferroelectric polarization occurs along the b-direction, in agreement with experiment.

0 citations0 reviewsNo code linkedOpen access