Person:

Olivares-Amaya, Roberto

Excess charge stabilization of molecules in metallic environments is of particular importance for fields such as molecular electronics and surface chemistry. We study the energetics of benzene and its anion between two metallic plates. We observe that orientational effects are important at small inter-plate separation. This leads to benzene oriented perpendicular to the gates being more stable than the parallel case due to induced dipole effects. We find that the benzene anion, known for being unstable in the gas-phase, is stabilized by the plates at zero bias and an inter-plate distance of 21 Å. We also observe the effect of benzene under a voltage bias generated by the plates; under a negative bias, the anion becomes destabilized. We use the electron localization function to analyze the changes in electron density due to the bias. These findings suggest that image effects such as those present in nanoscale devices, are able to stabilize excess charge and should be important to consider when modeling molecular transport junctions and charge-transfer effects.

Two new tools for the acceleration of computational chemistry codes using graphical processing units (GPUs) are presented. First, we propose a general black-box approach for the efficient GPU acceleration of matrix−matrix multiplications where the matrix size is too large for the whole computation to be held in the GPU’s onboard memory. Second, we show how to improve the accuracy of matrix multiplications when using only single-precision GPU devices by proposing a heterogeneous computing model, whereby single- and double-precision operations are evaluated in a mixed fashion on the GPU and central processing unit, respectively. The utility of the library is illustrated for quantum chemistry with application to the acceleration of resolution-of-the-identity second-order Møller−Plesset perturbation theory calculations for molecules, which we were previously unable to treat. In particular, for the 168-atom valinomycin molecule in a cc-pVDZ basis set, we observed speedups of 13.8, 7.8, and 10.1 times for single-, double- and mixed-precision general matrix multiply (SGEMM, DGEMM, and MGEMM), respectively. The corresponding errors in the correlation energy were reduced from −10.0 to −1.2 kcal mol(^{-1}) for SGEMM and MGEMM, respectively, while higher accuracy can be easily achieved with a different choice of cutoff parameter.

Time-dependent density functional theory (TDDFT) has recently been extended to describe many-body open quantum systems evolving under nonunitary dynamics according to a quantum master equation. In the master equation approach, electronic excitation spectra are broadened and shifted due to relaxation and dephasing of the electronic degrees of freedom by the surrounding environment. In this paper, we develop a formulation of TDDFT linear-response theory (LR-TDDFT) for many-body electronic systems evolving under a master equation, yielding broadened excitation spectra. This is done by mapping an interacting open quantum system onto a noninteracting open Kohn–Sham system yielding the correct nonequilibrium density evolution. A pseudoeigenvalue equation analogous to the Casida equations of the usual LR-TDDFT is derived for the Redfield master equation, yielding complex energies and Lamb shifts. As a simple demonstration, we calculate the spectrum of a (C^{2 +}) atom including natural linewidths, by treating the electromagnetic field vacuum as a photon bath. The performance of an adiabatic exchange-correlation kernel is analyzed and a first-order frequency-dependent correction to the bare Kohn–Sham linewidth based on the Görling–Levy perturbation theory is calculated.

Graphical processing units are now being used with dramatic effect to accelerate quantum chemistry calculations. However, early work exposed challenges involving memory bottlenecks and insufficient numerical precision. This research effort addresses those issues, proposing two new tools for accelerating matrix multiplications of arbitrary size where single-precision accuracy is not enough.

Modeling of electronic structure of molecules in electrostatic environments is of considerable relevance for surface-enhanced spectroscopy and molecular electronics. We have developed and implemented a novel approach to the molecular electronic structure in arbitrary electrostatic environments that is compatible with standard quantum chemical methods and can be applied to medium-sized and large molecules. The scheme denoted CheESE (chemistry in electrostatic environments) is based on the description of molecular electronic structure subject to a boundary condition on the system/environment interface. Thus, it is particularly suited to study molecules on metallic surfaces. The proposed model is capable of describing both electrostatic effects near nanostructured metallic surfaces and image-charge effects. We present an implementation of the CheESE model as a library module and show example applications to neutral and negatively charged molecules.