Linear-scaling time-dependent density-functional theory in the linear response formalism

A Resolution of Identity Technique to Speed up TDDFT with Hybrid Functionals: Implementation and Application to the Magic Cluster Series Au8n+4(SC6H5)4n+8 (n = 3-6)

The Resolution of Identity (RI) technique has been employed to speed up the use of hybrid exchange-correlation (xc) functionals at the TDDFT level using the Hybrid Diagonal Approximation. The RI has been implemented within the polTDDFT algorithm (a complex damped polarization method) in the AMS/ADF suite of programs. A speedup factor of 30 has been obtained with respect to a previous numerical implementation, albeit with the same level of accuracy. Comparison of TDDFT simulations with the experimental photoabsorption spectra of the cluster series Au8n+4(SR)4n+8(n = 3-6; R = C6H5) showed the excellent accuracy and efficiency of the method. Results were compared with those obtained via the more simplified and computationally cheaper TDDFT+TB and sTDDFT methods. The present method represents an accurate as well as computationally affordable approach to predict photoabsorption spectra of complex species, realizing an optimal compromise between accuracy and computational efficiency, and is suitable for applications to large metal clusters with sizes up to several hundreds of atoms.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods

Quantum algorithms for simulating electronic ground states are slower than popular classical mean-field algorithms such as Hartree-Fock and density functional theory but offer higher accuracy. Accordingly, quantum computers have been predominantly regarded as competitors to only the most accurate and costly classical methods for treating electron correlation. However, here we tighten bounds showing that certain first-quantized quantum algorithms enable exact time evolution of electronic systems with exponentially less space and polynomially fewer operations in basis set size than conventional real-time time-dependent Hartree-Fock and density functional theory. Although the need to sample observables in the quantum algorithm reduces the speedup, we show that one can estimate all elements of the k-particle reduced density matrix with a number of samples scaling only polylogarithmically in basis set size. We also introduce a more efficient quantum algorithm for first-quantized mean-field state preparation that is likely cheaper than the cost of time evolution. We conclude that quantum speedup is most pronounced for finite-temperature simulations and suggest several practically important electron dynamics problems with potential quantum advantage.

Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules

We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one s-type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra. Benchmarked on organic molecules and compared to standard TDDFT, TDDFT-ris has an average energy error of only 0.06 eV and yields absorption spectra in close agreement with TDDFT. Thus, TDDFT-ris enables simulation of realistic absorption spectra in large molecules that would be inaccessible from standard TDDFT.

Accelerating Linear-Response Time-Dependent Hybrid Density Functional Theory with Low-Rank Decomposition Techniques in the Plane-Wave Basis

We present an efficient low-rank implementation of linear-response time-dependent density functional theory for hybrid functionals (hybrid-LR-TDDFT) within the plane-wave pseudopotential framework. The adaptively compressed exchange (ACE) operator and the natural transition orbitals (NTOs) are introduced to build the low-rank representation of the nonlocal exchange operator in the hybrid-LR-TDDFT Hamiltonian. Numerical tests demonstrate that the ACE approximation significantly reduces the computational cost of applying the nonlocal exchange operator without loss of accuracy, and the NTO approximation can further accelerate the hybrid-LR-TDDFT calculations by introducing an NTO cutoff parameter. This new method enables us to effectively study the excitonic properties of two-dimensional MoS₂ consisting of 216 atoms and ∼1900 electrons with range-separated hybrid functionals on a single graphics processing unit.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.

Orbital-Free Methods for Plasmonics: Linear Response

Plasmonic systems, such as metal nanoparticles, are widely used in different application areas, going from biology to photovoltaics.The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties.When the characteristic sizes/distances reach a few nanometers, non-local and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate. Methods based on the Time-Dependent Density-Functional Theory (TD-DFT) represent the current reference for the description of quantum effects. However, TD-DFT is based on knowledge of all occupied orbitals whose calculation is computationally prohibitive to model large plasmonic systems of interest for applications.On the other hand, methods based on the Orbital-Free (OF) formulation of TD-DFT, can scale linearly with the system size.In this Review, OF methods ranging from semiclassical models to the quantum hydrodynamic theory, will be derived from the linear response TD-DFT, so that the key approximations and properties of each method can be clearly highlighted. The accuracy of the various approximations will be then validated for the linear optical properties of jellium nanoparticles, the most relevant model system in plasmonics. OF methods can describe the collective excitations in plasmonic systems with great accuracy andwithout system-tuned parameters. The accuracy on these methods depends only on the accuracy on the (universal) kinetic energy functional of the ground-state electronic density. Current approximations and future development directions will be indicated.

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.

Plasmon Couplings from Subsystem Time-Dependent Density Functional Theory

Many applications in plasmonics are related to the coupling between metallic nanoparticles (MNPs) or between an emitter and a MNP. The theoretical analysis of such a coupling is thus of fundamental importance to analyze the plasmonic behavior and to design new systems. While classical methods neglect quantum and spill-out effects, time-dependent density functional theory (TD-DFT) considers all of them and with Kohn-Sham orbitals delocalized over the whole system. Thus, within TD-DFT, no definite separation of the subsystems (the single MNP or the emitter) and their couplings is directly available. This important feature is obtained here using the subsystem formulation of TD-DFT, which has been originally developed in the context of weakly interacting organic molecules. In subsystem TD-DFT, interacting MNPs are treated independently, thus allowing us to compute the plasmon couplings directly from the subsystem TD-DFT transition densities. We show that subsystem TD-DFT, as well as a simplified version of it in which kinetic contributions are neglected, can reproduce the reference TD-DFT calculations for gap distances greater than about 6 Å or even smaller in the case of hybrid plasmonic systems (i.e., molecules interacting with MNPs). We also show that the subsystem TD-DFT can be also used as a tool to analyze the impact of charge-transfer effects.

Accurate and Efficient Computation of Optical Absorption Spectra of Molecular Crystals: The Case of the Polymorphs of ROY

When calculating the optical absorption spectra of molecular crystals from first principles, the influence of the crystalline environment on the excitations is of significant importance. For such systems, however, methods to describe the excitations accurately can be computationally prohibitive due to the relatively large system sizes involved. In this work, we demonstrate a method that allows optical absorption spectra to be computed both efficiently and at high accuracy. Our approach is based on the spectral warping method successfully applied to molecules in solvent. It involves calculating the absorption spectrum of a supercell of the full molecular crystal using semi-local time-dependent density functional theory (TDDFT), before warping the spectrum using a transformation derived from smaller-scale semi-local and hybrid TDDFT calculations on isolated dimers. We demonstrate the power of this method on three polymorphs of the well-known color polymorphic compound ROY and find that it outperforms both small-scale hybrid TDDFT dimer calculations and large-scale semi-local TDDFT supercell calculations, when compared to the experiment.

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.

Two-level iterative solver for linear response time-dependent density functional theory with plane wave basis set

We present a two-level iterative solver for linear response time-dependent density functional theory by combining two forms of the Casida equation in the Kohn-Sham orbital representation and in Hutter's formulation. This two-level iterative solver has been implemented with the plane wave pseudopotential method for excited-state simulations of molecular and low-dimensional solid materials. Numerical studies with the Davidson algorithm demonstrate that this two-level iterative solver yields excited-state properties for molecules (benzene C₆H₆ and fullerene C₆₀) and low-dimensional semiconductors [two-dimensional molybdenum disulfide MoS₂ monolayer and rutile titanium dioxide TiO₂(110) surface] with significantly reduced computational cost and storage requirement compared with standard iterative algorithms. We apply our approach to investigate the photoinduced charge separation of methanol molecules adsorption on the rutile TiO₂(110) surface from the exciton perspective and validate that the photogenerated hole can be captured by methanol molecules.

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.

An efficient hybrid scheme for time dependent density functional theory

A hybrid approach able to perform Time Dependent Density Functional Theory (TDDFT) simulations with the same accuracy as that of hybrid exchange-correlation (xc-) functionals but at a fraction of the computational cost is developed, implemented, and validated. The scheme, which we name Hybrid Diagonal Approximation (HDA), consists in employing in the response function a hybrid xc-functional (containing a fraction of the non-local Hartree-Fock exchange) only for the diagonal elements of the omega matrix, while the adiabatic local density approximation is employed for the off-diagonal terms. HDA is especially (but not exclusively) advantageous when using Slater type orbital basis sets and allows one to employ them in a uniquely efficient way, as we demonstrate here by implementing HDA in a local version of the Amsterdam Density Functional code. The new protocol is tested on NH₃, C₆H₆, and the [Au₂₅(SCH₃)₁₈]^- cluster as prototypical cases ranging from small molecules to ligand-protected metal clusters, finding excellent agreement with respect to both full kernel TDDFT simulations and experimental data. Additionally, a specific comparison test between full kernel and HDA is considered at the Casida level on seven other molecular species, which further confirm the accuracy of the approach for all investigated systems. For the [Au₂₅(SCH₃)₁₈]^- cluster, a speedup by a factor of seven is obtained with respect to the full kernel. The HDA, therefore, promises to provide a quantitative description of the optical properties of medium-sized systems (nanoclusters) at an affordable cost, thanks to its computational efficiency, especially in combination with a complex polarization algorithm version of TDDFT.

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.

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.

Relativistic time-dependent density functional theories

The foundations, formalisms, technicalities, and practicalities of relativistic time-dependent density functional theories (R-TD-DFT) for spinor excited states of molecular systems containing heavy elements are critically reviewed.

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.

Time-Dependent Coupled Perturbed Hartree-Fock and Density-Functional-Theory Approach for Calculating Frequency-Dependent (Hyper)Polarizabilities with Nonorthogonal Localized Molecular Orbitals

The time-dependent coupled perturbed Hartree-Fock/density-functional-theory (TDHF/TDDFT) approach has been reformulated based on nonorthogonal localized molecular orbitals (NOLMOs). Based on the NOLMO Fock equation, we have derived the corresponding NOLMO-TDHF/TDDFT equations up to the third order, and the formula for the frequency-dependent (hyper)polarizabilities has been given. Our approach has been applied to calculate both static and dynamic (hyper)polarizabilities of molecules varying from small molecules to large molecules. The NOLMO-TDHF/TDDFT approach can reproduce the reference canonical molecular orbital (CMO) results for all of our testing calculations. With the help of ongoing development of optimized local virtual molecular orbitals, the NOLMO-TDHF/TDDFT approach would be a very efficient method for large system calculations and tp achieve linear scaling.

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.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

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.

The potential of imogolite nanotubes as (co-)photocatalysts: a linear-scaling density functional theory study

We report a linear-scaling density functional theory (DFT) study of the structure, wall-polarization absolute band-alignment and optical absorption of several, recently synthesized, open-ended imogolite (Imo) nanotubes (NTs), namely single-walled (SW) aluminosilicate (AlSi), SW aluminogermanate (AlGe), SW methylated aluminosilicate (AlSi-Me), and double-walled (DW) AlGe NTs. Simulations with three different semi-local and dispersion-corrected DFT-functionals reveal that the NT wall-polarization can be increased by nearly a factor of four going from SW-AlSi-Me to DW-AlGe. Absolute vacuum alignment of the NT electronic bands and comparison with those of rutile and anatase TiO2 suggest that the NTs may exhibit marked propensity to both photo-reduction and hole-scavenging. Characterization of the NTs' band-separation and optical properties reveal the occurrence of (near-)UV inside-outside charge-transfer excitations, which may be effective for electron-hole separation and enhanced photocatalytic activity. Finally, the effects of the NTs' wall-polarization on the absolute alignment of electron and hole acceptor states of interacting water (H2O) molecules are quantified and discussed.

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

We present a solution of the full time-dependent density-functional theory (TDDFT) eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspaces with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate gradient algorithm that is very memory-efficient. The algorithm is validated on a small test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of bacteriochlorophyll in solution. The largest systems considered in this work are of the same order of magnitude as a variety of widely studied pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.

Linear scaling density matrix real time TDDFT: Propagator unitarity and matrix truncation

Real time, density matrix based, time dependent density functional theory (TDDFT) proceeds through the propagation of the density matrix, as opposed to the Kohn-Sham orbitals. It is possible to reduce the computational workload by imposing spatial cutoff radii on sparse matrices, and the propagation of the density matrix in this manner provides direct access to the optical response of very large systems, which would be otherwise impractical to obtain using the standard formulations of TDDFT. Following a brief summary of our implementation, along with several benchmark tests illustrating the validity of the method, we present an exploration of the factors affecting the accuracy of the approach. In particular, we investigate the effect of basis set size and matrix truncation, the key approximation used in achieving linear scaling, on the propagator unitarity and optical spectra. Finally, we illustrate that, with an appropriate density matrix truncation range applied, the computational load scales linearly with the system size and discuss the limitations of the approach.

Coupled-Perturbed SCF Approach for Calculating Static Polarizabilities and Hyperpolarizabilities with Nonorthogonal Localized Molecular Orbitals

Coupled-perturbed self-consistent-field (CPSCF) approach has been broadly used for polarizabilities and hyperpolarizabilities computation. To extend this application to large systems, we have reformulated the CPSCF equations with nonorthogonal localized molecular orbitals (NOLMOs). NOLMOs are the most localized representation of electronic degrees of freedom. Methods based on NOLMOs are potentially ideal for investigating large systems. In atomic orbital representation, with a static external electric field added, the wave function and SCF operator of unperturbed NOLMO-SCF wave function/orbitals are expanded to different orders of perturbations. We have derived the corresponding equations up to the third order, which are significantly different from those of a conventional CPSCF method because of the release of the orthogonal restrictions on MOs. The solution to these equations has been implemented. Several chemical systems are used to verify our method. This work represents the first step toward efficient calculations of molecular response and excitation properties with NOLMOs.

Identifying and tracing potential energy surfaces of electronic excitations with specific character via their transition origins: application to oxirane

We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a wide range of molecular geometries.

Toward Ab Initio Optical Spectroscopy of the Fenna-Matthews-Olson Complex

We present progress toward a first-principles parametrization of the Hamiltonian of the Fenna-Matthews-Olson pigment-protein complex, a molecule that has become key to understanding the role of quantum dynamics in photosynthetic exciton energy transfer. To this end, we have performed fully quantum mechanical calculations on each of the seven bacteriochlorophyll pigments that make up the complex, including a significant proportion of their protein environment (more than 2000 atoms), using linear-scaling density functional theory exploiting a recent development for the computation of excited states. Local pigment transition energies and interpigment coupling between optical transitions have been calculated and are in good agreement with the literature consensus. Comparisons between simulated and experimental optical spectra point toward future work that may help to elucidate important design principles in these nanoscale devices.