Solvation Dynamics Research Articles

Symmetry-Breaking Charge-Separation in a Subphthalocyanine Dimer Resolved by Two-Dimensional Electronic Spectroscopy.

Understanding the role of structural and environmental dynamics in the excited state properties of strongly coupled chromophores is of paramount importance in molecular photonics. Ultrafast, coherent, and multidimensional spectroscopies have been utilized to investigate such dynamics in the simplest model system, the molecular dimer. Here, we present a half-broadband two-dimensional electronic spectroscopy (HB2DES) study of the previously reported ultrafast symmetry-breaking charge separation (SB-CS) in the subphthalocyanine oxo-bridged homodimer μ-OSubPc2. Electronic structure calculations and 2D cross-peaks reveal the dimer's excitonic structure, while ultrafast evolution of the multidimensional spectra unveils subtle features of structural relaxation, solvation dynamics, and inhomogeneous broadening in the SB-CS. Analysis of coherently excited vibrational motions reveals dimer-specific low-frequency Raman active modes coupled to higher-frequency vibrations localized on the SubPc cores. Finally, beatmap amplitude distributions characteristic of excitonic dimers with multiple bright states are reported and analyzed.

FULLY ANALYTICAL SOLUTION FOR THE FIRST SPECTRAL MOMENT OF A FLUOROPHORE IN SOLUTION

One of the common methods for assessing the solvation dynamics of a fluorescent system in a polar medium is the analysis of the first spectral moment behavior. So there is a need to build a model that takes into account important physical processes and at the same time is not demanding in terms of calculations. In this paper, an approach is proposed that considers the parameters of the exciting laser pulse and intramolecular vibrational relaxation for constructing the solvation dynamics. A method for building an analytical solution for the dynamics of the first spectral moment is presented. This model is free of numerical methods using, and it is a simple construction built by basics functions and an additional error function. The influence of all key parameters on the system relaxation is illustrated. It is shown that delocalization in time of the pump pulse suppresses the oscillatory behavior of the first moment and hides information about the damped intramolecular mode. The similar effect is also observed with the diffusion and inertial processes contribution insreasion to the reorganization energy at early stages of relaxation. The presented model can be generalized to take into account quantum transitions between vibrational sublevels of high-frequency intramolecular vibrations.

Solvation Structure and Dynamics of the Thiocyanate Anion in mixed N,N-Dimethylformamide-Water Solvents: A Molecular Dynamics Approach.

The solvation structure and dynamics of the thiocyanate anion at infinite dilution in mixed N, N-Dimethylformamide (DMF)-water liquid solvents was studied using classical molecular dynamics simulation techniques. The results obtained have indicated a preferential solvation of the thiocyanate anions by the water molecules, due to strong hydrogen bonding interactions between the anion and water molecules. A first hydration shell at short intermolecular distances is formed around the SCN- anion consisting mainly by water molecules, followed by a second shell consisting by both DMF and water molecules. The strong interactions between the thiocyanate anion and water molecules are further reflected upon the calculated intermittent residence lifetimes of water and DMF in the first and second solvation shells. The dependence of the reorientational relaxation times of the thiocyanate anion upon the mole fraction of DMF in the mixtures has been found to be in good agreement with experiment, revealing strong concentration effects upon these relaxation phenomena. An appreciable solvent composition effect upon the low frequency intermolecular vibrations, due to the anion-water interactions, has also been revealed by calculating the atomic velocity correlation functions and corresponding spectral densities of the anion.

Interplay of ethaline and water dynamics in a hydrated eutectic solvent: Deuteron and oxygen magnetic resonance studies of aqueous ethaline.

For many technological processes, the impact of water addition on the properties of deep eutectic solvents is of central importance. In this context, the impact of hydration on the reorientational dynamics of the deep eutectic solvent (DES) ethaline, a 2:1 molar mixture of ethylene glycol and choline chloride, was studied. Its overall response was explored by means of shear mechanical rheology. To achieve component-selective insights into the dynamics of this material, isotope-edited deuteron and oxygen spin-lattice and spin-spin relaxometry, as well as stimulated-echo spectroscopy, were applied and yielded motional correlation times from above room temperature down to the highly viscous regime. For all temperatures, the cholinium anion was found to reorient about two times slower than ethylene glycol, while the water and the ethylene glycol molecules display very similar mobilities. While hydration enhances the component dynamics with respect to that of dry ethaline, the present findings reveal that it does not detectably increase the heterogeneity of the solvent. Merely, the time scale similarity that is found for the hydrogen bond donor and the water molecules over a particularly wide temperature range impressively attests to the stability of the native solvent structure in the "water-in-DES" regime.

Unraveling Hydration Shell Dynamics and Viscosity Effects Around Cyanamide Probes via 2D IR Spectroscopy.

Hydration dynamics and solvent viscosity play critical roles in the structure and function of biomolecules. An overwhelming body of evidence suggests that protein and membrane fluctuations are closely linked to solvent fluctuations. While extensive research exists on the use of vibrational probes to detect local interactions and solvent dynamics, fewer studies have explored how the behavior of these reporters changes in response to bulk viscosity. To address this gap, two-dimensional infrared spectroscopy (2D IR) was employed in this study to investigate the ultrafast hydration dynamics around a cyanamide (NCN) probe attached to a nucleoside, deoxycytidine, in aqueous solutions with varying glycerol content. The use of a small vibrational probe on a targeted nucleic acid offers the potential to capture more localized hydration dynamics than alternative methods. The time scales for the frequency correlation decays were found to increase linearly with bulk viscosity, ranging from 0.9 to 11.4 ps over viscosities of 0.96-49.1 cP. Additionally, molecular dynamics (MD) simulations were performed to model the local hydration dynamics around the NCN probe. Interestingly, increasing the glycerol content did not significantly alter the hydration of the deoxycytidine. The MD simulations further suggested that the NCN probe's frequency fluctuations were primarily influenced by the dynamics of water in the second solvation shell. Cage correlation functions, which measure the movement of water molecules in and out of the second solvation shell, exhibited decays that closely matched those of the frequency-fluctuation correlation function (FFCF). These findings offer new insights into hydration dynamics and the impact of viscosity on biological systems.

Effects of Cosolvent and Nonsolvating Solvent on the Structural Dynamics of Organic Electrolytes in Sodium-Ion Batteries.

In sodium-ion batteries (SIBs), the performance of a single solvent often does not meet actual requirements and a cosolvent or nonsolvating solvent is needed. However, the effect of these electrolyte additives on the solvation structure and dynamics of Na+ in SIBs is yet to be fully understood. Herein, electrolyte structural dynamics are examined for NaPF6 in dimethyl carbonate (DMC) with 1,1,2,2-tetrafluoro-2,2,2-trifluoroethoxy ethane (HFE) as the nonsolvating solvent or propylene carbonate (PC) as the cosolvent using steady-state and time-resolved infrared (IR) spectroscopies. Molecular dynamics simulations show that the solvation size of Na+ decreases with a loosened structure upon adding the nonsolvating solvent, whereas its first solvation shell becomes denser and the second one becomes softened with decreased participation of the PF6- anion upon adding the cosolvent. While a decreased participation of DMC in the solvation layer of Na+ is suggested by linear IR results, an increased structural inhomogeneity (and hence overall more dynamical fluctuations) is found for Na+-coordinated DMC upon adding both additives by two-dimensional (2D) IR spectroscopy where the carbonyl stretch is used as a probe. The Na+/DMC/PC complex is found to structurally evolve slower than both Na+/DMC and Na+/DMC/HFE complexes on the picosecond time scales by spectral diffusion dynamics extracted from 2D IR diagonal signals. Frontier orbital theory calculations also indicate that both solvent additives are beneficial to increasing the stability of PF6-. The results obtained in this work provide important insights into the roles played by the solvent additives in influencing the solvation structures and dynamics of Na+, which are critical for understanding the Na+ transfer mechanism in SIBs.

Role of Hydration and Amino Acid Interactions on the Ion Permeation Mechanism in the hNaV1.5 Channel.

This study explores the ion selectivity and conduction mechanisms of the hNaV1.5 sodium channel using classical molecular dynamics simulations under an externally applied electric field. Our findings reveal distinct conduction mechanisms for Na+ and K+, primarily driven by differences in their hydration states when they diffuse close to the channel's selective filter (DEKA) and extracellular ring (EEDD). The Na+ ions undergo partial dehydration in the EEDD region, followed by a rehydration step in the DEKA ring, resulting in longer retention times and a deeper free energy minimum compared to K+. Conversely, the K+ ions exhibit a continuous dehydration process, facilitating a smoother passage through these key regions. These results indicate that ion selectivity and conductance are primarily governed by solvation dynamics, which, in turn, depend on the interactions with key charged residues within the channel. Additionally, we show that the delicate energetic balance between the interactions of the ions with the protein residues and with their solvation shells during the dehydration and rehydration processes is not properly captured by the force field. As a consequence, the selectivity of the channel is not well described, indicating that more accurate computational models must be applied to simulate ion conduction through eukaryotic NaV channels.

Multi-Site Red-Edge Excitation Shift Reveals the Residue-Specific Solvation Dynamics during the Native to Amyloid-like Transition of an Amyloidogenic Protein.

Changes in water-protein interactions are crucial for proteins to achieve functional and nonfunctional conformations during structural transitions by modulating local stability. Amyloid-like protein aggregates in deteriorating neurons are hallmarks of neurodegenerative disorders. These aggregates form through significant structural changes, transitioning from functional native conformations to supramolecular cross-β-sheet structures via misfolded and oligomeric intermediates in a multistep process. However, the site-specific dynamics of water molecules from the native to misfolded conformations and further to oligomeric and compact amyloid structures remain poorly understood. In this study, we used the fluorescence method known as red-edge excitation shift (REES) to investigate the solvation dynamics at specific sites in various equilibrium conformations en route to the misfolding and aggregation of the functional domain of the TDP-43 protein (TDP-43tRRM). We generated three single tryptophan-single cysteine mutants of TDP-43tRRM, with the cysteines at different positions and tryptophan at a fixed position. Each sole cysteine was fluorescently labeled and used as a site-specific fluorophore along with the single tryptophan, creating four monitorable sites for REES studies. By investigating the site-specific extent of REES, we developed a residue-specific solvation dynamics map of TDP-43tRRM during its misfolding and aggregation. Our observations revealed that solvation dynamics progressively became more rigid and heterogeneous to varying extents at different sites during the transition from native to amyloid-like conformations.

Solvation Dynamics and Microheterogeneity in Deep Eutectic Solvents.

Deep eutectic solvents have attracted considerable attention due to their unique properties and their potential to replace conventional solvents in diverse applications, such as catalysis, energy storage, and green chemistry. However, despite their broad use, the microscopic mechanisms governing solvation dynamics and the role of hydrogen bonding in deep eutectic solvents remain insufficiently understood. In this article, we present our contributions toward unravelling the micro heterogeneity within deep eutectic solvents by combining vibrational Stark spectroscopy and two-dimensional infrared spectroscopy with molecular dynamics simulations. Our findings demonstrate how the composition, constituents, and addition of water significantly influence the heterogeneous hydrogen bonding network and solvent dynamics within these systems. These insights provide valuable guidance for the design of next-generation solvents tailored to specific applications. By integrating experimental and computational approaches, this work sheds light on the intricate relationship between solvation dynamics and nanostructure in deep eutectic solvents, ultimately paving the way for innovative advances in solvent design.

Effect of Binary Solvent Mixtures on Ion-Pair Formation and Solvation Dynamics

This study determines the impact of binary solvent mixtures, particularly water and Dimethylformamide (DMF), on the formation of the ion-pair and dynamic solvation. The ion interaction dynamics of halides and nitrates with M²⁺ and M³⁺ have been assessed through conductivity studies and also ion-pair stability in addition to assessing solvation thermodynamics by varying the temperature. Experiments along with theoretical understanding lead the paper to highlight the critical modulatory effect of solvent polarity, hydrogen bonding, and viscosity of the medium affecting ion-ion interactions and finally, structures of solvation.

Study on the formation mechanism and effective manipulation of polymorphs and solvates in Osimertinib-Caffeic acid multi-component crystal with distinct properties.

Investigating the formation mechanism and effective manipulation of multi-component crystal polymorphs is crucial for facilitating industrial drug development. Herein, five novel Osimertinib-caffeic acid forms were first strategically tailored by varying solvent selection. Theoretical analysis demonstrated this polymorphism is correlated with multiple hydrogen bond donors-acceptors within multi-component system, which provides manipulation space for reconfiguration of intermolecular interactions and structural competition, while solvent further induced or involved in hydrogen-bonded rearrangements. Molecular dynamics simulation and solvent characterization revealed strong solute-solvent interaction and hydrogen bond donor propensity promoted the generation of hydrate. Small molecular size or large cavity of crystal facilitated solvent entry, resulting in the generation of channel-type solvates. Higher solvation free energy implied faster reaction rate and poorer solvent removal ability for solvates. Thermal analysis confirmed removal of the solvent occupying crystal channels for precursor solvates led to polymorph formation while maintaining the structure. Properties assessments revealed multi-component crystals significantly enhanced the physicochemical properties of Osimertinib. Polymorphs and hydrate exhibited remarkable distinctions in solubility, dissolution rate and physical stability. Nanoindentation and tableting tests confirmed distinct mechanical properties of different forms. Overall, this exploration has the potential to reshape the landscape of multi-component crystal polymorph development, paving the way for manipulating intermolecular interactions.

Bioinspired Molecular Electrets: Structure-Function Relationships When Dipoles Appear Everywhere

The importance of charge transfer (CT) and charge transport cannot be overstated (Phys. Chem. Chem. Phys. 2020, 22, 21583-21629). The shift to low-dielectric-constant carbon materials and devices, however, leads to the emergence of dominating electrostatic effects from localized dipole-generated fields. Resorting to molecular electrets allows examining the multifaceted nature of such effects. We develop bioinspired molecular electrets with enhanced CT capabilities that comprise anthranilamide residues. Even a single electret amino-acid residue can rectify CT (Angew. Chem. Int. Ed. 2018, 57, 12365-12369; J. Am. Chem. Soc. 2014, 136, 12966-12973). The molecular macrodipoles and the CT properties of these electrets strongly depend on their structure and conformational flexibility. Recently, we discovered an emergence of picosecond dipole transients (two-to-three times larger than the average dipole values) that originate from the solvent dynamics (J. Am. Chem. Soc. 2024, 146, 5162-5172). Further investigation reveals that these dipole transients appear universal and are not limited to large polar macromolecular structures. Such unexpectedly large dipole jumps provide underlying indication for the gating trends observed in picosecond and sub-picosecond CT kinetics (J. Am. Chem. Soc. 2016, 138, 12826-12832; Proc. Natl. Acad. Sci. USA 2021, 118, e2026462118). These unprecedented insights about the enormous picosecond dipole dynamics underlines new structure-function relationships essential for organic electronics and photonics, as well as for energy science and engineering.

(Invited) Towards AI2 Electrochemistry(AI2 = AI × ab initio)

It is known that electrode materials undergo dynamic structural changes at in-situ/in-operando conditions. Yet, the majority of computational studies only consider the static structures of electrode materials. When the materials are submerged in liquid solution, dynamic solvation effects are often completely ignored, or treated with dielectric continuum models, often lacking validation. The situations are about to change. Thanks to the latest development of in-situ experimental techniques and state-of-the-art computational methods, dynamics of electrode materials has recently drawn more and more attentions in many research areas. In this talk, I will present our recent progress on modeling dynamic catalysis and electrochemistry using ab initio molecular dynamics (AIMD). The high computational cost of AIMD however limits its application to small model systems consisting of hundreds of atoms at timescale of tens of ps. While, the latest development of AI accelerated AIMD (AI2MD) significantly increases the size and timescale, showing great promise for in situ modeling of realistic electrochemical systems.

Testing the prediction capability of trained models: sildenafil citrate solubility in propylene glycol and 1-propanol mixtures at different temperatures

ABSTRACT This study investigates the thermodynamic characteristics, equilibrium solubility, and solvation dynamics of sildenafil citrate within binary solvent systems composed of propylene glycol and 1-propanol. Utilizing both determination and computational approaches, the solubility patterns of sildenafil citrate were quantitatively assessed, revealing direct proportionality with increments in temperature and the concentration of propylene glycol. A comparative analysis using three distinct mathematical models to represent the solid–liquid phase equilibrium was conducted. Among these, the Jouyban-Acree-van’t Hoff model yielded the most comprehensive correlation by covering the effects of both temperature and solvent composition variations on the solubility, with a mean relative deviation (MRD%) of 6.4%, suggesting a robust alignment between the modelled and experimental solubility data. Previously trained models were tested to predict the solubility of the drug in mono- and mixed-solvents at various temperatures where accurate predictions were obtained with the MRD% of 19.1% (for mono-solvents), 8.1 and 14.4% (for mixed solvent systems).

Lagrangian formulation of nuclear-electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems.

We present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear-electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

The Role of Water Networks in Phosphodiesterase Inhibitor Dissociation and Kinetic Selectivity.

In search of new opportunities to develop Trypanosoma brucei phosphodiesterase B1 (TbrPDEB1) inhibitors that have selectivity over the off-target human PDE4 (hPDE4), different stages of a fragment-growing campaign were studied using a variety of biochemical, structural, thermodynamic, and kinetic binding assays. Remarkable differences in binding kinetics were identified and this kinetic selectivity was explored with computational methods, including molecular dynamics and interaction fingerprint analyses. These studies indicate that a key hydrogen bond between GlnQ.50 and the inhibitors is exposed to a water channel in TbrPDEB1, leading to fast unbinding. This water channel is not present in hPDE4, leading to inhibitors with a longer residence time. The computer-aided drug design protocols were applied to a recently disclosed TbrPDEB1 inhibitor with a different scaffold and our results confirm that shielding this key hydrogen bond through disruption of the water channel represents a viable design strategy to develop more selective inhibitors of TbrPDEB1. Our work shows how computational protocols can be used to understand the contribution of solvent dynamics to inhibitor binding, and our results can be applied in the design of selective inhibitors for homologous PDEs found in related parasites.

Superior Photostability of the Unnatural Base 6-Amino-5-nitropyridin-2-ol: A Case Study Using Ultrafast Broadband Fluorescence, Transient Absorption, and Theoretical Computation.

6-Amino-5-nitropyridin-2-ol (Z), a nitroaromatic compound and a base for Hachimoji nucleic acids, holds significant potential in expanding the genetic alphabet, as well as in synthetic biology and biotechnology. Despite its promising applications, the spectral characterization and photoinduced properties of Z have remained largely unexplored until now. This study presents a comprehensive investigation into its excited state dynamics in various solvents, utilizing state-of-the-art ultrafast broadband time-resolved fluorescence and transient absorption spectroscopy, complemented by computational methods. The acquired results provide direct experimental evidence that, upon photoexcitation, Z emits prompt fluorescence from a nearly planar structure in its excited state, independent of solvent properties. This state deactivates nonradiatively within sub-picoseconds through internal conversion with a unitary yield, primarily mediated by the rotation of the nitro group. This unusually rapid deactivation pathway entirely excludes the involvement of long-lived nπ* states, triplet states, and photoproducts, which are commonly observed in most nitroaromatic compounds and natural DNA and RNA bases. Our findings underscore that Z, as an unnatural base, exhibits superior photostability compared to canonical natural bases. This provides valuable insights into the photodynamics of nitroaromatic compounds, which is beneficial for strategic substitution design in environmental and biological applications.

Impact of ethanol on the flotation efficiency of imidazolium ionic liquids as collectors: Insights from dynamic surface tension and solvation analysis

Impact of ethanol on the flotation efficiency of imidazolium ionic liquids as collectors: Insights from dynamic surface tension and solvation analysis

Anomalous Water Penetration in Al3+ Dissolution.

The physicochemical characterization of trivalent ions is limited due to a lack of accurate force fields. By leveraging the latest machine learning force field to model aqueous AlCl3, we discover that upon dissolution of Al3+, water molecules beyond the second hydration shell are involved in the hydration process. A combination of scissoring of coordinating water is followed by synchronized secondary motion of water in the second solvation shell due to hydrogen bonding. Consequently, the water beyond the second solvation penetrates through the second solvation shell and coordinates to the Al3+. Our study reveals a novel microscopic understanding of solvation dynamics for the trivalent ion.

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R2) and a low mean squared error. Specifically, the achieved R2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. The present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.