Search

A detailed analysis of the electronic structure and charge distribution around the trigonal site of Li-doped polyacetylene is reported using finite chain and periodic descriptions of the polymer. Atoms-in-molecules (AIM) analysis is done to characterize the nature of the bond between Li and the polymer backbone through the location of the bond critical points and computation of the total charge on the atomic basins around the doping site. We find that the Li atom donates practically one electron to the π-system, in accordance with the classical Su-Schriffer-Heeger model. However, despite that the Li atom is equidistant from the three closest C atoms in the geometric soliton, a single Li-C bond critical point is found. The AIM quantitative analysis of the electronic density reveals that the Li(+) ion is immersed into the polymer π-cloud in a way that resembles a metallic bonding interaction.

The structure and cohesive energy of crystalline urea have been investigated at the ab initio level of calculation. The performance of different Hamiltonians in dealing with a hydrogen-bonded molecular crystal as crystalline urea is assessed. Detailed calculations carried out by adopting both HF and some of the most popular DFT methods in solid-state chemistry are reported. Local, gradient-corrected, and hybrid functionals have been adopted: SVWN, PW91, PBE, B3LYP, and PBE0. First, a 6-31G(d,p) basis set has been adopted, and then the basis set dependence of computed results has been investigated at the B3LYP level. All calculations were carried out by using a development version of the periodic ab initio code CRYSTAL06, which allows full optimization of lattice parameters and atomic coordinates. With the 6-31G(d,p) basis set, structural features are well reproduced by hybrid methods and GGA. LDA gives lattice parameters and hydrogen-bond distances that are too small relative to experiment, while at the HF level the opposite trend is observed. Results show that hybrid methods are more accurate than HF and both LDA and GGA functionals, with a trend in the computed properties similar to that of hydrogen-bonded molecular complexes. When BSSE and ZPE are taken into account, all methods, except LDA, give computed cohesive energies that are underestimated with respect to the experimental sublimation enthalpy. Dispersion energy, not properly taken into account by DFT methods, plays a crucial role. Such a deficiency also affects dramatically the computed crystalline structure, especially when large basis sets are adopted. We show that this is an artifact due to the BSSE. Indeed, with small basis sets the BSSE gives an extra-binding that compensates for the missing dispersion forces, thus yielding structures in fortuitous agreement with experiment.

Details on the mechanism of HF catalyzed isobutylene-isobutane alkylation were investigated. On the basis of available experimental data and high-level quantum chemical calculations, a detailed reaction mechanism is proposed taking into account solvation effects of the medium. On the basis of our computational results, we explain why the density of the liquid media and stirring rates are the most important parameters to achieve maximum yield of alkylate, in agreement with experimental findings. The ab initio Car-Parrinello molecular dynamics calculations show that isobutylene is irreversibly protonated in the liquid HF medium at higher densities, leading to the ion pair formation, which is shown to be a minimum on the potential energy surface after optimization using periodic boundary conditions. The HF medium solvates preferentially the fluoride anion, which is found as solvated [FHF](-) or solvated F(-.)(HF)(3). On the other hand, the tert-butyl cation is weakly solvated, where the closest HF molecules appear at a distance of about 2.9 Angstrom with the fluorine termination of an HF chain.

The vibrational spectrum of Mg(3)Al(2)Si(3)O(12) pyrope is calculated at the Gamma point by using the periodic ab initio CRYSTAL program that adopts an all-electron Gaussian-type basis set and the B3LYP Hamiltonian. The full set of frequencies (17 IR active, 25 RAMAN active, 55 silent modes) is calculated. The effect of the basis set and of the computational parameters on the calculated frequencies is discussed. It is shown that the mean absolute difference with respect to the experimental IR and RAMAN data is as small as 6 and 8 cm(-1), respectively. The IR and RAMAN modes are fully characterized by various tools such as isotopic substitution, direct inspection of the eigenvectors, and graphical representation. The present calculation permits to clarify some of the assignment and interpretation problems raised by experiment and previous simulations with force fields.