Review of recent results on excursion set models

Many images consist of two or more 'phases', where a phase is a collection of homogeneous zones. For example, the phases may represent the presence of different sulphides in an ore sample. Frequently, these phases exhibit very little structure, though all connected components of a given phase may be similar in some sense. As a consequence, random set models are commonly used to model such images. The Boolean model and models derived from the Boolean model are often chosen. An alternative approach to modelling such images is to use the excursion sets of random fields to model each phase. In this paper, the properties of excursion sets will be firstly discussed in terms of modelling binary images. Ways of extending these models to multi-phase images will then be explored. A desirable feature of any model is to be able to fit it to data reasonably well. Different methods for fitting random set models based on excursion sets will be presented and some of the difficulties with these methods will be discussed.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

Tetrachromacy, oil droplets and bird plumage colours

There is a growing body of data on avian eyes, including measurements of visual pigment and oil droplet spectral absorption, and of receptor densities and their distributions across the retina. These data are sufficient to predict psychophysical colour discrimination thresholds for light-adapted eyes, and hence provide a basis for relating eye design to visual needs. We examine the advantages of coloured oil droplets, UV vision and tetrachromacy for discriminating a diverse set of avian plumage spectra under natural illumination. Discriminability is enhanced both by tetrachromacy and coloured oil droplets. Oil droplets may also improve colour constancy. Comparison of the performance of a pigeon's eye, where the shortest wavelength receptor peak is at 410 nm, with that of the passerine Leiothrix, where the ultraviolet-sensitive peak is at 365 nm, generally shows a small advantage to the latter, but this advantage depends critically on the noise level in the sensitivity mechanism and on the set of spectra being viewed.

Structural basis of autoregulation of phenylalanine hydroxylase

Phenylalanine hydroxylase converts phenylalanine to tyrosine, a rate-limiting step in phenylalanine catabolism and protein and neurotransmitter biosynthesis. It is tightly regulated by the substrates phenylalanine and tetrahydrobiopterin and by phosphorylation. We present the crystal structures of dephosphorylated and phosphorylated forms of a dimeric enzyme with catalytic and regulatory properties of the wild-type protein. The structures reveal a catalytic domain flexibly linked to a regulatory domain. The latter consists of an N-terminal autoregulatory sequence (containing Ser 16, which is the site of phosphorylation) that extends over the active site pocket, and an alpha-beta sandwich core that is, unexpectedly, structurally related to both pterin dehydratase and the regulatory domains of metabolic enzymes. Phosphorylation has no major structural effects in the absence of phenylalanine, suggesting that phenylalanine and phosphorylation act in concert to activate the enzyme through a combination of intrasteric and possibly allosteric mechanisms.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

An "openable," high-strength gradient set for orthopedic MRI

A novel three-axis gradient set and RF resonator for orthopedic MRT has been designed and constructed. The set is openable and may be wrapped around injured joints. The design methodology used was the minimization of magnetic field spherical harmonics by simulated annealing. Splitting of the longitudinal coil presents the major design challenge to a fully openable gradient set and in order to efficiently design such coils, we have developed a new fast algorithm for determining the magnetic field spherical harmonics generated by an are of multiturn wire. The algorithm allows a realistic impression of the effect of split longitudinal designs. A prototype set was constructed based on the new designs and tested in a 2-T clinical research system. The set generated 12 mT/m/A with a linear region of 12 cm and a switching time of 100 mu s, conforming closely with theoretical predictions. Preliminary images from the set are presented. (C) 1999 Academic Press.

Designs for an asymmetric gradient set and a compact superconducting magnet for neural magnetic resonance imaging

Imaging of the head and neck is the most commonly performed clinical magnetic resonance imaging (MRI) examination [R. G. Evans and J. R. G. Evans, AJR 157, 603 (1991)]. This is usually undertaken in a generalist MRI instrument containing superconducting magnet system capable of imaging all organs. These generalist instruments are large, typically having a bore of 0.9-1.0 m and a length of 1.7-2.5 m and therefore are expensive to site, somewhat claustrophobic to the patient, and offer little access by attending physicians. In this article, we present the design of a compact, superconducting MRI magnet for head and neck imaging that is less than 0.8 m in length and discuss in detail the design of an asymmetric gradient coil set, tailored to the magnet profile. In particular, the introduction of a radio-frequency FM modulation scheme in concert with a gradient sequence allows the epoch of the linear region of the gradient set to be much closer to the end of the gradient structure than was previously possible. Images from a prototype gradient set demonstrate the effectiveness of the designs. (C) 1999 American Institute of Physics. [S0034-6748(99)04910-2].

Approach to designing rotating drum bioreactors for solid-state fermentation on the basis of dimensionless design factors

The development of large-scale solid-stale fermentation (SSF) processes is hampered by the lack of simple tools for the design of SSF bioreactors. The use of semifundamental mathematical models to design and operate SSF bioreactors can be complex. In this work, dimensionless design factors are used to predict the effects of scale and of operational variables on the performance of rotating drum bioreactors. The dimensionless design factor (DDF) is a ratio of the rate of heat generation to the rate of heat removal at the time of peak heat production. It can be used to predict maximum temperatures reached within the substrate bed for given operational variables. Alternatively, given the maximum temperature that can be tolerated during the fermentation, it can be used to explore the combinations of operating variables that prevent that temperature from being exceeded. Comparison of the predictions of the DDF approach with literature data for operation of rotating drums suggests that the DDF is a useful tool. The DDF approach was used to explore the consequences of three scale-up strategies on the required air flow rates and maximum temperatures achieved in the substrate bed as the bioreactor size was increased on the basis of geometric similarity. The first of these strategies was to maintain the superficial flow rate of the process air through the drum constant. The second was to maintain the ratio of volumes of air per volume of bioreactor constant. The third strategy was to adjust the air flow rate with increase in scale in such a manner as to maintain constant the maximum temperature attained in the substrate bed during the fermentation. (C) 2000 John Wiley & Sons, Inc.