Efficient D-π-π-A-Type Dye Sensitizer Based on a Benzothiadiazole Moiety: A Computational Study.

The design of highly efficient sensitizers is one of the most significant areas in dye-sensitized solar cell (DSSC) research. We studied a series of benzothiadiazole-based D-π-π-A organic dyes, putting emphasis on the influence of the donor moiety on the DSSC's efficiency. Using (linear-response time-dependent) density functional theory ((TD)DFT)) with the CAM-B3LYP functional, different donor groups were characterized in terms of electronic absorption spectra and key photovoltaic parameters. As a reference, a dye was considered that had a benzothiadiazole fragment linked via thiophene rings to a diphenylamine donor and a cyanoacrylic-acid acceptor. The different systems were first studied in terms of individual performance parameters, which eventually aggregated into power conversion efficiency. Only the amino-substituted species showed a modest increase, whereas the dimethylamino case showed a decrease.

Site-Selective Real-Time Observation of Bimolecular Electron Transfer in a Photocatalytic System Using L-Edge X-Ray Absorption Spectroscopy*.

Time-resolved X-ray absorption spectroscopy has been utilized to monitor the bimolecular electron transfer in a photocatalytic water splitting system. This has been possible by uniting the local probe and element specific character of X-ray transitions with insights from high-level ab initio calculations. The specific target has been a heteroleptic [IrIII (ppy)2 (bpy)]+ photosensitizer, in combination with triethylamine as a sacrificial reductant and Fe3(CO)12 as a water reduction catalyst. The relevant molecular transitions have been characterized via high-resolution Ir L-edge X-ray absorption spectroscopy on the picosecond time scale and restricted active space self-consistent field calculations. The presented methods and results will enhance our understanding of functionally relevant bimolecular electron transfer reactions and thus will pave the road to rational optimization of photocatalytic performance.

Strong Exciton-Vibrational Coupling in Molecular Assemblies. Dynamics Using the Polaron Transformation in HEOM Space.

In Frenkel exciton dynamics of aggregated molecules, the polaron transformation (PT) technique leads to decoupling of diagonal elements in the subspace of excited electronic states from vibrations. In this article we describe for the first time how PT becomes applicable in the framework of the "Hierarchical Equations of Motion" (HEOM) approach for treatment of open quantum systems. We extend the concept of formulating operators in HEOM space by deriving hierarchical equations of PT which lead to a shift in the excited state potential energy surface to compensate its displacement. While the assumption of thermal equilibration of the vibrational oscillators, introduced by PT, results in a stationary state in a monomer, in a dimer under the same assumption nonequilibrium dynamics appears because of the interplay of the transfer process and vibrational equilibration. Both vertical transitions generating a vibrationally hot state and initially equilibrated vibrational oscillators evolve toward the same stationary asymptotic state associated with polaron formation. The effect of PT on the dynamics of this process depends on initial excitation and basis representation of the electronic system. The developed approach facilitates a generic formulation of quantum master equations involving perturbative treatment of polaron dynamics.

Intraband dynamics and exciton trapping in the LH2 complex of Rhodopseudomonas acidophila.

Over the last several decades, the light-harvesting protein complexes of purple bacteria have been among the most popular model systems for energy transport in excitonic systems in the weak and intermediate intermolecular coupling regime. Despite this extensive body of scientific work, significant questions regarding the excitonic states and the photo-induced dynamics remain. Here, we address the low-temperature electronic structure and excitation dynamics in the light-harvesting complex 2 of Rhodopseudomonas acidophila by two-dimensional electronic spectroscopy. We find that, although at cryogenic temperature energy relaxation is very rapid, exciton mobility is limited over a significant range of excitation energies. This points to the presence of a sub-200 fs, spatially local energy-relaxation mechanism and suggests that local trapping might contribute substantially more in cryogenic experiments than under physiological conditions where the thermal energy is comparable to or larger than the static disorder.

Ground- and Excited-State Properties of Iron(II) Complexes Linked to Organic Chromophores.

Two new bichromophoric complexes, [Fe(bim-ant)2]2+ and [Fe(bim-pyr)2]2+ ([H2-bim]2+ = 1,1'-(pyridine-2,6-diyl)bis(3-methyl-1H-imidazol-3-ium); ant = 9-anthracenyl; pyr = 1-pyrenyl), are investigated to explore the possibility of tuning the excited-state behavior in photoactive iron(II) complexes to design substitutes for noble-metal compounds. The ground-state properties of both complexes are characterized thoroughly by electrochemical methods and optical absorption spectroscopy, complemented by time-dependent density functional theory calculations. The excited states are investigated by static and time-resolved luminescence and femtosecond transient absorption spectroscopy. Both complexes exhibit room temperature luminescence, which originates from singlet states dominated by the chromophore (1Chrom). In the cationic pro-ligands and in the iron(II) complexes, the emission is shifted to red by up to 110 nm (5780 cm-1). This offers the possibility of tuning the organic chromophore emission by metal-ion coordination. The fluorescence lifetimes of the complexes are in the nanosecond range, while triplet metal-to-ligand charge-transfer (3MLCT) lifetimes are around 14 ps. An antenna effect as in ruthenium(II) polypyridine complexes connected to an organic chromophore is found in the form of an internal conversion within 3.4 ns from the 1Chrom to the 1MLCT states. Because no singlet oxygen forms from triplet oxygen in the presence of the iron(II) complexes and light, efficient intersystem crossing to the triplet state of the organic chromophore (3Chrom) is not promoted in the iron(II) complexes.

Molecular level picture of the interplay between pH and phosphate binding at the goethite-water interface.

The soil pH plays a substantial role in controlling phosphorus (P) adsorption and mobilization. These processes are strongly affected by the phosphate interaction strength with P-fixing soil minerals such as goethite. The target of the current contribution is to draw a molecular level picture of the interplay between pH and phosphate binding at the goethite-water interface via a joint experimental-theoretical approach. Periodic density functional theory (DFT) calculations were carried out to provide a molecular level understanding of the pH dependence of P adsorption. To validate the modeling approach, adsorption experiments of phosphate at goethite were performed in the pH range of 4-12. There was agreement between experiments and simulations in the description of the adsorption behavior by two pH-dependent successive stages. The adsorption increases along the pH change from 4 to 8. A further increase of pH leads to a decrease of adsorption. By comparing with literature data it is concluded that the first stage will be observed only if there is no significant change of the surface charge at low pH. Moreover, the molecular modeling results point to the abundance of the monodentate (M) binding motif at both extremely low and high pH ranges. Otherwise, the bidentate (B) one is predominant along the intermediate pH range. These observations could resolve the existing debate about the assignment of phosphate-goethite binding motifs. Furthermore, the results point to a decrease of pH upon phosphate sorption due to an induced acidification of soil solution. The present joint experimental-theoretical approach provides a better understanding and description of the existing phosphate sorption experiments and highlights new findings at the atomistic/molecular scale.

The effect of N-heterocyclic carbene units on the absorption spectra of Fe(II) complexes: a challenge for theory.

The absorption spectra of five Fe(ii) homoleptic and heteroleptic complexes containing strong sigma-donating N-heterocyclic carbene (NHC) and polypyridyl ligands have been theoretically characterized using a tuned range-separation functional. From a benchmark comparison of the obtained results against other functionals and a multiconfigurational reference, it is concluded that none of the methods is completely satisfactory to describe the absorption spectra. As a compromise using 20% exact exchange, the electronic excited states underlying the absorption spectra are analyzed. The low-lying energy band of all the compounds shows predominant metal-to-ligand charge transfer (MLCT) character while the triplet excited states have metal-centered (MC) nature, which becomes more pronounced with increasing the number of NHC-donor groups. Excited MC states with partial charge transfer to the NHC-donor groups are higher in energy than comparable states without these contributions. The presence of the low-lying MC states prevents the formation of long-lived MLCT states.

Exciton transfer using rates extracted from the "hierarchical equations of motion".

Frenkel exciton population dynamics of an excitonic dimer is studied by comparing the results from a quantum master equation involving rates from second-order perturbative treatment with respect to the excitonic coupling with the non-perturbative results from "Hierarchical Equations of Motion" (HEOM). By formulating generic Liouville-space expressions for the rates, we can choose to evaluate them either via HEOM propagations or by applying the cumulant expansion. The coupling of electronic transitions to bath modes is modeled either as overdamped oscillators for the description of thermal bath components or as underdamped oscillators to account for intramolecular vibrations. Cases of initial nonequilibrium and equilibrium vibrations are discussed. In the case of HEOM, initial equilibration enters via a polaron transformation. Pointing out the differences between the nonequilibrium and equilibrium approach in the context of the projection operator formalism, we identify a further description, where the transfer dynamics is driven only by fluctuations without involvement of dissipation. Despite this approximation, this approach can also yield meaningful results in certain parameter regimes. While for the chosen model, HEOM has no technical advantage for evaluation of the rate expressions compared to cumulant expansion, there are situations where only evaluation with HEOM is applicable. For instance, a separation of reference and interaction Hamiltonian via a polaron transformation to account for the interplay between Coulomb coupling and vibrational oscillations of the bath at the level of a second-order treatment can be adjusted for a treatment with HEOM.

Ab Initio Molecular Dynamics Simulations of the Interaction between Organic Phosphates and Goethite.

Today's fertilizers rely heavily on mining phosphorus (P) rocks. These rocks are known to become exhausted in near future, and therefore effective P use is crucial to avoid food shortage. A substantial amount of P from fertilizers gets adsorbed onto soil minerals to become unavailable to plants. Understanding P interaction with these minerals would help efforts that improve P efficiency. To this end, we performed a molecular level analysis of the interaction of common organic P compounds (glycerolphosphate (GP) and inositol hexaphosphate (IHP)) with the abundant soil mineral (goethite) in presence of water. Molecular dynamics simulations are performed for goethite-IHP/GP-water complexes using the multiscale quantum mechanics/molecular mechanics method. Results show that GP forms monodentate (M) and bidentate mononuclear (B) motifs with B being more stable than M. IHP interacts through multiple phosphate groups with the 3M motif being most stable. The order of goethite-IHP/GP interaction energies is GP M < GP B < IHP M < IHP 3M. Water is important in these interactions as multiple proton transfers occur and hydrogen bonds are formed between goethite-IHP/GP complexes and water. We also present theoretically calculated infrared spectra which match reasonably well with frequencies reported in literature.

Infrared Laser Excitation Controlled Reaction Acceleration in the Electronic Ground State.

Propelling a ground state reaction by mode-specific vibrational excitation via infrared (IR) light offers a novel route to carry out ground state chemistry. Here, we describe the acceleration of a bimolecular alcoholysis reaction as a paradigm for IR light-driven ground state reactions. Instead of resorting to coherent control, IR light is used for direct or indirect vibrational excitation of the reaction coordinate (RC) overcoming the activation energy and promoting the ground state reaction with negligible heating of the sample. Thus, knowledge of the RC is crucial to pick the reaction accelerating vibrations. Alternatively, upon mapping the reaction accelerating vibrations an image of the RC can be reconstructed. We discuss the concept of RCs and examine strategies to use vibrational energy relaxation pathways to single out vibrations belonging to the RC. The influence of the solvent interaction and limitations due to conformational heterogeneity are considered. We provide an application example generating microstructures of polymers and address the use for chemical synthesis in general.

QM/MM simulations of organic phosphorus adsorption at the diaspore-water interface.

Phosphorus (P) immobilization and thus its availability for plants are mainly affected by the strong interaction of phosphates with soil components especially soil mineral surfaces. The related reactions have been studied extensively via sorption experiments especially by carrying out adsorption of ortho-phosphates onto Fe-oxide surfaces. But a molecular-level understanding of the P-binding mechanisms at the mineral-water interface is still lacking, especially for forest eco-systems. Therefore, the current contribution provides an investigation of the molecular binding mechanisms for two abundant phosphates in forest soils, inositol hexaphosphate (IHP) and glycerolphosphate (GP), at the diaspore mineral surface. Here a hybrid electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) based molecular dynamics simulation has been applied to explore the diaspore-IHP/GP-water interactions. The results provide evidence for the formation of different P-diaspore binding motifs involving monodentate (M) and bidentate (B) for GP and two (2M) as well as three (3M) monodentates for IHP. The interaction energy results indicated the abundance of the GP B motif compared to the M one. The IHP 3M motif has a higher total interaction energy compared to its 2M motif, but exhibits a lower interaction energy per bond. Compared to GP, IHP exhibited stronger interaction with the surface as well as with water. Water was found to play an important role in controlling these diaspore-IHP/GP-water interactions. The interfacial water molecules form moderately strong H-bonds (HBs) with GP and IHP as well as with the diaspore surface. For all the diaspore-IHP/GP-water complexes, the interaction of water with the diaspore exceeds that with the studied phosphates. Furthermore, some water molecules form covalent bonds with diaspore Al atoms while others dissociate at the surface to protons and hydroxyl groups leading to proton transfer processes. Finally, the current results confirm the previous experimental conclusions indicating the importance of the number of phosphate groups, HBs, and proton transfers in controlling the P-binding at soil mineral surfaces.

Infrared spectroscopic characterization of phosphate binding at the goethite-water interface.

The interaction between phosphates and soil mineral surfaces, such as Fe- and Al-(oxyhydr)oxides, plays a crucial role in the immobilization of P and thus its availability for plants. The reactions of phosphates with Fe-hydroxides and especially goethite have been studied extensively. But a molecular-level picture of the phosphate binding mechanisms at the goethite-water interface is still lacking. Therefore, in the current contribution we have explored the molecular binding mechanisms for the adsorbed phosphate at the goethite-water interface by performing sorption kinetics experiments for orthophosphate and characterizing the adsorbed species by FT-IR spectroscopy. In parallel, periodic DFT calculations have been performed to explore the interaction mechanisms and to assign the IR spectra for monodentate (M) and bidentate (B) orthophosphate complexes at two different goethite surface planes (010 and 100) in the presence of water. In general, our interaction energy results give evidence that the mono-protonated B phosphate complex is favored to be formed at the goethite-water interface, although the M motif could exist as a minor fraction. Moreover, it was found that water plays an important role in controlling the phosphate adsorption process at the goethite surfaces. The interfacial water molecules form H-bonds (HBs) with the phosphate as well as with the goethite surface atoms. Furthermore, some water molecules form covalent bonds with goethite's Fe atoms while others dissociate at the surface to protons and hydroxyl groups. The present theoretical assignment of IR spectra introduces a benchmark for characterizing experimental IR data for the adsorbed KH2PO4 species at the goethite-water interface. In particular, the IR spectra of the mono-protonated (2O + 1Fe) B complex at the 010 goethite surface plane and the M complex at the 100 goethite surface plane were found to be consistent with the experimental data. In order to explore the role of different abundances of surface planes and binding motifs, IR spectra obtained from weighted averages have been analyzed. The results confirmed the conclusions drawn from interaction energy calculations.

On computing spectral densities from classical, semiclassical, and quantum simulations.

The Caldeira-Leggett model provides a compact characterization of a thermal environment in terms of a spectral density function, which has led to a variety of numerically exact quantum methods for reduced density matrix propagation. Since spectral densities are often computed from classical molecular dynamics simulations, we investigate in this paper whether quantum effects should be accounted for in the calculations. Therefore, we reformulate the recently developed Fourier method for spectral density calculations from semiclassical simulations which approximately allow for quantum effects. We propose two possible protocols based on either correlation functions or expectation values. These protocols are tested on a generic Calderra-Leggett model for the linearized semiclassical initial-value representation (LSC-IVR), the thawed Gaussian wave packet dynamics (TGWD), and hybrid schemes combining the two with the more accurate Herman-Kluk formula. Surprisingly, spectral densities from the LSC-IVR method, which treats the dynamics completely classically, are found to be extremely accurate, even in the quantum regime, where this method does not give a correct description of the correlation functions and expectation values. In contrast, the TGWD method turns out as too inaccurate for spectral density calculations, and the hybrid schemes perform well only if the system is close to the classical regime. This implies that, if the bath has a Caldeira-Leggett form, spectral densities are insensitive to quantum effects and any effort to approximately account for them rather leads to errors. Hence, in this case, spectral densities can be computed from classical simulations and used in a reduced quantum simulation as well.

Influence of metal ions on glyphosate detection by FMOC-Cl.

Glyphosate (GLP, N-(phosphonomethyl)glycine) is the most important broadband herbicide in the world, but discussions are controversial regarding its environmental behaviour and distribution. Residue analyses in a variety of environmental samples are commonly conducted by HPLC-MS where GLP needs to be derivatised with 9-fluoromethoxycarnonyl chloride (FMOC-Cl). Since this derivatisation reaction was suspected to be inhibited by metal ions in the sample matrix, the present study provides a comprehensive experimental study of the effect of metal ions (Al3+, Ca2+, Cd2+, Co2+, Cu2+, Fe2+, Fe3+, Mg2+, Mn2+, Zn2+) on derivatisation and GLP recovery. Results show that some metals (Cd2+, Co2+, Cu2+, Mn2+ and Zn2+) decreased the GLP recovery down to 19 to 59%. Complementary, quantum chemical modelling of 1:1 GLP-metal complexes as well as their reactivity with respect to FMOC-Cl was performed. Here, a decrease in reactivity of FMOC-Cl towards GLP-metal complexes is observed; i.e. the reaction is non-spontaneous in contrast to the free GLP case. The present results are in accord with previous studies and provide an explanation that full GLP recovery in different matrices was never reached. Remedy strategies to compensate for the inhibition effect are explored such as pH adjustment to acidic or alkaline conditions or addition of ethylenediaminetetraacetic acid (EDTA). In general, our results question the use of internal isotopic labelled standards (ILS) since this presupposes the presence of the analyte and the ILS in the same (free) form.

Erratum: Ultrafast Spin-State Dynamics in Transition-Metal Complexes Triggered by Soft-X-Ray Light [Phys. Rev. Lett. 118, 023001 (2017)].

This corrects the article DOI: 10.1103/PhysRevLett.118.023001.

Unravelling the nature of glyphosate binding to goethite surfaces by ab initio molecular dynamics simulations.

Investigation of the interaction between glyphosate (GLP) and soil minerals is essential for understanding GLP's fate in the environment. Whereas GLP-goethite binding has been discussed extensively, the impact of water as well as of different goethite surface planes has not been studied yet. In this contribution, periodic density functional theory-based molecular dynamics simulations are applied to explore possible binding mechanisms for GLP with three goethite surface planes (010, 001, and 100) in the presence of water. The investigation included several binding motifs of monodentate (M) and bidentate (B) type. It was found that the binding stability increases in the order M@001 < M@010 < (2O + 2Fe) B@100 < M@100 < (1O + 2Fe) B@001 < (2O + 1Fe) B@010. This behavior has been traced to the presence of intramolecular H-bonds (HBs) in GLP as well as intermolecular HBs between GLP and water, GLP and goethite, and water and goethite. These interactions are accompanied by proton transfer from GLP to water and to goethite, and from water to goethite as well as water dissociation at the goethite surface. Further, it was observed that the OH- species can replace the adsorbed GLP at the goethite surface, which could explain the well-known drastic drop in GLP adsorption at high pH. The present results highlight the role of water in the GLP-goethite interaction and provide a molecular level perspective on available experimental data.

Photoelectron shake-ups as a probe of molecular symmetry: 4d XPS analysis of I3- in solution.

A combination of multi-configurational restricted active space calculations with a Dyson orbital formalism has been applied for accurate simulations of 4d photo-electron spectra of the I3- molecular ion. The analysis based on the occupation numbers of natural orbitals allowed to predict and rationalize the spectral fingerprints of solvent-induced nuclear asymmetry. In particular, it demonstrates how the nuclear asymmetry directly causes an increase of shake-up intensity. The relative intensity of shake-up and main features of the I 4d XPS spectrum could therefore serve as a simplified experimental observable of structural asymmetry, complementary to changes in the shape of the main spectral features.

Using fluorescence detected two-dimensional spectroscopy to investigate initial exciton delocalization between coupled chromophores.

Förster theory describes electronic exciton energy migration in molecular assemblies as an incoherent hopping process between donor and acceptor molecules. The rate is expressed in terms of the overlap integral between donor fluorescence and acceptor absorption spectra. Typical time scales for systems like photosynthetic antennae are on the order of a few picoseconds. Prior to transfer, it is assumed that the initially excited donor molecule has equilibrated with respect to the local environment. However, upon excitation and during the equilibration phase, the state of the system needs to be described by the full density matrix, including coherences between donor and acceptor states. While being intuitively clear, addressing this regime experimentally has been a challenge until the recently reported advances in fluorescence detected two-dimensional spectroscopy. Here, we demonstrate using fourth order perturbation theory the conditions for the presence of donor-acceptor coherence induced cross-peaks at zero waiting time between the first and the second pair of pulses. The approach is illustrated for a heterodimer model which facilitates an analytical solution.

Simulating vibronic spectra via Matsubara-like dynamics: Coping with the sign problem.

Measuring vibronic spectra probes dynamical processes in molecular systems. When interpreted via suitable theoretical tools, the experimental data provides comprehensive information about the system in question. For complex many-body problems, such an approach usually requires the formulation of proper classical-like approximations, which is particularly challenging if multiple electronic states are involved. In this work, we express the imaginary-time shifted time correlation function and, thus, the vibronic spectrum in terms of the so-called Matsubara dynamics, which combines quantum statistics and classical-like dynamics. By applying the Matsubara approximation in the adiabatic limit, we derive a formal generalization of the existing Matsubara dynamics formalism to multiple potential energy surfaces (PESs), which, however, does not feature all the defining properties of its single-PES counterpart though suffering equally from the sign problem. The mathematical analysis for two shifted harmonic oscillators suggests a new modified method to practically simulate the standard correlation function via Matsubara-like dynamics. Importantly, this modified method samples the thermal Wigner function without suffering from the sign problem and yields an accurate approximation to the vibronic absorption spectrum, not only for the harmonic system but also for the anharmonic one.

Quasi-classical approaches to vibronic spectra revisited.

The framework to approach quasi-classical dynamics in the electronic ground state is well established and is based on the Kubo-transformed time correlation function (TCF), being the most classical-like quantum TCF. Here we discuss whether the choice of the Kubo-transformed TCF as a starting point for simulating vibronic spectra is as unambiguous as it is for vibrational ones. Employing imaginary-time path integral techniques in combination with the interaction representation allowed us to formulate a method for simulating vibronic spectra in the adiabatic regime that takes nuclear quantum effects and dynamics on multiple potential energy surfaces into account. Further, a generalized quantum TCF is proposed that contains many well-established TCFs, including the Kubo one, as particular cases. Importantly, it also provides a framework to construct new quantum TCFs. Applying the developed methodology to the generalized TCF leads to a plethora of simulation protocols, which are based on the well-known TCFs as well as on new ones. Their performance is investigated on 1D anharmonic model systems at finite temperatures. It is shown that the protocols based on the new TCFs may lead to superior results with respect to those based on the common ones. The strategies to find the optimal approach are discussed.