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.

The calculation of force constants and normal co-ordinates II: methyl fluoride

The mathematical difficulties which can arise in the force constant refinement procedure for calculating force constants and normal co-ordinates are described and discussed. The method has been applied to the methyl fluoride molecule, using an electronic computer. The best values of the twelve force constants in the most general harmonic potential field were obtained to fit twenty-two independently observed experimental data, these being the six vibration frequencies, three Coriolis zeta constants and two centrifugal stretching constants DJ and DJK, for both CH3F and CD3F. The calculations have been repeated both with and without anharmonicity corrections to the vibration frequencies. All the experimental data were weighted according to the reliability of the observations, and the corresponding standard errors and correlation coefficients of the force constants have been deduced. The final force constants are discussed briefly, and compared with previous treatments, particularly with a recent Urey-Bradley treatment for this molecule.

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.

Diagonalization and Jordan Normal Form – motivation through Maple

Following an introduction to the diagonalization of matrices, one of the more difficult topics for students to grasp in linear algebra is the concept of Jordan normal form. In this note, we show how the important notions of diagonalization and Jordan normal form can be introduced and developed through the use of the computer algebra package Maple®.

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 co-ordinates--VII Normal co-ordinates in the E species of allene

Force constant and normal co-ordinate calculations are reported for the E species vibrations of the allene molecule. Data on the fundamental vibration frequencies of allene-h4, allene-d4 and allene-1.1-d2 and on the five experimentally determined Coriolis zeta constants of C3H4 and C3D4, were used in a force constant refinement procedure. Allowing for product and sum rules this gives 21 independent data which were used to refine to the most general harmonic force field (10 parameters) with one constraint (in the absence of any constraints the refinement was not satisfactory). The results have been used to calculate the complete ζz Coriolis interaction matrix for the allene-1.1-d2 molecule, and hence to calculate the expected rotational structure of the perpendicular bending vibrations of this molecule; the good agreement obtained with the observed spectra is a check on our results.

Local mode and normal mode interpretations of the CH and CD stretching vibrational manifolds of 1,1-dideuteroethene

The =CH2 AND =CD2 stretching vibrational overtones of H2C=CD2 have been studied up to V= 6 and V= 3, respectively. We report their interpretation in terms of a transition from normal to local modes, involving Fermi resonance with the C=C stretching and CH2 scissoring vibrations. We discuss the alternative representation of the vibrational Hamiltonian matrix in local mode and normal mode basis functions, and conclude that the normal mode basis offers greater flexibility in representing small anharmonic couplings with other modes.

Effect of type and timing of oil supplements to sows during pregnancy on the growth performance and endocrine profile of low and normal birth weight offspring

Eighty-eight multiparous sows were used to evaluate whether type and timing of oil supplementation during gestation influences the incidence of low birth weight (LBW). Sows were allocated (eight per treatment) commercial sow pellets (3 kg/d; control diet) or an experimental diet consisting of control diet plus 10 % extra energy in the form of excess pellets, palm oil, olive oil (OO), sunflower oil (SO) or fish oil; experimental diets were fed during either the first half (G1) or second half (G2) of gestation. Growth performance and endocrine profile of LBW ( < 1·09 kg) and normal birth weight (NBW; 1·46–1·64 kg) offspring were compared. Maternal dietary supplementation altered the distribution curve for piglet birth weight. SOG1 sows had a greater proportion of LBW piglets (P < 0·05), whilst it was reduced in the OOG1 group (P < 0·05). Growth rate of LBW piglets was lower compared with their NBW siblings (P < 0·05) when dietary supplementation was offered in G2 but were similar for G1. At birth, LBW offspring of supplemented animals possessed more fat compared with the control group (P < 0·05); LBW offspring of control animals exhibited a more rapid decline in fat free mass/kg prior to weaning. Plasma metabolites and insulin concentrations were influenced by maternal diet and birth weight. In conclusion, maternal dietary supplementation altered the distribution of piglet birth weights and improved the energy status of LBW piglets. Supplementation with MUFA during G1 reduced the incidence of LBW, whereas PUFA had the reverse effect.

Rates of production of acetate, propionate, and butyrate in the rumen of lactating dairy cows given normal and low-roughage diets

Five lactating dairy cows with a permanent cannula in the rumen were given ( kg DM/d) a normal diet (7.8 concentrates, 5.1 hay) or a low-roughage (LR) diet (11.5 concentrates, 1.2 hay) in two meals daily in a two-period crossover design. Milk fat (g/kg) was severely reduced on diet LR. To measure rates of production of individual volatile fatty acids (VFA) in the rumen, 0.5 mCi 1-C-14-acetic acid, 2-C-14-propionic acid, or 1-C-14-n-butyric acid were infused into the rumen for 22 h at intervals of 2 to 6 d; rumen samples were taken over the last 12 h. To measure rumen volume, we infused Cr-EDTA into the rumen continuously, and polyethylene glycol was injected 2 h before the morning feed. Results were very variable, so volumes measured by rumen emptying were used instead. Net production of propionic acid more than doubled on LR, but acetate and butyrate production was only numerically lower. Net production rates pooled across both diets were significantly related to concentrations for each VFA. Molar proportions of net production were only slightly higher than molar proportions of concentrations for acetate and propionate but were lower for butyrate. The net energy value (MJ/d) of production of the three VFA increased from 89.5 on normal to 109.1 on LR, equivalent to 55 and 64% of digestible energy, respectively. Fully interchanging, three-pool models of VFA C fluxes are presented. It is concluded that net production rates of VFA can be measured in non-steady states without the need to measure rumen volumes.

The F test for comparing two normal variances: correct and incorrect calculation of the two-sided p-value?

This article illustrates that not all statistical software packages are correctly calculating a p-value for the classical F test comparison of two independent Normal variances. This is illustrated with a simple example, and the reasons why are discussed. Eight different software packages are considered.

Estimation following selection of the largest of two normal means

This paper considers the problem of estimation when one of a number of populations, assumed normal with known common variance, is selected on the basis of it having the largest observed mean. Conditional on selection of the population, the observed mean is a biased estimate of the true mean. This problem arises in the analysis of clinical trials in which selection is made between a number of experimental treatments that are compared with each other either with or without an additional control treatment. Attempts to obtain approximately unbiased estimates in this setting have been proposed by Shen [2001. An improved method of evaluating drug effect in a multiple dose clinical trial. Statist. Medicine 20, 1913–1929] and Stallard and Todd [2005. Point estimates and confidence regions for sequential trials involving selection. J. Statist. Plann. Inference 135, 402–419]. This paper explores the problem in the simple setting in which two experimental treatments are compared in a single analysis. It is shown that in this case the estimate of Stallard and Todd is the maximum-likelihood estimate (m.l.e.), and this is compared with the estimate proposed by Shen. In particular, it is shown that the m.l.e. has infinite expectation whatever the true value of the mean being estimated. We show that there is no conditionally unbiased estimator, and propose a new family of approximately conditionally unbiased estimators, comparing these with the estimators suggested by Shen.

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.