Source author record

Jinlong Yang

Jinlong Yang 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

67works
15topics
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

67 published item(s)

preprint2026arXiv

DCVD: Dual-Channel Cross-Modal Fusion for Joint Vulnerability Detection and Localization

Software vulnerability detection plays a critical role in ensuring system security, where real-world auditing requires not only determining whether a function is vulnerable but also pinpointing the specific lines responsible. However, existing approaches either rely on a single information source -- sequential, structural, or semantic -- failing to jointly exploit the complementary strengths across modalities, or treat statement-level localization merely as a byproduct of function-level detection without explicit line-level supervision. To address these limitations, we propose DCVD (Dual-Channel Cross-Modal Vulnerability Detection), a unified framework that performs joint function-level detection and statement-level localization. DCVD extracts control-dependency and semantic features through two parallel branches and integrates them via contrastive alignment coupled with bidirectional cross-attention, effectively bridging the cross-modal representation gap. It further introduces explicit supervision signals at both the function and statement levels, enabling collaborative optimization across the two granularities. Extensive experiments on a large-scale real-world vulnerability benchmark demonstrate that DCVD consistently outperforms state-of-the-art methods on both function-level detection and statement-level localization. Our code is available at https://github.com/vinsontang1/DCVD.

preprint2026arXiv

VulTriage: Triple-Path Context Augmentation for LLM-Based Vulnerability Detection

Automated vulnerability detection is a fundamental task in software security, yet existing learning-based methods still struggle to capture the structural dependencies, domain-specific vulnerability knowledge, and complex program semantics required for accurate detection. Recent Large Language Models (LLMs) have shown strong code understanding ability, but directly prompting them with raw source code often leads to missed vulnerabilities or false alarms, especially when vulnerable and benign functions differ only in subtle semantic details. To address this, we propose VulTriage, a triple-path context augmentation framework for LLM-based vulnerability detection. VulTriage enhances the LLM input through three complementary paths: a Control Path that extracts and verbalizes AST, CFG, and DFG information to expose control and data dependencies; a Knowledge Path that retrieves relevant CWE-derived vulnerability patterns and examples through hybrid dense--sparse retrieval; and a Semantic Path that summarizes the functional behavior of the code before the final judgment. These contexts are integrated into a unified instruction to guide the LLM toward more reliable vulnerability reasoning. Experiments on the PrimeVul pair test set show that VulTriage achieves state-of-the-art performance, outperforming existing deep learning and LLM-based baselines on key pair-wise and classification metrics. Further ablation studies verify the effectiveness of each path, and additional experiments on the Kotlin dataset demonstrate the generalization ability of VulTriage under low-resource and class-imbalanced settings. Our code is available at https://github.com/vinsontang1/VulTriage

preprint2024arXiv

Complex-valued K-means clustering of interpolative separable density fitting algorithm for large-scale hybrid functional enabled \textit{ab initio} molecular dynamics simulations within plane waves

K-means clustering, as a classic unsupervised machine learning algorithm, is the key step to select the interpolation sampling points in interpolative separable density fitting (ISDF) decomposition. Real-valued K-means clustering for accelerating the ISDF decomposition has been demonstrated for large-scale hybrid functional enabled \textit{ab initio} molecular dynamics (hybrid AIMD) simulations within plane-wave basis sets where the Kohn-Sham orbitals are real-valued. However, it is unclear whether such K-means clustering works for complex-valued Kohn-Sham orbitals. Here, we apply the K-means clustering into hybrid AIMD simulations for complex-valued Kohn-Sham orbitals and use an improved weight function defined as the sum of the square modulus of complex-valued Kohn-Sham orbitals in K-means clustering. Numerical results demonstrate that this improved weight function in K-means clustering algorithm yields smoother and more delocalized interpolation sampling points, resulting in smoother energy potential, smaller energy drift and longer time steps for hybrid AIMD simulations compared to the previous weight function used in the real-valued K-means algorithm. In particular, we find that this improved algorithm can obtain more accurate oxygen-oxygen radial distribution functions in liquid water molecules and more accurate power spectrum in crystal silicon dioxide compared to the previous K-means algorithm. Finally, we describe a massively parallel implementation of this ISDF decomposition to accelerate large-scale complex-valued hybrid AIMD simulations containing thousands of atoms (2,744 atoms), which can scale up to 5,504 CPU cores on modern supercomputers.

preprint2023arXiv

Quantum circuit matrix product state ansatz for large-scale simulations of molecules

As in the density matrix renormalization group (DMRG) method, approximating many-body wave function of electrons using a matrix product state (MPS) is a promising way to solve electronic structure problems. The expressibility of an MPS is determined by the size of the matrices or in other words the bond dimension, which unfortunately should be very large in many cases. In this study, we propose to calculate the ground state energies of molecular systems by variationally optimizing quantum circuit MPS (QCMPS) with a relatively small number of qubits. It is demonstrated that with carefully chosen circuit structure and orbital localization scheme, QCMPS can reach a similar accuracy as that achieved in DMRG with an exponentially large bond dimension. QCMPS simulation of a linear molecule with 50 orbitals can reach the chemical accuracy using only 6 qubits at a moderate circuit depth. These results suggest that QCMPS is a promising wave function ansatz in the variational quantum eigensolver algorithm for molecular systems.

preprint2023arXiv

Quantum Neural Network Inspired Hardware Adaptable Ansatz for Efficient Quantum Simulation of Chemical Systems

The variational quantum eigensolver is a promising way to solve the Schrödinger equation on a noisy intermediate-scale quantum (NISQ) computer, while its success relies on a well-designed wavefunction ansatz. Compared to physically motivated ansatzes, hardware heuristic ansatzes usually lead to a shallower circuit, but it may still be too deep for an NISQ device. Inspired by the quantum neural network, we propose a new hardware heuristic ansatz where the circuit depth can be significantly reduced by introducing ancilla qubits, which makes a practical simulation of a chemical reaction with more than 20 atoms feasible on a currently available quantum computer. More importantly, the expressibility of this new ansatz can be improved by increasing either the depth or the width of the circuit, which makes it adaptable to different hardware environments. These results open a new avenue to develop practical applications of quantum computation in the NISQ era.

preprint2022arXiv

Divide-and-conquer variational quantum algorithms for large-scale electronic structure simulations

Exploring the potential application of quantum computers in material design and drug discovery has attracted a lot of interest in the age of quantum computing. However, the quantum resource requirement for solving practical electronic structure problems are far beyond the capacity of near-term quantum devices. In this work, we integrate the divide-and-conquer (DC) approaches into the variational quantum eigensolver (VQE) for large-scale quantum computational chemistry simulations. Two popular divide-and-conquer schemes, including many-body expansion~(MBE) fragmentation theory and density matrix embedding theory~(DMET), are employed to divide complicated problems into many small parts that are easy to implement on near-term quantum computers. Pilot applications of these methods to systems consisting of tens of atoms are performed with adaptive VQE algorithms. This work should encourage further studies of using the philosophy of DC to solve electronic structure problems on quantum computers.

preprint2022arXiv

Exploring accurate potential energy surfaces via integrating variational quantum eigensovler with machine learning

The potential energy surface (PES) is crucial for interpreting a variety of chemical reaction processes. However, predicting accurate PESs with high-level electronic structure methods is a challenging task due to the high computational cost. As an appealing application of quantum computing, we show in this work that variational quantum algorithms can be integrated with machine learning (ML) techniques as a promising scheme for exploring accurate PESs. Different from using a ML model to represent the potential energy, we encode the molecular geometry information into a deep neural network (DNN) for representing parameters of the variational quantum eigensolver (VQE), leaving the PES to the wave function ansatz. Once the DNN model is trained, the variational optimization procedure that hinders the application of the VQE to complex systems is avoided and thus the evaluation of PESs is significantly accelerated. Numerical results demonstrate that a simple DNN model is able to reproduce accurate PESs for small molecules.

preprint2022arXiv

Ferromagnetic Negative Charge-Transfer Insulator: from Theoretical Proposal to Material Realization

Here we propose another type of ferromagnetic semiconductors: ferromagnetic negative charge-transfer insulator (FNCTI). In FNCTI, the negative charge-transfer states strongly enhance the ferromagnetic (FM) exchange interactions and the orbital hybridization gap permits the magnetic molecular orbitals as the underlying magnetic units rather than local atomic orbitals. Thus the FM exchange interactions are rather strong and decay slowly due to the large spearding of magnetic molecular orbitals. This is distinct from the superexchange mechanism where FM exchange interactions are quite weak as summarized in the well-known Goodenough-Kanamori-Anderson semi-empirical rules.Through first-principle calculations with the hybrid functional, PbO-type CrAs monolayer is mapped out to be a FNCTI, which possesses a band gap $\sim$ 0.35 eV, FM nearest-/next-nearest-neighbor exchange coupling strength $\sim$ 57/40 meV, and a high $T_c$ $\sim$ 1500 K respectively. It is believed that the existence of FNCTI validates the long-pending hypothesis by D. I. Khomskii and G. A. Sawatzky in 1997 [Solid State Commun. 102, 87 (1997)].

preprint2022arXiv

gDNA: Towards Generative Detailed Neural Avatars

To make 3D human avatars widely available, we must be able to generate a variety of 3D virtual humans with varied identities and shapes in arbitrary poses. This task is challenging due to the diversity of clothed body shapes, their complex articulations, and the resulting rich, yet stochastic geometric detail in clothing. Hence, current methods to represent 3D people do not provide a full generative model of people in clothing. In this paper, we propose a novel method that learns to generate detailed 3D shapes of people in a variety of garments with corresponding skinning weights. Specifically, we devise a multi-subject forward skinning module that is learned from only a few posed, un-rigged scans per subject. To capture the stochastic nature of high-frequency details in garments, we leverage an adversarial loss formulation that encourages the model to capture the underlying statistics. We provide empirical evidence that this leads to realistic generation of local details such as wrinkles. We show that our model is able to generate natural human avatars wearing diverse and detailed clothing. Furthermore, we show that our method can be used on the task of fitting human models to raw scans, outperforming the previous state-of-the-art.

preprint2022arXiv

ICON: Implicit Clothed humans Obtained from Normals

Current methods for learning realistic and animatable 3D clothed avatars need either posed 3D scans or 2D images with carefully controlled user poses. In contrast, our goal is to learn an avatar from only 2D images of people in unconstrained poses. Given a set of images, our method estimates a detailed 3D surface from each image and then combines these into an animatable avatar. Implicit functions are well suited to the first task, as they can capture details like hair and clothes. Current methods, however, are not robust to varied human poses and often produce 3D surfaces with broken or disembodied limbs, missing details, or non-human shapes. The problem is that these methods use global feature encoders that are sensitive to global pose. To address this, we propose ICON ("Implicit Clothed humans Obtained from Normals"), which, instead, uses local features. ICON has two main modules, both of which exploit the SMPL(-X) body model. First, ICON infers detailed clothed-human normals (front/back) conditioned on the SMPL(-X) normals. Second, a visibility-aware implicit surface regressor produces an iso-surface of a human occupancy field. Importantly, at inference time, a feedback loop alternates between refining the SMPL(-X) mesh using the inferred clothed normals and then refining the normals. Given multiple reconstructed frames of a subject in varied poses, we use SCANimate to produce an animatable avatar from them. Evaluation on the AGORA and CAPE datasets shows that ICON outperforms the state of the art in reconstruction, even with heavily limited training data. Additionally, it is much more robust to out-of-distribution samples, e.g., in-the-wild poses/images and out-of-frame cropping. ICON takes a step towards robust 3D clothed human reconstruction from in-the-wild images. This enables creating avatars directly from video with personalized and natural pose-dependent cloth deformation.

preprint2022arXiv

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. We describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

preprint2022arXiv

Neural Point-based Shape Modeling of Humans in Challenging Clothing

Parametric 3D body models like SMPL only represent minimally-clothed people and are hard to extend to clothing because they have a fixed mesh topology and resolution. To address these limitations, recent work uses implicit surfaces or point clouds to model clothed bodies. While not limited by topology, such methods still struggle to model clothing that deviates significantly from the body, such as skirts and dresses. This is because they rely on the body to canonicalize the clothed surface by reposing it to a reference shape. Unfortunately, this process is poorly defined when clothing is far from the body. Additionally, they use linear blend skinning to pose the body and the skinning weights are tied to the underlying body parts. In contrast, we model the clothing deformation in a local coordinate space without canonicalization. We also relax the skinning weights to let multiple body parts influence the surface. Specifically, we extend point-based methods with a coarse stage, that replaces canonicalization with a learned pose-independent "coarse shape" that can capture the rough surface geometry of clothing like skirts. We then refine this using a network that infers the linear blend skinning weights and pose dependent displacements from the coarse representation. The approach works well for garments that both conform to, and deviate from, the body. We demonstrate the usefulness of our approach by learning person-specific avatars from examples and then show how they can be animated in new poses and motions. We also show that the method can learn directly from raw scans with missing data, greatly simplifying the process of creating realistic avatars. Code is available for research purposes at {\small\url{https://qianlim.github.io/SkiRT}}.

preprint2022arXiv

Q$^2$Chemistry: A quantum computation platform for quantum chemistry

Quantum computer provides new opportunities for quantum chemistry. In this article, we present a versatile, extensible, and efficient software package, named Q$^2$Chemistry, for developing quantum algorithms and quantum inspired classical algorithms in the field of quantum chemistry. In Q$^2$Chemistry, wave function and Hamiltonian can be conveniently mapped into the qubit space, then quantum circuits can be generated according to a specific quantum algorithm already implemented in the package or newly developed by the users. The generated circuits can be dispatched to either a physical quantum computer, if available, or to the internal virtual quantum computer realized by simulating quantum circuit on classical supercomputers. As demonstrated by our benchmark simulations with up to 72 qubit, Q$^2$Chemistry achieves excellent performance in simulating medium scale quantum circuits. Application of Q$^2$Chemistry to simulate molecules and periodic systems are given with performance analysis.

preprint2022arXiv

Reducing circuit depth in adaptive variational quantum algorithms via effective Hamiltonian theories

Electronic structure simulation is an anticipated application for quantum computers. Due to high-dimensional quantum entanglement in strongly correlated systems, the quantum resources required to perform such simulations are far beyond the capacity of current quantum devices. To reduce the quantum circuit complexity, it has been suggested to incorporate a part of the electronic correlation into an effective Hamiltonian, which is often obtained from a similarity transformation of the electronic Hamiltonian. In this work, we introduce a new transformation in the form of a product of a linear combination of excitation operators to construct the effective Hamiltonian with finite terms. To demonstrate its accuracy, we also consider an equivalent adaptive variational algorithm with this transformation and show that it can obtain an accurate ground state wave function. The effective Hamiltonian defined with this new transformation is incorporated into the adaptive variational quantum algorithms to maintain constant-size quantum circuits. The new computational scheme is assessed by performing numerical simulations for small molecules. Chemical accuracy is achieved with a much shallower circuit depth.

preprint2022arXiv

Two-Dimensional Semiconducting Metal Organic Frameworks with Auxetic Effect, Room Temperature Ferrimagnetism, Chiral Ferroelectricity, Bipolar Spin Polarization and Topological Nodal Lines/Points

Two-dimensional (2D) semiconductors integrated with two or more functions are the cornerstone for constructing multifunctional nanodevices, but remain largely limited. Here, by tuning the spin state of organic linkers and the symmetry/topology of crystal lattice, we predict a class of unprecedented multifunctional semiconductors in 2D Cr(II) five-membered heterocyclic metal organic frameworks that simultaneously possess auxetic effect, room temperature ferrimagnetism, chiral ferroe-lectricity, electrically reversible spin polarization and topological nodal lines/points. Taking 2D Cr(TDZ)$_2$ (TDZ=1.2.5-thiadiazole) as an exemplification, the auxetic effect is produced by the anti-tetra-chiral lattice structure. The high temperature ferrimagnetism originates from the strong d-p direct magnetic exchange interaction between Cr cations and TDZ doublet radical anions. Meanwhile, the clockwise-counterclockwise alignment of TDZ' dipoles results in unique 2D chiral ferroelectricity with atomic-scale vortex-antivortex states. 2D Cr(TDZ)$_2$ is an intrinsic bipolar magnetic semiconductor where half-metallic conduction with switchable spin-polarization direction can be induced by applying a gate voltage. Besides, the symmetry of the little group C$_4$ of lattice structure endows 2D Cr(TDZ)$_2$ with topological nodal lines and a quadratic nodal point in the Brillouin zone near the Fermi level.

preprint2021arXiv

Two-Dimensional Bipolar Magnetic Semiconductor with High Curie Temperature and Electrically Controllable Spin Polarization Realized in Exfoliated Cr(pyrazine)$_2$ Monolayer

Exploring two-dimensional (2D) magnetic semiconductors with room temperature magnetic ordering and electrically controllable spin polarization is a highly desirable but challenging task for nanospintronics. Here, through first principles calculations, we propose to realize such a material by exfoliating the recently synthesized organometallic layered crystal Li$_{0.7}$[Cr(pyz)$_2$]Cl$_{0.7}$0.25$\cdot$(THF) (pyz = pyrazine, THF = tetrahydrofuran) [Science 370, 587 (2020)]. The feasibility of exfoliation is confirmed by the rather low exfoliation energy of 0.27 J/m$^2$, even smaller than that of graphite. In exfoliated Cr(pyz)$_2$ monolayer, each pyrazine ring grabs one electron from the Cr atom to become a radical anion, then a strong $d$-$p$ direct exchange magnetic interaction emerges between Cr cations and pyrazine radicals, resulting in room temperature ferrimagnetism with a Curie temperature of 342 K. Moreover, Cr(pyz)$_2$ monolayer is revealed to be an intrinsic bipolar magnetic semiconductor where electrical doping can induce half-metallic conduction with controllable spin-polarization direction.

preprint2020arXiv

Can High-Temperature Reactions Be Described by a Minimum Energy Path Model? Steric Hindrance Matters

High-temperature reactions widely exist in nature. However, they are difficult to be characterized either experimentally or computationally. The routinely used minimum energy path (MEP) model in computational modeling of chemical reactions is not justified to describe high-temperature reactions since high-energy structures are actively involved there. In this study, using CH4 decomposition on the Cu(111) surface as an example, we systematically compare MEP results with those obtained by explicitly sampling all relevant structures via ab initio molecular dynamics (AIMD) simulations at different temperatures. Interestingly, we find that, for reactions protected by a strong steric hindrance effect, the MEP is still effectively followed even at a temperature close to the Cu melting point. In contrast, without such a protection, the flexibility of surface Cu atoms can lead to a significant free energy barrier reduction at a high temperature. Accordingly, some conclusions about graphene growth mechanisms based on MEP calculations should be revisited. Physical insights provided by this study can deepen our understanding on high-temperature surface reactions.

preprint2020arXiv

Doping dependence of electronic structure of infinite-layer NdNiO2

We investigate the electronic structure of nickelate superconductor NdNiO2 upon hole doping, by means of density-functional theory and dynamical mean-field theory. We demonstrate the strong intrinsic hybridization between strongly correlated states formed by Ni-3dx2-y2 orbital and itinerant electrons due to Nd-5d and Ni-3dz2 orbitals, producing a valence-fluctuating correlated metal as the normal state of hole-doped NdNiO2. The Hund's rule appears to play a dominating role on multi-orbital physics in the lightly doped compound, while its effect is gradually reduced by increasing the doping level. Crucially, the hole-doping leads to intricate effects on Ni-3d orbitals, such as a non-monotonic change of electron occupation in lightly doped level, and a flipping orbital configuration in the overdoped regime. Additionaly, we also map out the topology of Fermi surface at different doping levels. These findings render a preferred window to peek into electron pairing and superconductivity.

preprint2020arXiv

Extreme-Scale Density Functional Theory High Performance Computing of DGDFT for Tens of Thousands of Atoms using Millions of Cores on Sunway TaihuLight

High performance computing (HPC) is a powerful tool to accelerate the Kohn-Sham density functional theory (KS-DFT) calculations on modern heterogeneous supercomputers. Here, we describe a massively extreme-scale parallel and portable implementation of discontinuous Galerkin density functional theory (DGDFT) method on the Sunway TaihuLight supercomputer. The DGDFT method uses the adaptive local basis (ALB) functions generated on-the-fly during the self-consistent field (SCF) iteration to solve the KS equations with the high precision comparable to that of plane-wave basis set. In particular, the DGDFT method adopts a two-level parallelization strategy that makes use of different types of data distribution, task scheduling, and data communication schemes, and combines with the feature of master-slave multi-thread heterogeneous parallelism of SW26010 processor, resulting in extreme-scale HPC KS-DFT calculations on the Sunway TaihuLight supercomputer. We show that the DGDFT method can scale up to 8,519,680 processing cores (131,072 core groups) on the Sunway TaihuLight supercomputer for investigating the electronic structures of two-dimensional (2D) metallic graphene systems containing tens of thousands of carbon atoms.

preprint2020arXiv

Interpolative separable density fitting decomposition for accelerating Hartree-Fock exchange calculations within numerical atomic orbitals

The high cost associated with the evaluation of Hartree-Fock exchange (HFX) makes hybrid functionals computationally challenging for large systems. In this work, we present an efficient way to accelerate HFX calculations with numerical atomic basis sets. Our approach is based on the recently proposed interpolative separable density fitting (ISDF) decomposition to construct a low rank approximation of HFX matrix, which avoids explicit calculations of the electron repulsion integrals (ERIs) and significantly reduces the computational cost. We implement the ISDF method for hybrid functional (PBE0) calculations in the HONPAS package. We take benzene and polycyclic aromatic hydrocarbons molecules as examples and demonstrate that hybrid functionals with ISDF yields quite promising results at a significantly reduced computational cost. Especially, the ISDF approach reduces the total cost for evaluating HFX matrix by nearly 2 orders of magnitude compared to conventional approaches of direct evaluation of ERIs.

preprint2020arXiv

Learning to Dress 3D People in Generative Clothing

Three-dimensional human body models are widely used in the analysis of human pose and motion. Existing models, however, are learned from minimally-clothed 3D scans and thus do not generalize to the complexity of dressed people in common images and videos. Additionally, current models lack the expressive power needed to represent the complex non-linear geometry of pose-dependent clothing shapes. To address this, we learn a generative 3D mesh model of clothed people from 3D scans with varying pose and clothing. Specifically, we train a conditional Mesh-VAE-GAN to learn the clothing deformation from the SMPL body model, making clothing an additional term in SMPL. Our model is conditioned on both pose and clothing type, giving the ability to draw samples of clothing to dress different body shapes in a variety of styles and poses. To preserve wrinkle detail, our Mesh-VAE-GAN extends patchwise discriminators to 3D meshes. Our model, named CAPE, represents global shape and fine local structure, effectively extending the SMPL body model to clothing. To our knowledge, this is the first generative model that directly dresses 3D human body meshes and generalizes to different poses. The model, code and data are available for research purposes at https://cape.is.tue.mpg.de.

preprint2020arXiv

Simulating periodic systems on quantum computer

The variational quantum eigensolver (VQE) is one of the most appealing quantum algorithms to simulate electronic structure properties of molecules on near-term noisy intermediate-scale quantum devices. In this work, we generalize the VQE algorithm for simulating extended systems. However, the numerical study of an one-dimensional (1D) infinite hydrogen chain using existing VQE algorithms shows a remarkable deviation of the ground state energy with respect to the exact full configuration interaction (FCI) result. Here, we present two schemes to improve the accuracy of quantum simulations for extended systems. The first one is a modified VQE algorithm, which introduces an unitary transformation of Hartree-Fock orbitals to avoid the complex Hamiltonian. The second one is a Post-VQE approach combining VQE with the quantum subspace expansion approach (VQE/QSE). Numerical benchmark calculations demonstrate that both of two schemes provide an accurate enough description of the potential energy curve of the 1D hydrogen chain. In addition, excited states computed with the VQE/QSE approach also agree very well with FCI results.

preprint2020arXiv

The dynamic parallel distribution algorithm for hybrid density-functional calculations in HONPAS package

This work presents a dynamic parallel distribution scheme for the Hartree-Fock exchange~(HFX) calculations based on the real-space NAO2GTO framework. The most time-consuming electron repulsion integrals~(ERIs) calculation is perfectly load-balanced with 2-level master-worker dynamic parallel scheme, the density matrix and the HFX matrix are both stored in the sparse format, the network communication time is minimized via only communicating the index of the batched ERIs and the final sparse matrix form of the HFX matrix. The performance of this dynamic scalable distributed algorithm has been demonstrated by several examples of large scale hybrid density-functional calculations on Tianhe-2 supercomputers, including both molecular and solid states systems with multiple dimensions, and illustrates good scalability.

preprint2020arXiv

The Moving-Grid Effect in the Harmonic Vibrational Frequency Calculations with Numeric Atom-Centered Orbitals

When using atom-centered integration grids, the portion of the grid that belongs to a certain atom also moves when this atom is displaced. In the paper, we investigate the moving-grid effect in the calculation of the harmonic vibrational frequencies when using all-electron full-potential numeric atomic-centered orbitals as the basis set. We find that, unlike the first order derivative (i.e., forces), the moving-grid effect plays an essential role for the second order derivatives (i.e., vibrational frequencies). Further analysis reveals that predominantly diagonal force constant terms are affected, which can be bypassed efficiently by invoking translational symmetry. Our approaches have been demonstrated in both finite (molecules) and extended (periodic) systems.

preprint2020arXiv

The static parallel distribution algorithms for hybrid density-functional calculations in HONPAS package

Hybrid density-functional calculation is one of the most commonly adopted electronic structure theory used in computational chemistry and materials science because of its balance between accuracy and computational cost. Recently, we have developed a novel scheme called NAO2GTO to achieve linear scaling (Order-N) calculations for hybrid density-functionals. In our scheme, the most time-consuming step is the calculation of the electron repulsion integrals (ERIs) part. So how to create an even distribution of these ERIs in parallel implementation is an issue of particular importance. Here, we present two static scalable distributed algorithms for the ERIs computation. Firstly, the ERIs are distributed over ERIs shell pairs. Secondly, the ERIs is distributed over ERIs shell quartets. In both algorithms, the calculation of ERIs is independent of each other, so the communication time is minimized. We show our speedup results to demonstrate the performance of these static parallel distributed algorithms in the Hefei Order-N packages for \textit{ab initio} simulations (HONPAS).

preprint2019arXiv

Electronic and Magnetic Structure of Infinite-layer $\textrm{NdNiO}_2$: Trace of Antiferromagnetic Metal

The recent discovery of Sr-doped infinite-layer nickelate $\textrm{NdNiO}_2$ [D. Li et al. Nature 572, 624 (2019)] offers an exciting platform for investigating unconventional superconductivity in nickelatebased compounds. In this work, we present a first-principles calculations for the electronic and magnetic properties of undoped parent $\textrm{NdNiO}_2$. Intriguingly, we found that: 1) the paramagnetic phase has complex Fermi pockets with 3D characters near the Fermi level; 2) by including electronelectron interactions, 3d-electrons of Ni tend to form $(π, π, π)$ antiferromagnetic ordering at low temperatures; 3) with moderate interaction strength, 5d-electrons of Nd contribute small Fermi pockets that could weaken the magnetic order akin to the self-doping effect. Our results provide a plausible interpretation for the experimentally observed resistivity minimum and Hall coefficient drop. Moreover, we elucidate that antiferromagnetic ordering in $\textrm{NdNiO}_2$ is relatively weak, arising from the small exchange coupling between 3d-electrons of Niand also hybridization with 5d-electrons of Nd.

preprint2016arXiv

A many-body GW+BSE investigation of electronic and optical properties of C2N

A newly synthesized layered material C2N was investigated based on many- body perturbation theory using GW plus Bethe-Salpeter equation approach. The electronic band gap was determined to be ranging from 3.75 to 1.89 eV from monolayer to bulk. Significant GW quasiparticle corrections, of more than 0.9 eV, to the Kohn-Sham band gaps from the local density approximation (LDA) calculations are found. Strong excitonic effects play a crucial role in optical properties. We found large binding energies of greater than 0.6 eV for bound excitons in few-layer C2N, while it is only 0.04 eV in bulk C2N. All the structures exhibit strong and broad optical absorption in the visible light region, which makes C2N a promising candidate for solar energy conversion, such as photocatalytic water splitting.

preprint2016arXiv

Two-Dimensional Boron Oxides with Dirac Loop and Strongly anisotropic Carrier Mobility

Recently, two-dimensional boron sheets have attracted a lot of attentions owing to their structural polymorphs and outstanding properties. And, due to chemical complexity and electron deficiency of B atoms, the 2D boron sheets are easy affected by the environment. So, exploring novel 2D boron oxides gets highly needed. In this study, we theoretically explored 2D boron oxides structures and their basic properties. We found 2D boron oxides are metals or semimetals, when oxygen concentration is low. More interesting, the B6O1 exhibits Dirac Loop near the Fermi level and the maximum Fermi velocity is estimated as high as 0.61*10E6 m/s, which much close to that in graphene. In addition, when the oxygen concentration is one forth, the most stable B4O1 get a semiconductor with a direct band gap of 1.24 eV and a strong anisotropic carrier mobility. Such huge differences of carrier mobility between two directions have not been reported in other 2D materials. Besides, when oxygen concentra-tion is higher than one forth, the 2D boron sheets will get unstable based our results. We also proposed a method how to synthesize these systems. Our work provide a basic knowledge of 2D boron oxides and may lead to novel 2D boron oxides discovered. The unique electronic properties of the B6O1 and B4O1 render them promising 2D materials for applications in high-performance nanodevices.

preprint2016arXiv

Water on Silicene: Hydrogen Bond Autocatalysis Induced Physisorption-Chemisorption-Dissociation Transition

A single water molecule has nothing special. However, macroscopic water displays many anomalous properties at the interface, such as a high surface tension, hydrophobicity and hydrophillicity. Although the underlying mechanism is still elusive, hydrogen bond is expected to have played an important role. An interesting question is if the few-water molecule clusters will be qualitatively different from a single molecule. Using adsorption behavior as an example, by carefully choosing two-dimensional silicene as the substrate material, we demonstrate that water monomer, dimer and trimer show contrasting properties. The additional water molecules in dimer and trimer induce a transition from physisorption to chemisorption then to dissociation on silicene. Such a hydrogen bond autocatalytic effect is expected to have a broad application potential in silicene-based water molecule sensor and metal-free catalyst for water dissociation.

preprint2016arXiv

δ-phosphorene: a two dimensional material with high negative Poisson's ratio

As a basic mechanical parameter, Poisson's ratio (ν) measures the mechanical responses of solids against external loads. In rare cases, materials have a negative Poisson's ratio (NPR), and present an interesting auxetic effect. That is, when a material is stretched in one direction, it will expand in the perpendicular direction. To design modern nanoscale electromechanical devices with special functions, two dimensional (2D) auxetic materials are highly desirable. In this work, based on first principles calculations, we rediscover the previously proposed δ-phosphorene (δ-P) nanosheets [Jie Guan et al., Phys. Rev. Lett. 2014, 113, 046804] are good auxetic materials with a high NPR. The results show that the Young's modulus and Poisson's ratio of δ-P are all anisotropic. The NPR value along the grooved direction is up to -0.267, which is much higher than the recently reported 2D auxetic materials. The auxetic effect of δ-P originated from its puckered structure is robust and insensitive to the number of layers due to weak interlayer interactions. Moreover, δ-P possesses good flexibility because of its relatively small Young's modulus and high critical crack strain. If δ-P can be synthesized, these extraordinary properties would endow it great potential in designing low dimensional electromechanical devices.

preprint2015arXiv

Few-Layer C2N: A Promising Metal-free Photocatalyst for Water Splitting

Successful synthesis of the nitrogenated holey two-dimensional structures C2N (Nat. Commun. 2015, 6, 6486) using simply wet-chemical reaction offer a cost-effective way to generate other 2D materials with novel optical and electronic properties. On basis of the density functional theory calculations, we investigate electronic properties of monolayer and multilayer C2N. We find that few-layer C2N have a direct bandgap and the direct bandgap of the system can vary from 2.47 eV for monolayer to 1.84 eV for a five-layer. Besides, for the few-layer C2N, appropriate band gap, band edge alignments, and strong visible-light absorption demonstrate it may be a potential metal-free visible-light driven photocatalyst for water splitting.

preprint2015arXiv

First-principles study of two-dimensional van der Waals heterojunctions

Research on graphene and other two-dimensional (2D) materials, such as silicene, germanene, phosphorene, hexagonal boron nitride (h-BN), graphitic carbon nitride (g-C3N4), graphitic zinc oxide (g-ZnO) and molybdenum disulphide (MoS2), has recently received considerable interest owing to their outstanding properties and wide applications. Looking beyond this field, combining the electronic structures of 2D materials in ultrathin van der Waals heterojunctions has also emerged to widely study theoretically and experimentally to explore some new properties and potential applications beyond their single components. Here, this article reviews our recent theoretical studies on the structural, electronic, electrical and optical properties of 2D van der Waals heterojunctions using density functional theory calculations, including the Graphene/Silicene, Graphene/Phosphorene, Graphene/g-ZnO, Graphene/MoS2 and g-C3N4/MoS2 heterojunctions. Our theoretical simulations, designs and calculations show that novel 2D van der Waals heterojunctions provide a promising future for electronic, electrochemical, photovoltaic, photoresponsive and memory devices in the experiments.

preprint2015arXiv

Is phosphorene with intrinsic defect still an ideal anode material?

The diffusion of Li in electrode materials is key factor to charging/discharging rate capacity of Li-ion battery (LIB). Recently, two-dimensional phosphorene has been proposed as a very promising electrode material due to its ultrafast and directional lithium diffusion, as well as large energy capacity. Here, on the basis of density functional theory, we report that the intrinsic point defects, including vacancy and stone-wales defects, will block the directional ultrafast diffusion of lithium in phosphorene. On the defect-free phosphorene, diffusion of Li along the zigzag lattice direction is 1.6 billion times faster than along the armchair lattice direction, and 260 times faster than that in graphite. By introducing intrinsic vacancy and stone-wales defect, the diffusion energy barrier of Li along zigzag lattice direction increases sharply to the range of 0.17 ~ 0.49 eV, which block the ultrafast migration of lithium along the zigzag lattice direction. Meanwhile, the open circuit voltage increases with the emergence of defects, which is not suitable for anode materials. In addition, the formation energies of defects in phosphorene are much lower than those in graphene and silicene sheet, therefore, it is highly demanded to generate defect-free phosphorene for LIB applications.

preprint2015arXiv

Multilayer C2N: Effect of Stacking Order and Number of Layers on Bandgap and Its Controlled Electronic Properties by External Electric Field

Successful synthesis of the nitrogenated holey two-dimensional structures C2N (Nat. Commun. 2015, 6, 1-7) using simply wet-chemical reaction offer a cost-effective way to generate other 2D materials with novel optical and electronic properties. Using the few-layer C2N as models, we have performed an ab initio study of electronic properties of layered C2N. Band gaps of this system exhibit monotone decreasing as the number of layers increase. And a direct-gap to indirect-gap transition at the bulk C2N. Besides, when we apply an out-of-plane electric field on few-layer C2N, the band gap of multilayer C2N will be decreased as the electric field increased and a semiconductor-semimetal transition will happen for five-layer C2N under an appropriate electric field, whereas the band gap of monolayer C2N is unchanged under electric field. Owing to their tunable bandgaps in a wide range, layers C2N will have tremendous opportunities to be applied in nanoscale electronic and optoelectronic devices.

preprint2014arXiv

Bilayer Graphene Growth via a Penetration Mechanism

From both fundamental and technical points of view, a precise control of the layer number of graphene samples is very important. To reach this goal, atomic scale mechanisms of multilayer graphene growth on metal surfaces should be understood. Although it is a geometrically favorable pathway to transport carbon species to interface and then form a new graphene layer there, penetration of a graphene overlayer is not a chemically straightforward process. In this study, the possibility of different active species to penetrate a graphene overlayer on Cu(111) surface is investigated based on first principles calculations. It is found that carbon atom penetration can be realized via an atom exchange process, which leads to a new graphene growth mechanism. Based on this result, a bilayer graphene growth protocol is proposed to obtain high quality samples. Such a penetration possibility also provides a great flexibility for designed growth of graphene nanostructures.

preprint2014arXiv

Defect in Phosphorene

Defects are inevitably present in materials and always can affect their properties. Here, first-principles calculations are performed to systematically study the stability, structural and electronic properties of ten kinds of point defects in semiconducting phosphorene, including the Stone-Wales (SW-1 and SW-2) defect, single (SV59 and SV5566) and double vacancy (DV585-1, DV585-2, DV555777-1, DV555777-2, DV555777-3 and DV4104) defects. We find that these defects are all much easily created in phosphorene with higher areal density compared with graphene and silicene. They are easy distinguish each other and correlate with their defective atomic structures with simulated scanning tunneling microscopy images at positive bias. The SW, DV585-1, DV555777 and DV4104 defects have little effect on phosphorene's electronic properties and defective phosphorene monolayers still show semiconducting with similar band gap values to perfect phosphorene. The SV59 and DV585-2 defects can introduce unoccupied localized states into phosphorene's fundamental band gap. Specifically, the SV59 and 5566 defects can induce hole doping in phosphorene, and only the stable SV59 defect can result in local magnetic moments in phosphorene different from all other defects.

preprint2014arXiv

Electrical Control of Carriers Spin Orientation in FeVTiSi Heusler Alloy

The direct control of carriers spin by electric field under room temperature is one of the most important challenges in the field of spintronics. For this purpose, we here propose a quaternary Heusler alloy FeVTiSi. Based on first principles calculations, FeVTiSi alloy is found to be an intrinsic bipolar magnetic semiconductor in which the valence band and conduction band approach the Fermi level through opposite spin channels. Thus FeVTiSi alloy can conduct completely spin-polarized currents with tunable spin-polarization direction simply by applying a gate voltage. Furthermore, by Monte Carlo simulations based on the classical Heisenberg Hamiltonian, the Curie temperature of FeVTiSi alloy is predicted to be as high as 1293 K, far above the room temperature. The bipolar magnetic semiconducting character and the high Curie temperature endow the FeVTiSi alloy great potentials in developing electrically controllable spintronic devices working at room temperature.

preprint2014arXiv

Electronic Structure of Large-Scale Graphene Nanoflakes

With the help of the recently developed SIESTA-PEXSI method [J. Phys.: Condens. Matter \textbf{26}, 305503 (2014)], we perform Kohn-Sham density functional theory (DFT) calculations to study the stability and electronic structure of hexagonal graphene nanoflakes (GNFs) with up to 11,700 atoms. We find the electronic properties of GNFs, including their cohesive energy, HOMO-LUMO energy gap, edge states and aromaticity, depend sensitively on the type of edges (ACGNFs and ZZGNFs), size and the number of electrons. We observe that, due to the edge-induced strain effect in ACGNFs, large-scale ACGNFs' cohesive energy decreases as their size increases. This trend does not hold for ZZGNFs due to the presence of many edge states in ZZGNFs. We find that the energy gaps $E_g$ of GNFs all decay with respect to $1/L$, where $L$ is the size of the GNF, in a linear fashion. But as their size increases, ZZGNFs exhibit more localized edge states. We believe the presence of these states makes their gap decrease more rapidly. In particular, when $L$ is larger than 6.40 $nm$, we find that ZZGNFs exhibit metallic characteristics. Furthermore, we find that the aromatic structures of GNFs appear to depend only on whether the system has $4N$ or $4N+2$ electrons, where $N$ is an integer.

preprint2014arXiv

Electronic Structures of Hybrid Graphene/Phosphorene Nanocomposite

Combining the electronic structures of two-dimensional monolayers in ultrathin hybrid nanocomposites is expected to display new properties beyond their simplex components. Here, first-principles calculations are performed to study the structural, electronic and optical properties of hybrid graphene and phosphorene nanocomposite. It turns out that weak van der Waals interactions dominate between graphene and phosphorene with their intrinsic electronic properties preserved. Hybrid graphene and phosphorene nanocomposite shows tunable band gaps at graphene's Dirac point and a transition from hole doing to electron doing for graphene as the interfacial distance decreases. Charge transfer between graphene to phosphorene induces interfacial electron-hole pairs in hybrid graphene and phosphorene nanocomposite with enhanced visible light response.

preprint2014arXiv

First-Principles Study of Hybrid Graphene and MoS$_2$ Nanocomposites

Combining the electronic properties of graphene and molybdenum disulphide (MoS$_2$) monolayers in two-dimensional (2D) ultrathin hybrid nanocomposites have been synthesized experimentally to create excellent electronic, electrochemical, photovoltaic, photoresponsive and memory devices. Here, first-principles calculations are performed to investigate the electronic, electrical and optical properties in hybrid G/MoS$_2$ and G/MoS$_2$/G nanocomposites. It turns out that weak van der Waals interactions dominate between graphene and MoS$_2$ with their intrinsic electronic properties preserved. Interestingly, tunable p-type doping of graphene is very easy to achieve by applying electric fields perpendicular to hybrid G/MoS$_2$ and G/MoS$_2$/G nanocomposites, because electrons can easily transfer from the Dirac point of graphene to the conduction band of MoS$_2$ due to the work function of graphene close to the electronic affinity of MoS$_2$. Vertical electric fields can generate strong p-type but weak n-type doping of graphene, inducing electron-hole pairs in hybrid G/MoS$_2$/G sandwiched nanocomposites. Moreover, improved optical properties in hybrid G/MoS$_2$ and G/MoS$_2$/G nanocomposites are also expected with potential photovoltaic and photoresponsive applications.

preprint2014arXiv

First-principles Study of the Interactions of Electron Donor and Acceptor Molecules with Phosphorene

Density functional theory calculations have been carried out to investigate single-layer phosphorene functionalized with two kinds of organic molecules, i.e. an electrophilic molecule tetracyano-p-quinodimethane (TCNQ) as electron acceptor and a nucleophilic molecule tetrathia-fulvalene (TTF) as electron donor. The TCNQ molecule introduces shallow acceptor states in the gap of phosphorene close to the valence band edge (VBE), which makes the doped system a p-type semiconductor. However, when the TTF molecule is adsorbed on the phosphorene, the occupied molecular states introduced into the gap are of deep donor states so that effective n-doping for transport cannot be realized. This disadvantageous situation can be amended by applying an external electric field perpendicular to the phosphorene surface with direction from the phosphorene to the TTF molecule, under which the TTF-introduced donor states move closer to conduction band edge (CBE) of the phosphorene and then the TTF-doped phosphorene system becomes an n-type semiconductor. The effective bipolar doping of single-layer phosphorene via molecular adsorption predicted above, especially n-doping against its native p-doping propensity, would broaden the way to the application of this new type of two-dimensional material in nanoelectronic and optoelectronic devices.

preprint2014arXiv

Nano-scale displacement sensing based on Van der Waals interaction

We propose the nano-scale displacement sensor with high resolution for weak-force systems could be realized based on vertical stacked two-dimensional (2D) atomic corrugated layer materials bound through Van der Waals (VdW) interaction. Using first-principles calculations, we found the electronic structure of bi-layer blue phosphorus (BLBP) varies appreciably to both the lateral and vertical interlayer displacement. The variation of electronic structure due to the lateral displacement is attributed to the changing of the interlayer distance dz led by atomic layer corrugation, which is in a uniform picture with vertical displacement. Despite different stacking configurations, the change of in-direct band gap is proportional to dz-2. This stacking configuration independent dz-2 law is found also works for other graphene-like corrugated bi-layer materials, for example MoS2. By measuring the tunable electronic structure using absorption spectroscopy, the nano-scale displacement could be detected. BLBP represents a large family of bi-layer 2D atomic corrugated materials for which the electronic structure is sensitive to the interlayer vertical and lateral displacement, thus could be used for nano-scale displacement sensor. Since this kind of sensor is established on atomic layers coupled through VdW interaction, it provides unique applications in measurements of nano-scale displacement induced by tiny external force.

preprint2013arXiv

Control of Spin in La(Mn,Zn)AsO Alloy by Carrier Doping

The control of spin without magnetic field is one of challenges in developing spintronic devices. In an attempt to solve this problem, we proposed a novel hypothetic LaMn0.5Zn0.5AsO alloy from two experimentally synthesized rare earth element transition metal arsenide oxides, i.e. LaMnAsO and LaZnAsO. On the basis of the first-principles calculations with strong-correlated correction, we found that the LaMn0.5Zn0.5AsO alloy is an antiferromagnetic semiconductor at ground state, while bipolar magnetic semiconductor at ferromagnetic state. Both electron and hole doping in the LaMn0.5Zn0.5AsO alloy induces the transition from antiferromagnetic to ferromagnetic, as well as semiconductor to half metal. In particular, the spin-polarization direction is switchable depending on the doped carrier's type. As carrier doping can be realized easily in experiment by applying a gate voltage, the LaMn0.5Zn0.5AsO alloy stands for a promising spintronic material to generate and control the spin-polarized carriers with electric field.

preprint2013arXiv

Orientation Sensitive Nonlinear Growth of Graphene: A Geometry-determined Epitaxial Growth Mechanism

Although the corresponding carbon-metal interactions can be very different, a similar nonlinear growth behavior of graphene has been observed for different metal substrates. To understand this interesting experimental observation, a multiscale $\lq\lq$standing-on-the-front" kinetic Monte Carlo study is performed. An extraordinary robust geometry effect is identified, which solely determines the growth kinetics and makes the details of carbon-metal interaction not relevant at all. Based on such a geometry-determined mechanism, epitaxial growth behavior of graphene can be easily predicted in many cases. As an example, an orientation-sensitive growth kinetics of graphene on Ir(111) surface has been studied. Our results demonstrate that lattice mismatch pattern at the atomic level plays an important role for macroscopic epitaxial growth.

preprint2012arXiv

A First-principles Prediction of Two-Dimensional Superconductivity in Pristine B2C Single layer

Based on first-principles lattice dynamics and electron-phonon coupling calculations, B2C sheet is predicted to be a two-dimensional (2D) phonon-mediated superconductor with a relatively high transition temperature (Tc). The electron-phonon coupling parameter calculated is 0.92, and it is mainly contributed by low frequency out-of-plane phonon modes and electronic states with a π character. When the Coulomb pseudopotential is set to 0.10, the estimated temperature Tc is 19.2 K. To be best of our knowledge, B2C is the first pristine 2D superconductor with a Tc higher than the boiling point of liquid helium.

preprint2012arXiv

Bipolar Magnetic Semiconductors: A New Class of Spintronics Materials

Electrical control of spin polarization is very desirable in spintronics, since electric field can be easily applied locally in contrast with magnetic field. Here, we propose a new concept of bipolar magnetic semiconductor (BMS) in which completely spin-polarized currents with reversible spin polarization can be created and controlled simply by applying a gate voltage. This is a result of the unique electronic structure of BMS, where the valence and conduction bands possess opposite spin polarization when approaching the Fermi level. Our band structure and spin-polarized electronic transport calculations on semi-hydrogenated single-walled carbon nanotubes confirm the existence of BMS materials and demonstrate the electrical control of spin-polarization in them.

preprint2012arXiv

Hydrogenated Bilayer Wurtzite SiC Nanofilms: A Two-Dimensional Bipolar Magnetic Semiconductor Material

Recently, a new kind of spintronics materials, bipolar magnetic semiconductor (BMS), has been proposed. The spin polarization of BMS can be conveniently controlled by a gate voltage, which makes it very attractive in device engineering. Now, the main challenge is finding more BMS materials. In this article, we propose that hydrogenated wurtzite SiC nanofilm is a two-dimensional BMS material. Its BMS character is very robust under the effect of strain, substrate, or even a strong electric field. The proposed two-dimensional BMS material paves the way to use this promising new material in an integrated circuit.

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.

preprint2011arXiv

Atomistic Mechanisms of Nonlinear Graphene Growth on Ir Surface

As a two-dimensional material, graphene can be naturally obtained via epitaxial growth on a suitable substrate. Growth condition optimization usually requires an atomistic level understanding of the growth mechanism. In this article, we perform a mechanistic study about graphene growth on Ir(111) surface by combining first principles calculations and kinetic Monte Carlo (kMC) simulations. Small carbon clusters on the Ir surface are checked first. On terraces, arching chain configurations are favorable in energy and they are also of relatively high mobilities. At steps, some magic two-dimensional compact structures are identified, which show clear relevance to the nucleation process. Attachment of carbon species to a graphene edge is then studied. Due to the effect of substrate, at some edge sites, atomic carbon attachment becomes thermodynamically unfavorable. Graphene growth at these difficult sites has to proceed via cluster attachment, which is the growth rate determining step. Based on such an inhomogeneous growth picture, kMC simulations are made possible by successfully separating different timescales, and they well reproduce the experimentally observed nonlinear kinetics. Different growth rates and nonlinear behaviors are predicted for different graphene orientations, which is consistent with available experimental results. Importantly, as a phenomenon originated from lattice mismatch, inhomogeneity revealed in this case is expected to be quite universal and it should also make important roles in many other hetero-epitaxial systems.

preprint2011arXiv

Chemical Vapor Deposition Growth of Graphene using Other Hydrocarbon Sources

Graphene has attracted a lot of research interests due to its exotic properties and a wide spectrum of potential applications. Chemical vapor deposition (CVD) from gaseous hydrocarbon sources has shown great promises for large-scale graphene growth. However, high growth temperature, typically 1000°C, is required for such growth. Here we demonstrate a revised CVD route to grow graphene on Cu foils at low temperature, adopting solid and liquid hydrocarbon feedstocks. For solid PMMA and polystyrene precursors, centimeter-scale monolayer graphene films are synthesized at a growth temperature down to 400°C. When benzene is used as the hydrocarbon source, monolayer graphene flakes with excellent quality are achieved at a growth temperature as low as 300°C. The successful low-temperature growth can be qualitatively understood from the first principles calculations. Our work might pave a way to undemanding route for economical and convenient graphene growth.

preprint2011arXiv

CO2 dissociation activated through electron attachment on reduced rutile TiO2(110)-1x1 surface

Converting CO$_2$ to useful compounds through the solar photocatalytic reduction has been one of the most promising strategies for artificial carbon recycling. The highly relevant photocatalytic substrate for CO$_2$ conversion has been the popular TiO$_2$ surfaces. However, the lack of accurate fundamental parameters that determine the CO$_2$ reduction on TiO$_2$ has limited our ability to control these complicated photocatalysis processes. We have systematically studied the reduction of CO2 at specific sites of the rutile TiO$_2$(110)-1x1 surface using scanning tunneling microscopy at 80 K. The dissociation of CO2 molecules is found to be activated by one electron attachment process and its energy threshold, corresponding to the CO$_2^{\dot-}$/CO$_2$ redox potential, is unambiguously determined to be 2.3 eV higher than the onset of the TiO$_2$ conduction band. The dissociation rate as a function of electron injection energy is also provided. Such information can be used as practical guidelines for the design of effective catalysts for CO$_2$ photoreduction.

preprint2011arXiv

Diamondization of Graphene and Graphene-BN Bilayers: Chemical Functionalization and Electronic Structure Engineering

In this article, based on first-principles calculations, we systematically study functionalization induced diamonization of graphene bilayer and graphene-BN hybrid bilayer. With single-side functionalization, the diamondized structures are magnetic semiconductor. Interestingly, if both sides of the bilayer are functionalized, diamondization becomes spontaneous without a barrier. On the other hand, when the bottom layer of the bilayer graphene is replaced by a single hexagonal BN layer, the diamondized structure becomes nonmagnetic metal. The tunable electronic and magnetic properties pave new avenues to construct graphene-based electronics and spintronics devices.

preprint2011arXiv

Evidence of Photocatalytic Dissociation of Water on TiO2 with Atomic Resolution

Photocatalytic water splitting reaction on TiO2 surface is one of the fundamental issues that bears significant implication in hydrogen energy technology and has been extensively studied. However, the existence of the very first reaction step, the direct photo-dissociation of water, has been disregarded. Here, we provide unambiguously experimental evidence to demonstrate that adsorbed water molecules on reduced rutile TiO2(110)-1\times1 surface can be dissociated under UV irradiation using low temperature scanning tunneling microscopy. It is identified that a water molecule at fivefold coordinated Ti (Ti5c) site can be photocatalytically dissociated, resulting in a hydroxyl at Ti5c and another hydroxyl at bridge oxygen row. Our findings reveal a missing link in the photocatalytic water splitting reaction chain, which greatly contribute to the detailed understanding of underlying mechanism.

preprint2011arXiv

First-Principles Thermodynamics of Graphene Growth on Cu Surface

Chemical vapor deposition (CVD) is an important method to synthesis grapheme on a substract. Recently, Cu becomes the most popular CVD substrate for graphene growth. Here, we combine electronic structure calculation, molecular dynamics simulation, and thermodynamics analysis to study the graphene growth process on Cu surface. As a fundamentally important but previously overlooked fact, we find that carbon atoms are thermodynamically unfavorable on Cu surface under typical experimental conditions. The active species for graphene growth should thus mainly be CHx instead of atomic carbon. Based on this new picture, the nucleation behavior can be understood, which explains many experimental observations and also provides us a guide to improve graphene sample quality.

preprint2011arXiv

Oxygen molecule dissociation on carbon nanostructures with different types of nitrogen doping

Energy barrier of oxygen molecule dissociation on carbon nanotube or graphene with different types of nitrogen doping is investigated using density functional theory. The results show that the energy barriers can be reduced efficiently by all types of nitrogen doping in both carbon nanotubes and graphene. Graphite-like nitrogen and Stone-Wales defect nitrogen decrease the energy barrier more efficiently than pyridine-like nitrogen, and a dissociation barrier lower than 0.2 eV can be obtained. Higher nitrogen concentration reduces the energy barrier much more efficiently for graphite-like nitrogen. These observations are closely related to partial occupation of π* orbitals and change of work functions. Our results thus provide useful insights into the oxygen reduction reactions.

preprint2011arXiv

Quantum Size Effect and Electronic Stability of Freestanding Metal Atom Wires

Using DFT calculations, we present a thorough study of the quantum size effects on the stability of freestanding metal atom wires. Our systems include Na, Ag, Au, In, Ga and Pb atom wires, i.e. $s$, $sd$, and $sp$ electron prototypes. We found that the total energy always oscillates with the wire length, which clearly indicates the existence of preferred lengths. Increasing the length, the s-system exhibits even-odd oscillations following a $a/x +b/x^2$ decay law in the stability, which can be attributed to electron band filling and quantum confinement along the wire. The $sd$-system exhibits a similar oscillation pattern, even in the presence of $sd$ hybridization. In $sp$-system, the energy oscillations are beyond the simple even-odd nature, likely due to unpaired p orbitals and the corresponding nontrival band filling. Our findings clearly demonstrate that electronic contribution is quite critical to the stability of freestanding wires, and this stability may be important even when wires are deposited on substrates or strained. This study sheds light on the underlying formation mechanism of metal atom wires.

preprint2011arXiv

Reversible Transition Between Thermodynamically Stable Phases with Low Density of Oxygen Vacancies on SrTiO$_3$(110) Surface

The surface reconstruction of SrTiO$_3$(110) is studied with scanning tunneling microscopy and density functional theory (DFT) calculations. The reversible phase transition between (4$\times$1) and (5$\times$1) is controlled by adjusting the surface metal concentration [Sr] or [Ti]. Resolving the atomic structures of the surface, DFT calculations verify that the phase stability changes upon the chemical potential of Sr or Ti. Particularly, the density of oxygen vacancies is low on the thermodynamically stabilized SrTiO$_3$(110) surface.

preprint2011arXiv

Structure of graphene oxide: thermodynamics versus kinetics

Graphene oxide (GO) is an important intermediate to prepare graphene and it is also a versatile material with various applications. However, despite its importance, the detailed structure of GO is still unclear. For example, previous theoretical studies based on energetics have suggested that hydroxyl chain is an important structural motif of GO, which, however, is found to be contrary to nuclear magnetic resonance (NMR) experiment. In this study, we check both thermodynamic and kinetic aspects missed previously. First principles thermodynamics gives a free energy based stability ordering similar to that based on energetics, and hydroxyl chain is thus thermodynamically still favorable. At the same time, by checking the calculated vibrational frequencies, we note that hydroxyl chain structure is also inconsistent with infrared experiment. Therefore, kinetics during GO synthesis is expected to make an important role in GO structure. Transition state calculations predict large energy barriers between local minima, which suggests that experimentally obtained GO has a kinetically constrained structure.

preprint2011arXiv

Tunable band gap of graphane nanoribbons under uniaxial elastic strain: a first-principles study

In this Letter, we investigate the strain-induced band-gap modulation of both armchair and zigzag graphane nanoribbons based on the first-principles calculations. Within the elastic range, the band gap changes linearly with the uniaxial strain, where the band-gap of graphane nanoribbonsis more sensitive to compressive than tensile deformation. The band-gap deformation mainly originates from the shift of valence band edge of graphane nanoribbons under stain. Our results imply the great potential of graphane nanoribbons in the pressure sensor and optical electronics applications.

preprint2010arXiv

First Principles NMR Signatures of Graphene Oxide

Nuclear magnetic resonance (NMR) has been widely used in the graphene oxide (GO) structure study. However, the detailed relationship between its spectroscopic features and the GO structural configuration has remained elusive. Based on first principles $^{13}$C chemical shift calculations using the gauge including projector augmented waves (GIPAW) method, we provide a spectrum-structure connection. Chemical shift of carbon is found to be very sensitive to atomic environment, even with an identical oxidation group. Factors determining the chemical shifts for epoxy and hydroxy groups have been discussed. GO structures previously reported in the literature have been checked from the NMR point of view. The energetically favorable hydroxy chain structure is not expected to be widely existed in real GO samples according to our NMR simulations. The epoxy pair we proposed previously is also supported by chemical shift calculations.

preprint2010arXiv

Nearly Free Electron States in Graphane Nanoribbon Superlattice

Nearly free electron (NFE) state has been widely studied in low dimensional systems. Based on first-principles calculations, we identify two types of NFE states in graphane nanoribbon superlattice, similar to those of graphene nanoribbons and boron nitride nanoribbons. Effect of electron doping on the NFE states in graphane nanoribbon superlattice has been studied, and it is possible to open a vacuum transport channel via electron doping.

preprint2010arXiv

Nearly Free Electron States in Graphene Nanoribbon Superlattices

Nearly free electron (NFE) state is an important kind of unoccupied state in low dimensional systems. Although it is intensively studied, a clear picture on its physical origin and its response behavior to external perturbations is still not available. Our systematic first-principles study based on graphene nanoribbon superlattices suggests that there are actually two kinds of NFE states, which can be understood by a simple Kronig-Penney potential model. An atom-scattering-free NFE transport channel can be obtained via electron doping, which may be used as a conceptually new field effect transistor.

preprint2010arXiv

Rotation and dissociation dynamics of a single O2 molecule on the Pt(111) surface determined from a first principles study

The STM induced rotation and dissociation dynamics of a single oxygen molecule on the Pt(111) surface have been finally determined by first principles calculations together with a newly developed statistical model for inelastic electron tunneling. Several long-standing puzzles associated with these dynamic processes in this classic system have been fully resolved. It is found that the unexpected low energy barrier of the O2 rotation is originated from an ingenious pathway, while the prior occupation of the metastable hcp-hollow site after the O2 dissociation can be attributed to a dynamic process of surface accommodation. The experimentally observed non-integer power-law dependence of the rotation rate as a function of the current can be perfectly explained by taking into account the randomness of multi-electron inelastic tunneling processes.

preprint2009arXiv

Universal $\sqrt{2}\times\sqrt{2}$ structure and short-range charge order at the surfaces of BaFe$_{2-x}$Co$_{x}$As$_{2}$ compounds with various Co doping levels

The structure and electronic order at the cleaved (001) surfaces of the newly-discovered pnictide superconductors BaFe$_{2-x}$Co$_{x}$As$_{2}$ with x ranging from 0 to 0.32 are systematically investigated by scanning tunneling microscopy. A $\sqrt{2}\times\sqrt{2}$ surface structure is revealed for all the compounds, and is identified to be Ba layer with half Ba atoms lifted-off by combination with theoretical simulation. A universal short-range charge order is observed at this $\sqrt{2}\times\sqrt{2}$ surface associated with an energy gap of about 30 meV for all the compounds.

preprint2006arXiv

Electronic Structure of Bilayer Graphene: A Real-space Green's Function Study

In this paper, a real-space analytical expression for the free Green's function (propagator) of bilayer graphene is derived based on the effective-mass approximation. Green's function displays highly spatial anisotropy with three-fold rotational symmetry. The calculated local density of states (LDOS) of a perfect bilayer graphene produces the main features of the observed scanning tunneling microscopy (STM) images of graphite at low bias voltage. Some predicted features of the LDOS can be verified by STM measurements. In addition, we also calculate the LDOS of bilayer graphene with vacancies by using the multiple-scattering theory (scatterings are localized around the vacancy of bilayer graphene). We observe that the interference patterns are determined mainly by the intrinsic properties of the propagator and the symmetry of the vacancies.