A. Ghysels

Theoretical and experimental IR data are combined to gain insight into the methanol and ethanol conversion over an acidic H-SAPO-34 catalyst. The theoretical simulations use a large finite cluster and the initial physisorption energy of both alcohols is calculated. Dispersive contributions turn out to be vital and ethanol adsorbs stronger than methanol with approximately 14 kJ mol(-1). Calculated IR spectra of the alcohols and of formed aromatic cations upon conversion are also analyzed and support the peak assignment of the experimental in situ DRIFT spectra, in particular for the C-H and C=C regions. Theoretical IR spectra of the gas phase compounds are compared with those of the molecules loaded in a SAPO cluster and the observed shifts of the peak positions are discussed. To get a better understanding of these framework-guest interactions, a new theoretical procedure is proposed based on a normal mode analysis. A cumulative overlap function is defined and enables the characterization of individual peaks as well as induced frequency shifts upon adsorption. (C) 2010 Elsevier B. V. All rights reserved.

In an earlier work, the authors developed a new method, the mobile block Hessian (MBH) approach, to accurately calculate vibrational modes for partially optimized molecular structures [ J. Chem. Phys. 2007, 126 (22), 224102.]. It is based on the introduction of blocks, consisting of groups of atoms, that can move as rigid bodies. The internal geometry of the blocks need not correspond to an overall optimization state of the total molecular structure. The standard MBH approach considers free blocks with six degrees of freedom. In the extended MBH approach introduced herein, the blocks can be connected by one or two adjoining atoms, which further reduces the number of degrees of freedom. The new approach paves the way for the normal-mode analysis of biomolecules such as proteins. It rests on the hypothesis that low-frequency modes of proteins can be described as pure rigid-body motions of blocks of consecutive amino acid residues. The method is validated for a series of small molecules and further applied to alanine dipeptide as a prototype to describe vibrational interactions between two peptide units; to crambin, a small protein with 46 amino acid residues; and to ICE/caspase-1, which contains 518 amino acid residues.

Silicalite-1 zeolite exhibits a characteristic pentasil framework vibration around 540−550 cm−1. In the initial stages of zeolite synthesis, however, this band is observed at much higher wavenumbers: literature shows this vibration to depend on particle size and to shift over 100 cm−1 with increasing condensation. In this work, the pentasil vibration frequency was derived from theoretical molecular dynamics simulations to obtain the correct IR band assignments for important nanoparticles. The IR spectroscopic fingerprint of oligomeric five-ring containing precursors proposed in the literature was computed and compared with experimental data. Our theoretical results show that, while isolated five-membered rings show characteristic vibrational bands around 650 cm−1, the combination of five-membered rings in the full MFI-type structure readily generates the bathochromic shift to the typical pentasil vibration around 550 cm−1. As opposed to what was previously believed, the IR band does not shift gradually as nanoparticle size increases, but it is highly dependent on the specific way structural units are added. The most important feature is the appearance of an additional band when double five-membered rings are included, which allows for a clear distinction between the key stages of early zeolite nucleation. Furthermore, the combination of the simulated spectra with the experimental observation of this spectral feature in nanoparticles extracted from silicalite-1 clear solutions supports their structured nature. The theoretical insights on the dependency of pentasil vibrations with the degree of condensation offer valuable support toward future investigations on the genesis of a zeolite crystal.

Partial optimization is a useful technique to reduce the computational load in simulations of extended systems. In such nonequilibrium structures, the accurate calculation of localized vibrational modes can be troublesome, since the standard normal mode analysis becomes inappropriate. In a previous paper [ A. Ghysels et al., J. Chem. Phys. 126, 224102 (2007) ], the mobile block Hessian (MBH) approach was presented to deal with the vibrational analysis in partially optimized systems. In the MBH model, the nonoptimized regions of the system are represented by one or several blocks, which can move as rigid bodies with respect to the atoms of the optimized region. In this way unphysical imaginary frequencies are avoided and the translational/rotational invariance of the potential energy surface is fully respected. In this paper we focus on issues concerning the practical numerical implementation of the MBH model. The MBH normal mode equations are worked out for several coordinate choices. The introduction of a consistent group-theoretical notation facilitates the treatment of both the case of a single block and the case of multiple blocks. Special attention is paid to the formulation in terms of Cartesian variables, in order to provide a link with the standard output of common molecular modeling programs.