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.

Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach

Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.

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.

Conformational preferences of the flexible galactofuranose sugar in gas-phase

In contrast with the predominant pyranose form of galactose, galactofuranose is known to be highly flexible. Such flexibility poses a remarkable challenge in terms of structural studies, thus hindering the in depth understanding of the structure/function relationship in this rare sugar. A thorough computational study based on molecular dynamics and density functional theory supported by vibrational spectroscopy in the gas phase was carried out to provide a better understanding of the instrinsic conformational preferences of galactofuranose. Based on energetic and spectroscopic criteria, we report a subtantially reduced conformational landscape: methyl α-D-galactofuranose adopts E₂/¹E conformations and methyl β-D-galactofuranose adopts ¹T₂/¹E conformations.

Vibrational circular dichroism spectra of natural products by means of the nuclear velocity perturbation theory

We present the application of the recently implemented nuclear velocity perturbation theory, using the combined Gaussian and plane waves approach in CP2K, to the vibrational circular dichroism (VCD) spectra of a set of natural products. Even though the calculations were carried out for isolated molecules in the gas-phase limit, neglecting inter-molecular interactions and anharmonic effects, the match between simulated and experimental spectra is reasonable. We also study the influence of different density functionals on the conformational search and the resulting VCD spectra via group coupling matrices (GCMs). The GCM analysis reveals that the VCD signal can in some cases arise from moieties which are close to each other and in other cases from moieties far from each other. Differences in spectra obtained using different exchange-correlation density functionals can be attributed to interaction terms between different moieties in the molecules changing their sign.

How Crystal Symmetry Dictates Non-Local Vibrational Circular Dichroism in the Solid State

Solid-State Vibrational Circular Dichroism (VCD) can be used to determine the absolute structure of chiral crystals, but its interpretation remains a challenge in modern spectroscopy. In this work, we investigate the effect of a twofold screw axis on the solid-state VCD spectrum in a combined experimental and theoretical analysis of P2₁ crystals of (S)-(+)-1-indanol. Even though the space group is achiral, a single proper symmetry operation has an important impact on the VCD spectrum, which reflects the supramolecular chirality of the crystal. Distinguishing between contributions originating from molecular chirality and from chiral crystal packing, we find that while IR absorption hardly depends on the symmetry of the space group, the situation is different for VCD, where completely new non-local patterns emerge. Understanding the two underlying mechanisms, namely gauge transport and direct coupling, will help to use VCD to distinguish polymorphic forms.

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.

Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method

The paper illustrates the Activity Weighted Velocities (AWV) methodology to compute Vibrational Circular Dichroism (VCD) anharmonic spectra from Density Functional Theory (DFT) molecular dynamics. AWV calculates the spectra by the Fourier Transform of the time correlation functions of velocities, weighted by specific observables: the Atomic Polar Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows to correctly reproduce the experimental spectra for systems in the gas and liquid phases, both in the case of weakly and strongly interacting systems. The comparison with the experimental spectra is striking especially in the fingerprint region, as demonstrated by the three benchmark systems discussed: (1S)-Fenchone in the gas phase, (S)-(−)-Propylene oxide in the liquid phase, and (R)-(−)-2-butanol in the liquid phase. The time evolution of APTs and AATs can be adequately described by a linear combination of the tensors of a small set of appropriate reference structures, strongly reducing the computational cost without compromising accuracy. Additionally, AWV allows the partition of the spectral signal in its molecular components without any expensive postprocessing and any localization of the charge density or the wave function.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

Spectroscopic Signatures of Hydrogen-Bonding Motifs in Protonic Ionic Liquid Systems: Insights from Diethylammonium Nitrate in the Solid State

Diethylammonium nitrate, [N0 0 2 2][NO3], and its perdeuterated analogue, [N D D 2 2] [NO3], were structurally characterized and studied by infrared, Raman, and inelastic neutron scattering (INS) spectroscopy. Using these experimental data along with state-of-the-art computational materials modeling, we report unambiguous spectroscopic signatures of hydrogen-bonding interactions between the two counterions. An exhaustive assignment of the spectral features observed with each technique has been provided, and a number of distinct modes related to NH···O dynamics have been identified. We put a particular emphasis on a detailed interpretation of the high-resolution, broadband INS experiments. In particular, the INS data highlight the importance of conformational degrees of freedom within the alkyl chains, a ubiquitous feature of ionic liquid (IL) systems. These findings also enable an in-depth physicochemical understanding of protonic IL systems, a first and necessary step to the tailoring of hydrogen-bonding networks in this important class of materials.

The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

We present a computational study of vibrational circular dichroism (VCD) in solutions of (S)-lactic acid, relying on ab initio molecular dynamics (AIMD) and full solvation with bulk water. We discuss the effect of the hydrogen bond network on the aggregation behaviour of the acid: while aggregates of the solute represent conditions encountered in a weakly interacting solvent, the presence of water drastically interferes with the clusters - more strongly than originally anticipated. For both scenarios we computed the VCD spectra by means of nuclear velocity perturbation theory (NVPT). The comparison with experimental data allows us to establish a VCD-structure relationship that includes the solvent network around the chiral solute. We suggest that fundamental modes with strong polarisation such as the carbonyl stretching vibration can borrow VCD from the chirally restructured solvent cage, which extends the common explanatory models of VCD generation in aqueous solution.

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.