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.

A design method for superconducting MRI magnets with ferromagnetic material

The emphasis of this work is on the optimal design of MRI magnets with both superconducting coils and ferromagnetic rings. The work is directed to the automated design of MRI magnet systems containing superconducting wire and both `cold' and `warm' iron. Details of the optimization procedure are given and the results show combined superconducting and iron material MRI magnets with excellent field characteristics. Strong, homogeneous central magnetic fields are produced with little stray or external field leakage. The field calculations are performed using a semi-analytical method for both current coil and iron material sources. Design examples for symmetric, open and asymmetric clinical MRI magnets containing both superconducting coils and ferromagnetic material are presented.

Light and electron microscopic analysis of aquaporin 1-like-immunoreactive amacrine cells in the rat retina

Aquaporin 1 (AQP1; also known as CHIP, a channel-forming integral membrane protein of 28 kDa) is the first protein to be shown to function as a water channel and has been recently shown to be present in the rat retina. We previously showed (Kim et al. [1998] Neurosci Lett 244:52-54) that AQP1-like immunoreactivity is present in a certain population of amacrine cells in the rat retina. This study was conducted to characterize these cells in more detail, With immunocytochemistry using specific antisera against AQP1, whole-mount preparations and 50-mum-thick vibratome sections were examined by light and electron microscopy. These cells were a class of amacrine cells, which had symmetric bistratified dendritic trees ramified in stratum 2 and in the border of strata 3 and 4 of the inner plexiform layer (IPL). Their dendritic field diameters ranged from 90 to 230 mum. Double labeling with antisera against AQP1 and gamma-aminobutyric acid or glycine demonstrated that these AQP1-like-immunoreactive amacrine cells were immunoreactive for glycine. Their most frequent synaptic input was from other amacrine cell processes in both sublaminae a and b of the IPL, followed by a few cone bipolar cells. Their primary targets were other amacrine cells and ganglion cells in both sublaminae a and b of the IPL. In addition, synaptic output Onto bipolar cells was rarely observed in sublamina b of the IPL. Thus, the AQP1 antibody labels a class of glycinergic amacrine cells with small to medium-sized dendritic fields in the rat retina. (C) 2002 Wiley-Liss, Inc.

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.

Method for determining the shear stress in cylindrical systems

We develop a method for determining the elements of the pressure tensor at a radius r in a cylindrically symmetric system, analogous to the so-called method of planes used in planar systems [B. D. Todd, Denis J. Evans, and Peter J. Daivis, Phys. Rev. E 52, 1627 (1995)]. We demonstrate its application in determining the radial shear stress dependence during molecular dynamics simulations of the forced flow of methane in cylindrical silica mesopores. Such expressions are useful for the examination of constitutive relations in the context of transport in confined systems.

Mapping the royal road and other hierarchical functions

In this paper we present a technique for visualising hierarchical and symmetric, multimodal fitness functions that have been investigated in the evolutionary computation literature. The focus of this technique is on landscapes in moderate-dimensional, binary spaces (i.e., fitness functions defined over {0, 1}(n), for n less than or equal to 16). The visualisation approach involves an unfolding of the hyperspace into a two-dimensional graph, whose layout represents the topology of the space using a recursive relationship, and whose shading defines the shape of the cost surface defined on the space. Using this technique we present case-study explorations of three fitness functions: royal road, hierarchical-if-and-only-if (H-IFF), and hierarchically decomposable functions (HDF). The visualisation approach provides an insight into the properties of these functions, particularly with respect to the size and shape of the basins of attraction around each of the local optima.

Diálogos sobre a alfabetização, a leitura e a escrita no Programa Provinha Brasil

Esta pesquisa documental analisa as concepções de alfabetização, leitura e escrita subjacentes à Provinha Brasil no período 2008-2012 e o panorama em que esse programa de avaliação é produzido. Parte do referencial bakhtiniano e do conceito de alfabetização de Gontijo (2008, 2013). Ao tomar a Provinha como gênero do discurso, discute os elos precedentes dentro do contexto de produção dessa avaliação, a autoria do Programa e seus principais destinatários. Constata que a Provinha é criada como resposta às demandas de avaliação da alfabetização provenientes de organismos internacionais como o Banco Mundial e a Organização das Nações Unidas para a Educação (Unesco). A avaliação é elaborada pelo Instituto Nacional de Estudos e Pesquisas Educacionais Anísio Teixeira (Inep) como órgão que coordena as avaliações no País, em colaboração com pesquisadores de universidades e de organizações da sociedade civil, para demonstrar confiabilidade científica aliada à participação democrática no processo de produção. Seus principais destinatários são gestores de Secretarias de Educação e professores. Aos primeiros, cabe aderir ao programa de avaliação e tomar medidas administrativas para sua operacionalização nas redes. Os docentes têm o papel central de seguir as orientações do material e reorganizar sua prática em função de melhorias nos desempenhos das crianças no teste. Estas, por sua vez, são desconsideradas como sujeitos de dizeres e é legitimado um discurso homogeneizador sobre seu desenvolvimento. A partir dos testes aplicados e das matrizes de referência e seus eixos, a pesquisa analisa como a diferenciação teórica entre alfabetização e letramento se concretiza na organização das provas. A alfabetização, entendida como apropriação do sistema de escrita, é avaliada no primeiro eixo do teste principalmente como identificação de unidades menores da língua, como letras, sílabas e fonemas. As habilidades de leitura, ligadas ao letramento como concebido nos pressupostos do programa, são aferidas ora como decodificação de palavras e frases descontextualizadas, ora como apreensão de significado predeterminado do texto. A escrita somente é avaliada no ano de 2008 e por meio de itens que solicitavam codificação de palavras e frases ditadas pelo aplicador. Desse modo, a Provinha Brasil contribui para a subtração das potencialidades políticas e transformadoras do aprendizado da língua materna no País.

Estudo da evolução da paisagem da área do município de Linhares (ES) nos anos de 1985 e de 2013/2014, por meio das métricas

Trata a presente pesquisa sobre o estudo da evolução da área da paisagem que compõe o município de Linhares (ES), nos anos de 1985 e 2013/2014, que constitui o maior município do Estado do Espírito Santo. Foi realizada, por meio de processamento digital de imagens de satélites LANDSAT 5 e 8, a classificação de uso e ocupação da terra da área em estudo. Além disso usou-se as imagens do satélite RapdEye para acurácia da classificação digital das imagens. A partir dos resultados levantados de uso e ocupação foram definidas as matrizes da paisagem para 1985 e 2013/2014, bem como avaliadas as manchas que compõem a matriz. Foram aplicadas as métricas da paisagem utilizando a ferramenta de estatística Fragstats, possibilitando o cálculo dos Índices de Paisagem afim de avaliar a evolução qualitativa e quantitativa da paisagem do município de Linhares.

Avaliação da adsorção do herbicida 2,4-D em carvão ativado em pó utilizando água com diferentes qualidades

O 2,4-diclorofenoxiacético (2,4-D) é um dos herbicidas mais consumidos no Brasil e é preferencialmente usado devido a sua boa seletividade e baixo custo. Possui alta toxidade e baixa biodegradabilidade, oferecendo risco à saúde humana e ao meio ambiente, podendo ser encontrado em solos, águas superficiais e subterrâneas. Estudos mostram que o tratamento convencional da água possui baixa eficácia na remoção de microcontaminantes, com isso várias técnicas têm sido utilizadas na remoção de compostos em água, como a adsorção por carvão ativado. Apresenta-se a adsorção em carvão ativado tem se demonstrado como tecnologia eficiente na remoção de diversos contaminantes, dentre eles os agrotóxicos. Assim, o presente trabalho objetivou avaliar a adsorção do 2,4-D por três carvões ativados em pó (CAP) em água ultrapura e em água bruta do Rio Santa Maria da Vitória. A quantificação do herbicida foi analisada por cromatografia líquida de alta eficiência, após concentração da amostra pelo método de extração em fase sólida. Os ensaios de adsorção foram realizados com carvões ativados derivados da casca de coco (CAP-01), pinus (CAP-02) e palha de café (CAP-03), que foram caracterizados e avaliados na sua capacidade de remoção do 2,4-D nas duas matrizes de água. Dois modelos de isoterma de adsorção, Langmuir e Freundlich, foram aplicados para descrever os dados de adsorção, que indicaram o CAP-02 como o carvão que apresentou a melhor capacidade de adsorção do 2,4-D entre os carvões estudados, tanto em água ultrapura quanto em água bruta. Nos ensaios realizados em água bruta, houve redução da adsorção do 2,4-D para as três amostras de CAP, quando comparado com os ensaios realizados em água ultrapura, indicando interferência de compostos, como a matéria orgânica, no processo de adsorção.

Automatic prebent customized prosthesis for pectus excavatum minimally invasive surgery correction

Pectus excavatum is the most common deformity of the thorax. A minimally invasive surgical correction is commonly carried out to remodel the anterior chest wall by using an intrathoracic convex prosthesis in the substernal position. The process of prosthesis modeling and bending still remains an area of improvement. The authors developed a new system, i3DExcavatum, which can automatically model and bend the bar preoperatively based on a thoracic CT scan. This article presents a comparison between automatic and manual bending. The i3DExcavatum was used to personalize prostheses for 41 patients who underwent pectus excavatum surgical correction between 2007 and 2012. Regarding the anatomical variations, the soft-tissue thicknesses external to the ribs show that both symmetric and asymmetric patients always have asymmetric variations, by comparing the patients’ sides. It highlighted that the prosthesis bar should be modeled according to each patient’s rib positions and dimensions. The average differences between the skin and costal line curvature lengths were 84 ± 4 mm and 96 ± 11 mm, for male and female patients, respectively. On the other hand, the i3DExcavatum ensured a smooth curvature of the surgical prosthesis and was capable of predicting and simulating a virtual shape and size of the bar for asymmetric and symmetric patients. In conclusion, the i3DExcavatum allows preoperative personalization according to the thoracic morphology of each patient. It reduces surgery time and minimizes the margin error introduced by the manually bent bar, which only uses a template that copies the chest wall curvature.