Solvent Effects and Aggregation Phenomena Studied by Vibrational Optical Activity and Molecular Dynamics: The Case of Pantolactone

Stereochemistry of natural products from vibrational circular dichroism

Secondary metabolites from land and marine (micro)organisms have been at the focus of the drug discovery process for many years. One of the reasons for this success is nature's incredible ability to create intricate molecular scaffolds. Such structural richness, however, makes the structural elucidation, and the absolute configuration assignment in particular, a challenging process. Vibrational circular dichroism (VCD) has emerged as one of the most reliable and versatile methods to unambiguously assign both the absolute configuration and conformations of chiral molecules in solution. Although VCD is no longer a curiosity in the field of molecular spectroscopy after 50 years since its first report, it is still underutilized by natural product chemists worldwide for varying reasons. Herein, we highlight the evolution of the application of VCD to natural product chemistry, focusing on its strengths as well as points that still need improvement. General guidelines for the correct application of VCD to stereochemical studies are also provided.

Testing the Simplified Molecular Dynamics Approach to Improve the Reproduction of ECD Spectra and Monitor Aggregation

A simplified molecular-dynamics-based electronic circular dichroism (ECD) approach was tested on three condensed derivatives with limited conformational flexibility and an isochroman-2H-chromene hybrid, the ECD spectra of which could not be precisely reproduced by the conventional ECD calculation protocol. Application of explicit solvent molecules at the molecular mechanics (MD) level in the dynamics simulations and subsequent TDDFT-ECD calculation for the unoptimized MD structures was able to improve the agreements between experimental and computed spectra. Since enhancements were achieved even for molecules with limited conformational flexibility, deformations caused by the solvent molecules and multitudes of conformers produced with unoptimized geometries seem to be key factors for better agreement. The MD approach could confirm that aggregation of the phenanthrene natural product luzulin A had a significant contribution to a specific wavelength range of the experimental ECD. The MD approach has proved that dimer formation occurred in solution and this was responsible for the anomalous ECD spectrum. The scope and limitations of the method have also been discussed.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Conformational distributions of tetrahydro-2-turoic acid in water at different pH values by their IR and vibrational circular dichroism spectra

Infrared (IR) and vibrational circular dichroism (VCD) spectra of tetrahydro-2-furoic acid (THFA) in aqueous solutions under several different pH conditions were recorded. To interpret the IR and VCD spectra of THFA obtained in highly acidic and basic aqueous solutions, extensive and systematic conformational searches were conducted to acquire the low-energy minima for both the neutral and deprotonated forms of THFA species, as well as their hydrated clusters. This was accomplished by using the conformer-rotamer ensemble sampling tool (CREST) with an implicit solvation model for water. The CREST candidates were further optimized at the B3LYP-D3BJ/def2-TZVP level of theory. The simulated VCD spectra of the neutral THFA conformers in the polarizable continuum model (PCM) of water alone exhibit little agreement with the experimental data under highly acidic conditions. Applying the clusters-in-a-liquid solvation model by considering the monohydrate THFA conformers in the PCM of water, significantly improved agreement with the experimental data. Similarly, the deprotonated THFA species solvated with one to four explicit water molecules in the PCM of water were considered. While the IR and VCD spectra of the deprotonated THFA monohydrate conformers offer the best agreement with the experimental data, other larger hydrated clusters, particularly the dihydrates, also contribute to the spectra. Through the synergistic combined experimental and theoretical approach used in the study, comprehensive conformational distributions of the predominant THFA species across various pH conditions were extracted.

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.

Pitfalls in the Optimization of Conformer Populations to Maximize the Similarity between Predicted and Experimental Chiroptical Spectra

The conformational populations of pantolactone, epichlorohydrin, and N-acetyl-tryptophan methyl ester were investigated by using similarity analysis between their calculated and experimental chiroptical spectra. By performing the analysis on pantolactone using two different chiroptical methods, namely, vibrational circular dichroism and Raman optical activity, it was found that the optimal sets of conformers do not match between the two methods, indicating that the conformational populations obtained by optimizing the similarity between calculated and experimental spectra are unlikely to be more accurate than energy-based Boltzmann populations. Also, it was found for pantolactone, epichlorohydrin, and N-acetyl-tryptophan methyl ester that the similarity between calculated and experimental spectra would often not vary significantly if each of the populated conformers was discarded, one at a time. This observation indicates that more than one set of conformers can provide acceptable similarity between the predicted and experimental spectra. Therefore, the correct set of conformers cannot be accurately determined by similarity analysis.

Influence of microsolvation on vibrational circular dichroism spectra in dimethyl sulfoxide solvent: A Bottom-Up approach using Quantum cluster growth

Chiroptical spectroscopic measurements serve as routine methods to assign the absolute configuration of chiral compounds and interpret their conformational behavior in solution. One common challenge is the use of strongly hydrogen-bonding solvents, which can significantly bias the conformational ensemble and affect the vibrational circular dichroism (VCD) active bands in solution. One such solvent is dimethyl sulfoxide (DMSO)-an excellent solvent for stubborn compounds-that must be explicitly considered in VCD analysis. Explicit consideration of solvent remains a critical challenge in chiroptical spectroscopy due to the need to explore solute-solvent conformational space and the computational expense in modeling these clusters. Interested in the recent development of the Quantum Cluster Growth (QCG) program by the Grimme lab, we set out to model and interpret previously reported VCD spectra for several molecules using their efficient program. Our purposes are two-fold: (1) to investigate the applicability of the QCG program to the problem of reproducing VCD spectra in DMSO solvent and (2) to identify limitations in using this approach. We find that we can conveniently model and analyze the VCD spectra of investigated molecules in DMSO. However, the final set of conformers used for VCD calculations are functional dependent and different sets of conformers can provide satisfactory quantitative agreement between experimental and predicted VCD spectra. We hope that this study provides guidance for future chiroptical studies in the challenging DMSO solvent.

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.

Infrared and vibrational circular dichroism spectra of methyl β-D-glucopyranose in water: The application of the quantum cluster growth and clusters-in-a-liquid solvation models

The infrared (IR) and vibrational circular dichroism (VCD) spectra of methyl β-D-glucopyranose in water were measured. Both implicit and explicit solvation models were utilized to explain the observed spectra. The vast body of existing experimental and theoretical data suggested that about eight explicit water molecules are needed to account for the solvent effects, supported by the current Quantum Cluster Growth (QCG) analysis. Extensive manual and systematic conformational searches of the molecular target and its water clusters were carried out by using a recently developed conformational searching tool, conformer-rotamer ensemble sampling tool (CREST), and the microsolvation model in the associated QCG code. The Boltzmann averaged IR and VCD spectra of the methyl β-D-glucopyranose-(water)n (n = 8) conformers in the PCM of water provide better agreement with the experimental ones than those with n = 0, 1, and 2. The explicit solvation with eight water molecules was shown to greatly modify the conformational preference of methyl β-D-glucopyranose from its monomeric form. Further analyses show that the result is consistent with the existence of long-lived methyl β-D-glucopyranose monohydrates with the additional explicit water effects being accounted for with the quantum mechanical treatment of the other seven close-by water molecules in the PCM of water.

Comparison of different density functional theory methods for the calculation of vibrational circular dichroism spectra

The determination of the absolute configuration (AC) of an organic molecule is still a challenging task for which the combination of spectroscopic with quantum-mechanical methods has become a promising approach. In this study, we investigated the accuracy of DFT methods (480 overall combinations of 15 functionals, 16 basis sets, and 2 solvation models) to calculate the VCD spectra of six chiral organic molecules in order to benchmark their capability to facilitate the determination of the AC.

Conformational analysis of amphetamine and methamphetamine: a comprehensive approach by vibrational and chiroptical spectroscopy

After cannabis, the most commonly used illicit substance worldwide is amphetamine and its derivatives, such as methamphetamine, with an ever-increasing number of synthetic modifications. Thus, fast and reliable methods are needed to identify them according to their spectral patterns and structures. Here, we have investigated the use of molecular spectroscopy methods to describe the 3D structures of these substances in a solution that models the physiological environment. The substances were analyzed by Raman and infrared (IR) absorption spectroscopy and by chiroptical methods, vibrational circular dichroism (VCD) and Raman optical activity (ROA). The obtained experimental data were supported by three different computational approaches based on density functional theory (DFT) and molecular dynamics (MD). Successful interpretation relies on good agreement between experimental and predicted spectra. The determination of the conformer populations of the studied molecules was based on maximizing the similarity overlap of weighted conformer spectra by a global minimization algorithm. Very good agreement was obtained between the experimental spectra and optimized-population weighted spectra from MD, providing a detailed insight into the structure of the molecules and their interaction with the solvent. The relative population of three amphetamine and six methamphetamine conformers was determined and is consistent with a previous NMR study. However, this work shows that only a few isolated conformers are not sufficient for the successful interpretation of the spectra, but the entire conformational space needs to be sampled appropriately and explicit interaction with the solvent needs to be included.

Chirality transfer observed in Raman optical activity spectra

Chirality transfer (also called induced chirality) is a phenomenon present in chiroptical spectra that manifests itself as a new band or bands of an achiral molecule interacting with a chiral one. In the Raman optical activity (ROA) spectroscopy, the bands of achiral solvents have been recently observed, but the latest papers have shown that they corresponded to the new ECD-Raman (eCP-Raman) effect. Here, we show an unambiguous example of chirality transfer observed in the ROA spectra. The spectra registered for the (1:1) mixtures of achiral benzonitrile with the enantiomers of 2,2,2-trifluoro-1-phenylethanol, 1-phenylethanol, and 1-phenylethylamine exhibited the v(CN) vibration band at about 2230 cm^-1. The ROA measurements were repeated several times to ensure the reliability of the phenomenon. Calculations revealed the CN···HO or CN···HNH hydrogen bond formation accompanied by the π···π or CH···π interactions. The interaction strength was shown to be an important factor for the pronouncement of the ROA chirality transfer effect.

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.

Dynamic Stereochemistry of a Biphenyl-Bisprolineamide Model Catalyst and its Imidazolidinone Intermediates

In this study, we characterize the dynamic stereochemistry of a biphenyl-2,2'-bis(proline amide) catalyst in chloroform and DMSO as representative weakly and strongly hydrogen bonding solvents. Using vibrational circular dichroism (VCD) spectroscopy and density functional theory (DFT) based spectra calculations, we show that the preferred axial stereochemistry of the catalyst is determined by solute-solvent interactions. Explicitly considering solvation with DMSO molecules is found to be essential to correctly predict the conformational preferences of the catalyst. Furthermore, we investigate the stereochemistry of the corresponding enamines and imidazolidinones that are formed upon reaction with isovaleraldehyde. The enamines are found to rapidly convert to endo-imidazolidinones and the thermodynamically favored exo-imidazolidinones are formed only slowly. The present study demonstrates that the stereochemistry of these imidazolidinones can be deduced directly from the VCD spectra analysis without any further detailed analysis of NMR spectra. Hence, we herein exemplify the use of VCD spectroscopy for an in situ characterization of intermediates relevant in asymmetric catalysts.

Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database

A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.

VCD spectroscopy reveals conformational changes of chiral crown ethers upon complexation of potassium and ammonium cations

Two chiral derivatives of 18-crown-6, namely the host molecules 2,3-diphenyl- and 2-phenyl-18c6, serve as model systems to investigate whether VCD spectroscopy can be used to monitor conformational changes occurring upon complexation of guests. Host-guest complexes of both crown ethers were prepared by addition of KNO₃. The more bulky 2,3-diphenyl-18c6 is found to undergo major conformational changes upon encapsulation of K⁺, which are revealed as characteristic changes of the VCD spectral signatures. In contrast, while 2-phenyl-18c6 also incorporates K⁺ into the macrocycle, strong conformational changes are not occurring and thus spectral changes are negligible. With an octyl ammonium cation as guest molecule, 2,3-diphenyl-18c6 shows the same conformational and spectral changes that were observed for K⁺-complexes. In addition, the asymmetric NH₃-deformation modes are found to gain VCD intensity through an induced VCD process. An analysis of the vibrational spectra enables a differentiation of VCD active and inactive guest modes: There appears to be a correlation between the symmetry of the vibrational mode and the induced VCD intensity. While this finding makes the host-guest complexes interesting systems for future theoretical studies on the origin of induced VCD signatures, the observations described in this study demonstrate that VCD spectroscopy is indeed a suitable technique for the characterization of supramolecular host-guest complexes.

Raman Optical Activity of N-Acetyl-L-Cysteine in water and in methanol: the "clusters-in-a-liquid" model and ab initio molecular dynamics simulations

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, i.e. the DFT based clusters-in-a-liquid solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations , which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.

The solid state VCD of a novel N-acylhydrazone trifluoroacetate

A novel N-acylhydrazone with pharmaceutical importance was subject of structural and IR/VCD investigations in the solid state. In the crystal structure, dimers of anion-cation pairs are stabilized by H-bonding and ionic interactions. Some less common interaction types, like C=N···C-NH₃⁺ (σ-hole) interactions, hydrazone-aromatic interactions and dispersive contacts of the CF₃ groups are also present in the crystal. Satisfactory reproduction of the solid state IR and VCD spectra required that quantum-chemical calculations be done on a tetramer (four cation-anion pairs) cut out from the crystal structure, exhibiting key intermolecular interactions. Ten DFT functionals were assessed as to the agreement between the calculated and experimental spectra. Various approaches to scaling of the calculated frequencies were applied. The best results were yielded with individual (optimized) frequency scaling factors (FSFs) and band half-widths at half maximum-(HWHM) for four separate spectral subregions. The best matching between the experimental and theoretical spectra (according to SimIR, SimVCD and SimVDF indices) was found for the B3PW91 functional, however, a few other functionals follow closely in the ranking. Based on the quantum chemical calculations, spectral assignments have been made.

Helicity control of a perfluorinated carbon chain within a chiral supramolecular cage monitored by VCD

Confinement within supramolecular systems is the leading technology to finely tune guest functional properties. In this communication we report the synthesis of a chiral supramolecular cage able to bias the helicity of a perfluorinated carbon chain hosted within the cage. We monitor the phenomenon of chiral induction by Vibrational Circular Dichroism (VCD) experiments complemented by DFT calculations over the possible conformers.

Use of Raman and Raman optical activity to extract atomistic details of saccharides in aqueous solution

Sugars are crucial components in biosystems and industrial applications. In aqueous environments, the natural state of short saccharides or charged glycosaminoglycans is floating and wiggling in solution. Therefore, tools to characterize their structure in a native aqueous environment are crucial but not always available. Here, we show that a combination of Raman/ROA and, on occasions, NMR experiments with Molecular Dynamics (MD) and Quantum Mechanics (QM) is a viable method to gain insights into structural features of sugars in solutions. Combining these methods provides information about accessible ring puckering conformers and their proportions. It also provides information about the conformation of the linkage between the sugar monomers, i.e., glycosidic bonds, allowing for identifying significantly accessible conformers and their relative abundance. For mixtures of sugar moieties, this method enables the deconvolution of the Raman/ROA spectra to find the actual amounts of its molecular constituents, serving as an effective analytical technique. For example, it allows calculating anomeric ratios for reducing sugars and analyzing more complex sugar mixtures to elucidate their real content. Altogether, we show that combining Raman/ROA spectroscopies with simulations is a versatile method applicable to saccharides. It allows for accessing many features with precision comparable to other methods routinely used for this task, making it a viable alternative. Furthermore, we prove that the proposed technique can scale up by studying the complicated raffinose trisaccharide, and therefore, we expect its wide adoption to characterize sugar structural features in solution.

VCD spectroscopy distinguishes the enamine and iminium ion of a 1,1’-binaphthyl azepine

We present a VCD spectroscopic characterization of a chiral 1,1’-binaphthyl azepine catalyst and show that the VCD spectra of an in-situ generated enamine and an ex-situ prepared iminium ion are...

Chiral analysis of pantolactone with molecular rotational resonance spectroscopy

A molecular rotational resonance spectroscopy method for measuring the enantiomeric excess of pantolactone, an intermediate in the synthesis of panthenol and pantothenic acid, is presented. The enantiomers are distinguished via complexation with a small chiral tag molecule, which produces diastereomeric complexes in the pulsed jet expansion used to inject the sample into the spectrometer. These complexes have distinct moments of inertia, so their spectra are resolved by MRR spectroscopy. Quantitative enantiomeric excess (EE) measurements are made by taking the ratio of normalized complex signal levels when a chiral tag sample of high, known EE is used, while the absolute configuration of the sample can be determined from electronic structure calculations of the complex geometries. These measurements can be performed without the need for reference samples with known enantiopurity. Two instruments were used in the analysis. A broadband, chirped-pulse spectrometer is used to perform structural characterization of the complexes. The broadband spectrometer is also used to determine the EE; however, this approach requires relatively long measurement times. A targeted MRR spectrometer is also used to demonstrate EE analysis with approximately 15-min sample-to-sample cycle time. The quantitative accuracy of the method is demonstrated by comparison with chiral gas chromatography and through the measurement of a series of reference samples prepared from mixtures of (R)-pantolactone and (S)-pantolactone samples of known EE.

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.

How important are the intermolecular hydrogen bonding interactions in methanol solvent for interpreting the chiroptical properties?

Two crispine A analogs and tetrahydrofuro[2,3-b]furan-3,3a(6aH)-diol, endowed with hydroxyl groups that can participate in intramolecular hydrogen bonding, have been synthesized and experimental vibrational circular dichroism (VCD) spectra and optical rotatory dispersion (ORD) data have been measured in CD₃OD/CH₃OH solvents. The absolute configurations (ACs) of these compounds have been determined using their synthetic schemes, supplemented wherever possible with X-ray diffraction data. The ACs are also analyzed with quantum chemical (QC) calculations of VCD and ORD utilizing implicit solvation as well as explicit solvation models, with the later employing classical molecular dynamics (MD) simulations. It is found that VCD calculations with implicit solvation model are adequate for determining the ACs, despite propensity of studied compounds for intermolecular hydrogen bonding between solute and solvent molecules. This observation is important because time-consuming MD simulations may not be necessary in the type of situations studied here. Additionally, it is found that the QC predicted VCD spectra provided enough diastereomer discrimination for determining the correct AC of studied compounds independently. The same observation did not apply to ORD.

Absolute Configuration and Conformation of (-)-R-t-Butylphenylphosphinoamidate: Chiroptical Spectroscopy and X-ray Analysis

The absolute configuration and conformations of (−)-tert-butylphenylphosphinoamidate were determined using three different chiroptical spectroscopic methods, namely vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and optical rotatory dispersion (ORD). In each of the spectroscopic methods used, experimental data for the (−)-enantiomer of tert-butylphenylphosphinoamidate were measured in the solution phase. Using the concentration-dependent experimental infrared spectra, the existence of dimers in the solution was investigated, and the monomer–dimer equilibrium constant was determined. Concomitant quantum mechanical predictions of the VCD, ECD, and ORD for monomeric tert-butylphenylphosphinoamidate were carried out using density functional theory (DFT) calculations using the B3LYP functional and the 6-31G(d), 6-311G(2d,2p) and aug-cc-pVDZ basis sets. Similar predictions for dimeric tert-butylphenylphosphinoamidate were also obtained using the B3LYP/6-31G(d) method. A comparison of theoretically predicted data with the corresponding experimental data led to the elucidation of the absolute configuration as (−)-(R)-tert-butylphenylphosphinoamidate with one predominant conformation in the solution. This conclusion was independently supported by X-ray analysis of the complex with (+)-R-2,2′-dihydroxy-1,1′-binaphthol ((+)-R- BINOL).