Hydronium Ion Transport across the Liquid/Liquid Interface Assisted by a Phase-Transfer Catalyst: Structure and Thermodynamics Using Molecular Dynamics Simulation

Molecular dynamics simulations are used to examine the thermodynamic and structural aspects of the transfer of the classical hydronium ion (H3O+) across a water/1,2-dichloroethane (DCE) interface assisted by the phase-transfer catalyst (PTC) tetrakis(pentafluorophenyl) borate anion (TPFB-). The free energy of transfer from water to DCE of the H3O+-TPFB- ion pair is calculated to be 6 ± 1 kcal/mol, significantly less than that of the free hydronium ion (17 ± 1 kcal/mol). The ion pair is relatively stable at the interface and in the organic phase when it is accompanied by three water molecules with a small barrier to dissociation that supports its utility as a PTC. An examination of the hydration structure that accompanies the transfer of the ion pair shows that the ion pair, like the free hydronium ion, is transferred with the assistance of a finger-like water structure.

Optimized Collective Variable for Collapse Transition in Linear Hydrophobic Polymers: Importance of Hydration Water and End-to-End Distance

Choosing an appropriate collective variable (CV) for any biomolecular process is a challenging task. Researchers are developing methods to solve this issue using a variety of methodologies, most recently using machine learning (ML) methods. In this work, we investigate the mechanism of collapse transition across various lengths of polymer systems through adaptively sampled multiple short trajectories utilizing the Time Lagged Independent Component Analysis (TICA) framework. From TICA analysis, it is revealed that the radius of gyration (Rg) and end-to-end distance serve as good order parameters (OPs) for these systems describing overall energy landscapes. Markov state model (MSM) and mean first passage time (MFPT) analysis suggest that hydration water (Nw) plays a determining role in dictating the time scale and barrier for the collapsed transition for the C40 system. P-fold analysis on identifying transition state ensembles (TSE) identified by committor analysis also strengthens the role of Nw in such a transition. TICA, MSM, and committor analyses on the collapse transition for C45 reveal similarities with C40 systems in different aspects. Furthermore, we propose a pipeline integrating XGBoost regression along with an interpretable ML model, Shapley Additive exPlanation (SHAP) to precisely elucidate the contribution of each OP locally at the TSE. Through this approach, we observe that the collapse transition is primarily driven by Nw for both polymer systems. A carefully designed protocol for the collapsed transition of C60 systems indirectly reiterates the above result. Overall, our results suggest that while the end-to-end distance should be considered for better resolution of metastable states in the landscape, Nw is the crucial coordinate to be used in enhanced sampling for the exploration of actual collapse transitions for linear hydrophobic polymer systems. The Python code for analyzing the contribution of different OPs in the TSE using an ML-aided protocol is available on GitHub (https://github.com/saikat-ai/linear_polymer_project).

Unraveling the excitonics of light emission from metal-halide perovskite quantum dots

Metal halide semicondictor perovskites have been under intense investigation for their promise in light absorptive applications like photovoltaics. They have more recently experienced interest for their promise in light emissive applications. A key aspect of perovskites is their glassy, ionic lattice that exhibits dynamical disorder. One possible result of this dynamical disorder is their strong coupling between electronic and lattice degrees of freedom which may confer remarkable properties for light emission such as defect tolerance. How does the system, comprised of excitons, couple to the bath, comprised of lattice modes? How does this system-bath interaction give rise to novel light emissive properties and how do these properties give insight into the nature of these materials? We review recent work from this group in which time-resolved photoluminescence spectroscopy is used to reveal such insights. Based upon a fast time resolution of 3 ps, energy resolution, and temperature dependence, a wide variety of insights are gleaned. These insights include: lattice contributions to the emission linewidths, multiexciton formation, hot carrier cooling, excitonic fine structure, single dot superradiance, and a breakdown of the Condon approximation, all due to complex structural dynamics in these materials.

On the nature of initial solvation in bulk polar liquids: Gaussian or exponential?

Measurement of time evolution of fluorescence of a probe solute has been a quintessential technique to quantify how dipolar solvent molecules dynamically minimize the free energy of an electronically excited probe. During such solvation dynamics in bulk liquids, a substantial part of relaxation was shown to complete within sub-100 fs from time-gated fluorescence measurements, as also predicted by molecular dynamics simulation studies. However, equivalent quantification of solvation timescales by femtosecond pump-probe and broadband fluorescence measurements revealed an exponential nature of this initial relaxation having quite different timescales. Here, we set out to unveil the reason behind these puzzling contradictions. We introduce a method for estimating probe wavelength-dependent instrument response and demonstrate that the observation of the Gaussian vs exponential nature of initial relaxation is indeed dependent on the method of data analysis. These findings call for further experimental investigation and parallel development of theoretical models to elucidate the molecular-level mechanism accounting for different types of early time solvation.

Nanosecond solvation dynamics in a polymer electrolyte for lithium batteries

Solvation dynamics critically affect charge transport. Spectroscopic experiments and computer simulations show that these dynamics in aqueous systems occur on a picosecond timescale. In the case of organic electrolytes, however, conflicting values ranging from 1 to several 100 picoseconds have been reported. We resolve this conflict by studying mixtures of an organic polymer and a lithium salt. Lithium ions coordinate with multiple polymer chains, resulting in temporary crosslinks. Relaxation of these crosslinks, detected by quasielastic neutron scattering, are directly related to solvation dynamics. Simulations reveal a broad spectrum of relaxation times. The average timescale for solvation dynamics in both experiment and simulation is one nanosecond. We present the direct measurement of ultraslow dynamics of solvation shell break-up in an electrolyte.

Experimental and Theoretical Rotational Diffusion Studies of 7DM4M1M1,8, N-2(1H)-one and 7A4T2H1B-2-one in Series of Alcohol Solvents: Stoke's-Einstein-Debye and Alavi-Waldeck Models

Rotational diffusion studies of two solutes 7-(dimethylamino)-4-methoxy-1-methyl-1,8-naphthyridin-2(1H)-one (7DM4M1M1,8, N-2(1H)-one) and 7-amino-4-(trifluoromethyl)-2H-1-benzopyran-2-one (7A4T2H1B-2-one) having equal volumes but different chemical natures are studied in series of alcohol solvents at 303 K using steady-state methods. HOMO-LUMO, Electron density, Molecular electrostatic potential (MEP), etc., are obtained from computational calculations using Gaussian 09 software. Rotational reorientation times of 7DM4M1M1,8, N-2(1H)-one solute molecule is found to be less than 7A4T2H1B-2-one solute molecule indicates it rotates slowly in chosen solvents. Stoke's-Einstein-Debye (SED) model with stick boundary conditions for the 7A4T2H1B-2-one solute molecule is modeled to describe mechanical friction. Polar solutes along with mechanical friction also experience dielectric friction. Both these frictions being non-separable, the Alavi-Waldeck (AW) model is studied for dielectric friction contribution to the total friction solute experiences in solvents. AW model effectively explains the observed dielectric friction in alcohol solvents.

Solvation dynamics on the diffusion timescale elucidated using energy-represented dynamics theory

Photoexcitation of a solute alters the solute-solvent interaction, resulting in the nonequilibrium relaxation of the solvation structure, often called a dynamic Stokes shift or solvation dynamics. Thanks to the local nature of the solute-solvent interaction, the characteristics of the local solvent environment dissolving the solute can be captured by the observation of this process. Recently, we derived the energy-represented Smoluchowski-Vlasov (ERSV) equation, a diffusion equation for molecular liquids, which can be used to analyze the solvation dynamics on the diffusion timescale. This equation expresses the time development for the solvent distribution on the solute-solvent pair interaction energy (energy coordinate). Since the energy coordinate can effectively treat the solvent flexibility in addition to the position and orientation, the ERSV equation can be utilized in various solvent systems. Here, we apply the ERSV equation to the solvation dynamics of 6-propionyl-2-dimethylamino naphthalene (Prodan) in water and different alcohol solvents (methanol, ethanol, and 1-propanol) for clarifying the differences of the relaxation processes among these solvents. Prodan is a solvent-sensitive fluorescent probe and is thus widely utilized for investigating heterogeneous environments. On the long timescale, the ERSV equation satisfactorily reproduces the relaxation time correlation functions obtained from the molecular dynamics (MD) simulations for these solvents. We reveal that the relaxation time coefficient on the diffusion timescale linearly correlates with the inverse of the translational diffusion coefficients for the alcohol solvents because of the Prodan-solvent energy distributions among the alcohols. In the case of water, the time coefficient deviates from the linear relationship for the alcohols due to the difference in the extent of importance of the collective motion between the water and alcohol solvents.

Spectroscopic view on the interaction between the psoralen derivative amotosalen and DNA

Psoralens are eponymous for PUVA (psoralen plus UV-A radiation) therapy, which inter alia can be used to treat various skin diseases. Based on the same underlying mechanism of action, the synthetic psoralen amotosalen (AMO) is utilized in the pathogen reduction technology of the INTERCEPT® Blood System to inactivate pathogens in plasma and platelet components. The photophysical behavior of AMO in the absence of DNA is remarkably similar to that of the recently studied psoralen 4'-aminomethyl-4,5',8-trimethylpsoralen (AMT). By means of steady-state and time-resolved spectroscopy, intercalation and photochemistry of AMO and synthetic DNA were studied. AMO intercalates with a higher affinity into A,T-only DNA (KD = 8.9 × 10-5 M) than into G,C-only DNA (KD = 6.9 × 10-4 M). AMO covalently photobinds to A,T-only DNA with a reaction quantum yield of ΦR = 0.11. Like AMT, it does not photoreact following intercalation into G,C-only DNA. Femto- and nanosecond transient absorption spectroscopy reveals the characteristic pattern of photobinding to A,T-only DNA. For AMO and G,C-only DNA, signatures of a photoinduced electron transfer are recorded.

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this work, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Caught in the act: real-time observation of the solvent response that promotes excited-state proton transfer in pyranine

Photo-induced excited-state proton transfer (ESPT) reactions are of central importance in many biological and chemical processes. Identifying mechanistic details of the solvent reorganizations that facilitate proton transfer however, is challenging for current experimental and theoretical approaches. Using optical pump THz probe (OPTP) spectroscopy and molecular dynamics simulations, we were able to elucidate the ultrafast changes in the solvation environment for three derivatives of pyranine: the photoacid HPTS, the methoxy derivative MPTS, and the photobase OPTS. Experimentally, we find damped oscillations in the THz signal at short times and our simulations enable their assignment to vibrational energy transfer beatings between the photoexcited chromophore and nearby solvent molecules. The simulations of HPTS reveal strikingly efficient sub-ps energy transfer into a particular solvent mode, that is active near 4 THz, and which can provide the requisite energy required for solvent reorganization promoting proton transfer. Similar oscillations are present in the THz signal for all three derivatives, however the signal is damped rapidly for HPTS (within 0.4 ps) and more slowly for MPTS (within 1.4 ps) and OPTS (within 2.0 ps). For HPTS, we also characterize an additional phonon-like propagation of the proton into the bulk with a 140 ps period and an 83 ps damping time. Thermalization of the solvent occurs on a time scale exceeding 120 ps.

A Thermally Activated Delayed Fluorescence Emitter Investigated by Time-Resolved Near-Infrared Spectroscopy

Emitters for organic light-emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) require small singlet (S₁ )-triplet (T₁ ) energy gaps as well as fast intersystem crossing (ISC) transitions. These transitions can be mediated by vibronic mixing with higher excited states S_n and T_n (n=2, 3, 4, …). For a prototypical TADF emitter consisting of a triarylamine and a dicyanobenzene moiety (TAA-DCN) it is shown that these higher states can be located energetically by time-resolved near-infrared (NIR) spectroscopy.

The photophysics of 2-cyanoindole probed by femtosecond spectroscopy

The photophysics of 2-cyanoindole (2-CI) in solution (water, 2,2,2-trifluoroethanol, acetonitrile‚ and tetrahydrofuran) was investigated by steady-state as well as time resolved fluorescence and absorption spectroscopy. The fluorescence quantum yield of 2-cyanoindole is strongly sensitive to the solvent. In water the quantum yield is as low as 4.4 × 10–4. In tetrahydrofuran, it amounts to a yield of 0.057. For 2-CI dissolved in water, a bi-exponential fluorescence decay with time constants of ∼1 ps and ∼8 ps is observed. For short wavelength excitation (266 nm) the initial fluorescence anisotropy is close to zero. For excitation with 310 nm it amounts to 0.2. In water, femtosecond transient absorption reveals that the fluorescence decay is solely due to internal conversion to the ground state. In aprotic solvents, the fluorescence decay takes much longer (acetonitrile: ∼900 ps, tetrahydrofuran: ∼2.6 ns) and intersystem crossing contributes.

Graphical abstract

Modular synthesis, host-guest complexation and solvation-controlled relaxation of nanohoops with donor-acceptor structures

Carbon nanohoops with donor-acceptor (D-A) structures are attractive electronic materials and biological fluorophores, but their synthesis is usually challenging. Moreover, the preparation of D-A nanohoop fluorophores exhibiting high fluorescence quantum yields beyond 500 nm remains a key challenge. This study presents a modular synthetic approach based on an efficient metal-free cyclocondensation reaction that readily produced nine congeners with D-A or donor-acceptor-donor' (D-A-D') structures, one of which is water-soluble. The tailored molecular design of nanohoops enabled a systematic and detailed study of their host-guest complexation with fullerene, optical properties, and charge transfer (CT) dynamics using X-ray crystallography, fluorescence titration, steady and ultrafast transient absorption spectroscopy, and theoretical calculations. The findings revealed intriguing physical properties associated with D-A motifs, such as tight binding with fullerene, moderate fluorescence quantum yields (37-67%) beyond 540 nm, and unique solvation-controlled CT relaxation of D-A-D' nanohoops, where two CT states (D-A and A-D') can be effectively tuned by solvation, resulting in dramatically changed relaxation pathways in different solvents.

Explicit Projection of Stokes Shifts onto Solvent Motion in an Aqueous Liquid and Linear Response Theory

We investigate the molecular origin of the fluorescence Stokes shift in an aqueous liquid. By examining the speed of energy change, the solvation response function is explicitly projected onto the translational and rotational motions of water molecules for both nonequilibrium relaxation and equilibrium fluctuations. Molecular dynamics simulations of a tryptophan solution show that these two processes have highly consistent dynamics, not only for the total response function but also for the decomposed components in terms of specific molecular movements. We found that the rotational mode governs the relaxation of the Stokes shift, whereas the translational mode contributes non-negligibly with slower dynamics. This consistency implies the similarity of the underlying translational and rotational movements of water molecules as the system is far away from and at equilibrium, supporting the validity of the linear response theory at the molecular level. The decomposition methodology is also applicable to a rigid solvent.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

3-aminoquinoline: a turn-on fluorescent probe for preferential solvation in binary solvent mixtures

3-Aminoquinoline (3AQ) has been used as a fluorescent probe for preferential solvation in hexane-ethanol solvent mixtures. Results of the present experiment have been put into context by comparison with prior observations with 5-aminoquinoline (5AQ) as the probe. 3AQ exhibits a relatively small change of dipole moment (Δμ= 2.2 D) upon photoexcitation, compared to 5AQ (Δμ= 6.1D), which might appear to be a hindrance in the way of its use as a solvation probe. Indeed, the values of parameters like spectral shifts are smaller for the present experiment with 3AQ. At the smallest concentration of alcohol used, its local mole fraction around the probe is significantly lower than in the previous experiments with 5AQ. However, these apparent disadvantages are outweighed by the significant increase in fluorescence intensity and lifetime observed with increasing concentration of ethanol in the solvent mixture, as opposed to the drastic fluorescence quenching that occurs for 5AQ. This is a marked advantage in the use of 3AQ in studies like the present one. The local mole fraction of ethanol and preferential solvation index experienced by 3AQ are in line with those reported for 5AQ. The disadvantage of the smaller magnitude of Δμpersists in the time resolved fluorescence experiments, for solvent mixtures with very low ethanol content. Negligible wavelength dependence of fluorescence transients of 3AQ is observed forx_p= 0.002,. However, this effect is outweighed at higher alcohol concentrations, for which nanosecond dynamics of preferential solvation is observed.

Choline chloride and ethylene glycol based deep eutectic solvent (DES) versus hydroxyl functionalized room temperature ionic liquids (RTILs): assessing the differences in microscopic behaviour between the DES and RTILs

With the aim of understanding the differences in the behavior of deep eutectic solvents (DESs) and room temperature ionic liquids (RTILs) in terms of their structure, dynamics, and intra- and intermolecular interactions, three different ILs and one DES having similar functionalities (hydroxyl) have been investigated by using both ensembled average and single-molecule spectroscopic techniques. Specifically, for this purpose, a choline chloride based DES (ethaline) and three hydroxyl functionalized ILs (1-(2-hydroxyethyl)-3-imidazolium bis(trifluoromethanesulfonyl)imide ([OHEMIM][NTF₂]), N-(2-hydroxyl ethyl)-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([OHEMPy][NTf₂]), and N-(2-hydroxyethyl)-N,N-dimethylpropan-1-aminium bis(trifluoromethanesulfonyl)imide ([OHC3CH][NTf₂])) are employed and investigated by EPR, time-resolved fluorescence, NMR and FCS studies. Estimation of polarity through EPR spectroscopy has revealed that the hydroxyl ILs employed in these studies are hyper-polar (close to water) in nature, whereas the polarity of the DES is found to be close to those of aliphatic polyhydroxy-alcohols. Interestingly, both time-resolved fluorescence anisotropy and FCS studies on these systems have suggested that the hydroxyl ILs are more dynamically heterogeneous than the DES. More interestingly, PFG-NMR measurements have indicated that the fluid structure of ethaline is relatively more associated as compared to those of the ILs despite the fact that all the cations have the same hydroxyl functionalities. All these investigations have essentially demonstrated that, despite having similar functionalities, both the DES and hydroxyl ILs employed in the present study exhibit microscopic behaviours that are significantly different from each other, indicating the interplay of various intermolecular interactions within the constituent species in governing the behaviours of these solvent systems.

Selective desolvation in two-step nucleation mechanism steers crystal structure formation

The two-step nucleation (TSN) theory and crystal structure prediction (CSP) techniques are two disjointed yet popular methods to predict nucleation rate and crystal structure, respectively. The TSN theory is a well-established mechanism to describe the nucleation of a wide range of crystalline materials in different solvents. However, it has never been expanded to predict the crystal structure or polymorphism. On the contrary, the existing CSP techniques only empirically account for the solvent effects. As a result, the TSN theory and CSP techniques continue to evolve as separate methods to predict two essential attributes of nucleation - rate and structure. Here we bridge this gap and show for the first time how a crystal structure is formed within the framework of TSN theory. A sequential desolvation mechanism is proposed in TSN, where the first step involves partial desolvation to form dense clusters followed by selective desolvation of functional groups directing the formation of crystal structure. We investigate the effect of the specific interaction on the degree of solvation around different functional groups of glutamic acid molecules using molecular simulations. The simulated energy landscape and activation barriers at increasing supersaturations suggest sequential and selective desolvation. We validate computationally and experimentally that the crystal structure formation and polymorph selection are due to a previously unrecognized consequence of supersaturation-driven asymmetric desolvation of molecules.

Solvent-Dependent Excited-State Evolution of Prodan Dyes

Excited-state character and dynamics of two 6-(dimethylamino)-2-acylnaphthalene dyes (Prodan and Badan-SCH₂CH₂OH) were studied by picosecond time-resolved IR spectroscopy (TRIR) in solvents of different polarity and relaxation times: hexane, CD₃OD, and glycerol-d₈. In all these solvents, near-UV excitation initially produced the same S₁(ππ*) excited state characterized by a broad TRIR signal. A very fast decay (3, ∼100 ps) followed in hexane, whereas conversion to a distinct IR spectrum with a ν(C═O) band downshifted by 76 cm^-1 occurred in polar/H-bonding solvents, slowing down on going from CD₃OD (1, 23 ps) to glycerol-d₈ (5.5, 51, 330 ps). The final relaxed excited state was assigned as planar Me₂N → C═O intramolecular charge transfer S₁(ICT) by comparing experimental and TDDFT-calculated spectra. TRIR conversion kinetics are comparable to those of early stages of multiexponential fluorescence decay and dynamic fluorescence red-shift. This work presents a strong evidence that Prodan-type dyes undergo solvation-driven charge separation in their S₁ state, which is responsible for the dynamic fluorescence Stokes shift observed in polar/H-bonding solvents. The time evolution of the optically prepared S₁(ππ*) state to the S₁(ICT) final state reflects environment relaxation and solvation dynamics. This finding rationalizes the widespread use of Prodan-type dyes as probes of environment dynamics and polarity.

Solvent Effects on Ultrafast Photochemical Pathways

Photochemical reactions are increasingly being used for chemical and materials synthesis, for example, in photoredox catalysis, and generally involve photoexcitation of molecular chromophores dissolved in a liquid solvent. The choice of solvent influences the outcomes of the photochemistry because solute-solvent interactions modify the energies of and crossings between electronic states of the chromophores, and they affect the evolving structures of the photoexcited molecules. Ultrafast laser spectroscopy methods with femtosecond to picosecond time resolution can resolve the dynamics of these photoexcited molecules as they undergo structural and electronic changes, relax back to the ground state, dissipate their excess internal energy to the surrounding solvent, or undergo photochemical reactions. In this Account, we illustrate how experimental studies using ultrafast lasers can reveal the influences that different solvents or cosolutes exert on the photoinduced nonadiabatic dynamics of internal conversion and intersystem crossing in nonradiative relaxation pathways. Although the environment surrounding a solute molecule is rapidly changing, with fluctuations in the coordination to neighboring solvent molecules occurring on femtosecond or picosecond time scales, we show that it is possible to photoexcite selectively only those molecular chromophores transiently experiencing specific solute-solvent interactions such as intermolecular hydrogen bonding.The effects of different solvation environments on the photodynamics are illustrated using four selected examples of photochemical processes in which the solvent has a marked effect on the outcomes. We first consider two aromatic carbonyl compounds, benzophenone and acetophenone, which are known to undergo fast intersystem crossing to populate the first excited triplet state on time scales of a few picoseconds. We show that the nonadiabatic excited-state dynamics are modified by transient hydrogen bonding of the carbonyl group to a protic solvent or by coordination to a metal cation cosolute. We then examine how different solvents modify the competition between two alternative relaxation pathways in a photoexcited UVA-sunscreen molecule, diethylamino hydroxybenzoyl hexyl benzoate (DHHB). This relaxation back to the ground electronic state is an essential part of the effective operation of the sunscreen compound, but the dynamics are sensitive to the surrounding environment. Finally, we consider how solvents of different polarity affect the energies and lifetimes of excited states with locally excited or charge-transfer character in heterocyclic organic compounds used as excited-state electron donors for photoredox catalysis. With these and other examples, we seek to develop a molecular level understanding of how the choice of solution environment might be used to control the outcomes of photochemical reactions.

Ultrafast Rotational and Translational Energy Relaxation in Neat Liquids

The excess energy flow pathways during rotational and translational relaxation induced by rotational or translational excitation of a single molecule of and within each of four different neat liquids (H₂O, MeOH, CCl₄, and CH₄) are studied using classical molecular dynamics simulations and energy flux analysis. For all four liquids, the relaxation processes for both types of excitation are ultrafast, but the energy flow is significantly faster for the polar, hydrogen-bonded (H-bonded) liquids H₂O and MeOH. Whereas the majority of the initial excess energy is transferred into hindered rotations (librations) for rotational excitation in the H-bonded liquids, an almost equal efficiency for transfer to translational and rotational motions is observed in the nonpolar, non-H-bonded liquids CCl₄ and CH₄. For translational excitation, transfer to translational motions dominates for all liquids. In general, the energy flows are quite local; i.e., more than 70% of the energy flows directly to the first solvent shell molecules, reaching almost 100% for CCl₄ and CH₄. Finally, the determined validity of linear response theory for these nonequilibrium relaxation processes is quite solvent-dependent, with the deviation from linear response most marked for rotational excitation and for the nonpolar liquids.

The Role of Hydrogen Bonding in the Preferential Solvation of 5-Aminoquinoline in Binary Solvent Mixtures

5-Aminoquinoline (5AQ) has been used as a fluorescent probe of preferential solvation (PS) in binary solvent mixtures in which the nonpolar component is diethyl ether and the polar component is protic (methanol) or aprotic (acetonitrile). Hence, the roles of solvent polarity and solute-solvent hydrogen bonding have been delineated. Positive deviations of spectral shifts from a linear dependence on the concentration of the polar component, signifying PS, are markedly more pronounced in case of the protic solvent. Solvation dynamics on a nanosecond time scale mark the formation of the solvation shell around the fluorescent probe. Time-resolved area-normalized emission spectra indicate the occurrence of the continuous solvation of the excited state when the polar component is acetonitrile. In contrast, two distinct states were observed when the polar component was methanol, the second state being the hydrogen bonded one. Translational diffusion is the rate-determining step for formation of the solvation shell. The time constant associated with it has been estimated from rise times observed in fluorescence transients monitored at the red end of the fluorescence spectra and also from the time evolution of the spectral width of time-resolved emission spectra.

Intramolecular Charge Transfer of 1-Aminoanthraquinone and Ultrafast Solvation Dynamics of Dimethylsulfoxide

The intramolecular charge transfer (ICT) of 1-aminoanthraquinone (AAQ) in the excited state strongly depends on its solvent properties, and the twisted geometry of its amino group has been recommended for the twisted ICT (TICT) state by recent theoretical works. We report the transient Raman spectra of AAQ in a dimethylsulfoxide (DMSO) solution by femtosecond stimulated Raman spectroscopy to provide clear experimental evidence for the TICT state of AAQ. The ultrafast (~110 fs) TICT dynamics of AAQ were observed from the major vibrational modes of AAQ including the ν_C-N + δ_CH and ν_C=O modes. The coherent oscillations in the vibrational bands of AAQ strongly coupled to the nuclear coordinate for the TICT process have been observed, which showed its anharmonic coupling to the low frequency out of the plane deformation modes. The vibrational mode of solvent DMSO, ν_S=O showed a decrease in intensity, especially in the hydrogen-bonded species of DMSO, which clearly shows that the solvation dynamics of DMSO, including hydrogen bonding, are crucial to understanding the reaction dynamics of AAQ with the ultrafast structural changes accompanying the TICT.

From fluorogens to fluorophores by elucidation and suppression of ultrafast excited state processes of a Schiff base

Strategies have been explored for developing strongly fluorescent species out of a weakly fluorescent Schiff base, 2-(((pyridin-2-ylmethyl)imino)methyl)phenol (salampy). The locally excited enolic state of salampy undergoes an intramolecular proton transfer with a time constant of ca. 200 fs. The emissive cis keto state thus formed decays completely within 50 ps. Its fast decay and miniscule fluorescence quantum yield are attributed to efficient non-radiative channels associated with conformational relaxation. The anionic form, salampy^-, has a significantly longer fluorescence lifetime of 800 ps. Its emissive state evolves in tens of picoseconds, from the locally excited state, by solvent and conformational relaxation. Both the neutral and anionic forms have a fluorescence lifetime of about 6 ns at 77 K, a temperature at which all activated nonradiative channels are blocked. This lifetime is similar to that obtained at room temperature, upon rigidification of the anion by complexation with Zn²⁺. Two such complexes have been studied. The first is binuclear, with acetate bridge between the two Zn²⁺ ions. The second, with ClO₄^- as the counterion, is mononuclear with two salampy ligands ligating the metal ion. Unlike a previous report on a different Schiff base, in which the ligands are π-stacked in its dimeric Zn²⁺ complex, no additional nonradiative deactivation pathway opens up in the Zn complexes of salampy, which are devoid of such stacking. The complex of salampy with Al³⁺ has an even longer fluorescence lifetime of 9 ns, indicating a greater degree of rigidification and consequent suppression of nonradiative processes.

Optical Kerr Effect of Liquid Acetonitrile Probed by Femtosecond Time-Resolved X-ray Liquidography

Optical Kerr effect (OKE) spectroscopy is a method that measures the time-dependent change of the birefringence induced by an optical laser pulse using another optical laser pulse and has been used often to study the ultrafast dynamics of molecular liquids. Here we demonstrate an alternative method, femtosecond time-resolved X-ray liquidography (fs-TRXL), where the microscopic structural motions related to the OKE response can be monitored using a different type of probe, i.e., X-ray solution scattering. By applying fs-TRXL to acetonitrile and a dye solution in acetonitrile, we demonstrate that different types of molecular motions around photoaligned molecules can be resolved selectively, even without any theoretical modeling, based on the anisotropy of two-dimensional scattering patterns and extra structural information contained in the q-space scattering data. Specifically, the dynamics of reorientational (libration and orientational diffusion) and translational (interaction-induced motion) motions are captured separately by anisotropic and isotropic scattering signals, respectively. Furthermore, the two different types of reorientational motions are distinguished from each other by their own characteristic scattering patterns and time scales. The measured time-resolved scattering signals are in excellent agreement with the simulated scattering signals based on a molecular dynamics simulation for plausible molecular configurations, providing the detailed structural description of the OKE response in liquid acetonitrile.

Synthesis and Photophysics of Water-Soluble Psoralens with Red-Shifted Absorption

8-Methoxypsoralen (8-MOP) serves as a PUVA (psoralen + UV-A) agent in the treatment of certain skin diseases. Derivatives of 8-MOP with cationic aromatic substituents at the five positions were synthesized and characterized by steady-state, femtosecond and nanosecond spectroscopy as well as cyclic voltammetry. The aromatic substituents' positive charge increases the water solubility and the affinity toward intercalation into DNA. The aromatic substituents were supposed to lower the psoralen S₁ energy and thereby suppress a photo-induced electron transfer (PET) with guanine-bearing DNA. Such a suppression of this PET is expected to increase the propensity of psoralens to photo-addition to DNA. For derivatives bearing methylpyridinium residues, femtosecond spectroscopy revealed an intramolecular PET occurring on the picosecond time scale. This PET precludes the population of the triplet state. As triplet states are the precursor state for the photo-addition to DNA, their intermolecular PET renders these derivatives ineffective in terms of PUVA. For two derivatives bearing trimethylphenylammonium moieties, such an intramolecular PET does not occur and the triplet state is populated. Surprisingly, these compounds also exhibit no PUVA activity. Based on these findings, implications for further optimization of PUVA agents are discussed.

Acridones: Strongly Emissive HIGHrISC Fluorophores

An acridone derivative (N-methyl-difluoro-acridone, NMA-dF) is characterized with respect to its utility as an emitter in organic light emitting diodes (OLEDs). Using steady-state and time-resolved spectroscopy as well as quantum chemistry, its ability to convert singlet and triplet excitons into light was scrutinized. NMA-dF emits in the deep blue range of the visible spectrum. Its fluorescence emission occurs with quantum yields close to 1 and a radiative rate constant of ≈5 × 10⁸ s^-1. So, it processes singlet excitons very efficiently. Using 1,4-dichlorobenzene as a sensitizer, it is shown that NMA-dF also converts triplet excitons into light. With the aid of quantum chemistry, this is related to a reverse intersystem crossing starting from a higher triplet state (HIGHrISC).

Effects of solvation shell relaxation on chain association mechanisms in poly(3-hexylthiophene) solutions

Using poly(3-hexylthiophene) (P3HT) as a model conjugated polymer and atomistic molecular dynamics simulations with carefully verified force fields, we performed in-depth investigations of solvation shell properties of P3HT chains (15 repeating units per chain) in two representative groups of non-polar (or aprotic) organic solvents (better solvents: ortho-dichlorobenzene, bromobenzene, and chlorobenzene; poorer solvents: chloroform, para-xylene, and toluene). We demonstrated that solvation shell relaxation properties in P3HT solutions dictate the formation of regular π-π associations and, hence, crystallinity through the initial chain association and subsequent chain sliding. In contrast, the mean features of polymer-solvent interactions, including solvation free energy and radial distribution function, present little or no difference for all solvent media investigated. Better-solvent media were revealed to bear relatively large values of the first solvation shell relaxation time (τ1 ≫ 100 ps) as well as larger ratios of relaxation times for the first two solvation shells (τ1/τ2 > 2), and vice versa for poorer-solvent media (τ1 ≪ 100 ps and τ1/τ2 < 2). The linear hexyl side-chain unit was noted to substantially enlarge both quantities while notably reducing the solvation free energy as well. As discussed herein, these findings shed new light on the mechanistic features by which solvent quality impacts the degree of π-π association crucial for modern applications with crystalline conjugated polymers.

Strong Electronic Coupling-Induced Ultrafast Charge Transfer in Donor-Pyrene-Acceptor Systems

In this study, we decipher the charge transfer (CT) processes in donor-pyrene-acceptor (DPA) molecules via various time-resolved spectroscopic measurements. It has been challenging to unravel the ultrafast CT dynamics in DPA molecules because they exhibit an initial CT emission in the same spectral range as the locally excited (LE) emission. However, we finally observed the CT rate of ∼200 fs in DPA molecules from the time-resolved fluorescence anisotropy decay profiles. Our measurements allow us to suggest that the LE and CT states of DPA systems have isoenergetic potential surfaces and that the introduction of the acceptor to the pyrene moiety gives rise to strong electronic coupling between the LE and CT states. Therefore, we determined that this solvent-independent ultrafast CT occurs through the adiabatic potential energy surface and that the CT characteristics are enhanced in DPA compared to the donor-pyrene-donor system.

Solvated Electrons in Solids-Ferroelectric Large Polarons in Lead Halide Perovskites

Solvation plays a pivotal role in chemistry and biology. A solid-state analogy of solvation is polaron formation, but the magnitude of Coulomb screening is typically an order of magnitude weaker than that of solvation in aqueous solutions. Here, we describe a new class of polarons, the ferroelectric large polaron, proposed initially by Miyata and Zhu in 2018 (Miyata, K.; Zhu, X.-Y. Ferroelectric Large Polarons. Nat. Mater. 2018, 17 (5), 379-381). This type of polaron allows efficient Coulomb screening of an electron or hole by extended ordering of dipoles from symmetry-broken unit cells. The local ordering is reflected in the ferroelectric-like THz dielectric responses of lead halide perovskites (LHPs) and may be partially responsible for their exceptional optoelectronic performances. Despite the likely absence of long-range ferroelectricity in LHPs, a charge carrier may be localized to and/or induce the formation of nanoscale domain boundaries of locally ordered dipoles. Based on the known planar nature of energetically favorable domain boundaries in ferroelectric materials, we propose that a ferroelectric polaron localizes to planar boundaries of transient polar nanodomains. This proposal is supported by dynamic simulations showing sheet-like transient electron or hole wave functions in LHPs. Thus, the Belgian-waffle-shaped ferroelectric polaron in the three-dimensional LHP crystal structure is a large polaron in two dimensions and a small polaron in the perpendicular direction. The ferroelectric large polaron may form in other crystalline solids characterized by dynamic symmetry breaking and polar fluctuations. We suggest that the ability to form ferroelectric large polarons can be a general principle for the efficient screening of charge carriers from scattering with other charge carriers, with charged defects and with longitudinal optical phonons, thus contributing to enhanced optoelectronic properties.

Implications of relaxation dynamics of collapsed conjugated polymeric nanoparticles for light-harvesting applications

The mechanism of the formation of nanoparticles (collapsed state) from the extended state of polymers and their ultrafast excited state relaxation dynamics are illustrated.

Solvent Effect on Supramolecular Self-Assembly of Chlorophylls a on Chemically Reduced Graphene Oxide

Solvent plays an important role in the surface interaction of molecules. In this study, we use "chlorophyll a", an archetypical molecule, to investigate its supramolecular self-assembly with chemically reduced graphene oxide in three different types of solvents: polar protic, polar aprotic, and non-polar. It was observed that only a polar protic solvent that can donate protons facilitates the hydrogen bonding between chlorophyll a and chemically reduced graphene oxide nanosheets in a hybrid system. The formation of hydrogen bonds further initiates the other non-covalent interactions such as π-π stacking and hydrophobic interaction, which altogether play a key driving force for supramolecular self-assembly of chlorophylls on chemically reduced graphene oxides. The experimental results are strongly supported by density functional theory calculations, which show robust electron coupling between chlorophylls and chemically reduced graphene oxide.

Solvation Dynamics in Water. 4. On the Initial Regime of Solvation Relaxation

It is shown, by means of numerical and analytic work, that initial molecular momenta play little significant role in the initial fast solvation relaxation that follows electronic excitation of, and charge creation for, a standard model system of a solute in water. Instead, the nonequilibrium dynamics are predominantly described by noninertial "steering" by the torques directly generated by the newly created charge distribution. It is this process that largely overcomes inertia and drives the relaxation dynamics on a time scale of a few tens of femtoseconds in the key initial regime of the dynamics. These results are discussed in the context of commonly employed descriptions such as inertial, Gaussian, and underdamped dynamical behavior.

Time-Resolved Probing of the Nonequilibrium Structural Solvation Dynamics by the Time-Dependent Stokes Shift

The time-dependent fluorescence Stokes shift monitors the relaxation of the polarization of a polar solvent in the surroundings of a photoexcited solute molecule but also the structural variation of the solute following photoexcitation and the subsequent molecular charge redistribution. Here, we formulate a simple nonequilibrium quantum theory of solvation for an explicitly time-dependent continuous solvent. The time-dependent solvent induces nonequilibrium fluctuations on the solvent dynamics which are directly reflected in different time components in the time-dependent Stokes shift. We illustrate the structural dynamics in the presence of an explicitly time-dependent solvent by the example of a dynamically shrinking solute which leads to a bimodal Stokes shift. Interestingly, both contributions are mutually coupled. Furthermore, we can explain a prominent long-tail decay of the Stokes shift associated with slow structural dynamical variations.

Direct Observation of Different One- and Two-Photon Fluorescent States in a Pyrrolo[3,2-b]pyrrole Fluorophore

Two-photon fluorophores are frequently employed to obtain superior spatial resolution in optical microscopy applications. To guide the rational design of these molecules, a detailed understanding of their excited-state deactivation pathways after two-photon excitation is beneficial, especially to assess the often-assumed presumption that the one- and two-photon excited-state dynamics are similar after excitation. Here, we showcase the breakdown of this assumption for one- and two-photon excitation of a centrosymmetric pyrrolo[3,2-b]pyrrole chromophore by combining time-resolved fluorescence and broadband femtosecond transient absorption spectroscopy. Compared to one-photon excitation, where radiative decay dominates the photodynamics, two-photon excitation leads to dynamics arising from increased nonradiative decay pathways. These different photodynamics are manifest to different quantum yields, thus highlighting the types of time-resolved studies described here to be valuable guideposts in the design of two-photon fluorophores for imaging applications.

To unravel the connection between the non-equilibrium and equilibrium solvation dynamics of tryptophan: success and failure of the linear response theory of fluorescence Stokes shift

The connections between the non-equilibrium solvation dynamics upon optical transitions and the system's equilibrium fluctuations are explored in aqueous liquid. Linear response theory correlates time-dependent fluorescence with the equilibrium time correlation functions. In the previous work [T. Li, J. Chem. Theory Comput., 2017, 13, 1867], Stokes shift was explicitly decomposed into the contributions of various order time correlation functions on the excited state surface. Gaussian fluctuations of the solute-solvent interactions validate linear response theory. Correspondingly, the deviation of the Gaussian statistics causes the inefficiency of linear response evaluation. The above mechanism is thoroughly tested in this study. By employing molecular simulations, multiple non-equilibrium processes, not necessarily initiated from the ground state equilibrium minimum, were examined for tryptophan. Both the success and failure of linear response theory are found for this simple system and the mechanism is analyzed. These observations, assisted by the width dynamics, the initial state linear response approach, and the variation of the solvation structures, integrally verify the virtue of the excited state Gaussian statistics on the dynamics of Stokes shift.

Self-thermophoresis at the nanoscale using light induced solvation dynamics

Downsizing microswimmers to the nanoscale, and using light as an externally controlled fuel, are two important goals within the field of active matter. Here we demonstrate using all-atom molecular dynamics simulations that solvation relaxation, the solvent dynamics induced after visible light electronic excitation of a fluorophore, can be used to propel nanoparticles immersed in polar solvents. We show that fullerenes functionalized with fluorophore molecules in liquid water exhibit substantial enhanced mobility under external excitation, with a propulsion speed proportional to the power dissipated into the system. We show that the propulsion mechanism is quantitatively consistent with a molecular scale instance of self-thermophoresis. Strategies to direct the motion of functionalized fullerenes in a given direction using confined environments are also discussed.

Lennard-Jones Parameters Determined to Reproduce the Solubility of NaCl and KCl in SPC/E, TIP3P, and TIP4P/2005 Water

Most classical nonpolarizable ion potential models underestimate the solubility values of NaCl and KCl in water significantly. We determine Lennard-Jones parameters of Na⁺, K⁺, and Cl^- that reproduce the solubility as well as the hydration free energy in dilute aqueous solutions for three water potential models, SPC/E, TIP3P, and TIP4P/2005. The ion-oxygen distance in the solution and the cation-anion distance in salt are also considered in the parametrization. In addition to the target properties, the hydration enthalpy, hydration entropy, self-diffusion coefficient, coordination number, lattice energy, enthalpy of solution, density, viscosity, and number of contact ion pairs are calculated for comparison with 17 frequently used or recently developed ion potential models. The overall performance of each ion model is represented by a global score using a scheme that was originally developed for comparison of water potential models. The global score is better for our models than for the other 17 models not only because of the quite good prediction for the solubility but also because of the relatively small deviation from the experimental value for many of the other properties.