Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Size and power properties of the AH test of evolutionary change

We investigate the size and power properties of the AH test of evolutionary change. This involves examining whether the size results are sensitive to both the number of individual frequencies estimated and the spectral shape adopted under the null hypothesis. The power tests examine whether the test has good power to detect shifts in both spectral position (variance) and spectral shape (autocovariance structure).

Safety and efficacy of low anterior en bloc resection as part of cytoreductive surgery for patients with ovarian cancer

Objective. To examine the feasibility and safety of a low anterior resection of the rectosigmoid plus adjacent pelvic tumour as part of primary cytoreduction for ovarian cancer. Methods. This study included 65 consecutive patients with primary ovarian cancer who had debulking surgery from 1996 through 2000. All patients underwent an en bloc resection of ovarian cancer and a rectosigmoid resection followed by an end-to-end anastomosis. Parameters for safety and efficacy were considered as primary statistical endpoints for the aim of this analysis. Results. Postoperative residual tumour was nil, 1 cm in 14, 34, and 14 patients, respectively. The median postoperative hospital stay was 11 days (range, 6 to 50 days). Intraoperative complications included an injury to the urinary bladder in one patient. Postoperative complications included wound complications (n=14, 21.5%), septicemia (n=9, 13.8%), cardiac complications (n=7, 10.8%), thromboembolic complications (n=5, 7.7%) ileus (n=2, 3.1%) anastomotic leak (n=2, 3.1%) and fistula (n=1, 1.5%). Reasons for a reoperation during the same admission included repair of an anastomotic leak (n=1), postoperative hemorrhage (n=1), and wound debridement (n=1). Wound complications, septicemia, and anastomotic leak formation were more frequent in patients who had a serum albumin level of less than or equal to 30 g/L preoperatively. There was one surgically related mortality in a patient who died from a cerebral vascular accident 2 days postoperatively. Conclusions. An en bloc resection as part of primary cytoreductive surgery for ovarian cancer is effective and its morbidity is acceptably low. (C) 2001 Academic Press.

Quantification of sulfur and sulfur-containing compounds in wastewaters by means of a combination of liquid chromatographic methods

Low-micromolar concentrations of sulfite, thiosulfate and sulfide, present in synthetic wastewater or anaerobic digester effluent, were quantified by means of derivatization with monobromobimane, followed by HPLC separation with fluorescence detection. The concentration of elemental sulfur was determined, after its extraction with chloroform from the derivatized sample, by HPLC with UV detection. Recoveries of sulfide (both matrices), and of thiosulfate and sulfite (synthetic wastewater) were between 98 and 103%. The in-run RSDs on separate derivatizations were 13 and 19% for sulfite (two tests), between 1.5 and 6.6% for thiosulfate (two tests) and between 4.1 and 7.7% for sulfide (three tests). Response factors for derivatives of sulfide and thiosulfate, but not sulfite, were steady over a 13-month period during which 730 samples were analysed. Dithionate and tetrathionate did not seem to be detectable with this method. The distinctness of the elemental sulfur and the derivatizing-agent peaks was improved considerably by detecting elution at 297 instead of 263 nm. (C) 2002 Elsevier Science B.V. All rights reserved.

Integrability and exact spectrum of a pairing model for nucleons

A pairing model for nucleons, introduced by Richardson in 1966, which describes proton-neutron pairing as well as proton-proton and neutron-neutron pairing, is re-examined in the context of the quantum inverse scattering method. Specifically, this shows that the model is integrable by enabling the explicit construction of the conserved operators. We determine the eigenvalues of these operators in terms of the Bethe ansatz, which in turn leads to an expression for the energy eigenvalues of the Hamiltonian.

An analysis of occupational outcomes for Indigenous and Asian employees in Australia

This paper examines occupational performance in Australia across three racial groups in Australia: Indigenous Australians; Asian people, defined as all those whose language spoken at home was either Chinese, Vietnamese or other forms of a South-east or East Asian language; and white people, defined as the residual category. The paper has as its starting point, observed differences in occupational attainment among the three groups in Australia and sets out to account for these observed differences on the basis of both race and non-racial attributes such as, age, education and area of residence.

Carbene and nitrene rearrangements: A theoretical study of cyclic allenes and carbenes, carbodiimides, and azirines

B3LYP/6-31G(d) calculations of structures, energies, and infrared spectra of several rearrangement products of (hetero)aromatic nitrenes and carbenes are reported. 3-Isoquinolylnitrene 36 ring closes to the azirine 37 prior to ring expansion to the potentially stable but unobserved seven-membered-ring carbodiimide 38 and diazacycloheptatrienylidene C-s-39S. A new, stable cycloheptatrienylidene, C-s-19S, is located on the naphthylcarbene energy surface. 4-Quinolylnitrene undergoes reaction via the azirine 50 in solution, but ring expansion to the stable seven-membered-ring ketenimine 47 under Ar matrix photolysis conditions. There is excellent agreement between calculated infrared spectra of 1,5-diazacyclohepta-1,2,4,6-tetraene 54 (obtained by photolysis of 4-pyridyl azide), 1-azacyclohepta-1,2,4,6-tetraene 5, 1-azacyclohepta-1,3,5,6-tetraene 55, and 1-azacyclohepta-1,3,4,6-tetraene 56 and the available experimental data.

Automatic construction of explicit R matrices for the one-parameter families of irreducible typical highest weight (0(m)vertical bar alpha(n)) representations of U-q[gl(m vertical bar n)]

We detail the automatic construction of R matrices corresponding to (the tensor products of) the (O-m\alpha(n)) families of highest-weight representations of the quantum superalgebras Uq[gl(m\n)]. These representations are irreducible, contain a free complex parameter a, and are 2(mn)-dimensional. Our R matrices are actually (sparse) rank 4 tensors, containing a total of 2(4mn) components, each of which is in general an algebraic expression in the two complex variables q and a. Although the constructions are straightforward, we describe them in full here, to fill a perceived gap in the literature. As the algorithms are generally impracticable for manual calculation, we have implemented the entire process in MATHEMATICA; illustrating our results with U-q [gl(3\1)]. (C) 2002 Published by Elsevier Science B.V.

From molecular complexes to zwitterions and final products. Reactions between C3O2 and amines

The formation of molecular complexes (prereactive intermediates) between C3O2 and amines (ammonia, dimethylamine, trimethylamine, and 4-(dimethylamino)pyridine) as well as the subsequent transformation of the complexes into C3O2-amine zwitterions in cryogenic matrixes (ca. 40 K) has been observed. In the case of dimethylamine, the formation of tetramethylmalonamide has also been documented. Calculations using density functional theory (B3LYP/6-31G(2d, p)) are used to assign all above species and are in excellent agreement with the IR spectra.

Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

The residual effects of cooperative-learning experiences: A two-year follow-up

The author investigated how training in small-group and interpersonal behaviors affected children's behavior and interactions as they worked in small groups 2 years later. The authors assigned 52 fifth graders, who had been trained 2 years previously in cooperative group behaviors, to the trained condition and 36 fifth graders, who had not previously been trained, to the untrained condition. Both were reconstituted from the pool of students who had participated previously in group activities. The results showed a residual training effect, with the children in the trained groups being more cooperative and helpful than their untrained peers.

Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Comparison of an ELISA and an LC/MS/MS method for measuring tacrolimus concentrations and making dosage decisions in transplant recipients

This study compared an enzyme-linked immunosorbent assay (ELISA) to a liquid chromatography-tandem mass spectrometry (LC/MS/MS) technique for measurement of tacrolimus concentrations in adult kidney and liver transplant recipients, and investigated how assay choice influenced pharmacokinetic parameter estimates and drug dosage decisions. Tacrolimus concentrations measured by both ELISA and LC/MS/MS from 29 kidney (n = 98 samples) and 27 liver (n = 97 samples) transplant recipients were used to evaluate the performance of these methods in the clinical setting. Tacrolimus concentrations measured by the two techniques were compared via regression analysis. Population pharmacokinetic models were developed independently using ELISA and LC/MS/MS data from 76 kidney recipients. Derived kinetic parameters were used to formulate typical dosing regimens for concentration targeting. Dosage recommendations for the two assays were compared. The relation between LC/MS/MS and ELISA measurements was best described by the regression equation ELISA = 1.02 . (LC/MS/MS) + 0.14 in kidney recipients, and ELISA = 1.12 . (LC/MS/MS) - 0.87 in liver recipients. ELISA displayed less accuracy than LC/MS/MS at lower tacrolimus concentrations. Population pharmacokinetic models based on ELISA and LC/MS/MS data were similar with residual random errors of 4.1 ng/mL and 3.7 ng/mL, respectively. Assay choice gave rise to dosage prediction differences ranging from 0% to 30%. ELISA measurements of tacrolimus are not automatically interchangeable with LC/MS/MS values. Assay differences were greatest in adult liver recipients, probably reflecting periods of liver dysfunction and impaired biliary secretion of metabolites. While the majority of data collected in this study suggested assay differences in adult kidney recipients were minimal, findings of ELISA dosage underpredictions of up to 25% in the long term must be investigated further.

Levene tests of homogeneity of variance for general block and treatment designs

This article develops a weighted least squares version of Levene's test of homogeneity of variance for a general design, available both for univariate and multivariate situations. When the design is balanced, the univariate and two common multivariate test statistics turn out to be proportional to the corresponding ordinary least squares test statistics obtained from an analysis of variance of the absolute values of the standardized mean-based residuals from the original analysis of the data. The constant of proportionality is simply a design-dependent multiplier (which does not necessarily tend to unity). Explicit results are presented for randomized block and Latin square designs and are illustrated for factorial treatment designs and split-plot experiments. The distribution of the univariate test statistic is close to a standard F-distribution, although it can be slightly underdispersed. For a complex design, the test assesses homogeneity of variance across blocks, treatments, or treatment factors and offers an objective interpretation of residual plot.

Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.