Visible-Light-Activated Carbon Monoxide Release from Porphyrin-Flavonol Hybrids.

We report on porphyrin-flavonol hybrids consisting of a porphyrin antenna and four covalently bound 3-hydroxyflavone (flavonol) groups, which act as highly efficient photoactivatable carbon monoxide (CO)-releasing molecules (photoCORMs). These bichromophoric systems enable activation of the UV-absorbing flavonol chromophore by visible light up to 650 nm and offer precise spatial and temporal control of CO administration. The physicochemical properties of the porphyrin antenna system can also be tuned by inserting a metal cation. Our computational study revealed that the process occurs via endergonic triplet-triplet energy transfer from porphyrin to flavonol and may become feasible thanks to flavonol energy stabilization upon intramolecular proton transfer. This mechanism was also indirectly supported by steady-state and transient absorption spectroscopy techniques. Additionally, the porphyrin-flavonol hybrids were found to be biologically benign. With four flavonol CO donors attached to a single porphyrin chromophore, high CO release yields, excellent uncaging cross sections, low toxicity, and CO therapeutic properties, these photoCORMs offer exceptional potential for their further development and future biological and medical applications.

How Does Mg2+(aq) Interact with ATP(aq)? Biomolecular Structure through the Lens of Liquid-Jet Photoemission Spectroscopy.

Liquid-jet photoemission spectroscopy (LJ-PES) allows for a direct probing of electronic structure in aqueous solutions. We show the applicability of the approach to biomolecules in a complex environment, exploring site-specific information on the interaction of adenosine triphosphate in the aqueous phase (ATP(aq)) with magnesium (Mg2+(aq)), highlighting the synergy brought about by the simultaneous analysis of different regions in the photoelectron spectrum. In particular, we demonstrate intermolecular Coulombic decay (ICD) spectroscopy as a new and powerful addition to the arsenal of techniques for biomolecular structure investigation. We apply LJ-PES assisted by electronic-structure calculations to study ATP(aq) solutions with and without dissolved Mg2+. Valence photoelectron data reveal spectral changes in the phosphate and adenine features of ATP(aq) due to interactions with the divalent cation. Chemical shifts in Mg 2p, Mg 2s, P 2p, and P 2s core-level spectra as a function of the Mg2+/ATP concentration ratio are correlated to the formation of [Mg(ATP) 2]6-(aq), [MgATP]2-(aq), and [Mg2ATP](aq) complexes, demonstrating the element sensitivity of the technique to Mg2+-phosphate interactions. The most direct probe of the intermolecular interactions between ATP(aq) and Mg2+(aq) is delivered by the emerging ICD electrons following ionization of Mg 1s electrons. ICD spectra are shown to sensitively probe ligand exchange in the Mg2+-ATP(aq) coordination environment. In addition, we report and compare P 2s data from ATP(aq) and adenosine mono- and diphosphate (AMP(aq) and ADP(aq), respectively) solutions, probing the electronic structure of the phosphate chain and the local environment of individual phosphate units in ATP(aq). Our results provide a comprehensive view of the electronic structure of ATP(aq) and Mg2+-ATP(aq) complexes relevant to phosphorylation and dephosphorylation reactions that are central to bioenergetics in living organisms.

Spin-Vibronic Coupling Controls the Intersystem Crossing of Iodine-Substituted BODIPY Triplet Chromophores.

4,4-Difluoro-4-borata-3a-azonia-4a-aza-s-indacene (BODIPY) dyes are extensively used in various applications of their triplet states, ranging from photoredox catalysis, through triplet sensitization to photodynamic therapy. However, the rational design of BODIPY triplet chromophores by ab initio modelling is limited by their strong interactions of spin, electronic and vibrational dynamics. In particular, spin-vibronic coupling is often overlooked when estimating intersystem crossing (ISC) rates. In this study, a combined experimental and theoretical approach using spin-vibronic coupling to correctly describe ISC in BODIPY dyes was developed. For this purpose, seven π-extended BODIPY derivatives with iodine atoms in different positions were examined. It was found that the heavy-atom effect of iodine atoms is site specific, causing high triplet yields in only some positions. This site-specific ISC was explained by El-Sayed rules, so both the contribution and character of the molecular orbitals involved in the excitation must be considered when predicting the ISC rates. Overall, the rational design of BODIPY triplet chromophores requires using (i) the high-quality electronic structure theory, including both static and dynamical correlations; and (ii) the two-component wave function Hamiltonian, and rationalizing; and (iii) ISC based on the character of the molecular orbitals of heavy atoms involved in the excitation, expanding El-Sayed rules beyond their traditional applications.

Fast and accurate excited states predictions: machine learning and diabatization.

The efficiency of machine learning algorithms for electronically excited states is far behind ground-state applications. One of the underlying problems is the insufficient smoothness of the fitted potential energy surfaces and other properties in the vicinity of state crossings and conical intersections, which is a prerequisite for an efficient regression. Smooth surfaces can be obtained by switching to the diabatic basis. However, diabatization itself is still an outstanding problem. We overcome these limitations by solving both problems at once. We use a machine learning approach combining clustering and regression techniques to correct for the deficiencies of property-based diabatization which, in return, provides us with smooth surfaces that can be easily fitted. Our approach extends the applicability of property-based diabatization to multidimensional systems. We utilize the proposed diabatization scheme to achieve higher prediction accuracy for adiabatic states and we show its performance by reconstructing global potential energy surfaces of excited states of nitrosyl fluoride and formaldehyde. While the proposed methodology is independent of the specific property-based diabatization and regression algorithm, we show its performance for kernel ridge regression and a very simple diabatization based on transition multipoles. Compared to most other algorithms based on machine learning, our approach needs only a small amount of training data.

Liquid-jet photoemission spectroscopy as a structural tool: site-specific acid-base chemistry of vitamin C.

Liquid-jet photoemission spectroscopy (LJ-PES) directly probes the electronic structure of solutes and solvents. It also emerges as a novel tool to explore chemical structure in aqueous solutions, yet the scope of the approach has to be examined. Here, we present a pH-dependent liquid-jet photoelectron spectroscopic investigation of ascorbic acid (vitamin C). We combine core-level photoelectron spectroscopy and ab initio calculations, allowing us to site-specifically explore the acid-base chemistry of the biomolecule. For the first time, we demonstrate the capability of the method to simultaneously assign two deprotonation sites within the molecule. We show that a large change in chemical shift appears even for atoms distant several bonds from the chemically modified group. Furthermore, we present a highly efficient and accurate computational protocol based on a single structure using the maximum-overlap method for modeling core-level photoelectron spectra in aqueous environments. This work poses a broader question: to what extent can LJ-PES complement established structural techniques such as nuclear magnetic resonance? Answering this question is highly relevant in view of the large number of incorrect molecular structures published.

Predicting the photodynamics of cyclobutanone triggered by a laser pulse at 200 nm and its MeV-UED signals-A trajectory surface hopping and XMS-CASPT2 perspective.

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born-Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S1, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function-ΔPDF(r)] for each (interpolated) pathway as well as a time-resolved signal [ΔPDF(r,t) and ΔI/I(s,t)] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both ΔPDF(r,t) and ΔI/I(s,t)] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.

Photocatalytic Generation of Trifluoromethyl Nitrene for Alkene Aziridination.

N-Trifluoromethylated organics may be applied in drug design, agrochemical synthesis, and materials science, among other areas. Yet, despite recent advances in the synthesis of aliphatic, cyclic and heterocyclic N-trifluoromethyl compounds, no strategy based on trifluoromethyl nitrene has hitherto been explored. Here we describe the formation of triplet trifluoromethyl nitrene from azidotrifluoromethane, a stable and safe-to-use precursor, by visible light photocatalysis. The addition of CF3 N to alkenes via biradical intermediates afforded previously unknown aziridines substituted with trifluoromethyl group on the nitrogen atom. The obtained aziridines were converted into either N-trifluoromethylimidazolines, via formal [3+2] cycloaddition with nitriles, mediated by a Lewis acid, or into N-trifluoromethylaldimines, via ring opening and aryl group migration mediated by a strong Brønsted acid. Our findings open new opportunities for the development of novel classes of N-CF3 compounds with possible applications in the life sciences.

Spin-Vibronic Control of Intersystem Crossing in Iodine-Substituted Heptamethine Cyanines.

Spin-orbit coupling between electronic states of different multiplicity can be strongly coupled to molecular vibrations, and this interaction is becoming recognized as an important mechanism for controlling the course of photochemical reactions. Here, we show that the involvement of spin-vibronic coupling is essential for understanding the photophysics and photochemistry of heptamethine cyanines (Cy7), bearing iodine as a heavy atom in the C3' position of the chain and/or a 3H-indolium core, as potential triplet sensitizers and singlet oxygen producers in methanol and aqueous solutions. The sensitization efficiency was found to be an order of magnitude higher for the chain-substituted than the 3H-indolium core-substituted derivatives. Our ab initio calculations demonstrate that while all optimal structures of Cy7 are characterized by negligible spin-orbit coupling (tenths of cm-1) with no dependence on the position of the substituent, molecular vibrations lead to its significant increase (tens of cm-1 for the chain-substituted cyanines), which allowed us to interpret the observed position dependence.

Excitation and fragmentation of the dielectric gas C4F7N: Electrons vs photons.

C4F7N is a promising candidate for the replacement of sulfur hexafluoride as an insulating medium, and it is important to understand the chemical changes initiated in the molecule by collision with free electrons, specifically the formation of neutral fragments. The first step of neutral fragmentation is electronic excitation, yet neither the absorption spectrum in the vacuum ultraviolet (VUV) region nor the electron energy loss spectrum have previously been reported. Here, we experimentally probed the excited states by VUV photoabsorption spectroscopy and electron energy loss spectroscopy (EELS). We found that the distribution of states populated upon electron impact with low-energy electrons is significantly different from that following photoabsorption. This difference was confirmed and interpreted with ab initio modeling of both VUV and EELS spectra. We propose here a new computational protocol for the simulation of EELS spectra combining the Born approximation with approximate forms of correlated wave functions, which allows us to calculate the (usually very expensive) scattering cross sections at a cost similar to the calculation of oscillator strengths. Finally, we perform semi-classical non-adiabatic dynamics simulations to investigate the possible neutral fragments of the molecule formed through electron-induced neutral dissociation. We show that the product distribution is highly non-statistical.

Probing aqueous ions with non-local Auger relaxation.

Non-local analogues of Auger decay are increasingly recognized as important relaxation processes in the condensed phase. Here, we explore non-local autoionization, specifically Intermolecular Coulombic Decay (ICD), of a series of aqueous-phase isoelectronic cations following 1s core-level ionization. In particular, we focus on Na+, Mg2+, and Al3+ ions. We unambiguously identify the ICD contribution to the K-edge Auger spectrum. The different strength of the ion-water interactions is manifested by varying intensities of the respective signals: the ICD signal intensity is greatest for the Al3+ case, weaker for Mg2+, and absent for weakly-solvent-bound Na+. With the assistance of ab initio calculations and molecular dynamics simulations, we provide a microscopic understanding of the non-local decay processes. We assign the ICD signals to decay processes ending in two-hole states, delocalized between the central ion and neighbouring water. Importantly, these processes are shown to be highly selective with respect to the promoted water solvent ionization channels. Furthermore, using a core-hole-clock analysis, the associated ICD timescales are estimated to be around 76 fs for Mg2+ and 34 fs for Al3+. Building on these results, we argue that Auger and ICD spectroscopy represents a unique tool for the exploration of intra- and inter-molecular structure in the liquid phase, simultaneously providing both structural and electronic information.

Electron-induced fragmentation of water droplets: Simulation study.

The transport of free electrons in a water environment is still poorly understood. We show that additional insight can be brought about by investigating fragmentation patterns of finite-size particles upon electron impact ionization. We have developed a composite protocol aiming to simulate fragmentation of water clusters by electrons with kinetic energies in the range of up to 100 eV. The ionization events for atomistically described molecular clusters are identified by a kinetic Monte Carlo procedure. We subsequently model the fragmentation with classical molecular dynamics simulations, calibrated by non-adiabatic quantum mechanics/molecular mechanics simulations of the ionization process. We consider one-electron ionizations, energy transfer via electronic excitation events, elastic scattering, and also the autoionization events through intermolecular Coulombic decay. The simulations reveal that larger water clusters are often ionized repeatedly, which is the cause of substantial fragmentation. After losing most of its energy, low-energy electrons further contribute to fragmentation by electronic excitations. The simultaneous measurement of cluster size distribution before and after the ionization represents a sensitive measure of the energy transferred into the system by an incident electron.

Bimolecular reactions on sticky and slippery clusters: Electron-induced reactions of hydrogen peroxide.

Nanoparticles can serve as an efficient reaction environment for bimolecular reactions as the reactants concentrate either inside the nanoparticle or on the surface of the nanoparticle. The reaction rate is then controlled by the rate of formation of the reaction pairs. We demonstrate this concept on the example of electron-induced reactions in hydrogen peroxide. We consider two types of nanoparticle environments: solid argon particles, only weakly interacting with the hydrogen peroxide reactant, and ice particles with a much stronger interaction. The formation of hydrogen peroxide dimers is investigated via classical molecular dynamics (MD) simulations on a microsecond timescale. With a modified force field for hydrogen peroxide, we found out a fast formation and stabilization of the hydrogen peroxide dimer for argon nanoparticles, while the reaction pair was formed reversibly at a much slower rate on the water nanoparticles. We have further investigated the electron-induced reactions using non-adiabatic ab initio MD simulations, identifying the possible reaction products upon the ionization or electron attachment. The major reaction path in all cases corresponded to a proton transfer. The computational findings are supported by mass spectrometry experiments, where large ArM and (H2O)M nanoparticles are generated, and several hydrogen peroxide molecules are embedded on these nanoparticles in a pickup process. Subsequently, the nanoparticles are ionized either positively by 70 eV electrons or negatively by electron attachment at electron energies below 5 eV. The recorded mass spectra demonstrate the efficient coagulation of H2O2 on ArM, while it is quite limited on (H2O)M.

Energy partitioning and spin-orbit effects in the photodissociation of higher chloroalkanes.

We investigate the photodissociation dynamics of the C-Cl bond in chloroalkanes CH3Cl, n-C3H7Cl, i-C3H7Cl, n-C5H11Cl, combining velocity map imaging (VMI) experiment and direct ab initio dynamical simulations. The Cl fragment kinetic energy distributions (KEDs) from the VMI experiment exhibit a single peak with maximum close to 0.8 eV, irrespective of the alkyl chain length and C-Cl bond position. In contrary to CH3Cl, where less than 10% of the available energy is deposited into the internal excitation of the CH3 fragment, for all higher chloroalkanes around 40% to 60% of the available energy goes into the alkyl fragment excitation. We apply the classical hard spheres and spectator model to explain the energy partitioning, and compare the classical approach with direct ab initio dynamics simulations. The alkyl chain appears to be a soft, energy absorbing unit. We further investigate the role of the spin-orbit effects on the excitation and dynamics. Combining our experimental data with theory allows us to derive the probability of the direct absorption into the triplet electronic state as well as the probabilities for intersystem crossing. The results indicate an increasing direct absorption into the triplet state with increasing alkyl chain length.

Temperature evolution in IR action spectroscopy experiments with sodium doped water clusters.

The combination of supersonic expansions with IR action spectroscopy techniques is the basis of many successful approaches to study cluster structure and dynamics. The effects of temperature and temperature evolution are important with regard to both the cluster synthesis and the cluster dynamics upon IR excitation. In the past the combination of the sodium doping technique with IR excitation enhanced near threshold photoionization has been successfully applied to study neutral, especially water clusters. In this work we follow an overall examination approach for inspecting the interplay of cluster temperature and cluster structure in the initial cooling process and in the IR excitation induced heating of the clusters. In molecular simulations, we study the temperature dependent photoionization spectra of the sodium doped clusters and the evaporative cooling process by water molecule ejection at the cluster surface. We present a comprehensive analysis that provides constraints for the temperature evolution from the nozzle to cluster detection in the mass spectrometer. We attribute the IR action effect to the strong temperature dependence of sodium solvation in the IR excited clusters and we discuss the effects of geometry changes during the IR multi-photon absorption process with regard to application prospects of the method.

Electron attachment to microhydrated 4-nitro- and 4-bromo-thiophenol.

We investigate the effect of microhydration on electron attachment to thiophenols with halogen (Br) and nitro (NO2) functional groups in the para position. We focus on the formation of anions upon the attachment of low-energy electrons with energies below 8 eV to heterogeneous clusters of the thiophenols with water. For nitro-thiophenol (NTP), the primary reaction channel observed is the associative electron attachment, irrespective of the microhydration. On the other hand, bromothiophenol (BTP) fragments significantly upon the electron attachment, producing Br- and (BTP-H)- anions. Microhydration suppresses fragmentation of both molecules, however in bromothiophenol, the Br- channel remains intense and Br(H2O)n- hydrated fragment clusters are observed. The results are supported by the reaction energetics obtained from ab initio calculations. Different dissociation dynamics of NTP and BTP can be related to different products of their plasmon induced reactions on Au nanoparticles. Computational modeling of the simplified BTP(H2O) system indicates that the electron attachment products reflect the structure of neutral precursor clusters - the anion dissociation dynamics is controlled by the hydration site.

Following in Emil Fischer's Footsteps: A Site-Selective Probe of Glucose Acid-Base Chemistry.

Liquid-jet photoelectron spectroscopy was applied to determine the first acid dissociation constant (pKa) of aqueous-phase glucose while simultaneously identifying the spectroscopic signature of the respective deprotonation site. Valence spectra from solutions at pH values below and above the first pKa reveal a change in glucose's lowest ionization energy upon the deprotonation of neutral glucose and the subsequent emergence of its anionic counterpart. Site-specific insights into the solution-pH-dependent molecular structure changes are also shown to be accessible via C 1s photoelectron spectroscopy. The spectra reveal a considerably lower C 1s binding energy of the carbon site associated with the deprotonated hydroxyl group. The occurrence of photoelectron spectral fingerprints of cyclic and linear glucose prior to and upon deprotonation are also discussed. The experimental data are interpreted with the aid of electronic structure calculations. Our findings highlight the potential of liquid-jet photoelectron spectroscopy to act as a site-selective probe of the molecular structures that underpin the acid-base chemistry of polyprotic systems with relevance to environmental chemistry and biochemistry.

Deciphering the Structure-Property Relations in Substituted Heptamethine Cyanines.

Heptamethine cyanines (Cy7) are fluorophores essential for modern bioimaging techniques and chemistry. Here, we systematically evaluated the photochemical and photophysical properties of a library of Cy7 derivatives containing diverse substituents in different positions of the heptamethine chain. A single substitution allows modulation of their absorption maxima in the range of 693-805 nm and photophysical properties, such as quantum yields of singlet-oxygen formation, decomposition, and fluorescence or affinity to singlet oxygen, within 2-3 orders of magnitude. The same substituent in different positions of the chain often exhibits distinctly contradictory effects, demonstrating that both the type and position of the substituent are pivotal for the design of Cy7-based applications. The combination of experimental results with quantum-chemical calculations provides insights into the structure-property relationship, the elucidation of which will accelerate the development of cyanines with properties tailored for specific applications, such as fluorescent probes and sensors, photouncaging, photodynamic therapy, or singlet-oxygen detection.

Wavelength-Dependent Photochemistry and Biological Relevance of a Bilirubin Dipyrrinone Subunit.

Phototherapy is a standard treatment for severe neonatal jaundice to remove toxic bilirubin from the blood. Here, the wavelength-dependent photochemistry of vinylneoxanthobilirubic acid methyl ester, a simplified model of a bilirubin dipyrrinone subunit responsible for a lumirubin-like structural rearrangement, was thoroughly investigated by liquid chromatography and mass and absorption spectroscopies, with the application of a multivariate curve resolution analysis method supplemented with quantum chemical calculations. Irradiation of the model chromophore leads to reversible Z → E photoisomerization followed by reversible photocyclization to a seven-membered ring system (formed as a mixture of diastereomers). Both the isomerization processes are efficient (ΦZE ∼ ΦEZ ∼ 0.16) when irradiated in the wavelength range of 360-410 nm, whereas the E-isomer cyclization (Φc = 0.006-0.008) and cycloreversion (Φ-c = 0.002-0.004) reactions are significantly less efficient. The quantum yields of all processes were found to depend strongly on the wavelength of irradiation, especially when lower energy photons were used. Upon irradiation in the tail of the absorption bands (490 nm), both the isomers exhibit more efficient photoisomerization (ΦZE ∼ ΦEZ ∼ 0.30) and cyclization (Φc = ∼0.07). In addition, the isomeric bilirubin dipyrrinone subunits were found to possess important antioxidant activities while being substantially less toxic than bilirubin.

Solvation energies of ions with ensemble cluster-continuum approach.

Solvation free energies can be advantageously estimated by cluster-continuum approaches. They proved useful especially for systems with high charge density. However, the clusters are assumed to be single minimum rigid species. It is an invalid condition for larger clusters and it complicates the assessment of convergence with the system size. We present a new variant of the cluster-continuum approach, "Ensemble Cluster-Continuum" scheme, where the single minima problem is circumvented by a thermodynamic cycle based on vertical quantities (ionization energies, electron affinities). Solvation free energies are calculated for a charged-neutralized system and solvation correction for the vertical quantities is estimated for an ensemble of structures from molecular dynamics simulation. We test the scheme on a set of various types of anions and cations, we study the convergence of the cluster-continuum model and assess various types of errors. The quantitative data depend on the applied continuum solvation model yet the convergence is analogous. We argue that the assessment of convergence provides a measure of the reliability of the calculated solvation energies.

Ionization energies in solution with the QM:QM approach.

We discuss a fragment-based QM:QM scheme as a practical way to access the energetics of vertical electronic processes in the condensed phase. In the QM:QM scheme, we decompose the large molecular system into small fragments, which interact solely electrostatically. The energies of the fragments are calculated in a self-consistent field generated by the other fragments and the total energy of the system is calculated as a sum of the fragment energies. We show on two test cases (cytosine and a sodium cation) that the method allows one to accurately simulate the shift of vertical ionization energies (VIE) while going from the gas phase to the bulk. For both examples, the predicted solvent shifts and peak widths estimated at the DFT level agree well with the experimental observations. We argue that the QM:QM approach is more suitable than either an electrostatic embedding based QM/MM approach, a full quantum description at the DFT level with a generally used functional or a combination of both. We also discuss the potential scope of the applicability for other electronic processes such as Auger decay.