Convergence of Computed Aqueous Absorption Spectra with Explicit Quantum Mechanical Solvent.

For reliable condensed phase simulations, an accurate model that includes both short- and long-range interactions is required. Short- and long-range interactions can be particularly strong in aqueous solution, where hydrogen-bonding may play a large role at short range and polarization may play a large role at long range. Although short-range solute-solvent interactions such as charge transfer, hydrogen bonding, and solute-solvent polarization can be taken into account with a quantum mechanical (QM) treatment of the solvent, it is unclear how much QM solvent is necessary to accurately model interactions with different solutes. In this work, we investigate the effect of explicit QM solvent on absorption spectra computed for a series of solutes with decreasing polarity. By adjusting the boundary between QM and classical molecular mechanical solvent to include up to 400 QM water molecules, convergence of the calculated absorption spectra with respect to the size of the QM region is achieved. We find that the rate of convergence does not correlate with solute polarity when excitation energies are calculated using time-dependent density functional theory with a range-separated hybrid functional, but it does correlate with solute polarity when using configuration interaction singles. We also find that larger basis sets converge the computed spectrum with fewer QM solvent molecules. To optimize the computational cost with respect to convergence, we test a mixed basis set with more basis functions for atoms of the chromophore and the solvent molecules that are nearest to it and fewer basis functions for the atoms of the remaining solvent molecules in the QM region. Our results show that using a mixed basis set is a potentially effective way to significantly lower the computational cost while reproducing the results computed with larger basis sets.

Convergence of Excitation Energies in Mixed Quantum and Classical Solvent: Comparison of Continuum and Point Charge Models.

Mixed quantum mechanical (QM)/classical methods provide a computationally efficient approach to modeling both ground and excited states in the condensed phase. To accurately model short-range interactions, some amount of the environment can be included in the QM region, whereas a classical model can treat long-range interactions to maintain computational affordability. The best computational protocol for these mixed QM/classical methods can be determined by examining convergence of molecular properties. Here, we compare molecular mechanical (MM) fixed point charges to a polarizable continuum model (PCM) for computing electronic excitations in solution. We computed the excitation energy of three pairs of neutral/anionic molecules in aqueous solvent, including up to 250 water molecules in the QM region. Interestingly, the convergence is similar for MM point charges and a PCM, with convergence achieved when at least one full solvation shell is treated with QM. Although the van der Waals (VDW) definition of the PCM cavity is adequate with small amounts of QM solvent, larger QM solvent layers had gaps in the VDW PCM cavity, leading to asymptotically incorrect excitation energies. Given that the VDW cavity leads to unphysical solute-solvent interactions, we advise using a solvent-excluded surface cavity for QM/PCM calculations that include QM solvent.

Multistate Density Functional Theory for Effective Diabatic Electronic Coupling.

Multistate density functional theory (MSDFT) is presented to estimate the effective transfer integral associated with electron and hole transfer reactions. In this approach, the charge-localized diabatic states are defined by block localization of Kohn-Sham orbitals, which constrain the electron density for each diabatic state in orbital space. This differs from the procedure used in constrained density functional theory that partitions the density within specific spatial regions. For a series of model systems, the computed transfer integrals are consistent with experimental data and show the expected exponential attenuation with the donor-acceptor separation. The present method can be used to model charge transfer reactions including processes involving coupled electron and proton transfer.

Peak-shifting in real-time time-dependent density functional theory.

In recent years, the development and application of real-time time-dependent density functional theory (RT-TDDFT) has gained momentum as a computationally efficient method for modeling electron dynamics and properties that require going beyond a linear response of the electron density. However, the RT-TDDFT method within the adiabatic approximation can unphysically shift absorption peaks throughout the electron dynamics. Here, we investigate the origin of these time-dependent resonances observed in RT-TDDFT spectra. Using both exact exchange and hybrid exchange-correlation approximate functionals, adiabatic RT-TDDFT gives time-dependent absorption spectra in which the peaks shift in energy as populations of the excited states fluctuate, while exact wave function methods yield peaks that are constant in energy but vary in intensity. The magnitude of the RT-TDDFT peak shift depends on the frequency and intensity of the applied field, in line with previous studies, but it oscillates as a function of time-dependent molecular orbital populations, consistent with a time-dependent superposition electron density. For the first time, we provide a rationale for the direction and magnitude of the time-dependent peak shifts based on the molecular electronic structure. For three small molecules, H2, HeH(+), and LiH, we give contrasting examples of peak-shifting to both higher and lower energies. The shifting is explained as coupled one-electron transitions to a higher and a lower lying state. Whether the peak shifts to higher or lower energies depends on the relative energetics of these one-electron transitions.

Two-electron Rabi oscillations in real-time time-dependent density-functional theory.

We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S0 state and the doubly-excited S2 state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation.

Explicit polarization: a quantum mechanical framework for developing next generation force fields.

Conspectus Molecular mechanical force fields have been successfully used to model condensed-phase and biological systems for a half century. By means of careful parametrization, such classical force fields can be used to provide useful interpretations of experimental findings and predictions of certain properties. Yet, there is a need to further improve computational accuracy for the quantitative prediction of biomolecular interactions and to model properties that depend on the wave functions and not just the energy terms. A new strategy called explicit polarization (X-Pol) has been developed to construct the potential energy surface and wave functions for macromolecular and liquid-phase simulations on the basis of quantum mechanics rather than only using quantum mechanical results to fit analytic force fields. In this spirit, this approach is called a quantum mechanical force field (QMFF). X-Pol is a general fragment method for electronic structure calculations based on the partition of a condensed-phase or macromolecular system into subsystems ("fragments") to achieve computational efficiency. Here, intrafragment energy and the mutual electronic polarization of interfragment interactions are treated explicitly using quantum mechanics. X-Pol can be used as a general, multilevel electronic structure model for macromolecular systems, and it can also serve as a new-generation force field. As a quantum chemical model, a variational many-body (VMB) expansion approach is used to systematically improve interfragment interactions, including exchange repulsion, charge delocalization, dispersion, and other correlation energies. As a quantum mechanical force field, these energy terms are approximated by empirical functions in the spirit of conventional molecular mechanics. This Account first reviews the formulation of X-Pol, in the full variationally correct version, in the faster embedded version, and with systematic many-body improvements. We discuss illustrative examples involving water clusters (which show the power of two-body corrections), ethylmethylimidazolium acetate ionic liquids (which reveal that the amount of charge transfer between anion and cation is much smaller than what has been assumed in some classical simulations), and a solvated protein in aqueous solution (which shows that the average charge distribution of carbonyl groups along the polypeptide chain depends strongly on their position in the sequence, whereas they are fixed in most classical force fields). The development of QMFFs also offers an opportunity to extend the accuracy of biochemical simulations to areas where classical force fields are often insufficient, especially in the areas of spectroscopy, reactivity, and enzyme catalysis.

The quantum coherent mechanism for singlet fission: experiment and theory.

The absorption of one photon by a semiconductor material usually creates one electron-hole pair. However, this general rule breaks down in a few organic semiconductors, such as pentacene and tetracene, where one photon absorption may result in two electron-hole pairs. This process, where a singlet exciton transforms to two triplet excitons, can have quantum yields as high as 200%. Singlet fission may be useful to solar cell technologies to increase the power conversion efficiency beyond the so-called Shockley-Queisser limit. Through time-resolved two-photon photoemission (TR-2PPE) spectroscopy in crystalline pentacene and tetracene, our lab has recently provided the first spectroscopic signatures in singlet fission of a critical intermediate known as the multiexciton state (also called a correlated triplet pair). More importantly, we found that population of the multiexciton state rises at the same time as the singlet state on the ultrafast time scale upon photoexcitation. This observation does not fit with the traditional view of singlet fission involving the incoherent conversion of a singlet to a triplet pair. However, it provides an experimental foundation for a quantum coherent mechanism in which the electronic coupling creates a quantum superposition of the singlet and the multiexciton state immediately after optical excitation. In this Account, we review key experimental findings from TR-2PPE experiments and present a theoretical analysis of the quantum coherent mechanism based on electronic structural and density matrix calculations for crystalline tetracene lattices. Using multistate density functional theory, we find that the direct electronic coupling between singlet and multiexciton states is too weak to explain the experimental observation. Instead, indirect coupling via charge transfer intermediate states is two orders of magnitude stronger, and dominates the dynamics for ultrafast multiexciton formation. Density matrix calculation for the crystalline tetracene lattice satisfactorily accounts for the experimental observations. It also reveals the critical roles of the charge transfer states and the high dephasing rates in ensuring the ultrafast formation of multiexciton states. In addition, we address the origins of microscopic relaxation and dephasing rates, and adopt these rates in a quantum master equation description. We show the need to take the theoretical effort one step further in the near future by combining high-level electronic structure calculations with accurate quantum relaxation dynamics for large systems.

The Third Dimension of a More O'Ferrall-Jencks Diagram for Hydrogen Atom Transfer in the Isoelectronic Hydrogen Exchange Reactions of (PhX)(2)H(•) with X = O, NH, and CH(2).

A critical element in theoretical characterization of the mechanism of proton-coupled electron transfer (PCET) reactions, including hydrogen atom transfer (HAT), is the formulation of the electron and proton localized diabatic states, based on which a More O'Ferrall-Jencks diagram can be represented to determine the step-wise and concerted nature of the reaction. Although the More O'Ferrall-Jencks diabatic states have often been used empirically to develop theoretical models for PCET reactions, the potential energy surfaces for these states have never been determined directly based on first principles calculations using electronic structure theory. The difficulty is due to a lack of practical method to constrain electron and proton localized diabatic states in wave function or density functional theory calculations. Employing a multistate density functional theory (MSDFT), in which the electron and proton localized diabatic configurations are constructed through block-localization of Kohn-Sham orbitals, we show that distinction between concerted proton-electron transfer (CPET) and HAT, which are not distinguishable experimentally from phenomenological kinetic data, can be made by examining the third dimension of a More O'Ferrall-Jencks diagram that includes both the ground and excited state potential surfaces. In addition, we formulate a pair of effective two-state valence bond models to represent the CPET and HAT mechanisms. We found that the lower energy of the CPET and HAT effective diabatic states at the intersection point can be used as an energetic criterion to distinguish the two mechanisms. In the isoelectronic series of hydrogen exchange reaction in (PhX)(2)H(•), where X = O, NH, and CH(2), there is a continuous transition from a CPET mechanism for the phenoxy radical-phenol pair to a HAT process for benzyl radical and toluene, while the reaction between PhNH(2) and PhNH(•) has a mechanism intermediate of CPET and HAT. The electronically nonadiabatic nature of the CPET mechanism in the phenol system can be attributed to the overlap interactions between the ground and excited state surfaces, resulting in roughly orthogonal minimum energy paths on the adiabatic ground and excited state potential energy surfaces. On the other hand, the minimum energy path on the adiabatic ground state for the HAT mechanism coincides with that on the excited state, producing a large electronic coupling that separates the two surfaces by more than 120 kcal/mol.

Formyloxyl radical-gold nanoparticle binding: a theoretical study.

The citrate reduction method is one of the simplest and most common methods used in the synthesis of gold nanoparticles. It has been thought that citrate acts as both a reducing agent for the gold salt and as the capping agent. However, it has recently been reported using density functional theory (DFT) that electron density builds up on uncomplexed apex gold atoms and the binding of formate (the simplest carboxylate and a model for citrate) becomes unfavorable after two additions, limiting citrate's utility as a capping agent. In this study, Au(20)-formyloxyl radical interactions are investigated using DFT at the BP86/DZ level of theory to model neutral carboxylate-gold nanoparticle binding (corresponding to carboxylates interacting with a partially oxidized gold nanoparticle). Binding energies are refined using a TZP basis set. It is found that the incremental binding energies of formyloxyl radicals remain highly favorable through eight additions (the highest number tested). The addition of one formyloxyl radical is 56 kJ/mol less than the addition of one formate but becomes 210 kJ/mol more favorable for the second addition. The range of binding energies through the eight additions is 154-331 kJ/mol. Furthermore, after the third addition, the most favorable geometries feature distortion of the gold tetrahedron. These results suggest that oxidized species formed in the citrate reduction method are likely capping agents and that binding of these ligands may affect the properties of the nanoparticles through distortion of the gold structure.

Origin of intense chiroptical effects in undecagold subnanometer particles.

Time-dependent density functional theory (TDDFT) calculations are employed to examine the optical absorption and circular dichroism (CD) spectra of undecagold Au(11)L(4)X(2)(+) (X = Cl, Br) complexes and their Au(2)X(2)L precursors, where L is either 2,2'-bis(diphenylphosphino)-1,1'-binaphtyl (BINAP) or 1,4-diphosphino-1,3-butadiene (dpb). These systems exhibit intense and mirror-image Cotton effects in their CD spectra. Experimental peak positions are well reproduced in the calculations. The low energy peaks of Au(11)L(4)X(2)(+) arise primarily from transitions between delocalized metal superatom orbitals. Bidentate phosphine ligands have both a structural and electronic impact on the system. The lowest energy structure of Au(11)L(4)X(2)(+) has a chiral C(2) geometry, whereas monodentate phosphine ligands lead to a C(1) structure. In addition, the chiral core structure of Au(11)L(4)X(2)(+) is not sufficient to explain the strong Cotton effects, and the intensity of the CD spectrum is increased by the presence of the bidentate phosphine ligands.