DFTB+, a software package for efficient approximate density functional theory based atomistic simulations.

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.

Determining GaN Nanowire Polarity and its Influence on Light Emission in the Scanning Electron Microscope.

The crystal polarity of noncentrosymmetric wurtzite GaN nanowires is determined nondestructively in the scanning electron microscope using electron backscatter diffraction (EBSD). The impact of the nanowire polarity on light emission is then investigated using cathodoluminescence (CL) spectroscopy. EBSD can determine polarity of noncentrosymmetric crystals by interrogating differences in the intensity distribution of bands of the EBSD pattern associated with semipolar planes. Experimental EBSD patterns from an array of GaN nanowires are compared with theoretical patterns produced using dynamical electron simulations to reveal whether they are Ga- or N-polar or, as in several cases, of mixed polarity. CL spectroscopy demonstrates the effect of the polarity on light emission, with spectra obtained from nanowires of known polarity revealing a small but measurable shift (≈28 meV) in the GaN near band edge emission energy between those with Ga and N polarity. We attributed this energy shift to a difference in impurity incorporation in nanowires of different crystal polarity. This approach can be employed to nondestructively identify polarity in a wide range of noncentrosymmetric nanoscale material systems and provide direct comparison with their luminescence.

Plasmon modes in single gold nanodiscs.

Optical properties of single gold nanodiscs were studied by scanning near-field optical microscopy. Near-field transmission spectra of a single nanodisc exhibited multiple plasmon resonances in the visible to near-infrared region. Near-field transmission images observed at these resonance wavelengths show wavy spatial features depending on the wavelength of observation. To clarify physical pictures of the images, theoretical simulations based on spatial correlation between electromagnetic fundamental modes inside and outside of the disc were performed. Simulated images reproduced the observed spatial structures excited in the disc. Mode-analysis of the simulated images indicates that the spatial features observed in the transmission images originate mainly from a few fundamental plasmon modes of the disc.

Formation of Helices in Graphene Nanoribbons under Torsion.

We use objective boundary conditions and self-consistent charge density-functional-based tight-binding to simulate at the atomistic scale the formation of helices in narrow graphene nanoribbons with armchair edges terminated with fluorine and hydrogen. We interpret the microscopic data using an inextensible, unshearable elastic rod model, which considers both bending and torsional strains. When fitted to the atomistic data, the simple rod model uses closed-form solutions for a cubic equation to predict the strain energy and morphology at a given twist angle and the crossover point between pure torsion and a helix. Our modeling and simulation bring key insights into the origin of the helical graphene morphologies stored inside of carbon nanotubes. They can be useful for designing chiral nanoribbons with tailored properties.

Ewald summation on a helix: A route to self-consistent charge density-functional based tight-binding objective molecular dynamics.

We explore the generalization to the helical case of the classical Ewald method, the harbinger of all modern self-consistent treatments of waves in crystals, including ab initio electronic structure methods. Ewald-like formulas that do not rely on a unit cell with translational symmetry prove to be numerically tractable and able to provide the crucial component needed for coupling objective molecular dynamics with the self-consistent charge density-functional based tight-binding treatment of the inter-atomic interactions. The robustness of the method in addressing complex hetero-nuclear nano- and bio-systems is demonstrated with illustrative simulations on a helical boron nitride nanotube, a screw dislocated zinc oxide nanowire, and an ideal DNA molecule.

Rapid nondestructive analysis of threading dislocations in wurtzite materials using the scanning electron microscope.

We describe the use of electron channeling contrast imaging in the scanning electron microscope to rapidly and reliably image and identify threading dislocations (TDs) in materials with the wurtzite crystal structure. In electron channeling contrast imaging, vertical TDs are revealed as spots with black-white contrast. We have developed a simple geometric procedure which exploits the differences observed in the direction of this black-white contrast for screw, edge, and mixed dislocations for two electron channeling contrast images acquired from two symmetrically equivalent crystal planes whose g vectors are at 120° to each other. Our approach allows unambiguous identification of all TDs without the need to compare results with dynamical simulations of channeling contrast.

Geometrical Mie theory for resonances in nanoparticles of any shape.

We give a geometrical theory of resonances in Maxwell's equations that generalizes the Mie formulae for spheres to all scattering channels of any dielectric or metallic particle without sharp edges. We show that the electromagnetic response of a particle is given by a set of modes of internal and scattered fields that are coupled pairwise on the surface of the particle and reveal that resonances in nanoparticles and excess noise in macroscopic cavities have the same origin. We give examples of two types of optical resonances: those in which a single pair of internal and scattered modes become strongly aligned in the sense defined in this paper, and those resulting from constructive interference of many pairs of weakly aligned modes, an effect relevant for sensing. This approach calculates resonances for every significant mode of particles, demonstrating that modes can be either bright or dark depending on the incident field. Using this extra mode information we then outline how excitation can be optimized. Finally, we apply this theory to gold particles with shapes often used in experiments, demonstrating effects including a Fano-like resonance.

DFTB+, a sparse matrix-based implementation of the DFTB method.

A new Fortran 95 implementation of the DFTB (density functional-based tight binding) method has been developed, where the sparsity of the DFTB system of equations has been exploited. Conventional dense algebra is used only to evaluate the eigenproblems of the system and long-range Coulombic terms, but drop-in O(N) or O(N2) modules are planned to replace the small code sections that these entail. The developed sparse storage structure is discussed in detail, and a short overview of other features of the new code is given.

Self-interaction and strong correlation in DFTB.

The density functional based tight-binding (DFTB) method can benefit substantially from a number of developments in density functional theory (DFT) while also providing a simple analytical proving ground for new extensions. This contribution begins by demonstrating the variational nature of charge-self-consistent DFTB (SCC-DFTB), proving the presence of a defined ground-state in this class of methods. Because the ground state of the SCC-DFTB method itself can be qualitatively incorrect for some systems, suitable forms of the recent LDA+U functionals for SCC-DFTB are also presented. This leads to both a new semilocal self-interaction correction scheme and a new physical argument for the choice of parameters in the LDA+U method. The locality of these corrections can only partly repair the HOMO-LUMO gap and chemical potential discontinuity, hence a novel method for introducing this further physics into the method is also presented, leading to exact derivative discontinuities in this theory at low computational cost. The prototypical system NiO is used as an illustration for these developments.

Initial steps toward automating the fitting of DFTB Erep(r).

The most time-consuming part of developing new parametrizations for the density functional based tight-binding (DFTB) method consists of producing accurate and transferable repulsive pair potentials. In the conventional approach to repulsive parametrization, every possible diatomic combination of the elements covered by the parametrization must be individually hand-constructed. We present an initial attempt to automate some of this time-consuming process. We consider a simple genetic algorithm-based approach to the fitting problem.

An efficient LDA+U based tight binding approach.

The functionals usually applied in DFT calculations have deficiencies in describing systems with strongly localized electrons such as transition metals or rare earth (RE) compounds. In this work, we present the self-consistent charge density based functional tight binding (SCC-DFTB) calculation scheme including LDA+U like potentials and apply it for the simulation of RE-doped GaN. DFTB parameters for the simulation of GaN and a selection of rare earth ions, where the f electrons were explicitly included in the valence, have been created. The results of the simulations were tested against experimental data (where present) and against various more sophisticated but computationally more costly DFT calculations. Our approach is found to correctly reproduce the geometry and the energetic of the studied systems.