Source author record

Jijun Zhao

Jijun Zhao 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

20works
7topics
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

20 published item(s)

preprint2026arXiv

Efficient Screening of Organic Singlet Fission Molecules Using Graph Neural Networks

Singlet fission (SF) provides a promising strategy for surpassing the Shockley-Queisser limit in photovoltaics. However, the identification of efficient SF materials is hindered by the limited availability of suitable molecular candidates and the high computational costs associated with conventional quantum-chemical methods for excited states. In this study, we introduce a high-throughput screening framework that integrates a graph neural network (GNN) with multi-level validation to accelerate the discovery of SF-active molecules. Trained on a previously reported FORMED database, the GNN achieves state-of-the-art accuracy in predicting SF-relevant excited-state properties, demonstrating a mean absolute error of about 0.1 eV for S1, T1, and T2 excitation energies. This capability facilitates the efficient screening of over 20 million molecular structures from both OE62 and QO2Mol databases. Our framework significantly reduces the computational demand associated with Time-Dependent Density Functional Theory validation by four orders of magnitude and identifies 180 potential SF molecules along with more than 1000 conformers. Subsequent assessments regarding synthetic accessibility, GW approximation and Bethe-Salpeter equation calculations further highlight a subset of experimentally feasible candidates among these SF candidates. The approach presented herein exemplifies an effective strategy for accelerating the discovery of functional molecules with optoelectronic applications.

preprint2026arXiv

Ir3Ge20: A 3n-Connected Cloverleaf-Shaped Supercluster

Group 14 Zintl clusters are promising molecular building blocks for nanoscale architecture. Endohedral variants, which encapsulate d/f-block metals within p-block semimetal cages, provide insights into intermetallic bonding and compound formation. In this study, experimental photoelectron spectroscopy and first-principles calculations were used to investigate the Ir-doped germanium cluster species. A new C2v symmetric building block, IrGe12, was identified, serving as the basis for designing the supercluster Ir3Ge20 with a cloverleaf-shaped, D3h symmetric architecture. This structure consists of three interconnected aromatic IrGe12 units linked by Ge-Ge sigma bonds, forming shielding cones. Ir3Ge20 follows the 5n rule with 100 valence electrons, featuring a core-Ir3 unit with a d10 closed-shell configuration sharing electrons with the Ge20 skeleton. The stability, chemical bonding, and aromaticity of Ir3Ge20 were confirmed, demonstrating a novel approach to precise atom manipulation in cluster-based materials and devices.

preprint2022arXiv

Eliminating edge electronic and phonon states of phosphorene nanoribbon by unique edge reconstruction

Edge termination plays a vital role in determining the properties of 2D materials. By performing compelling ab initio simulations, a lowest-energy U-edge [ZZ(U)] reconstruction is revealed in the bilayer phosphorene. Such reconstruction reduces 60% edge energy compared with the pristine one and occurs almost without energy barrier, implying it should be the dominating edge in reality. The electronic band structure of phosphorene nanoribbon with such reconstruction resembles that of intrinsic 2D layer, exhibiting nearly edgeless band characteristics. Although ZZ(U) changes the topology of phosphorene nanoribbon (PNR), simulated TEM, STEM and STM images indicates it is very hard to be identified. One possible identify method is IR/Raman analyses because ZZ(U) edge alters vibrational modes dramatically. Beyond, it also increases the thermal conductivity of PNR 1.4 and 2.3 times than the pristine and Klein edges.

preprint2022arXiv

Mechanical Anisotropy and Multiple Direction-Dependent Dirac States in the Synthesized Ag3C20 Monolayer

Recently, a 2D orthorhombic silver-organic framework, Ag3C20 monolayer, was synthesized by assembling organic molecules linked with multiple aryl-metal bonds. Herein, via first-principles study, we demonstrate that owing to the unique bonding feature, Ag3C20 monolayer not only exhibits strong mechanical anisotropy, but also possesses various tunable direction-dependent Dirac states. Around the Fermi level blow, the intrinsic Dirac points form two antiparallel quasi type-III nodal lines protected by mirror symmetry, which can further evolve into hybrid nodal loops under tiny strains. Intriguingly, near the Fermi level above, a special semi-Dirac state can emerge under a critical strain by merging two type-I Dirac cone, which harbors direction-dependent strongly localized fermions, normal massive carries, and ultrafast Dirac fermions at the same time. These findings suggest that the mechanically sensitive Ag3C20 monolayer is a promising 2D material to realize the interesting Dirac physics and highly anisotropic multiple carries transport.

preprint2020arXiv

Enhanced Ferromagnetism of CrI3 Bilayer by Self-Intercalation

Two-dimensional (2D) ferromagnets with high Curie temperature have long been the pursuit for electronic and spintronic applications. CrI3 is a rising star of intrinsic 2D ferromagnets, however, it suffers from weak exchange coupling. Here we propose a general strategy of self-intercalation to achieve enhanced ferromagnetism in bilayer CrI3. We showed that filling either Cr or I atoms into the van der Waals gap of stacked and twisted CrI3 bilayers can induce the double exchange effect and significantly strengthen the interlayer ferromagnetic coupling. According to our first-principles calculations, the intercalated native atoms act as covalent bridge between two CrI3 layers and lead to discrepant oxidation states for the Cr atoms. These theoretical results offer a facile route to achieve high-Curie-temperature 2D magnets for device implementation.

preprint2020arXiv

Experimental Realization of Two-Dimensional Buckled Lieb lattice

Two-dimensional (2D) materials with a Lieb lattice can host exotic electronic band structures. Such a system does not exist in nature, and it is also difficult to obtain in the laboratory due to its structural instability. Here, we experimentally realized a 2D system composed of a tin overlayer on an aluminum substrate by molecular beam epitaxy. The specific arrangement of Sn atoms on the Al(100) surface, which benefits from favorable interface interactions, forms a stabilized buckled Lieb lattice. Our theoretical calculations indicate a partially broken nodal line loop protected by its mirror reflection symmetry and a topologically nontrivial insulating state with a spin-orbital coupling (SOC) effect in the band structure of this Lieb lattice. The electronic structure of this system has also been experimentally characterized by scanning tunnelling spectroscopy and angle-resolved photoemmision spectroscopy. Our work provides an appealing method for constructing 2D quantum materials based on the Lieb lattice.

preprint2020arXiv

High-Curie-temperature ferromagnetism in bilayer CrI3 on bulk semiconducting substrates

Two-dimensional (2D) ferromagnetic (FM) semiconductors with high Curie temperature have long been pursued for electronic and spintronic applications. Here we provide a general strategy to achieve robust FM state in bilayer CrI3 of the monoclinic stacking, which intrinsically has interlayer antiferromagnetic (AFM) order and weak in-plane FM coupling. We showed that the proximity effect from bulk semiconducting substrates induces electronic doping and significantly increases the FM nearest-neighbor exchange for bilayer CrI3, leading to the AFM-to-FM transition for the interlayer spin configuration as well as enhanced intralayer FM coupling. By first-principles calculations and Monte Carlo simulations, bulk and 2D semiconductors providing different interaction strengths from strong covalent bonding to weak van der Waals (vdW) interaction with CrI3 are compared to thoroughly address the substrate effect on magnetic behavior and Curie temperature of bilayer CrI3. These theoretical results offer a facile route for direct synthesis of 2D ferromagnets on proper semiconducting substrates to achieve high Curie temperature for device implementation.

preprint2020arXiv

Selective C-C Coupling by Spatially Confined Dimeric Metal Centers

Direct conversion of carbon dioxide (CO2) to high-energy fuels and high-value chemicals is a fascinating sustainable strategy. For most of the current electrocatalysts for CO2 reduction, however, multi-carbon products are inhibited by large overpotentials and low selectivity. For practical applications, there remains a big gap of knowledge in proper manipulation of the C-C coupling process. Herein, we exploit dispersed 3d transition metal dimers as spatially confined dual reaction centers for selective reduction of CO2 to liquid fuels. Various nitrogenated holey carbon monolayers are shown to be promising templates to stabilize these metal dimers and dictate their electronic structures, allowing precise control of the catalytic activity and product selectivity. By comprehensive first-principles calculations, we screen the suitable transition metal dimers that universally have high activity for ethanol (C2H5OH). Furthermore, remarkable selectivity for C2H5OH against other C1 and C2 products is found for Fe2 dimer anchored on C2N monolayer. The correlation between the activity and d band center of the supported metal dimer as well as the role of electronic coupling between the metal dimer and the carbon substrates are thoroughly elucidated.

preprint2020arXiv

Ultra-high Hydrogen Storage Capacity of Holey Graphyne

Holey graphyne (HGY), a novel 2D single-crystalline carbon allotrope, was synthesized most recently by Castro-Stephens coupling reaction. The natural existing uniform periodic holes in the 2D carbon-carbon network demonstrate its tremendous potential application in the area of energy storage. Herein, we conducted density functional theory calculation to predict the hydrogen storage capacity of HGY sheet. It's found the Li-decorated single-layer HGY can serve as a promising candidate for hydrogen storage. Our numerical calculations demonstrate that Li atoms can bind strongly to the HGY sheet without the formation of Li clusters, and each Li atom can anchor four H2 molecules with the average adsorption energy about -0.22 eV/H2. The largest hydrogen storage capacity of the doped HGY sheet can arrive as high as 12.8 wt%, this value largely surpasses the target of the U. S. Department of Energy (9 wt%), showing the Li/HGY complex is an ideal hydrogen storage material at ambient conditions. In addition, we investigate the polarization mechanism of the storage media and and find that the polarization stemed from both the electric field induced by the ionic Li decorated on the HGY and the weak polarized hydrogen molecules dominated the H2 adsorption process.

preprint2019arXiv

Physical properties and device applications of graphene oxide

Graphene oxide (GO), the functionalized graphene with oxygenated groups (mainly epoxy and hydroxyl), has attracted resurgent interests in the past decade owing to its large surface area, superior physical and chemical properties, and easy composition with other materials via surface functional groups. Usually, GO is used as an important raw material for mass production of graphene via reduction. However, under different conditions, the coverage, types, and arrangements of oxygen-containing groups in GO can be varied, which give rise to excellent and controllable physical properties, such as tunable electronic and mechanical properties depending closely on oxidation degree, suppressed thermal conductivity, optical transparency and fluorescence, and nonlinear optical properties. Based on these outstanding properties, many electronic, optical, optoelectronic, and thermoelectric devices with high performance can be achieved on the basis of GO. Here we present a comprehensive review on recent progress of GO, focusing on the atomic structures, fundamental physical properties, and related device applications, including transparent and flexible conductors, field-effect transistors, electrical and optical sensors, fluorescence quenchers, optical limiters and absorbers, surface enhanced Raman scattering detectors, solar cells, light-emitting diodes, and thermal rectifiers.

preprint2018arXiv

Understanding the mechanical properties of reduced activation steels

Reduced activation ferritic/martensitic (RAFM) steels are structural materials with potential application in Generation-IV fission and fusion reactors. We use density-functional theory to scrutinize the micro-mechanical properties of the main alloy phases of three RAFM steels based on the body-centered cubic FeCrWVMn solid solution. We assess the lattice parameters and elastic properties of ferromagnetic $α$-Fe and Fe$_{91}$Cr$_{9}$, which are the main building blocks of the RAFM steels, and present a detailed analysis of the calculated alloying effects of V, Cr, Mn, and W on the mechanical properties of Fe$_{91}$Cr$_{9}$. The composition dependence of the elastic parameters is decomposed into electronic and volumetric contributions and studied for alloying levels that cover the typical intervals in RAFM steels. A linear superposition of the individual solute effects on the properties of Fe$_{91}$Cr$_{9}$ is shown to provide an excellent approximation for the \emph{ab initio} values obtained for the RAFM steels. The intrinsic ductility is evaluated through Rice's phenomenological theory using the surface and unstable stacking fault energies, and the predictions are contrasted with those obtained by empirical criteria. Alloying with V or W is found to enhance the ductility, whereas additional Cr or Mn turns the RAFM base alloys more brittle.

preprint2016arXiv

Alloying effect on the ideal tensile strength of ferromagnetic and paramagnetic bcc iron

Using \emph{ab initio} alloy theory formulated within the exact muffin-tin orbitals theory in combination with the coherent potential approximation, we investigate the ideal tensile strength (ITS) in the $[001]$ direction of bcc ferro-/ferrimagnetic (FFM) and paramagnetic (PM) Fe$_{1-x}M_{x}$ ($M=$ Al, V, Cr, Mn, Co, or Ni) random alloys. The ITS of ferromagnetic (FM) Fe is calculated to be $12.6$\,GPa, in agreement with available data, while the PM phase turns out to posses a significantly lower value of $0.7\,$GPa. Alloyed to the FM matrix, we predict that V, Cr, and Co increase the ITS of Fe, while Al and Ni decrease it. Manganese yields a weak non-monotonic alloying behavior. In comparison to FM Fe, the alloying effect of Al and Co to PM Fe is reversed and the relative magnitude of the ITS can be altered more strongly for any of the solutes. All considered binaries are intrinsically brittle and fail by cleavage of the $(001)$ planes under uniaxial tensile loading in both magnetic phases. We show that the previously established ITS model based on structural energy differences proves successful in the PM Fe-alloys but is of limited use in the case of the FFM Fe-based alloys. The different performance is attributed to the specific interplay between magnetism and volume change in response to uniaxial tension. We establish a strong correlation between the compositional effect on the ITS and the one on the shear elastic constant $C'$ for the PM alloys and briefly discuss the relation between hardenability and the ITS.

preprint2016arXiv

Band alignment of two-dimensional lateral heterostructures

Recent experimental synthesis of two-dimensional (2D) heterostructures opens a door to new opportunities in tailoring the electronic properties for novel 2D devices. Here, we show that a wide range of lateral 2D heterostructures could have a prominent advantage over the traditional three-dimensional (3D) heterostructures, because their band alignments are insensitive to the interfacial conditions. They should be at the Schottky-Mott limits for semiconductor-metal junctions and at the Anderson limits for semiconductor junctions, respectively. This fundamental difference from the 3D heterostructures is rooted in the fact that, in the asymptotic limit of large distance, the effect of the interfacial dipole vanishes for 2D systems. Due to the slow decay of the dipole field and the dependence on the vacuum thickness, however, studies based on first-principles calculations often failed to reach such a conclusion. Taking graphene/hexagonal-BN and MoS2/WS2 lateral heterostructures as the respective prototypes, we show that the converged junction width can be order of magnitude longer than that for 3D junctions. The present results provide vital guidance to high-quality transport devices wherever a lateral 2D heterostructure is involved.

preprint2016arXiv

Tensile strain-induced softening of iron at high temperature

In weakly ferromagnetic materials, already small changes in the atomic configuration triggered by temperature or chemistry can alter the magnetic interactions responsible for the non-random atomic-spin orientation. Different magnetic states, in turn, can give rise to substantially different macroscopic properties. A classical example is iron, which exhibits a great variety of properties as one gradually removes the magnetic long-range order by raising the temperature towards and beyond its Curie point of $T_{\text{C}}^{0}=1043$\,K. Using first-principles theory, here we demonstrate that uniaxial tensile strain can also destabilize the magnetic order in iron and eventually lead to a ferromagnetic to paramagnetic transition at temperatures far below $T_{\text{C}}^{0}$. In consequence, the intrinsic strength of the ideal single-crystal body-centered cubic iron dramatically weakens above a critical temperature of $\sim 500$\,K. The discovered strain-induced magneto-mechanical softening provides a plausible atomic-level mechanism behind the observed drop of the measured strength of Fe whiskers around $300-500$\,K. Alloying additions which have the capability to partially restore the magnetic order in the strained Fe lattice, push the critical temperature for the strength-softening scenario towards the magnetic transition temperature of the undeformed lattice. This can result in a surprisingly large alloying-driven strengthening effect at high temperature as illustrated here in the case of Fe-Co alloy.

preprint2015arXiv

Point defects in epitaxial silicene on Ag(111) surface

Silicene, a counterpart of graphene, has achieved rapid development due to its exotic electronic properties and excellent compatibility with the mature silicon-based semiconductor technology. Its low room-temperature mobility of about 100 cm2V-1s-1, however, inhibits device applications such as in field-effect transistors. Generally, defects and grain boundaries would act as scattering centers and thus reduce the carrier mobility. In this paper, the morphologies of various point defects in epitaxial silicene on Ag(111) surfaces have been systematically investigated using first-principles calculations combined with experimental scanning tunneling microscope (STM) observations. The STM signatures for various defects in epitaxial silicene on Ag(111) surface are identified. In particular, the formation energies of point defects in Ag(111)-supported silicene sheets show an interesting dependence on the superstructures, which, in turn, may have implications for controlling the defect density during the synthesis of silicene. Through estimating the concentrations of various point defects in different silicene superstructures, the mystery of the defective appearance of v13*v13 silicene in experiments is revealed, and 4*4 silicene sheet is thought to be the most suitable structure for future device applications.

preprint2015arXiv

Stacking fault energy of face-centered cubic metals: thermodynamic and ab initio approaches

The formation energy of the interface between face-centered cubic (fcc) and hexagonal close packed (hcp) structures is a key parameter in determining the stacking fault energy (SFE) of fcc metals and alloys using thermodynamic calculations. Often the contribution of the planar fault energy to the SFE has the same order of magnitude as the bulk part, and thus the lack of a precise information about it can become the limiting factor in thermodynamic predictions. Here, we differentiate between the actual interfacial energy for the coherent fcc(111)/hcp(0001) interface and the "pseudo-interfacial energy" that enters the thermodynamic expression for the SFE. Using first-principles calculations, we determine the coherent and pseudo- interfacial energies for six elemental metals (Al, Ni, Cu, Ag, Pt, and Au) and for three paramagnetic Fe-Cr-Ni alloys. Our results show that the two interfacial energies significantly differ from each other. We observe a strong chemistry dependence of both interfacial energies. The calculated pseudo-interfacial energies for the Fe-Cr-Ni steels agree well with the available literature data.

preprint2014arXiv

Anomalous ideal tensile strength of ferromagnetic Fe and Fe-rich alloys

Within the same failure mode, iron has the lowest ideal tensile strength among the transition metals crystallizing in the body-centered cubic structure. Here, we demonstrate that this anomalously low strength of Fe originates partly from magnetism and is reflected in unexpected alloying effects in dilute Fe(\emph{M}) (\emph{M} = Al, V, Cr, Mn, Co, Ni) binaries. We employ the structural energy difference and the magnetic pressure to disentangle the magnetic effect on the ideal tensile strength from the chemical effect. We find that the investigated solutes strongly alter the magnetic response of the Fe host from the weak towards a stronger ferromagnetic behavior, which is explained based on single-particle band energies.

preprint2014arXiv

Tuning the Band Gap in Silicene by Oxidation

Silicene monolayers grown on Ag(111) surfaces demonstrate a band gap that is tunable by oxygen adatoms from semimetallic to semiconducting type. By using low-temperature scanning tunneling microscopy, it is found that the adsorption configurations and amounts of oxygen adatoms on the silicene surface are critical for band-gap engineering, which is dominated by different buckled structures in R13xR13, 4x4, and 2R3x2R3 silicene layers. The Si-O-Si bonds are the most energy-favored species formed on R13xR13, 4x4, and 2R3x2R3 structures under oxidation, which is verified by in-situ Raman spectroscopy as well as first-principles calculations. The silicene monolayers retain their structures when fully covered by oxygen adatoms. Our work demonstrates the feasibility of tuning the band gap of silicene with oxygen adatoms, which, in turn, expands the base of available two-dimensional electronic materials for devices with properties that is hardly achieved with graphene oxide.

preprint2013arXiv

Ideal strength of random alloys from first-principles theory

The all-electron exact muffin-tin orbitals method in combination with the coherent-potential appproximation has been employed to investigate the ideal tensile strengths of elemental V, Mo solids and V- and Mo-based random solid solutions. The present ideal tensile strengths, calculated assuming isotropic Poisson contraction, are 16.1, 26.7 and 37.6 GPa for bcc V in the [001], [111] and [110] directions, respectively, and 26.7 GPa for bcc Mo in the [001] direction, which are all in good agreement with the available theoretical data. When a few percent Tc is introduced in Mo, it is found that the ideal strength decreases in the [001] direction. For the V-based alloys, Cr increases and Ti decreases the ideal tensile strength in all principal directions. Adding the same concentration of Cr and Ti to V leads to ternary alloys with similar ideal strength values as that of pure V. The alloying effects on the ideal strength is explained using the electronic band structure.

preprint2011arXiv

Graphene Nucleation on Transition Metal Surface: Structure Transformation and Role of the Metal Step Edge

The nucleation of graphene on a transition metal (TM) surface, either on a terrace or near a step edge, is systematically explored using density functional theory (DFT) calculations and applying the two-dimensional (2D) crystal nucleation theory. Careful optimization of the supported carbon clusters, CN (with size N ranging from 1 to 24), on the Ni(111) surface indicates a ground state structure transformation from a one-dimensional (1D) C chain to a two-dimensional (2D) sp2 C network at N ~ 10-12. Furthermore, the crucial parameters controlling graphene growth on the metal surface, nucleation barrier, nucleus size, and the nucleation rate on a terrace or near a step edge, are calculated. In agreement with numerous experimental observations, our analysis shows that graphene nucleation near a metal step edge is superior to that on a terrace. Based on our analysis, we propose the use of seeded graphene to synthesize high-quality graphene in large area.