An experimental and finite element poroelastic creep response analysis of an intervertebral hydrogel disc model in axial compression

A hydrogel intervertebral disc (lVD) model consisting of an inner nucleus core and an outer anulus ring was manufactured from 30 and 35% by weight Poly(vinyl alcohol) hydrogel (PVA-H) concentrations and subjected to axial compression in between saturated porous endplates at 200 N for 11 h, 30 min. Repeat experiments (n = 4) on different samples (N = 2) show good reproducibility of fluid loss and axial deformation. An axisymmetric nonlinear poroelastic finite element model with variable permeability was developed using commercial finite element software to compare axial deformation and predicted fluid loss with experimental data. The FE predictions indicate differential fluid loss similar to that of biological IVDs, with the nucleus losing more water than the anulus, and there is overall good agreement between experimental and finite element predicted fluid loss. The stress distribution pattern indicates important similarities with the biological lVD that includes stress transference from the nucleus to the anulus upon sustained loading and renders it suitable as a model that can be used in future studies to better understand the role of fluid and stress in biological IVDs. (C) 2005 Springer Science + Business Media, Inc.

Excitation spectra of the spin-1/2 triangular-lattice Heisenberg antiferromagnet

We use series expansion methods to calculate the dispersion relation of the one-magnon excitations for the spin-(1)/(2) triangular-lattice nearest-neighbor Heisenberg antiferromagnet above a three-sublattice ordered ground state. Several striking features are observed compared to the classical (large-S) spin-wave spectra. Whereas, at low energies the dispersion is only weakly renormalized by quantum fluctuations, significant anomalies are observed at high energies. In particular, we find rotonlike minima at special wave vectors and strong downward renormalization in large parts of the Brillouin zone, leading to very flat or dispersionless modes. We present detailed comparison of our calculated excitation energies in the Brillouin zone with the spin-wave dispersion to order 1/S calculated recently by Starykh, Chubukov, and Abanov [Phys. Rev. B74, 180403(R) (2006)]. We find many common features but also some quantitative and qualitative differences. We show that at temperatures as low as 0.1J the thermally excited rotons make a significant contribution to the entropy. Consequently, unlike for the square lattice model, a nonlinear sigma model description of the finite-temperature properties is only applicable at temperatures < 0.1J. Finally, we review recent NMR measurements on the organic compound kappa-(BEDT-TTF)(2)Cu-2(CN)(3). We argue that these are inconsistent with long-range order and a description of the low-energy excitations in terms of interacting magnons, and that therefore a Heisenberg model with only nearest-neighbor exchange does not offer an adequate description of this material.

The lattice solid model to simulate the physics of rocks and earthquakes: Incorporation of friction

The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity-Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick-slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics. (C) 1999 Academic Press.

Characterizing mixing and measurement in quantum mechanics

What fundamental constraints characterize the relationship between a mixture rho = Sigma (i)p(i)rho (i) of quantum states, the states rho (i) being mixed, and the probabilities p(i)? What fundamental constraints characterize the relationship between prior and posterior states in a quantum measurement? In this paper we show that then are many surprisingly strong constraints on these mixing and measurement processes that can be expressed simply in terms of the eigenvalues of the quantum states involved. These constraints capture in a succinct fashion what it means to say that a quantum measurement acquires information about the system being measured, and considerably simplify the proofs of many results about entanglement transformation.

On minimal models of process systems

Minimal representations are known to have no redundant elements, and are therefore of great importance. Based on the notions of performance and size indices and measures for process systems, the paper proposes conditions for a process model being minimal in a set of functionally equivalent models with respect to a size norm. Generalized versions of known procedures to obtain minimal process models for a given modelling goal, model reduction based on sensitivity analysis and incremental model building are proposed and discussed. The notions and procedures are illustrated and compared on a simple example, that of a simple nonlinear fermentation process with different modelling goals and on a case study of a heat exchanger modelling. (C) 2004 Elsevier Ltd. All rights reserved.

Evaluation of different models to describe egg and pupal development of Xyleborus fornicatus Eichh. (Coleoptera: Scolytidae), the shot-hole borer of tea in Sri Lanka

Abstract Development data of eggs and pupae of Xyleborus fornicatus Eichh. (Coleoptera: Scolytidae), the shot-hole borer of tea in Sri Lanka, at constant temperatures were used to evaluate a linear and seven nonlinear models for insect development. Model evaluation was based on fit to data (residual sum of squares and coefficient of determination or coefficient of nonlinear regression), number of measurable parameters, the biological value of the fitted coefficients and accuracy in the estimation of thresholds. Of the nonlinear models, the Lactin model fitted experimental data well and along with the linear model, can be used to describe the temperature-dependent development of this species.

Investigating the performance analysis of EASI algorithm and EKENS algorithm in nonlinear model

This paper investigates the performance of EASI algorithm and the proposed EKENS algorithm for linear and nonlinear mixtures. The proposed EKENS algorithm is based on the modified equivariant algorithm and kernel density estimation. Theory and characteristic of both the algorithms are discussed for blind source separation model. The separation structure of nonlinear mixtures is based on a nonlinear stage followed by a linear stage. Simulations with artificial and natural data demonstrate the feasibility and good performance of the proposed EKENS algorithm.

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.

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.

Quantum transport and integrability of the Anderson model for a quantum dot with multiple leads

We show that an Anderson Hamiltonian describing a quantum dot connected to multiple leads is integrable. A general expression for the nonlinear conductance is obtained by combining the Bethe ansatz exact solution with Landauer-Buttiker theory. In the Kondo regime, a closed form expression is given for the matrix conductance at zero temperature and when all the leads are close to the symmetric point. A bias-induced splitting of the Kondo resonance is possible for three or more leads. Specifically, for N leads, with each at a different chemical potential, there can be N-1 Kondo peaks in the conductance.

Nonlinear quantum optical computing via measurement

We show how the measurement induced model of quantum computation proposed by Raussendorf and Briegel ( 2001, Phys. Rev. Letts., 86, 5188) can be adapted to a nonlinear optical interaction. This optical implementation requires a Kerr nonlinearity, a single photon source, a single photon detector and fast feed forward. Although nondeterministic optical quantum information proposals such as that suggested by KLM ( 2001, Nature, 409, 46) do not require a Kerr nonlinearity they do require complex reconfigurable optical networks. The proposal in this paper has the benefit of a single static optical layout with fixed device parameters, where the algorithm is defined by the final measurement procedure.

Calibration of the channel that determines the omega-hydroxylation regiospecificity of cytochrome P4504A1 - Catalytic oxidation of 12-halododecanoic acids

The fatty acid omega-hydroxylation regiospecificity of CYP4 enzymes may result from presentation of the terminal carbon to the oxidizing species via a narrow channel that restricts access to the other carbon atoms. To test this hypothesis, the oxidation of 12-iodo-, 12-bromo-, and 12-chlorododecanoic acids by recombinant CYP4A1 has been examined. Although all three 12-halododecanoic acids bind to CYP4A1 with similar dissociation constants, the 12-chloro and 12-bromo fatty acids are oxidized to 12-hydroxydodecanoic acid and 12-oxododecanoic acid, whereas the 12-iodo analogue is very poorly oxidized. Incubations in (H2O)-O-18 show that the 12-hydroxydodecanoic acid oxygen derives from water, whereas that in the aldehyde derives from O-2. The alcohol thus arises from oxidation of the halide to an oxohalonium species that is hydrolyzed by water, whereas the aldehyde arises by a conventional carbon hydroxylation-elimination mechanism. No irreversible inactivation of CYP4A1 is observed during 12-halododecanoic acid oxidation. Control experiments show that CYP2E1, which has an omega-1 regiospecificity, primarily oxidizes 12-halododecanoic acids to the omega-aldehyde rather than alcohol product. Incubation of CYP4A1 with 12,12-[H-2](2)-12-chlorododecanoic acid causes a 2-3-fold increase in halogen versus carbon oxidation. The fact that the order of substrate oxidation (Br > Cl >> I) approximates the inverse of the intrinsic oxidizability of the halogen atoms is consistent with presentation of the halide terminus via a channel that accommodates the chloride and bromide but not iodide atoms, which implies an effective channel diameter greater than 3.90 angstrom but smaller than 4.30 angstrom.

Vitamin B-12 release from P(HEMA-co-THFMA) in water and SBF: A model drug release study

A model drug release study on the ingress of water and Kokubo simulated body fluid (SBF) into poly(2-hydroxyethyl methacrylate) (THFMA) and its copolymers with tetrahydrofurfuryl methacrylate (THFMA) loaded with vitamin B-12 was undertaken over the temperature range 298-318 K. The polymers were studied as cylinders and were loaded with either 5 or 10 wt-% of the drug. The drug release from the polymers was found to follow a Fickian diffusion mechanism in the early stages of the drug release, with higher normalized release rates at higher temperatures and higher drug loadings. The normalized release rates were also found to be higher for the SBF solution than for water. The copolymer composition was found to have a significant effect on the rate of release of the drug, with the rate falling rapidly between HEMA mole fractions of 1.0 and 0.8, but for lower mole fractions of HEMA the normalized release rate decreased more slowly. This behaviour followed the trend found for the changes in the equilibrium penetrant contents for the copolymers.

General foundation for the nonlinear ponderomotive four-force in laser-plasma interactions

The interaction of electromagnetic radiation with plasmas is studied in relativistic four-vector formalism. A gauge and Lorentz invariant ponderomotive four-force is derived from the time dependent nonlinear three-force of Hora (1985). This four-force, due to its Lorentz invariance, contains new magnetic field terms. A new gauge and Lorentz invariant model of the response of plasma to electromagnetic radiation is then devised. An expression for the dispersion relation is obtained from this model. It is then proved that the magnetic permeability of plasma is unity for a general reference frame. This is an important result since it has been previously assumed in many plasma models.

Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project

To foster ongoing international cooperation beyond ACES (APEC Cooperation for Earthquake Simulation) on the simulation of solid earth phenomena, agreement was reached to work towards establishment of a frontier international research institute for simulating the solid earth: iSERVO = International Solid Earth Research Virtual Observatory institute (http://www.iservo.edu.au). This paper outlines a key Australian contribution towards the iSERVO institute seed project, this is the construction of: (1) a typical intraplate fault system model using practical fault system data of South Australia (i.e., SA interacting fault model), which includes data management and editing, geometrical modeling and mesh generation; and (2) a finite-element based software tool, which is built on our long-term and ongoing effort to develop the R-minimum strategy based finite-element computational algorithm and software tool for modelling three-dimensional nonlinear frictional contact behavior between multiple deformable bodies with the arbitrarily-shaped contact element strategy. A numerical simulation of the SA fault system is carried out using this software tool to demonstrate its capability and our efforts towards seeding the iSERVO Institute.