Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.

Tissue engineering for bone regeneration using differentiated alveolar bone cells in collagen scaffolds

Regeneration of osseous defects by a tissue-engineering approach provides a novel means of treatment utilizing cell biology, materials science, and molecular biology. In this study the concept of tissue engineering was tested with collagen type I matrices seeded with cells with osteogenic potential and implanted into sites where osseous damage had occurred. Explant cultures of cells from human alveolar bone and gingiva were established. When seeded into a three-dimensional type I collagen-based scaffold, the bone-derived cells maintained their osteoblastic phenotype as monitored by mRNA and protein levels of the bone-related proteins including bone sialoprotein, osteocalcin, osteopontin, bone morphogenetic proteins 2 and 4, and alkaline phosphatase. These in vitro-developed matrices were implanted into critical-size bone defects in skulls of immunodeficient (SCID) mice. Wound healing was monitored for up to 4 weeks. When measured by microdensitometry the bone density within defects filled with osteoblast-derived matrix was significantly higher compared with defects filled with either collagen scaffold alone or collagen scaffold impregnated with gingival fibroblasts. New bone formation was found at all the sites treated with the osteoblast-derived matrix at 28 days, whereas no obvious new bone formation was identified at the same time point in the control groups. In situ hybridization for the human-specific Alu gene sequence indicated that the newly formed bone tissue resulted from both transplanted human osteoblasts and endogenous mesenchymal stem cells. The results indicate that cells derived from human alveolar bone can be incorporated into bioengineered scaffolds and synthesize a matrix, which on implantation can induce new bone formation.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

Ultrastructure of developing tooth plates in the Australian lungfish, Neoceratodus forsteri (Osteichthyes : Dipnoi)

While the lungfish dentition is partially understood as far as morphology and light microscopic structure is concerned, the ultrastructure is not. Each tooth plate is associated with a dental lamina that develops from the inner layer of endodermal cells that form the oral epithelium. Dentines, bone and cartilage of the jaws differentiate from mesenchyme cells aggregating beneath the oral endothelium. Enamel, in the developing and in the mature form, has similarities to that of other early vertebrates, but unusual characters appear as development proceeds. Ameloblasts are capable of secreting enamel, and, with mononuclear osteoclasts, of remodelling the bone below the tooth plate. The forms of dentine, all based largely on an extracellular matrix of collagen and mineralised with biological apatite, differ from each other and from the underlying bone in the ultrastructure of associated cells and in the mineralised extracellular matrices produced. Cell processes emerging from the odontoblasts and from the osteoblasts vary in length, degree of branching and of anastomoses between the processes, although all of the cell types have large amounts of rough endoplasmic reticulum. Mineralisation of the extracellular matrices varies among the enamel, dentines and bone in the tooth plate. In addition, the development of the hard tissues of the tooth plates indicates that many of the similarities in fine structure of the dentition in lungfish, to tissues in other fish and amphibia, apparent early in development, disappear as the dentition matures. (C) 2003 Elsevier Ltd. All rights reserved.

N-aminopyrroledione-hydrazonoketene-pyrazolium oxide-pyrazolone rearrangements and pyrazolone tautomerism

Flash vacuum thermolysis (FVT) of 1-(dimethylamino)pyrrole-2,3-diones 5 causes extrusion of CO with formation of transient hydrazonoketenes 7. The transient ketenes 7 are observable in the form of weak bands at 2130 (7a) or 2115 cm(-1) (7b) in the Ar matrix IR spectra resulting from either FVT or photolysis of either 5 or 1,1- dimethylpyrazolium-5- oxides 8, and these absorptions are in excellent agreement with B3LYP/6-31G* frequency calculations. Under FVT conditions the ketenes 7 cyclize to pyrazolium oxides 8, which undergo 1,4-migration of a methyl group to yield 1,4-dimethyl-3-phenylpyrazole-5(4H)-one 9a and 1,4,4-trimethyl-3-phenylpyrazole-5(4H)-one 9b. All three tautomers of 9a have been characterized, viz. the CH form 9a (most stable form in the gas phase, the solid state and solvents of low polarity), the OH form 9a' (metastable solid at room temperature) and the NH form 9a (stable in aprotic dipolar solvents). The isomeric 1,4-dimethyl-5-phenylpyrazole-3(2H)-one 12 tautomerizes to the 3-hydroxypyrazole 12'. The crystal structure of the hydrochloride 14 of 9a'/9a is reported, representing the first structurally characterised example of a protonated 5-hydroxypyrazole.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.