Comparison of Different Three-site Interaction Potentials for Liquid Acetonitrile

Entropy Drives Accelerated Ion Diffusion upon Carbon Dioxide Expansion of Electrolytes

Carbon dioxide-expanded liquids, organic solvents with high concentrations of soluble carbon dioxide (CO2) at mild pressures, have gained attention as green catalytic media due to their improved properties over traditional solvents. More recently, carbon dioxide-expanded electrolytes (CXEs) have demonstrated improved reaction rates in the electrochemical reduction of CO2, by increasing the rate of delivery of CO2 to the electrode while maintaining facile charge transport. However, recent studies indicate that the limiting behavior of CXEs at higher CO2 pressures is a decline in solution conductivity due to reduced polarity, leading to poorer charge screening and greater ion pairing. In this article, we employ molecular dynamics simulations to investigate the energetic driving forces behind the diffusive properties of an acetonitrile and tetrapropylammonium hexafluorophosphate (TPrAPF6) CXE with increasing CO2 concentration. Our results indicate that entropy drives solvent and electrolyte diffusion with increasing CO2 pressure. The activation energy of ion diffusion increases with higher concentrations of CO2, indicating that increasing the temperature may improve solution conductivity in these systems. This trend in the activation energies is traced to stronger cation-anion Coulombic interactions due to weaker solvent screening at high CO2 concentrations, suggesting that the choice of ion may provide a route to diminish this effect.

Atomic-scale observation of solvent reorganization influencing photoinduced structural dynamics in a copper complex photosensitizer

Photochemical reactions in solution are governed by a complex interplay between transient intramolecular electronic and nuclear structural changes and accompanying solvent rearrangements. State-of-the-art time-resolved X-ray solution scattering has emerged in the last decade as a powerful technique to observe solute and solvent motions in real time. However, disentangling solute and solvent dynamics and how they mutually influence each other remains challenging. Here, we simultaneously measure femtosecond X-ray emission and scattering to track both the intramolecular and solvation structural dynamics following photoexcitation of a solvated copper photosensitizer. Quantitative analysis assisted by molecular dynamics simulations reveals a two-step ligand flattening strongly coupled to the solvent reorganization, which conventional optical methods could not discern. First, a ballistic flattening triggers coherent motions of surrounding acetonitrile molecules. In turn, the approach of acetonitrile molecules to the copper atom mediates the decay of intramolecular coherent vibrations and induces a further ligand flattening. These direct structural insights reveal that photoinduced solute and solvent motions can be intimately intertwined, explaining how the key initial steps of light harvesting are affected by the solvent on the atomic time and length scale. Ultimately, this work takes a step forward in understanding the microscopic mechanisms of the bidirectional influence between transient solvent reorganization and photoinduced solute structural dynamics.

Polarizable Force Field for Acetonitrile Based on the Single-Center Multipole Expansion

A polarizable potential function describing the interaction between acetonitrile molecules is introduced. The molecules are described as rigid and linear, with three mass sites corresponding to the CH₃ group (methyl, Me), the central carbon atom (C), and the nitrogen atom (N). The electrostatic interaction is represented using a single-center multipole expansion as has been done previously for H₂O [Wikfeldt et al., Phys. Chem. Chem. Phys. 15, 16542 (2013)], by including multipole moments from dipole up to and including hexadecapole, as well as anisotropic dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole polarizability tensors. The model is free of point charges. The non-electrostatic part is described in a pair-wise fashion by a Born-Mayer repulsion and damped dispersion attraction. The potential function is parameterized to fit the interaction energy of small (CH₃CN)_n, n = 2-6, clusters calculated using the PBE0 hybrid functional with an additional atomic many-body dispersion contribution. The parameterized potential function is found to compare well with results of the electronic structure calculations of dissociation curves for different dimer orientations and cohesive properties (the equilibrium volume, cohesive energy, and the bulk modulus) of the α-phase of acetonitrile crystal. The average value of the molecular dipole moment obtained in the α-phase is 5.53 D, corresponding to ca. 40% increase as compared to the dipole moment of an isolated acetonitrile molecule, 3.92 D. The calculated densities of solid and liquid acetonitrile turn out to be 8-10% higher than experimental values. This appears to be caused by an overestimate of the atomic many-body dispersion interaction in the density functional calculations used as input in the parametrization of the potential function.

Enhancing the Accuracy of Ab Initio Molecular Dynamics by Fine Tuning of Effective Two-Body Interactions: Acetonitrile as a Test Case

Grimme's dispersion-corrected density functional theory (DFT-D) methods have emerged among the most practical approaches to perform accurate quantum mechanical calculations on molecular systems ranging from small clusters to microscopic and mesoscopic samples, i.e., including hundreds or thousands of molecules. Moreover, DFT-D functionals can be easily integrated into popular ab initio molecular dynamics (MD) software packages to carry out first-principles condensed-phase simulations at an affordable computational cost. Here, starting from the well-established D3 version of the dispersion-correction term, we present a simple protocol to improve the accurate description of the intermolecular interactions of molecular clusters of growing size, considering acetonitrile as a test case. Optimization of the interaction energy was performed with reference to diffusion quantum Monte Carlo calculations, successfully reaching the same inherent accuracy of the latter (statistical error of ∼0.1 kcal/mol per molecule). The refined DFT-D3 model was then used to perform ab initio MD simulations of liquid acetonitrile, again showing significant improvements toward available experimental data with respect to the default correction.

Solvation of anthraquinone and TEMPO redox-active species in acetonitrile using a polarizable force field

Redox-active molecules are of interest in many fields, such as medicine, catalysis, or energy storage. In particular, in supercapacitor applications, they can be grafted to ionic liquids to form so-called biredox ionic liquids. To completely understand the structural and transport properties of such systems, an insight at the molecular scale is often required, but few force fields are developed ad hoc for these molecules. Moreover, they do not include polarization effects, which can lead to inaccurate solvation and dynamical properties. In this work, we developed polarizable force fields for redox-active species anthraquinone (AQ) and 2,2,6,6-tetra-methylpiperidinyl-1-oxyl (TEMPO) in their oxidized and reduced states as well as for acetonitrile. We validate the structural properties of AQ, AQ^•-, AQ^2-, TEMPO^•, and TEMPO⁺ in acetonitrile against density functional theory-based molecular dynamics simulations and we study the solvation of these redox molecules in acetonitrile. This work is a first step toward the characterization of the role played by AQ and TEMPO in electrochemical and catalytic devices.

Constrained Molecular Dynamic Simulation of the Potential Mean Force of Lithium Bromide Ion Pairs in Acetonitrile

Molecular dynamic simulations of Li+, and Br− ions in acetonitrile were carried out. The simulated structural properties were compared to experimental data. The solvent potentials of Li+-Br−, Li+-Li+, and Br−-Br− were evaluated using constrained molecular dynamics (CMD) simulations, to determine the solvent contribution to the total force acting on the solute and estimate the liquid arrangements according to the potential of mean force (PMF) values. The PMF of friction kernels and the passage across the Li+-Br− barrier was studied using the Grote–Hynes theory. The union-separation development happens in the polarization confining system.

On the interplay of solvent and conformational effects in simulated excited-state dynamics of a copper phenanthroline photosensitizer

QM/MM direct dynamics simulations in acetonitrile reveal the interplay between solvent and conformational effects in the photoinduced ultrafast flattening of a copper photosensitizer.

Molecular Insights into Dipole Relaxation Processes in Water-Lysine Mixtures

Dielectric spectroscopy is a robust method to investigate relaxations of molecular dipoles. It is particularly useful for studies of biological solutions because of the potential of this method to cover a broad range of dynamical time scales typical for such systems. However, this technique does not provide any information about the nature of the molecular motions, which leads to a certain underemployment of dielectric spectroscopy for gaining microscopic understanding of material properties. For such detailed understanding, computer simulations are valuable tools because they can provide information about the nature of molecular motions observed by, for example, dielectric spectroscopy and to further complement them with structural information. In this work, we acquire information about the nature of dipole relaxation, in n-lysine solutions by means of molecular dynamics simulations. Our results indicate that the experimentally observed main relaxation process of n-lysine has different origins for the single monomer and the polypeptide chains. The relaxation of 1-lysine is due to the motions of whole molecules, whereas the experimentally observed relaxation of 3-lysine and 4-lysine is due to the motions of the residues, which, in turn, are promoted by water relaxation. Furthermore, we propose a new structural model of the lysine amino acids, which can quantitatively account for the experimental dielectric relaxation data. Hydrogen bonding and the structure of water are also discussed in terms of their influence on relaxation processes.

Excited-state solvation structure of transition metal complexes from molecular dynamics simulations and assessment of partial atomic charge methods

Excited-state solvation structure (radial distribution function) of transition metal complexes by classical and mixed quantum-classical (QM/MM) molecular dynamics simulations.

Structures, intermolecular interactions, and chemical hardness of binary water-organic solvents: a molecular dynamics study

The evolution of structural properties, thermodynamics and averaged (dynamic) total hardness values as a function of the composition of binary water-organic solvents, was rationalized in view of the intermolecular interactions. The organic solvents considered were ethanol, acetonitrile, and isopropanol at 0.25, 0.5, 0.75, and 1 mass fractions, and the results were obtained using molecular dynamics simulations. The site-to-site radial distribution functions reveal a well-defined peak for the first coordination shell in all solvents. A characteristic peak of the second coordination shell exists in aqueous mixtures of acetonitrile, whereas in the water-alcohol solvents, a second peak develops with the increase in alcohol content. From the computed coordination numbers, averaged hydrogen bonds and their lifetimes, we found that water mixed with acetonitrile largely preserves its structural features and promotes the acetonitrile structuring. Both the water and alcohol structures in their mixtures are disturbed and form hydrogen bonds between molecules of different kinds. The dynamic hardness values are obtained as the average over the total hardness values of 1200 snapshots per solvent type, extracted from the equilibrium dynamics. The dynamic hardness profile has a non-linear evolution with the liquid compositions, similarly to the thermodynamic properties of these non-ideal solvents. Graphical abstract Computed dynamic total hardness, as a function of the cosolvent mass fraction for water-ethanol (EtOH), water-isopropanol (2PrOH) and water-acetonitrile (AN).

On the nature of solvatochromic effect: The riboflavin absorption spectrum as a case study

We present here the calculation of the absorption spectrum of riboflavin in acetonitrile and dimethyl sulfoxide using a hybrid quantum/classical approach, namely the perturbed matrix method, based on quantum mechanical calculations and molecular dynamics simulations. The calculated spectra are compared to the absorption spectrum of riboflavin previously calculated in water and to the experimental spectra obtained in all three solvents. The experimentally observed variations in the absorption spectra upon change of the solvent environment are well reproduced by the calculated spectra. In addition, the nature of the excited states of riboflavin interacting with different solvents is investigated, showing that environment effects determine a recombination of the gas-phase electronic states and that such a recombination is strongly affected by the polarity of the solvent inducing significant changes in the absorption spectra.

Solvent-induced conformational changes in cyclic peptides: a vibrational circular dichroism study

The three-dimensional structure of a peptide is strongly influenced by its solvent environment. In the present study, we study three cyclic tetrapeptides which serve as model peptides for β-turns. They are of the general structure cyclo(Boc-Cys-Pro-X-Cys-OMe) with the amino acid X being either glycine (1), or L- or D-leucine (L- or D-2). Using vibrational circular dichroism (VCD) spectroscopy, we confirm previous NMR results which showed that D-2 adopts predominantly a βII turn structure in apolar and polar solvents. Our results for L-2 indicate a preference for a βI structure over βII. With increasing solvent polarity, the preference for 1 is shifted from βII towards βI. This conformational change goes along with the breaking of an intramolecular hydrogen bond which stabilizes the βII conformation. Instead, a hydrogen bond with a solvent molecule can stabilize the βI turn conformation.

Dissimilar dynamics of coupled water vibrations

Dissimilar dynamics of coupled stretch vibrations of a water molecule are revealed by two-dimensional IR correlation spectroscopy. These are caused by essentially non-Gaussian fluctuations of the electric field exerted by the environment on the individual OH stretch vibrations. Non-Gaussian statistics of the individual site frequency fluctuations results in distinctively different dephasing of the symmetric and asymmetric eigenmodes. This phenomenon can only be described if the assumption of Gaussian dynamics in the traditional theories is abandoned.

Exploring the Dynamics of Calix[4]pyrrole: Effect of Solvent and Fluorine Substitution

Molecular dynamics simulations show that calix[4]pyrrole (CP) and octafluorocalix[4]pyrrole (8F-CP) are extremely flexible molecules. CP mainly adopts the 1,3-alternate conformation in all the solvents, although the percentage of alternative conformations increases in polar solvents, especially those with good hydrogen-bonding acceptor properties. However, in the case of 8F-CP, the cone conformation is the most populated in some solvents. Transitions between conformers are common and fast, and both CP and 8F-CP can adopt the cone conformation needed for optimum interaction with anions more easily than would be predicted on the basis of previous gas-phase calculations. Furthermore, the present studies show that when a fluoride anion is specifically placed initially in close proximity to CP and 8F-CP in their respective 1,3-alternate conformations, an extremely fast change to the cone conformation is observed in both cases. The results suggest that preorganization does not represent a major impediment to anion-binding for either CP or 8F-CP, and that ion-induced conformational changes can follow different mechanisms depending on the solvent and the chemical substituents present on the calix[4]pyrrole beta-pyrrolic positions.

New six-site acetonitrile model for simulations of liquid acetonitrile and its aqueous mixtures

A new six site flexible acetonitrile molecular model is developed. The AMBER force field was used for description of intramolecular parameters, the atomic charges were calculated from a high level ab initio theory and finally the Lennard-Jones parameters were tuned to fit the experimental density and evaporation heat. The obtained in this way model reproduces correctly densities of water-acetonitrile mixtures as well as provides qualitative description of the dielectric permittivity and self-diffusion coefficients.

1H and 13C-NMR and molecular dynamics studies of cyclosporin a interacting with magnesium(II) or cerium(III) in acetonitrile. Conformational changes and cis-trans conversion of peptide bonds

Cyclosporin A (CsA) is an important drug used to prevent graft rejection in organ transplantations. Its immunosuppressive activity is related to the inhibition of T-cell activation through binding with the proteins Cyclophilin (Cyp) and, subsequently, Calcineurin (CN). In the complex with its target (Cyp), CsA adopts a conformation with all trans peptide bonds and this feature is very important for its pharmacological action. Unfortunately, CsA can cause several side effects, and it can favor the excretion of calcium and magnesium. To evaluate the possible role of conformational effects induced by these two metal ions in the action mechanism of CsA, its complexes with Mg(II) and Ce(III) (the latter as a paramagnetic probe for calcium) have been examined by two-dimensional NMR and relaxation rate analysis. The conformations of the two complexes and of the free form have been determined by restrained molecular dynamics calculations based on the experimentally obtained metal-proton and interproton distances. The findings here ratify the formation of 1:1 complexes of CsA with both Mg(II) and Ce(III), with metal coordination taking place on carbonyl oxygens and substantially altering the peptide structure with respect to the free form, although the residues involved and the resulting conformational changes, including cis-trans conversion of peptide bonds, are different for the two metals.

A transferable electrostatic map for solvation effects on amide I vibrations and its application to linear and two-dimensional spectroscopy

A method for modeling infrared solvent shifts using the electrostatic field generated by the solvent is presented. The method is applied to the amide I vibration of N-methyl acetamide. Using ab initio calculations the fundamental frequency, anharmonicity, and the transition dipoles between the three lowest vibrational states are parametrized in terms of the electrostatic field. The generated map, which takes into account the electric field and its gradients at four molecular positions, is tested in a number of common solvents. Agreement of solvent shift and linewidths with experimental Fourier transform infrared (FTIR) data is found to within seven and four wave numbers, respectively, for polar solvents. This shows that in these solvents electrostatic contributions dominate solvation effects and the map is transferable between these types of solvents. The effect of motional narrowing arising from the fast solvent fluctuations is found to be substantial for the FTIR spectra. Also the two-dimensional infrared (2DIR) spectra, simulated using the constructed map, reproduce experimental results very well. The effect of anharmonicity fluctuations on the 2DIR spectra was found to be negligible.

Effect of the Range of Interactions on the Properties of Fluids. 2. Structure and Phase Behavior of Acetonitrile, Hydrogen Fluoride, and Formic Acid

To complete the study on the effect of the long-range part of Coulombic interactions on properties of complex polar and associating fluids, we have investigated in detail three compounds with extreme features: acetonitrile for its unusually large dipole moment, hydrogen fluoride with very strong hydrogen bonding, and formic acid for its potential formation of different n-mers in liquid and gaseous phases. The effect of the long-range Coulombic interactions on both the structure and thermodynamics of the homogeneous phase, and on the vapor-liquid equilibria has been examined using the same decomposition of realistic potential models into a short-range part and a residual part as in the previous paper [Kettler, M.; et al. J. Phys. Chem. B 2002, 106, 7537-7546]. The present results fully confirm the previous findings that the properties of polar and associating systems are determined primarily by the short-range interactions regardless of their nature, i.e., contributions arising from the long-range interactions constitute only a small portion of the total properties, and thus that the short-range potential counterpart of full realistic models can be used as a convenient reference for a successful perturbation expansion.

Collective contributions to the dielectric relaxation of hydrogen-bonded liquids

Dielectric relaxation times are often interpreted in terms of the reorientation of dipolar species or aggregates. The relevant time correlation function contains, however, cross terms between dipole moments of different particles. In the static case, these cross terms are accounted for by the Kirkwood factor g(K). Theories and molecular dynamics simulations suggest that such cross correlations may also affect the time-dependent properties, as reflected in the dielectric spectra. We present an experimental method for detecting effects of such cross correlations in dielectric spectra by a comparative analysis of dielectric and magnetic relaxation data. We demonstrate that such collective contributions can substantially affect dielectric relaxation. Experiments for n-pentanol (g(K)=3.06 at 298 K) and 2,2-dimethyl-3-ethyl-pentane-3-ol (g(K)=0.59) and their solutions in carbon tetrachloride show that in systems with g(K)>1, the cross correlations slow down dielectric relaxation. In systems with g(K)<1, dielectric relaxation is enhanced. The results conform to theoretical predictions by Madden and Kivelson [Adv. Chem. Phys. 56, 467 (1984)] and to results of molecular dynamics simulations. The relaxation enhancement by cross terms in the case of g(K)<1 is difficult to rationalize by conventional models of dielectric relaxation.

Molecular dynamics simulations of acetonitrile/dimethyl sulfoxide liquid mixtures

Binary liquid mixtures of dimethyl sulfoxide and acetonitrile at the three molar fractions 0.25, 0.50, and 0.75 have been investigated by molecular dynamics computer simulations. Thermodynamic states corresponding to liquid-vapor coexistence at a temperature of 298 K were considered. Intermolecular interactions were described by potential models of the site-site (12-6) Lennard-Jones plus Coulomb type that have been developed for the description of the pure liquids. Dimethyl sulfoxide has been represented by four interactions sites and acetonitrile by a three- as well as a six-site potential model. We have calculated thermodynamic properties and the intermolecular pair distribution functions. Intermolecular interaction energies indicate deviations from the behavior of ideal mixtures. The local mole fraction analysis demonstrates that dimethyl sulfoxide is preferentially solvated by acetonitrile and that the first solvation shell surrounding acetonitrile molecules is significantly enriched by dimethyl sulfoxide. The nonideal behavior in the mixtures is not affected by the choice of the three- or the six-site potential model for acetonitrile. Orientational correlations of dipole vectors within the first solvation shells indicate that the relative molecular orientations found in pure acetonitrile and dimethyl sulfoxide are maintained in the mixtures. Parallel and antiparallel dipole-dipole configurations determine first shell acetonitrile-dimethyl sulfoxide configurations. Dynamical features of the mixtures are discussed in terms of diffusion constants and orientational correlation times as obtained from the time correlation functions for linear velocities and molecular dipole moments, respectively. Computed relaxation times indicate faster reorientational motion for dimethyl sulfoxide if acetonitrile is added. In contrast, the orientational dynamics of acetonitrile becomes stronger correlated upon dilution with dimethyl sulfoxide. The diffusion coefficients for both compounds follow this tendency.

Theoretical study of anion binding to calix[4]pyrrole: the effects of solvent, fluorine substitution, cosolute, and water traces

The binding of different anions to calix[4]pyrrole has been studied by means of molecular dynamics coupled to thermodynamic integration calculations. The effect of different apolar solvents, octafluoro substitution, and the change in binding free energy derived from the presence of cosolute and water traces (the hydrated salt used to introduce the anion in the solution) were examined. Calculations allow us to rationalize the differential binding of ions to calix[4]pyrrole and octafluorocalix[4]pyrrole as well as to predict the behavior in new solvents for which experimental data are not available yet. It is found that both calix[4]pyrrole and octafluorocalix[4]pyrrole have a dramatic preference for F- in the gas phase and pure aprotic solvents, but the situation can change dramatically in protic solvents or in the presence of the hydrated cation which is used as cosolute of the anion. Overall, our results provide interesting clues for a better understanding of the process detected experimentally as "binding".