We have carried out a thorough benchmark of the frozen density-embedding (FDE) method for calculating hole transfer couplings. We have considered 10 exchange-correlation functionals, 3 nonadditive kinetic energy functionals, and 3 basis sets. Overall, we conclude that with a 7% mean relative unsigned error, the PBE and PW91 functionals coupled with the PW91k nonadditive kinetic energy functional and a TZP basis set constitute the most stable and accurate levels of theory for hole transfer coupling calculations. The FDE-ET method is found to be an excellent tool for computing diabatic couplings for hole transfer reactions.

No abstract

Exploiting the machinery of Constrained Density Functional Theory (CDFT), we propose a variational method for calculating low-lying excited states of molecular systems. We dub this method eXcited CDFT (XCDFT). Excited states are obtained by self-consistently constraining a user-defined population of electrons, N, in the virtual space of a reference set of occupied orbitals. By imposing this population to be N = 1.0, we computed the first excited state of 15 molecules from a test set. Our results show that XCDFT achieves an accuracy in the predicted excitation energy only slightly worse than linear-response time-dependent DFT (TDDFT), but without incurring into problems of variational collapse typical of the more commonly adopted ΔSCF method. In addition, we selected a few challenging processes to test the limits of applicability of XCDFT. We find that in contrast to TDDFT, XCDFT is capable of reproducing energy surfaces featuring conical intersections (azobenzene and H) with correct topology and correct overall energetics also away from the intersection. Venturing to condensed-phase systems, XCDFT reproduces the TDDFT solvatochromic shift of benzaldehyde when it is embedded by a cluster of water molecules. Thus, we find XCDFT to be a competitive method among single-reference methods for computations of excited states in terms of time to solution, rate of convergence, and accuracy of the result.

In the past two decades, many research groups worldwide have tried to understand and categorize simple regimes in the charge transfer of such biological systems as DNA. Theoretically speaking, the lack of exact theories for electron-nuclear dynamics on one side, and poor quality of the parameters needed by model Hamiltonians and nonadiabatic dynamics alike (such as couplings and site energies) on the other, are the two main difficulties for an appropriate description of the charge transfer phenomena. In this work, we present an application of a previously benchmarked and linear-scaling subsystem DFT method for the calculation of couplings, site energies and superexchange decay factors (β ) of several biological donor-acceptor dyads, as well as double stranded DNA oligomers comprised of up to 5 base pairs. The calculations are all-electron, and provide a clear view of the role of the environment on superexchange couplings in DNA -they follow experimental trends and confirm previous semiempirical calculations. The subsystem DFT method is proven to be an excellent tool for long-range, bridge-mediated coupling and site energy calculations of embedded molecular systems.

In the present work, we revise and extend the Characteristic Crystallographic Element (CCE) norm, an algorithm used to simultaneously detect radial and orientational similarity of computer-generated structures with respect to specific reference crystals and local symmetries. Based on the identification of point group symmetry elements, the CCE descriptor is able to gauge local structure with high precision and finely distinguish between competing morphologies. As test cases we use computer-generated monomeric and polymer systems of spherical particles interacting with the hard-sphere and square-well attractive potentials. We demonstrate that the CCE norm is able to detect and differentiate, between others, among: hexagonal close packed (HCP), face centered cubic (FCC), hexagonal (HEX) and body centered cubic (BCC) crystals as well as non-crystallographic fivefold (FIV) local symmetry in bulk 3-D systems; triangular (TRI), square (SQU) and honeycomb (HON) crystals, as well as pentagonal (PEN) local symmetry in thin films of one-layer thickness (2-D systems). The descriptor is general and can be applied to identify the symmetry elements of any point group for arbitrary atomic or particulate system in two or three dimensions, in the bulk or under confinement.