Time-dependent Kohn-Sham electron dynamics coupled with nonequilibrium plasmonic response via atomistic electromagnetic model.

Computational modeling of plasmon-mediated molecular photophysical and photochemical behaviors can help us better understand and tune the bound molecular properties and reactivity and make better decisions to design and control nanostructures. However, computational investigations of coupled plasmon-molecule systems are challenging due to the lack of accurate and efficient protocols to simulate these systems. Here, we present a hybrid scheme by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency dependent fluctuating charge (TD-ωFQ) model. At first, we transform ωFQ in the frequency-domain, an atomistic electromagnetic model for the plasmonic response of plasmonic metal nanoparticles (PMNPs), into the time-domain and derive its equation-of-motion formulation. The TD-ωFQ introduces the nonequilibrium plasmonic response of PMNPs and atomistic interactions to the electronic excitation of the quantum mechanical (QM) region. Then, we combine TD-ωFQ with RT-TDDFT. The derived RT-TDDFT/TD-ωFQ scheme allows us to effectively simulate the plasmon-mediated "real-time" electronic dynamics and even the coupled electron-nuclear dynamics by combining them with the nuclear dynamics approaches. As a first application of the RT-TDDFT/TD-ωFQ method, we study the nonradiative decay rate and plasmon-enhanced absorption spectra of two small molecules in the proximity of sodium MNPs. Thanks to the atomistic nature of the ωFQ model, the edge effect of MNP on absorption enhancement has also been investigated and unveiled.

Real-Time Simulation of Ultrafast Electronic Dynamics of Nanoscale Systems Involving an Organic Molecule and a Nanoparticle Dimer.

Understanding and predicting the behavior of nanomaterials composed of plasmons interacting with quantum emitters at ultrafast timescales is crucial for the better manipulation of light at the nanoscale and advancing technologies like ultrafast communication and computing. Here we perform a simulation of the "real-time" electronic dynamics of a coupled molecule-metal nanoparticle dimer interacting with an ultrashort resonant laser pulse by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency-dependent fluctuating charge (TD-ωFQ) model, an atomistic electromagnetic (AEM) model for the dynamic plasmonic response of nanoparticles. It is shown that the induced dipoles evolve from an exponential decay pattern to a beat pattern with an increase in coupling strength, which is altered by changing the molecular orientation relative to the dimer axis. It is further shown that in the strong coupling regime, both the excited molecule and the plasmon relax rapidly due to the molecule-plasmon interaction, and the efficient coherent energy exchange between the interacting molecule and plasmon modes occurs on a femtosecond (fs) timescale. This work provides guidance on manipulating light-matter interaction and studying molecular plasmonics at extremely fast timescales.

Identification of the interchromophore interaction in the electronic absorption and circular dichroism spectra of bis-phenanthrenes.

We characterize the low-lying excited electronic states of a series of bis-phenanthrenes using our newly developed diabatic scheme called the fragment particle-hole density (FPHD) method and calculate both the electronic absorption and circular dichroism (ECD) spectra using the time-dependent density functional theory (TDDFT) and the FPHD-based exciton model which couples intrachromophore local excitations (LEs) and the interchromophore charge-transfer excitations (CTEs). TDDFT treats each bis-phenanthrene as a single molecule while the mixed LE-CTE exciton model partitions the molecule into two phenanthrene-based aromatic moieties, and then applies the electronic coupling between the various quasi-diabatic states to cover the interactions. It is found that TDDFT and the mixed LE-CTE model reproduce all experimentally observed trends in the spectral profiles, and the hybridization between LE and CTE states is displayed differently in absorption and ECD spectral intensities, as it usually decreases the absorption maxima and affects the positive/negative extrema of the ECD irregularly. By comparing the results yielded by the LE-CTE model with and without the LE-CTE coupling, we identify the contribution of CTE on the main dipole-allowed transitions.

Electronic Couplings for Singlet Fission Processes Based on the Fragment Particle-Hole Densities.

A new diabatization scheme is proposed to calculate the electronic couplings for the singlet fission process in multichromophoric systems. In this approach, a robust descriptor that treats single and multiple excitations on an equal footing is adopted to quantify the localization degree of the particle and hole densities of the electronic states. By maximally localizing the particles and holes in terms of predefined molecular fragments, quasi-diabatic states with well-defined characters (locally excited, charge transfer, correlated triplet pair, etc.) can be automatically constructed as the linear combinations of the adiabatic ones, and the electronic couplings can be directly obtained. This approach is very general in that it applies to electronic states with various spin multiplicities and can be combined with various kinds of preliminary electronic structure calculations. Due to the high numerical efficiency, it is able to manipulate more than 100 electronic states in diabatization. The applications to the tetracene dimer and trimer reveal that high-lying multiply excited charge transfer states have significant influences on both the formation and separation of the correlated triplet pair and can even enlarge the coupling for the latter process by 1 order of magnitude.

Analytical derivative couplings within the framework of time-dependent density functional theory coupled with conductor-like polarizable continuum model: Formalism, implementation, and applications.

The nonadiabatic phenomena, which are characterized by a strong coupling between electronic and nuclear motions, are ubiquitous. The nonadiabatic effect of the studied system can be significantly affected by the surrounding environment, such as solvents, in which such nonadiabatic process takes place. It is essential to develop the theoretical models to simulate these processes while accurately modeling the solvent environment. The time-dependent density functional theory (TDDFT) is currently the most efficient approach to describe the electronic structures and dynamics of complex systems, while the polarizable continuum model (PCM) represents one of the most successful examples among continuum solvation models. Here, we formulate the first-order derivative couplings (DCs) between the ground and excited states as well as between two excited states by utilizing time-independent equation of motion formalism within the framework of both linear response and spin flip formulations of TDDFT/CPCM (the conductor-like PCM), and implement the analytical DCs into the Q-CHEM electronic structure software package. The analytic implementation is validated by the comparison of the analytical and finite-difference results, and reproducing geometric phase effect in the protonated formaldimine test case. Taking 4-(N,N-dimethylamino)benzonitrile and uracil in the gas phase and solution as an example, we demonstrate that the solvent effect is essential not only for the excitation energies of the low-lying excited-states but also for the DCs between these states. Finally, we calculate the internal conversion rate of benzophenone in a solvent with DC being used. The current implementation of analytical DCs together with the existing analytical gradient and Hessian of TDDFT/PCM excited states allows one to study the nonadiabatic effects of relatively large systems in solutions with low computational cost.

Understanding the mechanism of plasmon-driven water splitting: hot electron injection and a near field enhancement effect.

Utilizing plasmon-generated hot carriers to drive chemical reactions has currently become an active area of research in solar photocatalysis at the nanoscale. However, the mechanism underlying exact transfer and the generation dynamics of hot carriers, and the strategies used to further improve the quantum efficiency of the photocatalytic reaction still deserve further investigation. In this work, we perform a nonadiabatic excited-state dynamics study to depict the correlation between the reaction rate of plasmon-driven water splitting (PDWS) and the sizes of gold particles, the incident light frequency and intensity, and the near-field spatial distribution. Four model systems, H2O and Au20@H2O separately interacting with the laser field and the near field generated by the Au nanoparticle (NP) with a few nanometers in size, have been investigated. Our simulated results clearly unveil the mechanism of PDWS and hot-electron injection in a Schottky-free junction: the electrons populated on the antibonding orbitals of H2O are mandatory to drive the OH bond breaking and the strong orbital hybridization between Au20 and H2O creates the conditions for direct electron injection. We further find that the linear dependence of the reaction rate and the field amplitude only holds at a relatively weak field and it breaks down when the second OH bond begins to dissociate and field-induced water fragmentation occurs at a very intensive field, and that with the guarantee of electron injection, the water splitting rate increases with an increase in the NP size. This study will be helpful for further improving the efficiency of photochemical reactions involving plasmon-generated hot carriers and expanding the applications of hot carriers in a variety of chemical reactions.

Mechanism and Origins of Stereoinduction in an Asymmetric Friedel-Crafts Alkylation Reaction of Chalcone Catalyzed by Chiral N, N'-Dioxide-Sc(III) Complex.

The mechanism and selectivity of the asymmetric Friedel-Crafts (F-C) alkylation reaction between indole and chalcone catalyzed by chiral N, N'-dioxide-Sc(III) complexes were investigated at the M06/6-311+G(d,p)//M06/[LANL2DZ,6-31G(d)](SMD,CH2Cl2) level. The reaction occurred via a three-step mechanism: (i) the C3-Cß bond formation by interacting the most mucleophilic C3 center of indole with the most electrophilic Cß center of chalcone; (ii) the abstraction of the proton at the C3 atom of indole by counterion OTf-; (iii) proton transfer from HOTf to the Cα atom of chalcone, generating the F-C alkylation product. The reaction preferred to occur along the favorable re-face attack pathway, producing the dominant R-product. The turnover frequency (TOF) of catalysis was predicted to be 1.59 × 10-7 s-1, with a rate constant of K( T) = 1.58 × 10-7 exp(-29057/ RT) dm6·mol-2·s-1 over the temperature range of 248-368 K. Activation strain model (ASM) and energy decomposition analysis (EDA), as well as noncovalent interaction (NCI) analysis, for the stereocontrolling transition state revealed that the substituent attached to the N atom of the amide subunits as well as the amino acid backbone of ligand played important roles in chiral inductivity. The benzyl group with structural flexibility tended to form strong π-π stacking with substrate as well as the terminal phenyl group of chalcone, stabilizing re-face attack transition state.