Dissoziation schwacher Säuren über den gesamten Molenbruchbereich

Exploring the Influence of the Phosphorus-Heteroatom Substitution in Nicotine on Its Electronic and Vibrational Spectroscopic Properties

The influence of phosphorus substitution of nitrogen in heterocyclic compounds on the vibrational spectroscopy as well as frontier molecular orbitals are analyzed. Nicotine with two nitrogen atoms in its structure is taken as the sample system to be studied computationally. By replacing the nitrogen atom in one or both rings of this molecule with phosphorus, three nicotine derivatives are created. The vibrational circular dichroism and infrared spectra of these four molecules in their monomer state, as well as the assemblies up to trimers are determined. The aforementioned spectra are calculated using static quantum chemical calculations employing a cluster-weighted approach. The calculated gas phase spectra of nicotine are compared to their respective experimental spectra. It is observed that the nicotine derivatives with phosphorus in the methylpyrrolidine ring have considerably different gas phase and bulk phase vibrational circular dichroism spectra when compared to nicotine. The phosphorus substitution reduces the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital as well as altering the polarizability and reactivity of the investigated molecules.

Efficient Prediction of Mole Fraction Related Vibrational Frequency Shifts

While so far it has been possible to calculate vibrational spectra of mixtures at a particular composition, we present here a novel cluster approach for a fast and robust calculation of mole fraction dependent infrared and vibrational circular dichroism spectra at the example of acetonitrile/(R)-butan-2-ol mixtures. By assigning weights to a limited number of quantum chemically calculated clusters, vibrational spectra can be obtained at any desired composition by a weighted average of the single cluster spectra. In this way, peak positions carrying information about intermolecular interactions can be predicted. We show that mole fraction dependent peak shifts can be accurately modeled and, that experimentally recorded infrared spectra can be reproduced with high accuracy over the entire mixing range. Because only a very limited number of clusters is required, the presented approach is a valuable and computationally efficient tool to access mole fraction dependent spectra of mixtures on a routine basis.

Cluster Analysis in Liquids: A Novel Tool in TRAVIS

We present a novel cluster analysis implemented in our open-source software TRAVIS and its application to realistic and complex chemical systems. The underlying algorithm is exclusively based on atom distances. Using a two-dimensional model system, we first introduce different cluster analysis functions and their application to single snapshots and trajectories including periodicity and temporal propagation. Using molecular dynamics simulations of pure water with varying system size, we show that our cluster analysis is size-independent. Furthermore, we observe a similar clustering behavior of pure water in classical and ab initio molecular dynamics simulations, showing that our cluster analysis is universal. In order to emphasize the application to more complex systems and mixtures, we additionally apply the cluster analysis to ab initio molecular dynamics simulations of the [C₂C₁Im][OAc] ionic liquid and its mixture with water. Using that, we show that our cluster analysis is able to analyze the clustering of the individual components in a mixture as well as the clustering of the ionic liquid with water.

Conformer Weighting and Differently Sized Cluster Weighting for Nicotine and Its Phosphorus Derivatives

Weighting methods applied to systems with many conformers have been broadly employed to calculate thermodynamic properties, structural characteristics, and populations. To better understand and test the sensitivity of conventional weighting methods, the conformational distributions of nicotine and its phosphorus-substituted derivatives are investigated. The weighting schemes used for this are all based on Boltzmann statistics. Classical Boltzmann factors based on the electronic energy and the Gibbs free energy are calculated at different quantum chemical levels of theory and compared to cluster weights obtained by the quantum cluster equilibrium method. Furthermore, the influence of the modified rigid-rotor-harmonic-oscillator (mRRHO) approximation on the cluster weights is investigated. The substitution of the nitrogen atom in the methylpyrrolidine ring by a phosphorus atom results in more monomer conformers and clusters being populated. The conformational distribution of the monomers remained stable at different levels of theory and weighting methods. However, going to dimers and trimers, we observe a significant influence of the level of theory, weighting method, and mRRHO cutoff on the populations of these clusters. We show that mRRHO cutoff values of 50 and 100 cm^-1 yield similar results, which is why 50 cm^-1 is recommended as a robust choice. Furthermore, we observe that the global minimum for ΔE₀ and ΔG varies in a few cases and that the global minimum is not always the dominantly occupied structure.

Uncertainty quantification of phase transition quantities from cluster weighting calculations

In this work, we investigate how uncertainties in experimental input data influence the results of quantum cluster equilibrium calculations. In particular, we focus on the calculation of vaporization enthalpies and entropies of seven organic liquids, compare two computational approaches for their calculation and investigate how these properties are affected by changes in the experimental input data. It is observed that the vaporization enthalpies and entropies show a smooth dependence on changes in the reference density and boiling point. The reference density is found to have only a small influence on the vaporization thermodynamics, whereas the boiling point has a large influence on the vaporization enthalpy but only a small influence on the vaporization entropy. Furthermore we employed the Gauss--Hermite estimator in order to quantify the uncertainty in the thermodynamic functions that stems from inaccuracies in the experimental reference data at the example of the vaporization enthalpy of (\textit{R})-butan-2-ol. We quantify the uncertainty as 30.95~$\cdot$10$^{-3}$~kJ~mol$^{-1}$. Additionally we compare the convergence behaviour and computational effort of the Gauss--Hermite estimator with the Monte Carlo approach and show the superiority of the former. By this, we present how uncertainty quantification can be applied to examples from theoretical chemistry.

The Ionic Product of Water in the Eye of the Quantum Cluster Equilibrium

The theoretical description of water properties continues to be a challenge. Using quantum cluster equilibrium (QCE) theory, we combine state-of-the-art quantum chemistry and statistical thermodynamic methods with the almost historical Clausius-Clapeyron relation to study water self-dissociation and the thermodynamics of vaporization. We pay particular attention to the treatment of internal rotations and their impact on the investigated properties by employing the modified rigid-rotor-harmonic-oscillator (mRRHO) approach. We also study a novel QCE parameter-optimization procedure. Both the ionic product and the vaporization enthalpy yield an astonishing agreement with experimental reference data. A significant influence of the mRRHO approach is observed for cluster populations and, consequently, for the ionic product. Thermodynamic properties are less affected by the treatment of these low-frequency modes.

Hydrogen Bonding and Vaporization Thermodynamics in Hexafluoroisopropanol-Acetone and -Methanol Mixtures. A Joined Cluster Analysis and Molecular Dynamic Study

Binary mixtures of hexafluoroisopropanol with either methanol or acetone are analyzed via classical molecular dynamics simulations and quantum cluster equilibrium calculations. In particular, their populations and thermodynamic properties are investigated with the binary quantum cluster equilibrium method, using our in‐house code peacemaker 2.8, upgraded with temperature‐dependent parameters. A novel approach, where the final density from classical molecular dynamics, has been used to generate the necessary reference isobars. The hydrogen bond network in both type of mixtures at molar fraction of hexafluoroisopropanol of 0.2, 0.5, and 0.8 respectively is investigated via the molecular dynamics trajectories and the cluster results. In particular, the populations show that mixed clusters are preferred in both systems even at 0.2 molar fractions of hexafluoroisopropanol. Enthalpies and entropies of vaporization are calculated for the neat and mixed systems and found to be in good agreement with experimental values.

Calculation of improved enthalpy and entropy of vaporization by a modified partition function in quantum cluster equilibrium theory

In this work, we present an altered partition function that leads to an improved calculation of the enthalpy and entropy of vaporization in the framework of quantum cluster equilibrium theory. The changes are based on a previously suggested modification [S. Grimme, Chem. Eur. J. 18, 9955-9964 (2012)] of the molecular entropy calculation in the gas phase. Here, the low energy vibrational frequencies in the vibrational partition function are treated as hindered rotations instead of vibrations. The new scheme is tested on a set of nine organic solvents for the calculation of the enthalpy and entropy of vaporization. The enthalpies and entropies of vaporization show improvements from 6.5 error to 3.3 kJ mol^-1 deviation to experiment and from 28.4 error to 13.5 J mol^-1 K^-1 deviation to experiment, respectively. The effect of the corrected partition function is visible in the different populations of clusters, which become physically more meaningful in that larger clusters are higher populated in the liquid phase and the gas phase is mainly populated by the monomers. Furthermore, the corrected partition function also overcomes technical difficulties and leads to an increased stability of the calculations in regard to the size of the cluster set.

Water in Protic Ionic Liquid Electrolytes: From Solvent Separated Ion Pairs to Water Clusters

The large electrochemical and cycling stability of "water-in-salt" systems have rendered promising prospective electrolytes for batteries. The impact of addition of water on the properties of ionic liquids has already been addressed in several publications. In this contribution, we focus on the changes in the state of water. Therefore, we investigated the protic ionic liquid N-butyl-pyrrolidinium bis(trifluoromethanesulfonyl)imide with varying water content at different temperatures with the aid of molecular dynamics simulations. It is revealed that at very low concentrations, the water is well dispersed and best characterized as shared solvent molecules. At higher concentrations, the water forms larger aggregates and is increasingly approaching a bulk-like state. While the librational and rotational dynamics of the water molecules become faster with increasing concentration, the translational dynamics are found to become slower. Further, all dynamics are found to be faster if the temperature increases. The trends of these findings are well in line with the experimental measured conductivities.

Activity coefficients of binary methanol alcohol mixtures from cluster weighting

The hydrogen bond network of different small alcohols is investigated via cluster analysis. Methanol/alcohol mixtures are studied with increasing chain length and branching of the molecule. Those changes can play an important role in different fields, including solvent and metal extraction. The extended tight binding method GFN2-xTB allows the evaluation and geometry optimization of thousands of clusters built via a genetic algorithm. Interaction energies and geometries are evaluated and discussed for the neat systems. Thermodynamic properties, such as vaporization enthalpies and activity coefficients, are calculated with the binary quantum cluster equilibrium (bQCE) approach using our in-house code peacemaker 2.8. Combined distribution functions of the distances against the angles of the hydrogen bonds are evaluated for neat and mixed clusters and weighted by the equilibrium populations achieved from bQCE calculations.