DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.

Time delay between singly and doubly ionizing wavepackets in laser-driven helium

We present calculations of the time delay between single and double ionization of helium, obtained from full-dimensionality numerical integrations of the helium-laser Schroedinger equation. The notion of a quantum mechanical time delay is defined in terms of the interval between correlated bursts of single and double ionization. Calculations are performed at 390 and 780 nm in laser intensities that range from 2 X 10^14 to 14 X 10^14 W /cm^2. We find results consistent with the rescattering model of double ionization but supporting its classical interpretation only at 780 nm.

Local and semilocal density functional computations for crystals of 1-alkyl-3-methyl-imidazolium salts

The accuracy and reliability of popular density functional approximations for the compounds giving origin to room temperature ionic liquids have been assessed by computing the T=0 K crystal structure of several 1-alkyl-3-methyl-imidazolium salts. Two prototypical exchange-correlation approximations have been considered, i.e., the local density approximation (LDA) and one gradient corrected scheme [PBE-GGA, Phys. Rev. Lett. 77, 3865 (1996)]. Comparison with low-temperature x-ray diffraction data shows that the equilibrium volume predicted by either approximations is affected by large errors, nearly equal in magnitude (~10%), and of opposite sign. In both cases the error can be traced to a poor description of the intermolecular interactions, while the intramolecular structure is fairly well reproduced by LDA and PBE-GGA. The PBE-GGA optimization of atomic positions within the experimental unit cell provides results in good agreement with the x-ray structure. The correct system volume can also be restored by supplementing PBE-GGA with empirical dispersion terms reproducing the r-6 attractive tail of the van der Waals interactions.

On the origin of the forward peak and backward oscillations in the F + H2(v = 0) to HF(v' = 2) + H reaction

We present a semiclassical complex angular momentum (CAM) analysis of the forward scattering peak which occurs at a translational collision energy around 32 meV in the quantum mechanical calculations for the F + H2(v = 0, j = 0) ? HF(v' = 2, j' = 0) + H reaction on the Stark–Werner potential energy surface. The semiclassical CAM theory is modified to cover the forward and backward scattering angles. The peak is shown to result from constructive/destructive interference of the two Regge states associated with two resonances, one in the transition state region and the other in the exit channel van der Waals well. In addition, we demonstrate that the oscillations in the energy dependence of the backward differential cross section are caused by the interference between the direct backward scattering and the decay of the two resonance complexes returning to the backward direction after one full rotation.

Isolated vibrational wavepackets in D-2(+): Defining superposition conditions and wavepacket distinguishability

Tunnel ionization of room-temperature D-2 in an ultrashort (12 femtosecond) near infrared (800 nm) pump laser pulse excites a vibrational wavepacket in the D-2(+) ions; a rotational wavepacket is also excited in residual D-2 molecules. Both wavepacket types are collapsed a variable time later by an ultrashort probe pulse. We isolate the vibrational wavepacket and quantify its evolution dynamics through theoretical comparison. Requirements for quantum computation (initial coherence and quantum state retrieval) are studied using this well-defined (small number of initial states at room temperature, initial wavepacket spatially localized) single-electron molecular prototype by temporally stretching the pump and probe pulses.

Absolute single and multiple charge exchange cross sections for highly charged C, O, and Ne ions on H2O, CO, and CO2

Reported herein are measured absolute single, double, and triple charge exchange (CE) cross sections for the highly charged ions (HCIs) Cq+ (q=5,6), Oq+ (q=6,7,8), and Neq+ (q=7,8) colliding with the molecular species H2O, CO, and CO2. Present data can be applied to interpreting observations of x-ray emissions from comets as they interact with the solar wind. As such, the ion impact energies of 7.0q keV (1.62–3.06 keV/amu) are representative of the fast solar wind, and data at 1.5q keV for O6+ (0.56 keV/amu) on CO and CO2 and 3.5q keV for O5+ (1.09 keV/amu) on CO provide checks of the energy dependence of the cross sections at intermediate and typical slow solar wind velocities. The HCIs are generated within a 14 GHz electron cyclotron resonance ion source. Absolute CE measurements are made using a retarding potential energy analyzer, with measurement of the target gas cell pressure and incident and final ion currents. Trends in the cross sections are discussed in light of the classical overbarrier model (OBM), extended OBM, and with recent results of the classical trajectory Monte Carlo theory.