Estimation of formamide harmonic and anharmonic modes in the Kohn-Sham limit using the polarization consistent basis sets

Application of the FFLUX Force Field to Molecular Crystals: A Study of Formamide

In this work, we present the first application of the quantum chemical topology force field FFLUX to the solid state. FFLUX utilizes Gaussian process regression machine learning models trained on data from the interacting quantum atom partitioning scheme to predict atomic energies and flexible multipole moments that change with geometry. Here, the ambient (α) and high-pressure (β) polymorphs of formamide are used as test systems and optimized using FFLUX. Optimizing the structures with increasing multipolar ranks indicates that the lattice parameters of the α phase differ by less than 5% to the experimental structure when multipole moments up to the quadrupole are used. These differences are found to be in line with the dispersion-corrected density functional theory. Lattice dynamics calculations are also found to be possible using FFLUX, yielding harmonic phonon spectra comparable to dispersion-corrected DFT while enabling larger supercells to be considered than is typically possible with first-principles calculations. These promising results indicate that FFLUX can be used to accurately determine properties of molecular solids that are difficult to access using DFT, including the structural dynamics, free energies, and properties at finite temperature.

Basis Set Extrapolation of Vibrational Frequencies

We investigate three different approaches for extrapolating harmonic vibrational frequencies to the complete basis set limit, by direct extrapolation of the frequencies, by calculation of frequencies based on extrapolated Hessians, or based on the Hessian from optimization of the extrapolated energy surface. For regular molecules, the three extrapolation procedures yield essentially identical results, but for loosely bound complexes, the frequencies derived from extrapolated Hessians yield unpredictable behavior. None of the basis set extrapolations, however, provide any significant improvement over the results upon which the extrapolation is based.

Effect of unsaturated substituents in the reaction of Criegee intermediates with water vapor

Criegee intermediates (CIs), formed in the reactions of unsaturated hydrocarbons with ozone, are very reactive carbonyl oxides and have recently been suggested as important oxidants in the atmosphere. In this work, we studied the substituent effect on the water monomer and dimer reaction with CIs which include up to three carbon atoms at the QCISD(T)/CBS//B3LYP/6-311+G(2d,2p) level. Our calculation showed that for saturated CIs with a hydrogen atom on the same side as the terminal oxygen atom, the reaction with water vapor would likely dominate the removal processes of these CIs in the atmosphere. On the other hand, for unsaturated CIs, the reactivity toward water vapor decreases compared to the saturated species allowing them to survive in humid atmospheric environments. We also evaluated the kinetic isotope effect in the reaction between CI and water vapor by performing calculations with deuterated water. We found that tunneling is not important and the kinetic isotope effect mainly comes from the difference in the zero point energy between water and deuterated water.

Method and basis set dependence of the NICS indexes of aromaticity for benzene

The role of theory level in prediction of benzene magnetic indexes of aromaticity is analysed and compared with calculated nuclear magnetic shieldings of ³ He used as NMR probe. Three closely related nucleus-independent chemical shift (NICS) based indexes were calculated for benzene at SCF-HF, MP2, and DFT levels of theory and the impact of basis set on these quantities was studied. The changes of benzene NICS(0), NICS(1), and NICS(1)zz parameters calculated using SCF-HF, MP2 and several density functionals were within 1 to 3 ppm. Similar deviations between magnetic indexes of aromaticity were observed for values calculated with selected basis sets. Only very small effect of polar solvent on benzene aromaticity was predicted. The ³ He nuclear magnetic isotropic shielding (σ) and its zz-components (σ_zz ) of helium atom approaching the centre of benzene ring from above produced similar curves versus benzene-He distance to NICS parameters calculated for similarly moving Bq ghost atom. We also propose an experimental verification of NICS calculations by designing the ³ He NMR measurement for benzene saturated with helium gas or in low temperature matrices.

Calculation of 15N NMR chemical shifts: Recent advances and perspectives

Recent advances in computation of ¹⁵N NMR chemical shifts are reviewed, concentrating mainly on practical aspects of computational protocols and accuracy factors. The review includes the discussion of the level of theory, the choice of density functionals and basis sets together with taking into account solvent effects, rovibrational corrections and relativistic effects. Computational aspects of ¹⁵N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.

DFT studies on the structural and vibrational properties of polyenes

Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C = C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C = C) and ν(C-C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon-carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C = C units with both ends substituted by tert-butyl groups. The corrected ν(C = C) and ν(C-C) wavenumbers for all-trans molecules were used to estimate the presence of 9 - 12 C = C units in all-trans polyene pigment in red coral.

Predicting the structure and vibrational frequencies of ethylene using harmonic and anharmonic approaches at the Kohn-Sham complete basis set limit

In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested.

The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies.

Graphical Abstract

Estimating the carbonyl anharmonic vibrational frequency from affordable harmonic frequency calculations

A linear correlation between harmonic and anharmonic frequencies of water calculated at B3LYP level of theory was observed with a number of basis sets. Similar relationships were found in both the gas phase and solution for several small molecules. The best correlation was found for C = O stretch mode in formaldehyde, formamide and N-methylacetamide. The average difference between B3LYP harmonic and anharmonic ν(C = O) frequencies calculated with several basis sets in these molecules was 30 cm(-1). The ad hoc correction of -30 cm(-1), added to harmonic frequencies of two different carbonyl groups present in a structure of a larger molecule was tested as a fast way of predicting anharmonic frequencies without elaborated calculations. The proposed approach was tested successfully on a larger molecule of E and Z isomers of N-acetyl-α,β-dehydrophenylalanine N',N'-dimethylamide [Ac-(E/Z)-ΔPhe-NMe(2)] and the estimated anharmonic ν(C = O) frequencies were close to directly calculated results.