The calculation of force constants and normal coordinates IV: XH4 and XH3 molecules

Normal coordinate calculations of XH4 and XH3 molecules are reviewed and discussed. It is shown that for most of these molecules the true values of the force constants in the most General Harmonic Force Field can be uniquely determined only by making use of vibration-rotation interaction constants. It is emphasized that without these extra data the GFF is not determined. The results are compared with various model force fields for these molecules.

The calculation of accurate normal coordinates I: general theory, and application to methane

The calculation of accurate and reliable vibrational potential functions and normal co-ordinates is discussed, for such simple polyatomic molecules as it may be possible. Such calculations should be corrected for the effects of anharmonicity and of resonance interactions between the vibrational states, and should be fitted to all the available information on all isotopic species: particularly the vibrational frequencies, Coriolis zeta constants and centrifugal distortion constants. The difficulties of making these corrections, and of making use of the observed data are reviewed. A programme for the Ferranti Mercury Computer is described by means of which harmonic vibration frequencies and normal co-ordinate vectors, zeta factors and centrifugal distortion constants can be calculated, from a given force field and from given G-matrix elements, etc. The programme has been used on up to 5 × 5 secular equations for which a single calculation and output of results takes approximately l min; it can readily be extended to larger determinants. The best methods of using such a programme and the possibility of reversing the direction of calculation are discussed. The methods are applied to calculating the best possible vibrational potential function for the methane molecule, making use of all the observed data.

Normal coordinates for the planar vibrations of benzene

Force constants and normal coordinates have been recalculated for all of the in-plane vibrations of benzene, making use of the recently observed data on one of the Coriolis constants in the E2g species from the work of Callomon et al. The extent to which the force field is uniquely determined by the data is reviewed for each symmetry species in turn, and the results of a force constant refinement calculation are reported in which a modified valency force field was used based on the hybrid orbital model. The results show considerable differences from Whiffen's normal coordinates for benzene, but somewhat smaller differences from Scherer's recent calculations.

Normal coordinates in the methyl halides

The potential‐energy functions found by Chang for the methyl halides have been put into valence‐type form and revised to eliminate inconsistencies and to accord with the true (nontetrahedral) geometry and the normal frequencies (corrected for Fermi resonance and anharmonicity). The resulting valence‐type force constants and normal coordinates are given for light (CH3) and heavy (CD3) chlorides, bromides, and iodides.

The calculation of force constants and normal coordinates III: constrained force fields for the methyl halides

A method is discussed for imposing any desired constraint on the force field obtained in a force constant refinement calculation. The application of this method to force constant refinement calculations for the methyl halide molecules is reported. All available data on the vibration frequencies, Coriolis interaction constants and centrifugal stretching constants of CH3X and CD3X molecules were used in the refinements, but despite this apparent abundance of data it was found that constraints were necessary in order to obtain a unique solution to the force field. The results of unconstrained calculations, and of three different constrained calculations, are reported in this paper. The constrained models reported are a Urey—Bradley force field, a modified valence force field, and a constraint based on orbital-following bond-hybridization arguments developed in the following paper. The results are discussed, and compared with previous results for these molecules. The third of the above models is found to reproduce the observed data better than either of the first two, and additional reasons are given for preferring this solution to the force field for the methyl halide molecules.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Investigation of short-time isomerization dynamics in p-nitroazobenzene from resonance Raman intensity analysis

Resonance Raman (RR) spectra are presented for p-nitroazobenzene dissolved in chloroform using 18 excitation Wavelengths, covering the region of (1)(n --> pi*) electronic transition. Raman intensities are observed for various totally symmetric fundamentals, namely, C-C, C-N, N=N, and N-O stretching vibrations, indicating that upon photoexcitation the excited-state evolution occurs along all of these vibrational coordinates. For a few fundamentals, interestingly, in p-nitroazobenzene, it is observed that the RR intensities decrease near the maxima of the resonant electronic (1)(n --> pi*) transition. This is attributed to the interference from preresonant scattering due to the strongly allowed (1)(pi --> pi*) electronic transition. The electronic absorption spectrum and the absolute Raman cross section for the nine Franck-Condon active fundamentals of p-nitroazobenzene have been successfully modeled using Heller's time-dependent formalism for Raman scattering. This employs harmonic description of the lowest energy (1)(n --> pi*) potential energy surface. The short-time isomerization dynamics is then examined from a priori knowledge of the ground-state normal mode descriptions of p-nitroazobenzene to convert the wave packet motion in dimensionless normal coordinates to internal coordinates. It is observed that within 20 fs after photoexcitation in p-nitroazobenzene, the N=N and C-N stretching vibrations undergo significant changes and the unsubstituted phenyl ring and the nitro stretching vibrations are also distorted considerably.

Investigation of molecular polarizabilities and derivatives in halomethanes

The research undertaken was to obtain absolute Raman intensities for the symmetric stretching vibrations of the methyl halides, CH3X with (X=F, CI, Br), by experiment and theory. The intensities were experimentally measured using the Ar+ ion gas laser as excitation source, a Spex 14018 double monochromator and a RCA C-31034 photomultiplier tube as detector. These intensities arise from changes in the derivative of the polarizability (8 a'), with respect to vibration along a normal coordinate (8qi). It was intended that these derivatives obtained with respect to normal coordinates would be converted to derivatives with respect to internal coordinates, for a quantitative comparison with theory. Theoretical numerical polarizability derivatives for the stretching vibrations are obtained using the following procedure. A vibration was simulated in the molecule by increasi.ng and decreasing the respective bond by the amount ±o.oosA for the C-H bonds and ±o.oIA for the C-X (X=F, CI, Br) bond. The derivative was obtained by taking the difference in the polarizability for the equilibrium geometry and the geometry when a particular bond is changed. This difference, when divided by the amount of change in each bond and the number of bonds present results in the derivative of the polarizability with respect to internal coordinate i.e., !1u/!1r. These derivatives were obtained by two methods: I} ab initio molecular orbital calculation and 2} theory of atoms in molecules (AIM) analysis. Due to errors in the experimental setup only a qualitative analysis of the results was undertaken relative to the theory. Theoretically it is predicted that the symmetric carbonhalogen stretch vibrations are more intense than the respective carbon-hydrogen stretch, but only for the methyl chloride and bromide. The carbon fluorine stretch is less intense than the carbon-hydrogen stretch, a fact which is attributed to the small size and high electronegativity of the fluorine atom. The experimental observations are seen to agree qualitatively with the theory results. It is hoped that when the experiment is repeated, a quantitative comparison can be made. The analysis by the theory of atoms in molecules, along with providing polarizabilities and polarizability derivatives, gives additional information outlined below. The theory provides a pictorial description of the main factors contributing to the molecular polarizability and polarizability derivative. These contributions are from the charge transfer and atomic dipole terms i.e., transfer of charge from one atom to another and the reorganization of atomic electronic charge distribution due to presence of an electric field. The linear relationship between polarizability and molecular volume was also observed.

The electronic spectra and the determination of the excited state geometry of 3A" propynal /|nby Chhiu-Tsu Lin. -- 260 St. Catharines, Ont. : [s. n.],

Impurity free eluission spectra of HCCCHO and DCCCHO have been rephotographed using the electronic-energy-exchange method with benzene as a carrier gas. The near ultraviolet spectra of ReeCHO and DCCCHO were photographed in a sorption under conditions of high resolution with absorption path lengths up to 100 meters. The emission and absorption spectra of Propynal resulting from 3 n 1 t 1\ - A excitation has been reanalyzed in som.e detail. Botrl of the eH out-of-plane wagging modes were found to have negative anharmonicity. A barrier height of 56.8/0.0 cm- 1 and a nonplanar oft , , equilibrium angle of 17 3 /30 are calculated for the V 10/ lJ 11 modes. The in-plane and out-of-plane v1. brational modes in the 3A." and 1a~. ' elec ronic states of Propynal were subjected to a normal coordinate treatment in the approximat :on of tIle Urey-Bradley force field. From the relative oscillator strengths of the trans1·t1·0ns connect i ng t he v ibrat1•0n1ess lA' , state and t,he V1· bron1·C 3· if levels of the A state, the differences in equilibrium configuration were evaluated from an approximate Franck-Condon analysis based on the ground state normal coordinates. As this treatment gave 512 possible geometrical structures for the upper state, it 4 was necessary to resort to a comparison of the observed and calculated moments of inertia along with chemical intuition to isolate the structure. A test of the correctness of the calculated structure change and the vibrational assignment was raade by evaluating the intensities of the inplane and out-oi-plane fundarnental, sequence, and cross sequellce transitions y the exact Franck-Condon method.

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

Relations among several nuclear and electronic density functional reactivity indexes

A set of connections among several nuclear and electronic indexes of reactivity in the framework of the conceptual Density Functional Theory by using an expansion ofthe energy functional in terms of the total number of electrons and the normal coordinates within a canonical ensemble was derived. The relations obtained provided explicit links between important quantities related to the chemical reactivity of a system. This paper particularly demonstrates that the derivative of the electronic energy with respect to the external potential of a system in its equilibrium geometry was equal to the negative of the nuclear repulsion derivative with respect to the external potential

Desenvolupament, implementació i aplicació de noves metodologies per a l'estudi de la reactivitat química i de les propietats òptiques no lineals

El coneixement de la superfície d'energia potencial (PES) ha estat essencial en el món de la química teòrica per tal de discutir tant la reactivitat química com l'estructura i l'espectroscòpia molecular. En el camp de la reactivitat química es hem proposat continuar amb el desenvolupament de nova metodologia dins el marc de la teoria del funcional de la densitat conceptual. En particular aquesta tesis es centrarà en els següents punts: a) El nombre i la naturalesa dels seus punts estacionaris del PES poden sofrir canvis radicals modificant el nivell de càlcul utilitzats, de tal manera que per estar segurs de la seva naturalesa cal anar a nivells de càlcul molt elevats. La duresa és una mesura de la resistència d'un sistema químic a canviar la seva configuració electrònica, i segons el principi de màxima duresa on hi hagi un mínim o un màxim d'energia trobarem un màxim o un mínim de duresa, respectivament. A l'escollir tot un conjunt de reaccions problemàtiques des del punt de vista de presència de punts estacionaris erronis, hem observat que els perfils de duresa són més independents de la base i del mètode utilitzats, a més a més sempre presenten el perfil correcte. b) Hem desenvolupat noves expressions basades en les integracions dels kernels de duresa per tal de determinar la duresa global d'una molècula de manera més precisa que la utilitzada habitualment que està basada en el càlcul numèric de la derivada segona de l'energia respecte al número d'electrons. c) Hem estudiat la validesa del principis de màxima duresa i de mínima polaritzabiliat en les vibracions asimètriques en sistemes aromàtics. Hem trobat que per aquests sistemes alguns modes vibracionals incompleixen aquests principis i hem analitzat la relació d'aquest l'incompliment amb l'efecte de l'acoblament pseudo-Jahn-Teller. A més a més, hem postulat tot un conjunt de regles molt senzilles que ens permetien deduir si una molècula compliria o no aquests principis sense la realització de cap càlcul previ. Tota aquesta informació ha estat essencial per poder determinar exactament quines són les causes del compliment o l'incompliment del MHP i MPP. d) Finalment, hem realitzat una expansió de l'energia funcional en termes del nombre d'electrons i de les coordenades normals dintre del conjunt canònic. En la comparació d'aquesta expansió amb l'expansió de l'energia del nombre d'electrons i del potencial extern hem pogut recuperar d'una altra forma diferent tot un conjunt de relacions ja conegudes entre alguns coneguts descriptors de reactivitat del funcional de la densitat i en poden establir tot un conjunt de noves relacions i de nous descriptors. Dins del marc de les propietats moleculars es proposa generalitzar i millorar la metodologia pel càlcul de la contribució vibracional (Pvib) a les propietats òptiques no lineals (NLO). Tot i que la Pvib no s'ha tingut en compte en la majoria dels estudis teòrics publicats de les propietats NLO, recentment s'ha comprovat que la Pvib de diversos polímers orgànics amb altes propietats òptiques no lineals és fins i tot més gran que la contribució electrònica. Per tant, tenir en compte la Pvib és essencial en el disseny dels nous materials òptics no lineals utilitzats en el camp de la informàtica, les telecomunicacions i la tecnologia làser. Les principals línies d'aquesta tesis sobre aquest tema són: a) Hem calculat per primera vegada els termes d'alt ordre de Pvib de diversos polímers orgànics amb l'objectiu d'avaluar la seva importància i la convergència de les sèries de Taylor que defineixen aquestes contribucions vibracionals. b) Hem avaluat les contribucions electròniques i vibracionals per una sèrie de molècules orgàniques representatives utilitzant diferents metodologies, per tal de poder de determinar quina és la manera més senzilla per poder calcular les propietats NLO amb una precisió semiquantitativa.

Vibrational intensities VI: ethylene and its deutero-isotopes

Absolute intensity measurements have been made on the fundamental vibrations of ethylene and four of its deuteroisotopes. The bands were pressure broadened with nitrogen at 50 atmos, and the intensities were determined by the method of Wilson and Wells except that the observed optical density was integrated against logv rather than v. Normal coordinates have been calculated, and the intensities have been interpreted in terms of quantities (∂p/∂Si) giving the change in dipole moment with respect to each internal symmetry coordinate. Data from the different isotopic species have been used to eliminate ambiguities in the interpretation. Effective bond moments are calculated for each symmetry coordinate.

Vibrational intensities VII: ethane and ethane-d6

Absolute intensity measurements have been made on the fundamental vibrations of C2H6 and C2D6, using the extrapolation method of Wilson and Wells and using nitrogen at pressures up to 50 atmospheres to broaden the bands. The absorption coefficient was integrated against the logarithm of the frequency. Normal coordinates were calculated from the potential function of Hansen and Dennison, and were used to interpret the results in terms of quantities (∂p/∂Si) giving the change of dipole moment with respect to the symmetry coordinates Si. Consistency of data between the isotopes was used both to eliminate ambiguities in the interpretation, and as a criterion in separating overlapping pairs of absorption bands. The results have been interpreted in terms of bond effective moments.