Differentiation of solvatomorphs of active pharmaceutical ingredients (API) by solid-state vibrational circular dichroism (VCD)

Here, we present the new application of solid-state Vibrational Circular Dichroism (VCD) spectroscopy to differentiate several dutasteride (DS) solvatomorphs - the model active pharmaceutical ingredient (API). Several crystalline DS hydrochloride hydrates solvated with methanol, ethanol, acetonitrile, acetone, and acetic acid were prepared. In contrast to almost identical IR spectra, the VCD ones were very sensitive to changes in the sample composition. We marked significant differences in the shape of VCD spectra of studied DS solvatomorphs, DS hydrates, and DS polymorphic forms. Our findings, supported by DFT calculations, show that VCD spectroscopy has the pronounced ability to distinguish their crystal arrangements. We believe that this contribution will extend the use of VCD in the pharmaceutical industry for developing and designing new chiral drug products for the identification, description, and in-depth probing of several pharmaceutical solvatomorphs in the future.

VCD spectra of chiral naphthalene-1-carboxamides in the solid-state

The VCD spectra of chiral 2,3-dihydro-1H-benzo[de]isoquinolin-1-one (8-substituted naphthalene-1-carboxamide, BIQ) were studied in KBr pellets. The X-ray diffractometry revealed that the Me, Ph, and pClPh BIQs crystalize in the monoclinic P21, while nBu, pMePh, and oMeOPh BIQs in the orthorhombic P212121 space group. Only the Me-BIQ crystal exhibits the presence of cyclic amide dimers, while the others contain chains of the amid group hydrogen bonds. For all BIQs, except pMePh, the most intense IR band in the 1750-1550 cm-1 region is located at ca. 1680 cm-1 and is accompanied by two weak ones at ca. 1618 and 1590 cm-1. For the pMePh derivative, four almost equally intense IR bands at 1662, 1639, 1614, and 1588 cm-1 are observed. This region of the IR spectra of BIQs, but pMePh, is well reproduced by calculations based on BIQ monomers. On the other hand, the complex IR pattern of pMePh is computationally reproduced when larger crystal fragments, like octamers, are considered. Registration of the VCD spectra enabled recognizing the complexity of IR contours at ca. 1680 cm-1 by the corresponding VCD motives. For (i) Me, Ph and pClPh (R)-enantiomers, two (+)(-) bands were distinguished and for (ii) nBu and pMePh ones, one VCD band with right-side asymmetry was found. For (iii) oMeOPh the VCD pattern cannot be unambiguously assigned. Thus, the VCD spectra in the ν(C=O) range diverse the studied compounds. Among the set of molecules, pMePh has exceptional crystal geometry. Therefore, its most intense ν(C=O) band position and shape can be connected with the geometry of the hydrogen bonds, interactions, and crystal packing. Interpretation of the VCD spectra is based on linear and packed BIQ octamers. This cluster model can reproduce the main features of the solid-state VCD of BIQs.

Rearrangement of thiazolidine derivatives - a synthesis of a chiral fused oxathiane-γ-lactam bicyclic system

Reaction of L-cysteine with carbonyl compounds leads to thiazolidine derivatives which undergo a stereoselective conversion to two types of chiral bicyclic products bearing two or three stereogenic centers, including the first fused oxathiane-γ-lactam system.

UV-vis and ECD spectroelectrochemistry of atropisomeric naphthalenediimide derivative

The UV-vis and ECD spectroelectrochemistry (SEC) of a chiral binaphthalenylamine derivative of the N-butyl naphthalenediimide (NDIB-NH₂) enantiomers were applied to measure UV-vis and ECD spectra of NDIB-NH₂ radicals and dianion formed in the reduction and oxidation processes observed in cyclic voltammetry (CV). The CV curves and EPR spectroelectrochemistry enabled us to establish conditions at which a radical-anion [NDIB-NH₂]̇^.-, a dianion [NDIB-NH₂]^2-, and a radical-cation [NDIB-NH₂]̇^.+ are formed. The DFT restricted open-shell CAM-B3LYP-D3/def2TZVP/PCM calculations demonstrated that in the radical-anion [NDIB-NH₂]̇^.-, spin is spread over the NDI system while in the radical-cation [NDIB-NH₂]̇⁺ it is spread over the aminonaphthalene moiety. The UV-vis spectra of radical-anion and dianion show the most significant changes in the 400-800 nm range. In that range, the ECD spectra varied with the change of electrode potential more than the UV-vis did and enabled the identification of a new ECD band of [NDIB-NH₂]̇^.- at ca. 400 nm hidden in the background in the UV spectra at -1000 mV. A broad structured ECD pattern with a maximum at ca. 530 nm was observed for [NDIB-NH₂]̇^.- (-1000 mV), while a single smooth ECD band of [NDIB-NH₂]^2- was located at 520 nm (-1750 mV). For the first time, an isosbestic point (455 nm) was found in ECD spectroelectrochemical measurements for the radical-cation [NDIB-NH₂]̇^.+ in equilibrium with the NDIB-NH₂ neutral form. The TD-DFT CAM-B3LYP-D3/6-31G** calculations combined with the hybrid (explicit combined with implicit) solvation model fairly well reproduced the UV-vis and ECD SEC of neutral and redox forms of NDIB-NH₂ but the ECD spectrum of [NDIB-NH₂]̇^.+ above 390 nm.

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.

New chiral ECD-Raman spectroscopy of atropisomeric naphthalenediimides

In this study, we found that a recently discovered ECD-Raman effect dominated over the natural Raman optical activity in a series of atropisomeric naphthalenediimides, and we investigated the kind of information about the molecular structure that could be obtained from the spectra. The ECD-Raman effect is polarised Raman scattering modulated by electronic circular dichroism. We showed that the spectra significantly depended on the substitution of the solute and/or the change of the solvent. Moreover, the spectra could be well-predicted by the theory, thus providing an interesting tool to monitor the chirality of the binaphthyl compounds.

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.

On Raman optical activity sign-switching between the ground and excited states leading to an unusual resonance ROA induced chirality

Raman optical activity (ROA) spectra recorded for a chiral naphthalene diimide derivative (nBu-NDI–BINAM) dissolved in a series of solvents exhibit strong solute to solvent induced chirality with: (1) dominating bands of solvents, (2) nBu-NDI–BINAM resonance ROA (RROA) bands which are barely visible, (3) monosignate RROA Solvent spectra with an unexpected sign concordant with that of the ECD band of the resonant electronic state, (4) bisignate RROA bands for a few solvents, and (5) superposition of non-resonant and resonant ROA bands of the chiral solvents. The unusual ROA enhancement was explained in terms of resonance energy transfer with resonant Raman emission. The surprising RROA sign-switching was found to be due to specific conformational equilibria where one solute conformer dominates in the ground and the other in the first excited singlet state, and, the signs of the related ECD bands of these two conformers are opposite.

Unusual solute to solvent induced chirality in ROA comes from specific conformer equilibria in the ground and the excited states.

Enantioselective sensing of (S)-Thalidomide in blood plasma with a chiral naphthalene diimide derivative

Fast, simple in use and highly effective voltammetric enantiosensor dedicated for determination of thalidomide (TD) enantiomers (especially towards the toxic (S)-enantiomer) in blood plasma is still desirable. Here we have proven that newly synthesized chiral naphthalene diimide (NDI) derivatives are excellent electroactive materials for TD enantiosensors. The recognition process relies on the specific interaction between the chiral NDI receptor and the thalidomide enantiomer of the opposite configuration. This unique specific interaction between (S)-thalidomide and (R)-NDI derivative counterparts, evident in the DPV voltammograms, was confirmed by molecular modeling. The demonstrated voltammetric enantiosensors are characterized by the low detection limit at the level of μg·L^-1, wide analytical range from 5·10^-4 - 10 mg·L^-1, high selectivity and long lifetime. The results of the recovery rates showed a very good degree of accuracy towards the determination of (S)-thalidomide in the blood samples, so it can be successfully used in the analysis of clinical samples.

Chiral Lanthanide Complexes with l- and d-Alanine: An X-ray and Vibrational Circular Dichroism Study

A whole series of [Ln(H2O)4(Ala)2]26+ dimeric cationic lanthanide complexes with both L- and D-alanine enantiomers was synthesized. The single-crystal X-ray diffraction at 100 and 292 K shows the formation of two types of dimers (I and II) in crystals. Between the dimer centers, the alanine molecules behave as bridging (μ2-O,O'-) and chelating bridging (μ2-O,O,O'-) ligands. The first type of bridge is present in dimers I, while both bridge forms can be observed in dimers II. The IR and vibrational circular dichroism (VCD) spectra of all L- and D-alanine complexes were registered in the 1750-1250 cm-1 range as KBr pellets. Despite all the studied complexes are exhibiting similar crystal structures, the spectra reveal correlations or trends with the Ln-O1 distances which exemplify the lanthanide contraction effect in the IR spectra. This is especially true for the positions and intensities of some IR bands. Unexpectedly, the ν(C=O) VCD bands are quite intense and their composed shapes reveal the inequivalence of the C=O vibrators in the unit cell which vary with the lanthanide. Unlike in the IR spectra, the ν(C=O) VCD band positions are only weakly correlated with the change of Ln and the VCD intensities at most show some trends. Nevertheless, this is the first observation of the lanthanide contraction effect in the VCD spectra. Generally, for the heavier lanthanides (Ln: Dy-Lu), the VCD band maxima are very close to each other and the mirror reflection of the band of two enantiomers is usually better than that of the lighter Lns. DFT calculations show that the higher the multiplicity the higher the stability of the system. Actually, the molecular geometry in crystals (at 100 K) is well predicted based on the highest-spin structures. Also, the simulated IR and VCD spectra strongly depend on the Ln electron configuration but the best overall agreement was reached for the Lu complex, which is the only system with a fully filled f-shell.

Revision of the Structure of N, O-Diacetylsolasodine. Unusual Epimerization at the Spiro Carbon Atom during Acetylation of Solasodine

The steroidal alkaloid solasodine (1) undergoes inversion of configuration at the C-22 spiro atom when treated with acetic anhydride-pyridine at ambient temperature. The basic solvolysis of the N, O-diacetyl derivative (2) reverses the reaction, yielding the starting solasodine (1). The mechanisms of both processes (acetylation and deacetylation) were studied in terms of possible reaction pathways.

Solvation of 2-(hydroxymethyl)-2,5,7,8-tetramethyl-chroman-6-ol revealed by circular dichroism: a case of chromane helicity rule breaking

The primary goal of this work is to clarify why 2-(hydroxymethyl)-2,5,7,8-tetramethyl-chroman-6-ol {(S)-TMChM} deviates from the chromane helicity rule under solvent change. The rule, applicable to determining the absolute configuration of molecules containing the chromane chromophore, binds the sign of the 1Lb Cotton effect (CE) with the helicity of the dihydropyran ring. In case of TMChM, however, this CE exhibits extreme solvent dependence: it is negative in non-coordinating solvents and positive in coordinating ones, irrespective of the helicity of the heterocyclic ring. TD-DFT calculations using PCM and hybrid solvation models were performed to explain origin of this phenomenon. It turned out that the 1Lb CE sign directly depends on the position of the phenolic OH group at carbon atom C6 (OHC6). In the absence of interactions with solvents (as in CCl4 or nC6H14) or when a solvent plays proton donor role (as in CHCl3), the OHC6 lies in the phenyl plane and the 1Lb CE sign follows the P/M helicity rule. In contrast, in proton acceptor solvents, like DMSO, CH3OH or CH3CN, the OHC6 group is deflected from the phenyl plane, and the 1Lb CE sign of individual (S)-TMChM conformers depends on the sector in which the OHC6 is located. Thus, in solution, the 1Lb CE sign is an average over different orientations of the OHC6 group and can be positive (as in DMSO and CH3OH) or negative (as in CH3CN) which means that it does not follow the chromane helicity rule. The impact of OHC6 on the 1Lb CE sign and thus the conclusions for the stereochemistry of chromans are demonstrated here for the first time. Additionally, a comparison of experimental and simulated ECD spectra, supported by VCD data, allowed to determine the geometry of intermolecular clusters formed in different solvents.

Exploring the Sponge Consortium Plakortis symbiotica-Xestospongia deweerdtae as a Potential Source of Antimicrobial Compounds and Probing the Pharmacophore for Antituberculosis Activity of Smenothiazole A by Diverted Total Synthesis

Fractionation of the bioactive CHCl₃-MeOH (1:1) extracts obtained from two collections of the sponge consortium Plakortis symbiotica-Xestospongia deweerdtae from Puerto Rico provided two new plakinidone analogues, designated as plakinidone B (2) and plakinidone C (3), as well as the known plakinidone (1), plakortolide F (4), and smenothiazole A (5). The structures of 1-5 were characterized on the basis of 1D and 2D NMR spectroscopic, IR, UV, and HRMS analysis. The absolute configurations of plakinidones 2 and 3 were established through chemical correlation methods, VCD/ECD experiments, and spectroscopic data comparisons. When assayed in vitro against Mycobacterium tuberculosis H₃₇Rv, none of the plakinidones 1-3 displayed significant activity, whereas smenothiazole A (5) was the most active compound, exhibiting an MIC value of 4.1 μg/mL. Synthesis and subsequent biological screening of 8, a dechlorinated version of smenothiazole A, revealed that the chlorine atom in 5 is indispensable for anti-TB activity.

Modelling the matrix shift on the vibrational frequency of ThO by DFT-D3 calculations

Benchmark calculations with a goal to find dispersion-corrected DFT-D3 methods suitable for a reliable estimation of matrix shifts on the vibrational frequency were carried out on the ThO molecule in three rare gas (Rg = Ne, Ar, and Kr) matrices. The matrices were modelled by the explicit approach, in which a single and a double shell of Rg atoms around ThO was considered. The selection of exchange-correlation functionals was based on test calculations on triatomic ThO⋯Rg models. The B3LYP, PBE0, CAM-B3LYP, and LC-ωPBE functionals were found to be the best suited for the estimation of matrix shifts. The single shell of Rg's around ThO accounted for a major part of the shifts; the addition of a second Rg shell resulted only in a minor improvement. Continuum solvation models considerably overestimated the effect of Rg matrices both when the whole matrix was treated by the model and when the first shell was treated explicitly and the rest with a continuum solvation model.

Reassignment of the absolute configuration of plakinidone from the sponge consortium Plakortis halichondrioides-Xestospongia deweerdtae using a combination of synthesis and a chiroptical approach

Recent work by Wu et al. in connection with the first synthesis of the marine natural product plakinidone revealed that the most salient feature of its purported structure, a six-membered perlactone moiety, was in fact a five-membered lactone, i.e. a 3-methyl-4-hydroxy-2(5H)-furanone or tetronic acid ring. With the planar structure of plakinidone confidently revised, we undertook a new investigation to unambiguously establish its absolute configuration. Upon preparing two stable derivatives 1 and 5 from a sample of naturally occurring plakinidone extracted from the sponge association Plakortis halichondriodes-Xetospongia deweerdtae, the absolute configuration was assigned by synthesis and vibrational and electronic circular dichroism (VCD and ECD) measurements in combination with density functional theory calculations at the B3LYP/aug-cc-pVDZ/PCM(CH3CN) level of theory. Our combined efforts and the agreement between the experimental and calculated VCD/ECD spectra of 1 revealed that the absolute configuration of plakinidone was in fact (11S,17R) and not the formerly reported (11S,17S) diastereomer assigned by Wu et al.. Therefore, we propose that natural plakinidone is accurately represented by structure 12.

Oxidation of Olefins with Benzeneseleninic Anhydride in the Presence of TMSOTf

A new oxidizing system for olefins, consisting of benzeneseleninic anhydride and trimethylsilyl triflate, was studied. The highly reactive benzeneseleninyl cation is presumably formed under these conditions. It has been shown that different products are formed with this species depending on the specific structure of olefin. The 1,1-disubstituted olefins afforded mostly α,β-unsaturated carbonyl compounds. The sterically encumbered tri- or tetrasubstituted olefins yielded 1,2- or 1,4-dihydroxylated products, presumably via four-membered cyclic intermediates.

On stability, chirality measures, and theoretical VCD spectra of the chiral C58X2 fullerenes (X = N, B)

The stability of all 23 C(58)N(2) and C(58)B(2) heterofullerenes in the singlet and triplet states was determined at the B3LYP/6-31G** level. In equilibrium mixture the achiral (1,4) C(58)N(2) isomer would be populated in ca. 95.8%, the chiral (1,16) one in ca. 3.3%, and the achiral (1,4) C(58)B(2) in 100%, whereas all triplet state isomers are less stable. Fourteen out of 23 C(58)X(2) are chiral. Four different chirality measures were calculated by our own CHIMEA program: pure geometrical, labeled, mass, and charge. Intercorrelations between the measures for all chiral compounds indicate that the pure geometrical chirality measure is unstable and should not be used in QSAR predictions of the other molecular properties, while the labeled and mass-weighted ones are promising QSAR descriptors. For each chiral C(58)N(2) molecule, some very strong VCD bands, of intensity comparable with that in the IR spectra, can serve in identification and characterization of the isomers.

Phenylisoserine in the gas-phase and water: Ab initio studies on neutral and zwitterion conformers

The conformational landscape of phenylisoserine (PhIS) was studied. Trial structures were generated by allowing for all combinations of single-bond rotamers. Based on the B3LYP/aug-cc-pVDZ calculations 54 conformers were found to be stable in the gas phase. The six most stable conformers were further optimized at the B3LYP/aug-cc-pVTZ and MP2/aug-cc-pVDZ levels for which characteristic intramolecular hydrogen bond types were classified. To estimate the influence of water on PhIS conformation, the IEF-PCM/B3LYP/aug-cc-pVDZ calculations were carried out and showed 51 neutral and six zwitterionic conformers to be stable in water solution. According to DFT calculations, the conformer equilibrium in the gas phase is dominated by one conformer, whereas the MP2 calculations suggest three PhIS structures to be significantly populated. Comparison of DFT and MP2 energies of all 57 structures stable in water indicates that, in practice, one zwitterionic and one neutral conformer determine the equilibrium in water. Based on the AIM calculations, we found that for the neutral conformers in vacuum and in water, d(H...B) is linearly correlated with Laplacian at the H-bond critical point.

Figure

Density Functional Theory Study on Vibrational Circular Dichroism as a Tool for Analysis of Intermolecular Systems: (1:1) Cysteine−Water Complex Conformations

This paper presents a discussion of the interaction energies for selected conformers of chiral l-cysteine and their (1:1) complexes with water at the B3LYP/aug-cc-pVDZ level. From among more than forty calculated 1:1 complexes three groups of complexes were singled out and examined by the B3LYP/aug-cc-pVDZ calculated vibrational circular dichroism (VCD) spectra. On the basis of analysis of the nu(OmicronEta) and nu(NuEta) and beta(OH2) and beta(NH2) ranges, the VCD spectra were found to be sensitive to conformational changes and water arrangement in cysteine complexes, and to be especially useful for discriminating between different chiral forms of intermolecular hydrogen-bonding complexes. In particular, we show that the VCD modes of an achiral water molecule after complex formation acquire significant rotational strengths whose signs change in line with the geometry of the complex. Moreover, for some water arrangements the VCD spectra can be sensitive to water-wagging conformers and, in temperatures low enough, the intensive nu(OmicronEtaWfree) and beta(H2O) VCD bands may be sufficiently separated to be splitted into pair of oppositely directed bands.

Ar-matrix IR spectra of 5-halouracils interpreted by means of DFT calculations

The infrared low-temperature Ar-matrix spectra of 5-halouracils and unsubstituted uracil were measured and interpreted in terms of the spectra calculated at the DFT/B3PW91/6-311G level followed by a potential energy distribution (PED) analysis. For the PED analysis, the sets of halouracil mode definitions were constructed so that dissimilarities in the interpretations of the different spectra were minimized. Anharmonic frequency calculations enabled more light to be shed on the Fermi resonance (FR) phenomena occurring in the nu(C=O) stretching vibrations region. For each halouracil vibrational spectrum, several FRs manifest themselves in the nu(C=O) stretching vibrations region. We show that the most frequent components participating in these resonances are the nu(C(4)=O(10)) frequency, a beta(N-H) mode frequency, and a beta(C=O) mode frequency. The experimental nu(N-H) frequencies are reproduced by the calculated anharmonic frequencies better than by the scaled harmonic ones, and the nu(C=O) frequencies respond in the opposite manner. The experimental frequencies located below 1500 cm(-1) are reproduced equally well by the two kinds of calculations.

Theoretical Prediction and the First IR Matrix Observation of SeveralL-Cysteine Molecule Conformers

For the first time the argon-matrix low-temperature IR spectra of cysteine are recorded. They reveal a quite complicated spectral pattern, which can also be reproduced in the N2 matrix. Assignment of the observed spectra is undertaken on the basis of comparison of the experimental and calculated B3LYP/aug-cc-pVDZ anharmonic IR spectra. The presence of at least three, and possibly even six or more, cysteine conformers with and without intramolecular hydrogen bonding is confirmed. On the basis of the calculated vibrational circular dichroism spectra, we predict this technique to be more distinctive for conformers than IR absorption is.

An ab Initio Study on the Allene−Isocyanic Acid and Ketene−Vinylimine [2 + 2] Cycloaddition Reaction Paths

The reaction paths of [2 + 2] cycloadditions of allene (H2C=C=CH2) to isocyanic acid (HN=C=O) and ketene (H2C=C=O) to vinylimine (H2C=C=NH), leading to all the possible 14 four-membered ring molecules, were investigated by the MP2/aug-cc-pVDZ method. In the two considered reactions, the 2-azetidinone (beta-lactam) ring compounds were predicted to be the most stable thermodynamically in the absence of an environment. Although 4-methylene-2-azetidinone is the most stable product of the ketene-vinylimine cycloaddition, its activation barrier is higher than that for 4-methylene-2-iminooxetane by ca. 6 kcal/mol. Therefore, the latter product can be obtained owing to kinetic control. The activation barriers in the allene-isocyanic acid reactions are quite high, 50-70 kcal/mol, whereas in the course of the ketene-vinylimine cycloaddition they are equal to ca. 30-55 kcal/mol. All the reactions studied were found to be concerted and mostly asynchronous. Simulation of the solvent environment (toluene, tetrahydrofuran, acetonitrile, and water) by using Tomasi's polarized continuum model with the integral equation formalism (IEF-PCM) method showed the allene-isocyanic reactions remained concerted, yet the activation barriers were somewhat higher than those in the gas phase, whereas the ketene-vinylimine reactions became stepwise. The larger the solvent dielectric constant, the lower the activation barriers found. The lowest-energy pathways in the gas phase and in solvent were confirmed by intrinsic reaction coordinate (IRC) calculations. The atoms in molecules (AIM) analysis of the electron density distribution in the transition-state (TS) structures allowed us to distinguish pericyclic from pseudopericyclic from nonplanar-pseudopericyclic types of reactions.

DFT study of vibrational circular dichroismspectra ofd-lactic acid–water complexes

This paper presents a discussion of the interaction energies, conformations, vibrational absorption (VA, harmonic and anharmonic) and vibrational circular dichroism (VCD) spectra for conformers of monomeric chiral d(-)-lactic acid and their complexes with water at the DFT(B3LYP)/aug-cc-pVDZ and DFT(B3LYP)/aug-cc-pVTZ levels. A detailed analysis has been performed principally for the two most stable complexes with water, differing by lactic acid conformation. The VCD spectra were found to be sensitive to conformational changes of both free and complexed molecules, and to be especially useful for discriminating between different chiral forms of intermolecular hydrogen bonding complexes. In particular, we show that the VCD modes of an achiral water molecule after complex formation acquire significant rotational strengths whose signs change in line with the geometry of the complex. Using the theoretical prediction, we demonstrate that the VCD technique can be used as a powerful tool for structural investigation of intermolecular interactions of chiral molecules and can yield information complementary to data obtained through other molecular spectroscopy methods.

Interactions in open-shell clusters: ab initio study of pre-reactive complex O(3P) + HCl

Van der Waals interactions between the ground-state triplet O(3P) atom and the closed-shell HCl molecule are investigated in the pre-reactive region. Three adiabatic (two of A'' symmetry and one of A' symmetry) and four non-relativistic diabatic potential energy surfaces are obtained by combining a restricted open-shell coupled cluster approach with the multireference configuration interaction method. The lower A'' adiabatic potential surface has a single minimum (D(e) = 589 cm(-1)) for a linear O...HCl configuration. The upper A'' potential has a weak (D(e) = 65 cm(-1)) minimum for a linear HCl...O configuration. The A' adiabatic potential has a weak (124 cm(-1)) T-shaped minimum. Adiabatic potentials intersect once in the O...HCl linear configuration and twice in the linear HCl...O geometry. The role of electrostatic interactions in shaping these potentials is discussed. The effects of spin-orbit coupling on this interaction are also investigated assuming a constant value of the SO parameter.

Reaction Paths of the [2 + 2] Cycloaddition of XCY Molecules (X, Y = S or O or CH2). Ab Initio Study

The reaction paths of [2+2] cycloaddition of the X=C=Y cumulenes were modeled at the MP2/aug-cc-pVDZ level. Cycloadditions of allene and CO2, CS2, or OCS lead in part to the same four-membered products as dimerizations of either ketene or thioketene or addition of ketene and thioketene, respectively. All the reactions studied are concerted and mostly asynchronous. The majority of the allene cycloadditions studied are endoergic and proceed with much higher activation barriers than do the alternative (thio)ketene additions. In comparison with the energy of the substrates, the four-membered cycles incorporating S-atoms are stabilized more than the analogous structures with O-atoms built into the rings. There are also some products that are thermodynamically disfavored, yet seem to be obtainable thanks to a relatively low barrier of the reaction. The AIM analysis of the electron density distribution in the transition state structures allowed distinguishing pericyclic from pseudopericyclic and nonplanar-pseudopericyclic types of reaction.