Solving Chemical Absorption Equilibria using Free Energy and Quantum Chemistry Calculations: Methodology, Limitations, and New Open-Source Software.

We developed an open-source chemical reaction equilibrium solver in Python (CASpy, https://github.com/omoultosEthTuDelft/CASpy) to compute the concentration of species in any reactive liquid-phase absorption system. We derived an expression for a mole fraction-based equilibrium constant as a function of excess chemical potential, standard ideal gas chemical potential, temperature, and volume.

Vapor pressures and vapor phase compositions of choline chloride urea and choline chloride ethylene glycol deep eutectic solvents from molecular simulation.

Despite the widespread acknowledgment that deep eutectic solvents (DESs) have negligible vapor pressures, very few studies in which the vapor pressures of these solvents are measured or computed are available. Similarly, the vapor phase composition is known for only a few DESs. In this study, for the first time, the vapor pressures and vapor phase compositions of choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs are computed using Monte Carlo simulations.

New Features of the Open Source Monte Carlo Software Brick-CFCMC: Thermodynamic Integration and Hybrid Trial Moves.

We present several new major features added to the Monte Carlo (MC) simulation code Brick-CFCMC for phase- and reaction equilibria calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one is thermodynamic integration for the computation of excess chemical potentials (μ).

Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution.

We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li, NH , F, and Cl, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules.

Solvation of the Guanidinium Ion in Pure Aqueous Environments: A Theoretical Study from an "Ab Initio"-Based Polarizable Force Field.

We report simulation results regarding the hydration process of the guanidinium cation in water droplets and in bulk liquid water, at a low concentration of 0.03 M, performed using a polarizable approach to model both water/water and ion/water interactions. In line with earlier theoretical studies, our simulations show a preferential orientation of guanidinium at water-vacuum interfaces, i.

Organic ion association in aqueous phase and ab initio-based force fields: The case of carboxylate/ammonium salts.

We performed molecular dynamics simulations of carboxylate/methylated ammonium ion pairs solvated in bulk water and of carboxylate/methylated ammonium salt solutions at ambient conditions using an ab initio-based polarizable force field whose parameters are assigned to reproduce only high end quantum computations, at the Møller-Plesset second-order perturbation theory/complete basis set limit level, regarding single ions and ion pairs as isolated and micro-hydrated in gas phase. Our results agree with the available experimental results regarding carboxylate/ammonium salt solutions. For instance, our force field approach predicts the percentage of acetate associated with ammonium ions in CHCOO/CHNH solutions at the 0.

Extrapolating Single Organic Ion Solvation Thermochemistry from Simulated Water Nanodroplets.

We compute the ion/water interaction energies of methylated ammonium cations and alkylated carboxylate anions solvated in large nanodroplets of 10 000 water molecules using 10 ns molecular dynamics simulations and an all-atom polarizable force-field approach. Together with our earlier results concerning the solvation of these organic ions in nanodroplets whose molecular sizes range from 50 to 1000, these new data allow us to discuss the reliability of extrapolating absolute single-ion bulk solvation energies from small ion/water droplets using common power-law functions of cluster size. We show that reliable estimates of these energies can be extrapolated from a small data set comprising the results of three droplets whose sizes are between 100 and 1000 using a basic power-law function of droplet size.

Simulated Solvation of Organic Ions II: Study of Linear Alkylated Carboxylate Ions in Water Nanodrops and in Liquid Water. Propensity for Air/Water Interface and Convergence to Bulk Solvation Properties.

We investigated the solvation of carboxylate ions from formate to hexanoate, in droplets of 50 to 1000 water molecules and neat water, by computations using standard molecular dynamics and sophisticated polarizable models. The carboxylate ions from methanoate to hexanoate show strong propensity for the air/water interface in small droplets. Only the ions larger than propanoate retain propensity for the interface in larger droplets, where their enthalpic stabilization by ion/water dispersion is reduced there by 3 kcal mol(-1) per CH2 group.

Physisorbed H2@Cu(100) surface: potential and spectroscopy.

Using an embedding approach, a 2-D potential energy function has been calculated to describe the physisorption interaction of H2 with a Cu(100) surface. For this purpose, a cluster model of the system calculated with highly correlated wavefunctions is combined with a periodic Density-Functional-Theory method using van der Waals-DF2 functional. Rotational and vibrational energy levels of physisorbed H2, as well as D2 and HD, are calculated using the 2D embedding corrected potential energy function.

Simulated solvation of organic ions: protonated methylamines in water nanodroplets. Convergence toward bulk properties and the absolute proton solvation enthalpy.

We applied an alternative, purely theoretical route to estimate thermodynamical properties of organic ions in bulk solution. The method performs a large ensemble of simulations of ions solvated in water nanodroplets of different sizes, using a polarizable molecular dynamics approach. We consider protonated ammonia and methylamines, and K(+) for comparison, solvated in droplets of 50-1000 water molecules.

Structural and atoms-in-molecules analysis of hydrogen-bond network around nitroxides in liquid water.

In this study, we investigated the hydrogen-bond network patterns involving the NO moieties of five small nitroxides in liquid water by analyzing nanosecond scale molecular dynamics trajectories. To this end, we implemented two types of hydrogen-bond definitions, based on electronic structure, using Bader's atoms-in-molecules analysis and based on geometric criteria. In each definition framework, the nitroxide/water hydrogen-bond networks appear very variable from a nitroxide to another.

Structure and spectromagnetic properties of the superoxide radical adduct of DMPO in water: elucidation by theoretical investigations.

In the field of spin trapping chemistry, the design of more efficient radical traps can be assisted by the development of theoretical methods able to give a quantitative evaluation of the electron paramagnetic resonance (EPR) spectrum features of the spin-adduct radical, even before initiating the experimental work. The superoxide radical adduct of the 5,5-dimethyl-1-pyrroline-N-oxide nitrone (DMPO-OOH) has been reported in a huge number of papers devoted to the study of the oxidative stress. Here, we present for the first time the theoretical study of DMPO-OOH in an explicit water solution, based on the combined QM/MM//MD protocol we recently proposed, featuring a full coupling between the solute and all the explicit water molecules.

Viscosity of the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquid from equilibrium and nonequilibrium molecular dynamics.

This work deals with the determination, using equilibrium and nonequilibrium molecular dynamics, of the viscosity of an ionic liquid: 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([emim][Tf2N]). A first method consists in computing the shear viscosity using the Green-Kubo formalism from the pressure tensor correlation function obtained from equilibrium simulations. On the other hand, the Newtonian viscosity can be extracted from the shear rate dependence of the viscosity in nonequilibrium simulations using the mode-coupling theory and the standard Cross and Carreau equations.

Further insights into the environmental effects on the computed hyperfine coupling constants of nitroxides in aqueous solution.

We investigated the main two factors influencing the mean hyperfine coupling constants of small nitroxide radicals in aqueous solution, i.e., the out-of-plane displacement of their nitrogen atom and the environmental effects (solvent effects), by means of the approach we previously developed and fine-tuned to study the solvation of the dimethyl nitroxide radical.

Prediction of nitroxide hyperfine coupling constants in solution from combined nanosecond scale simulations and quantum computations.

We present a combined theoretical approach based on analyzing molecular dynamics trajectories (at the nanosecond scale) generated by use of classical polarizable force fields and on quantum calculations to compute averaged hyperfine coupling constants. That method is used to estimate the constant of a prototypical nitroxide: the dimethylnitroxide. The molecule is embedded during the simulations in a cubic box containing about 500 water molecules and the molecular dynamics is generated using periodic conditions.