Search

Optimizing density-functional theory (DFT) and G0W0 calculations present coupled problems as orbitals from DFT are needed as G0W0 starting points. Applied to 341 two-dimensional (2D) materials, we demonstrate that CAM-B3LYP provides minimal changes in bandgap (e.g., mean absolute deviation of 0.23 eV) when used to start G0W0 calculations, compared to traditional functionals such as PBE, PBE0, and HSE06 (1.07 eV, 1.48 eV, and 1.51 eV, respectively). CAM-B3LYP also delivers the smallest changes in orbital representation. These and other results indicate the suitability of CAM-B3LYP as a density-functional approach for modelling 2D materials, as well as for use in optimizing G0W0 calculations. Our findings parallel well established features of applications to molecules, as well as for spectroscopic applications involving 3D materials.

Plastic scintillators that can discriminate between gamma rays, fast neutrons, and thermal neutrons were synthesized and characterized while considering the balance between processing and performance at the kilogram scale. These trade-offs were necessitated by the inclusion of 0.1 wt. % lithium-6 to enable detection of thermal neutrons. The synthesis and processing of these plastic scintillators on the kilogram scale required consideration of many factors. First, a comonomer (methacrylic acid) was used to solubilize salts of lithium-6, which allow for a thermal-neutron capture reaction that produces scintillation light following energy transfer. Second, scintillation performance and processability were considered because the increasing content of the comonomer resulted in a sharp decrease in the light output. The use of small amounts of comonomer (less than or equal to 3 wt. %) resulted in better performance but required high processing temperatures. At large scales, these high temperatures could initiate an exothermic polymerization that results in premature curing and/or defects. The deleterious effects of the comonomer may be mitigated by using m-terphenyl as a primary dye rather than 2,5-diphenyloxazole (PPO), which has been traditionally used in organic scintillators. Finally, the curing environment was controlled to avoid defects like cracking and discoloration while maintaining solubility of dopants during curing. For scintillators that were produced from kilogram-scale batches of precursors, the effective attenuation of scintillation light was characterized.

We discuss here new, enabling technologies for future photon-based neutrino detectors. These technologies touch nearly every aspect of such detectors: new scintillating materials, new methods of loading isotopes, new photon sensors and collectors, new approaches to simulation and analysis, and new front-end electronics and DAQ ideas. Of particular interest are technologies that enable broad physics programs in hybrid Cherenkov/scintillation detectors, such as slow fluors, water-based liquid scintillator, and spectral sorting of photons. Several new large-scale detector ideas are also discussed, including hybrid detectors like Theia, ArTEMIS, and generic slow-fluor detectors, as well as the very different SLIPs and LiquidO approaches to instrumenting photon-based detectors. A program of demonstrators for future detectors, including ANNIE, Eos, and NuDOT are also discussed.

This work describes a class of liquid scintillators that contain mostly water (>50 wt. % of the entire composition) and can discriminate between interactions induced by neutrons and gamma rays. By balancing the interface interactions between the components of the formulation, these scintillators form emulsions that can be thermodynamically stable. This approach, which considers a quantity known as the hydrophilic-lipophilic difference, requires consideration of the salinity and temperature as well as characterization of the surfactants and oil phase. Emulsions comprised of water and various oils were characterized first. Then, the effect of scintillating dyes in the oil phase was considered, followed by the construction of partial phase diagrams of the emulsions. For transparent oil-in-water emulsions with a single phase, the scintillation light yield and properties of pulse-shape discrimination were measured. The best performing scintillators contained 33 wt. % of a scintillating oil phase and exhibited a light yield that was as high as 18% of the light yield of a commercially available liquid scintillator that does not contain water (EJ-309). These water-based liquid scintillators exhibited a figure of merit of neutron/gamma ray discrimination as high as 1.79 at about 1500 keVee.

Density functionals with asymptotic corrections to the long-range potential provide entry-level methods for calculations on molecules that can sustain charge transfer, but similar applications in Materials Science are rare. We describe an implementation of the CAM-B3LYP range-separated functional within the Vienna Ab-initio Simulation Package (VASP) framework, together with its analytical functional derivatives. Results obtained for eight representative materials: aluminum, diamond, graphene, silicon, NaCl, MgO, 2D h-BN and 3D h-BN, indicate that CAM-B3LYP predictions embody mean-absolute deviations (MAD) compared to HSE06 that are reduced by a factor of 6 for lattice parameters, 4 for quasiparticle band gaps, 3 for the lowest optical excitation energies, and 6 for exciton binding energies. Further, CAM-B3LYP appears competitive compared to ab initio G0W0 and Bethe-Salpeter equation (BSE) approaches. The CAM-B3LYP implementation in VASP was verified by comparison of optimized geometries and reaction energies for isolated molecules taken from the ACCDB database, evaluated in large periodic unit cells, to analogous results obtained using Gaussian basis sets. Using standard GW pseudopotentials and energy cutoffs for the plane-wave calculations and the aug-cc-pV5Z basis set for the atomic-basis ones, the MAD in energy for 1738 chemical reactions was 0.34 kcal mol-1, whilst for 480 unique bond lengths this was 0.0036 {\AA}; these values reduced to 0.28 kcal mol-1 (largest error 0.94 kcal mol-1) and 0.0009 {\AA} by increasing the plane-wave cuttoff energy to 850 eV., Comment: 2 Figures, 7 tables

Extensive photochemical and spectroscopic properties of the $V_B^-$ defect in hexagonal boron nitride are calculated, concluding that the observed photoemission associated with recently observed optically-detected magnetic resonance is most likely of (1)3E" to (1)3A2' origin. Rapid intersystem crossing from the defect's triplet to singlet manifolds explains the observed short excited-state lifetime and very low quantum yield. New experimental results reveal smaller intrinsic spectral bandwidths than previously recognized, interpreted in terms spectral narrowing and zero-phonon-line shifting induced by the Jahn-Teller effect. Different types of computational methods are applied to map out the complex triplet and singlet defect manifolds, including the doubly ionised formulation of the equation-of-motion coupled-cluster theory that is designed to deal with the open-shell nature of defect states, and mixed quantum-mechanics/molecular-mechanics schemes enabling 5763-atom simulations. Two other energetically feasible spectral assignments from amongst the singlet and triplet manifolds are considered, but ruled out based on inappropriate photochemical properties.

The $V_N N_B$ defect in hexagonal boron nitride (h-BN), comprising a nitrogen vacancy adjacent to a nitrogen-for-boron substitution, is modelled in regard to its possible usefulness in a nanophotonics device. The modelling is done on both a simple model compound and on a 2D periodic representation of the defect, considering its magnetic and spectroscopic properties. The electronic distribution in $V_N N_B$ excited states is very open-shell in nature, and to deal with this two new computational methods are developed: one allows standard density-functional theory (DFT) calculations to be employed to evaluate state energies, the other introduces techniques needed to apply the VASP computational package to these and many other problems involving excited states. Also of general use, results from DFT calculations are then calibrated against those from ab initio methods, seeking robust computational schemes. These innovations allow 45 electronic states of the defect in its neutral, +1 and -1 charged forms to be considered. The charged forms of the defect are predicted to display properties of potential interest to nanophotonics.

Determination of the chemical and spectroscopic natures of defects in materials such as hexagonal boron nitride (h-BN) remains a serious challenge for both experiment and theory. To establish basics needs for reliable calculations, we consider a model defect $V_N N_B$ in h-BN in which a boron-for-nitrogen substitution is accompanied by a nitrogen vacancy, examining its lowest-energy transition, (1)2B1 to (1)2A1. This provides a relatively simple test system as open-shell and charge-transfer effects, that are difficult to model and can dominate defect spectroscopy, are believed to be small. We establish calculation convergence with respect to sample size using both cluster and 2D-periodic models, convergence with respect to numerical issues such as use of plane-wave or Gaussian-basis-set expansions, and convergence with respect to the treatment of electron correlation. The results strongly suggest that poor performance of computational methods for defects of other natures arise through intrinsic methodological shortcomings.

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN- defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being "molecular" and "materials" in nature and hence bridge the realms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials use only providing further improvements.

Single photon emitters (SPEs) in hexagonal boron nitride (hBN) have garnered significant attention over the last few years due to their superior optical properties. However, despite the vast range of experimental results and theoretical calculations, the defect structure responsible for the observed emission has remained elusive. Here, by controlling the incorporation of impurities into hBN and by comparing various synthesis methods, we provide direct evidence that the visible SPEs are carbon related. Room temperature optically detected magnetic resonance (ODMR) is demonstrated on ensembles of these defects. We also perform ion implantation experiments and confirm that only carbon implantation creates SPEs in the visible spectral range. Computational analysis of hundreds of potential carbon-based defect transitions suggest that the emission results from the negatively charged VBCN- defect, which experiences long-range out-of-plane deformations and is environmentally sensitive. Our results resolve a long-standing debate about the origin of single emitters at the visible range in hBN and will be key to deterministic engineering of these defects for quantum photonic devices.

Crystal defects in the two-dimensional insulator hexagonal boron nitride (hBN) can host localised electronic states that are candidates for applications in quantum technology, yet the precise chemical and structural nature of the defects measured in experiments remains unknown. Recently, optically detected magnetic resonance (ODMR) has been measured from defects in hBN for the first time, providing rich information. In one case a tentative assignment has been made to the VB- defect. Here, we report density-functional theory (DFT) calculations of the energetics, zero field splitting parameter, D, and hyperfine coupling tensor, A, that confirm this assignment. We also show that a second ODMR measurement is consistent with the VN- defect. In both cases, the location of the defect at the edge or centre of layers is important. This advances our atomic-scale understanding of crystal defects in hBN, which is necessary for harnessing the potential of these systems for quantum technology applications., Comment: 17 pages, 4 figures, 1 table

The screening of novel materials is an important topic in the field of materials science. Although traditional computational modeling, especially first-principles approaches, is a very useful and accurate tool to predict the properties of novel materials, it still demands extensive and expensive state-of-the-art computational resources. Additionally, they can be often extremely time consuming. We describe a time and resource-efficient machine learning approach to create a large dataset of structural properties of van der Waals layered structures. In particular, we focus on the interlayer energy and the elastic constant of layered materials composed of two different 2-dimensional (2D) structures, that are important for novel solid lubricant and super-lubricant materials. We show that machine learning models can recapitulate results of computationally expansive approaches (i.e. density functional theory) with high accuracy.

This report summarizes progress made in understanding properties such as zero-phonon-line energies, emission and absorption polarizations, electron-phonon couplings, strain tuning and hyperfine coupling of single photon emitters in hexagonal boron nitride. The primary aims of this research are to discover the chemical nature of the emitting centres and to facilitate deployment in device applications. Critical analyses of the experimental literature and data interpretation, as well as theoretical approaches used to predict properties, are made. In particular, computational and theoretical limitations and challenges are discussed, with a range of suggestions made to overcome these limitations, striving to achieve realistic predictions concerning the nature of emitting centers. A symbiotic relationship is required in which calculations focus on properties that can easily be measured, whilst experiments deliver results in a form facilitating mass-produced calculations., Comment: 43 pages, 7 figures

We show how van der Waals (vdW) forces outcompete covalent and ionic forces to control ferroelectric ordering in CuInP2S6 nanoflakes as well as in CuInP2S6 and CuBiP2Se6 crystals. While the self-assembly of these 2D layered materials is clearly controlled by vdW effects, this result indicates that the internal layer structure is also similarly controlled. Using up to 14 first-principles computational methods, we predict that the bilayers of both materials should be antiferroelectric. However, antiferroelectric nanoflakes and bulk materials are shown to embody two fundamentally different types of inter-layer interactions, with vdW forces strongly favouring one and strongly disfavouring the other compared to ferroelectric ordering. Strong specific vdW interactions involving the Cu atoms control this effect. Thickness-dependent significant cancellation of these two large opposing vdW contributions results in a small net effect that interacts with weak ionic contributions to control ferroelectric ordering.

Following the recent demonstration that van der Waals forces control the ferroelectric ordering of layers within nanoflakes and bulk samples of CuBiP2Se6 and CuInP2S6, it is demonstrated that they also control the internal geometrical structure of isolated monolayers of these materials. This internal structure involves large displacements of the copper atoms, either normal to the layer plane or else within the plane, that change its ligation environment. In both cases, the van der Waals dispersion force out-competes traditional bonding effects to control structure. However, we find that the aspects of the dispersion force giving rise to each effect are uncorrelated: long range effects control inter-layer ferroelectric ordering whereas short-range effects control internal layer structure. These conclusions are drawn considering predicted properties of monolayers, bilayers, and bulk materials obtained using 14 density-functional-theory based methods. While the different methods used often predict starkly different quantitative results, they concur as to the basic nature of ABP2X6 materials. Of the methods used, only the PBE-D3 and optPBEvdW methods were found to predict a wide range of observed properties without serious disparity. Finding optimal computational methods remains a significant challenge for which the unusual multi-scale nature of the van der Waals interactions in ABP2X6 materials provides demanding criteria.

Germanene has emerged as a novel two-dimensional material with various interesting properties and applications. Here we report the possibility of superconductivity in a stable potassium intercalated germanene compound, KGe$_2$, with a transition temperature $T_c \sim 11$ K, and an electron-phonon coupling of 1.9. Applying a 5\% tensile strain, which reduces the buckling height by 4.5\%, leads to the reduction of the electron-phonon coupling by 11\% and a slight increase in $T_c \sim 12$ K. That is, strong electron-phonon coupling results from the buckled structure of the germanene layers. Despite being an intercalated van der Waals material similar to intercalated graphite superconductors, it does not possess an occupied interlayer state.

Key properties of nine possible defect sites in hexagonal boron nitride (h-BN) are predicted using density-functional theory and are corrected by applying results from high-level ab initio calculations. Observed h-BN electron-paramagnetic resonance signals at 22.4, 20.83, and 352.70 MHz are assigned to VN, CN, and VNO2B, respectively, while the observed photoemission at 1.95 eV is assigned to VNCB. Detailed consideration of the available excited states, allowed spin-orbit couplings, zero-field splitting, and optical transitions are made for the two related defects VNCB and VBCN. VNCB is proposed for realizing long-lived quantum memory in h-BN. VBCN is predicted to have a triplet ground state, implying that spin initialization by optical means is feasible and suitable optical excitations are identified, making this defect of interest for possible quantum-qubit operations.

Defect states in 2D materials present many possible uses but both experimental and computational characterization of their spectroscopic properties is difficult. We provide and compare results from 13 DFT and ab initio computational methods for up to 25 excited states of a paradigm system, the VNCB defect in hexagonal boron nitride (h-BN). Studied include include: (i) potentially catastrophic effects for computational methods arising from the multi-reference nature of the closed-shell and open-shell states of the defect, which intrinsically involves broken chemical bonds, (ii) differing results from DFT and time-dependent DFT (TDDFT) calculations, (iii) comparison of cluster models to periodic-slab models of the defect, (iv) the starkly differing effects of nuclear relaxation on the various electronic states as broken bonds try to heal that control the widths of photoabsorption and photoemission spectra, (v) the effect of zero-point energy and entropy on free-energy differences, (vi) defect-localized and conduction/valence band transition natures, and (vii) strategies needed to ensure that the lowest-energy state of a defect can be computationally identified., Comment: https://pubs.acs.org/doi/pdf/10.1021/acs.jctc.7b01072

In 2012, Bjorkman et al. posed the question "Are we van der Waals ready?" [J. Phys.: Condens. Matter, 2012, 24, 424218] about the ability of ab initio modelling to reproduce van der Waals (vdW) dispersion forces in layered materials. The answer at that time was no, however. Here we report on a new generation of vdW dispersion models and show that one, fractionally-ionic atom (FIA) theory, offers close to quantitative predictions for layered structures. Furthermore, it does so from a qualitatively correct picture of dispersion forces. Other methods, such as D3 and optB88vdW also work well, albeit with some exceptions. We thus argue that we are nearly vdW ready, and that some modern dispersion methods are accurate enough to be used for nanomaterial prediction, albeit with some caution required.

Hexagonal boron nitride (hBN) has recently emerged as a fascinating platform for room-temperature quantum photonics due to the discovery of robust visible light single-photon emitters. In order to utilize these emitters, it is necessary to have a clear understanding of their atomic structure and the associated excitation processes that give rise to this single photon emission. Here we perform density-functional theory (DFT) and constrained DFT calculations for a range of hBN point defects in order to identify potential emission candidates. By applying a number of criteria on the electronic structure of the ground state and the atomic structure of the excited states of the considered defects, and then calculating the Huang-Rhys (HR) factor, we find that the CBVN defect, in which a carbon atom substitutes a boron atom and the opposite nitrogen atom is removed, is a potential emission source with a HR factor of 1.66, in good agreement with the experimental HR factor. We calculate the photoluminescence (PL) line shape for this defect and find that it reproduces a number of key features in the the experimental PL lineshape.