Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: Obtaining efficiency and accuracy with in situ optimised local orbitals

Can domain-based local pair natural orbitals approaches accurately predict phosphorescence energies?

Since the discovery of the peculiar conducting and optical properties of aromatics, many efforts have been made to characterize and predict their phosphorescence. This physical process is exploited in modern Organic Emitting Light Diodes (OLEDs), and it is also one of the processes decreasing the efficiency of Dye-sensitized solar cells (DSSCs). Herein, we propose a computational strategy for the accurate calculation of singlet-triplet gaps of aromatic compounds, which provides results that are in excellent agreement with available experimental data. Our approach relies on the domain-based local pair natural orbital (DLPNO) variant of the "gold standard" CCSD(T) method. The convergence of our results with respect to the key technical parameters of the calculation, such as the basis set used, the approximations employed in the perturbative triples correction, and the dimension of the PNOs space, was thoroughly discussed.

Transition-Based Constrained DFT for the Robust and Reliable Treatment of Excitations in Supramolecular Systems

Despite the variety of available computational approaches, state-of-the-art methods for calculating excitation energies, such as time-dependent density functional theory (TDDFT), are computationally demanding and thus limited to moderate system sizes. Here, we introduce a new variation of constrained DFT (CDFT), wherein the constraint corresponds to a particular transition (T), or a combination of transitions, between occupied and virtual orbitals, rather than a region of the simulation space as in traditional CDFT. We compare T-CDFT with TDDFT and ΔSCF results for the low-lying excited states (S1 and T1) of a set of gas-phase acene molecules and OLED emitters and with reference results from the literature. At the PBE level of theory, T-CDFT outperforms ΔSCF for both classes of molecules, while also proving to be more robust. For the local excitations seen in the acenes, T-CDFT and TDDFT perform equally well. For the charge transfer (CT)-like excitations seen in the OLED molecules, T-CDFT also performs well, in contrast to the severe energy underestimation seen with TDDFT. In other words, T-CDFT is equally applicable to both local excitations and CT states, providing more reliable excitation energies at a much lower computational cost than TDDFT cost. T-CDFT is designed for large systems and has been implemented in the linear-scaling BigDFT code. It is therefore ideally suited for exploring the effects of explicit environments on excitation energies, paving the way for future simulations of excited states in complex realistic morphologies, such as those which occur in OLED materials.

Efficiently Computing Excitations of Complex Systems: Linear-Scaling Time-Dependent Embedded Mean-Field Theory in Implicit Solvent

Quantum embedding schemes have the potential to significantly reduce the computational cost of first-principles calculations while maintaining accuracy, particularly for calculations of electronic excitations in complex systems. In this work, I combine time-dependent embedded mean field theory (TD-EMFT) with linear-scaling density functional theory and implicit solvation models, extending previous work within the ONETEP code. This provides a way to perform multilevel calculations of electronic excitations on very large systems, where long-range environmental effects, both quantum and classical in nature, are important. I demonstrate the power of this method by performing simulations on a variety of systems, including a molecular dimer, a chromophore in solution, and a doped molecular crystal. This work paves the way for high accuracy calculations to be performed on large-scale systems that were previously beyond the reach of quantum embedding schemes.

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Rapid Prediction of Ultraviolet-Visible Spectra from Conventional (Non-Time-Dependent) Density Functional Theory Calculations

We present an approximate approach for the simulation of UV/vis spectra using conventional [non-time-dependent (non-TD)] DFT computations. It uses Kohn-Sham orbitals and orbital energies to estimate both the excitation energies and the associated oscillator strengths. For a wide range of systems from small molecules to large molecular dyes used in electrochromic and solar-cell applications, reasonable UV/vis spectra are generated, each with just two conventional DFT computations. The accuracy is generally comparable to what one would expect from TD-DFT calculations. In comparison to TD-DFT, the protocol of the present study provides an intuitive and notably more rapid means for simulating electronic absorption properties. It enables efficient screening of materials for a wide range of relevant applications.

Machine Learning for Absorption Cross Sections

We present a machine learning (ML) method to accelerate the nuclear ensemble approach (NEA) for computing absorption cross sections. ML-NEA is used to calculate cross sections on vast ensembles of nuclear geometries to reduce the error due to insufficient statistical sampling. The electronic properties-excitation energies and oscillator strengths-are calculated with a reference electronic structure method only for a relatively few points in the ensemble. The KREG model (kernel-ridge-regression-based ML combined with the RE descriptor) as implemented in MLatom is used to predict these properties for the remaining tens of thousands of points in the ensemble without incurring much of additional computational cost. We demonstrate for two examples, benzene and a 9-dicyanomethylene derivative of acridine, that ML-NEA can produce statistically converged cross sections even for very challenging cases and even with as few as several hundreds of training points.

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.

Neutral excitation density-functional theory: an efficient and variational first-principles method for simulating neutral excitations in molecules

We introduce neutral excitation density-functional theory (XDFT), a computationally light, generally applicable, first-principles technique for calculating neutral electronic excitations. The concept is to generalise constrained density functional theory to free it from any assumptions about the spatial confinement of electrons and holes, but to maintain all the advantages of a variational method. The task of calculating the lowest excited state of a given symmetry is thereby simplified to one of performing a simple, low-cost sequence of coupled DFT calculations. We demonstrate the efficacy of the method by calculating the lowest single-particle singlet and triplet excitation energies in the well-known Thiel molecular test set, with results which are in good agreement with linear-response time-dependent density functional theory (LR-TDDFT). Furthermore, we show that XDFT can successfully capture two-electron excitations, in principle, offering a flexible approach to target specific effects beyond state-of-the-art adiabatic-kernel LR-TDDFT. Overall the method makes optical gaps and electron-hole binding energies readily accessible at a computational cost and scaling comparable to that of standard density functional theory. Owing to its multiple qualities beneficial to high-throughput studies where the optical gap is of particular interest; namely broad applicability, low computational demand, and ease of implementation and automation, XDFT presents as a viable candidate for research within materials discovery and informatics frameworks.

The ONETEP linear-scaling density functional theory program

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.

Computation of Electromagnetic Properties of Molecular Ensembles

We outline a methodology for efficiently computing the electromagnetic response of molecular ensembles. The methodology is based on the link that we establish between quantum-chemical simulations and the transfer matrix (T-matrix) approach, a common tool in physics and engineering. We exemplify and analyze the accuracy of the methodology by using the time-dependent Hartree-Fock theory simulation data of a single chiral molecule to compute the T-matrix of a cross-like arrangement of four copies of the molecule, and then computing the circular dichroism of the cross. The results are in very good agreement with full quantum-mechanical calculations on the cross. Importantly, the choice of computing circular dichroism is arbitrary: Any kind of electromagnetic response of an object can be computed from its T-matrix. We also show, by means of another example, how the methodology can be used to predict experimental measurements on a molecular material of macroscopic dimensions. This is possible because, once the T-matrices of the individual components of an ensemble are known, the electromagnetic response of the ensemble can be efficiently computed. This holds for arbitrary arrangements of a large number of molecules, as well as for periodic or aperiodic molecular arrays. We identify areas of research for further improving the accuracy of the method, as well as new fundamental and technological research avenues based on the use of the T-matrices of molecules and molecular ensembles for quantifying their degrees of symmetry breaking. We provide T-matrix-based formulas for computing traditional chiro-optical properties like (oriented) circular dichroism, and also for quantifying electromagnetic duality and electromagnetic chirality. The formulas are valid for light-matter interactions of arbitrarily-high multipolar orders.

Predicting solvatochromic shifts and colours of a solvated organic dye: The example of nile red

The solvatochromic shift, as well as the change in colour of the simple organic dye nile red, is studied in two polar and two non-polar solvents in the context of large-scale time-dependent density-functional theory (TDDFT) calculations treating large parts of the solvent environment from first principles. We show that an explicit solvent representation is vital to resolve absorption peak shifts between nile red in n-hexane and toluene, as well as acetone and ethanol. The origin of the failure of implicit solvent models for these solvents is identified as being due to the strong solute-solvent interactions in form of π-stacking and hydrogen bonding in the case of toluene and ethanol. We furthermore demonstrate that the failures of the computationally inexpensive Perdew-Burke-Ernzerhof (PBE) functional in describing some features of the excited state potential energy surface of the S₁ state of nile red can be corrected for in a straightforward fashion, relying only on a small number of calculations making use of more sophisticated range-separated hybrid functionals. The resulting solvatochromic shifts and predicted colours are in excellent agreement with experiment, showing the computational approach outlined in this work to yield very robust predictions of optical properties of dyes in solution.

Distance dependent photoacoustics revealed through DNA nanostructures

Molecular rulers that rely on the Förster resonance energy transfer (FRET) mechanism are widely used to investigate dynamic molecular processes that occur on the nanometer scale. However, the capabilities of these fluorescence molecular rulers are fundamentally limited to shallow imaging depths by light scattering in biological samples. Photoacoustic tomography (PAT) has recently emerged as a high resolution modality for in vivo imaging, coupling optical excitation with ultrasound detection. In this paper, we report the capability of PAT to probe distance-dependent FRET at centimeter depths. Using DNA nanotechnology we created several nanostructures with precisely positioned fluorophore-quencher pairs over a range of nanoscale separation distances. PAT of the DNA nanostructures showed distance-dependent photoacoustic signal enhancement and demonstrated the ability of PAT to reveal the FRET process deep within tissue mimicking phantoms. Further, we experimentally validated these DNA nanostructures as a novel and biocompatible strategy to augment the intrinsic photoacoustic signal generation capabilities of small molecule fluorescent dyes.

Unravelling the Roles of Size, Ligands, and Pressure in the Piezochromic Properties of CdS Nanocrystals

Understanding the effects of pressure-induced deformations on the optoelectronic properties of nanomaterials is important not only from the fundamental point of view but also for potential applications such as stress sensors and electromechanical devices. Here, we describe the novel insights into these piezochromic effects gained from using a linear-scaling density functional theory framework and an electronic enthalpy scheme, which allow us to accurately characterize the electronic structure of CdS nanocrystals with a zincblende-like core of experimentally relevant size. In particular, we focus on unravelling the complex interplay of size and surface (phenyl) ligands with pressure. We show that pressure-induced deformations are not simple isotropic scaling of the original structures and that the change in HOMO-LUMO gap with pressure results from two competing factors: (i) a bulk-like linear increase due to compression, which is offset by (ii) distortions and disorder and, to a lesser extent, orbital hybridization induced by ligands affecting the frontier orbitals. Moreover, we observe that the main peak in the optical absorption spectra is systematically red-shifted or blue-shifted, as pressure is increased up to 5 GPa, depending on the presence or absence of phenyl ligands. These heavily hybridize the frontier orbitals, causing a reduction in overlap and oscillator strength, so that at zero pressure, the lowest energy transition involves deeper hole orbitals than in the case of hydrogen-capped nanocrystals; the application of pressure induces greater delocalization over the whole nanocrystals bringing the frontier hole orbitals into play and resulting in an unexpected red shift for the phenyl-capped nanocrystals, in part caused by distortions. In response to a growing interest in relatively small nanocrystals that can be difficult to accurately characterize with experimental techniques, this work exemplifies the detailed understanding of structure-property relationships under pressure that can be obtained for realistic nanocrystals with state-of-the-art first-principles methods and used for the characterization and design of devices based on these and similar nanomaterials.

Linear-scaling density functional theory using the projector augmented wave method

Quantum mechanical simulation of realistic models of nanostructured systems, such as nanocrystals and crystalline interfaces, demands computational methods combining high-accuracy with low-order scaling with system size. Blöchl's projector augmented wave (PAW) approach enables all-electron (AE) calculations with the efficiency and systematic accuracy of plane-wave pseudopotential calculations. Meanwhile, linear-scaling (LS) approaches to density functional theory (DFT) allow for simulation of thousands of atoms in feasible computational effort. This article describes an adaptation of PAW for use in the LS-DFT framework provided by the ONETEP LS-DFT package. ONETEP uses optimisation of the density matrix through in situ-optimised local orbitals rather than the direct calculation of eigenstates as in traditional PAW approaches. The method is shown to be comparably accurate to both PAW and AE approaches and to exhibit improved convergence properties compared to norm-conserving pseudopotential methods.

Applications of large-scale density functional theory in biology

Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.