LOBSTER: Local orbital projections,atomic charges,and chemical-bonding analysis from projector-augmented-wave-based density-functional theory

We present an update on recently developed methodology and functionality in the computer program Local Orbital Basis Suite Toward Electronic-Structure Reconstruction (LOBSTER) for chemical-bonding analysis in periodic systems. LOBSTER is based on an analytic projection from projector-augmented wave (PAW) density-functional theory (DFT) computations (Maintz et al., J. Comput. Chem. 2013 , 34, 2557), reconstructing chemical information in terms of local, auxiliary atomic orbitals and thereby opening the output of PAW-based DFT codes to chemical interpretation. We demonstrate how LOBSTER has been improved by taking into account time-reversal symmetry, thereby speeding up the DFT and LOBSTER calculations by a factor of 2. Over the recent years, the functionalities have also been continually expanded, including accurate projected densities of states (DOSs), crystal orbital Hamilton population (COHP) analysis, atomic and orbital charges, gross populations, and the recently introduced k -dependent COHP. The software is offered free-of-charge for non-commercial research.     

The 5v 3 bands of 18O enriched ozone: line positions of 16O16O18O, 16O18O16O, 16O18O18O and 18O16O18O

The four 5v ₃ bands of ¹⁸O enriched ozone have been observed and analysed for the first time. Two species (¹⁶O¹⁸O¹⁶O and ¹⁸O¹⁶O¹⁸O) belong to the C_2v symmetry group and two other (¹⁸O¹⁸O¹⁶O and ¹⁶O¹⁶O¹⁸O) to the C_s symmetry group. They have been recorded at a resolution of 0.008 cm^?1 with a pathlength of 32.16 m. Despite the very weak absorptions observed, almost 250 energy levels have been derived for each of the 4 species, with J ? 35 and K _a ? 13, and suitable sets of Hamiltonian parameters have been determined. For 3 species it has been necessary to account for the resonance between the (005) and (311) states to correctly reproduce the spectra observed. These resonances, anharmonic for C_2v, and hybrid (both anhar-monic and Coriolis) for C_s symmetry confirm the accidentally extremely strong coupling between the (005) and (311) states for ¹⁶O₃, due in that case to the very close distance between unperturbed energy levels. This work also confirms the excellent prediction of band centres of these four species derived from the recently determined isotopically invariant molecular potential function.     

On the DFT Ground State of Crystalline Bromine and Iodine

We report on an erroneous ground state within common density functional theory (DFT) methods for the solid elements bromine and iodine. Phonon computations at the GGA level for both molecular crystals yield imaginary vibrational modes, erroneously indicating dynamic instability—that fact alone could easily pass as a computational artefact, but these imaginary modes lead to energetically more favorable and dynamically stable structures, made up of infinite monoatomic chains. In contrast, meta‐GGA and hybrid functionals yield the correct energetic order for bromine, while for iodine, most global hybrids do not improve the GGA result significantly. The qualitatively correct answer, in both cases, is given by the long‐range corrected hybrid LC‐ωPBE, the Minnesota functionals M06L and M06, and by periodic Hartree–Fock and MP2 theory. This poor performance of economic DFT functionals should be kept in mind, for example, during global structure optimizations of systems with significant contributions from halogen bonds.     

Hierarchically Structured Allotropes of Phosphorus from Data-Driven Exploration

The discovery of materials is increasingly guided by quantum-mechanical crystal-structure prediction, but the structural complexity in bulk and nanoscale materials remains a bottleneck. Here we demonstrate how data-driven approaches can vastly accelerate the search for complex structures, combining a machine-learning (ML) model for the potential-energy surface with efficient, fragment-based searching. We use the characteristic building units observed in Hittorf's and fibrous phosphorus to seed stochastic (“random”) structure searches over hundreds of thousands of runs. Our study identifies a family of hierarchically structured allotropes based on a P8 cage as principal building unit, including one-dimensional (1D) single and double helix structures, nanowires, and two-dimensional (2D) phosphorene allotropes with square-lattice and kagome topologies. These findings yield new insight into the intriguingly diverse structural chemistry of phosphorus, and they provide an example for how ML methods may, in the long run, be expected to accelerate the discovery of hierarchical nanostructures.     

Hydrogen-bonding networks from first-principles: exploring the guanidine crystal

Hydrogen bonding is among the most important interactions in molecular crystals, and examples are abundant. As a consequence of such interactions, many molecules crystallize in complex but intriguing structures, in contrast to the relatively simple packing principles of metallic or ionic solids. In this work, we present a computational approach based on plane-wave density-functional theory (DFT) and supercell techniques, aiming to understand and quantify hydrogen-bonded networks in the solid state and in two-, one-, and zero-dimensional fragments derived from the molecular crystal. With such methodology at hand, we investigate guanidine, a fitting example of a molecular crystal and an important compound for inorganic and organic chemistry alike. On the basis of our computations, we discuss the initially proposed layered structure of guanidine and identify both stabilizing and destabilizing cooperative interactions in the three crystalline dimensions.     

LOBSTER: A tool to extract chemical bonding from plane‐wave based DFT

The computer program LOBSTER (Local Orbital Basis Suite Towards Electronic‐Structure Reconstruction) enables chemical‐bonding analysis based on periodic plane‐wave (PAW) density‐functional theory (DFT) output and is applicable to a wide range of first‐principles simulations in solid‐state and materials chemistry. LOBSTER incorporates analytic projection routines described previously in this very journal [J. Comput. Chem. 2013 , 34, 2557] and offers improved functionality. It calculates, among others, atom‐projected densities of states (pDOS), projected crystal orbital Hamilton population (pCOHP) curves, and the recently introduced bond‐weighted distribution function (BWDF). The software is offered free‐of‐charge for non‐commercial research. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

From Atomistic Surface Chemistry to Nanocrystals of Functional Chalcogenides

Synthesis and utilization of nanocrystals are highly active fields of current research, but they require a thorough understanding of the underlying crystal surfaces. In this Minireview, we span the arc from surfaces to free nanocrystals, and onward to their chemical synthesis, using as examples lead selenide (PbSe), tin telluride (SnTe), and their direct chemical relatives. Besides experimental insights, we highlight the increasingly influential role played by quantum‐chemical simulations of surfaces and nanocrystals. What can theory do today, or possibly tomorrow; where are its limits? Answering these questions, and skillfully linking them to experiments, could open up new atomistically (that is, chemically) guided perspectives for nanosynthesis.