A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

A note on extension variances in IR2

Matheron (1971) proposed an approximation of the extension variance in IR. We propose in this note an extension of this formula in IR2, based on a MacLaurin formula. Its application is shown in an example, the estimation of the maximum depressional storage of a soil surface.

Kinetic analysis of vascular marker distribution in perfused rat livers after regeneration following partial hepatectomy

Background/Aims: Liver clearance models are based on information (or assumptions) on solute distribution kinetics within the microvasculatory system, The aim was to study albumin distribution kinetics in regenerated livers and in livers of normal adult rats, Methods: A novel mathematical model was used to evaluate the distribution space and the transit time dispersion of albumin in livers following regeneration after a two-thirds hepatectomy compared to livers of normal adult rats. Outflow curves of albumin measured after bolus injection in single-pass perfused rat livers were analyzed by correcting for the influence of catheters and fitting a long-tailed function to the data. Results: The curves were well described by the proposed model. The distribution volume and the transit time dispersion of albumin observed in the partial hepatectomy group were not significantly different from livers of normal adult rats. Conclusions: These findings suggest that the distribution space and the transit time dispersion of albumin (CV2) is relatively constant irrespective of the presence of rapid and extensive repair. This invariance of CV2 implies, as a first approximation, a similar degree of intrasinusoidal mixing, The finding that a sum of two (instead of one) inverse Gaussian densities is an appropriate empirical function to describe the outflow curve of vascular indicators has consequences for an improved prediction of hepatic solute extraction.

Two-level atom in a squeezed vacuum with finite bandwidth

The resonance fluorescence of a two-level atom driven by a coherent laser field and damped by a finite bandwidth squeezed vacuum is analysed. We extend the Yeoman and Barnett technique to a non-zero detuning of the driving field from the atomic resonance and discuss the role of squeezing bandwidth and the detuning in the level shifts, widths and intensities of the spectral lines. The approach is valid for arbitrary values of the Rabi frequency and detuning but for the squeezing bandwidths larger than the natural linewidth in order to satisfy the Markoff approximation. The narrowing of the spectral lines is interpreted in terms of the quadrature-noise spectrum. We find that, depending on the Rabi frequency, detuning and the squeezing phase, different factors contribute to the line narrowing. For a strong resonant driving field there is no squeezing in the emitted field and the fluorescence spectrum exactly reveals the noise spectrum. In this case the narrowing of the spectral lines arises from the noise reduction in the input squeezed vacuum. For a weak or detuned driving field the fluorescence exhibits a large squeezing and, as a consequence, the spectral lines have narrowed linewidths. Moreover, the fluorescence spectrum can be asymmetric about the central frequency despite the symmetrical distribution of the noise. The asymmetry arises from the absorption of photons by the squeezed vacuum which reduces the spontaneous emission. For an appropriate choice of the detuning some of the spectral lines can vanish despite that there is no population trapping. Again this process can be interpreted as arising from the absorption of photons by the squeezed vacuum. When the absorption is large it may compensate the spontaneous emission resulting in the vanishing of the fluorescence lines.

Epidermal iontophoresis: I. Development of the ionic mobility-pore model

Purpose, An integrated ionic mobility-pore model for epidermal iontophoresis is developed from theoretical considerations using both the free volume and pore restriction forms of the model for a range of solute radii (r(j)) approaching the pore radii (r(p)) as well as approximation of the pore restriction form for r(j)/r(p) < 0.4. In this model, we defined the determinants for iontophoresis as solute size (defined by MV, MW or radius), solute mobility, solute shape, solute charge, the Debye layer thickness, total current applied, solute concentration, fraction ionized, presence of extraneous ions (defined by solvent conductivity), epidermal permselectivity, partitioning rates to account for interaction of unionized and ionized lipophilic solutes with the wall of the pore and electroosmosis. Methods, The ionic mobility-pore model was developed from theoretical considerations to include each of the determinants of iontophoretic transport. The model was then used to reexamine iontophoretic flux conductivity and iontophoretic flux-fraction ionized literature data on the determinants of iontophoretic flux. Results. The ionic mobility-pore model was found to be consistent with existing experimental data and determinants defining iontophoretic transport. However, the predicted effects of solute size on iontophoresis are more consistent with the pore-restriction than free volume form of the model. A reanalysis of iontophoretic flux-conductivity data confirmed the model's prediction that, in the absence of significant electroosmosis, the reciprocal of flux is linearly related to either donor or receptor solution conductivity. Significant interaction with the pore walls, as described by the model, accounted for the reported pH dependence of the iontophoretic transport for a range of ionizable solutes. Conclusions. The ionic mobility-pore iontophoretic model developed enables a range of determinants of iontophoresis to be described in a single unifying equation which recognises a range of determinants of iontophoretic flux.

Theory of multidimensional parametric band-gap simultons

Multidimensional spatiotemporal parametric simultons (simultaneous solitary waves) are possible in a nonlinear chi((2)) medium with a Bragg grating structure, where large effective dispersion occurs near two resonant band gaps for the carrier and second-harmonic field, respectively. The enhanced dispersion allows much reduced interaction lengths, as compared to bulk medium parametric simultons. The nonlinear parametric band-gap medium permits higher-dimensional stationary waves to form. In addition, solitons can occur with lower input powers than conventional nonlinear Schrodinger equation gap solitons. In this paper, the equations for electromagnetic propagation in a grating structure with a parametric nonlinearity are derived from Maxwell's equation using a coupled mode Hamiltonian analysis in one, two, and three spatial dimensions. Simultaneous solitary wave solutions are proved to exist by reducing the equations to the coupled equations describing a nonlinear parametric waveguide, using the effective-mass approximation (EMA). Exact one-dimensional numerical solutions in agreement with the EMA solutions are also given. Direct numerical simulations show that the solutions have similar types of stability properties to the bulk case, providing the carrier waves are tuned to the two Bragg resonances, and the pulses have a width in frequency space less than the band gap. In summary, these equations describe a physically accessible localized nonlinear wave that is stable in up to 3 + 1 dimensions. Possible applications include photonic logic and switching devices. [S1063-651X(98)06109-1].

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Modeling the optical constants of GaP, InP, and InAs

An extension of the Adachi model with the adjustable broadening function, instead of the Lorentzian one, is employed to model the optical constants of GaP, InP, and InAs. Adjustable broadening is modeled by replacing the damping constant with the frequency-dependent expression. The improved flexibility of the model enables achieving an excellent agreement with the experimental data. The relative rms errors obtained for the refractive index equal 1.2% for GaP, 1.0% for InP, and 1.6% for InAs. (C) 1999 American Institute of Physics. [S0021-8979(99)05807-7].

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Modeling the optical constants of AlxGa1-xAs alloys

The extension of Adachi's model with a Gaussian-like broadening function, in place of Lorentzian, is used to model the optical dielectric function of the alloy AlxGa1-xAs. Gaussian-like broadening is accomplished by replacing the damping constant in the Lorentzian line shape with a frequency dependent expression. In this way, the comparative simplicity of the analytic formulas of the model is preserved, while the accuracy becomes comparable to that of more intricate models, and/or models with significantly more parameters. The employed model accurately describes the optical dielectric function in the spectral range from 1.5 to 6.0 eV within the entire alloy composition range. The relative rms error obtained for the refractive index is below 2.2% for all compositions. (C) 1999 American Institute of Physics. [S0021-8979(99)00512-5].

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Wall material properties of yeast cells. Part II. Analysis

In the preceding paper (Part I) force-deformation data were measured with the compression experiment in conjunction with the initial radial stretch ratio and the initial wall-thickness to cell-radius ratio for baker's yeast (Saccharomyces cerevisiae). In this paper, these data have been analysed with the mechanical model of Smith et al. (Smith, Moxham & Middelberg (1998) Chemical Engineering Science, 53, 3913-3922) with the wall constitutive behaviour defined a priori as incompressible and linear-elastic. This analysis determined the mean Young's modulus ((E) over bar), mean maximum von Mises stress-at-failure (<(sigma)over bar>(VM,f)) and mean maximum von Mises strain-at failure (<(epsilon)over bar>(VM,f)) to be (E) over bar = 150 +/- 15 MPa, <(sigma)over bar>(VM,f) = 70 +/- 4 MPa and <(epsilon)over bar>(VM,f) = 0.75 +/- 0.08, respectively. The mean Young's modulus was not dependent (P greater than or equal to 0.05) on external osmotic pressure (0-0.8 MPa) nor compression rate (1.03-7.68 mu m/s) suggesting the incompressible linear-elastic relationship is representative of the actual cell-wall constitutive behaviour. Hydraulic conductivities were also determined and were comparable to other similar cell types (0-2.5 mu m/MPa s). The hydraulic conductivity distribution was not dependent on external osmotic pressure (0-0.8 MPa) nor compression rate (1.03-7.68 mu m/s) suggesting inclusion of cell-wall permeability in the mechanical model is justified. <(epsilon)over bar>(VM,f) was independent of cell diameter and to a first-approximation unaffected (P greater than or equal to 0.01) by external osmotic pressure and compression rate, thus providing a reasonable failure criterion. This criterion states that the cell-wall material will break when the strain exceeds <(epsilon)over bar>(VM,f) = 0.75 +/- 0.08. Variability in overall cell strength during compression was shown to be primarily due to biological variability in the maximum von Mises strain-at-failure. These data represent the first estimates of cell-wall material properties for yeast and the first fundamental analysis of cell-compression data. They are essential for describing cell-disruption at the fundamental level of fluid-cell interactions in general bioprocesses. They also provide valuable new measurements for yeast-cell physiologists. (C) 2000 Elsevier Science Ltd. All rights reserved.

Kinetics of propranolol uptake in muscle, skin, and fat of the isolated perfused rat hindlimb

The tissue distribution kinetics of a highly bound solute, propranolol, was investigated in a heterogeneous organ, the isolated perfused limb, using the impulse-response technique and destructive sampling. The propranolol concentration in muscle, skin, and fat as well as in outflow perfusate was measured up to 30 min after injection. The resulting data were analysed assuming (1) vascular, muscle, skin and fat compartments as well mixed (compartmental model) and (2) using a distributed-in-space model which accounts for the noninstantaneous intravascular mixing and tissue distribution processes but consists only of a vascular and extravascular phase (two-phase model). The compartmental model adequately described propranolol concentration-time data in the three tissue compartments and the outflow concentration-time curve (except of the early mixing phase). In contrast, the two-phase model better described the outflow concentration-time curve but is limited in accounting only for the distribution kinetics in the dominant tissue, the muscle. The two-phase model well described the time course of propranolol concentration in muscle tissue, with parameter estimates similar to those obtained with the compartmental model. The results suggest, first that the uptake kinetics of propranolol into skin and fat cannot be analysed on the basis of outflow data alone and, second that the assumption of well-mixed compartments is a valid approximation from a practical point of view las, e.g., in physiological based pharmacokinetic modelling). The steady-state distribution volumes of skin and fat were only 16 and 4%, respectively, of that of muscle tissue (16.7 ml), with higher partition coefficient in fat (6.36) than in skin (2.64) and muscle (2.79. (C) 2000 Elsevier Science B.V. All rights reserved.

Absorptive quantum measurements via coherently coupled quantum dots

We propose an absorptive measurement scheme via coupled quantum dots based on studies of the quantum dynamics of coherently coupled dots. The system is described through a Markov master equation that is related to a measurable quantity, the current. We analyse the measurement configuration and calculate the correlations and noise spectra beyond the adiabatic approximation.

Adsorption hysteresis and criticality in regular mesoporous materials

In the present work, various theories predicting the critical diameter for the absence of capillary condensation and hysteresis are applied to experimental adsorption isotherms of vapors on regular mesoporous materials. Among the various theories studied, the tensile strength approximation proposed by the authors was found to be the most successful. Reversibility of nitrogen adsorption at 77.4 K was studied on pure MCM-41 of various pore sizes, as well as mixtures of pure MCM-41 samples in a 1:1 ratio. The results of PSD and hysteresis on MCM-41 mixtures are close to that expected from studies of the pure materials. The estimates of hysteresis critical temperature and diameter of MCM-41, HMS, FSM and KIT materials are also provided.