Ab initio calculation of the structural and electronic properties of carbon and boron nitride using ultrasoft pseudopotentials

Experimental and Computational Insights into the Catalytic Mechanism of Y1-xBaxCoO3-δ Perovskite Oxides with a Controlled Crystal Structure

The crystal structure of Co-based perovskite oxides (ACoO3) can be controlled by adjusting the A-site elements. In this study, we synthesized Y1-xBaxCoO3-δ (x = 0, 0.5, and 1.0) via a coprecipitation method and investigated their CO oxidation performances. YCoO3 (x = 0; cubic perovskite oxide; Pbnm) shows a higher catalytic performance than Y0.5Ba0.5CoO2.72 (x = 0.5; A-site-ordered double perovskite oxide; P4/nmm), which exhibits high oxygen nonstoichiometric properties, and BaCoO3 (x = 1.0; hexagonal perovskite oxide; P63/mmc), which contains high-valent Co4+ species. To elucidate the reaction mechanism, we conducted isotopic experiments with CO and 18O2. The CO oxidation reaction on YCoO3 proceeds via the Langmuir-Hinshelwood mechanism, which is a surface reaction of CO and O2 gas that does not utilize lattice oxygen. Because of the significantly smaller specific surface area of YCoO3 compared with that of the reference Pt/Al2O3, the bulk features of the crystal structures affect the catalytic reaction. When density functional theory is applied, YCoO3 clearly exhibits semiconducting properties in the ground state with the diamagnetic t2g6eg0 states, which can translate to a magnetic t2g5eg1 configuration upon excitation by a relatively low energy of 0.64 eV. We propose that the unique nature of YCoO3 activates oxygen in the gas phase, thereby enabling the smooth oxidation of CO. This study demonstrates that the bulk properties originating from the crystal structure contribute to the catalytic activity and reaction mechanism.

Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces

Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.

Toward accurate modeling of structure and energetics of bulk hexagonal boron nitride

Materials that exhibit both strong covalent and weak van der Waals interactions pose a considerable challenge to many computational methods, such as DFT. This makes assessing the accuracy of calculated properties, such as exfoliation energies in layered materials like hexagonal boron nitride (h-BN) problematic, when experimental data are not available. In this paper, we investigate the accuracy of equilibrium lattice constants and exfoliation energy calculation for various DFT-based computational approaches in bulk h-BN. We contrast these results with available experiments and reference fixed-node diffusion quantum Monte Carlo (QMC) results. From our reference QMC calculation, we obtained an exfoliation energy of - 33 ± 2 meV/atom (-0.38 ± 0.02 J/m2 ).

Machine Learning-Assisted Computational Screening of Adhesive Molecules Derived from Dihydroxyphenyl Alanine

Marine mussels adhere to virtually any surface via 3,4-dihydroxyphenyl-L-alanines (L-DOPA), an amino acid largely contained in their foot proteins. The biofriendly, water-repellent, and strong adhesion of L-DOPA are unparalleled by any synthetic adhesive. Inspired by this, we computationally designed diverse derivatives of DOPA and studied their potential as adhesives or coating materials. We used first-principles calculations to investigate the adsorption of the DOPA derivatives on graphite. The presence of an electron-withdrawing group, such as nitrogen dioxide, strengthens the adsorption by increasing the π-π interaction between DOPA and graphite. To quantify the distribution of electron charge and to gain insights into the charge distribution at interfaces, we performed Bader charge analysis and examined charge density difference plots. We developed a quantitative structure-property relationship (QSPR) model using an artificial neural network (ANN) to predict the adsorption energy. Using the three-dimensional and quantum mechanical electrostatic potential of a molecule as a descriptor, the present quantum NN model shows promising performance as a predictive QSPR model.

Cathodoluminescence spectroscopy of monolayer hexagonal boron nitride

Cathodoluminescence (CL) spectroscopy is a suitable technique for studying the luminescent properties of optoelectronic materials because CL has no limitation on the excitable bandgap energy and eliminates ambiguous signals due to simple light scattering and resonant Raman scattering potentially involved in the photoluminescence spectra. However, direct CL measurements of atomically thin two-dimensional materials have been difficult due to the small excitation volume that interacts with high-energy electron beams. Herein, distinct CL signals from a monolayer hexagonal BN (hBN), namely mBN, epitaxial film grown on a graphite substrate are shown by using a CL system capable of large-area and surface-sensitive excitation. Spatially resolved CL spectra at 13 K exhibited a predominant 5.5-eV emission band, which has been ascribed to originate from multilayered aggregates of hBN, markedly at thicker areas formed on the step edges of the substrate. Conversely, a faint peak at 6.04 ± 0.01 eV was routinely observed from atomically flat areas, which is assigned as being due to the recombination of phonon-assisted direct excitons of mBN. The CL results support the transition from indirect bandgap in bulk hBN to direct bandgap in mBN. The results also encourage one to elucidate emission properties of other low-dimensional materials by using the present CL configuration.

Adsorption behavior of O2 on U-Nb surface: a first-principles study

CONTEXT

Oxygen is a chemically active gas. Metal uranium will be rapidly oxidized when exposed to oxygen directly, forming a complex and uncompacted oxide layer on the surface. The U-O2 reaction system is a very complex metal oxidation system. The surface oxidation reaction of uranium-niobium alloy is more complex, and it is also the main corrosion form of uranium-niobium alloy. Exploring the microscopic mechanism of surface oxidation corrosion of uranium-niobium alloy is of great significance for understanding the surface corrosion phenomenon in practical applications. The adsorption of O2 on U-Nb (5at.%) surface was investigated by first-principles calculation using a periodic slab model within the density functional theory (DFT). The effect of different levels of Nb doping on O adsorption was investigated, and the stability of O-U and O-Nb adsorption was studied in conjunction with the Bader charge distribution, differential charge density map, electron density of states, and surface work function. The results showed that physical adsorption happens when O2 is vertically adsorbed at the top site of Nb atom, whereas dissociative adsorption happened in all other configurations considered here. The interaction between adsorbed O2 and U-Nb surface was essentially the hybridization between the O/2p orbital electrons and the U/6d, U/5f, Nb/4p, Nb/5s, and Nb/4d orbital electrons to form a relatively stable ionic bond. When Nb was doped into the second layer of the surface, the most stable O2 adsorption configuration was top-site horizontal adsorption. The adsorption energy was - 22.38 eV, which was more negative than the adsorption energy of - 21.38 eV for the first doped layer horizontal adsorption at hollow-site. The interaction was essentially the hybridization between O/2p orbital electrons and U/6d and U/5f electrons and between O/2s orbital electrons and U/6p orbital electrons to form a relatively stable ionic bond.

METHOD

In this paper, the adsorption and dissociation of oxygen molecules on the surface of U-Nb system were systematically studied by first-principles calculations, in order to explore the surface oxidation corrosion mechanism of uranium-niobium alloy from the micro level. The VASP program uses PAW pseudopotential, which belongs to the density functional method, and uses DSP.

Hexagonal Boron Nitride for Photonic Device Applications: A Review

Hexagonal boron nitride (hBN) has emerged as a key two-dimensional material. Its importance is linked to that of graphene because it provides an ideal substrate for graphene with minimal lattice mismatch and maintains its high carrier mobility. Moreover, hBN has unique properties in the deep ultraviolet (DUV) and infrared (IR) wavelength bands owing to its indirect bandgap structure and hyperbolic phonon polaritons (HPPs). This review examines the physical properties and applications of hBN-based photonic devices that operate in these bands. A brief background on BN is provided, and the theoretical background of the intrinsic nature of the indirect bandgap structure and HPPs is discussed. Subsequently, the development of DUV-based light-emitting diodes and photodetectors based on hBN's bandgap in the DUV wavelength band is reviewed. Thereafter, IR absorbers/emitters, hyperlenses, and surface-enhanced IR absorption microscopy applications using HPPs in the IR wavelength band are examined. Finally, future challenges related to hBN fabrication using chemical vapor deposition and techniques for transferring hBN to a substrate are discussed. Emerging techniques to control HPPs are also examined. This review is intended to assist researchers in both industry and academia in the design and development of unique hBN-based photonic devices operating in the DUV and IR wavelength regions.

Efficient fabrication of single-wall carbon nanotube nanoreactors by defect-induced cutting

Single-wall carbon nanotubes (SWCNTs) with ultra-thin channels are considered promising nanoreactors for confined catalysis, chemical reactions, and drug delivery. The fabrication of SWCNT nanoreactors by cutting usually suffers from low efficiency and poor controllability. Here we develop a defect-induced gas etching method to efficiently cut SWCNTs and to obtain nanoreactors with ultrasmall confined space. H₂ plasma treatment was performed to generate defects in the walls of SWCNTs, then H₂O vapor was used as a "knife" to cut SWCNTs at the defect sites, and short cut-SWCNTs with an average length of 175 nm were controllably obtained with a high yield of 75% under optimized conditions. WO₃@SWCNT derivatives with different morphologies were synthesized using short cut-SWCNTs as nanoreactors. The radiation resistance of WO₃@SWCNT hybrids improved obviously, thus providing a platform for the synthesis of novel SWCNT-based derivatives with fascinating properties.

Predicting the structures and vibrational spectra of molecular crystals containing large molecules with the generalized energy-based fragmentation approach

The generalized energy-based fragmentation (GEBF) approach under periodic boundary conditions (PBCs) has been developed to facilitate calculations of molecular crystals containing large molecules. The PBC-GEBF approach can help predict structures and properties of molecular crystals at different theory levels by performing molecular quantum chemistry calculations on a series of non-periodic subsystems constructed from the studied systems. A more rigorous formula of the forces on translational vectors of molecular crystals was proposed and implemented, enabling more reliable predictions of crystal structures. Our benchmark results on several typical molecular crystals show that the PBC-GEBF approach could reproduce the forces on atoms and the translational vectors and the optimized crystal structures from the corresponding conventional periodic methods. The improved PBC-GEBF approach is then applied to predict the crystal structures and vibrational spectra of two molecular crystals containing large molecules. The PBC-GEBF approach can provide a satisfactory description on the crystal structure of a molecular crystal containing 312 atoms in a unit cell at density-fitting second-order Møller-Plesset perturbation theory and density functional theory (DFT) levels and the infrared vibrational spectra of another molecular crystal containing 864 atoms in a unit cell at the DFT level. The PBC-GEBF approach is expected to be a promising theoretical tool for electronic structure calculations on molecular crystals containing large molecules.

Room-Temperature Deep-UV Photoluminescence from Low-Dimensional Hexagonal Boron Nitride Prepared Using a Facile Synthesis

Identification and evaluation of defect levels in low-dimensional materials is an important aspect in quantum science. In this article, we report a facile synthesis method of low-dimensional hexagonal boron nitride (h-BN) and study light emission characteristics due to the defects. The thermal annealing procedure is optimized to obtain clean multilayered h-BN as revealed by transmission electron microscopy. UV-vis spectroscopy shows the optical energy gap of 5.28 eV, which is comparable to the reported energy gap for exfoliated, clean h-BN samples. X-ray photoelectron spectroscopy reveals the location of the valence band edge at 2 eV. The optimized synthesis route of h-BN generates two kinds of defects, which are characterized using room-temperature photoluminescence (PL) measurements. The defects emit light at 4.18 eV [deep-UV (DUV)] and 3.44 eV (UV) photons. The intensity of PL has an oscillatory dependence on the excitation energy for the defect emitting DUV light. A series of spectral lines are observed with the energy ranging between 2.56 and 3.44 eV. The average peak-to-peak energy separation is about 125 meV. The locations of the spectral lines can be modeled using Franck-Condon-type transition and associated with displaced harmonic oscillator approximation. Our facile route gives an easier approach to prepare clean h-BN, which is essential for classical two-dimensional material-based electronics and single-photon-based quantum devices.

Engineering van der Waals Materials for Advanced Metaphotonics

The outstanding chemical and physical properties of 2D materials, together with their atomically thin nature, make them ideal candidates for metaphotonic device integration and construction, which requires deep subwavelength light-matter interaction to achieve optical functionalities beyond conventional optical phenomena observed in naturally available materials. In addition to their intrinsic properties, the possibility to further manipulate the properties of 2D materials via chemical or physical engineering dramatically enhances their capability, evoking new science on light-matter interaction, leading to leaped performance of existing functional devices and giving birth to new metaphotonic devices that were unattainable previously. Comprehensive understanding of the intrinsic properties of 2D materials, approaches and capabilities for chemical and physical engineering methods, the resulting property modifications and novel functionalities, and applications of metaphotonic devices are provided in this review. Through reviewing the detailed progress in each aspect and the state-of-the-art achievement, insightful analyses of the outstanding challenges and future directions are elucidated in this cross-disciplinary comprehensive review with the aim to provide an overall development picture in the field of 2D material metaphotonics and promote rapid progress in this fast emerging and prosperous field.

Polytypes of sp2-Bonded Boron Nitride

The sp2-bonded layered compound boron nitride (BN) exists in more than a handful of different polytypes (i.e., different layer stacking sequences) with similar formation energies, which makes obtaining a pure monotype of single crystals extremely tricky. The co-existence of polytypes in a similar crystal leads to the formation of many interfaces and structural defects having a deleterious influence on the internal quantum efficiency of the light emission and on charge carrier mobility. However, despite this, lasing operation was reported at 215 nm, which has shifted interest in sp2- bonded BN from basic science laboratories to optoelectronic and electrical device applications. Here, we describe some of the known physical properties of a variety of BN polytypes and their performances for deep ultraviolet emission in the specific case of second harmonic generation of light.

Structures and Electromagnetic Properties of Boron Nitride Nanoribbons Doped with Transition Metals

Inspired by the recent discovery of the Ti-doped BN nanocages, here we report the design of novel BN nanoribbons (BNNRs) doped with fourth-row transition metals (Sc-Cu) and the prediction of their structural and electromagnetic properties. First-principles calculations and ab initio molecular dynamics simulations show that Ti-doped BNNR possesses both thermodynamic and kinetic stability at high temperatures for synthesis of BN materials. Metal doping may make the nonmagnetic pristine BNNR ferromagnetic or antiferromagnetic, depending on the metal. The doping with all considered metals reduces substantially the band gap of pristine BNNR. For example, Sc-doped BNNR is ferromagnetic with an indirect band gap of 1.18 eV, while V-doped nanoribbon is antiferromagnetic with a direct gap of 2.50 eV. Remarkably, the carrier mobility in both materials is significantly enhanced compared to the pristine BNNR. Our findings suggest that doping with different metals may endow BNNRs with versatile electronic and magnetic properties.

Structure, Properties and Applications of Two-Dimensional Hexagonal Boron Nitride

Hexagonal boron nitride (h-BN) has emerged as a strong candidate for two-dimensional (2D) material owing to its exciting optoelectrical properties combined with mechanical robustness, thermal stability, and chemical inertness. Super-thin h-BN layers have gained significant attention from the scientific community for many applications, including nanoelectronics, photonics, biomedical, anti-corrosion, and catalysis, among others. This review provides a systematic elaboration of the structural, electrical, mechanical, optical, and thermal properties of h-BN followed by a comprehensive account of state-of-the-art synthesis strategies for 2D h-BN, including chemical exfoliation, chemical, and physical vapor deposition, and other methods that have been successfully developed in recent years. It further elaborates a wide variety of processing routes developed for doping, substitution, functionalization, and combination with other materials to form heterostructures. Based on the extraordinary properties and thermal-mechanical-chemical stability of 2D h-BN, various potential applications of these structures are described.

High-efficient ab initio Bayesian active learning method and applications in prediction of two-dimensional functional materials

Beyond the conventional trial-and-error method, machine learning offers a great opportunity to accelerate the discovery of functional materials, but still often suffers from difficulties such as limited materials data and the unbalanced distribution of target properties. Here, we propose the ab initio Bayesian active learning method that combines active learning and high-throughput ab initio calculations to accelerate the prediction of desired functional materials with ultrahigh efficiency and accuracy. We apply it as an instance to a large family (3119) of two-dimensional hexagonal binary compounds with unbalanced materials properties, and accurately screen out the materials with maximal electric polarization and proper photovoltaic band gaps, respectively, whereas the computational costs are significantly reduced by only calculating a few tenths of the possible candidates in comparison with a random search. This approach shows the enormous advantages for the cases with unbalanced distribution of target properties. It can be readily applied to seek a broad range of advanced materials.

Lattice dynamics in the conformational environment of multilayered hexagonal boron nitride (h-BN) results in peculiar infrared optical responses

Stacking mismatches in hexagonal boron nitride (h-BN) nanostructures affect their photonic, mechanical, and thermal properties. To access information about the stacked configuration of layered ensembles, highly sophisticated techniques like X-ray photoemission spectroscopy or electron microscopy are necessary. Here, instead, by taking advantage of the geometrical and chemical nature of h-BN, we show how simple structural models, based on shortened interplanar distances, can produce effective charge densities. Accounting these in the non-analytical part of the lattice dynamical description makes it possible to access information about the composition of differently stacked variants in experimental samples characterized by infrared spectroscopy. The results are obtained by density functional theory and confirmed by various functionals and pseudopotential approximations. Even though the method is shown using h-BN, the conclusions are more general and show how effective dielectric models can be considered as valuable theoretical pathways for the vibrational structure of any layered material.

A possible unaccounted source of atmospheric sulfate formation: amine-promoted hydrolysis and non-radical oxidation of sulfur dioxide

Numerous field and laboratory studies have shown that amines, especially dimethylamine (DMA), are crucial to atmospheric particulate nucleation. However, the molecular mechanism by which amines lead to atmospheric particulate formation is still not fully understood. Herein, we show that DMA molecules can also promote the conversion of atmospheric SO2 to sulfate. Based on ab initio simulations, we find that in the presence of DMA, the originally endothermic and kinetically unfavourable hydrolysis reaction between gaseous SO2 and water vapour can become both exothermic and kinetically favourable. The resulting product, bisulfite NH2(CH3)2 +·HSO3 -, can be readily oxidized by ozone under ambient conditions. Kinetic analysis suggests that the hydrolysis rate of SO2 and DMA with water vapour becomes highly competitive with and comparable to the rate of the reaction between SO2 and OH·, especially under the conditions of heavily polluted air and high humidity. We also find that the oxidants NO2 and N2O5 (whose role in sulfate formation is still under debate) appear to play a much less significant role than ozone in the aqueous oxidation reaction of SO2. The newly identified oxidation mechanism of SO2 promoted by both DMA and O3 provides another important new source of sulfate formation in the atmosphere.

Prediction of a novel robust superconducting state in TaS2 under high pressure

A novel superconducting I4/mmm phase has been predicted in TaS₂ under high pressure, illustrating an unusual superconductor–metal–superconductor transition.

Multiscale molecular thermodynamics of graphene-oxide liquid-phase exfoliation

Liquid-phase exfoliation is one of the most feasible methods for mass-production of two-dimensional (2D) nanomaterials such as graphene, graphene-oxide (GO), etc. Assessing requirements for successful exfoliation necessitates molecular-level thermodynamic analysis that can provide quantitative measures such as free energy changes. Here we explain this methodology and apply it to the production of GO that is used as a precursor for graphene synthesis and as an ultrathin substrate for many applications. Three different routes to GO exfoliation are studied, namely parallel and perpendicular to the GO surface as well as exfoliation via edge bending, using multi-scale combination of density functional, force field, and continuum approaches. Detailed analysis of free energy variations reveals relative feasibility of different exfoliation mechanisms and their dependence on system size and surface coverage. The methodology is general and can be applied to liquid-phase exfoliation of other 2D nanomaterials.

Novel carbon polymorphs with cumulative double bonds in three-dimensional sp-sp2 hybrid framework

A conspicuous amount of theoretical study has been published on the properties of carbon allotropes with alternate single and triple bonds, (-C[triple bond, length as m-dash]C-)n. However, theoretical characterizations of carbon allotropes with cumulative double bonds ([double bond, length as m-dash]C[double bond, length as m-dash]C[double bond, length as m-dash])n is almost non-existent in literature. Based upon first-principles calculations, two new three-dimensional (3D) microporous carbon allotropes consisting of whorl chains connected by cumulative double bonds in a sp-sp2 hybrid framework have been proposed in this study. One of these structures, namely, Trig-C9 was obtained by an evolutionary particle swarm structural search, while the other structure, denoted as Trig-C15, was obtained by inserting double bonds into Trig-C9. Both the 3D sp-sp2 hybridized carbons have a trigonal structure with 9 and 15 atoms in the hexagonal primitive cells. The calculated results demonstrate that these polymorphs are thermodynamically, mechanically, and dynamically feasible. Trig-C9 and Trig-C15 are indirect semiconductors with band gaps of 2.70 eV and 1.25 eV, respectively. Their unique frameworks render them mechanical ductility and significant elastic anisotropy. These results open up new horizons for the exploration of new carbon phases with unique structural, mechanical, and electronic properties.

Role of zero-point effects in stabilizing the ground state structure of bulk Fe2P

Structural stability of Fe₂P is investigated in detail using first-principles calculations based on density functional theory. While the orthorhombic C23 phase is found to be energetically more stable, the experiments suggest it to be hexagonal C22 phase. In the present study, we show that in order to obtain the correct ground state structure of Fe₂P from the first-principles based methods it is utmost necessary to consider the zero-point effects such as zero-point vibrations and spin fluctuations. This study demonstrates an exceptional case where a bulk material is stabilized by quantum effects, which are usually important in low-dimensional materials. Our results also indicate the possibility of magnetic field induced structural quantum phase transition in Fe₂P, which should form the basis for further theoretical and experimental efforts.

A theoretical study of formaldehyde adsorption and decomposition on a WC (0001) surface

A lot of research attention has been paid to designing and exploring efficient adsorbents for HCHO adsorption and decomposition. Herein, we have demonstrated a highly active material, WC, for HCHO adsorption through first-principles calculations. Due to the exposed three-coordinated W atoms (W_3c) of the WC (0001) surface, HCHO molecules can be settled on the WC (0001) surface through newly formed O_F–W_3c and/or C_F–W_3c bonds, forming different adsorption configurations. When settled on the WC (0001) surface, the molecular configuration of the HCHO molecule and the corresponding C_F–H_F and C_F–O_F bond lengths would be greatly changed. Due to the enlarged C_F–H_F and C_F–O_F bond lengths, the adsorbed HCHO molecules tend to dissociate through two possible pathways; these are the two-step C_F–H_F bond scission or the one-step C_F–O_F bond scission. The former results in two H adatoms and a CO molecule chemisorbed to the surface and the latter produces an O adatom and a CH₂ group on the surface. Further Cl-NEB calculations demonstrate that the pre-adsorbed O atom has little influence on the first C_F–H_F bond scission and the C_F–O_F bond scission, while promoting the second C_F–H_F bond scission. Considering the dissociative products, H and CH₂ have the potential to couple into a CH₃ group (or even a CH₄ molecule) and two CH₂ groups may couple into a C₂H₄ molecule. In the end, we propose that OH ions may couple with the dissociative products of HCHO, so an alkali solution could be used to post-process the WC (0001) surface to restore its surface active sites. These results demonstrated the potential of WC in HCHO sensing and abatement.

WC is a material capable of HCHO adsorption and dissociation, indicating its potential application in HCHO sensing and elimination.

Resonant X-ray Emission of Hexagonal Boron Nitride

The electronic structure of hexagonal boron nitride (h-BN) is explored using measurements of x-ray absorption and resonant inelastic x-ray scattering (RIXS) at the nitrogen K edge (1s) in tandem with calculations using many-body perturbation theory within the GW and Bethe-Salpeter equation (BSE) approximations. Our calculations include the effects of lattice disorder from phonons activated thermally and from zero point energy. They highlight the influence of disorder on near-edge x-ray spectra.

A Density Functional Theory Study on the Effect of Lattice Impurities on the Electronic Structures and Reactivity of Fluorite

Fluorite (CaF2), a halogen elemental mineral, always co-exists with other minerals. The Ca element in fluorite is often replaced by rare earth elements (REEs), such as cerium (Ce) and yttrium (Y). In this work, the electronic structures of fluorite crystals containing REE (Ce, Th, U, and Y) impurities were studied by density functional theory (DFT). The calculated results showed that the presence of impurities increased the lattice parameter of fluorite. The impurities caused the Fermi level to shift towards the high energy direction, making the fluorite accept electrons more easily. The impurities except Y led to the occurrence of an impurities state in the valence band. The Mullinken population values of F–REE bonds were larger than that of F–Ca and F–F bonds, and F–Y bonds had the largest population value. Analysis of the frontier molecular orbital showed that the impurities contributed greatly to the lowest unoccupied molecular orbital (LUMO). The interaction between oleic acid and impurities-bearing fluorite were discussed. The results suggested that the incorporation of impurities would enhanced the reactivity of fluorite with oleic acid.

Novel hybrid C/BN two-dimensional heterostructures

Here we present an investigation of new quasi-two-dimensional heterostructures based on the alternation of bounded carbon and boron nitride layers (C/BN). We carried out a theoretical study of the atomic structure, stability and electronic properties of the proposed heterostructures. Such ultrathin quasi-two-dimensional C/BN films can be synthesized by means of chemically induced phase transition by connection of the layers of multilayered h-BN/graphene van der Waals heterostructures, which is indicated by the negative phase transition pressure in the calculated phase diagrams (P, T) of the films. It was shown that the band gap value of the C/BN films spans the infrared and visible spectrum. We hope that the proposed films and fabrication method can be considered as a possible route to obtain nanostructures with a controllable band gap in wide energy range. This makes these materials potentially suitable for a variety of applications, including photovoltaics, photoelectronics and more.

Computational Search for Novel Hard Chromium-Based Materials

Nitrides, carbides, and borides of transition metals are an attractive class of hard materials. Our recent preliminary explorations of the binary chemical compounds indicated that chromium-based materials are among the hardest transition metal compounds. Motivated by this, here we explore in detail the binary Cr-B, Cr-C, and Cr-N systems using global optimization techniques. Calculated enthalpy of formation and hardness of predicted materials were used for Pareto optimization to define the hardest materials with the lowest energy. Our calculations recover all numerous known stable compounds (except Cr₂₃C₆ with its large unit cell) and discover a novel stable phase Pmn2₁-Cr₂C. We resolve the structure of Cr₂N and find it to be of anti-CaCl₂ type (space group Pnnm). Many of these phases possess remarkable hardness, but only CrB₄ is superhard (Vickers hardness 48 GPa). Among chromium compounds, borides generally possess the highest hardnesses and greatest stability. Under pressure, we predict stabilization of a layered TMDC-like phase of Cr₂N, a WC-type phase of CrN, and a new compound CrN₄. Nitrogen-rich chromium nitride CrN₄ is a high-energy-density material featuring polymeric nitrogen chains. In the presence of metal atoms (e.g., Cr), polymerization of nitrogen takes place at much lower pressures; CrN₄ becomes stable at ∼15 GPa (cf. 110 GPa for synthesis of pure polymeric nitrogen).

Ti3BN monolayer: the MXene-like material predicted by first-principles calculations

An MXene-like Ti₃BN monolayer whose electronic properties could be modulated has been predicted following the strategy of “atomic transmutation”.

C 20 - T carbon: a novel superhard sp (3) carbon allotrope with large cavities

Through first-principles calculations, we predict a new superhard carbon allotrope named C 20  -  T, which possesses a cubic T symmetry with space group No.198(P213). This new carbon allotrope has an all-sp (3) hybridized bonding network with 20 atoms in its primitive unit cell. The dynamic, mechanical, and thermal stabilities of this new carbon phase at zero pressure are confirmed by using a variety of state-of-the-art theoretical calculations. Interestingly, despite the fact that C 20  -  T carbon has a porous structure with large cavities, our calculations identify its superhard properties with the Vickers hardness of 72.76 Gpa. The ideal tensile and shear strength of C 20  -  T carbon are calculated to be 71.1 and 55.2 GPa respectively, comparable to that of c-BN. Electronic band calculations reveal that this new carbon allotrope is a transparent insulator with an indirect band gap of 5.44 eV. These results broaden our understanding of superhard carbon allotropes.

Adsorption and dynamics of group IV, V atoms and molecular oxygen on semiconductor group IV (0 0 1) surfaces

In this review we address (1) the co-adsorption of group V (As, Sb, Bi) atoms and molecular oxygen on the Si(0 0 1) surface and (2) the adsorption and dynamics of Sb, Bi, Si and Ge ad-dimers on the Si(0 0 1) and Ge(0 0 1) surfaces. The adsorption and diffusion processes of group IV and V atoms on the (0 0 1) surfaces of group IV semiconductor surfaces have been studied using multi-configuration self-consistent field methods and density functional theory calculations. Results obtained by various types of first-principle total energy calculations are mutually compared and discussed. Our results demonstrate the capability of these quantum chemistry methods to provide relevant and reliable information on the interaction between adsorbate and semiconductor surfaces.

H18 Carbon: A New Metallic Phase with sp2-sp3 Hybridized Bonding Network

Design and synthesis of three-dimensional metallic carbons are currently one of the hot issues in contemporary condensed matter physics because of their fascinating properties. Here, based on first-principles calculations, we discover a novel stable metallic carbon allotrope (termed H₁₈ carbon) in () symmetry with a mixed sp²-sp³ hybridized bonding network. The dynamical stability of H₁₈ carbon is verified by phonon mode analysis and molecular dynamics simulations, and its mechanical stability is analyzed by elastic constants, bulk modulus, and shear modulus. By simulating the x-ray diffraction patterns, we propose that H₁₈ carbon would be one of the unidentified carbon phases observed in recent detonation experiments. The analysis of the band structure and density of states reveal that this new carbon phase has a metallic feature mainly due to the C atoms with sp² hybridization. This novel 3D metallic carbon phase is anticipated to be useful for practical applications such as electronic and mechanical devices.

Temperature dependence of the electronic structure of semiconductors and insulators

The renormalization of electronic eigenenergies due to electron-phonon coupling (temperature dependence and zero-point motion effect) is sizable in many materials with light atoms. This effect, often neglected in ab initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the recent progresses in this field and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a slowly converging q-point integration. For non-zero Born effective charges, we show that a divergence appears in the electron-phonon matrix elements at q → Γ, leading to a divergence of the adiabatic renormalization at band extrema. This problem is exacerbated by the slow convergence of Born effective charges with electronic wavevector sampling, which leaves residual Born effective charges in ab initio calculations on materials that are physically devoid of such charges. Here, we propose a solution that improves this convergence. However, for materials where Born effective charges are physically non-zero, the divergence of the renormalization indicates a breakdown of the adiabatic harmonic approximation, which we assess here by switching to the non-adiabatic harmonic approximation. Also, we study the convergence behavior of the renormalization and develop reliable extrapolation schemes to obtain the converged results. Finally, the adiabatic and non-adiabatic theories, with corrections for the slow Born effective charge convergence problem (and the associated divergence) are applied to the study of five semiconductors and insulators: α-AlN, β-AlN, BN, diamond, and silicon. For these five materials, we present the zero-point renormalization, temperature dependence, phonon-induced lifetime broadening, and the renormalized electronic band structure.

Temperature dependent structural properties and bending rigidity of pristine and defective hexagonal boron nitride

Structural and thermodynamical properties of monolayer pristine and defective boron nitride sheets (h-BN) have been investigated in a wide temperature range by carrying out atomistic simulations using a tuned Tersoff-type inter-atomic empirical potential. The temperature dependence of lattice parameter, radial distribution function, specific heat at constant volume, linear thermal expansion coefficient and the height correlation function of the thermally excited ripples on pristine as well as defective h-BN sheet have been investigated. Specific heat shows considerable increase beyond the Dulong-Petit limit at high temperatures, which is interpreted as a signature of strong anharmonicity present in h-BN. Analysis of the height fluctuations, ⟨h2⟩, shows that the bending rigidity and variance of height fluctuations are strongly temperature dependent and this is explained using the continuum theory of membranes. A detailed study of the height-height correlation function shows deviation from the prediction of harmonic theory of membranes as a consequence of the strong anharmonicity in h-BN. It is also seen that the variance of the height fluctuations increases with defect concentration.

van der Waals forces in density functional theory: a review of the vdW-DF method

A density functional theory (DFT) that accounts for van der Waals (vdW) interactions in condensed matter, materials physics, chemistry, and biology is reviewed. The insights that led to the construction of the Rutgers-Chalmers van der Waals density functional (vdW-DF) are presented with the aim of giving a historical perspective, while also emphasizing more recent efforts which have sought to improve its accuracy. In addition to technical details, we discuss a range of recent applications that illustrate the necessity of including dispersion interactions in DFT. This review highlights the value of the vdW-DF method as a general-purpose method, not only for dispersion bound systems, but also in densely packed systems where these types of interactions are traditionally thought to be negligible.

Direction-controlled chemical doping for reversible G-phonon mixing in ABC trilayer graphene

Not only the apparent atomic arrangement but the charge distribution also defines the crystalline symmetry that dictates the electronic and vibrational structures. In this work, we report reversible and direction-controlled chemical doping that modifies the inversion symmetry of AB-bilayer and ABC-trilayer graphene. For the "top-down" and "bottom-up" hole injection into graphene sheets, we employed molecular adsorption of electronegative I2 and annealing-induced interfacial hole doping, respectively. The chemical breakdown of the inversion symmetry led to the mixing of the G phonons, Raman active Eg and Raman-inactive Eu modes, which was manifested as the two split G peaks, G(-) and G(+). The broken inversion symmetry could be recovered by removing the hole dopants by simple rinsing or interfacial molecular replacement. Alternatively, the symmetry could be regained by double-side charge injection, which eliminated G(-) and formed an additional peak, G(o), originating from the barely doped interior layer. Chemical modification of crystalline symmetry as demonstrated in the current study can be applied to other low dimensional crystals in tuning their various material properties.

Pressure evolution of the potential barriers for transformations of layered BN to dense structures

The energy barrier and stacking way from layered BN to dense phase under pressure.

Radiation-Induced Nucleation of Diamond from Amorphous Carbon: Effect of Hydrogen

Electron irradiation of anthracite functionalized by dodecyl groups leads to recrystallization of the carbon network into diamonds. The diamonds range in size from ∼2 to ∼10 nm and exhibit {111} spacing of 2.1 Å. A bulk process consistent with bias-enhanced nucleation is proposed in which the dodecyl group provides hydrogen during electron irradiation. Recrystallization into diamond occurs in the hydrogenated graphitic subsurface layers. Unfunctionalized anthracite could not be converted into diamond during electron irradiation. The dependence of the phase transition pressure on cluster size was estimated, and it was found that diamond particles with a radius up to 20 nm could be formed.

Phase diagram of quasi-two-dimensional carbon, from graphene to diamond

We explore how a few-layer graphene can undergo phase transformation into thin diamond film under reduced or no pressure, if the process is facilitated by hydrogenation of the surfaces. Such a "chemically induced phase transition" is inherently nanoscale phenomenon, when the surface conditions directly affect thermodynamics, and the transition pressure depends greatly on film thickness. For the first time we obtain, by ab initio computations of the Gibbs free energy, a phase diagram (P, T, h) of quasi-two-dimensional carbon-diamond film versus multilayered graphene. It describes accurately the role of film thickness h and shows the feasibility of creating novel quasi-two-dimensional materials. Further, the role of finite diameter of graphene flakes and possible formation of the diamond films with the (110) surface are described as well.