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.

A General Iron-Catalyzed Decarboxylative Oxygenation of Aliphatic Carboxylic Acids.

We report an iron-catalyzed decarboxylative C(sp3)-O bond-forming reaction under mild, base-free conditions with visible light irradiation. The transformation uses readily available and structurally diverse carboxylic acids, iron photocatalyst, and 2,2,6,6-tetramethylpiperidine 1-oxyl (TEMPO) derivatives as oxygenation reagents. The process exhibits a broad scope in acids possessing a wide range of stereoelectronic properties and functional groups. The developed reaction was applied to late-stage oxygenation of a series of bio-active molecules. The reaction leverages the ability of iron complexes to generate carbon-centered radicals directly from carboxylic acids by photoinduced carboxylate-to-iron charge transfer. Kinetic, electrochemical, EPR, UV/Vis, HRMS, and DFT studies revealed that TEMPO has a triple role in the reaction: as an oxygenation reagent, an oxidant to turn over the Fe-catalyst, and an internal base for the carboxylic acid deprotonation. The obtained TEMPO adducts represent versatile synthetic intermediates that were further engaged in C-C and C-heteroatom bond-forming reactions using commercial organo-photocatalysts and nucleophilic reagents.

Selective Chirality Transfer to the Bis(trifluoromethylsulfonyl)imide Anion of an Ionic Liquid.

Chirality transfer from the chiral molecule (R)-1,2-propylene oxide to the achiral anion of the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquid is observed. The chiral probe selectively affects one part of the binary ionic liquid, i. e., it has previously been shown experimentally and theoretically that this particular imidazolium cation can be affected by chirality transfer, but in the present system chirality is almost exclusively transferred to the anion and not to both parts of the solvent (anion and cation). This observation is of high relevance because of its selectivity and because anion effects are usually much more important in ionic liquid research than cation effects. From ab initio molecular dynamics simulations, a conformational analysis and dissected vibrational circular dichroism spectra are obtained to study the chirality transfer. While in the neat ionic liquid two mirror imaged trans conformers of the anion occur almost equally, we observe an excess of one of these conformers in the presence of the chiral solute, causing optical activity of the anion. Although the cis conformers are not tremendously affected by the chirality transfer, they gain in total population when (R)-1,2-propylene oxide is dissolved in the ionic liquid.

Solvent Dependency of Catalyst-Substrate Aggregation Through π-π Stacking in Photoredox Catalysis.

Assemblies of photoredox catalysts and their target substrates prior to photoexcitation is a phenomenon naïvely overlooked by the majority of synthetic chemists, but can have profound influences on reactivity and selectivity in photocatalytic reactions. In this study, we determine the aggregation states of triarylamine radical cationic photocatalysts with various target arene substrates in different solvents by specifically parameterized polarizable molecular dynamics simulations. A π-stacking interaction previously implicated by more expensive, less-representative quantum calculations is confirmed. Critically, this study presents new insights on: i) the ability of solvents (MeCN vs DMF) to make or break a photocatalytic reaction by promoting (MeCN) or demoting (DMF) its catalyst-substrate assemblies, which is a determining factor for reactivity, ii) the average "lifetimes" of assemblies in solution from a dynamic simulation. We find that both in the ground state and the photoexcited state, the cationic radical assemblies remain intact for periods often higher than 60 ps, rendering them ideally suitable to undergo intra-assembly electron transfer reactions upon photoexcitation. Such aspects have not addressed by previous studies on synthetic photocatalytic reactions involving non-covalent assemblies.

Stronger Together! Mechanistic Investigation into Synergistic Effects during Homogeneous Carbon Dioxide Hydrogenation Using a Heterobimetallic Catalyst.

A series of group 6 heterobimetallic complexes [M0;IrIII] (M = Cr, Mo, W) were synthesized and fully characterized, and the catalytic behavior was studied. The heterobimetallic complex [Mo0;IrIII] (C1) was by far the most active and has shown a considerable synergistic effect, with both metals actively participating in homogeneous carbon dioxide hydrogenation, leading to formate salts. Based on theoretical calculations, the synergistic interaction is due to Pauli repulsion, lowering the transition state and thus enabling higher catalytic activity. The mechanism of both the hydrogenation itself and the synergistic interaction was studied by NMR spectroscopy, kinetic measurements, and theoretical calculations. The homogeneous nature of the reaction was proven using in situ high-pressure (HP) NMR experiments. The same experiments also showed that the octahedral Mo(CO)3P3 moiety of the complex is stable under the reaction conditions. The hydride complex is the resting state because the hydride transfer is the rate-determining step. This is supported by kinetic measurements, in situ HP NMR experiments, and theoretical calculations and is in contrast to the monometallic IrIII counterpart of C1.

CONAN─Novel Tool to Create and Analyze Liquids in Confined Space.

Modeling of complex liquids at solid surfaces and in confinement is gaining attention due to an increase in computer power and advancement of simulation techniques. Therefore, tools to set up structures and for analysis are needed. In this paper, we present CONAN─a Python code designed to facilitate the study of liquids interacting with solid structures, such as walls or pores. Among other things, the program provides the option to generate a variety of different structures, including carbon walls and nanotubes and their boron nitride analogs, as well as the ability to analyze various structural properties of confined and interfacial liquids. In the case of the ionic liquid 1-butyl-3-methylimidazolium acetate in carbon nanotubes of different sizes, we demonstrate the abilities of our tool. The average density within the confinement highly depends on the carbon nanotube size, and it is generally lower than the density of the bulk liquid. The arrangement of the individual species within the tube also depends on size, with radial layers forming within the tubular confinement. The density is largely increased in the respective layers, while it is drastically reduced between the layers.

Locality in amino-acid based imidazolium ionic liquids.

Several amino-acid based imidazolium ILs are investigated through the use of ab initio molecular dynamics (AIMD), which includes full polarization. The electric dipole moment distribution and polarization is used as a means of characterizing and understanding these complex systems. Various charge scheme methods were analyzed (Wannier function, Blöchl, Löwdin and Mulliken charge schemes and Voronoi tessellation) to determine their ability to predict dipole moments. These results and the following comparison of methods further deepen the knowledge of polarization by highlighting the importance of the anion and cation separately on polarizability contribution and the need to select a suitable method to predict these. The angular probability distribution is utilized to measure the degree of locality in monopole-dipole electrostatic interactions, which showed no preferential alignment over 700 pm. In addition, the IR and Raman spectra from Voronoi tessellation of [C2C1Im][ala] were analyzed. In these, the strongest signalling peaks showed consistency with experiment and the ability to differentiate between anion and cation components of the IL.

Ring-Opening Reaction of 1-Phospha-2-Azanorbornenes via P-N Bond Cleavage and Reversibility Studies.

The reactive P-N bond in 1-phospha-2-azanorbornenes is readily cleaved by simple alcohols to afford P-chiral 2,3-dihydrophosphole derivatives as a racemic mixture. The isolation of the products was not possible due to the reversibility of the reaction, which could, however, be stopped by sulfurization of the phosphorus atom, thus efficiently blocking the lone pair of electrons, as exemplified for 6b yielding structurally characterized 8b. Additionally, the influence of the substituent in the α position to the phosphorus atom (H, Ph, 2-py, CN) on the reversibility of the reaction was studied. Extensive theoretical calculations for understanding the ring-closing mechanism suggested that a multi-step reaction with one or more intermediates was most probable.

Synergistic Catalysis in Heterobimetallic Complexes for Homogeneous Carbon Dioxide Hydrogenation.

Two heterobimetallic Mo,M' complexes (M' = IrIII, RhIII) were synthesized and fully characterized. Their catalytic activity in homogeneous carbon dioxide hydrogenation to formate was studied. A pronounced synergistic effect between the two metals was found, most notably between Mo and Ir, leading to a fourfold increase in activity compared with a binary mixture of the two monometallic counterparts. This synergism can be attributed to spatial proximity of the two metals rather than electronic interactions. To further understand the nature of this interaction, the mechanism of the CO2 hydrogenation to formate by a monometallic IrIII catalyst was studied using computational and spectroscopic methods. The resting state of the reaction was found to be the metal-base adduct, whereas the rate-determining step is the inner-sphere hydride transfer to CO2. Based on these findings, the synergism in the heterobimetallic complex is beneficial in this key step, most likely by further activating the CO2.

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.

Recognition in Chiral Ionic Liquids: The Achiral Cation Makes the Difference!

By simulating butan-2-ol dissolved in the chiral ionic liquid 1-ethyl-3-methylimidazolium (S)-alaninate, we investigate the chiral recognition of butan-2-ol in the ionic liquid. The hydrogen bonding between the chiral anion and both enantiomers of butan-2-ol is similar; however, both chiral molecules (anion and alcohol) induce an asymmetry in the achiral cation which leads to a more favorable environment for the alcohol in the heterochiral system as compared to the homochiral system and hence provides an energetic stabilization of the former.

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 [C2C1Im][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.

Charge transfer and polarisability in ionic liquids: a case study.

The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C4C1im][NTf2] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.

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 ΔE0 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 thermodynamic functions that stems from inaccuracies in the experimental reference data for the example of the vaporization enthalpy of (R)-butan-2-ol. We quantify the uncertainty as 30.95 · 10-3 kJ mol-1. In addition, we compare the convergence behavior and computational effort of the Gauss-Hermite estimator with the Monte Carlo approach and show the superiority of the former. Using this study, 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.

Molecular level insight into the solvation of cellulose in deep eutectic solvents.

Deep eutectic solvents as sustainable and new-generation solvents show potential in the field of cellulose dissolution. Although these novel materials are tested for numerous industrial, environmental, and medical applications, little is known about the structural features of cellulose interacting with deep eutectic solvents. In this work, the interplay of cellulose is studied in two deep eutectic solvents: choline acetate mixed with urea and choline chloride mixed with urea using classical molecular dynamics simulations. Dissolution of cellulose in the studied liquids was not observed to be in agreement with experimental work from the literature. However, a slight swelling in the chloride, as compared to the acetate-based solvent, is apparent. A possible rationale might be found in the stronger hydrogen bonding of the chloride anion compared to the acetate anion with the hydrogen atoms of the cellulose. Moreover, chloride approaches the outer glucose units comparatively more, which could be interpreted as the onset of entering and thus dissolving the cellulose as was previously observed. Specific hydrogen bonds between all units are analyzed and discussed in detail.

How is CO2 absorbed into a deep eutectic solvent?

Deep eutectic solvents show great potential as CO2 absorbents, which is highly desirable for the sustainable development of CO2 reduction and prevention of global climate changes. Ab initio molecular dynamics simulations in the isothermal-isobaric ensemble at pressures of 1 MPa and 5 MPa and at the corresponding experimental density are carried out to investigate the CO2 absorption in choline chloride: ethylene glycol deep eutectic solvent. Based on the structural analysis, there is a strong anion and hydrogen bond donor effect and a minor cation effect on CO2 solvation in the solvent. Instead of cooperation, a competition between the anion and the hydrogen bond donor (ethylene glycol) for the interaction with CO2 is indicated. While at a lower pressure, the ethylene glycol-CO2 interaction dominates, at a higher pressure, it is the chloride-CO2 interaction. Thus, it is possible to use the same advantages within the deep eutectic solvent as the CO2 absorbent as in ionic liquids, but in the hydrogen bond, a donor can be exploited.

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.

Understanding the Complex Surface Interplay for Extraction: A Molecular Dynamics Study.

By means of classical molecular dynamics simulation the interfacial properties of methanol and n-dodecane, which are two potential candidate solvents for use in non-aqueous liquid-liquid extraction, were assessed. The question of how the interface changes depending on the concentration of extractant (tri-n-butyl phosphate) and salt (LiCl) is addressed. Two different models to represent systems were used to evaluate how LiCl and tri-n-butyl phosphate affect mutual miscibility, and how the last-named behaves depending on the chemical environment. Tri-n-butyl phosphate increases the mutual solubility of the solvents, whereas LiCl counteracts it. The extractant was found to be mostly adsorbed on the interface between the solvents, and therefore the structural features of the adsorption were investigated. Adsorption of tri-n-butyl phosphate changes depending on its concentration and the presence of LiCl. It exhibits a preferential orientation in which the butyl chains point at the n-dodecane phase and the phosphate group points at the methanol phase. For high concentrations of tri-n-butyl phosphate, its molecular orientation is preserved by diffusion of the excess molecules into both the methanol and n-dodecane phases. However, LiCl hinders the diffusion into the methanol phase, and thus increases the concentration of tri-n-butyl phosphate at the interface and forces a rearrangement with subsequent loss of orientation.