Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Gas permeation in silicon-oxide/polymer (SiOx/PET) barrier films: role of the oxide lattice, nano-defects and macro-defects

We propose a model for permeation in oxide coated gas barrier films. The model accounts for diffusion through the amorphous oxide lattice, nano-defects within the lattice, and macro-defects. The presence of nano-defects indicate the oxide layer is more similar to a nano-porous solid (such as zeolite) than silica glass with respect to permeation properties. This explains why the permeability of oxide coated polymers is much greater, and the activation energy of permeation much lower, than values expected for polymers coated with glass. We have used the model to interpret permeability and activation energies measured for the inert gases (He, Ne and Ar) in evaporated SiOx films of varying thickness (13-70 nm) coated on a polymer substrate. Atomic force and scanning electron microscopy were used to study the structure of the oxide layer. Although no defects could be detected by microscopy, the permeation data indicate that macro-defects (>1 nm), nano-defects (0.3-0.4 nm) and the lattice interstices (<0.3 nm) all contribute to the total permeation. (C) 2002 Elsevier Science B.V. All rights reserved.

Quasimode theory of quantum optical processes in photonic band gap materials

This paper deals with atomic systems coupled to a structured reservoir of quantum EM field modes, with particular relevance to atoms interacting with the field in photonic band gap materials. The case of high Q cavities has been treated elsewhere using Fano diagonalization based on a quasimode approach, showing that the cavity quasimodes are responsible for pseudomodes introduced to treat non-Markovian behaviour. The paper considers a simple model of a photonic band gap case, where the spatially dependent permittivity consists of a constant term plus a small spatially periodic term that leads to a narrow band gap in the spectrum of mode frequencies. Most treatments of photonic band gap materials are based on the true modes, obtained numerically by solving the Helmholtz equation for the actual spatially periodic permittivity. Here the field modes are first treated in terms of a simpler quasimode approach, in which the quasimodes are plane waves associated with the constant permittivity term. Couplings between the quasimodes occur owing to the small periodic term in the permittivity, with selection rules for the coupled modes being related to the reciprocal lattice vectors. This produces a field Hamiltonian in quasimode form. A matrix diagonalization method may be applied to relate true mode annihilation operators to those for quasimodes. The atomic transitions are coupled to all the quasimodes, and the true mode atom-EM field coupling constants (one-photon Rabi frequencies) are related to those for the quasimodes and also expressions are obtained for the true mode density. The results for the one-photon Rabi frequencies differ from those assumed in other work. Expressions for atomic decay rates are obtained using the Fermi Golden rule, although these are valid only well away from the band gaps.

Dynamic Behavior of the Jahn-Teller distorted Cu(H2O)(6)(2+) ion in Cu2+ doped Cs-2[Zn(H2O)(6)](ZrF6)(2) and the crystal structure of the host lattice

The temperature dependence of the X- and Q-band EPR spectra of Cs-2[Zn(H2O)(6)](ZrF6)(2) containing similar to1% Cu2+ is reported. All three molecular g-values vary with temperature, and their behavior is interpreted using a model in which the potential surface of the Jahn-Teller distorted Cu(H2O)(6)(2+) ion is perturbed by an orthorhombic strain induced by interactions with the surrounding lattice. The strain parameters are significantly smaller than those reported previously for the Cu(H2O)(6)(2+) ion in similar lattices. The temperature dependence of the two higher g-values suggests that in the present compound the lattice interactions change slightly with temperature. The crystal structure of the Cs-2[Zn(H2O)(6)](ZrF6)(2) host is reported, and the geometry of the Zn(H2O)(6)(2+) ion is correlated with lattice strain parameters derived from the EPR spectrum of the guest Cu2+ complex.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Analysis of adsorption data of graphitized thermal carbon black with a DFT-lattice gas theory

In this paper we analyzed the adsorption of gases and vapors on graphitised thermal carbon black by using a modified DFT-lattice theory, in which we assume that the behavior of the first layer in the adsorption film is different from those of second and higher layers. The effects of various parameters on the topology of the adsorption isotherm were first investigated, and the model was then applied in the analysis of adsorption data of numerous substances on carbon black. We have found that the first layer in the adsorption film behaves differently from the second and higher layers in such a way that the adsorbate-adsorbate interaction energy in the first layer is less than that of second and higher layers, and the same is observed for the partition function. Furthermore, the adsorbate-adsorbate and adsorbate-adsorbent interaction energies obtained from the fitting are consistently lower than the corresponding values obtained from the viscosity data and calculated from the Lorentz-Berthelot rule, respectively.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Responsiveness, affinity constants and beta-adrenoceptor reserves for isoprenaline on aortae from normo-, pre and hypertensive rats

The objective of this study was to determine the responsiveness, affinity constants and beta-adrenoceptor reserves for isoprenaline on the isolated aorta in the maturation of normotensive and hypertensive rats. The effects of a very slowly reversible antagonist, bromoacetylalprenololmenthane (BAAM), on the relaxant responses of the aortae of 5- and 14-week-old Wistar Kyoto normotensive rats (WKY) and spontaneously hypertensive rats (SHRs) to isoprenaline were determined. Five-week-old SHRs are pre-hypertensive and the aortic rings are less responsive to isoprenaline than age-matched WKY (pD(2) values: WKY, 8.40; SHRs, 8.03). Similar relaxant responses to forskolin were obtained on the aortae of 5- and 14-week-old WKY and SHRs. The K-A value for isoprenaline at the aortic beta(2)-adrenoceptors of the 5-week-old WKY was 2.1 x 10(-7) M, and similar values were obtained on the aortae of 5-week-old SHR and 14-week-old WKY and SHRs. In the maturation of the WKY aortae from 5 to 14 weeks, there was a reduction in the maximum response, a major loss of sensitivity and a loss of 2-adrenoceptor reserve for isoprenaline. On 5-week-old SHR aorta, the sensitivity to isoprenaline was 2.5-fold lower, and the beta(2)-adrenoceptor reserve was less than on age-matched WKY. In the development of hypertension on the SHR aorta from 5 to 14 weeks, there was a reduction in the maximum response to isoprenaline. At 14 weeks, the sensitivity and the 2-adrenoceptor reserve to isoprenaline were similar, but the maximum responses were lower on the SHR than WKY. As there are differences in pre-hypertensive SHR and age-matched WKY aortic responses to isoprenaline, it is no longer valid to consider that the loss of responsiveness to isoprenaline in hypertension is solely owing to the hypertension. There are no changes in affinity, but major changes in the sensitivity, maximum responses and aortic beta(2)-adrenoceptor reserves to isoprenaline in the maturation of normotensive and pre-hypertensive aortae.

Effect of Co addition on the lattice parameter, electrical conductivity and sintering of gadolinia-doped ceria

Co-sintering aid has been added to Ce1.9Gd0.1O1.95 (CGO) by treating a commercial powder with Co(NO3)(2) (COCGO), X-ray diffraction (XRD) measurements of lattice parameter indicated that the Co was located on the CGO particle surface after calcination at 650 degreesC. After heat treatment at temperatures above 650 degreesC, the room temperature lattice parameter of CGO was found to increase, indicating redistribution of the Gd. Compared to CGO, the lattice parameter of CGO + 2 cation% Co (2CoCGO) was lower for a given temperature (650-1100 degreesC), A.C. impedance revealed that the lattice conductivity of 2CoCGO was enhanced when densified at lower temperatures, Transmission electron microscopy (TEM) showed that, even after sintering for 4 h at 980 degreesC, most of the Co was located at grain boundaries. (C) 2002 Published by Elsevier Science B.V.

Structural variations and formation constants of first-row transition metal complexes of biologically active aroylhydrazones

Iron chelators of the 2-pyridinecarbaldehyde isonicotinoylhydrazone (HPCIH) class show high potential for the treatment of iron overload diseases. In the present study, selected first-row transition metal (from Mn to Zn) complexes with HPCIH and 2-pyridinecarbaldehyde (4'-aminobenzoyl)hydrazone (HPCAH) were synthesised and characterised. Crystallography reveals that HPCAH exclusively forms bis complexes with divalent transition metals, with each ligand coordinating meridionally through its pyridine-N, imine-N and carbonyl-O atoms, forming distorted octahedral cis-MN4O2 complexes. Complexes of HPCIH were more varied and unpredictable, with metal/ligand ratios of 1:1, 1:2, 2:2 and 3:2 obtained with different metal ions. The isonicotinoyl ring N-atom in HPCIH was found to be an effective ligand, and this resulted in the varied metal/ligand ratios observed. The formation constants of divalent metal complexes with HPCIH were determined by potentiometric titrations and the values obtained were consistent with similar tridentate ligands and with the Irving-Williams order. ((C) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2003).

The folding of knotted proteins: insights from lattice simulations

We carry out systematic Monte Carlo simulations of Go lattice proteins to investigate and compare the folding processes of two model proteins whose native structures differ from each other due to the presence of a trefoil knot located near the terminus of one of the protein chains. We show that the folding time of the knotted fold is larger than that of the unknotted protein and that this difference in folding time is particularly striking in the temperature region below the optimal folding temperature. Both proteins display similar folding transition temperatures, which is indicative of similar thermal stabilities. By using the folding probability reaction coordinate as an estimator of folding progression we have found out that the formation of the knot is mainly a late folding event in our shallow knot system.

Biodegradation rate constants in different NF/UF fractions of cork processing wastewaters

Cork processing wastewater is an aqueous complex mixture of organic compounds that have been extracted from cork planks during the boiling process. These compounds, such as polysaccharides and polyphenols, have different biodegradability rates, which depend not only on the natureof the compound but also on the size of the compound. The aim of this study is to determine the biochemical oxygen demands (BOD) and biodegradationrate constants (k) for different cork wastewater fractions with different organic matter characteristics. These wastewater fractions were obtained using membrane separation processes, namely nanofiltration (NF) and ultrafiltration (UF). The nanofiltration and ultrafiltration membranes molecular weight cut-offs (MWCO) ranged from 0.125 to 91 kDa. The results obtained showed that the biodegradation rate constant for the cork processing wastewater was around 0.3 d(-1) and the k values for the permeates varied between 0.27-0.72 d(-1), being the lower values observed for permeates generated by the membranes with higher MWCO and the higher values observed for the permeates generated by the membranes with lower MWCO. These higher k values indicate that the biodegradable organic matter that is permeated by the membranes with tighter MWCO is more readily biodegradated.