Feature tracking on the unit-sphere

Techniques used in a previous study of the objective identification and tracking of meteorological features in model data are extended to the unit sphere. An alternative feature detection scheme is described based on cubic interpolation for the sphere and local maximization. The extension of the tracking technique, used in the previous study, to the unit sphere is described. An example of the application of these techniques to a global relative vorticity field from a model integration are presented and discussed.

A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

Wavelets and the generalization of the variogram

The experimental variogram computed in the usual way by the method of moments and the Haar wavelet transform are similar in that they filter data and yield informative summaries that may be interpreted. The variogram filters out constant values; wavelets can filter variation at several spatial scales and thereby provide a richer repertoire for analysis and demand no assumptions other than that of finite variance. This paper compares the two functions, identifying that part of the Haar wavelet transform that gives it its advantages. It goes on to show that the generalized variogram of order k=1, 2, and 3 filters linear, quadratic, and cubic polynomials from the data, respectively, which correspond with more complex wavelets in Daubechies's family. The additional filter coefficients of the latter can reveal features of the data that are not evident in its usual form. Three examples in which data recorded at regular intervals on transects are analyzed illustrate the extended form of the variogram. The apparent periodicity of gilgais in Australia seems to be accentuated as filter coefficients are added, but otherwise the analysis provides no new insight. Analysis of hyerpsectral data with a strong linear trend showed that the wavelet-based variograms filtered it out. Adding filter coefficients in the analysis of the topsoil across the Jurassic scarplands of England changed the upper bound of the variogram; it then resembled the within-class variogram computed by the method of moments. To elucidate these results, we simulated several series of data to represent a random process with values fluctuating about a mean, data with long-range linear trend, data with local trend, and data with stepped transitions. The results suggest that the wavelet variogram can filter out the effects of long-range trend, but not local trend, and of transitions from one class to another, as across boundaries.

An X-ray micro-tomography system optimised for the low-dose study of living organisms

An X-ray micro-tomography system has been designed that is dedicated to the low-dose imaging of radiation sensitive living organisms and has been used to image the early development of the first few days of plant development immediately after germination. The system is based on third-generation X-ray micro-tomography system and consists of an X-ray tube, two-dimensional X-ray detector and a mechanical sample manipulation stage. The X-ray source is a 50 kVp X-ray tube with a silver target with a filter to centre the X-ray spectrum on 22 keV.A 100 mm diameter X-ray image intensifier (XRII) is used to collect the two-dimensional projection images. The rotation tomography table incorporates a linear translation mechanism to eliminate ring artefact that is commonly associated with third-generation tomography systems' Developing maize seeds (Triticum aestivum) have been imaged using the system with a cubic voxel linear dimension of 100 mum, over a diameter of 25 mm and the root lengths and volumes measured. The X-ray dose to the plants was also assessed and found to have no effect on the plant root development. (C) 2003 Elsevier Science Ltd. All rights reserved.

Anharmonic force field of acetylene

The harmonic and anharmonic force field of acetylene has been determined in a least-squares calculation from recently determined data on the spectroscopic constants of various isotopic species (including the vibrational l-doubling constant). A general quadratic and cubic force field was used, but a constrained quartic force field containing only 8 of the 23 possible quartic constants. The results are discussed and compared with earlier work.

A refined quartic forcefield for acetylene: accurate calculation of the vibrational spectrum

It is now possible to calculate the nine-dimensional rovibrational wavefunctions of sequentially bonded four-atom molecules variationally without dynamical approximation. In the case of HCCH, the simplest such molecule, many hundreds of rovibrational (J = 0, 1, 2) levels can be converged to better than 1.5 cm −1. Variational calculations of this kind are used here systematically to refine the well-known quartic valence-coordinate forcefleld of Strey and Mills [J.Mol. Spectrosc.59, 103-115 (1976)] against experimental term values up to three C-H stretch quanta for the principal and two deuterated isotopomers, yielding a new surface that reproduces the energies of all the known Σ, Π, and Δ states of these species up to the energy of two C-H stretch quanta with an rms error of 3 cm−1 . The refined forcefield is used to study the resonances associated with the accidental degeneracies (ν2 + ν4 + ν5, ν3) and (ν2 + 2ν5, ν1) in the principal isotopomer, leading to a clarification of the assignment of she experimentally detected states in the 2ν3 and 3ν3, polyads, and to the finding that vibrational Coriolis (kinetic energy) terms, rather than quartic anharmonicities in the potential, are the primary cause of the resonant interactions. Using a new cubic ab initio electric dipole field to calculate IR absorption coefficients, 24 undetected Σ and Π states of 1H12C12C1H and 5 undetected Σ states of D12C12CD are identified as candidates for experimental study, and their calculated energies and assignments are given.

Excited vibrational state microwave spectra, structure, and force-field of trifluorosilane

The J = 2−1 microwave spectrum of six isotopic species of HSiF3 has been observed and assigned in excited states of five of the six fundamental vibrations. The assignment is based on relative intensities, double resonance experiments, and trial anharmonic force constant calculations. Analysis of the spectra leads to experimental values for five of the constants, all three l-doubling constants qt, one Fermi resonance constant φ233, and one zeta constant. The harmonic force field has been refined to all the available data on vibration wavenumbers, centrifugal distortion constants, and zeta constants. The cubic anharmonic force field has been refined to the data on and qt constants, using two models: a valence force model with two cubic force constants for SiH and SiF stretching, and a more sophisticated model. With the help of these calculations, the following equilibrium structure has been determined: re(SiH) = 1.4468(±5) Å, re(SiF) = 1.5624(±1) Å, HSiF = 110.64(±3)°,

The anharmonic force field of HCN and HCP

The quadratic, cubic, and quartic force field of HCN has been calculated by a least squares refinement to fit the most recent observed data on the vibration-rotation constants of HCN, DCN and H13CN. All of the observed parameters are fitted within their standard errors of observation. The corresponding parameters for other isotopic species are calculated. For HCP and DCP the more limited data available have been fitted to an anharmonic force field using constraints based on comparison with HCN. Using this force field the zero-point rotational constants B0 have been corrected to obtain the equilibrium constants Be, and hence the equilibrium structure has been determined to be re(CH) = 1•0692(7)A, and re(CP) = 1•5398(2)A.

An anharmonic force field for SO3

An anharmonic force field for SO3 based on the valence force model has been investigated. The results of extending the model to include some further estimated cubic interaction potential constants have also been investigated. The phenomenological parameters calculated from both model force fields agree with those few values which have been experimentally determined. A calculation of the inertia defect has been made, and thus the value of C0 has been determined. The equilibrium structure has been determined to be: re = 1.4184 ± 0.0010 Å.

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.

Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.

Radon in homes and risk of lung cancer: collaborative analysis of individual data from 13 European case-control studies

Objective: To determine the risk of lung cancer associated with exposure at home to the radioactive disintegration products of naturally Occurring radon gas. Design: Collaborative analysis of individual data from 13 case-control studies of residential radon and lung cancer. Setting Nine European countries. Subjects 7148 cases Of lung cancer and 14 208 controls. Main outcome measures: Relative risks of lung cancer and radon gas concentrations in homes inhabited during the previous 5-34 years measured in becquerels (radon disintegrations per second) per cubic incite (Bq/m(3)) Of household air. Results: The mean measured radon concentration in homes of people in tire control group was 97 Bq/m(3), with 11% measuring > 200 and 4% measuring > 400 Bq/m(3). For cases of lung cancer the mean concentration was 104 Bq/m(3). The risk of lung cancer increased by 8.4% (95% confidence interval 3.0% to 15.8%) per 100 Bq/m(3) increase in measured radon (P = 0.0007). This corresponds to an increase of 16% (5% to 31%) per 100 Bq/m(3) increase in usual radon-that is, after correction for the dilution caused by random uncertainties in measuring radon concentrations. The dose-response relation seemed to be linear with no threshold and remained significant (P=0.04) in analyses limited to individuals from homes with measured radon < 200 Bq/m(3). The proportionate excess risk did not differ significantly with study, age, sex, or smoking. In the absence of other causes of death, the absolute risks of lung cancer by age 75 years at usual radon concentrations of 0, 100, and 400 Bq/m(3) would be about 0.4%, 0.5%, and 0.7%, respectively, for lifelong non-smokers, and about 25 times greater (10%, 12%, and 16%) for cigarette smokers. Conclusions: Collectively, though not separately, these studies show appreciable hazards from residential radon, particularly for smokers and recent ex-smokers, and indicate that it is responsible for about 2% of all deaths from cancer in Europe.

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.

Carbon-decorated FePt nanoparticles

FePt magnetic nanoparticles are an important candidate material for many future magnetic applications. FePt exists as two main phases, that is, a disordered face-centered cubic (fcc) structure, which is generally prepared by chemical methods at low temperatures, and the high-temperature chemically ordered face-centered tetragonal (fct) structure. The fee FePt, with low coercivity but associated with superparamagnetic properties, may find applications as a magnetic fluid or as a nanoscale carrier for chemical or biochemical species in biomedical areas, while fct FePt is proposed for use in ultrahigh-density magnetic recording applications. However, for both of these applications an enhancement of the intrinsically weak magnetic properties, the avoidance of magnetic interferences from neighbor particles, and the improved stability of the small magnetic body remain key practical issues. We report a simple synthetic method for producing FePt nanoparticles that involves hydrothermal treatment of Fe and Pt precursors in glucose followed by calcination at 900 degrees C. This new method produces thermally stable spheroidal graphite nanoparticles (large and fullerene-like) that encapsulate or decorate FePt particles of ca. 5 nm with no severe macroscopic particle coalescence. Also, a low coercivity of the material is recorded; indicative of small magnetic interference from neighboring carbon-coated particles. Thus, this simple synthetic method involves the use of a more environmentally acceptable glucose/aqueous phase to offer a protective coating for FePt nanoparticles. It is also believed that such a synthetic protocol can be readily extended to the preparation of other graphite-coated magnetic iron alloys of controlled size, stoichiometry, and physical properties.