Measuring entangled qutrits and their use for quantum bit commitment

We produce and holographically measure entangled qudits encoded in transverse spatial modes of single photons. With the novel use of a quantum state tomography method that only requires two-state superpositions, we achieve the most complete characterization of entangled qutrits to date. Ideally, entangled qutrits provide better security than qubits in quantum bit commitment: we model the sensitivity of this to mixture and show experimentally and theoretically that qutrits with even a small amount of decoherence cannot offer increased security over qubits.

Emergent anisotropy and flow alignment in viscous rock

A novel class of nonlinear, visco-elastic rheologies has recently been developed by MUHLHAUS et al. (2002a, b). The theory was originally developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of the layer surfaces or slip planes in the context of crystallographic slip is determined by the normal vector the so-called director of these surfaces. Here the model (MUHLHAUS et al., 2002a, b) is generalized to include thermal effects; it is shown that in 2-D steady states the director is given by the gradient of the flow potential. The model is applied to anisotropic simple shear where the directors are initially parallel to the shear direction. The relative effects of textural hardening and thermal softening are demonstrated. We then turn to natural convection and compare the time evolution and approximately steady states of isotropic and anisotropic convection for a Rayleigh number Ra=5.64x10(5) for aspect ratios of the experimental domain of 1 and 2, respectively. The isotropic case has a simple steady-state solution, whereas in the orthotropic convection model patterns evolve continuously in the core of the convection cell, which makes only a near-steady condition possible. This near-steady state condition shows well aligned boundary layers, and the number of convection cells which develop appears to be reduced in the orthotropic case. At the moderate Rayleigh numbers explored here we found only minor influences in the change from aspect ratio one to two in the model domain.

Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

Geopolymeric adsorbents from fly ash for dye removal from aqueous solution

Adsorbents from coal fly ash treated by a solid-state fusion method using NaOH were prepared. It was found that amorphous aluminosilicate, geopolymers would be formed. These fly ash-derived inorganic polymers were assessed as potential adsorbents for removal of some basic dyes, methylene blue and crystal violet, from aqueous solution. It was found that the adsorption capacity of the synthesised adsorbents depends on the preparation conditions such as NaOH:fly-ash ratio and fusion temperature with the optimal conditions being at 121 weight ratio of Na:fly-ash at 250-350 degrees C. The synthesised materials exhibit much higher adsorption capacity than fly ash itself and natural zeolite. The adsorption isotherm can be fitted by Langmuir and Freundlich models while the two-site Langmuir model producing the best results. It was also found that the fly ash derived geopolymeric adsorbents show higher adsorption capacity for crystal violet than methylene blue and the adsorption temperature influences the adsorption capacity. Kinetic studies show that the adsorption process follows the pseudo second-order kinetics. (c) 2006 Elsevier Inc. All rights reserved.

Hepatic structure-pharmacokinetic relationships: The hepatic disposition and metabolite kinetics of a homologous series of O-acyl derivatives of salicylic acid

1 The hepatic disposition and metabolite kinetics of a homologous series of O-acyl (acetyl, propionyl, butanoyl, pentanoyl, hexanoyl and octanoyl) esters of salicylic acid (C2SA, C3SA, C4SA, C5SA, C6SA and C8SA, respectively) was determined using a single-pass, in-sills rat liver preparation. 2 The hepatic venous outflow profiles for the parent esters and the generated metabolite, salicylic acid (SA) were analysed by HPLC. Non-parametric moments analysis was used to determine the area under the curve (AUC'), mean transit time (MTT) and normalized variance (CV2) for the parent esters and generated SA. 3 Pregenerated SA ([C-14]-salicylic acid) was injected into each liver with the parent ester to determine its distribution characteristics. 4 The overall recovery of ester plus metabolite was 89% of the ester dose injected and independent of the ester carbon number, suggesting that ester extraction was due to hepatic metabolism to salicylic acid. 5 The metabolite AUC' value increased directly with the lipophilicity of the parent ester (from 0.12 for C2SA to 0.95 for C8SA). By contrast, the parent AUC' decreased with the lipophilicity (from 0.85 for C2SA to zero for C8SA). The metabolite MTT value also showed a trend to increase with the lipophilicity of the parent ester (from 15.72 s for C3SA to 61.97 s for C8SA). However, the parent MTT value shows no significant change across the series. 6 The two-compartment dispersion model was used to derive the kinetic parameters for parent ester, pregenerated SA and generated SA. Consequently, these parameters were used to estimate the values of AUG', MITT and CV2 for the parent ester and metabolite. The moments values obtained using the two-compartment dispersion model show similar trends to the corresponding moments values obtained from the outflow profiles using a non-parametric approach. 7 The more lipophilic aspirin analogues are more confined to the portal circulation after oral administration than aspirin due to their more extensive hepatic elimination avoiding systemic prostacyclin inhibition. Given that aspirin's selectivity as an anti-thrombotic agent has been postulated to be due to selective anti-platelet effects in the portal circulation, the more lipophilic and highly extracted analogues are potentially more selective anti-thrombotic agents than aspirin.

Hepatic disposition and metabolite kinetics of a homologous series of diflunisal esters

The hepatic disposition and metabolite kinetics of a homologous series of diflunisal O-acyl esters (acetyl, butanoyl, pentanoyl, anti hexanoyl) were determined using a single-pass perfused in situ rat liver preparation. The experiments were conducted using 2% BSA Krebs-Henseleit buffer (pH 7.4), and perfusions were performed at 30 mL/min in each liver. O-Acyl esters of diflunisal and pregenerated diflunisal were injected separately into the portal vein. The venous outflow samples containing the esters and metabolite diflunisal were analyzed by high performance liquid chromatography (HPLC). The normalized outflow concentration-time profiles for each parent ester and the formed metabolite, diflunisal, were analyzed using statistical moments analysis and the two-compartment dispersion model. Data (presented as mean +/- standard error for triplicate experiments) was compared using ANOVA repeated measures, significance level P < 0.05. The hepatic availability (AUC'), the fraction of the injected dose recovered in the outflowing perfusate, for O-acetyldiflunisal (C2D = 0.21 +/- 0.03) was significantly lower than the other esters (0.34-0.38). However, R-N/f(u), the removal efficiency number R-N divided by the unbound fraction in perfusate f(u), which represents the removal efficiency of unbound ester by the liver, was significantly higher for the most lipophilic ester (O-hexanoyldiflunisal, C6D = 16.50 +/- 0.22) compared to the other members of the series (9.57 to 11.17). The most lipophilic ester, C6D, had the largest permeability surface area (PS) product (94.52 +/- 38.20 mt min-l g-l liver) and tissue distribution value VT (35.62 +/- 11.33 mL g(-1) liver) in this series. The MTT of these O-acyl esters of diflunisal were not significantly different from one another. However, the metabolite diflunisal MTTs tended to increase with the increase in the parent ester lipophilicity (11.41 +/- 2.19 s for C2D to 38.63 +/- 9.81 s for C6D). The two-compartment dispersion model equations adequately described the outflow profiles for the parent esters and the metabolite diflunisal formed from the O-acyl esters of diflunisal in the liver.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

Unitary R-matrices for topological quantum computing

The main problem with current approaches to quantum computing is the difficulty of establishing and maintaining entanglement. A Topological Quantum Computer (TQC) aims to overcome this by using different physical processes that are topological in nature and which are less susceptible to disturbance by the environment. In a (2+1)-dimensional system, pseudoparticles called anyons have statistics that fall somewhere between bosons and fermions. The exchange of two anyons, an effect called braiding from knot theory, can occur in two different ways. The quantum states corresponding to the two elementary braids constitute a two-state system allowing the definition of a computational basis. Quantum gates can be built up from patterns of braids and for quantum computing it is essential that the operator describing the braiding-the R-matrix-be described by a unitary operator. The physics of anyonic systems is governed by quantum groups, in particular the quasi-triangular Hopf algebras obtained from finite groups by the application of the Drinfeld quantum double construction. Their representation theory has been described in detail by Gould and Tsohantjis, and in this review article we relate the work of Gould to TQC schemes, particularly that of Kauffman.

Applying linear time-varying constraints to econometric models: With an application to demand systems

When linear equality constraints are invariant through time they can be incorporated into estimation by restricted least squares. If, however, the constraints are time-varying, this standard methodology cannot be applied. In this paper we show how to incorporate linear time-varying constraints into the estimation of econometric models. The method involves the augmentation of the observation equation of a state-space model prior to estimation by the Kalman filter. Numerical optimisation routines are used for the estimation. A simple example drawn from demand analysis is used to illustrate the method and its application.

Impulse-response studies on tracer doses of [C-14]lignocaine and its multiple metabolites in the perfused rat liver

The outflow-concentration-time profiles for lignocaine (lidocaine) and its metabolites have been measured after bolus impulse administration of [C-14]lignocaine into the perfused rat liver. Livers from female Sprague-Dawley rats were perfused in a once-through fashion with red-blood-cell-free Krebs-Henseleit buffer containing 0 or 2% bovine serum albumin. Perfusate flow rates of 20 and 30 mL min(-1) were used and both normal and retrograde flow directions were employed. Significant amounts of metabolite were detected in the effluent perfusate soon after lignocaine injection. The early appearance of metabolite contributed to bimodal outflow profiles observed for total C-14 radioactivity. The lignocaine outflow profiles were well characterized by the two-compartment dispersion model, with efflux rate << influx rate. The profiles for lignocaine metabolites were also characterized in terms of a simplified two-compartment dispersion model. Lignocaine was found to be extensively metabolized under the experimental conditions with the hepatic availability ranging between 0.09 and 0.18. Generally lignocaine and metabolite availability showed no significant change with alterations in perfusate flow rate from 20 to 30 mt min(-1) or protein content from 0 to 2%. A significant increase in lignocaine availability occurred when 1200 mu M unlabelled lignocaine was added to the perfusate. Solute mean transit times generally decreased with increasing flow rate and with increasing perfusate protein content. The results confirm that lignocaine pharmacokinetics in the liver closely follow the predictions of the well-stirred model. The increase in lignocaine availability when 1200 mu M unlabelled lignocaine was added to the perfusate is consistent with saturation of the hydroxylation metabolic pathways of lignocaine metabolism.

Structure-hepatic disposition relationships for cationic drugs in isolated perfused rat livers: Transmembrane exchange and cytoplasmic binding process

This work studied the structure-hepatic disposition relationships for cationic drugs of varying lipophilicity using a single-pass, in situ rat liver preparation. The lipophilicity among the cationic drugs studied in this work is in the following order: diltiazem. propranolol. labetalol. prazosin. antipyrine. atenolol. Parameters characterizing the hepatic distribution and elimination kinetics of the drugs were estimated using the multiple indicator dilution method. The kinetic model used to describe drug transport (the two-phase stochastic model) integrated cytoplasmic binding kinetics and belongs to the class of barrier-limited and space-distributed liver models. Hepatic extraction ratio (E) (0.30-0.92) increased with lipophilicity. The intracellular binding rate constant (k(on)) and the equilibrium amount ratios characterizing the slowly and rapidly equilibrating binding sites (K-S and K-R) increase with the lipophilicity of drug (k(on) : 0.05-0.35 s(-1); K-S : 0.61-16.67; K-R : 0.36-0.95), whereas the intracellular unbinding rate constant (k(off)) decreases with the lipophilicity of drug (0.081-0.021 s(-1)). The partition ratio of influx (k(in)) and efflux rate constant (k(out)), k(in)/k(out), increases with increasing pK(a) value of the drug [from 1.72 for antipyrine (pK(a) = 1.45) to 9.76 for propranolol (pK(a) = 9.45)], the differences in k(in/kout) for the different drugs mainly arising from ion trapping in the mitochondria and lysosomes. The value of intrinsic elimination clearance (CLint), permeation clearance (CLpT), and permeability-surface area product (PS) all increase with the lipophilicity of drug [CLint (ml . min(-1) . g(-1) of liver): 10.08-67.41; CLpT (ml . min(-1) . g(-1) of liver): 10.80-5.35; PS (ml . min(-1) . g(-1) of liver): 14.59-90.54]. It is concluded that cationic drug kinetics in the liver can be modeled using models that integrate the presence of cytoplasmic binding, a hepatocyte barrier, and a vascular transit density function.

Pharmacokinetics and pharmacodynamics of melphalan in isolated limb infusion for recurrent localized limb malignancy

Isolated limb infusion (ILI) is an attractive, less complex alternative to Isolated limb perfusion (ILP). It has a lower morbidity in treating localized recurrences and in transit metastases of the limb for tumours such as melanoma, Merkel cell tumour and Kaposi's sarcoma, allowing administration of high concentrations of cytotoxic agent to the affected limb under hypoxic conditions. Melphalan is the preferred cytotoxic agent for the treatment of melanoma by ILP or ILI. We report pharmacokinetic data from 12 patients treated by ILI for tumours of the limb in Brisbane. The kinetics of drug distribution in the limb was calculated using a two-compartment vascular model, where both tissue and infusate act as well-stirred compartments. Analysis of melphalan concentrations in the perfusate during ILI showed good agreement between the values measured and the concentrations predicted by the model. Recirculation and wash-out flow rates, tissue concentrations and the permeability surface area product (PS) were calculated. Correlations between the PS value and the drug concentrations In the perfusate and tissue were supported by the results. These data contribute to a better understanding of the distribution of melphalan during ILI in the limb, and offer the opportunity to optimize the drug regimen for patients undergoing ILI. (C) 2001 Lippincott Williams & Wilkins.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

A framework for extreme-event control in wastewater treatment

In this paper an approach to extreme event control in wastewater treatment plant operation by use of automatic supervisory control is discussed. The framework presented is based on the fact that different operational conditions manifest themselves as clusters in a multivariate measurement space. These clusters are identified and linked to specific and corresponding events by use of principal component analysis and fuzzy c-means clustering. A reduced system model is assigned to each type of extreme event and used to calculate appropriate local controller set points. In earlier work we have shown that this approach is applicable to wastewater treatment control using look-up tables to determine current set points. In this work we focus on the automatic determination of appropriate set points by use of steady state and dynamic predictions. The performance of a relatively simple steady-state supervisory controller is compared with that of a model predictive supervisory controller. Also, a look-up table approach is included in the comparison, as it provides a simple and robust alternative to the steady-state and model predictive controllers, The methodology is illustrated in a simulation study.

Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.