Magnitude of zero-point vibrational corrections to optical rotation in rigid organic molecules: a time-dependent density functional study

The zero-point vibrational corrections (ZPVCs) to the optical rotation of 22 rigid organic molecules have been calculated using time-dependent density functional theory with the B3LYP hybrid functional. We outline an implementation for calculating ZPVCs that can be applied with a variety of quantum chemistry programs and methods. It is shown that the ZPVCs to optical rotation have a wide range of values and can be quite significant depending on the molecule. On average, it has been determined that vibrational corrections can account for about 20% of the optical rotation for the equilibrium value. It is also concluded that vibrational effects alone cannot be the only factor in improving the calculated values of optical rotation with respect to experimental data measured in solution.     

Temperature dependence of the optical rotation in six bicyclic organic molecules calculated by vibrational averaging.

The vibrational corrections and the temperature dependence of the specific rotation of six rigid organic molecules (alpha-pinene, beta-pinene, cis-pinane, camphene, camphor, and fenchone) were calculated at three wavelengths using hybrid time-dependent density functional theory (TDDFT). A technique for calculating the temperature dependence of the vibrational average of a molecular property has been applied to obtain the specific rotation of the molecules as a function of temperature. For cases in which accurate equilibrium optical rotations can be obtained as a "base value," and for which there is little effect from solvation, accurate predictions of the trends in the temperature-dependence of the specific rotations can be calculated. For other cases, the method can be used to extract purely vibrational contributions to the overall temperature dependence of optical rotation.     

Gas-phase electron diffraction,vibrational spectroscopy,and quantum chemical studies of the molecular structure of 3,3-dimethyl-3-silathiane

The 3,3-dimethyl-3-silathiane molecule was studied by gas-phase electron diffraction and vibrational spectroscopy. The initial geometrical parameters and the force field were calculated by the B3LYP/6-311+G** method; the vibrational amplitudes of atomic pairs and vibrational corrections were calculated using the scaled B3LYP/6-311+G** force field. The molecular conformation was found to be a distorted chair with structural parameters close to the expected ones.     

Vibrational spectral analysis and first hyperpolarizability studies of 1-bromonaphthalene based on ab initio and DFT methods

In this work, the experimental and theoretical vibrational spectra of 1-bromonaphthalene (1-BN) were studied. FTIR and FT Raman spectra were recorded in the region 4000-400 cm(-1) and 3500-50 cm(-1), respectively. The structural and spectroscopic data of the molecule in the ground state were calculated by using ab initio Hartree-Fock and density functional method (B3LYP) with the 6-311++G(d,p) basis set. The vibrational frequencies were calculated and scaled values were compared with experimental FTIR and FT Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The optimized geometric parameters were calculated. The predicted first hyperpolarizability also shows that the molecule might have a reasonably good nonlinear optical (NLO) behaviour. The calculated HOMO-LUMO energy gap reveals that charge transfer occurs within the molecule.     

Equilibrium conformation of dimethylaminodichlorophosphine molecule: the use of the results of quantum-chemical calculations and spectroscopic measurements in the analysis of gas-phase electron diffraction data

Geometric parameters and force fields of two stable isomers of dimethylaminodichlorophosphine molecule, a gauche-conformer with C₁ symmetry (A) and anti-conformer with C _s symmetry (D), resulting from internal rotation about the P—N bond, were calculated in the RHF/6-31G* approximation. Using the scaled quantum-chemical force field for the most stable conformer A, the first reliable interpretation of the vibrational spectra of the light and perdeuterated isotopomers of dimethylaminodichlorophosphine was obtained. The root-mean-square vibrational amplitudes, harmonic and anharmonic vibrational corrections, and centrifugal distortion corrections were also calculated. Structural analysis of electron diffraction data was performed with consideration of nonlinear kinematic effects at the first-order level of perturbation theory. The experimental values of the equilibrium geometric parameters were estimated. The results obtained suggest a nonplanar equilibrium configuration of the amino group in the dimethylaminodichlorophosphine molecule.     

Comparison of time-dependent density-functional theory and coupled cluster theory for the calculation of the optical rotations of chiral molecules

A comparison of the abilities of time-dependent density-functional theory (TDDFT) and coupled cluster (CC) theory to reproduce experimental sodium D-line specific rotations for 13 conformationally rigid organic molecules is reported. The test set includes alkanes, alkenes, and ketones with known absolute configurations. TDDFT calculations make use of gauge-including atomic orbitals and give origin-independent specific rotations. CC rotations are computed using both the origin-independent dipole-velocity and origin-dependent dipole-length representations. The mean absolute deviations of calculated and experimental rotations are of comparable magnitudes for all three methods. The origin-independent DFT and CC methods give the same sign of [alpha]D for every molecule except norbornanone. For every large-rotation ketone and alkene for which DFT and CC yield the incorrect sign as compared to liquid-phase experimental data, the corresponding optical rotatory dispersion (ORD) curve is bisignate, suggesting that the two models cannot reliably reproduce the relative excitation energies and antagonistic rotational strengths of multiple competing electronic states that contribute to the total long-wavelength rotation. Several potential sources of error in the theoretical treatments are considered, including basis set incompleteness, vibrational and temperature effects, electron correlation, and solvent effects.     

Time-dependent density functional theory for calculating origin-independent optical rotation and rotatory strength tensors

An approach to calculate origin-independent electronic chiroptical property tensors using time-dependent density functional theory (TDDFT) and gauge-including atomic orbital (GIAO) basis sets is evaluated. Computations of origin-dependent optical rotation tensors and of rotatory strengths needed to simulate circular dichroism spectra are presented. The optical rotation tensor computations employ solutions of coupled perturbed Kohn-Sham equations for a dynamic electric field and a static magnetic field. Because the magnetic field is time independent, the GIAO treatment is somewhat simplified compared to a previously reported method, at some added computational cost if hybrid functionals are employed. GIAO rotatory strengths are also calculated, using transition density matrices from a standard TDDFT excitation energy module. A new implementation in the NWChem quantum chemistry package is employed for representative computations of origin-invariant chiroptical response tensors for methyloxirane, norbornenone, and the ketosteroid androstadienone. For the steroid molecule the vibrational structure of the CD spectrum is modeled explicitly by using calculated Franck-Condon factors. The agreement with experiment is favorable.     

Coupled cluster and density functional theory studies of the vibrational contribution to the optical rotation of (S)-propylene oxide

In a previous study (Chemical Physics Letters 2005, 401 , 385) we computed the optical rotatory dispersion of (S)-propylene oxide in gas phase and solution using the hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3. Even for the highly correlated CC3 model combined with a flexible basis set, the theoretical gas-phase specific rotation at 355 nm was found to be negative in contrast to the experimental result. We argued that vibrational contributions could be crucial for obtaining a complete understanding of the experimental result. Here, we show that this indeed is the case by using coupled cluster models and density functional theory methods to calculate the vibrational contributions to the gas-phase specific rotation at 355, 589.3, and 633 nm. While density functional theory (B3LYP and SAOP functionals) overestimates the specific rotation at 355 nm by approximately 1 order of magnitude and yields an incorrect sign at 589.3 and 633 nm, the coupled cluster results are in excellent agreement with the experimentally measured optical rotations. We find that all vibrational modes contribute significantly to the optical rotation and that temperature effects must be taken into account.     

FT-IR, UV-vis, 1H and 13C NMR spectra and the equilibrium structure of organic dye molecule disperse red 1 acrylate: a combined experimental and theoretical analysis

This study reports the characterization of disperse red 1 acrylate compound by spectral techniques and quantum chemical calculations. The spectroscopic properties were analyzed by FT-IR, UV-vis, (1)H NMR and (13)C NMR techniques. FT-IR spectrum in solid state was recorded in the region 4000-400 cm(-1). The UV-vis absorption spectrum of the compound that dissolved in methanol was recorded in the range of 200-800 nm. The (1)H and (13)C NMR spectra were recorded in CDCl(3) solution. The structural and spectroscopic data of the molecule in the ground state were calculated using density functional theory (DFT) employing B3LYP exchange correlation and the 6-311++G(d,p) basis set. The vibrational wavenumbers were calculated and scaled values were compared with experimental FT-IR spectrum. A satisfactory consistency between the experimental and theoretical spectra was obtained and it shows that the hybrid DFT method is very useful in predicting accurate vibrational structure, especially for high-frequency region. The complete assignments were performed on the basis of the experimental results and total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. Isotropic chemical shifts were calculated using the gauge-invariant atomic orbital (GIAO) method. A study on the electronic properties were performed by timedependent DFT (TD-DFT) and CIS(D) approach. To investigate non linear optical properties, the electric dipole moment μ, polarizability α, anisotropy of polarizability Δα and molecular first hyperpolarizability β were computed. The linear polarizabilities and first hyperpolarizabilities of the studied molecule indicate that the compound can be a good candidate of nonlinear optical materials.     

Zero-point vibrational effects on proton shieldings: functional-group contributions from ab initio calculations.

We investigate the effects of zero-point vibrational motion on the nuclear magnetic shielding constants of a large number of organic molecules. The vibrational corrections include anharmonic contributions from the potential energy surface and harmonic contributions from the curvature of the property surface. Particular attention is paid to vibrational corrections to hydrogen shielding constants where we show that vibrational corrections may be substantial, ranging from about +0.50 to -0.70 ppm, and thus demonstrating that ignoring these effects may give errors in the chemical shifts by more than 1 ppm in certain extreme cases. These effects can therefore not be neglected when comparing calculated results with experiment, not even for the chemical shifts. However, we also demonstrate that the vibrational corrections to the hydrogen shieldings are to a large extent transferable from one molecule to another. We have tabulated functional vibrational corrections to the hydrogen shieldings, based on results for more than 35 molecules. Unfortunately, no similar transferability has been observed for the vibrational corrections to shielding constants of other nuclei such as carbon, nitrogen, or oxygen.     

Dynamic (hyper)polarizabilities of the sulphur dioxide molecule: coupled cluster calculations including vibrational corrections

In this work we report results for dynamical (hyper)polarizabilities of the sulphur dioxide molecule with inclusion of vibrational corrections. The electronic contributions were computed analytically at the single and double coupled cluster level through response theories for the frequencies 0, 0.0239, 0.0428, 0.0656, 0.0720, and 0.0886 hartree. Contributions of the connected triple excitations to the dynamic electronic properties were also estimated through the multiplicative correction scheme. Vibrational corrections were calculated by means of the perturbation theoretical method. The results obtained show that the zero point vibrational correction is very small for all properties studied while the pure vibrational correction is relevant for the dc-Pockels effect, intensity dependent refractive index, and dc-Kerr effect. For these nonlinear optical processes, the pure vibrational corrections represent approximately 75%, 13%, and 6% of the corresponding electronic contributions for the higher frequencies quoted. The results presented for the polarizability are in good agreement with experimental values available in the literature. For the hyperpolarizabilities we have not obtained experimental results with precision sufficient for comparison.     

Temperature dependence of the optical rotation of fenchone calculated by vibrational averaging

On the basis of vibrational averaging, the temperature dependence of the optical rotation for fenchone has been calculated using TDDFT with the B3LYP hybrid functional at three wavelengths. The results show that very good agreement is obtained between theory and experiment. It is concluded that temperature-dependent vibrational effects are likely to account for much of the observed temperature dependence in optical rotation exhibited by rigid organic molecules in case there is only a weak temperature-dependent interaction with the solvent.     

Vibrational relaxation of oxygen in an argon cage

The vibrational relaxation of oxygen embedded in an argon cage through vibrational to local translation, rotation, and argon phonon modes has been studied using semiclassical procedures. The collision model is based on the trapped molecule undergoing the restricted motions (local translation and hindered rotation) in a cage formed by its twelve nearest argon neighbors in a face-centered-cubic structure. At 85 K in the liquid argon temperature range, the deexcitation probability of O(2)(v=1) is 5.8 x 10(-12) and the relaxation rate constant with the collision frequency from local translation is 23 s(-1). The rate constant decreases to 5.1 s(-1) at 50 K and to 0.016 s(-1) at 10 K in the solid argon temperature range. Transfer of the vibrational energy to local translation, rotation (both hindered and free), and argon phonon modes is the relaxation pathway for the trapped oxygen molecule.     

Solvent effects on Raman optical activity spectra calculated using the polarizable continuum model

The integral equation formulation of the polarizable continuum model (IEFPCM) has been extended to the calculation of solvent effects on vibrational Raman optical activity spectra. Gauge-origin independence of the differential scattering intensities of right and left circularly polarized light is ensured through the use of London atomic orbitals. Density functional theory (DFT) calculations have been carried out for bromochlorofluoromethane, methyloxirane, and epichlorhydrin. The results indicate that solvent effects on the ROA differential scattering intensities can be substantial, and vary in sign and magnitude for different vibrational modes. It is demonstrated that both direct and indirect effects are important in determining the total solvent effects on the ROA differential scattering intensities. Local field effects are shown to be in general small, whereas electronic nonequilibrium solvation has a profound effect on the calculated solvent effects compared to an equilibrium solvation model. For molecules with several conformations, the changes in the relative stability of the different conformers also lead to noticeable changes in the ROA spectra.     

Is styrene planar in liquid phases?

The proton NMR spectra of two (13)C-labeled isotopomers of styrene dissolved in two liquid crystalline solvents have been obtained and analyzed to yield four sets each of 24 dipolar couplings. These couplings were then used to investigate the structure of the ring and the ene fragments of the molecule, and the position of the maximum, phi(0), in the ring-ene bond rotational probability distribution. To do this, the effect on the dipolar couplings of small-amplitude vibrational motion was taken into account using vibrational wave functions calculated by molecular orbital and density functional methods. It is concluded that the NMR data are consistent with the ring fragment, averaged over the ring-ene rotation, planar, while the ene fragment is not. The value of phi(0) is found to be 18.0 degrees +/-0.2 degrees for the two solutions, compared with a value of 27 degrees calculated by the molecular method MP2/6-31G(*).     

Monohydrates of cuprous chloride and argentous chloride: H2O⋅⋅⋅CuCl and H2O⋅⋅⋅AgCl characterized by rotational spectroscopy and ab initio calculations

Pure rotational spectra of the ground vibrational states of ten isotopologues of each of H(2)O???CuCl and H(2)O???AgCl have been measured and analyzed to determine rotational constants and hyperfine coupling constants for each molecule. The molecular structure and spectroscopic parameters determined from the experimental data are presented alongside the results of calculations at the CCSD(T) level. Both experiment and theory are consistent with structures that are nonplanar at equilibrium. The heavy atoms are collinear while the local C(2) axis of the water molecule intersects the axis defined by the heavy atoms at an angle, φ = 40.9(13)° for Cu and φ = 37.4(16)° for Ag. In the zero-point state, each molecule is effectively planar, undergoing rapid inversion between two equivalent structures where φ has equal magnitude but opposite sign. The equilibrium geometry has C(s) symmetry, however. The ab initio calculations confirm that the timescale of this inversion is at least an order of magnitude faster than that of rotation of the molecule in the lowest rotational energy levels. The molecular geometries are rationalized using simple rules that invoke the electrostatic interactions within the complexes. Centrifugal distortion constants, Δ(J) and Δ(JK), nuclear quadrupole coupling constants, χ(aa)(Cu), χ(aa)(Cl), (χ(bb) - χ(cc))(Cu), and (χ(bb) - χ(cc))(Cl), and the nuclear spin-rotation constant of the copper atom, C(bb)(Cu)+C(cc)(Cu), are also presented.     

Chiroptical spectroscopic determination of molecular structures of chiral sulfinamides: t-butanesulfinamide

The absolute configuration of t-butanesulfinamide has been determined as (-)-(S) using three different chiroptical spectroscopic methods, namely, vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and optical rotatory dispersion (ORD). Furthermore, the predominant conformation of this molecule is determined to have S=O and NH2 groups staggered with respect to the three methyl groups and to have amine hydrogens in gauche orientation with respect to S=O. The quality of predictions obtained for vibrational properties, namely, vibrational absorption and VCD, is found to be satisfactory with the B3LYP functional and 6-31G* basis set. However, this basis set is found to be inadequate for obtaining reliable predictions of electronic properties, namely, electronic absorption and ECD, but a larger aug-cc-pVDZ basis set is found to provide satisfactory prediction of electronic properties. t-Butanesulfinamide serves as an example which invalidates the recommendation of using the 6-31G* basis set for molecules that exhibit the same sign for the long-wavelength ECD band and ORD. This molecule also emphasizes the importance of simultaneous investigation of ECD and ORD, and the use of multiple chiroptical spectroscopic methods, for reliable determination of molecular stereochemistry.     

Effects of complex formation on vibrational circular dichroism spectra

The determination of absolute configurations of chiral compounds using VCD is performed by comparing measured vibrational circular dichroism (VCD) spectra with calculated spectra. The process is based on two facts: the two enantiomers have rotational strengths of opposite sign, and the absolute configuration of the molecule used in the calculation is known. However, calculations on isolated molecules very often predict VCD intensities of very different magnitude or even different signs compared to the spectra measured in solution. Therefore, we have carefully analyzed what type of changes are induced by complexation of a solvent molecule to a solute. In the theoretical example of benzoyl-benzoic acid (in a particular chiral conformation) hydrogen bonded to the achiral NH3, we distinguish six cases, ranging from no or very small changes in the rotational strengths of solute modes (case A) to changes of sign of rotational strengths (case B), changes in magnitude (case C), nonzero rotational strengths for modes of the achiral solvent ("transfer of chirality", case D), large frequency shifts accompanied by giant enhancements of the IR and VCD intensities of modes involved in hydrogen bonding (case E), and emergence of new peaks (case F). In this work, all of these situations will be discussed and their origin will be elucidated. On the basis of our analysis, we advocate that codes for VCD rotational strength calculation should output for each mode i the angle xi(i) between the electric and magnetic transition dipole moments because only "robust modes" with xi far from 90 degrees should be used for the determination of the absolute configuration.     

Computational examination on monomeric,dimeric, trimeric structural and vibrational interactions,AIM, Hirshfeld,IGM and oxygenated solvent effect on optical properties for pyridine N-oxide

This work reports DFT (Density functional theory) and topological examination by means of AIM (Atom in molecule) theory, Laplacian electron density (ED), Electron –Localization function (ELF) and Hirshfeld surface were executed monomer, dimer and trimer structure of PNO molecule. These structures were optimized and Structural parameters like bond – length and Angles are compared by using B3LYP/6-311G++(d,p) basis set. Computed vibrational – frequencies and corresponding PED (Potential energy – Distribution) were also compared. Polarizability, Dipole moment and the electronic properties were calculated for all three cases of PNO molecule. Hirshfeld surface investigation has been executed to achieve the C–H?O/N–O?H type of hydrogen -bond intermolecular interactions. Moreover, Reactive site, inter/intra interactions are calculated, the optical behaviors for monomer structure of PNO were investigate by using UV–Vis spectroscopy in oxygenated solvents.