Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Bayesian analysis of the error correction model

This paper presents a method for estimating the posterior probability density of the cointegrating rank of a multivariate error correction model. A second contribution is the careful elicitation of the prior for the cointegrating vectors derived from a prior on the cointegrating space. This prior obtains naturally from treating the cointegrating space as the parameter of interest in inference and overcomes problems previously encountered in Bayesian cointegration analysis. Using this new prior and Laplace approximation, an estimator for the posterior probability of the rank is given. The approach performs well compared with information criteria in Monte Carlo experiments. (C) 2003 Elsevier B.V. All rights reserved.

Quantum phase-space analysis of the pendular cavity

We perform a quantum-mechanical analysis of the pendular cavity, using the positive-P representation, showing that the quantum state of the moving mirror, a macroscopic object, has noticeable effects on the dynamics. This system has previously been proposed as a candidate for the quantum-limited measurement of small displacements of the mirror due to radiation pressure, for the production of states with entanglement between the mirror and the field, and even for superposition states of the mirror. However, when we treat the oscillating mirror quantum mechanically, we find that it always oscillates, has no stationary steady state, and exhibits uncertainties in position and momentum which are typically larger than the mean values. This means that previous linearized fluctuation analyses which have been used to predict these highly quantum states are of limited use. We find that the achievable accuracy in measurement is fat, worse than the standard quantum limit due to thermal noise, which, for typical experimental parameters, is overwhelming even at 2 mK

Non-diagonal solutions of the reflection equation for the trigonometric A((1)) (n-1) vertex model

We obtain a class of non-diagonal solutions of the reflection equation for the trigonometric A(n-1)((1)) vertex model. The solutions can be expressed in terms of intertwinner matrix and its inverse, which intertwine two trigonometric R-matrices. In addition to a discrete (positive integer) parameter l, 1 less than or equal to l less than or equal to n, the solution contains n + 2 continuous boundary parameters.

Existence results for partial neutral functional differential equations with state-dependent delay

In this paper we study the existence of mild solutions for a class of first order abstract partial neutral differential equations with state-dependent delay. (C) 2008 Elsevier Ltd. All rights reserved.

Use of pleopod morphology to determine sexual dimorphism and maturity in hermit crabs: Isocheles sawayai as a model

In the Anomura, studies on growth patterns are infrequent, possibly because the heterogeneity of the group, especially in terms of morphology, makes it difficult to construct generalized growth models. Particularly hermit crabs are an interesting group to evaluate aspects of growth, because of their unique body. Isocheles sawayai, a hermit crab found only in the western Atlantic Ocean, poorly known with respect to its sexual dimorphism and maturity, was investigated here based on morphometry. Monthly collections (July 2001 through June 2003) were made from a shrimp fishing boat in the Caraguatatuba region on the northern coast of the state of SA o pound Paulo, Brazil. The specimens were measured and weighed, and had their sex checked. Throughout the sampling period, 374 specimens of I. sawayai were collected (11.23% nonovigerous females, 6.69% ovigerous females, 79.41% males and 2.67% intersexes). The size at which morphological sexual maturity was reached by both sexes ranged from 4.0 to 4.3 mm shield length, according to the relative growth and the size of the smallest ovigerous female. Sexual dimorphism was shown by males, which were significantly larger than females, and by differences in growth pattern between the sexes, especially for relationships that involved the pleopods, which is related to their different functions in males and females. The present study is one of the first to use pleopod morphometry to determine sexual maturity and dimorphism in hermit crabs, especially for species with intersexuality such as I. sawayai.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

Algebraic Bethe ansatz for the anisotropic supersymmetric U model

We present an algebraic Bethe ansatz for the anisotropic supersymmetric U model for correlated electrons on the unrestricted 4(L)-dimensional electronic Hilbert space x(n=l)(L)C(4)(where L is the lattice length). The supersymmetry algebra of the local Hamiltonian is the quantum superalgebra U-q[gl(2\1)] and the model contains two symmetry-preserving free real parameters; the quantization parameter q and the Hubbard interaction parameter U. The parameter U arises from the one-parameter family of inequivalent typical four-dimensional irreps of U-q[gl(2\1)]. Eigenstates of the model are determined by the algebraic Bethe ansatz on a one-dimensional periodic lattice.

Pure two-atom squeezed state in spontaneous emission from two non-identical atoms

We calculate the stationary state of the system of two non-identical two-level atoms driven by a finite-bandwidth two-mode squeezed vacuum. It is well known that two identical two-level atoms driven by a broadband squeezed vacuum may decay to a pure state, called the pure two-atom squeezed state, and that the presence of the antisymmetric state can change its purity. Here, we show that for small interatomic separations the stationary state of two non-identical atoms is not sensitive to the presence of the antisymmetric state and is the pure two-atom squeezed state. This effect is a consequence of the fact that in the system of two non-identical atoms the antisymmetric state is no longer the trapping state. We also calculate the squeezing properties of the emitted field and find that the squeezing spectrum of the output field may exhibit larger squeezing than that in the input squeezed vacuum. Moreover, we show that squeezing in the total field attains the optimum value which can ever be achieved in the field emitted by two atoms.

Influence of Estrogen Deficiency on Bone Around Osseointegrated Dental Implants: An Experimental Study in the Rat Jaw Model

Purpose: The aim of this study was to evaluate the influence of estrogen deficiency on bone around osseointegrated dental implants in a rat jaw model. Materials and Methods: This study used 16 female rats that had the first molars bilaterally extracted and were allowed to heal for 30 days before implant placement. Sixty days after implant placement, the animals were randomly subjected to sham surgery or ovariectomy (OVX). The animals were euthanized 90 days after OVX. Bone-to-implant contact, bone area fraction occupancy between implant threads, mineral density, turnover markers, and cells positive for tartrate-resistant acid phosphatase were assessed for the 2 groups. Results: The results showed that OVX group presented a decrease of systemic bone density, alterations in bone turnover markers, and an increase of cells positive for tartrate-resistant acid phosphatase compared with the sham-surgery group. However, no difference relative to bone-to-implant contact and bone area fraction occupancy was observed between groups. Conclusions: The findings of this study demonstrate that estrogen deficiency may not be considered a risk factor for osseointegrated implant failure in jaw bone. (C) 2011 American Association of Oral and Maxillofacial Surgeons J Oral Maxillofac Surg 69:1911-1918, 2011

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

An isolated perfused human placental lobule model for multiple indicator dilution studies

We describe a method for multiple indicator dilution studies in the isolated perfused human placental lobule developed to investigate the relationships between changes in pressure and flow and solute clearance. A peripheral lobule of a human placenta is perfused with a tissue culture-based medium and the perfusate oxygen tension, arterial and venous pressures, pH and perfusion temperature continuously monitored by a computerized system. Flow rates are readily changed. Bolus injections of vascular, extracellular and water space markers, and study compounds can be made into either maternal or fetal circulations, and precisely timed outflow fractions can be collected with computer-controlled fraction collectors, allowing simultaneous determination of concentration-time profiles of each marker. (C) 1997 Elsevier Science Inc.

Cell cycle model to describe animal cell size variation and lag between cell number and biomass dynamics

The use of cell numbers rather than mass to quantify the size of the biotic phase in animal cell cultures causes several problems. First, the cell size varies with growth conditions, thus yields expressed in terms of cell numbers cannot be used in the normal mass balance sense. Second, experience from microbial systems shows that cell number dynamics lag behind biomass dynamics. This work demonstrates that this lag phenomenon also occurs in animal cell culture. Both the lag phenomenon and the variation in cell size are explained using a simple model of the cell cycle. The basis for the model is that onset of DNA synthesis requires accumulation of G1 cyclins to a prescribed level. This requirement is translated into a requirement for a cell to reach a critical size before commencement of DNA synthesis. A slower gl-owing cell will spend more time in G1 before reaching the critical mass. In contrast, the period between onset of DNA synthesis and mitosis, tau(B), is fixed. The two parameters in the model, the critical size and tau(B), were determined from eight steady-state measurements of mean cell size in a continuous hybridoma culture. Using these parameters, it was possible to predict with reasonable accuracy the transient behavior in a separate shift-up culture, i.e., a culture where cells were transferred from a lean environment to a rich environment. The implications for analyzing experimental data for animal cell culture are discussed. (C) 1997 John Wiley & Sons, Inc.