Exact solution for the Bariev model with boundary fields

The Bariev model with open boundary conditions is introduced and analysed in detail in the framework of the Quantum Inverse Scattering Method. Two classes of independent boundary reflecting K-matrices leading to four different types of boundary fields are obtained by solving the reflection equations. The models are exactly solved by means of the algebraic nested Bethe ansatz method and the four sets or Bethe ansatz equations as well as their corresponding energy expressions are derived. (C) 2001 Elsevier Science B.V. All rights reserved.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

A semi-analytical solution for vibration of rectangular plates with abrupt thickness variation

A semi-analytical analysis of free vibration of plates with cross-sectional discontinuities due to abrupt changes in thickness is presented. A basic square element divided into suitable subdomains dependent upon the positions of these abrupt changes is used as the basic building element, Admissible functions that satisfy the essential or geometric boundary conditions are used to define the transverse deflection of each subdomain. Continuities in the displacement, slope, moment and higher derivatives between adjacent subdomains are enforced at the interconnecting edges. The resulting global energy functional from the proper assembly of the coupled strain and kinetic energy contributions of each subdomain is then minimized via the Ritz procedure to extract the frequencies and mode shapes. Contour plots of a range of new mode shapes are presented for the enhancement of understanding the dynamic behavior of this class of plates, (C) 2001 Elsevier Science Ltd, All rights reserved.

Predation of brown tree snakes (Boiga irregularis) in Australia

We conducted a literature review to address the potential for using a native, vertebrate predator of brown tree snakes (Boiga irregularis) as a biological control method on Guam. Both actual and potential predators were included in our review. We located two actual predators (red-bellied black snakes (Pseudechis porphyriacus) and cane toads (Bufo marinus)) and 55 potential predators of brown tree snakes. However, none of the native predators of brown tree snakes appear likely candidates as a biological control method on Guam due to their lack of selectivity in their feeding habits and unknown aspects of their natural history. (C) 2002 Elsevier Science Ltd. All rights reserved.

A high definition, finite difference time domain method

A high definition, finite difference time domain (HD-FDTD) method is presented in this paper. This new method allows the FDTD method to be efficiently applied over a very large frequency range including low frequencies, which are problematic for conventional FDTD methods. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and has been verified against analytical solutions within the frequency range 50 Hz-1 GHz. As an example of the lower frequency range, the method has been applied to the problem of induced eddy currents in the human body resulting from the pulsed magnetic field gradients of an MRI system. The new method only requires approximately 0.3% of the source period to obtain an accurate solution. (C) 2003 Elsevier Science Inc. All rights reserved.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.

Analysis of the free vibration of rectangular plates with central cut-outs using the discrete Ritz method

We present an analysis of the free vibration of plates with internal discontinuities due to central cut-outs. A numerical formulation for a basic L-shaped element which is divided into appropriate sub-domains that are dependent upon the location of the cut-out is used as the basic building element. Trial functions formed to satisfy certain boundary conditions are employed to define the transverse deflection of each sub-domain. Mathematical treatments in terms of the continuities in displacement, slope, moment, and higher derivatives between the adjacent sub-domains are enforced at the interconnecting edges. The energy functional results, from the proper assembly of the coupled strain and kinetic energy contributions of each sub-domain, are minimized via the Ritz procedure to extract the vibration frequencies and. mode shapes of the plates. The procedures are demonstrated by considering plates with central cut-outs that are subjected to two types of boundary conditions. (C) 2003 Elsevier Ltd. All rights reserved.

A convergent solution-phase synthesis of the macrocycle Ac-Phe-[Orn-Pro-D-Cha-Trp-Arg], a potent new antiinflammatory drug

Relatively few cyclic peptides have reached the pharmaceutical marketplace during the past decade, most produced through fermentation rather than made synthetically. Generally, this class of compounds is synthesized for research purposes on milligram scales by solid-phase methods, but if the potential of macrocyclic peptidomimetics is to be realized, low-cost larger scale solution-phase syntheses need to be devised and optimized to provide sufficient quantities for preclinical, clinical, and commercial uses. Here, we describe a cheap, medium-scale, solution-phase synthesis of the first reported highly potent, selective, and orally active antagonist of the human C5a receptor. This compound, Ac-Phe[Orn-Pro-D-Cha-Trp-Arg], known as 3D53, is a macrocyclic peptidomimetic of the human plasma protein C5a and displays excellent antiinflammatory activity in numerous animal models of human disease. In a convergent approach, two tripeptide fragments Ac-Phe-Orn-(Boc)-Pro-OH and H-D-Cha-Trp(For)-Arg-OEt were first prepared by high-yielding solution-phase couplings using a mixed anhydride method before coupling them to give a linear hexapeptide which, after deprotection, was obtained in 38% overall yield from the commercially available amino acids. Cyclization in solution using BOP reagent gave the antagonist in 33% yield (13% overall) after HPLC purification. Significant features of the synthesis were that the Arg side chain was left unprotected throughout, the component Boe-D-Cha-OH was obtained very efficiently via hydrogenation Of D-Phe with PtO2 in TFA/water, the tripeptides were coupled at the Pro-Cha junction to minimize racemization via the oxazolone pathway, and the entire synthesis was carried out without purification of any intermediates. The target cyclic product was purified (>97%) by reversed-phase HPLC. This convergent synthesis with minimal use of protecting groups allowed batches of 50100 g to be prepared efficiently in high yield using standard laboratory equipment. This type of procedure should be useful for making even larger quantities of this and other macrocyclic peptidomimetic drugs.

Exact solution of the XXZ Gaudin model with generic open boundaries

The XXZ Gaudin model with generic integrable boundaries specified by generic non-diagonal K-matrices is studied. The commuting families of Gaudin operators are diagonalized by the algebraic Bethe ansatz method. The eigenvalues and the corresponding Bethe ansatz equations are obtained. (C) 2004 Elsevier B.V. All rights reserved.

Semi-analytical solution for nonlinear vibration of laminated FGM plates with geometric imperfections

This paper investigates the nonlinear vibration of imperfect shear deformable laminated rectangular plates comprising a homogeneous substrate and two layers of functionally graded materials (FGMs). A theoretical formulation based on Reddy's higher-order shear deformation plate theory is presented in terms of deflection, mid-plane rotations, and the stress function. A semi-analytical method, which makes use of the one-dimensional differential quadrature method, the Galerkin technique, and an iteration process, is used to obtain the vibration frequencies for plates with various boundary conditions. Material properties are assumed to be temperature-dependent. Special attention is given to the effects of sine type imperfection, localized imperfection, and global imperfection on linear and nonlinear vibration behavior. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with graded silicon nitride/stainless steel layers. It is shown that the vibration frequencies are very much dependent on the vibration amplitude and the imperfection mode and its magnitude. While most of the imperfect laminated plates show the well-known hard-spring vibration, those with free edges can display soft-spring vibration behavior at certain imperfection levels. The influences of material composition, temperature-dependence of material properties and side-to-thickness ratio are also discussed. (C) 2004 Elsevier Ltd. All rights reserved.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Pseudo-ternary phase diagrams of aqueous mixtures of Quil A, cholesterol and phospholipid prepared by the lipid-film hydration method

Pseudo-ternary phase diagrams of the polar lipids Quil A, cholesterol (Chol) and phosphatidylcholine (PC) in aqueous mixtures prepared by the lipid film hydration method (where dried lipid film of phospholipids and cholesterol are hydrated by an aqueous solution of Quil A) were investigated in terms of the types of particulate structures formed therein. Negative staining transmission electron microscopy and polarized light microscopy were used to characterize the colloidal and coarse dispersed particles present in the systems. Pseudo-ternary phase diagrams were established for lipid mixtures hydrated in water and in Tris buffer (pH 7.4). The effect of equilibration time was also studied with respect to systems hydrated in water where the samples were stored for 2 months at 4degreesC. Depending on the mass ratio of Quil A, Chol and PC in the systems, various colloidal particles including ISCOM matrices, liposomes, ring-like micelles and worm-like micelles were observed. Other colloidal particles were also observed as minor structures in the presence of these predominant colloids including helices, layered structures and lamellae (hexagonal pattern of ring-like micelles). In terms of the conditions which appeared to promote the formation of ISCOM matrices, the area of the phase diagrams associated with systems containing these structures increased in the order: hydrated in water/short equilibration period < hydrated in buffer/short equilibration period < hydrated in water/prolonged equilibration period. ISCOM matrices appeared to form over time from samples, which initially contained a high concentration of ring-like micelles suggesting that these colloidal structures may be precursors to ISCOM matrix formation. Helices were also frequently found in samples containing ISCOM matrices as a minor colloidal structure. Equilibration time and presence of buffer salts also promoted the formation of liposomes in systems not containing Quil A. These parameters however, did not appear to significantly affect the occurrence and predominance of other structures present in the pseudo-binary systems containing Quil A. Pseudo-ternary phase diagrams of PC, Chol and Quil A are important to identify combinations which will produce different colloidal structures, particularly ISCOM matrices, by the method of lipid film hydration. Colloidal structures comprising these three components are readily prepared by hydration of dried lipid films and may have application in vaccine delivery where the functionality of ISCOMs has clearly been demonstrated. (C) 2003 Elsevier B.V. All rights reserved.