A new method for calculating optical absorption spectra within linear-scaling density-functional theory (LS-DFT) is presented, incorporating a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states. Three different schemes are compared and the most promising of these, based on the use of a projection operator, has been implemented in a fullyfunctional LS-DFT code. The method has been applied to the calculation of optical absorption spectra for the metal-free phthalocyanine molecule and the conjugated polymer poly(para-phenylene). Excellent agreement with results from a traditional DFT code is obtained.

We present a comparison of methods for treating the electrostatic interactions of finite, isolated systems within periodic boundary conditions (PBCs), within density functional theory (DFT), with particular emphasis on linear-scaling (LS) DFT. Often, PBCs are not physically realistic but are an unavoidable consequence of the choice of basis set and the efficacy of using Fourier transforms to compute the Hartree potential. In such cases the effects of PBCs on the calculations need to be avoided, so that the results obtained represent the open rather than the periodic boundary. The very large systems encountered in LS-DFT make the demands of the supercell approximation for isolated systems more difficult to manage, and we show cases where the open boundary (infinite cell) result cannot be obtained from extrapolation of calculations from periodic cells of increasing size. We discuss, implement, and test three very different approaches for overcoming or circumventing the effects of PBCs: truncation of the Coulomb interaction combined with padding of the simulation cell, approaches based on the minimum image convention, and the explicit use of open boundary conditions (OBCs). We have implemented these approaches in the ONETEP LS-DFT program and applied them to a range of systems, including a polar nanorod and a protein. We compare their accuracy, complexity, and rate of convergence with simulation cell size. We demonstrate that corrective approaches within PBCs can achieve the OBC result more efficiently and accurately than pure OBC approaches.

While monolayer forms of 2D materials are well-characterised both experimentally and theoretically, properties of bilayer heterostructures are not nearly so well-known. We employ high-accuracy linear-scaling DFT calculations utilising non-local van-der-Waals functionals to explore the possible constructions of the MoS2/MoSe2 interface. Utilising large supercells, we vary rotation, translation and separation of the layers without introducing unrealistic strain. The energy landscape shows very low variations under rotation, with no strongly preferred alignments. By unfolding the spectral function into the primitive cells, we show that the monolayers are more independent than in homo-bilayers, and that the electronic bandstructure of each layer is tunable through rotation, thus influencing hole effective masses.

We investigate a recently developed approach [P. L. Silvestrelli, Phys. Rev. Lett. 100, 053002 (2008); J. Phys. Chem. A 113, 5224 (2009)] that uses maximally localized Wannier functions to evaluate the van der Waals contribution to the total energy of a system calculated with density-functional theory. We test it on a set of atomic and molecular dimers of increasing complexity (argon, methane, ethene, benzene, phthalocyanine, and copper phthalocyanine) and demonstrate that the method, as originally proposed, has a number of shortcomings that hamper its predictive power. In order to overcome these problems, we have developed and implemented a number of improvements to the method and show that these modifications give rise to calculated binding energies and equilibrium geometries that are in closer agreement to results of quantum-chemical coupled-cluster calculations.