Co-reporter:Xiao-Qing Liang, Xiao-Jiao Deng, Sheng-Jie Lu, Xiao-Ming Huang, Ji-Jun Zhao, Hong-Guang Xu, Wei-Jun Zheng, and Xiao Cheng Zeng
The Journal of Physical Chemistry C March 30, 2017 Volume 121(Issue 12) pp:7037-7037
Publication Date(Web):March 3, 2017
DOI:10.1021/acs.jpcc.7b00943
We present a joint experimental and theoretical study on double iron atom doped germanium clusters, Fe2Gen–/0 (n = 3–12). The experimental photoelectron spectra of cluster anions are reasonably reproduced by theoretical simulations. The low-lying structures of the iron-doped semiconductor clusters are obtained by using an ab initio computation-based genetic-algorithm global optimization method. We find that the smaller-sized Fe2Gen– (n = 3–8) clusters adopt bipyramid-based geometries, while the larger ones (n ≥ 9) adopt polyhedral cagelike structures with one interior Fe atom. Interestingly, starting from n = 8, the most stable anionic clusters Fe2Gen– exhibit structures that are different from that of their neutral counterparts Fe2Gen. Robust ferromagnetic interaction is found between the two doped iron atoms in the neutral clusters Fe2Gen, while the total spin moment always remains at 4 μB for all the neutral double iron atom doped germanium clusters up to n = 12. This behavior is in stark contrast to the magnetic quenching behavior typically observed in germanium clusters doped with a single Fe atom.
Co-reporter:Si Zhou, Nanshu Liu, Zhiyu Wang, and Jijun Zhao
ACS Applied Materials & Interfaces July 12, 2017 Volume 9(Issue 27) pp:22578-22578
Publication Date(Web):June 16, 2017
DOI:10.1021/acsami.7b05755
Composites of transition metal and carbon-based materials are promising bifunctional catalysts for the oxygen reduction reaction (ORR) and oxygen evolution reaction (OER), and are widely used in rechargeable metal–air batteries. However, the mechanism of their enhanced bicatalytic activities remains elusive. Herein, we construct N-doped graphene supported by Co(111) and Fe(110) substrates as bifunctional catalysts for ORR and OER in alkaline media. First-principles calculations show that these heterostructures possess a large number of active sites for ORR and OER with overpotentials comparable to those of noble metal benchmark catalysts. The catalytic activity is modulated by the coupling strength between graphene and the metal substrates, as well as the charge distribution in the graphitic sheet, which is delicately mediated by N dopants. These theoretical results uncover the key parameters that govern the bicatalytic properties of hybrid materials and help prescribe the principles for designing multifunctional electrocatalysts of high performance.Keywords: graphene; overpotential; oxygen evolution; oxygen reduction; rechargeable metal−air battery; transition metal;
Co-reporter:Yuanchang Li;Yue Qi;Nannan Han;Zhepeng Zhang;Xiebo Zhou;Bing Deng;Qiucheng Li;Mengxi Liu;Zhongfan Liu;Yanfeng Zhang
ACS Nano February 28, 2017 Volume 11(Issue 2) pp:1807-1815
Publication Date(Web):January 21, 2017
DOI:10.1021/acsnano.6b07773
Hetero-epitaxial growth of hexagonal boron nitride (h-BN) from the edges of graphene domains or vice versa has been widely observed during synthesis of in-plane heterostructures of h-BN-G on Rh(111), Ir(111), and even Cu foil. We report that on a strongly coupled Re(0001) substrate via a similar two-step sequential growth strategy, h-BN preferably nucleated on the edges of Re(0001) steps rather than on the edges of existing graphene domains. Statistically, one-third of the domain boundaries of graphene and h-BN were patched seamlessly, and the others were characterized by obvious “defect lines” when the total coverage approached a full monolayer. This imperfect merging behavior can be explained by translational misalignment and lattice mismatch of the resulting separated component domains. According to density functional theory calculations, this coexisting patching and non-patching growth behavior was radically mediated by the strong adlayer–substrate (A–S) interactions, as well as the disparate formation energies of the attachment of B–N pairs or B–N lines along the edges of the Re(0001) steps versus the graphene domains. This work will be of fundamental significance for the controllable synthesis of in-plane heterostructures constructed from two-dimensional layered materials with consideration of A–S interactions.Keywords: edges of Re steps and graphene domains; graphene and hexagonal boron nitride heterostructures; preferable nucleation; ultra-high-vacuum scanning tunneling microscopy/spectroscopy;
Co-reporter:Nannan Han, Nan Gao, and Jijun Zhao
The Journal of Physical Chemistry C August 24, 2017 Volume 121(Issue 33) pp:17893-17893
Publication Date(Web):August 2, 2017
DOI:10.1021/acs.jpcc.7b04209
Blue phosphorene (blue P), a new two-dimensional allotrope of phosphorus, has attracted great attentions due to its high carrier mobility, suitable band gap and appreciable stability comparable to black phosphorene (black P). Motivated by recent experimental success in synthesizing monolayer blue P on Au(111) surface, here we investigate the nucleation mechanism and growth behavior of PN clusters on Au(111) substrate by ab initio calculations. During the initial growth stage, PN clusters transform from dispersed atoms to zigzag chain at N = 4, and further turn to ring-based one-dimensional chain at N = 11. This peculiar behavior is ascribed to the competition between the interaction among P atoms and the attraction on P atoms by Au substrate. On the basis of the interaction energies between black/blue phosphorene and Au(111) surface, we further propose that monolayer blue P can be synthesized on a chemical active metal substrate, while a relatively inert metal substrate would be beneficial to the growth of monolayer black P. Our theoretical findings offer experimentalists insightful guidance to fabricate monolayer blue P or black P by choosing appropriate substrate.
Co-reporter:Shaohong Liu;Zhiyu Wang;Si Zhou;Fengjiao Yu;Mengzhou Yu;Chang-Yang Chiang;Wuzong Zhou;Jieshan Qiu
Advanced Materials 2017 Volume 29(Issue 31) pp:
Publication Date(Web):2017/08/01
DOI:10.1002/adma.201700874
The oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) are cornerstone reactions for many renewable energy technologies. Developing cheap yet durable substitutes of precious-metal catalysts, especially the bifunctional electrocatalysts with high activity for both ORR and OER reactions and their streamlined coupling process, are highly desirable to reduce the processing cost and complexity of renewable energy systems. Here, a facile strategy is reported for synthesizing double-shelled hybrid nanocages with outer shells of Co-N-doped graphitic carbon (Co-NGC) and inner shells of N-doped microporous carbon (NC) by templating against core–shell metal–organic frameworks. The double-shelled NC@Co-NGC nanocages well integrate the high activity of Co-NGC shells into the robust NC hollow framework with enhanced diffusion kinetics, exhibiting superior electrocatalytic properties to Pt and RuO2 as a bifunctional electrocatalyst for ORR and OER, and hold a promise as efficient air electrode catalysts in Zn–air batteries. First-principles calculations reveal that the high catalytic activities of Co-NGC shells are due to the synergistic electron transfer and redistribution between the Co nanoparticles, the graphitic carbon, and the doped N species. Strong yet favorable adsorption of an OOH* intermediate on the high density of uncoordinated hollow-site C atoms with respect to the Co lattice in the Co-NGC structure is a vital rate-determining step to achieve excellent bifunctional electrocatalytic activity.
Co-reporter:Linwei Sai;Xue Wu;Xiaoming Huang
Journal of Cluster Science 2017 Volume 28( Issue 3) pp:1729-1737
Publication Date(Web):15 February 2017
DOI:10.1007/s10876-017-1181-5
Using our improved genetic algorithm combined with density functional theory calculations, we perform unbiased global search for the most stable structures of Sin clusters with n = 21–29. A new principal direction crossover operation is introduced to search the prolate structures on the potential energy surface. We confirmed previous results for most sized Sin clusters and find more stable structure at n = 23 and 26. In particular, we find a new endohedral cage structure at Si23, which is 0.4 eV lower in energy than previously reported prolate structure and exhibit a distinct peak on the second order difference of total cluster energy. We discuss the competition between prolate and spherical motifs in the medium-sized Sin clusters in light of binding energy and HOMO–LUMO gap. The transition of structural motif from prolate to near spherical in Sin clusters occurs at n = 26. The outer cages are composed of pentagon, hexagon, quadrangle and heptagon, with one to three interior atoms filled inside the cage. These results establish a more complete picture for the structural evolution of the medium-sized silicon clusters. Average bond lengths and Mayer bond order are also discussed.
Co-reporter:Pengbo Zhang, Jianhua Ding, Dan Sun, Jijun Zhao
Journal of Nuclear Materials 2017 Volume 492(Volume 492) pp:
Publication Date(Web):15 August 2017
DOI:10.1016/j.jnucmat.2017.05.022
•Interstitial noble gas atoms tend to stay together with each other by self-trapping.•H/He prefers to locate interstitial sites nearby Ne atom than other interstitials.•Noble gas atoms like Ne can act as a trapping site for H/He impurities in bcc Fe.We investigate the energetics and clustering trend of noble gas atoms (He, Ne, and Ar) in bcc Fe, and their interactions with vacancy or H/He impurities using first-principles calculations. We determine the formation energy of single and double noble gas atoms inside Fe host lattice as well as the resulted volume changes. The Ne/Ar formation energy is two and three times that of He. The attraction between Ne/Ar and vacancy is stronger than He-vacancy, indicating higher dissolution energy of Ne/Ar. The interstitial Ne-Ne/Ar-Ar pairs have stronger attractions (−1.91 eV/−1.40 eV) than He-He (−0.37 eV), forming stable <110> configurations. Such strong attraction means that He/Ne/Ar tend to aggregate, which can be well explained by the lower electron density induced by interstitial noble gas atoms and its strong repulsion with Fe atoms. Moreover, H/He energetically prefers to occupy the tetrahedral sites nearby Ne/Ar atom. The attraction energies of He-Ne/He-Ar pairs (−1.01 eV/-0.85 eV) are much stronger than those of H-Ne/H-Ar (−0.22 eV/−0.10 eV) and their charge density differences are discussed. The distinct attraction strengths by various noble gas atoms provide a preliminary explanation for the difference in irradiation effects on Fe solid by He, Ne, Ar, and He+H/Ne+He. These findings improve our understanding about the behavior of noble gas atoms and gas bubble formation in iron under irradiation.
Co-reporter:Si Zhou, Nanshu Liu, Jijun Zhao
Computational Materials Science 2017 Volume 130(Volume 130) pp:
Publication Date(Web):1 April 2017
DOI:10.1016/j.commatsci.2017.01.009
•Monolayer black and blue phosphorus quantum dots possess highly tunable energy gap.•Their frontier molecular orbitals straddle the water redox potentials.•They can harvest the full spectrum of visible light.•They are promising candidates of visible-light photocatalysts for water splitting.The discovery and synthesis of phosphorene stimulate the exploration of other two-dimensional (2D) phosphorus allotropes as well as phosphorus nanostructures. The recent successful preparation of phosphorene quantum dots further broadens the application areas of this novel 2D material. Herein we propose, for the first time, to sculpture monolayer black and blue phosphorus into quantum dots for photocatalytic solar water splitting. The size dependency of electronic, optical and photocatalytic properties of monolayer phosphorus quantum dots are systematically investigated. According to our first-principles calculations, these nanostructures possess highly tunable energy gap and work function; their frontier molecular orbitals straddle the water redox potentials with reasonable oxidizing and reducing power; they can harvest the full spectrum of visible light. Therefore, monolayer phosphorus quantum dots are promising candidates of visible-light photocatalysts for water splitting.Download high-res image (228KB)Download full-size image
Co-reporter:Linwei Sai;Xue Wu;Nan Gao;R. Bruce King
Nanoscale (2009-Present) 2017 vol. 9(Issue 37) pp:13905-13909
Publication Date(Web):2017/09/28
DOI:10.1039/C7NR02399E
Using a genetic algorithm combined with density functional theory calculations, we perform a global search for the lowest-energy structures of Bn clusters with n = 46, 48, 50. Competition among different structural motifs including a hollow cage, core–shell, bilayer, and quasi-planar, is investigated. For B46, a core–shell B4@B42 structure resembling the larger Bn clusters with n ≥ 68 is found to compete with a quasi-planar structure with a central hexagonal hole. A quasi-planar configuration with two connected hexagonal holes is most favorable for B50. More interestingly, an unprecedented bilayer structure is unveiled at B48, which can be extended to a two-dimensional bilayer phase exhibiting appreciable stability. Our results suggest alternatives to the cage motif as lower-energy Bn cluster structures with n > 50.
Co-reporter:Nannan Han;Hongsheng Liu;Junfeng Zhang;Junfeng Gao
Nanoscale (2009-Present) 2017 vol. 9(Issue 10) pp:3585-3592
Publication Date(Web):2017/03/09
DOI:10.1039/C6NR09962A
The in-plane combination of graphene (G) and hexagonal-boron nitride (h-BN) leads to lateral h-BN/G heterostructures, which are promising candidates for novel two-dimensional electronics. The quality of the interface between G and h-BN domains is crucial for the device performance. By comprehensive first-principles calculations, we explore the heteroepitaxial growth of graphene along the edge of an h-BN domain on a Cu(111) surface and compare it with that on a Cu(111) terrace. We find that the graphene nucleation site strongly depends on the chemical potential of carbon and predeposited h-BN coverage. Under the suitable carbon concentration and coverage of h-BN, graphene mainly grows along the h-BN edge, leading to a sharp and straight h-BN/G interface. Our results provide insightful knowledge to synthesize well-defined h-BN/G and other lateral heterostructures.
Co-reporter:Zhifeng Liu;Xiaojuan Liu
Nanoscale (2009-Present) 2017 vol. 9(Issue 47) pp:18781-18787
Publication Date(Web):2017/12/07
DOI:10.1039/C7NR06431D
Superhalogens, which have larger electron affinity than any halogen, play an important role in physical chemistry and materials design because of their applications in hydrogen storage and lithium-ion batteries. Inspired by the unique geometries and electronic properties of II–VI/III–V cage clusters, particularly the experimentally synthesized B12N12, we propose a core–shell structure model to design new superhalogens. The idea is assessed by conducting ab initio calculations on endohedral cage clusters X@B12N12 (X = F, Cl, Br) and other similar systems. With an exceptionally large electron affinity of 5.36 eV, the stable F@B12N12 cluster behaves as a novel superhalogen that can serve as a building block for Li salts and hyperhalogens. The findings highlight a new route for the discovery of superhalogens and provide useful building blocks for the bottom-up design of materials.
Co-reporter:Heng-Fu Lin;Li-Min Liu
Journal of Materials Chemistry C 2017 vol. 5(Issue 9) pp:2291-2300
Publication Date(Web):2017/03/02
DOI:10.1039/C7TC00013H
A recently discovered layered semiconducting material, namely, black phosphorus, has a bandgap that depends on the number of layers. It is thus feasible to create lateral heterostructures using the same material with different thicknesses. The structural and electronic properties of the lateral heterostructures of bilayer and monolayer phosphorene (BP/MP) are investigated by first-principles calculations. Hydrogen passivated heterostructures have much lower formation energy when compared to unpassivated heterostructures or phosphorene grain boundaries. The electronic band structures of the heterostructures with and without hydrogen passivation are greatly dependent on interface orientation and exhibit three distinct families of characteristics: direct bandgap semiconductor, indirect bandgap semiconductor and metal. Additionally, a type-I to type-II band alignment transition takes place when the ribbon widths of the BP and MP regions decrease, enabling continuous modulation of the band offset. The magnitude of the band offset can also be effectively tuned by changing the stacking type of BP and the interface orientation. These theoretical findings would be helpful in the design and optimization of BP/MP heterostructures for electronic and optoelectronic applications.
Co-reporter:Kai Cheng;Yu Guo;Nannan Han;Yan Su;Junfeng Zhang
Journal of Materials Chemistry C 2017 vol. 5(Issue 15) pp:3788-3795
Publication Date(Web):2017/04/13
DOI:10.1039/C7TC00595D
Lateral semiconductor/semiconductor heterostructures made up of two-dimensional (2D) monolayer or few-layer materials provide new opportunities for 2D devices. Herein, we propose four lateral heterostructures constructed by phosphorene-like monolayer group-IV monochalcogenides, including GeS/GeSe, SnS/GeSe, SnSe/GeS and GeS/SnS. Using first-principles calculations, we investigated the energetics and electronic properties of these lateral heterostructures. The band structures and formation energies from supercell calculations demonstrate that these heterostructures retain semiconducting behavior and can be easily synthesized in the laboratory. The band offsets of monolayer, bilayer and trilayer heterojunctions at the Anderson limit are calculated from the valence/conduction band edges with respect to the vacuum energy level for each individual component. Among them, some heterostructures belong to type II band alignment and are promising for a high-efficiency solar cell.
Co-reporter:Wei Gan, Nannan Han, Chao Yang, Peng Wu, Qin Liu, Wen Zhu, Shuangming Chen, Chuanqiang Wu, Muhammad Habib, Yuan Sang, Zahir Muhammad, Jijun Zhao, and Li Song
ACS Nano 2017 Volume 11(Issue 2) pp:
Publication Date(Web):January 13, 2017
DOI:10.1021/acsnano.6b06144
Here we demonstrate a ternary Cu2NiZn alloy substrate for controllably synthesizing monolayer graphene using a liquid carbon precursor cyclohexane via a facile CVD route. In contrast with elemental metal or bimetal substrates, the alloy-induced synergistic effects that provide an ideal metallic platform for much easier dehydrogenation of hydrocarbon molecules, more reasonable strength of adsorption energy of carbon monomer on surface and lower formation energies of carbon chains, largely renders the success growth of monolayer graphene with higher electrical mobility and lower defects. The growth mechanism is systemically investigated by our DFT calculations. This study provides a selective route for realizing high-quality graphene monolayer via a scalable synthetic method by using economic liquid carbon supplies and multialloy metal substrates.Keywords: alloy; chemical vapor deposition; density functional theory; graphene; growth mechanism;
Co-reporter:Nan Gao, Hongsheng Liu, Si Zhou, Yizhen Bai, and Jijun Zhao
The Journal of Physical Chemistry C 2017 Volume 121(Issue 9) pp:
Publication Date(Web):February 17, 2017
DOI:10.1021/acs.jpcc.7b00023
Beyond graphene, other group IV monolayers with honeycomb lattice, including silicene, germanene, and stanene, have attracted much attention due to their peculiar physical properties and potential applications in future electronic devices. However, since sp3 hybridization is more favorable than sp2 hybridization for Si, Ge, and Sn, these group IV monolayers have to be stabilized by metal surfaces during epitaxial synthesis. Using systematical first-principles calculations, here we investigate the interactions between these monolayers and various metal surfaces, i.e., Ag(111), Ir(111), Pt(111), Al(111), Au(111), and Cu(111). STM images, charge density difference, and partial density of states of these monolayer/metal systems have been calculated and discussed. In combination with the known experimental facts, we find that a moderate strength of interaction at 0.6–0.7 eV/atom is beneficial for the epitaxial growth of silicene and germanene without too much buckling or in-plane distortion. We further propose that the Al(111) substrate might be a good choice for synthesis of stanene with low-buckled structure.
Co-reporter:Jijun Zhao, Hongsheng Liu, Zhiming Yu, Ruge Quhe, Si Zhou, Yangyang Wang, Cheng Cheng Liu, Hongxia Zhong, Nannan Han, Jing Lu, Yugui Yao, Kehui Wu
Progress in Materials Science 2016 Volume 83() pp:24-151
Publication Date(Web):October 2016
DOI:10.1016/j.pmatsci.2016.04.001
Silicene, a silicon analogue of graphene, has attracted increasing attention during the past few years. As early as in 1994, the possibility of stage corrugation in the Si analogs of graphite had already been theoretically explored. But there were very few studies on silicene until 2009, when silicene with a low buckled structure was confirmed to be dynamically stable by ab initio calculations. In spite of the low buckled geometry, silicene shares most of the outstanding electronic properties of planar graphene (e.g., the “Dirac cone”, high Fermi velocity and carrier mobility). Compared with graphene, silicene has several prominent advantages: (1) a much stronger spin–orbit coupling, which may lead to a realization of quantum spin Hall effect in the experimentally accessible temperature, (2) a better tunability of the band gap, which is necessary for an effective field effect transistor (FET) operating at room temperature, (3) an easier valley polarization and more suitability for valleytronics study. From 2012, monolayer silicene sheets of different superstructures were successfully synthesized on various substrates, including Ag(1 1 1), Ir(1 1 1), ZrB2(0 0 0 1), ZrC(1 1 1) and MoS2 surfaces. Multilayer silicene sheets have also been grown on Ag(1 1 1) surface. The experimental successes have stimulated many efforts to explore the intrinsic properties as well as potential device applications of silicene, including quantum spin Hall effect, quantum anomalous Hall effect, quantum valley Hall effect, superconductivity, band engineering, magnetism, thermoelectric effect, gas sensor, tunneling FET, spin filter, and spin FET, etc. Recently, a silicene FET has been fabricated, which shows the expected ambipolar Dirac charge transport and paves the way towards silicene-based nanoelectronics. This comprehensive review covers all the important theoretical and experimental advances on silicene to date, from the basic theory of intrinsic properties, experimental synthesis and characterization, modulation of physical properties by modifications, and finally to device explorations.
Co-reporter:Junfeng Zhang, Weiyu Xie, Xiaohong Xu, Shengbai Zhang, and Jijun Zhao
Chemistry of Materials 2016 Volume 28(Issue 14) pp:5022
Publication Date(Web):July 5, 2016
DOI:10.1021/acs.chemmater.6b01764
The in-plane heterostructures composed of graphene and hexagonal boron nitride (G/BN), as the first kind of two-dimensional metal/semiconductor heterostructures of one-atom thickness, are attractive for both fundamental low-dimensional physics and nanoscale devices because of the tailorable electronic properties. The atomic structures and electronic properties of interfaces in lateral G/BN heterostructures are investigated by first-principles calculations. The symmetric armchair interfaces have a similar formation energy but a larger band gap compared with the nonsymmetric interfaces. G/BN heterostructures with zigzag-type interfaces constructed under the guide of Clar’s rule are found to possess a lower formation energy than those with abrupt interfaces and open a finite band gap. In addition to the zigzag and armchair interfaces, other misorientated interfaces with pentagon and heptagon rings are also stable with low formation energies of 4.4–6.8 eV/nm. These theoretical results are important to clarify the correlation between atomic structures and electronic properties of in-plane G/BN heterostructures and establish a fundamental picture for further theoretical studies and device design.
Co-reporter:Si Zhou and Jijun Zhao
Nanoscale 2016 vol. 8(Issue 16) pp:8910-8918
Publication Date(Web):30 Mar 2016
DOI:10.1039/C5NR08810K
Graphene, a superior 2D material with high carrier mobility, has limited application in electronic devices due to zero band gap. In this regard, boron and nitrogen atoms have been integrated into the graphene lattice to fabricate 2D semiconducting heterostructures. It is an intriguing question whether oxygen can, as a replacement of nitrogen, enter the sp2 honeycomb lattice and form stable B–C–O monolayer structures. Here we explore the atomic structures, energetic and thermodynamic stability, and electronic properties of various 2D B–C–O alloys using first-principles calculations. Our results show that oxygen can be stably incorporated into the graphene lattice by bonding with boron. The B and O species favor forming alternate patterns into the chain- or ring-like structures embedded in the pristine graphene regions. These B–C–O hybrid sheets can be either metals or semiconductors depending on the B:O ratio. The semiconducting (B2O)nCm and (B6O3)nCm phases exist under the B- and O-rich conditions, and possess a tunable band gap of 1.0–3.8 eV and high carrier mobility, retaining ∼1000 cm2 V−1 s−1 even for half coverage of B and O atoms. These B–C–O alloys form a new class of 2D materials that are promising candidates for high-speed electronic devices.
Co-reporter:Jinxuan Liu, Wencai Zhou, Jianxi Liu, Yamato Fujimori, Tomohiro Higashino, Hiroshi Imahori, Xue Jiang, Jijun Zhao, Tsuneaki Sakurai, Yusuke Hattori, Wakana Matsuda, Shu Seki, Suresh Kumar Garlapati, Subho Dasgupta, Engelbert Redel, Licheng Sun and Christof Wöll
Journal of Materials Chemistry A 2016 vol. 4(Issue 33) pp:12739-12747
Publication Date(Web):13 Jul 2016
DOI:10.1039/C6TA04898F
We demonstrate the fabrication of a new class of epitaxial porphyrin metal–organic framework thin films whose photophysical properties can be tuned by the introduction of electron-donating diphenylamine (DPA) groups into the porphyrin skeleton. The attachment of DPA groups results in strongly improved absorption characteristics, yielding the highest photocarrier generation efficiency reported so far. DFT calculations identify a shift of the charge localization pattern in the VBM (lowest unoccupied molecular orbital), confirming that the introduction of the DPA groups is the main reason for the shift of the optical absorption spectrum and the improved photocurrent generation.
Co-reporter:Si Zhou, Yu Guo and Jijun Zhao
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 15) pp:10607-10615
Publication Date(Web):17 Mar 2016
DOI:10.1039/C6CP01012A
The thermoelectric properties of two-dimensional (2D) materials are of great interest for both fundamental science and device applications. Graphene oxide (GO), whose physical properties are highly tailorable by chemical and structural modifications, is a potential 2D thermoelectric material. In this report, we pattern nanoroads on GO sheets with epoxide functionalization, and investigate their ballistic thermoelectric transport properties based on density functional theory and the nonequilibrium Green's function method. These graphene oxide nanoroads (GONRDs) are all semiconductors with their band gaps tunable by the road width, edge orientation, and the structure of the GO matrix. These nanostructures show appreciable electrical conductance at certain doping levels and enhanced thermopower of 127–287 μV K−1, yielding a power factor 4–22 times of the graphene value; meanwhile, the lattice thermal conductance is remarkably reduced to 15–22% of the graphene value; consequently, attaining the figure of merit of 0.05–0.75. Our theoretical results are not only helpful for understanding the thermoelectric properties of graphene and its derivatives, but also would guide the theoretical design and experimental fabrication of graphene-based thermoelectric devices of high performance.
Co-reporter:Tengfei Cao, Da Wang, Dong-Sheng Geng, Li-Min Liu and Jijun Zhao
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 10) pp:7156-7162
Publication Date(Web):27 Jan 2016
DOI:10.1039/C5CP06528C
The indirect bandgap character of silicon greatly limits its applications in electronic or optoelectronic devices, and direct bandgaps are highly desirable in all silicon allotropes. The successful synthesis of ultrathin or even monolayer silicon films experimentally has opened new opportunities to further modulate the electronic structure of silicon through external modulation. In this work, strain or electric field effects on the electronic structure of ultrathin silicon film (USF) are systematically explored. The results demonstrate that all USFs are indirect band-gap semiconductors; interestingly, tensile strain or electric field efficiently tunes the USFs into direct band gap semiconductors. The indirect to direct band gap transition in the USFs not only extends their light adsorption spectra into the visible light region but also greatly enhances the adsorption intensity. Because of this, strained USFs have great potential to be used as a high-performance photovoltaic material. Furthermore, the high stability, moderate band-gap and proper band edge positions demonstrate that monolayer and bilayer USFs can also be used as photocatalysts for water splitting.
Co-reporter:Xiaoqing Li, Stephan Schönecker, Jijun Zhao, Börje Johansson, Levente Vitos
Journal of Alloys and Compounds 2016 Volume 676() pp:565-574
Publication Date(Web):15 August 2016
DOI:10.1016/j.jallcom.2016.03.218
•Alloying effect on ideal tensile strength of bcc ferro-/paramagnetic Fe uncovered.•Fe and all considered Fe-based random alloys are intrinsically brittle.•Magnetism limits the applicability of an established ideal tensile strength model.•All solutes change weakly ferromagnetic Fe towards a stronger ferromagnetic behavior.•All considered solutes are expected to increase the hardenability of ferritic steel.Using 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) Fe1−xMx (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 Fe-alloys and briefly discuss the relation between hardenability and the ITS.
Co-reporter:Ruihuan Li, Pengbo Zhang, Xiaojie Li, Jianhua Ding, Yuanyuan Wang, Jijun Zhao, Levente Vitos
Computational Materials Science 2016 Volume 123() pp:85-92
Publication Date(Web):October 2016
DOI:10.1016/j.commatsci.2016.06.019
Density functional theory calculations have been performed to study the effects of alloying Cr and W on the stability and diffusivity of interstitial He impurity in body-centered cubic (bcc) Fe host lattice. The interaction between two close Cr/W atoms is repulsive. The relative stable position for an interstitial He remains the tetrahedral interstitial site in the presence Cr. Energetically, He prefers to locate far away from W inside Fe host lattice due to the strong repulsive interaction between He and W. On the other hand, the He migration barrier becomes lower in the presence of Cr (0.026 eV) and W (0.049 eV), as compared to 0.059 eV for pure Fe. Addition of Cr is benefit for He self-trapping, while W is against. The effective diffusivity of He decreases with increasing Cr and W concentrations. Moreover, the additions of Cr and W slightly hinder He being trapped by monovacancy.
Co-reporter:Lizhao Liu, Junfeng Zhang, Haili Gao, Lu Wang, Xue Jiang, Jijun Zhao
Computational Materials Science 2016 Volume 112(Part B) pp:527-546
Publication Date(Web):1 February 2016
DOI:10.1016/j.commatsci.2015.06.032
Graphene, a two-dimensional monolayer of carbon atoms in honeycomb structure, is a research hotspot in multidisciplinary due to its excellent physical properties. To further extend the applications of graphene, various strategies have been proposed to tailor its physical properties. Recently, our group has carried out systematically computational studies on modifying graphene, including hydrogenation, oxidation, and introduction of grain boundaries. Both the hydrogenation and oxidation will convert sp2 hybridized carbons into sp3 configurations, while formation of grain boundaries only makes the sp2 carbon bonds distorted. Employing density functional theory calculations, structures, physical properties and applications of these modified graphene were explored, such as structural phase diagram, mechanical and electronic properties, and photocatalytic applications. It turns out that many physical properties of graphene are tunable, endowing graphene promising applications in various fields. In this review article, we will generally summarize our recent works on the hydrogenated graphene, graphene oxide, and graphene grain boundaries.
Co-reporter:Xiaojie Li, Jie Fu, Ying Qin, Shengzhi Hao, Jijun Zhao
Computational Materials Science 2016 Volume 112(Part A) pp:75-79
Publication Date(Web):1 February 2016
DOI:10.1016/j.commatsci.2015.10.014
•Gupta potentials are developed for five HCP rare earth elements (Er, Dy, Gd, Tb and Lu).•Cohesive energy, elastic constants, lattice constants, phonon dispersion are reproduced.•Surface energies and formation energies of defects are computed.Empirical Gupta-type potentials based on the second-moment approximation of tight-binding model have been developed for five hexagonal-closed-packed (hcp) rare earth metals: Er, Dy, Gd, Tb, Lu. The potentials can reproduce experimental cohesive energies, lattice constants and elastic constants of Er, Dy, Gd, Tb and Lu. The calculated bulk moduli, shear moduli, sound velocities and Debye temperatures are in reasonable agreement with the measured values. Vacancy formation energies and surface energies of these metals are also calculated and compared with previous theoretical results. Our potentials are able to simulate the phonon dispersions of Tb, Er, Dy, Gd, and Lu solids and reproduce Tb’s experimental data in detail. The present Gupta potentials would be useful for future simulations of elementary metals of Er, Dy, Gd, Tb, Lu and their alloys.
Co-reporter:Tengfei Cao, Xibo Li, Limin Liu, Jijun Zhao
Computational Materials Science 2016 Volume 112(Part A) pp:297-303
Publication Date(Web):1 February 2016
DOI:10.1016/j.commatsci.2015.10.042
•The band gap of monolayer Black Phosphorous (M-BP) is only reduced by vertical electric field due to its puckered configuration.•The electron effective mass (EEM) in M-BP is anisotropic and is tunable under the electric field.•Upon applying in-plane strain, the band gap of M-BP is almost linearly modulated, with tensile (compression) strain enlarging (reducing) the band gap.•In-plane strain also effectively modulate the response of M-BP to the vertical electric field.Electric field or in-plane strain is used to tailor the electronic structures of monolayer Black Phosphorus (M-BP). Upon applying electric field, the band gap of M-BP is greatly reduced and insulator–metal transition happens under certain field intensity. The electric field impact on the electron effective mass (EEM) of M-BP is anisotropic. The EEM along armchair direction is increased and that in the zigzag direction is greatly reduced. Tensile strain under small magnitude enlarges the band gap of M-BP and starts to reduce it when the strain becomes relatively large. The anisotropic EEM in the M-BP can also be reversed by the tensile strain. Under tensile strain, the electronic structure of M-BP becomes to be more efficiently modulated by the electric field. Compression strain only reduces the band gap of M-BP and has little impact on the EEM. For the M-BP under compression strain, its electronic structure can hardly be altered by the electric field.
Co-reporter:Si Zhou and Jijun Zhao
The Journal of Physical Chemistry C 2016 Volume 120(Issue 38) pp:21691-21698
Publication Date(Web):September 1, 2016
DOI:10.1021/acs.jpcc.6b07651
Co-reporter:Zhaoyang Zheng, Yi-Yang Sun, Weiyu Xie, Jijun Zhao, and Shengbai Zhang
The Journal of Physical Chemistry C 2016 Volume 120(Issue 47) pp:
Publication Date(Web):November 7, 2016
DOI:10.1021/acs.jpcc.6b11150
Electronic applications require the ability to dope a material with a controllable amount of impurities. However, current understanding of the doping mechanism in colloidal–synthesized quantum dots (QDs) is still limited. This is in contrast with bulk semiconductors for which first-principles-based theories have been well established. Using prototype CdSe as an example, here we propose an atomistic theory for the doping of colloidal-synthesized QDs. The key in our theory is the evaluation of atomic chemical potential inside the solution, whose range can deviate considerably from the bulk value due to the presence of solvent. This theory, coupled to first-principles calculations and ab initio molecular dynamics, is able to explain the difference of doping limit in Mn (or Co)-doped CdSe QDs and their bulk counterparts. It also explains the doping behavior of a number of other 3d transition-metal impurities in CdSe QDs in contrast with the solid case.
Co-reporter:Yi Du;Jincheng Zhuang;Hongsheng Liu;Zhi Li;Jiaou Wang;Xun Xu;Haifeng Feng;Lan Chen;Kehui Wu;Xiaolin Wang;Shi Xue Dou
Science Advances 2016 Volume 2(Issue 7) pp:e1600067
Publication Date(Web):22 Jul 2016
DOI:10.1126/sciadv.1600067
Quasi-freestanding silicene with massless Dirac fermion characteristics has been successfully obtained by oxygen intercalation.
Co-reporter:Yingying Huang;Yan Su;Lu Wang;Xiao Cheng Zeng;Chongqin Zhu;Xiaoxiao Cao;Xue Jiang;Sheng Meng
Science Advances 2016 Volume 2(Issue 2) pp:e1501010
Publication Date(Web):12 Feb 2016
DOI:10.1126/sciadv.1501010
Researchers predict a new ice clathrate structure as the most stable ice polymorph with the lowest density in a negative-pressure region.
Co-reporter:Xue Jiang, Peng Wang and Jijun Zhao
Journal of Materials Chemistry A 2015 vol. 3(Issue 15) pp:7750-7758
Publication Date(Web):14 Jan 2015
DOI:10.1039/C4TA03438D
Since the graphene boom, great efforts have been devoted to two-dimensional (2D) monolayer materials with exciting possibilities for applications. Most known 2D materials are inorganic. Using the covalent triazine framework (CTF) as a representative, we explored 2D organic semiconductors using first-principles calculations. From a systematic study of the electronic band structures, work functions, CBM/VBM positions, and optical absorption spectra, we identified the CTF as a new class of 2D visible-light-driven organocatalyst for water splitting. Controllable construction of such CTFs from suitable organic subunits paves the way to correlate band alignment and chemical composites. In addition, multilayer CTFs have enhanced visible-light absorption compared to monolayer CTFs due to interlayer coupling. Our theoretical prediction not only has fulfilled the search for organic counterparts of inorganic photocatalysts for water splitting, but also would motivate scientists to further search for novel 2D organic materials with other technological applications.
Co-reporter:Jijun Zhao, Xiaoming Huang, Ruili Shi, Hongsheng Liu, Yan Su and R. Bruce King
Nanoscale 2015 vol. 7(Issue 37) pp:15086-15090
Publication Date(Web):25 Aug 2015
DOI:10.1039/C5NR04034E
Our ab initio global searches reveal the lowest-energy cage for B28, which is built from two B12 units and prevails over the competing structural isomers such as planar, bowl, and tube. This smallest boron cage extends the scope of all-boron fullerene and provides a new structural motif of boron clusters and nanostructures.
Co-reporter:Xue Jiang and Jijun Zhao
RSC Advances 2015 vol. 5(Issue 59) pp:48012-48023
Publication Date(Web):11 May 2015
DOI:10.1039/C5RA05852J
We investigated the high-pressure structures and properties of iron tetraborides (FeB4) using a combination of an ab initio high-throughput search and a particle-swarm optimization algorithm for crystal structure prediction. We found that, under compression, the boron sublattice in FeB4 from the buckled boron layer first polymerizes into B4 tetrahedral clusters and then forms cubo-octahedral B12 clusters. At 55 GPa, the orthorhombic crystal structure with a Pnnm space group (58-FeB4) transforms into a tetragonal I41/acd structure (142-FeB4), which is stable within a wide pressure range up to 695 GPa. Then, a cubic Imm phase (229-FeB4) emerges at higher pressures up to at least 1 TPa. The computed Vicker's hardnesses of 58-, 142-, and 229-FeB4 are 61.58, 47.44, and 50.87 GPa, respectively. All of them can be considered as superhard materials. Compared to the previously reported 58-FeB4 as a superhard superconductor, the B4 tetrahedral cluster-based 142-FeB4 is a superhard semiconductor with an indirect band gap of 1.34 eV. The pressure-induced metal-to-semiconductor transition can be related to a unique Fe–B–B three-center covalent bond. Moreover, 229-FeB4, which is composed of cubo-octahedral B12 clusters, is ferromagnetic with a magnetic moment of 0.929μB per Fe atom at ambient pressure. The magnetic moment will decrease rapidly with increasing pressure and be completely quenched as pressure exceeds 40 GPa. The pressure-induced evolution of boron cluster units not only adds new features to boron chemistry, but also gives rise to novel superhard semiconductors or ferromagnetic materials. Moreover, our results may inspire further experimental and theoretical interest in designing new materials using clusters as pseudo-atoms with expected properties.
Co-reporter:Hongsheng Liu, Nannan Han and Jijun Zhao
RSC Advances 2015 vol. 5(Issue 23) pp:17572-17581
Publication Date(Web):26 Jan 2015
DOI:10.1039/C4RA17320A
Monolayer transition metal dichalcogenides (TMDs) stand out in two-dimensional (2D) materials due to their potential applications in future microelectronic and optoelectronic devices. In experiments, field effect transistors (FET) based on the MoS2 monolayer are sensitive to environmental gases, especially O2. Thus, the oxidation of monolayer TMDs is a critical concern. By first-principles calculations, we reveal that a perfect single-layer sheet of TMDs stays intact when exposed in O2 due to the weak physical adsorption of O2. However, O2 can be chemically adsorbed onto the monolayer of TMDs (including MoS2, MoSe2, MoTe2, WS2, WSe2, and WTe2) with single vacancies of chalcogen, which are the most common defects in realistic TMD materials. The adsorption configurations and dissociation behavior of the O2 molecule at vacancy sites, as well as the possible diffusion behavior of oxygen adatoms on the TMD monolayer surface were explored. Oxidation significantly influenced the electronic properties of a defective MoS2 monolayer, while other defective TMD monolayers (especially MoTe2 and WTe2) suffered less from oxidation. Our theoretical results provide valuable atomistic insight into the oxidation of TMD monolayers and are useful for the future design of TMD-based 2D devices.
Co-reporter:Pengbo Zhang, Tingting Zou, Jijun Zhao
Journal of Nuclear Materials 2015 Volume 467(Part 1) pp:465-471
Publication Date(Web):December 2015
DOI:10.1016/j.jnucmat.2015.09.039
We investigated the atomistic mechanism of He–He and He–metal interactions in bcc transition metals (V, Nb, Ta, Cr, Mo, W, and Fe) using first-principles methods. We calculated formation energy and binding energy of He–He pair as function of distance within the host lattices. The strengths of He–He attraction in Cr, Mo, W, and Fe (0.37–1.11 eV) are significantly stronger than those in V, Nb, and Ta (0.06–0.17 eV). Such strong attractions mean that He atoms would spontaneously aggregate inside perfect Cr, Mo, W, and Fe host lattices in absence of defects like vacancies. The most stable configuration of He–He pair is <100> dumbbell in groups VB metals, whereas it adopts close <110> configuration in Cr, Mo, and Fe, and close <111> configuration in W. Overall speaking, the He–He equilibrium distances of 1.51–1.55 Å in the group VIB metals are shorter than 1.65–1.70 Å in the group VB metals. Moreover, the presence of interstitial He significantly facilitates vacancy formation and this effect is more pronounced in the group VIB metals. The present calculations help understand the He-metal/He–He interaction mechanism and make a prediction that He is easier to form He cluster and bubbles in the groups VIB metals and Fe.
Co-reporter:Zhaoyang Zheng, Xue Jiang, Jijun Zhao
Chemical Physics Letters 2015 Volume 628() pp:76-80
Publication Date(Web):16 May 2015
DOI:10.1016/j.cplett.2015.03.058
•Lattice constants of MOF-EMs can be reproduced well by optPBE-vdW scheme.•NHN are metallic and other MOF-EMs are semiconducting.•Impact sensitivities of MOF-EMs are discussed in terms of their band gaps.•Calculated bulk moduli of MOF-EMs range from 15.1 GPa (CuAN) to 35.0 GPa (NHN).The structural, electronic and elastic properties for metal-organic frameworks (MOFs) as energetic materials are investigated using non-local density functional theory with dispersion correction. The lattice constants of MOF-EMs are reproduced well by optPBE-vdW functional. The electronic structure analysis reveals that NHN is a metal, while the others are semiconductors or insulators with band gap from 0.1 eV to 4.7 eV. NHP, CHP, CHHP and CuAN are predicted to be magnetic. We also discuss the impact sensitivities of MOF-EMs in terms of their electronic structures. The calculated bulk modulus ranges from 15.1 GPa (CuAN) to 35.0 GPa (NHN).
Co-reporter:Guobao Li, Xin Yue, Gaixia Luo, Jijun Zhao
Computational Materials Science 2015 Volume 106() pp:15-22
Publication Date(Web):August 2015
DOI:10.1016/j.commatsci.2015.04.027
•Electrode potentials of NaMO2 compounds with Na concentration were computed.•Activation energy for Na ions diffusion in NaMO2 compounds were computed.•The relation between variation trend of AEV and TM was further confirmed.•O3 type structures with V and Co are of great advantage for Na ions diffusion.Recently, researches on sodium (Na) batteries get resurrected for large-scale applications under the background of limited lithium (Li) resource. Four layered sodium transition-metal (TM) oxides, NaMO2 (M = V, Cr, Co and Ni) with α-NaFeO2 type as potential cathode materials for sodium batteries, are systematically investigated by first-principles calculations. As two key properties for a successful cathode material, the electrode potential and activation energy (EA) of Na+ in their parent NaMO2 crystals (M = V, Cr, Co and Ni) are reported. Our results suggest that NaCrO2, NaCoO2 and NaNiO2 possess satisfying average electrode potential (AEP) and their potential platforms are close to that of LiCoO2 to a large extent. Meanwhile, EA of these materials are all acceptable for batteries. These first-principles results suggest that layered NaCrO2, NaCoO2 and NaNiO2 are viable candidates for large-scale applications.
Co-reporter:Ruihuan Li, Wenbo Li, Chong Zhang, Pengbo Zhang, Hongyu Fan, Dongping Liu, Levente Vitos, Jijun Zhao
Journal of Nuclear Materials 2015 Volume 457() pp:36-41
Publication Date(Web):February 2015
DOI:10.1016/j.jnucmat.2014.10.062
In fusion environment, large amounts of helium (He) atoms are produced by transmutation along with structural damage in the structural materials, causing material swelling and degrading of physical properties. To understand the microscopic mechanism of He trapping in vacancies and voids, we explored He–vacancy interactions in HenVam (Va for vacancy) clusters (n, m = 1–4) and multiple He trapping in a 7-atom void of silicon carbide (SiC) by first-principles calculations. The binding energy between He and the HenVam clusters increases with the number of vacancies, while the vacancy binding energy gradually increases with the number of He atoms. Furthermore, a small cavity of about 0.55 nm in diameter can accommodate up to 14 He atoms energetically and the corresponding internal pressure is estimated to be 2.5 GPa. The tendency of He trapping in small voids provides an explanation for the experimentally observed He bubble formation at vacancy defects in SiC materials.
Co-reporter:Wenbo Li, Jijun Zhao, Dejun Wang
Solid State Communications 2015 Volume 205() pp:28-32
Publication Date(Web):March 2015
DOI:10.1016/j.ssc.2014.12.020
Highlights•A new amorphous SiO2/SiC interface model is generated using the first-principles method.•Calculated valence- and conduction-band offsets of this interface are 2.84 eV and 2.76 eV, respectively.•Accurate charge transition levels of several defects at the SiO2/SiC interface are estimated.•Silicon interstitial in SiO2 and carbon dimer can account for the large interface states experimentally observed.A defect-free structural model of the amorphous SiO2/4H-SiC(0001) interface is presented through first-principle calculations. Following the potential lineup method, we first calculate the valence- and conduction-band offsets of this interface, which are in good agreement with the experimental values. Based on this interface model, we create several typical interface defects and estimate the accurate charge transition levels of these defects within the HSE06 hybrid functional scheme. The results indicate that the silicon interstitial in SiO2 and carbon dimers in both SiC and SiO2 are the possible candidates for the large interface states experimentally observed near the conduction band of 4H-SiC.
Co-reporter:Lizhao Liu, Xin Yue, Jijun Zhao, Qian Cheng, Jie Tang
Physica E: Low-dimensional Systems and Nanostructures 2015 Volume 69() pp:316-321
Publication Date(Web):May 2015
DOI:10.1016/j.physe.2015.02.006
•Graphene antidot lattices show a formation energy per length of 0.50 ~ 0.60 eV/nm.•Within a hole density of 10%, graphene antidot lattices are almost metallic.•An optimum pore diameter of ~0.86 nm was found for penetration of K and Cl atoms.Thermodynamic stabilities and electronic properties of graphene antidot lattices with hexagonal holes were examined using density functional theory calculations and several crucial factors related to the applications of supercapacitors were discussed. For the graphene antidot lattices with different hole sizes, the formation energy per edge length is about 0.50∼0.60 eV/nm, which is comparable to that of graphene nanoribbon edges. Within a hole density of 10%, the graphene antidot lattices can maintain the excellent electronic properties of perfect graphene due to negligible intervalley scattering. Further increasing the hole density will open a band gap. Taking the potassium chloride (KCl) electrolyte as an example, we further investigated the diffusion behaviors of potassium (K) and chlorine (Cl) atoms through the graphene antidot lattices. It was shown that K and Cl atoms can go through the holes with nearly no barrier at an appropriate hole size of 0.54 nm, which gives an optimum pore diameter of ∼0.86 nm. Therefore, the excellent graphene-like electronic properties and good penetrability for ions suggest promising applications of graphene antidot lattices in the field of supercapacitors.Compared with perfect graphene, porous graphene antidot lattices facilitate the diffusion of ions in passing through the graphene layers, which should be promising electrode materials for supercapacitors.
Co-reporter:Fen Li, Xue Jiang, Jijun Zhao, Shengbai Zhang
Nano Energy 2015 Volume 16() pp:488-515
Publication Date(Web):September 2015
DOI:10.1016/j.nanoen.2015.07.014
•The progresses of GO application in hydrogen storage and photocatalytic water splitting have been elaborately summarized.•The progresses of GO application in lithium batteries and supercapacitors have been systematically discussed.•We summarized the versatile application of GO in air and water purification.Graphene oxide (GO), the functionalized graphene with oxygen-containing chemical groups, has recently attracted resurgent interests because of its superior properties such as large surface area, mechanical stability, tunable electrical and optical properties. Moreover, the surface functional groups of hydroxyl, epoxy and carboxyl make GO an excellent candidate in coordinating with other materials or molecules. Owing to the expanded structural diversity and improved overall properties, GO and its composites hold great promise for versatile applications of energy storage/conversion and environment protection, including hydrogen storage materials, photocatalyst for water splitting, removal of air pollutants and water purification, as well as electrode materials for various lithium batteries and supercapacitors. In this review, we present an overview on the current successes, as well as the challenges, of the GO-based materials for energy and environmental applications.A schematic showing the GO/RGO-based hybrid materials for energy and environmental applications along with the SCI-indexed journal publications until now (January, 2015), searched from ISI by the keyword of graphene oxide combined with another one listed in the outer circle (left panel) or the X axis (right panel).Figure optionsDownload full-size imageDownload as PowerPoint slide
Co-reporter:Xiaoming Huang;Yan Su;Linwei Sai;Vijay Kumar
Journal of Cluster Science 2015 Volume 26( Issue 2) pp:389-409
Publication Date(Web):2015 March
DOI:10.1007/s10876-014-0829-7
The low-energy structures of PtnSnn (n = 1–10) and Pt3mSnm (m = 1–5) clusters have been determined using genetic algorithm incorporated with density functional theory. Platinum and tin atoms tend to mix with each other due to the energetically favorable Pt–Sn bonds. However, due to the larger atomic radius of Sn atoms, we find segregation of Sn atoms on the surface of PtnSnn clusters. This leaves one or two Pt atoms available for reaction and for larger clusters segregation of Sn could block the Pt sites. For Pt3mSnm clusters, Sn atoms are well separated in the cluster structures and prefer to form sharp vertices leaving triangular faces of three Pt atoms available for reactivity. The electronic properties such as highest occupied molecular orbital–lowest unoccupied molecular orbital gap, distribution of frontier orbitals, Mayer bond order, Mülliken atomic charge, and the density of states are discussed. Significant hybridization between the d orbitals of Pt and the p orbitals of Sn is revealed. These theoretical results provide the general trends for the structural and bonding characteristics of the Pt–Sn alloy clusters and help understand their catalytic behavior.
Co-reporter:Jijun Zhao, Xiaoming Huang, Peng Jin, Zhongfang Chen
Coordination Chemistry Reviews 2015 s 289–290() pp: 315-340
Publication Date(Web):
DOI:10.1016/j.ccr.2014.12.013
Co-reporter:Xiaoming Huang
The Journal of Physical Chemistry C 2015 Volume 119(Issue 20) pp:10987-10994
Publication Date(Web):December 29, 2014
DOI:10.1021/jp5112845
Vanadium-doped silicon cluster anions, V3Sin– (n = 3–14), have been generated by laser vaporization and investigated by anion photoelectron spectroscopy. The vertical detachment energies (VDEs) and adiabatic detachment energies (ADEs) of these clusters were obtained. Meanwhile, genetic algorithm (GA) combined with density functional theory (DFT) calculations are employed to determine their ground-state structures systematically. Excellent agreement is found between theory and experiment. Among the V3Sin– clusters, V3Si5–, V3Si9–, and V3Si12– are relatively more stable. Generally speaking, three V atoms prefer to stay close with others and form strong V–V bonds. Starting from V3Si11–, cage configurations with one interior V atom emerge.
Co-reporter:Lingli Tang, Ruili Shi, Yan Su, and Jijun Zhao
The Journal of Physical Chemistry A 2015 Volume 119(Issue 44) pp:10971-10979
Publication Date(Web):October 15, 2015
DOI:10.1021/acs.jpca.5b08073
In order to understand the cage fusion behavior during the nucleation processes of methane hydrate (MH), methane-encapsulated double-cage clusters (CH4)2(H2O)n (n = 30–43) and several multicage structures with three or more cages were studied employing DFT-D methods. We find that almost all the lowest-energy double-cage structures can be constructed by merging the most stable structures of the monocage clusters CH4(H2O)n (n = 18–24). Double-cage structures can achieve higher stability through sharing a hexagon than a pentagon, which may be applicable to larger fused cage clusters. The preference of hexagons during cage fusion should be favorable for the appearance of the cages including hexagons such as the 51262, 51264 cages during the MH nucleation process. The symmetric C–H stretching modes of methane molecules in the double-cage structures show a clear trend of red shift with increasing size of the composing monocages. Compared with the case of monocages, the stretching frequencies of methane molecules in double-cage structures shift slightly, indicating variation of monocage configuration when cage fusion occurs. The larger multicage structures are found to possess higher fusion energies through sharing more polygons. Their thermodynamic stabilities do not simply increase with the number of fused monocages and are affected by the spatial arrangement of the building cages.
Co-reporter:Yan Su, Xingfa Gao, Jijun Zhao
Carbon 2014 Volume 67() pp:146-155
Publication Date(Web):February 2014
DOI:10.1016/j.carbon.2013.09.073
We used density functional theory to study the reaction mechanisms of chemical reduction of graphene oxide (GO) by the sulfur-containing compounds HSO3− and H2SO3. We studied the reaction energy profiles for the following reactions: dehydroxylation of GO with one and two hydroxyl groups, de-epoxidation of GO with one or two epoxy groups and decarboxylation and decarbonylation of GO with carboxyl and carbonyl groups. We found that hydroxyl and epoxide groups could be easily reduced because of the lower energy barriers, whereas decarboxylation and decarbonylation reactions are not kinetically and thermodynamically easy because of the higher energy barriers. These reaction mechanisms at the atomistic level are not only supported by Chen’s experimental results [J. Phys. Chem. C 2010, 114, 19885], but are also beneficial for the development of new agents that could efficiently reduce GO.
Co-reporter:Jie Fu, Chong Zhang, Jijun Zhao
Computational Materials Science 2014 Volume 85() pp:142-146
Publication Date(Web):1 April 2014
DOI:10.1016/j.commatsci.2013.12.061
•Developed Gupta potentials for alkaline earth elements Ca and Sr of fcc phase.•Reproduced lattice constants, cohesive energy, and elastic constants.•Examined phonon dispersion, equation of state and fcc–bcc transformation.•Predicted surface energy, vacancy formation energy.Empirical Gupta-type potentials have been developed for alkaline earth metals, Ca and Sr. The potential parameters were fitted according to the experimental lattice constants, cohesive energies and elastic constants. Compared to experiments and DFT calculations, our potentials can reproduce the bulk modulus, shear modulus, and sound speed very well. More importantly, our potentials can also describe phonon dispersion relationships of Ca and Sr solids correctly. Energy profiles for fcc–bcc transformation along the Bain deformation path has been calculated, which yields the correct energy difference between fcc and bcc phases. Our potentials are also applied to the low-coordination environments, i.e., vacancy formation energies and surface energies of low-index surfaces, and a qualitative agreement with experiment and previous theoretical results is obtained.
Co-reporter:Jingcheng Xu, Jijun Zhao, Pavel Korzhavyi, Börje Johansson
Computational Materials Science 2014 Volume 84() pp:301-305
Publication Date(Web):March 2014
DOI:10.1016/j.commatsci.2013.12.032
•Disordered BCC Fe–Cr–W solid solutions are studied using ab initio calculations.•Cubic elastic constants are calculated as a function of composition in the range of 7.8–10 wt.% Cr, and 1–2 wt.% W.•The effect of alloying elements on the elastic properties of polycrystalline alloys is analyzed.Reduced activation ferritic/martensitic (RAFM) steels are considered as the primary candidate for the blanket module in International Thermonuclear Experimental Reactor (ITER), and the basic composition of RAFM steels is Fe–Cr–W alloy. Using the exact muffin-tin orbitals method (EMTO) combined with coherent potential approximation (CPA), we investigated composition dependence of elastic properties for Fe–Cr–W random alloys in the composition range of 7.8–10.0 wt.% of Cr and 1.0–2.0 wt.% of W. Bulk modulus, shear modulus, Young’s modulus, B/G ratio, and Poisson ratio are discussed as functions of ternary composition. The elastic moduli of Fe–Cr–W alloys increase with chromium content in the studied range of alloy compositions, the effect of tungsten is slight.
Co-reporter:Pengbo Zhang, Chong Zhang, Ruihuan Li, Jijun Zhao
Journal of Nuclear Materials 2014 Volume 444(1–3) pp:147-152
Publication Date(Web):January 2014
DOI:10.1016/j.jnucmat.2013.09.048
Co-reporter:Chong Zhang, Jie Fu, Ruihuan Li, Pengbo Zhang, Jijun Zhao, Chuang Dong
Journal of Nuclear Materials 2014 Volume 455(1–3) pp:354-359
Publication Date(Web):December 2014
DOI:10.1016/j.jnucmat.2014.07.011
Chinese low activation martensitic steel (CLAM) has been designed with decreased W content and increased Ta content to improve performance. We performed first-principles calculations to investigate the diffusion properties of solute element (Cr, W, Mn, V, Ta) and C diffusion with a nearby solute element inside bcc Fe. The self-diffusion coefficients and solute diffusion coefficients in Fe host were derived using the nine-frequency model. A relatively lower diffusivity was observed for W in paramagnetic state, implying enriched W concentration inside Fe host. The solute atom interacts strongly with C impurity, depending on the interatomic distance. According to our calculations, formation of Ta carbide precipitates is energetically preferred by trapping C impurity around Ta atom. Our theoretical results are helpful for investigating the evolution of microstructure of steels for engineering applications.
Co-reporter:Pengbo Zhang, Jijun Zhao
Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 2014 Volume 322() pp:34-40
Publication Date(Web):1 March 2014
DOI:10.1016/j.nimb.2013.12.030
We report an atomistic study of energetics and configurations of He–He pair in vacancy of bcc transition metals (V, Nb, Ta, Cr, Mo, W, and Fe) using first-principles methods. The results show He-vacancy attractions in the group VB metals are 2–3 times weaker than those in Cr, Mo, W and Fe. The 〈111〉 dumbbell configuration for He–He pair in vacancy is the most stable except for V. Calculated formation energies of He–He pair in vacancy of group VB metals (3.2–3.95 eV) are systematically lower than those of group VIB and VIIIB metals (4.25–5.67 eV), while He–He distance for V metal is greater than other metals. He–metal repulsion is stronger than He–He due to longer He–metal distance in metal vacancy. Calculated densities of states provide a reasonable explanation for most stable configuration; thus the stability of He–He pairs depends strongly on the electronic structure of the host and insignificantly on its atomic size. Moreover, the group-specific trends of He–He and self-interstitial configurations are discussed.
Co-reporter:Lizhao Liu;Feng Liu
Nano Research 2014 Volume 7( Issue 5) pp:626-657
Publication Date(Web):2014 May
DOI:10.1007/s12274-014-0431-1
Co-reporter:Fen Li, Junfeng Gao, Jian Zhang, Fen Xu, Jijun Zhao and Lixian Sun
Journal of Materials Chemistry A 2013 vol. 1(Issue 27) pp:8016-8022
Publication Date(Web):01 May 2013
DOI:10.1039/C3TA10800G
The incorporation of lithium amidoborane (LiAB) into graphene oxide (GO) and the dehydrogenation process of the GO–LiAB complex have been investigated for combining the chemical and physical hydrogen storage approaches. The obtained adsorption energy and minimum energy pathway (MEP) demonstrate that both of the two dominant groups of –O– and OH– contribute to the facile combination of GO and LiAB (GO3–LiAB). The GO3–LiAB complex has a better dehydrogenation performance than the pristine LiAB, which also indicates the feasibility of doping Li atoms on the GO surface. Atomic charge and bond length analyses match well with the MEP prediction. By dehydrogenation of the GO3–LiAB complex, we also achieve uniform metal doping on the GO surface for physisorption of H2. The GO–Li(n) products can store up to 5 wt% H2 and the GO–(Li3N3B3)(n) can still store 5 wt% H2. The dehydrogenation product of the GO3–LiAB complex has bridged the chemical and physical hydrogen storage approaches to move towards on-board hydrogen storage applications, which expands the scope for designing more efficient hydrogen storage materials.
Co-reporter:Junfeng Gao, Junfeng Zhang, Hongsheng Liu, Qinfang Zhang and Jijun Zhao
Nanoscale 2013 vol. 5(Issue 20) pp:9785-9792
Publication Date(Web):25 Jul 2013
DOI:10.1039/C3NR02826G
In the fabrication and processing of silicene monolayers, structural defects are almost inevitable. Using ab initio calculations, we systemically investigated the structures, formation energies, migration behaviors and electronic/magnetic properties of typical point defects in silicene, including the Stone-Wales (SW) defect, single and double vacancies (SVs and DVs), and adatoms. We found that SW can be effectively recovered by thermal annealing. SVs have much higher mobility than DVs and two SVs are very likely to coalesce into one DV to lower the energy. Existence of SW and DVs may induce small gaps in silicene, while the SV defect may transform semimetallic silicene into metallic. Adatoms are unexpectedly stable and can affect the electronic properties of silicene dramatically. Especially, Si adatoms as self-dopants in silicene sheets can induce long-range spin polarization as well as a remarkable band gap, thus achieving an all-silicon magnetic semiconductor. The present theoretical results provide valuable insights into identification of these defects in experiments and understanding their effects on the physical properties of silicene.
Co-reporter:Gaixia Luo, Lizhao Liu, Junfeng Zhang, Guobao Li, Baolin Wang, and Jijun Zhao
ACS Applied Materials & Interfaces 2013 Volume 5(Issue 21) pp:11184
Publication Date(Web):October 17, 2013
DOI:10.1021/am403427h
One great challenge for supercapacitor is to achieve high energy capacity and fast charge/discharge rates simultaneously. Porous graphene with large surface area is a promising candidate for electrode materials of supercapacitor. Using first-principles calculations and non-equilibrium Green’s function technique, we have explored the formation energies, mechanical properties, diffusion behaviors and electrical conductance of graphene sheets with various hole defects and/or nitrogen doping. Interestingly, graphene sheets with pyridinic-like holes (especially hexagonal holes) can be more easily doped with nitrogen and still retain the excellent mechanical properties of pristine graphene that is beneficial for the long cycle life. Porous graphene electrode with moderate hole diameter of 4.2–10 Å facilitates efficient access of electrolyte and exhibit excellent rate capability. In addition, doping with nitrogen as electron donors or proton attractors leads to charge accumulation and generates higher pseudocapacitance. Transmission coefficients of N-doped graphene sheets with pyridinic-like holes are only moderately reduced with regard to that of pristine graphene and are insensitive to the detailed geometry parameters. Overall, N-doped graphene with pyridinic-like holes exhibits exciting potentials for high performance energy storage in supercapacitor devices.Keywords: doping; graphene; hole; supercapacitors;
Co-reporter:Junfeng Zhang, Jijun Zhao
Carbon 2013 Volume 55() pp:151-159
Publication Date(Web):April 2013
DOI:10.1016/j.carbon.2012.12.021
The graphene grain boundaries with periodic length up to 18 Å have been studied using density functional theory. Atomic structures, thermodynamic stabilities and electronic properties of 40 grain boundaries with symmetric and nonsymmetric structures were investigated. According to the arrangements of pentagons and heptagons on the boundary, grain boundaries were cataloged into four classes. Some nonsymmetric grain boundaries constructed here have identical misorientation angles to the experimentally observed ones. The formation energies of grain boundaries can be correlated with the misorientation angle and inflection angle. Nonsymmetric grain boundaries possess comparable formation energies to their symmetric counterparts when the periodic length along the defect line is larger than 1 nm. Analysis of electronic density of states shows that the existence of a grain boundary usually increases the density of states near the Fermi level, whereas some symmetric grain boundaries can open a small band gap due to local sp2-to-sp3 rehybridization.
Co-reporter:Lizhao Liu, Junfeng Gao, Xu Guo and Jijun Zhao
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 40) pp:17134-17141
Publication Date(Web):14 Aug 2013
DOI:10.1039/C3CP52904E
Atomic structural models of zigzag-shaped carbon nanotubes (Z-CNTs) were constructed by periodically introducing pentagons and heptagons into pristine CNTs. In terms of formation energies, the Z-CNTs present comparable energetic stabilities to those of the pristine CNTs and are more stable than C60 fullerene. The mechanical properties of these Z-CNTs, including the Young's modulus, intrinsic strength and failure behaviour, were systematically investigated by first-principles computations. Compared with the pristine CNTs with an average Young's modulus of about 1.0 TPa, incorporation of pentagons and heptagons in the Z-CNTs will reduce the average Young's modulus to several hundreds of GPa. Moreover, the computational results also showed that under uniaxial tensile strain, the intrinsic strength and failure strain of the Z-CNTs are also lower than those of the pristine CNTs. Generally, the Young's modulus and intrinsic strength of the Z-CNTs are exponentially inverse to curvature, which can be expressed by simple formulae. In particular, the electronic properties of the armchair Z-CNTs can be tailored by uniaxial tensile strain. It was also found that through applying tensile strain, a semiconductor–metal or metal–semiconductor transition can be triggered. The localized–delocalized partial charge distribution near the Fermi energy for the strained Z-CNTs can explain the semiconductor–metal or metal–semiconductor transition. This significant electromechanical coupling effect suggests the Z-CNTs have potential applications in nanoscale electromechanical sensors and switches.
Co-reporter:Han Zhang, Zheng Duan, Xiaonan Zhang, Chao Liu, Junfeng Zhang and Jijun Zhao
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 28) pp:11794-11799
Publication Date(Web):15 May 2013
DOI:10.1039/C3CP44716B
We present a molecular dynamic simulation on the mechanical strength and fracture behavior of graphene grain boundaries (GBs). The intrinsic strength, critical failure strain, and failure mechanism of graphene GBs mainly rely on the temperature and inflection angle, whereas the Young's modulus does not vary significantly with either temperature or boundary configuration. The overall intrinsic strengths of inflected GBs can be correlated with infection angle by a linear term, which is irrelevant to the system temperature. The initial failure sites of GBs locate either on the boundary line or inside the domain at high temperature.
Co-reporter:Fen Li, Lixian Sun, Jijun Zhao, Fen Xu, Huai-Ying Zhou, Qing-Ming Zhang, Feng-Lei Huang
International Journal of Hydrogen Energy 2013 Volume 38(Issue 17) pp:6930-6937
Publication Date(Web):10 June 2013
DOI:10.1016/j.ijhydene.2013.03.076
•H2 generation behavior of Al16Mg and Al16Bi clusters was measured.•A feasible mechanism for hydrolysis reaction of Al–Bi alloy was proposed.•Bi in Al16Bi cluster can accelerate proton-transport in hydrolysis reaction.We have systematically investigated the hydrolysis mechanism of metal doped Al16M (M = Al, Mg and Bi) clusters with H2O molecules and proposed a reasonable elucidation for the experimentally observed fast H2 generation rate and high H2 yield in the Al–Bi based composite. Mg and Bi showed negative effect on the dissociation process of the first H2O molecule, but accelerated further H2 generation process. The investigation of persistent hydrolysis reactions demonstrated that the proton-transfer way makes the aluminum–water reaction a lasting process in the long-term H2 generation in existence of Bi atom, which explains not only the previously observed fast H2 generation rate but also high H2 yields in the Bi added Al powder. Our experimental results of hydrogen generation form Al–Bi (Mg) mixture and water are in good agreement with the theory prediction. The facilitated hydrolysis reaction in Al16Bi cluster is attributed to the weakened hydroxide adsorption with the presence of Bi in the aluminum cluster, which is the key factor to accelerate the proton-transfer process.
Co-reporter:Gaixia Luo, Jijun Zhao, Baolin Wang
Computational Materials Science 2013 Volume 68() pp:212-217
Publication Date(Web):February 2013
DOI:10.1016/j.commatsci.2012.10.027
Using first-principles computations, we investigated the effect of the graphite interlayer spacing and substitutional boron doping on the storage capacity of Li ions. We found that increasing the distance between graphite layers only moderately increases the capacity because of a combination of geometric and electronic effects. Doping with boron results in a noticeable increase in the saturation Li intercalation density by about 33.3% at boron contents of BC15 and BC7 with regard to the pristine graphite predicts that the maximum of saturation Li intercalation density locates at around 10 at% of boron content. The electronic structures of Li-intercalated graphite systems were analyzed to explain these effects.
Co-reporter:Chong Zhang, Pengbo Zhang, Ruihuan Li, Jijun Zhao, Chuang Dong
Journal of Nuclear Materials 2013 Volume 442(1–3) pp:370-376
Publication Date(Web):November 2013
DOI:10.1016/j.jnucmat.2013.08.032
Addition of trace amounts of Al, Si and Y into V–4Cr–4Ti alloy is beneficial for the mechanical properties under irradiation. It is thus important to investigate the influence of solute/trace elements on stabilities, energetics and diffusion behaviors of vacancy defects. We performed first-principles calculations to evaluate vacancy–solute/trace interaction inside dilute V–X (X = Ti, Cr, Al, Si, Y) and V–4Cr–4Ti–(Al, Si, Y) alloys. With addition of Si and Y, vacancy-based complexes tend to form near Ti–Si and Ti–Y pairs, while the effect of Al is negligible. Moreover, diffusion coefficients of solute/trace element in vanadium were derived using nine-frequency model. With high binding energy and low diffusion coefficient, Si atom is strongly attractive to vacancy in vanadium matrix. Our theoretical results suggest that the interactions between vacancy and solute/trace elements play some role in the evolution of microstructures inside vanadium alloys.
Co-reporter:Pengbo Zhang, Ruihuan Li, Jijun Zhao, Bin Wen
Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 2013 Volume 303() pp:75-80
Publication Date(Web):15 May 2013
DOI:10.1016/j.nimb.2012.10.043
Using first-principles calculations, we investigate the stability and energetics of H–He–vacancy complexes and He–He/He–H/H–H interactions in bulk vanadium to explore the synergetic effect of H and He impurities with vacancy. Inside vacancy space, He prefers to occupy octahedral site rather than vacancy center, different from of the cases of bcc Fe and W. Equilibrium He–H (2.33 Å) and H–H (2.45 Å) distances in the vacancy are longer than He–He distance (1.75 Å) and they exhibit a weak attraction. One He–vacancy complex can trap up to five H atoms and the stable configurations with different amount of trapped H atoms are discussed in details. If a He atom occupies vacancy center, formation of H2 molecule in the He–vacancy complex is almost impossible. Moreover, formation energy of a new vacancy around the H–He–vacancy complex does not remarkably reduce with increasing number of H atoms. We thus suggest that creation of a new vacancy around the He-complex via H aggregation is thermodynamically difficult.
Co-reporter:Ruihuan Li, Pengbo Zhang, Xiaoqing Li, Chong Zhang, Jijun Zhao
Journal of Nuclear Materials 2013 435(1–3) pp: 71-76
Publication Date(Web):
DOI:10.1016/j.jnucmat.2012.12.022
Co-reporter:Ruihuan Li, Pengbo Zhang, Chong Zhang, Xiaoming Huang, Jijun Zhao
Journal of Nuclear Materials 2013 440(1–3) pp: 557-561
Publication Date(Web):
DOI:10.1016/j.jnucmat.2013.03.068
Co-reporter:Hongsheng Liu, Junfeng Gao, and Jijun Zhao
The Journal of Physical Chemistry C 2013 Volume 117(Issue 20) pp:10353-10359
Publication Date(Web):April 30, 2013
DOI:10.1021/jp311836m
Silicene, a two-dimensional hexagonal lattice of silicon, has been synthesized recently and exhibits fascinating electronic properties that resemble graphene. The substrate effect on the electronic properties of silicene is important for the practical applications of silicene. First-principles calculations were performed for silicene on two kinds of representative inert substrates, that is, hexagonal boron nitride (h-BN) monolayer and SiC(0001) surface. The silicene–substrate interaction energies range in 0.067–0.089 eV per Si atom, belonging to typical van der Waals interaction. The characteristic Dirac cone is preserved for silicene on h-BN monolayer or hydrogenated Si-terminated SiC(0001) surface. On the other hand, the silicene becomes metallic when it is placed on a hydrogenated C-terminated SiC(0001) surface. This effect was explained by the work functions for silicene and the substrates. The present results provide some guidelines for selecting proper substrates for silicene in future microelectronic devices.
Co-reporter:Chong Zhang;Hua Tian;Chuanpu Hao;Qing Wang
Journal of Materials Science 2013 Volume 48( Issue 8) pp:3138-3146
Publication Date(Web):2013 April
DOI:10.1007/s10853-012-7091-x
Using the first-principles calculations, a cluster-plus-glue-atom model was employed to investigate the elastic and electronic properties of Ti–Mo–Nb alloys with cluster formula of [MoTi14] (glue atom)x (glue atom = Ti, Mo, Nb, x = 1 or 3) for a theoretical guidance in composition design of β titanium alloys. The bulk modulus, shear modulus, Young’s modulus, and Poisson ratio were evaluated from the calculated elastic constants using Voigt–Reuss–Hill average scheme on the periodic supercell model of cluster packing. The electronic properties of the Ti–Mo–Nb alloys were discussed by analyzing the electron density of state and Mulliken population. Meanwhile, we designed two series of Ti–Mo–Nb alloys, i.e., [MoTi14]X1 (X = Ti, Mo, Nb) and [YTi14]Nb3 (Y = Ti, Mo), and experimentally measured their mechanical properties. Our theoretical results (including mass density, Young’s modulus, ductility) based on our cluster packing model agreed well with the experimental data, especially for [TiMo14]X1 (X = Ti, Mo, Nb) alloy series. On the contrary, the random solid solution structures were mechanically unstable and the calculated values significantly deviated from the experiments. Based on the cluster-plus-glue-atom model, an Ashby map of E/ρ versus B/G was constructed and indicated the inverse correlation between stiffness and ductility, for which the random solid solution model was unable to reflect. The Mo/Ti = 1/14 rule derived from the cluster model may serve as an important guideline for composition design of Ti–Mo based systems to achieve low elastic modulus alloys with stable β phase.
Co-reporter:Junfeng Gao ; Jijun Zhao ;Feng Ding
Journal of the American Chemical Society 2012 Volume 134(Issue 14) pp:6204-6209
Publication Date(Web):March 14, 2012
DOI:10.1021/ja2104119
In vacuum, the bare zigzag (zz) edge of graphene is reconstructed into a line of pentagon–heptagon pairs, while the pristine armchair (ac) edge is retained. Our first-principle explorations of graphene edges on three metal surfaces [Cu(111), Co(111), and Ni(111)] indicate an opposite tendency, that is, the pristine zz edge is energetically favorable and the reconstructed ac edge with dangling C atoms is highly stable on Co(111) and Ni(111) surfaces. Insightful analysis shows that passivation of the graphene edge by metal surfaces is responsible for the dramatic differences. Beyond this, the unique edge configuration has a significant impact on the graphene CVD growth behavior.
Co-reporter:Lizhao Liu, Junfeng Zhang, Jijun Zhao and Feng Liu
Nanoscale 2012 vol. 4(Issue 19) pp:5910-5916
Publication Date(Web):30 Jul 2012
DOI:10.1039/C2NR31164J
The mechanical properties, including the Young's modulus and intrinsic strength, of graphene oxides are investigated by first-principles computations. Structural models of both ordered and amorphous graphene oxides are considered and compared. For the ordered graphene oxides, the Young's modulus is found to vary from 380 to 470 GPa as the coverage of oxygen groups changes, respectively. The corresponding variations in the Young's modulus of the amorphous graphene oxides with comparable coverage are smaller at 290–430 GPa. Similarly, the ordered graphene oxides also possess higher intrinsic strength compared with the amorphous ones. As coverage increases, both the Young's modulus and intrinsic strength decrease monotonically due to the breaking of the sp2 carbon network and lowering of the energetic stability for the ordered and amorphous graphene oxides. In addition, the band gap of the graphene oxide becomes narrower under uniaxial tensile strain, providing an efficient way to tune the electronic properties of graphene oxide-based materials.
Co-reporter:Lizhao Liu, Lu Wang, Junfeng Gao, Jijun Zhao, Xingfa Gao, Zhongfang Chen
Carbon 2012 Volume 50(Issue 4) pp:1690-1698
Publication Date(Web):April 2012
DOI:10.1016/j.carbon.2011.12.014
Based on the experimental observations, amorphous structural models of graphene oxides (GOs) were constructed and investigated by first-principles computations. Geometric structures, thermodynamic stabilities, and electron density of states of these amorphous GO models were examined and compared with the previously proposed ordered GO structures. The thermodynamically most favorable amorphous GO models always contain some locally ordered structures in the short range, due to a compromise of the formation of hydrogen bonds, the existence of dangling bonds, and the retention of the π bonds. Compared to the ordered counterparts, these amorphous GO structures possess good stability at low oxygen coverage. Varying the oxygen coverage and the ratio of epoxy and hydroxyl groups provides an efficient way to tune the electronic properties of the GO-based materials.
Co-reporter:Hua Tian, Yingmin Wang, Lei Liu, Xiaoqing Li, Chong Zhang, Jijun Zhao, Chuang Dong, Bin Wen
Journal of Non-Crystalline Solids 2012 Volume 358(Issue 15) pp:1730-1734
Publication Date(Web):1 August 2012
DOI:10.1016/j.jnoncrysol.2012.04.014
Ni–Ta bulk metallic glass (BMG) with compositions around Ni60Ta40 is a newly found binary BMG with high glass forming ability and extraordinary mechanical strength. Using ab initio molecular dynamics, the local atomic structure, elastic properties and electronic structures of Ni60Ta40 glass have been explored. The pair-correlation functions, coordination numbers, and chemical compositions of the most abundant local clusters have been analyzed. We demonstrated the existence of icosahedral Ni7Ta6 clusters as the major Ni-centered clusters, while the most popular Ta-centered cluster is Ta7Ni8. These findings agree with our previous cluster model of Ni–Ta binary BMG. The elastic moduli of Ni60Ta40 glass were also computed and the experimental Young's modulus is well reproduced. Analysis of electronic structures further revealed that the interaction between d electrons of Ni and Ta atoms is responsible for the experimentally observed ultrahigh mechanical strength for the Ni–Ta BMGs.Highlights► Using ab initio molecular dynamics, we study the Ni60Ta40 glass. ► The findings of the atomic structure agree with our previous cluster model. ► The experimental Young's modulus is well reproduced. ► We revealed the interaction between d electrons of Ni and Ta atoms.
Co-reporter:Linwei Sai, Lingli Tang, Xiaoming Huang, Guibin Chen, Jijun Zhao, Jun Wang
Chemical Physics Letters 2012 Volume 544() pp:7-12
Publication Date(Web):20 August 2012
DOI:10.1016/j.cplett.2012.06.050
Abstract
Using genetic algorithm combined with density functional theory calculations, we performed global search for the most stable structures of (WO3)n clusters for n = 2–12. Small (WO3)n clusters with n = 3 or 4 adopt ring-like configurations with W–O alternating arrangement. Starting from (WO3)8, the tungsten oxide clusters transform to symmetric spherical-like cages. The relative stability, HOMO–LUMO gap, electronic states of these (WO3)n clusters were discussed. Analysis of wavefunctions of frontier orbitals and electron density of states shows that the valence bands are dominated by the 2p electrons from oxygen and the conduction bands are mainly contributed by the 5d states from tungsten.
Co-reporter:Yuan Wen, Chunqiang Zhuang, Xue Jiang, Jijun Zhao, Xin Jiang
Diamond and Related Materials 2012 Volumes 27–28() pp:14-18
Publication Date(Web):July–August 2012
DOI:10.1016/j.diamond.2012.05.003
Based on first-principles calculations within a random solid solution model, we investigated the formation energies and mechanical properties (including Young's modulus, Vickers hardness, and ductility) as function of compositions on the ternary phase diagram of B–C–N crystals. According to our calculations, compositions in the C-rich (70–90%) areas possess relatively high Vickers hardness (up to 77.5 GPa), B-rich (30–35%) areas have better ductility, and the sp3 crystalline phase does not exist stably in the N-rich areas (15–35%). As for the formation energies, the most favorable area locates around the C-poor area (30–50%) in the diagram. These theoretical results provide useful insight for designing and synthesizing the ternary B–C–N crystal of desirable mechanical properties.Highlights► Searching potential superhard phases in large areas of the B–C–N diagram. ► Building the relationship of compositions and mechanical properties. ► Finding superhard areas in the B–C–N phase diagram.
Co-reporter:Xue Jiang, Chunqiang Zhuang, Jijun Zhao, Xin Jiang
Diamond and Related Materials 2012 Volume 23() pp:44-49
Publication Date(Web):March 2012
DOI:10.1016/j.diamond.2011.12.014
Based on ab initio molecular dynamics method, four amorphous CNx structures with different stoichiometries (CN0.47, CN0.67, CN0.92, and CN1.3) were generated within a 100-atom supercell. Characterizations of the pair correlation functions, bond length and the fraction of bond types of the amorphous carbon nitrides revealed that the N content in such structures plays a key role in determining the obtained structural networks. The calculations of the cohesive energy give direct evidence that experimentally synthesized amorphous carbon nitride films usually possess low nitrogen concentration. The competition relationship between sp2 and sp3 bonds was analyzed by electronic effect and electron density of states. The bonding states, especially the sp3, affect greatly the hardness of these a-CNx compounds, implying a possible way to enhance the hardness of films via tuning the bonding states.Highlights► We generated four with different stoichiometries a-CNx using MD method. ► The competition between sp2–sp3 of a-CNx is interpreted by the electronic effect. ► No superhard stoichiometry a-CNx was found. ► The hardness of the a-CNx was determined by the fraction of sp3 C–C bond.
Co-reporter:Liang Hong, Haoliang Wang, Jingxin Cheng, Xiaoming Huang, Linwei Sai, Jijun Zhao
Computational and Theoretical Chemistry 2012 Volume 993() pp:36-44
Publication Date(Web):1 August 2012
DOI:10.1016/j.comptc.2012.05.027
Genetic algorithm combined with first-principle calculations is used to globally search the potential energy surface of the most stable configurations of elementary Aum and Agn clusters, as well as AumAgn (5 ⩽ m + n ⩽ 12) binary clusters. The effects of size and composition (i.e., Au:Ag ratio) on the atomic structures, coordination numbers and electronic properties including the binding energies and formation energies of Au–Ag binary clusters are discussed. The critical Au:Ag ratios for the 2D–3D transition are obtained and it is found that Ag atoms sometimes play a more important role in determining the ground-state configuration of a Au–Ag bimetallic cluster. The electron density of states is further analyzed to explore the influence of Au and Ag atoms. Stronger s–d hybridization originated from relativistic effects of Au atom is observed in the planar structure with regard to the 3D structures.Graphical abstractHighlights► DFT-based genetic algorithm search most stable structures of Au–Ag clusters. ► Competition between planar and non-planar structures due to Au and Ag. ► Ag atoms prefer peripheral positions while Au atoms occupy central ones.
Co-reporter:Yuan Liu, Jijun Zhao, Jingcheng Xu
Computational and Theoretical Chemistry 2012 Volume 991() pp:165-173
Publication Date(Web):1 July 2012
DOI:10.1016/j.comptc.2012.04.016
Using molecular dynamic with consistent valence force field (CVFF), we simulated the dissociation behavior of carbon dioxide hydrate, along with methane and hydrogen hydrates for comparison. The detailed dissociation process is discussed in terms of structural snapshots, radial distribution functions, mean square displacements, and diffusion coefficient at different temperatures. With increasing temperature, the clathrate skeleton of water molecules is firstly distorted and damaged; then the encapsulated CO2 molecules escape and are distributed in the aqueous solution in the form of small CO2 bubbles; finally the clathrate skeleton of hydrate is completely destroyed and CO2 molecules aggregate into a large CO2 bubble. The dissociation of carbon dioxide hydrate is not only related to the overall occupancy, but also sensitive to the residence of specific cavities. Moreover, the diffusion barriers for guest molecules (CO2, CH4, H2) penetrating the water cages are obtained from ab initio calculations, which help illustrate the physical origin for the difference in the dissociation behavior of the corresponding hydrates, that is, CO2 and CH4 molecules have to diffuse after the water cages break while H2 can diffuse in the earlier stage due to low diffusion barrier.Graphical abstractHighlights► The dissociation process of carbon dioxide hydrate is revealed. ► The dissociation of sI CO2 hydrate is sensitive to the residence of specific cavities. ► The dissociation mechanism of gas hydrates are related to the diffusion barriers of guest molecule.
Co-reporter:Junfeng Zhang, Jijun Zhao, and Jianping Lu
ACS Nano 2012 Volume 6(Issue 3) pp:2704
Publication Date(Web):February 27, 2012
DOI:10.1021/nn3001356
As one-dimension line defects, grain boundaries (GBs) can affect many intrinsic properties of graphene. In this paper, the mechanical properties of 20 representative graphene grain boundaries were studied using density functional theory and molecular dynamics. With different arrangements of the pentagonal and heptagonal rings, the grain boundary may remain flat or become inflected up to 72°. For the flat GBs, the intrinsic tensile strength decreases linearly with the formation energy with a maximum value of 93 GPa, close to that of a perfect graphene. The intrinsic tensile strength of the inflected GBs is found to generally decrease with increasing inflection angle. Stone–Wales transformation is identified as the major failure mechanism of graphene GBs at high temperatures, whereas the initial fracture site can be either on the boundary line or inside the domain. These theoretical results constitute a useful picture of the grain boundary effect on the mechanical properties of polycrystalline graphene.Keywords: defects; grain boundary; graphene; intrinsic strength; tensile strain
Co-reporter:Pengbo Zhang, Jijun Zhao, Bin Wen
Journal of Nuclear Materials 2012 423(1–3) pp: 164-169
Publication Date(Web):
DOI:10.1016/j.jnucmat.2012.01.027
Co-reporter:Pengbo Zhang, Jijun Zhao, Bin Wen
Journal of Nuclear Materials 2012 429(1–3) pp: 216-220
Publication Date(Web):
DOI:10.1016/j.jnucmat.2012.06.018
Co-reporter:Lizhao Liu, Yan Su, Junfeng Gao, Jijun Zhao
Physica E: Low-dimensional Systems and Nanostructures 2012 Volume 46() pp:6-11
Publication Date(Web):September 2012
DOI:10.1016/j.physe.2012.08.011
First-principles calculations were carried out to explore the geometry, energetic stability and magnetic property of cobalt nanoclusters adsorbed on graphene. Three structural isomers of Co13 clusters (icosahedron, fcc, and hcp) were chosen as model systems for ferromagnetic nanoclusters. The adsorption energies for these three types of Co13 clusters supported on graphene range between 1.66 and 2.21 eV. Both binding energy and adsorption energy first decrease with reducing coverage and then converge to invariant values at the coverage of 21.7%, which is originated from the interaction between the Co13 cluster and its periodic images. Regardless the detailed cluster configuration, the cluster–graphene interaction reduces the magnetic moment of the Co13 cluster by ∼20% compared with its free-standing counterpart. Most importantly, induced magnetic moments of graphene due to Co13 adsorption were observed, which may offer a new opportunity to tune the magnetic properties of graphene for its device applications.Highlights► Chemical adsorption between the Co13 cluster and the graphene was predicted. ► Both the binding and adsorption energies converge to a constant for R<21.7%. ► Formation of CoC carbide reduces magnetic moment of the Co13 cluster by ∼20%. ► Induced magnetic moments of graphene for adsorption of Co13 clusters were found.
Co-reporter:Qinghong Yuan ; Junfeng Gao ; Haibo Shu ; Jijun Zhao ; Xiaoshuang Chen ;Feng Ding
Journal of the American Chemical Society 2011 Volume 134(Issue 6) pp:2970-2975
Publication Date(Web):November 12, 2011
DOI:10.1021/ja2050875
Ground-state structures of supported C clusters, CN (N = 16, ..., 26), on four selected transition metal surfaces [Rh(111), Ru(0001), Ni(111), and Cu(111)] are systematically explored by ab initio calculations. It is found that the core–shell structured C21, which is a fraction of C60 possessing three isolated pentagons and C3v symmetry, is a very stable magic cluster on all these metal surfaces. Comparison with experimental scanning tunneling microscopy images, dI/dV curves, and cluster heights proves that C21 is the experimentally observed dominating C precursor in graphene chemical vapor deposition (CVD) growth. The exceptional stability of the C21 cluster is attributed to its high symmetry, core–shell geometry, and strong binding between edge C atoms and the metal surfaces. Besides, the high barrier of two C21 clusters’ dimerization explains its temperature-dependent behavior in graphene CVD growth.
Co-reporter:Lizhao Liu, Lei Zhang, Haili Gao, Jijun Zhao
Carbon 2011 Volume 49(Issue 13) pp:4518-4523
Publication Date(Web):November 2011
DOI:10.1016/j.carbon.2011.06.062
Based on armchair carbon nanotubes (CNTs), we construct the structural models of symmetric armchair carbon nanotori of different tubular/radial diameters. Tight-binding (TB) calculations show that the energetic stabilities of these carbon nanotori rely on their symmetries and tubular diameters closely. Density functional theory (DFT) calculations on a carbon nanotorus with substitutional B/N doping reveal that B(N) dopant prefers to occupy heptagonal (pentagonal) sites in order to form a steady six π-electrons orbital. After B/N doping, the electron density of states (DOSs) near the Fermi energy are notablely enhanced. The pristine and doped carbon nanotori with diverse geometric and electronic properties provide new opportunities in the applications of nanotechnology.
Co-reporter:Lu Wang, Jijun Zhao, Lili Wang, Tianyin Yan, Yi-Yang Sun and Shengbai B. Zhang
Physical Chemistry Chemical Physics 2011 vol. 13(Issue 47) pp:21126-21131
Publication Date(Web):24 Oct 2011
DOI:10.1039/C1CP21778J
We propose titanium-decorated graphene oxide (Ti–GO) as an ideal sorbent for carbon monoxide (CO) capture and separation from gas mixtures. Based on first-principles calculations, Ti–GO exhibits a large binding energy of ∼70 kJ mol−1 for CO molecules, while the binding energies for other gases, such as N2, CO2, and CH4, are significantly smaller. The gas adsorption properties of Ti–GO are independent of the local GO structures once Ti atoms are anchored by the oxygen-containing groups on the GO surface. The strong interaction between CO molecule and Ti is a result of dative bonding, i.e., hybridization between an empty d orbital of Ti and an occupied p orbital of CO. Adsorption isotherms from grand canonical Monte Carlo simulations clearly demonstrate the strong selectivity of Ti–GO for CO adsorption in a mixture with other gas.
Co-reporter:Xue Jiang, Jijun Zhao, Xin Jiang
Computational Materials Science 2011 Volume 50(Issue 7) pp:2287-2290
Publication Date(Web):May 2011
DOI:10.1016/j.commatsci.2011.01.043
From a statistical manner, we collected and correlated experimental bulk (B), shear (G), Young’s modulus (E), and ductility (G/B) with Vickers hardness (Hv) for a number of covalent materials and fitted quantitative and simple Hv–G and Hv–E relationships. Using these experimental formulas and our first-principles calculations, we further predicted the microhardness of some novel potential hard/superhard covalent compounds (BC2N, AlMgB14, TiO2, ReC, and PtN2). It was found that none of them are superhard materials (Hv ⩾ 40 GPa) except BC2N. The present empirical formula builds up a bridge between Vickers hardness and first-principles calculations that is useful to evaluate and design promising hard/superhard materials.Research highlights► Elastic moduli are correlated with Vickers hardness for covalent materials. ► Linear relationships between shear/Young’s modulus and Vickers hardness are fitted. ► Hardness of some potential superhard compounds is predicted.
Co-reporter:Pengbo Zhang, Jijun Zhao, Ying Qin, Bin Wen
Journal of Nuclear Materials 2011 Volume 419(1–3) pp:1-8
Publication Date(Web):December 2011
DOI:10.1016/j.jnucmat.2011.08.023
We report the energetics, stability, and diffusion behavior of helium (He), vacancies (V), and helium–vacancy complex clusters HenVm (n, m = 0–4) in vanadium solid from first-principles calculations. For He, vacancy site is more energetically favorable than tetrahedral interstitial by ∼0.74 eV, while hydrogen always prefers to stay in tetrahedral sites in vanadium. He exhibits a low migration energy (0.06 eV) and can be easily trapped in vacancy. A nearly linear relationship between formation energy and the number of He or vacancy is obtained for He or vacancy clusters, and the weak binding energies of He clusters indicate that He clusters themselves are unstable. The binding energies and dissociation energies of He and vacancy to helium–vacancy complex clusters are computed and compared well with the experimental observation from helium desorption spectra. The cluster stability depends on He content. Finally, He diffusion coefficients are predicted to be (1.07–1.27) × 10−8 m2 s−1 at typical temperatures of 600–800 K. We thus propose that He aggregation via vacancy trapping should be the main mechanism for He bubble formation.Highlights► He prefers to stay in vacancy and has a fast diffusion rate in vanadium. ► He can stabilize helium-vacancy clusters by suppressing vacancy emission. ► He aggregation by vacancy is a main mechanism for He bubble formation. ► An optimal ratio of helium to vacancy for the most stable cluster is revealed. ► Theoretical results are confirmed by helium desorption spectra observation.
Co-reporter:Pengbo Zhang, Jijun Zhao, Ying Qin, Bin Wen
Journal of Nuclear Materials 2011 Volume 413(Issue 2) pp:90-94
Publication Date(Web):15 June 2011
DOI:10.1016/j.jnucmat.2011.03.031
The stability and migration behavior of helium and self defects in vanadium and V–4Cr–4Ti alloy are studied by first-principles calculations. The tetrahedral site is found as the most stable configuration for interstitial He, followed by the octahedral and substitutional sites. Among the self defects, the monovacancy has lower formation energy (1.71 eV for V and 2.14 eV for V–4Cr–4Ti alloy) than the self interstitial ones. The migration energies for He hopping between the tetrahedral sites are 0.06 and 0.09 eV for vanadium and V–4Cr–4Ti alloy, respectively. Our calculations reveal strong repulsion between two interstitial He atoms and strong attraction between He and vacancy, suggesting that vacancy acts as a trapping site for He impurity and a seed for further bubble formation.
Co-reporter:Xue Jiang, Chunqiang Zhuang, Xiaoqing Li, Lingwei Sai, Jijun Zhao, Xin Jiang
Diamond and Related Materials 2011 Volume 20(Issue 7) pp:891-895
Publication Date(Web):July 2011
DOI:10.1016/j.diamond.2011.05.001
Based on first-principles calculations, we present the distributions of mechanical properties and formation ability of amorphous BxCyNz solids on the ternary B–C–N phase diagram. Along the C–BN isoelectronic line, the formation energy shows symmetric distributions; the Young's modulus and ratio of bulk modulus and shear modulus (B/G) show zonal distributions. Amazingly, for some peculiar compositions (B: 13–17 at.%; C: 48–52 at.%; N: 33–35 at.%), B–C–N solids exhibit certain ductile characteristic that is comparable to metals. On the phase area (B: 15–30 at.%; C: 50–60 at.%; N: 20–30 at.%), B–C–N solids possess both excellent hardness and good formation ability. These theoretical results provide valuable guidance for intentionally synthesizing BxCyNz materials with desirable mechanical properties.Research highlights► B–C–N coatings as potential industrial applications. ► Relationship of mechanical property and composition has not been built. ► We studied the relationship. ► We evaluated the formation ability and mechanical property for different compositions.
Co-reporter:X.Q. Li, J.J. Zhao, J.C. Xu, X. Liu
Journal of Materials Science & Technology 2011 Volume 27(Issue 11) pp:1029-1033
Publication Date(Web):November 2011
DOI:10.1016/S1005-0302(11)60182-5
Co-reporter:Xiaoqing Li, Chong Zhang, Jijun Zhao, Börje Johansson
Computational Materials Science 2011 Volume 50(Issue 9) pp:2727-2731
Publication Date(Web):July 2011
DOI:10.1016/j.commatsci.2011.04.027
V–(4–5) wt.% Cr–(4–5) wt.% Ti alloys are important candidate structural materials for the first-wall and blanket in future fusion reactor. Thus it is necessary to study the fundamental mechanical properties and the irradiation effects of the V-based alloys. Within a random solid solution model, the elastic constants and ideal strength of the V–4Cr–4Ti and the V–5Cr–5Ti alloys were calculated and compared with those of pure V solid. According to the theoretical Cauchy pressure and the ratio of bulk modulus and shear modulus, both alloys exhibit good ductility. Within the 250-atom supercell, inclusion of one vacancy defect or one interstitial H (He) atom will further enhance the ductility of these alloys.Highlights► V–4Cr–4Ti and V–5Cr–5Ti alloys exhibit better ductile behavior with regard to pure V solid. ► One vacancy or interstitial can soften bcc V–4Cr–4Ti and V–5Cr–5Ti alloys and enhance their ductility.
Co-reporter:Lingli Tang, Linwei Sai, Jijun Zhao, Ruifeng Qiu
Computational and Theoretical Chemistry 2011 Volume 969(1–3) pp:35-43
Publication Date(Web):30 August 2011
DOI:10.1016/j.comptc.2011.05.009
Classical and heptagon-contained nonclassical fullerene isomers of C30–C40 and C50 are generated by Spiral algorithm and optimized using DFT methods. There are much larger number of nonclassical isomers than classical ones. Energies of isomers generally increase with the number of heptagons inside and those lowest-energy nonclassical isomers contain only one heptagon. For all fullerenes considered, there is no nonclassical isomer that energetically prevails over the lowest-energy classical counterpart. But the best nonclassical isomers are more stable than most classical ones. The validity of pentagon adjacency penalty rule (PAPR) and hexagon neighbor rule (HNR) is examined for the nonclassical fullerenes. PAPR is still roughly applicable while HNR almost breaks down. A new topological index U with better performance over PAPR and HNR is presented and its superiority is further confirmed on two subsets of C50 isomers. This U index can be considered as an extension of PAPR to nonclassical isomers with existence of heptagons.Graphical abstractThe best structure of the heptagon-contained nonclassical C34 fullerene.Highlights► Classical and heptagon-contained nonclassical fullerene isomers of C30-C40 and C50 are studied. ► Pentagon adjacency penalty rule (PAPR) is roughly applicable for nonclassical fullerenes with heptagons. ► Hexagon neighbor rule (HNR) almost breaks down for the heptagon-contained nonclassical fullerenes. ► A new topological index U with better performance over PAPR and HNR is presented.
Co-reporter:Fengyu Li;Lu Wang;John Rui-Hua Xie
Theoretical Chemistry Accounts 2011 Volume 130( Issue 2-3) pp:341-352
Publication Date(Web):2011 October
DOI:10.1007/s00214-011-0989-6
The geometric structures, stabilization energies, dipole moments, and vibrational frequencies of the neutral water clusters (H2O)n, with n = 1–10, were investigated using density functional theory along with a variety of exchange-correlation functionals (LDA with SVWN5 parameterization, GGA with BLYP, PW91, PBE, B3LYP, X3LYP, PBE0, PBE1W, M05-2X, M06-2X and M06-L parameterizations) as well as high-level ab initio MP2 and CCSD(T) methods. Using the MP2 and CCSD(T) results as benchmarks, the effects of exchange-correlation functionals and basis sets were carefully examined. Each functional has its advantage in certain aspects; for example, M05-2X and X3LYP yield better geometries, and the capability of these two functionals to distinguish the relative energies between isomers are more similar to MP2. The size of the split-valence basis set (6-31G or larger), diffuse functions on the oxygen atom, and d(p) polarization on the oxygen (hydrogen) atom are crucial for an accurate description of intermolecular interaction in water clusters. The 6-31+G(2d,p) basis set is thus recommended as a compromise between computational efficiency and accuracy for structural description. We further demonstrated that the numerical basis set, TNP, performs satisfactorily in describing structural parameters of water clusters.
Co-reporter:Haili Gao ; Lu Wang ; Jijun Zhao ; Feng Ding ;Jianping Lu
The Journal of Physical Chemistry C 2011 Volume 115(Issue 8) pp:3236-3242
Publication Date(Web):February 10, 2011
DOI:10.1021/jp1094454
The electronic states of partially hydrogenated graphene (HG) structures are studied by the density functional theory calculations. Several types of HG configurations, including randomly removing of H pair, randomly removing individual H atoms, and ordered H pairs removal, are investigated. We find that the configurations with randomly removing H pairs are most energetically favorable. More interestingly, the band gap for such configurations decrease with H concentration and approaches zero around 67% H coverage. The ability to continuously tune the band gap of hydrogenated graphene from 0 to 4.66 eV by different H coverage provides a new pathway for engineering the electronic structure of graphene materials and enhances their applications in electronics and photonics.
Co-reporter:Junfeng Gao ; Qinghong Yuan ; Hong Hu ; Jijun Zhao ;Feng Ding
The Journal of Physical Chemistry C 2011 Volume 115(Issue 36) pp:17695-17703
Publication Date(Web):August 11, 2011
DOI:10.1021/jp2051454
To understand the nucleation of carbon atoms to form graphene on transition metal substrates during chemical vapor deposition (CVD) synthesis, carbon clusters supported on Ni(111) surfaces, namely CN@Ni(111) (where N ≤ 24), were explored systematically using density functional theory (DFT) calculations. Very different from the freestanding C clusters, on a Ni(111) surface, the C chain configuration is superior to the C ring formation and dominates the ground state until N > 12. A ground state structural transition from a one-dimensional C chain to a two-dimensional sp2 network (or graphene island) occurs at N = 12. It is surprising that incorporating one to three 5-membered-rings (5MRs) or pentagons into a graphene island is required to achieve the energetically most stable structure. This deep insight into the supported C cluster formation is crucial for understanding the growth mechanism of graphene on a transition metal surfaces in CVD experiments and the experimental design of CVD graphene growth.
Co-reporter:Fuchao Jia, Chunqiang Zhuang, Changyu Guan, Jijun Zhao, Yizhen Bai, Xin Jiang
Vacuum 2011 Volume 85(Issue 9) pp:887-891
Publication Date(Web):25 February 2011
DOI:10.1016/j.vacuum.2011.01.004
Amorphous B–C–N films were fabricated on silicon (100) substrates using radio frequency reactive magnetron sputtering technique with variable N2/Ar flow ratios. The structures, chemical bonding and mechanical properties were characterized by Fourier transform infrared spectroscopy, X-ray photoelectron spectroscopy and nanoindentation. We found that the N concentration is insensitive to the increment of N2/Ar flow ratio while the B concentration decreases and C concentration increases. All B–C, C–N, and B–N bond contents increase as the N2/Ar flow ratio varies from 1/10 to 5/10. Further improving the N2/Ar flow ratio will promote N atoms prior to bonding with C, resulting in decreased B–C and increased C–N bond content, respectively. The changes of bond content lead to a shift of the main peaks of B1s, C1s, and N1s spectra toward higher binding energies. The hardness of B–C–N thin films is almost invariant with the N2/Ar flow ratio.Research highlights► Bond contents can give semi-quantitative estimation of bonding level. ► Bond contents can greatly affect mechanical properties of a material. ► B-C-N coatings were produced by magnetron sputtering technique. ► Relation of bond content and mechanical properties in B-C-N films was studied.
Co-reporter:Fen Li, Jijun Zhao, Börje Johansson, Lixian Sun
International Journal of Hydrogen Energy 2010 Volume 35(Issue 1) pp:266-271
Publication Date(Web):January 2010
DOI:10.1016/j.ijhydene.2009.10.061
We proposed a possible way of promoting the binding of H2 molecules on covalent organic frameworks crystals via substituting the bridge C2O2B rings with different metal-participated rings, which can naturally avoid the clustering of metal atoms. First-principles calculations on both crystalline phase and molecular fragments show that the H2 binding energy can be enhanced by a factor of four with regard to the undoped crystal, i.e. reaching about 10 kJ/mol. Grand canonical Monte Carlo simulations further confirm that such substitutional doping would improve the room temperature hydrogen storage capacity by a factor of two to three.
Co-reporter:Jijun Zhao, Chunqiang Zhuang, Xin Jiang
Diamond and Related Materials 2010 Volume 19(Issue 11) pp:1419-1422
Publication Date(Web):November 2010
DOI:10.1016/j.diamond.2010.08.010
Within a random solid solution model based on diamond lattice, the formation energies, lattice constants, and elastic properties of BC2N crystals were calculated using first-principles methods. We found that the degree of BN–C2 mixture defined by the ratio of bond contents is the key for the lattice parameter, relative stability, and physical properties in the BC2N solid solution. BC2N crystals with the most probable degree of mixture possess larger lattice parameters and less elastic moduli with respect to cubic-BN, and the high formation energy can be compensated by the contribution of configuration entropy under high-temperature conditions. The present model can explain the controversial experimental results on BC2N crystals synthesized under different conditions.► Random solid solution model for BC2N crystals. ► Degree of mixture as the key parameter. ► High formation energy can be compensated by the configurational entropy at high Tpatients.
Co-reporter:Lu Wang, Jijun Zhao, Fengyu Li, Zhongfang Chen
Chemical Physics Letters 2010 Volume 501(1–3) pp:16-19
Publication Date(Web):6 December 2010
DOI:10.1016/j.cplett.2010.10.052
Abstract
Unbiased search by first-principles simulated annealing revealed irregular cage configurations for medium-sized Bn clusters (namely, boron fullerene) with n = 32–56, which are more stable than the previously proposed symmetric cages. The stability of these irregular cages can be understood by the three-center bonds as well as the polygonal holes on the cage surface. The delocalized distribution of molecular orbitals as well as the negative nucleus-independent chemical shifts (NICS) values for each boron cage indicates strong aromaticity.
Co-reporter:Tao Xue, Jing Luo, Si Shen, Fengyu Li, Jijun Zhao
Chemical Physics Letters 2010 Volume 485(1–3) pp:26-30
Publication Date(Web):18 January 2010
DOI:10.1016/j.cplett.2009.12.019
Co-reporter:Junfeng Zhang, Fengyu Li, Xiangyang Miao, Jijun Zhao, Long Jing, Guohui Yang, Xiangfu Jia
Chemical Physics Letters 2010 Volume 492(1–3) pp:68-70
Publication Date(Web):26 May 2010
DOI:10.1016/j.cplett.2010.04.007
Abstract
The most stable structures of polyhydroxylated metallofullerene Gd@C82(OH)x (x = 1–24) were investigated using density functional theory calculations. After examining a number of structural isomers, some rules about the locations of hydroxyl groups on the endohedral fullerene cage of Gd@C82 were proposed. The hydroxylated carbon atoms tend to enclose two opposite six-membered rings for the most stable Gd@C82(OH)12 and seven six-membered rings for the most stable Gd@C82(OH)24, respectively. The hydroxyl groups prefer to locate on para-position of the six-membered ring or meta-position of the five-membered ring.
Co-reporter:Xue Jiang, Jijun Zhao, Chunqiang Zhuang, Bin Wen, Xin Jiang
Diamond and Related Materials 2010 Volume 19(Issue 1) pp:21-25
Publication Date(Web):January 2010
DOI:10.1016/j.diamond.2009.10.011
Mechanical and electronic properties of ultrathin hydrogenated nanodiamonds (with diameters from 0.71 nm to 1.4 nm) under uniaxial compression have been investigated by means of density functional theory calculations. The computed Young's moduli of nanodiamonds are lower than the bulk value and increase with size, which can be fitted to an empirical function of diameter. Similar to the bulk diamond, the HOMO–LUMO gaps of nanodiamond reduces under uniaxial strain, implying tunable electronic properties via mechanical deformations.
Co-reporter:Zhao Guo, Bin Lu, Xue Jiang, Jijun Zhao, Rui-Hua Xie
Physica E: Low-dimensional Systems and Nanostructures 2010 Volume 42(Issue 5) pp:1755-1762
Publication Date(Web):March 2010
DOI:10.1016/j.physe.2010.01.039
The lowest-energy structures of medium-sized Lin (n=20, 30, 40, 50) clusters are determined from simulated annealing technique followed by geometry optimization within the framework of density functional theory. The shapes of magic-number Li20 and Li40 clusters are nearly spherical while those of the other clusters are ellipsoid. The growth of Lin clusters is based on core of multi-layered pentagonal bipyramids with other atoms capped on the surface. The binding energies of the Lin clusters were computed and compared with experiments. At magic-number sizes (n=20, 40), Lin clusters possess relatively larger HOMO–LUMO gaps and higher ionization potentials, corresponding to the closure of electron shell. The molecular orbitals of the lithium clusters can be grouped into electron shells and their spatial distributions resemble the atomic orbitals. The average polarizability of the Li clusters reduces rapidly with cluster size and can be approximately described by a classical metallic sphere model. The optical absorption spectra of Lin clusters from time-dependent density functional theory calculations show giant resonance phenomenon, and resonance peak blueshifts with increasing cluster size.
Co-reporter:Baolin Wang, Xiaoqiu Wang and Jijun Zhao
The Journal of Physical Chemistry C 2010 Volume 114(Issue 13) pp:5741-5744
Publication Date(Web):March 5, 2010
DOI:10.1021/jp908472h
Very recently, mass spectroscopy of ZnO clusters revealed a hitherto unknown (ZnO)60 magic-number cluster with exceptional stability. Using first-principles approaches, we searched the most stable structures of medium-size (ZnO)n clusters by considering several possible structural motifs. Instead of the previously nominated nested cage for (ZnO)60, we found a sodalite structure via coalescence of (ZnO)12 cages, which was predicted to be a metastable phase in bulk ZnO solid. Due to the smaller influence of surface reconstruction, this sodalite motif is very competitive for larger (ZnO)n clusters up to n = 96.
Co-reporter:Jijun Zhao, Lu Wang, Fengyu Li, and Zhongfang Chen
The Journal of Physical Chemistry A 2010 Volume 114(Issue 37) pp:9969-9972
Publication Date(Web):August 10, 2010
DOI:10.1021/jp1018873
Unbiased search on the potential energy surface of medium-sized boron clusters, with B80, B74, and B68 as representatives, was carried out using simulated annealing incorporated with first-principles molecular dynamics. These boron clusters thermodynamically prefer the B12-centered core−shell structures, which resemble the fragment of bulk boron solids. Though these core−shell clusters lack a descriptive symmetry and may not be the true global minima, the core−shell B80 is about 25 meV/atom lower in energy than the widely assumed highly stable “magic” B80 fullerene. The electronic states and photoelectron spectra of these clusters are closely correlated to the structural motif, which may be helpful for detecting the cluster configurations in experiments.
Co-reporter:Lu Wang, Jijun Zhao, Fengyu Li, Haiping Fang and Jian Ping Lu
The Journal of Physical Chemistry C 2009 Volume 113(Issue 14) pp:5368-5375
Publication Date(Web):2017-2-22
DOI:10.1021/jp808873r
Water molecules confined inside a single-walled (6, 6) carbon nanotube were investigated using density functional theory. In this narrow-sized carbon nanotube (of about 0.8 nm in diameter), the encapsulated water molecules form chain-like configurations via hydrogen bonding. As compared to the water chains in vacuum, the intramolecular charge transfer in the encapsulated water chain is enhanced and the dipole moment is reduced due to the screening effect of the carbon nanotube. The tube−molecule interaction characterized by the coupling energy is about 0.28 eV per water molecule by local density approximation and 0.1 eV by general gradient approximation; the latter one is close to the results by empirical potentials. Weak coupling between the molecular orbitals of the encapsulated water molecules and the delocalized π electrons from the carbon nanotube was observed, implying that the tube−water interaction is not a simple effect of geometry confinement. Vibrational analysis revealed some unique hydrogen-bond modes for the water chains as well as red shift of the O−H stretching modes for the encapsulated water molecules with regard to the vacuum frequencies due to the tube−water interaction.
Co-reporter:Jianguang Wang, Lu Wang, Li Ma, Jijun Zhao, Baolin Wang, Guanghou Wang
Physica E: Low-dimensional Systems and Nanostructures 2009 Volume 41(Issue 5) pp:838-842
Publication Date(Web):March 2009
DOI:10.1016/j.physe.2008.12.018
Structural and electronic properties of single-walled TiO2 nanotubes were investigated using density-functional theory. The strain energies were computed as a function of tube diameter and the small-sized nanotubes are less stable than the large-diameter ones. All TiO2 nanotubes are semiconductors with indirect bandgaps independent of the tube chirality. The bandgaps of those nanotubes slightly increase with increase in tube diameter, and they are higher than that of the bulk TiO2 solid. The hydrogen-storage capacity for a selected TiO2 nanotube is 3.2 wt% with the binding energy of 0.05 eV per H2 molecule.
Co-reporter:Chunqiang Zhuang, Xue Jiang, Jijun Zhao, Bin Wen, Xin Jiang
Physica E: Low-dimensional Systems and Nanostructures 2009 Volume 41(Issue 8) pp:1427-1432
Publication Date(Web):August 2009
DOI:10.1016/j.physe.2009.04.011
The heat of formation and infrared spectra of hydrogenated nanodiamonds with different morphologies (up to 1.92 nm in diameter) have been investigated using density functional theory. Preferential growth orientation along 〈1 1 1〉 orientations corresponding to an octahedral shape was observed according to the computed heat of formation. The simulated infrared (IR) spectra show distinct dependence on both size and morphology of the nanodiamonds. The major IR peaks in the region of 900–1300 cm−1 shift toward high wavenumber as the nanocluster size increases, while the C–H vibrations in the range 2500–3500 cm−1 show an opposite trend. The present results would be useful for identifying the size and morphology of nanodiamonds in experimental IR spectroscopy.
Co-reporter:Chunqiang Zhuang, Jijun Zhao, Fuchao Jia, Changyu Guan, Zhanling Wu, Yizhen Bai, Xin Jiang
Surface and Coatings Technology 2009 204(5) pp: 713-717
Publication Date(Web):
DOI:10.1016/j.surfcoat.2009.09.031
Co-reporter:Lei Liu, Jianguo Du, Jijun Zhao, Hong Liu, Di Wu, Fuliang Zhao
Computer Physics Communications 2008 Volume 179(Issue 6) pp:417-423
Publication Date(Web):15 September 2008
DOI:10.1016/j.cpc.2008.04.003
The high-pressure and/or high-temperature behavior and α/βα/β phase transition of forsterite (Mg2SiO4) were investigated by first principle calculations: local density approximation (LDA), generalized gradient approximation (GGA) and the quasi-harmonic Debye model in consideration of the temperature effect. At zero pressure and zero temperature, the calculated values of lattice parameters and elastic constants are well concordant with experimental data. The calculated values of cell volume at different pressures and temperatures are compatible with experiments. The results calculated by LDA are more concordant with the experimental results than those by GGA for describing Mg2SiO4 crystals. The boundaries of the P–TP–T plots for α/βα/β phase transition are calculated to be 0.00779 GPa/K by LDA and 0.00772 GPa/K by GGA, suggesting the upper and lower limits of the phase boundary for α/βα/β forsterite. At temperature of less than 1000 K, the calculated results agree reasonably with experimental ones. Thus, the present results indicate that combination of first-principles and quasi-harmonic Debye model is an efficient approach to simulate the behavior of minerals at high pressure and/or high temperature.
Co-reporter:Li Ma, Jianguang Wang, Jijun Zhao, Guanghou Wang
Chemical Physics Letters 2008 Volume 452(1–3) pp:183-187
Publication Date(Web):4 February 2008
DOI:10.1016/j.cplett.2007.12.054
To discern the anisotropy of silicon nanowires (SiNWs) grown in different directions, the binding energy, heat of formation, and Young’s modulus of hydrogen-passivated SiNWs with various diameters and crystallographic orientations were calculated using all-electron density functional theory. In the size range studied, nanowires grown in the [1 1 0] direction are most stable while those in the [1 0 0] direction are energetically least favorable. Similar trend was observed in the computed Young’s modulus. With the same radius, the nanowire along the [1 1 0] direction possesses the highest Young’s modulus, while the [1 0 0] wire has the lowest value.The optimized SiNWS structures and Young’s modulus in four different growth directions as a function of wire mean diameter.
Co-reporter:Jijun Zhao, Hong Liu
Computational Materials Science 2008 Volume 42(Issue 4) pp:698-703
Publication Date(Web):June 2008
DOI:10.1016/j.commatsci.2007.10.008
1,1-Diamino-2,2-dinitroethylene (C2H4N4O4, FOX-7) is a newly found energetic materials with high performance and low sensitivity. We performed density functional theory (DFT) calculations within local density approximation (LDA) as well as generalized gradient approximation (GGA) to simulate the structural and electronic properties of FOX-7 crystal under high-pressure up to 4 GPa. Due to the insufficient treatment of weakly intermolecular interaction, LDA underestimates the lattice parameters and volume, while GGA overestimates them. It is interesting to find that a mixing of LDA and GGA results of lattice properties can lead to fairly good agreement with experiments. Examination of pressure-induced changes of molecular geometry shows that the C–N bonds are most sensitive under compression. The change of band gap as function of external pressure was also discussed.
Co-reporter:Bin Wen, Michael J. Bucknum, Jijun Zhao, Xu Guo, Tingju Li
Diamond and Related Materials 2008 Volume 17(7–10) pp:1353-1355
Publication Date(Web):July–October 2008
DOI:10.1016/j.diamond.2008.01.072
Diamond is the hardest material with superior mechanical stability. Due to the widespread application in high pressure research, stability limit of diamond at ultrahigh pressure is an importance issue for theoretical and experimental studies, which remains controversy between theoretical prediction and experimental observations. In this paper, we report first-principle calculations on the relative stability of diamond under non-hydrostatic compression and the computational results demonstrate that the cubic diamond would become unstable with respect to the lonsdalite under isotropic stress tensor with mean stress as low as 26 GPa, while agree well with experimental findings under similar non-hydrostatic loading conditions.
Co-reporter:Bin Wen, Jijun Zhao, Fudong Bai, Tingju Li
Intermetallics 2008 Volume 16(Issue 2) pp:333-339
Publication Date(Web):February 2008
DOI:10.1016/j.intermet.2007.11.003
The structural properties, heat of formation, elastic properties, and electronic structures of six kinds of binary Al–Ru intermetallics were studied using first-principle calculations. The structural properties of these intermetallic compounds agree well with previously reported experimental data except that the computed equilibrium lattice parameters of the orthorhombic Al5Ru2 crystal (Cmcm) are significantly different from the experimental data reported by Mi et al. [Intermetallics 2003;11:643]. With increasing Ru concentration, the bulk modulus of Al–Ru intermetallics increases linearly. Most of the Al–Ru alloys considered are conductors, while the Al2Ru alloy is an indirect band gap semiconductor with a band gap of 0.168 eV.
Co-reporter:Lu Wang, Jijun Zhao
Journal of Molecular Structure: THEOCHEM 2008 Volume 862(1–3) pp:133-137
Publication Date(Web):15 August 2008
DOI:10.1016/j.theochem.2008.05.011
Lowest-energy structures of InnPn (n = 1–12) clusters have been determined from a number of structural isomers using all-electron density functional calculations. The In–P alternating hollow cage-like structures emerge at n = 7. Size-dependent cluster properties such as binding energy, HOMO–LUMO gaps, electron affinities, and photoelectron spectra have been computed and discussed. The simulated electron affinities and photoelectron spectra agree reasonably with experiments. With exceptionally low electron affinity and large HOMO–LUMO gap, In3P3 was identified as a magic-numbered cluster of relatively higher stability, in agreement with experimental observations.
Co-reporter:Lu Wang ; Jijun Zhao ;Haiping Fang
The Journal of Physical Chemistry C 2008 Volume 112(Issue 31) pp:11779-11785
Publication Date(Web):July 10, 2008
DOI:10.1021/jp8048185
Encapsulation of (H2O)n clusters (n = 1−22) in fullerene cages of different diameters (0.73−1.41 nm) has been investigated using gradient-corrected density functional theory. A linear relationship between cavity volume and maximum number of the encapsulated water molecules has been obtained. The interaction between water molecules and the fullerene wall was identified as physisorption with an adsorption energy of about 1.1 kcal/mol per molecule. The equilibrium configurations of small confined water clusters (n < 12) roughly resemble those of gas-phase clusters, whereas larger water clusters tend to adopt cage-like configurations when they are encapsulated in fullerene cages of sufficiently large diameter (i.e., >1.4 nm). The dipole moments of water clusters in the confined phase are smaller than those in the gas phase due to the screening effect of the outer fullerene cage. These results might shed some light on the behavior of water clusters confined in the nonpolar cavities of biological interests.
Co-reporter:Jijun Zhao, Lu Wang, Jianming Jia, Xiaoshuang Chen, Xiaolan Zhou, Wei Lu
Chemical Physics Letters 2007 Volume 443(1–3) pp:29-33
Publication Date(Web):27 July 2007
DOI:10.1016/j.cplett.2007.06.055
The lowest-energy structures of AlnPn clusters up to n = 9 have been explored using all-electron density functional calculations with a gradient correction. For smaller AlnPn clusters with n = 1–4, we successfully reproduced the previously reported lowest-energy structures. Novel cage structures with Al–P alternating arrangement were observed for n ⩾ 5. The comparison of the lowest-energy structures of AlnPn clusters with those of Si2n, BnNn, and GanAsn clusters has been made. Size-dependent cluster properties such as binding energy, HOMO–LUMO gaps, electron affinities, and photoelectron spectra have been computed and compared with experiments.The lowest-energy structures of AlnPn clusters up to n = 9 are obtained from density functional theory calculations.
Co-reporter:Bin Wen, Jijun Zhao, Tingju Li
Chemical Physics Letters 2007 Volume 441(4–6) pp:318-321
Publication Date(Web):25 June 2007
DOI:10.1016/j.cplett.2007.05.054
Density functional theory calculations have been performed to compute the heat of formation of nanoscale graphene sheets and nanodiamonds with up to ∼450 carbon atoms. The accuracy of the computational scheme has been validated by calculating the heat of reaction of three precursor reactions between small sp2 and sp3 hydrocarbon molecules. By comparing the heat of formation for nanographites and nanodiamonds of different sizes, an sp3-to-sp2 crossover was estimated to occur at about n = 1060 (d ≃ 2.1 nm). The size-dependence of the structural and electronic properties of hydrogenated nanodiamonds has also been discussed.Nanodiamonds of octahedron shape with Td symmetry up to C455H196 are studied using density function theory.
Co-reporter:Dongxu Tian, Hualei Zhang, Jijun Zhao
Solid State Communications 2007 Volume 144(3–4) pp:174-179
Publication Date(Web):October 2007
DOI:10.1016/j.ssc.2007.05.020
Using a genetic algorithm followed by local optimization with density functional theory, the lowest-energy structures of Agn clusters in a size range of n=3–22n=3–22 were studied. The Agn (n=9–16n=9–16) clusters prefer compact structures of flat shape, while the Agn(n=19,21,22) clusters adopt amorphous packing based on a 13-atom icosahedral core. For Ag16, two competitive candidates for the lowest-energy structures, namely a hollow-cage structure and close-packed structures of flat shape, were found. Two competing candidates were found for Ag17 and Ag18: hollow-cage structures versus icosahedron-based compact structures. The lowest-energy structure of Ag20 is a highly symmetric tetrahedron with TdTd symmetry. These results are significantly different from those predicted in earlier works using empirical methods. The ionization potentials and electron affinities for the lowest-energy structures of Agn (n=3–22n=3–22) clusters were computed and compared with experimental values.
Co-reporter:R. Bruce King and Jijun Zhao
Chemical Communications 2006 (Issue 40) pp:4204-4205
Publication Date(Web):01 Sep 2006
DOI:10.1039/B607895H
The valence electrons in the recently reported icosahedral cluster [As@Ni12@As20]3− with a Russian matryoshka nesting doll structure can be partitioned so that the central As atom has the rare gas configuration, as As3−, and the intermediate Ni12 icosahedron receives 40 electrons from the lone pairs of the outer As20 dodecahedron to be isoelectronic with the Al13− jellium cluster found in molecular beam experiments.
Co-reporter:Jijun Zhao, Xiaoshuang Chen, John R.H. Xie
Analytica Chimica Acta 2006 Volume 568(1–2) pp:161-170
Publication Date(Web):24 May 2006
DOI:10.1016/j.aca.2006.02.006
Chemical doping of carbon nanotubes provides a variety of opportunities for tailoring the physical properties of carbon nanotubes. In this review, we discussed the optical properties of doped carbon nanotubes and the related applications as nanoscale photonic devices. The fundamental optical properties of carbon nanotubes with various chemical doping have been summarized. Novel optoelectronic and photonic devices based on doped carbon nanotubes, such as optical nonlinear materials, optical limiting devices, photovoltaic devices, etc., have been discussed.
Co-reporter:Jijun Zhao, Baolin Wang, Xiaolan Zhou, Xiaoshuang Chen, Wei Lu
Chemical Physics Letters 2006 Volume 422(1–3) pp:170-173
Publication Date(Web):28 April 2006
DOI:10.1016/j.cplett.2006.02.048
Abstract
The lowest-energy structures were obtained for GanNn clusters (n = 4–12) using gradient-corrected density functional theory. For each cluster size, a number of structural isomers were constructed and optimized. Crossover from ring to cage-like structures were found at Ga8N8. The observed cage-like configurations consisting of six-membered rings with gallium–nitrogen alternative arrangement can be viewed as embryo of wurtzite structure of gallium nitride crystal. The size-dependence of cluster properties such as binding energy, HOMO–LUMO gaps, ionization potentials, electron affinities, and Mulliken charges have been computed and analyzed. The bonding characteristics of the GanNn clusters were discussed.
Co-reporter:Dongxu Tian, Hualei Zhang, Jijun Zhao
Solid State Communications (October 2007) Volume 144(3–4) pp:174-179
Publication Date(Web):1 October 2007
DOI:10.1016/j.ssc.2007.05.020
Using a genetic algorithm followed by local optimization with density functional theory, the lowest-energy structures of Agn clusters in a size range of n=3–22 were studied. The Agn (n=9–16) clusters prefer compact structures of flat shape, while the Agn(n=19,21,22) clusters adopt amorphous packing based on a 13-atom icosahedral core. For Ag16, two competitive candidates for the lowest-energy structures, namely a hollow-cage structure and close-packed structures of flat shape, were found. Two competing candidates were found for Ag17 and Ag18: hollow-cage structures versus icosahedron-based compact structures. The lowest-energy structure of Ag20 is a highly symmetric tetrahedron with Td symmetry. These results are significantly different from those predicted in earlier works using empirical methods. The ionization potentials and electron affinities for the lowest-energy structures of Agn (n=3–22) clusters were computed and compared with experimental values.
Co-reporter:Yuanyuan Wang, Jijun Zhao, Chi Zhang
Fusion Engineering and Design (May 2017) Volume 118() pp:129-134
Publication Date(Web):1 May 2017
DOI:10.1016/j.fusengdes.2017.03.145
•The initial internal variable in the Anand model is modified by considering both temperature and irradiation dose.•The tensile stress-strain response is examined and analyzed under different temperatures and irradiation doses.•Yield strengths are predicted as functions of strain rate, temperature and irradiation dose.The viscoplastic equations with a modified initial internal variable are implemented into the finite element code to investigate stress-strain response and irradiation hardening of the materials under increased temperature and at different levels of irradiated dose. We applied this model to Mod 9Cr-1Mo steel. The predicted results are validated by the experimentally measured data. Furthermore, they show good agreement with the previous data from a constitutive crystal plasticity model in account of dislocation and interstitial loops. Three previous hardening models for predicting the yield strength of the material are discussed and compared with our simulation results.
Co-reporter:Yingying Huang, Chongqin Zhu, Lu Wang, Jijun Zhao, Xiao Cheng Zeng
Chemical Physics Letters (March 2017) Volume 671() pp:
Publication Date(Web):March 2017
DOI:10.1016/j.cplett.2017.01.035
•A new cubic ice clathrate (s-IV) with a record low mass density is predicted.•The s-IV can be fully stabilized with twelve CH4 molecules encapsulated in each large cavity.•In the P-T phase diagram of water, the s-IV is a more stable phase in certain range of negative pressures.Using extensive Monte Carlo packing algorithm and dispersion-corrected density functional theory optimization, we predict a new cubic crystalline phase of ice clathrate, named as s-IV, which is composed of eight large icosihexahedral cavities (12464418), eight intermediate dodecahedral cavities (6646), and sixteen small octahedral cavities (6246) per unit cell. Based on DFT calculations, we find that the s-IV ice clathrate with an extremely low mass density of 0.506 g/cm3. In the P-T phase diagram of water described by the TIP4P/2005 water model, the s-IV ice clathrate becomes a more stable ice polymorph in the negative-pressure region, e.g., below –3830 bar at 0 K, below –4882 bar at 115 K, and below –7292 bar at 200 K.
Co-reporter:Xue Jiang, Jawad Nisar, Biswarup Pathak, Jijun Zhao, Rajeev Ahuja
Journal of Catalysis (March 2013) Volume 299() pp:204-209
Publication Date(Web):1 March 2013
DOI:10.1016/j.jcat.2012.12.022
To elucidate the usage of graphene oxide (GO) as a photocatalysis material, we have studied the effect of epoxy and hydroxyl functionalization on the electronic structure, work function, CBM/VBM position, and optical absorption spectra of GO using density functional theory calculations. By varying the coverage and relative ratio of the surface epoxy (O) and hydroxyl (OH) groups, both band gap and work function of the GO materials can be tuned to meet the requirement of photocatalyst. Interestingly, the electronic structures of GO materials with 40–50% (33–67%) coverage and OH:O ratio of 2:1 (1:1) are suitable for both reduction and oxidation reactions for water splitting. Among of these systems, the GO composition with 50% coverage and OH:O (1:1) ratio can be very promising materials for visible-light-driven photocatalyst. Our results not only explain the recent experimental observations about 2-D graphene oxide as promising visible-light-driven photocatalyst materials but can also be very helpful in designing the optimal composition for higher performance.Graphical abstractSite levels of VBM and CBM for OH:O = 1 and OH:O = 2 graphene oxide with different coverage rate: (a), (b), (c), (d), and (e) represented C36O8H4, C24O8H4, and C24O12H6 with OH:O = 1, C20O6H4 and C16O6H4 with OH:O = 2, respectively. The dot lines are standard water redox potentials. The reference potential is the vacuum level.Our results not only explain the recent experimental observations that graphene oxide is a promising two-dimensional material for visible-light photocatalysis but can be very helpful in designing the optimal composition for higher performance.Download high-res image (129KB)Download full-size imageHighlights► Show GO is a promising material for visible-light-driven photocatalyst. ► Explain recent observations. ► Helpful in designing GO.
Co-reporter:Pengbo Zhang, Jijun Zhao, Ying Qin, Bin Wen
Fusion Engineering and Design (January 2011) Volume 86(Issue 1) pp:45-50
Publication Date(Web):1 January 2011
DOI:10.1016/j.fusengdes.2010.08.003
To understand the combined effect of plasma heating and neutron heating loadings, the distributions of temperature, stress, and strain in different two-dimensional first wall panel models under normal ITER operation condition were simulated using finite element method. The maximum temperature occurs at the Be armor, and reaches 461 °C. High thermal stresses (in the range of 80–200 MPa) are found at the interface between the Be armor and the CuCrZr layer. The maximum thermal stress reaches 324 MPa in the SS316L cooling tube (20 mm diameter), exceeding its yield strength and resulting in a maximum strain of about 1.7% at the tube inner surface. These simulation results are useful for the design and operation of ITER.
Co-reporter:Jie Fu, Xiaoqing Li, Börje Johansson, Jijun Zhao
Journal of Alloys and Compounds (25 May 2017) Volume 705() pp:369-375
Publication Date(Web):25 May 2017
DOI:10.1016/j.jallcom.2017.02.103
Co-reporter:Hongsheng Liu, Nannan Han, Jijun Zhao
Applied Surface Science (1 July 2017) Volume 409() pp:
Publication Date(Web):1 July 2017
DOI:10.1016/j.apsusc.2017.03.007
•Chemical potential phase diagrams of silicene/Ag(111) at varied temperatures.•The priorities of various silicene phases in experiments are explained.•A proper experimental condition to obtain homogeneous 4 × 4 silicene is recommended.Silicene, the single layer of silicon atoms arranged in a honeycomb lattice, has been synthesized in recent experiments and attracted significant attentions. Silicene is promising in future nanoelectronic devices due to its outstanding electronic properties. In experiments, however, different silicene superstructures coexist on Ag(111) substrate. For the device applications, homogenous silicene sheet with large scale and high quality is highly desired. Here, for the first time, we investigate both the temperature and the coverage effects on the thermal stability of epitaxial silicene on Ag(111) surface by ab initio molecular dynamics simulations. The relationship between the stability of various silicene superstructures and the growth conditions, including temperature and coverage of silicon atoms, is revealed by plotting the chemical potential phase diagram of silicene on Ag(111) surfaces at different temperatures. Our results are helpful for understanding the observed diversity of silicene phases on Ag(111) surfaces and provide some useful guidance for the synthesis of homogenous silicene phase in experiments.
Co-reporter:Lu Wang, Jijun Zhao, Lili Wang, Tianyin Yan, Yi-Yang Sun and Shengbai B. Zhang
Physical Chemistry Chemical Physics 2011 - vol. 13(Issue 47) pp:NaN21131-21131
Publication Date(Web):2011/10/24
DOI:10.1039/C1CP21778J
We propose titanium-decorated graphene oxide (Ti–GO) as an ideal sorbent for carbon monoxide (CO) capture and separation from gas mixtures. Based on first-principles calculations, Ti–GO exhibits a large binding energy of ∼70 kJ mol−1 for CO molecules, while the binding energies for other gases, such as N2, CO2, and CH4, are significantly smaller. The gas adsorption properties of Ti–GO are independent of the local GO structures once Ti atoms are anchored by the oxygen-containing groups on the GO surface. The strong interaction between CO molecule and Ti is a result of dative bonding, i.e., hybridization between an empty d orbital of Ti and an occupied p orbital of CO. Adsorption isotherms from grand canonical Monte Carlo simulations clearly demonstrate the strong selectivity of Ti–GO for CO adsorption in a mixture with other gas.
Co-reporter:Lizhao Liu, Junfeng Gao, Xu Guo and Jijun Zhao
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 40) pp:NaN17141-17141
Publication Date(Web):2013/08/14
DOI:10.1039/C3CP52904E
Atomic structural models of zigzag-shaped carbon nanotubes (Z-CNTs) were constructed by periodically introducing pentagons and heptagons into pristine CNTs. In terms of formation energies, the Z-CNTs present comparable energetic stabilities to those of the pristine CNTs and are more stable than C60 fullerene. The mechanical properties of these Z-CNTs, including the Young's modulus, intrinsic strength and failure behaviour, were systematically investigated by first-principles computations. Compared with the pristine CNTs with an average Young's modulus of about 1.0 TPa, incorporation of pentagons and heptagons in the Z-CNTs will reduce the average Young's modulus to several hundreds of GPa. Moreover, the computational results also showed that under uniaxial tensile strain, the intrinsic strength and failure strain of the Z-CNTs are also lower than those of the pristine CNTs. Generally, the Young's modulus and intrinsic strength of the Z-CNTs are exponentially inverse to curvature, which can be expressed by simple formulae. In particular, the electronic properties of the armchair Z-CNTs can be tailored by uniaxial tensile strain. It was also found that through applying tensile strain, a semiconductor–metal or metal–semiconductor transition can be triggered. The localized–delocalized partial charge distribution near the Fermi energy for the strained Z-CNTs can explain the semiconductor–metal or metal–semiconductor transition. This significant electromechanical coupling effect suggests the Z-CNTs have potential applications in nanoscale electromechanical sensors and switches.
Co-reporter:Han Zhang, Zheng Duan, Xiaonan Zhang, Chao Liu, Junfeng Zhang and Jijun Zhao
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 28) pp:NaN11799-11799
Publication Date(Web):2013/05/15
DOI:10.1039/C3CP44716B
We present a molecular dynamic simulation on the mechanical strength and fracture behavior of graphene grain boundaries (GBs). The intrinsic strength, critical failure strain, and failure mechanism of graphene GBs mainly rely on the temperature and inflection angle, whereas the Young's modulus does not vary significantly with either temperature or boundary configuration. The overall intrinsic strengths of inflected GBs can be correlated with infection angle by a linear term, which is irrelevant to the system temperature. The initial failure sites of GBs locate either on the boundary line or inside the domain at high temperature.
Co-reporter:Si Zhou, Yu Guo and Jijun Zhao
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 15) pp:NaN10615-10615
Publication Date(Web):2016/03/17
DOI:10.1039/C6CP01012A
The thermoelectric properties of two-dimensional (2D) materials are of great interest for both fundamental science and device applications. Graphene oxide (GO), whose physical properties are highly tailorable by chemical and structural modifications, is a potential 2D thermoelectric material. In this report, we pattern nanoroads on GO sheets with epoxide functionalization, and investigate their ballistic thermoelectric transport properties based on density functional theory and the nonequilibrium Green's function method. These graphene oxide nanoroads (GONRDs) are all semiconductors with their band gaps tunable by the road width, edge orientation, and the structure of the GO matrix. These nanostructures show appreciable electrical conductance at certain doping levels and enhanced thermopower of 127–287 μV K−1, yielding a power factor 4–22 times of the graphene value; meanwhile, the lattice thermal conductance is remarkably reduced to 15–22% of the graphene value; consequently, attaining the figure of merit of 0.05–0.75. Our theoretical results are not only helpful for understanding the thermoelectric properties of graphene and its derivatives, but also would guide the theoretical design and experimental fabrication of graphene-based thermoelectric devices of high performance.
Co-reporter:Tengfei Cao, Da Wang, Dong-Sheng Geng, Li-Min Liu and Jijun Zhao
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 10) pp:NaN7162-7162
Publication Date(Web):2016/01/27
DOI:10.1039/C5CP06528C
The indirect bandgap character of silicon greatly limits its applications in electronic or optoelectronic devices, and direct bandgaps are highly desirable in all silicon allotropes. The successful synthesis of ultrathin or even monolayer silicon films experimentally has opened new opportunities to further modulate the electronic structure of silicon through external modulation. In this work, strain or electric field effects on the electronic structure of ultrathin silicon film (USF) are systematically explored. The results demonstrate that all USFs are indirect band-gap semiconductors; interestingly, tensile strain or electric field efficiently tunes the USFs into direct band gap semiconductors. The indirect to direct band gap transition in the USFs not only extends their light adsorption spectra into the visible light region but also greatly enhances the adsorption intensity. Because of this, strained USFs have great potential to be used as a high-performance photovoltaic material. Furthermore, the high stability, moderate band-gap and proper band edge positions demonstrate that monolayer and bilayer USFs can also be used as photocatalysts for water splitting.
Co-reporter:Fen Li, Junfeng Gao, Jian Zhang, Fen Xu, Jijun Zhao and Lixian Sun
Journal of Materials Chemistry A 2013 - vol. 1(Issue 27) pp:NaN8022-8022
Publication Date(Web):2013/05/01
DOI:10.1039/C3TA10800G
The incorporation of lithium amidoborane (LiAB) into graphene oxide (GO) and the dehydrogenation process of the GO–LiAB complex have been investigated for combining the chemical and physical hydrogen storage approaches. The obtained adsorption energy and minimum energy pathway (MEP) demonstrate that both of the two dominant groups of –O– and OH– contribute to the facile combination of GO and LiAB (GO3–LiAB). The GO3–LiAB complex has a better dehydrogenation performance than the pristine LiAB, which also indicates the feasibility of doping Li atoms on the GO surface. Atomic charge and bond length analyses match well with the MEP prediction. By dehydrogenation of the GO3–LiAB complex, we also achieve uniform metal doping on the GO surface for physisorption of H2. The GO–Li(n) products can store up to 5 wt% H2 and the GO–(Li3N3B3)(n) can still store 5 wt% H2. The dehydrogenation product of the GO3–LiAB complex has bridged the chemical and physical hydrogen storage approaches to move towards on-board hydrogen storage applications, which expands the scope for designing more efficient hydrogen storage materials.
Co-reporter:Jinxuan Liu, Wencai Zhou, Jianxi Liu, Yamato Fujimori, Tomohiro Higashino, Hiroshi Imahori, Xue Jiang, Jijun Zhao, Tsuneaki Sakurai, Yusuke Hattori, Wakana Matsuda, Shu Seki, Suresh Kumar Garlapati, Subho Dasgupta, Engelbert Redel, Licheng Sun and Christof Wöll
Journal of Materials Chemistry A 2016 - vol. 4(Issue 33) pp:NaN12747-12747
Publication Date(Web):2016/07/13
DOI:10.1039/C6TA04898F
We demonstrate the fabrication of a new class of epitaxial porphyrin metal–organic framework thin films whose photophysical properties can be tuned by the introduction of electron-donating diphenylamine (DPA) groups into the porphyrin skeleton. The attachment of DPA groups results in strongly improved absorption characteristics, yielding the highest photocarrier generation efficiency reported so far. DFT calculations identify a shift of the charge localization pattern in the VBM (lowest unoccupied molecular orbital), confirming that the introduction of the DPA groups is the main reason for the shift of the optical absorption spectrum and the improved photocurrent generation.
Co-reporter:Xue Jiang, Peng Wang and Jijun Zhao
Journal of Materials Chemistry A 2015 - vol. 3(Issue 15) pp:NaN7758-7758
Publication Date(Web):2015/01/14
DOI:10.1039/C4TA03438D
Since the graphene boom, great efforts have been devoted to two-dimensional (2D) monolayer materials with exciting possibilities for applications. Most known 2D materials are inorganic. Using the covalent triazine framework (CTF) as a representative, we explored 2D organic semiconductors using first-principles calculations. From a systematic study of the electronic band structures, work functions, CBM/VBM positions, and optical absorption spectra, we identified the CTF as a new class of 2D visible-light-driven organocatalyst for water splitting. Controllable construction of such CTFs from suitable organic subunits paves the way to correlate band alignment and chemical composites. In addition, multilayer CTFs have enhanced visible-light absorption compared to monolayer CTFs due to interlayer coupling. Our theoretical prediction not only has fulfilled the search for organic counterparts of inorganic photocatalysts for water splitting, but also would motivate scientists to further search for novel 2D organic materials with other technological applications.
Co-reporter:Kai Cheng, Yu Guo, Nannan Han, Yan Su, Junfeng Zhang and Jijun Zhao
Journal of Materials Chemistry A 2017 - vol. 5(Issue 15) pp:NaN3795-3795
Publication Date(Web):2017/03/13
DOI:10.1039/C7TC00595D
Lateral semiconductor/semiconductor heterostructures made up of two-dimensional (2D) monolayer or few-layer materials provide new opportunities for 2D devices. Herein, we propose four lateral heterostructures constructed by phosphorene-like monolayer group-IV monochalcogenides, including GeS/GeSe, SnS/GeSe, SnSe/GeS and GeS/SnS. Using first-principles calculations, we investigated the energetics and electronic properties of these lateral heterostructures. The band structures and formation energies from supercell calculations demonstrate that these heterostructures retain semiconducting behavior and can be easily synthesized in the laboratory. The band offsets of monolayer, bilayer and trilayer heterojunctions at the Anderson limit are calculated from the valence/conduction band edges with respect to the vacuum energy level for each individual component. Among them, some heterostructures belong to type II band alignment and are promising for a high-efficiency solar cell.
Co-reporter:Fen Li, Yan Su and Jijun Zhao
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 36) pp:NaN25248-25248
Publication Date(Web):2016/08/12
DOI:10.1039/C6CP04071C
The advance of lithium sulfur batteries is now greatly restricted by the fast capacity fading induced by shuttle effect. Using first-principles calculations, various vacancies, N doping, and B,N co-doping in graphene sheets have been systematically explored for lithium polysufides entrapped in Li–S batteries. The Li⋯S, Li⋯C, Li⋯N and S⋯B bonds and Hirshfeld charges in the Li2S6 adsorbed defective graphene systems have been analyzed to understand the intrinsic mechanism of retaining lithium polysulfides in these systems. Total and local densities of states analyses elucidate the strongest adsorption sites among the N and B–N co-doped graphene systems. The overall electrochemical performance of Li–S batteries varies with the types of defects in graphene. Among the defective graphene systems, only the reconstructed pyrrole-like vacancy is effective for retaining lithium polysulfides. N doping induces a strong Li⋯N interaction in the defective graphene systems, in which the pyrrolic N rather than the pyridinic N plays a dominant role in trapping of lithium polysulfides. The shuttle effect can be further depressed via pyrrolic B,N co-doped defective graphene materials, especially the G–B–N-hex system with extremely strong adsorption of lithium polysulfides (4–5 eV), and simultaneous contribution from the strong Li⋯N and S⋯B interactions.
Co-reporter:Heng-Fu Lin, Li-Min Liu and Jijun Zhao
Journal of Materials Chemistry A 2017 - vol. 5(Issue 9) pp:NaN2300-2300
Publication Date(Web):2017/02/07
DOI:10.1039/C7TC00013H
A recently discovered layered semiconducting material, namely, black phosphorus, has a bandgap that depends on the number of layers. It is thus feasible to create lateral heterostructures using the same material with different thicknesses. The structural and electronic properties of the lateral heterostructures of bilayer and monolayer phosphorene (BP/MP) are investigated by first-principles calculations. Hydrogen passivated heterostructures have much lower formation energy when compared to unpassivated heterostructures or phosphorene grain boundaries. The electronic band structures of the heterostructures with and without hydrogen passivation are greatly dependent on interface orientation and exhibit three distinct families of characteristics: direct bandgap semiconductor, indirect bandgap semiconductor and metal. Additionally, a type-I to type-II band alignment transition takes place when the ribbon widths of the BP and MP regions decrease, enabling continuous modulation of the band offset. The magnitude of the band offset can also be effectively tuned by changing the stacking type of BP and the interface orientation. These theoretical findings would be helpful in the design and optimization of BP/MP heterostructures for electronic and optoelectronic applications.