This computational study characterises charge-transfer-to-solvent (CTTS) states of aqueous thiocyanate anion using equation-of-motion coupled-cluster methods combined with electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. Equilibrium sampling was carried out using classical molecular dynamics (MD) with standard force-fields and QM/MM ab initio molecular dynamics (AIMD) using density functional theory. The two calculations yield significantly different local structure around solvated SCN−
. Because of the diffuse character of CTTS states, they are very sensitive to the local structure of solvent around the solute and its dynamic fluctuations. Owing to this sensitivity, the spectra computed using MD and AIMD based snapshots differ considerably. This sensitivity suggests that the spectroscopy exploiting CTTS transitions can provide an experimental handle for assessing the quality of force-fields and density functionals. By combining CTTS-based spectroscopies with reliable theoretical modeling, detailed microscopic information of the solvent structure can be obtained. We present a robust computational protocol for modeling spectra of solvated anions and emphasise the use of an ab initio characterization of individual electronic transitions as CTTS or local excitations.
more »
« less
A generalized Rayleigh-Plesset equation for ions with solvent fluctuations
We introduce a mathematical modeling framework for the conformational dynamics of charged molecules (i.e., solutes) in an aqueous solvent (i.e., water or salted water). The solvent is treated as an incompressible fluid, and its fluctuating motion is described by the Stokes equation with the Landau–Lifschitz stochastic stress. The motion of the solute-solvent interface (i.e., the dielectric boundary) is determined by the fluid velocity together with the balance of the viscous force,hydrostatic pressure, surface tension, solute-solvent van der Waals interaction force, and electrostatic force. The electrostatic interactions are described by the dielectric Poisson–Boltzmann theory.Within such a framework, we derive a generalized Rayleigh–Plesset equation, a nonlinear stochastic ordinary differential equation (SODE), for the radius of a spherical charged molecule, such as anion. The spherical average of the stochastic stress leads to a multiplicative noise. We design and test numerical methods for solving the SODE and use the equation, together with explicit solvent molecular dynamics simulations, to study the effective radius of a single ion. Potentially, our general modeling framework can be used to efficiently determine the solute-solvent interfacial structures and predict the free energies of more complex molecular systems.
more »
« less
- Award ID(s):
- 1913144
- NSF-PAR ID:
- 10288376
- Date Published:
- Journal Name:
- SIAM journal on applied mathematics
- Volume:
- 81
- Issue:
- 3
- ISSN:
- 0036-1399
- Page Range / eLocation ID:
- 1098-1115
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)To help better interpret experimental measurement of nanoparticle size, it is important to understand how their diffusion depends on the physical and chemical features of surface ligands. In this study, explicit solvent molecular dynamics simulations are used to probe the effect of ligand charge and flexibility on the diffusion of small gold nanoparticles. The results suggest that despite a high bare charge (+18 e), cationic nanoparticles studied here have reduced diffusion constants compared to a hydrophobic gold nanoparticle by merely a modest amount. Increasing the ligand length by 10 CH2 units also has a limited impact on the diffusion constant. For the three particles studied here, the difference between estimated hydrodynamic radius and radius of gyration is on the order of one solvent layer (3–5 Å), confirming that the significant discrepancies found in the size of similar nanoparticles by recent transmission electron microscopy and dynamic light scattering measurements were due to aggregation under solution conditions. The limited impact of electrostatic friction on the diffusion of highly charged nanoparticles is found to be due to the strong anticorrelation between electrostatic and van der Waals forces between nanoparticle and environment, supporting the generality of recent observation for proteins by Matyushov and co-workers. Including the first shell of solvent molecules as part of the diffusing particle has a minor impact on the total force autocorrelation function but reduces the disparity in relaxation time between the total force and its electrostatic and van der Waals components.more » « less
-
Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O–H bonds pointing toward the solute. In contrast to the well accepted structure–function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute–solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.more » « less
-
null (Ed.)The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Potential applications and caveats to the usage of LD are included as a summary.more » « less
-
Electrolyte solutions play an important role in energy storage devices, whose performance relies heavily on the electrokinetic processes at sub-micron scales. Although fluctuations and stochastic features become more critical at small scales, the long-range Coulomb interactions pose a particular challenge for both theoretical analysis and simulation of fluid systems with fluctuating hydrodynamic and electrostatic interactions. Here, we present a theoretical framework based on the Landau–Lifshitz theory to derive closed-form expressions for fluctuation correlations in electrolyte solutions, indicating significantly different decorrelation processes of ionic concentration fluctuations from hydrodynamic fluctuations, which provides insights for understanding transport phenomena of coupled fluctuating hydrodynamics and electrokinetics. Furthermore, we simulate fluctuating electrokinetic systems using both molecular dynamics (MD) with explicit ions and mesoscopic charged dissipative particle dynamics (cDPD) with semi-implicit ions, from which we identify that the spatial probability density functions of local charge density follow a gamma distribution at sub-nanometre scale (i.e. $0.3\,{\rm nm}$ ) and converge to a Gaussian distribution above nanometre scales (i.e. $1.55\,{\rm nm}$ ), indicating the existence of a lower limit of length scale for mesoscale models using Gaussian fluctuations. The temporal correlation functions of both hydrodynamic and electrokinetic fluctuations are computed from all-atom MD and mesoscale cDPD simulations, showing good agreement with the theoretical predictions based on the linearized fluctuating hydrodynamics theory.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
