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.

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.

Coriolis perturbations in the infrared spectrum of ethylene

Rotational structure has been resolved and analyzed in two of the infrared‐active perpendicular bands of C2H4 vapor: the Type b fundamental band, ν10, at 826 cm—1, and the Type c fundamental band, ν7, at 949 cm—1. Many of the individual PP and RR branch lines have been observed. The analysis has been confined to values of the quantum number K≥3, for which energy levels ethylene shows no detectable deviations from a symmetric‐top rotational structure. The analysis reveals a Coriolis interaction between ν7 and ν10, and between ν4 and ν10, and values of the Coriolis constants ζ7,10z and ζ4,10y are obtained; these are related to normal coordinate calculations for the appropriate symmetry species, and force constants are derived to fit the observed zeta constants. The band center of ν10 has been revised from the original figure of 810 cm—1 to the new value, 826 cm—1, and the inactive frequency ν4 is estimated to lie at 1023±3 cm—1, in good agreement with the previous estimate of 1027 cm—1. The change in the value of ν10 leads to a suggested change in the value of the Raman‐active fundamental ν6 from 1236 to 1222 cm—1. New combination bands have been observed at 2174 cm—1, assigned as ν3+ν10; and at 2252 cm—1, assigned as ν4+ν6; also rotational structure has been resolved and analyzed in the ν6+ν10 band at 2048 cm—1. The new data obtained for the C2H4 molecule are summarized in Table XII, with all of the other data presently available on the vibrational and rotational constants.

The puckering coordinate in oxetane

The recently determined apucker(xz) = ∂Ixz/∂Qpucker constants for five isotopic species of oxetane have been interpreted quantitatively by allowing the ring-puckering normal coordinate to contain a small amount of methylene rocking and twisting motion. The observed data are most closely reproduced by assuming a methylene rock angle of about 3° and a twist angle for the -methylene groups of about 1° at the equilibrium conformation of the ring.

Advanced multivariate analysis to assess remediation of hydrocarbons in soils

Accurate monitoring of degradation levels in soils is essential in order to understand and achieve complete degradation of petroleum hydrocarbons in contaminated soils. We aimed to develop the use of multivariate methods for the monitoring of biodegradation of diesel in soils and to determine if diesel contaminated soils could be remediated to a chemical composition similar to that of an uncontaminated soil. An incubation experiment was set up with three contrasting soil types. Each soil was exposed to diesel at varying stages of degradation and then analysed for key hydrocarbons throughout 161 days of incubation. Hydrocarbon distributions were analysed by Principal Coordinate Analysis and similar samples grouped by cluster analysis. Variation and differences between samples were determined using permutational multivariate analysis of variance. It was found that all soils followed trajectories approaching the chemical composition of the unpolluted soil. Some contaminated soils were no longer significantly different to that of uncontaminated soil after 161 days of incubation. The use of cluster analysis allows the assignment of a percentage chemical similarity of a diesel contaminated soil to an uncontaminated soil sample. This will aid in the monitoring of hydrocarbon contaminated sites and the establishment of potential endpoints for successful remediation.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

General derivative relations for anharmonic force fields

The brace notation, introduced by Allen and Csaszar (1993, J. chem. Phys., 98, 2983), provides a simple and compact way to deal with derivatives of arbitrary non-tensorial quantities. One of its main advantages is that it builds the permutational symmetry of the derivatives directly into the formalism. The brace notation is applied to formulate the general nth-order Cartesian derivatives of internal coordinates, and to provide closed forms for general, nth-order transformation equations of anharmonic force fields, expressed as Taylor series, from internal to Cartesian or normal coordinate spaces.

High resolution photoacoustic overtone spectroscopy of monofluoroacetylene, HCCF, with a titanium:sapphire ring laser

Intracavity photoacoustic overtone spectrum of monofluoroacetylene, HCCF, has been recorded in the wave number region 10 750–14 500 cm−1 with a titanium:sapphire ring laser. The spectrum contains many dense vibration–rotation band systems which can be resolved with Doppler limited resolution. Altogether 58 individual overtone bands have been analyzed rotationally. Many of the observed bands show perturbations of which some have been attributed to anharmonic resonance interactions. A Fermi resonance model based on conventional rectilinear normal coordinate theory has been used to assign vibrationally bands from this work and from earlier studies. Many of the observed vibrational term values and rotational constants can be reproduced well with this model. The results show the importance of the Fermi resonance interactions at the high overtone excitations.

Quartic anharmonic resonances in acetylene molecules

Formulas are derived for the quartic anharmonic resonance coefficients observed to be important between C–H stretching and the combination of one quantum of C≡C stretching and two quanta of H–C≡C bending in a number of acetylene molecules. Examples of this resonance are ν3 with ν2+ν4+ν5 in 12C2H2, ν1 with ν2+2ν5 in 13C2H2, and ν1 with ν2+2ν4 in monofluoroacetylene and monochloroacetylene. The coefficients characterizing the resonances in these examples, which we denote K3,245, K1,255, and K1,244, arise from cubic and quartic terms in the anharmonic force field, in the normal coordinate representation, through second order and first order perturbation treatments respectively, where the second order resonances are calculated by a Van Vleck resonance formalism. The experimentally determined values of these coefficients are compared with values calculated from model anharmonic force fields.

Infrared intensities of furan, pyrrole and thiophene: beyond the double harmonic approximation

Infrared intensities of the fundamental, overtone and combination transitions in furan, pyrrole and thiophene have been calculated using the variational normal coordinate code MULTIMODE. We use pure vibrational wavefunctions, and quartic force fields and cubic dipole moment vector surfaces, generated by density functional theory. The results are compared graphically with second-order perturbation calculations and with relative intensities from experiment for furan and pyrrole.

The vibrations and tunnelling of malonaldehyde on a Moller-Plesset surface

The vibrations and tunnelling motion of malonaldehyde have been studied in their full dimensionality using an internal coordinate path Hamiltonian. In this representation there is one large amplitude internal coordinate s and 3N - 7 (=20) normal coordinates Q which are orthogonal to the large amplitude motion at all points. It is crucial that a high accuracy potential energy surface is used in order to obtain a good representation for the tunneling motion; we use a Moller-Plesset (MP2) surface. Our methodology is variational, that is we diagonalize a sufficiently large matrix in order to obtain the required vibrational levels, so an exact representation for the kinetic energy operator is used. In a harmonic valley representation (s, Q) complete convergence of the normal coordinate motions and the internal coordinate motions has been obtained; for the anharmonic valley in which we use two- and three-body terms in the surface (s, Q(1), Q(2)), we also obtain complete convergence. Our final computed stretching fundamentals are deficient because our potential energy surface is truncated at quartic terms in the normal coordinates, but our lower fundamentals are good.

Baroclinic instability on hot extrasolar planets

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value cal- culations are carried out – focusing on the freely-evolving, adiabatic situation. Using a high- resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot ex- trasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qual- itative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks – but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.

Combined affinity labelling and mass spectrometry analysis of differential cell surface protein expression in normal and prostate cancer cells

Differences in the expression of cell surface proteins between a normal prostate epithelial (1542-NP2TX) and a prostate cancer cell line (1542-CP3TX) derived from the same patient were investigated. A combination of affinity chromatographic purification of biotin-tagged surface proteins with mass spectrometry analysis identified 26 integral membrane proteins and 14 peripheral surface proteins. The findings confirm earlier reports of altered expression in prostate cancer for several cell surface proteins, including ALCAM/CD166, the Ephrin type A receptor, EGFR and the prostaglandin F2 receptor regulatory protein. In addition, several novel findings of differential expression were made, including the voltage-dependent anion selective channel proteins Porin 1 and 2, ecto-5'-nucleotidase (CD73) and Scavenger receptor B1. Cell surface protein expression changed both qualitatively and quantitatively when the cells were grown in the presence of either or both interferon INFalpha and INFgamma. Costimulation with type I and II interferons had additive or synergistic effects on the membrane density of several, mainly peripherally attached surface proteins. Concerted upregulation of surface exposed antigens may be of benefit in immuno-adjuvant-based treatment of interferon-responsive prostate cancer. In conclusion, this study demonstrates that differences in the expression of membrane proteins between normal and prostate cancer cells are reproducibly detectable following vectorial labelling with biotin, and that detailed analysis of extracellular-induced surface changes can be achieved by combining surface-specific labelling with high-resolution two-dimensional gel electrophoresis and mass spectrometry.

Combined affinity labelling and mass spectrometry analysis of differential cell surface protein expression in normal and prostate cancer cells

Differences in the expression of cell surface proteins between a normal prostate epithelial (1542-NP2TX) and a prostate cancer cell line (1542-CP3TX) derived from the same patient were investigated. A combination of affinity chromatographic purification of biotin-tagged surface proteins with mass spectrometry analysis identified 26 integral membrane proteins and 14 peripheral surface proteins. The findings confirm earlier reports of altered expression in prostate cancer for several cell surface proteins, including ALCAM/CD166, the Ephrin type A receptor, EGFR and the prostaglandin F2 receptor regulatory protein. In addition, several novel findings of differential expression were made, including the voltage-dependent anion selective channel proteins Porin 1 and 2, ecto-5'-nucleotidase (CD73) and Scavenger receptor B1. Cell surface protein expression changed both qualitatively and quantitatively when the cells were grown in the presence of either or both interferon INF alpha and INF gamma. Costimulation with type I and II interferons had additive or synergistic effects on the membrane density of several, mainly peripherally attached surface proteins. Concerted upregulation of surface exposed antigens may be of benefit in immuno-adjuvant-based treatment of interferon-responsive prostate cancer. In conclusion, this study demonstrates that differences in the expression of membrane proteins between normal and prostate cancer cells are reproducibly detectable following vectorial labelling with biotin, and that detailed analysis of extracellular-induced surface changes can be achieved by combining surface-specific labelling with high-resolution two-dimensional gel electrophoresis and mass spectrometry.