1
|
Okita K, Ito N, Morishita-Watanabe N, Umakoshi H, Kasahara K, Matubayasi N. Solvation dynamics on the diffusion timescale elucidated using energy-represented dynamics theory. Phys Chem Chem Phys 2024; 26:12852-12861. [PMID: 38623745 DOI: 10.1039/d4cp00235k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
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.
Collapse
Affiliation(s)
- Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Natsuumi Ito
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Nozomi Morishita-Watanabe
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Hiroshi Umakoshi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| |
Collapse
|
2
|
Okita K, Kasahara K, Matubayasi N. Diffusion theory of molecular liquids in the energy representation and application to solvation dynamics. J Chem Phys 2022; 157:244505. [PMID: 36586971 DOI: 10.1063/5.0125432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
The generalized Langevin equation (GLE) formalism is a useful theoretical fundament for analyzing dynamical phenomena rigorously. Despite the systematic formulation of dynamics theories with practical approximations, however, the applicability of GLE-based methods is still limited to simple polyatomic liquids due to the approximate treatment of molecular orientations involved in the static molecular liquid theory. Here, we propose an exact framework of dynamics based on the GLE formalism incorporating the energy representation theory of solution, an alternative static molecular liquid theory. A fundamental idea is the projection of the relative positions and orientations of solvents around a solute onto the solute-solvent interaction, namely the energy coordinate, enabling us to describe the dynamics on a one-dimensional coordinate. Introducing systematic approximations, such as the overdamped limit, leads to the molecular diffusion equation in the energy representation that is described in terms of the distribution function of solvents on the energy coordinate and the diffusion coefficients. The present theory is applied to the solvation dynamics triggered by the photoexcitation of benzonitrile. The long-time behavior of the solvation time correlation function is in good agreement with that obtained by the molecular dynamics simulation.
Collapse
Affiliation(s)
- Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
3
|
Yamaguchi T, Yoshida N. Solvation dynamics in electronically polarizable solvents: Theoretical treatment using solvent-polarizable three-dimensional reference interaction-site model theory combined with time-dependent density functional theory. J Chem Phys 2021; 154:044504. [PMID: 33514097 DOI: 10.1063/5.0036289] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The theory of solvation structure in an electronically polarizable solvent recently proposed by us, referred to as the "solvent-polarizable three-dimensional reference interaction-site model theory," is extended to dynamics in this study through the combination with time-dependent density functional theory. Test calculations are performed on model charge-transfer systems in water, and the effects of electronic polarizability on solvation dynamics are examined. The electronic polarizability slightly retards the solvation dynamics. This is ascribed to the decrease in the curvature of the nonequilibrium free energy profile along the solvation coordinate. The solvent relaxation is bimodal, and the faster and the slower modes are assigned to the reorientational and the translational modes, respectively, as was already reported by the surrogate theory combined with the site-site Smoluchowski-Vlasov equation. The relaxation path along the solvation coordinate is a little higher than the minimum free energy path because the translational mode is fixed in the time scale of the reorientational relaxation.
Collapse
Affiliation(s)
- Tsuyoshi Yamaguchi
- Graduate School of Engineering, Nagoya University, Chikusa, Nagoya 464-8603, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishiku, Fukuoka 819-0395, Japan
| |
Collapse
|
4
|
Kasahara K, Sato H. Dynamics theory for molecular liquids based on an interaction site model. Phys Chem Chem Phys 2017; 19:27917-27929. [DOI: 10.1039/c7cp05423h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Dynamics theories for molecular liquids based on an interaction site model have been developed over the past few decades and proved to be powerful tools to investigate various dynamical phenomena.
Collapse
Affiliation(s)
- Kento Kasahara
- Department of Molecular Engineering
- Kyoto University
- Japan
| | - Hirofumi Sato
- Department of Molecular Engineering and Elements Strategy for Catalysts and Batteries (ESICB)
- Kyoto University
- Japan
| |
Collapse
|
5
|
Masuda Y, Nakano T, Sugiyama M. First observation of ultrafast intramolecular proton transfer rate between electronic ground states in solution. J Phys Chem A 2012; 116:4485-94. [PMID: 22510164 DOI: 10.1021/jp2110874] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Despite the importance of ultrafast (time scale exceeding 10(-11) s) intramolecular proton transfer (PT) events between electronic ground states in solution, experimental determination of the rates of such reactions has not yet been accomplished because of the limitations of the utilized methods. The objective of this study was to evaluate the PT rates of intramolecular O···H···O hydrogen-bonded systems in solution through the (1)H spin-lattice relaxation times of the hydroxyl protons, induced by the (1)H-(17)O dipolar interactions (T(1dd)(OH)), taking into account the contribution of the OH reorientational motion to T(1dd)(OH). Solutions of the benzoic acid dimer (BA dimer), 1-benzoyl-6-hydroxy-6-phenylfulvene (Fulvene), and dibenzoylmethane (DBM) were chosen as test systems. For Fulvene in CCl(4), the PT time, τ(PT), was deduced to be 7 × 10(-11) s. In the case of the BA dimer in CCl(4), the τ(PT) value was considerably greater than the OH reorientational correlation time, τ(R(OH)) = 4.3 × 10(-11) s. In contrast, the experimental results for DBM in CCl(4) indicated that the proton is located about midway between the two oxygen atoms, that is, the PT potential energy surface is a single well or a double well with a PT barrier near or below the zero-point energy.
Collapse
Affiliation(s)
- Yuichi Masuda
- Department of Chemistry, Faculty of Science, Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan.
| | | | | |
Collapse
|
6
|
Nishiyama K, Takata K, Watanabe K, Shigematsu H. Solvation dynamics of coumarin 153 embedded in AOT+phenol organogels studied by time-resolved fluorescence spectroscopy. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.01.050] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
7
|
Aono S, Yamamoto T, Kato S. Solution reaction space Hamiltonian based on an electrostatic potential representation of solvent dynamics. J Chem Phys 2011; 134:144108. [DOI: 10.1063/1.3572057] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
8
|
Nishiyama K, Yamaguchi T, Hirata F. Solvation Dynamics in Polar Solvents Studied by Means of RISM/Mode-Coupling Theory. J Phys Chem B 2009; 113:2800-4. [DOI: 10.1021/jp809926g] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Katsura Nishiyama
- Faculty of Education, Shimane University, Matsue, Shimane 690-8504, Japan, Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Nagoya, Aichi 464-8603, Japan, and Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| | - Tsuyoshi Yamaguchi
- Faculty of Education, Shimane University, Matsue, Shimane 690-8504, Japan, Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Nagoya, Aichi 464-8603, Japan, and Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| | - Fumio Hirata
- Faculty of Education, Shimane University, Matsue, Shimane 690-8504, Japan, Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Nagoya, Aichi 464-8603, Japan, and Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| |
Collapse
|
9
|
Megerle U, Selmaier F, Lambert C, Riedle E, Lochbrunner S. Symmetry-dependent solvation of donor-substituted triarylboranes. Phys Chem Chem Phys 2008; 10:6245-51. [DOI: 10.1039/b806131a] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
10
|
Kobryn AE, Yamaguchi T, Hirata F. Site-site memory equation approach in study of density/pressure dependence of translational diffusion coefficient and rotational relaxation time of polar molecular solutions: acetonitrile in water, methanol in water, and methanol in acetonitrile. J Chem Phys 2007; 122:184511. [PMID: 15918733 DOI: 10.1063/1.1884998] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present results of the theoretical study and numerical calculation of the dynamics of molecular liquids based on the combination of the memory equation formalism and the reference interaction site model (RISM). Memory equations for the site-site intermediate scattering functions are studied in the mode-coupling approximation for the first-order memory kernels, while equilibrium properties such as site-site static structure factors are deduced from RISM. The results include the temperature-density (pressure) dependence of translational diffusion coefficients D and orientational relaxation times tau for acetonitrile in water, methanol in water, and methanol in acetonitrile--all in the limit of infinite dilution. Calculations are performed over the range of temperatures and densities employing the extended simple point charge model for water and optimized site-site potentials for acetonitrile and methanol. The theory is able to reproduce qualitatively all main features of temperature and density dependences of D and tau observed in real and computer experiments. In particular, anomalous behavior, i.e, the increase in mobility with density, is observed for D and tau of methanol in water, while acetonitrile in water and methanol in acetonitrile do not show deviations from the ordinary behavior. The variety exhibited by the different solute-solvent systems in the density dependence of the mobility is interpreted in terms of the two competing origins of friction, which interplay with each other as density increases: the collisional and dielectric frictions which, respectively, increase and decrease with increasing density.
Collapse
Affiliation(s)
- Alexander E Kobryn
- Department of Theoretical Study, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | | | | |
Collapse
|
11
|
Masuda Y, Shimizu C. Solvent Effect on Intramolecular Electron Transfer Rates of Mixed-Valence Biferrocene Monocation Derivatives. J Phys Chem A 2006; 110:7019-27. [PMID: 16737249 DOI: 10.1021/jp060087w] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Intramolecular electron transfer (ET) rates in various solvents of mixed-valence biferrocene monocation (Fe(II), Fe(III)) and the 1',1' ''-diiodo and 1',1' ''-diethyl derivatives (respectively abbreviated as BFC(+), I(2)BFC(+), and Et(2)BFC(+)) were determined by means of the spin-lattice relaxation times of the protons, taking into account the local magnetic field fluctuation caused by the electron hopping between the two ferrocene units. We also determined the ET rates of a mixed-valence diferrocenylacetylene monocation (DFA(+)) in order to examine the effect of the insertion of an acetylene bridge between the two ferrocene units. The insertion of the bridge decreased the ET rate, while the effect of substitution on the cyclopentadienyl rings on the rate was minor. The observed rates for each mixed-valence monocation in various solvents did not correlate with the reorganization energies, but we did find a significant contribution of the solvent dynamics. The observed rates were considerably higher than those expected on the basis of the Sumi-Marcus-Nalder model in which the solvents were regarded as dielectric continua. The slope of the logarithm plot of the pre-exponential factors in various solvents for each mixed-valence monocation versus the inverse of the longitudinal dielectric relaxation times of the solvents was significantly smaller than unity, and the slope for DFA(+) was larger than those for BFC(+), I(2)BFC(+), and Et(2)BFC(+). These results were ascribed to a partial contribution of the dielectric friction to the dynamics along the solvent coordinate; the extent of the contribution decreased with a reduction in the ET distance. For the dynamics along the solvent coordinate of the ET reactions in methanol, the observed rates indicated an important contribution by the minor dielectric relaxation components with faster relaxation times, rather than the major component with an extraordinarily long relaxation time.
Collapse
Affiliation(s)
- Y Masuda
- Department of Chemistry, Faculty of Science, Ochanomizu University, Bunkyo-ku, Tokyo.
| | | |
Collapse
|
12
|
Yamaguchi T, Matsuoka T, Koda S. Generalized Langevin theory on the dynamics of simple fluids under external fields. J Chem Phys 2005; 123:34504. [PMID: 16080741 DOI: 10.1063/1.1955455] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.
Collapse
Affiliation(s)
- T Yamaguchi
- Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Aichi, Japan.
| | | | | |
Collapse
|
13
|
Nishiyama K, Yamaguchi T, Hirata F, Okada T. Solvation dynamics in water investigated by RISM/mode-coupling theory. J Mol Liq 2005. [DOI: 10.1016/j.molliq.2004.10.009] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
14
|
Pressure dependence of diffusion coefficient and orientational relaxation time for acetonitrile and methanol in water: RISM/mode-coupling study. J Mol Liq 2005. [DOI: 10.1016/j.molliq.2004.10.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
15
|
González-Mozuelos P, Yeom MS, Olvera de la Cruz M. Molecular multivalent electrolytes: microstructure and screening lengths. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2005; 16:167-78. [PMID: 15729508 DOI: 10.1140/epje/e2005-00018-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2004] [Indexed: 05/23/2023]
Abstract
We study small rod-like molecular electrolytes solutions with their corresponding atomic counterions. The asymptotic length scales (decay length and wavelength) of the structural correlations are analyzed using the formalism of the dressed interaction site theory (DIST). The correlation functions are determined using the reference interaction site model equation complemented with a mixed approach in which the hypernetted-chain closure is used for the repulsive interactions, and the mean spherical approximation is used for the attractive interactions. The results from this scheme are in good agreement with the Monte Carlo computer simulations reported here. The asymptotic properties of the correlation functions of this molecular system are compared against those corresponding to two related simple (atomic) electrolyte models. The main conclusion is that the molecular structure of the ions lowers by two orders of magnitude the concentration at which the transition from monotonic to oscillatory decay occurs.
Collapse
Affiliation(s)
- P González-Mozuelos
- Departamento de Física, Centro de Investigación y de Estudios Avanzados del IPN, A.P. 14-740, 07000 México D.F., Mexico
| | | | | |
Collapse
|
16
|
Yamaguchi T, Matsuoka T, Koda S. Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute. J Chem Phys 2005; 122:14512. [PMID: 15638679 DOI: 10.1063/1.1828039] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.
Collapse
Affiliation(s)
- T Yamaguchi
- Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Nagoya, Aichi 464-8603, Japan.
| | | | | |
Collapse
|
17
|
Bedard-Hearn MJ, Larsen RE, Schwartz BJ. Hidden Breakdown of Linear Response: Projections of Molecular Motions in Nonequilibrium Simulations of Solvation Dynamics. J Phys Chem A 2003. [DOI: 10.1021/jp034383+] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|