DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

In two-dimensional (2D) borophene, the structural transition from triangular lattice to hexagonal lattice with an increase in vacancy concentration is a basic principle of constructing various borophene isomers. Here, by performing an extensive structural search of 4239 borophene isomers with both hexagonal holes (HHs) and large holes (LHs), we show that the structural transformation from triangular lattice to borophene with large holes is energetically more favorable. Borophene isomers with LHs are more stable than those with only HHs at high vacancy concentrations (>20%) and are just slightly less stable than those with only HHs at low vacancy concentrations. This discovery greatly expands the family of 2D borophene and opens a route for synthesizing new borophene isomers.

Low-dimensional carbon and boron nitride nanomaterials - hexagonal boron nitride, graphene, boron nitride nanotubes and carbon nanotubes - remain at the forefront of advanced materials research. Catalytic chemical vapour deposition has become an invaluable technique for reliably and cost-effectively synthesising these materials. In this review, we will emphasise how a synergy between experimental and theoretical methods has enhanced the understanding and optimisation of this synthetic technique. This review examines recent advances in the application of CVD to synthesising boron nitride and carbon nanomaterials and highlights where, in many cases, molecular simulations and quantum chemistry have provided key insights complementary to experimental investigation. This synergy is particularly prominent in the field of carbon nanotube and graphene CVD synthesis, and we propose here it will be the key to future advances in optimisation of CVD synthesis of boron nitride nanomaterials, boron nitride - carbon composite materials, and other nanomaterials generally.

In chemistry, theory of aromaticity or π bond resonance plays a central role in intuitively understanding the stability and properties of organic molecules. Here we present an analogue theory for σ bond resonance in flat boron materials, which allows us to determine the distribution of two-center two-electron and three-center two-electron bonds without quantum calculations. Based on this theory, three rules are proposed to draw the Kekulé-like bonding configurations for flat boron materials and to explore their properties intuitively. As an application of the theory, a simple explanation of why neutral borophene with ~1/9 hole has the highest stability and the effect of charge doping on borophene’s optimal hole concentration is provided with the assumption of σ and π orbital occupation balance. Like the aromaticity theory for carbon materials, this theory greatly deepens our understanding on boron materials and paves the way for the rational design of various boron-based materials.