The composite Euler method for stiff stochastic differential equations

In this paper we present the composite Euler method for the strong solution of stochastic differential equations driven by d-dimensional Wiener processes. This method is a combination of the semi-implicit Euler method and the implicit Euler method. At each step either the semi-implicit Euler method or the implicit Euler method is used in order to obtain better stability properties. We give criteria for selecting the semi-implicit Euler method or the implicit Euler method. For the linear test equation, the convergence properties of the composite Euler method depend on the criteria for selecting the methods. Numerical results suggest that the convergence properties of the composite Euler method applied to nonlinear SDEs is the same as those applied to linear equations. The stability properties of the composite Euler method are shown to be far superior to those of the Euler methods, and numerical results show that the composite Euler method is a very promising method. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

Robust parameter settings of evolutionary algorithms for the optimisation of agricultural systems models

Numerical optimisation methods are being more commonly applied to agricultural systems models, to identify the most profitable management strategies. The available optimisation algorithms are reviewed and compared, with literature and our studies identifying evolutionary algorithms (including genetic algorithms) as superior in this regard to simulated annealing, tabu search, hill-climbing, and direct-search methods. Results of a complex beef property optimisation, using a real-value genetic algorithm, are presented. The relative contributions of the range of operational options and parameters of this method are discussed, and general recommendations listed to assist practitioners applying evolutionary algorithms to the solution of agricultural systems. (C) 2001 Elsevier Science Ltd. All rights reserved.

Separation of mercury (II), lead (II), and silver (I) from aqueous solution using an aminopolycarboxylate silica

A quartz crystal microbalance modified by the attachment of silica particles derivatized with the aminopolycarboxylate ligand N-[(3-trimethoxysilyl)propyl]ethylenediamine-N,N',N'-triacetic acid has been employed to assess conditions under which mercury (II), lead (II), and silver (I) nitrates may be separated in aqueous solution. The separation protocol, which involved removal of Hg(II), as [HgI4](2-), and Pb(II) with H+ was successfully applied to a batchwise separation of the 3 metal ions.

Diverse solid-state and solution structures within a series of hexaamine dicopper(II) complexes

Mono- and dicopper(II) complexes of a series of potentially bridging hexaamine ligands have been prepared and characterized in the solid state by X-ray crystallography. The crystal structures of the following Cu-II complexes are reported: [Cu(HL3)](ClO4)(3), C11H31Cl3CuN6O12, monoclinic, P2(1)/n, a = 8.294(2) Angstrom, b = 18.364(3) Angstrom, c = 15.674(3) Angstrom, beta = 94.73(2)degrees, Z = 4; {[Cu-2(L-4)(CO3)](2)}(ClO4)(4). 4H(2)O, C40H100Cl4Cu4N12O26, triclinic, P (1) over bar, a = 9.4888(8) Angstrom, b=13.353(1) Angstrom,. c = 15.329(1) Angstrom, alpha = 111.250(7)degrees, beta = 90.068(8)degrees, gamma = 105.081(8)degrees, Z=1; [Cu-2(L-5)(OH2)(2)](ClO4)(4), C(13)H(36)Cl(4)Cu(2)Z(6)O(18), monoclinic, P2(1)/c, a = 7.225(2) Angstrom. b = 8.5555(5) Angstrom, c = 23.134(8) Angstrom, beta = 92.37(1)degrees, Z = 2; [Cu-2(L-6)(OH2)(2)](ClO4)(4). 3H(2)O, C14H44Cl4Cu2N6O21, monoclinic, P2(1)/a, a = 15.204(5) Angstrom, b = 7.6810(7) Angstrom, c = 29.370(1) Angstrom, beta = 100.42(2)degrees, Z = 4. Solution spectroscopic properties of the bimetallic complexes indicate that significant conformational changes occur upon dissolution, and this has been probed with EPR spectroscopy and molecular mechanics calculations.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Analysis of acousto-ultrasonic characteristics for contact-type transducers coupled to composite laminated plates

The acousto-ultrasonic (AU) input-output characteristics for contact-type transmitting and receiving transducers coupled to composite laminated plates are considered in this paper. Combining a multiple integral transform method, an ordinary discrete layer theory for the laminates and some simplifying assumptions for the electro-mechanical transduction behaviour of the transducers, an analytical solution is developed which can deal with all the wave processes involved in the AU measurement system, i.e, wave generation, wave propagation and wave reception. The spectral response of the normal contact pressure sensed by the receiving transducer due to an arbitrary input pulse excited by the transmitting transducer is obtained. To validate the new analytical-numerical spectral technique in the low-frequency regime, the results are compared with Mindlin plate theory solutions. Based on the analytical results, numerical calculations are carried out to investigate the influence of various external parameters such as frequency content of the input pulse, transmitter/receiver spacing and transducer aperture on the output of the measurement system. The results show that the presented analytical-numerical procedure is an effective tool for understanding the input-output characteristics of the AU technique for laminated plates. (C) 2001 Elsevier Science Ltd. All rights reserved.

Solution structures by H-1 NMR of the novel cyclic trypsin inhibitor SFTI-1 from sunflower seeds and an acyclic permutant

SFTI-1 is a recently discovered cyclic peptide trypsin inhibitor from sunflower seeds comprising 14 amino acid residues. It is the most potent known Bowman-Birk inhibitor and the only naturally occurring cyclic one. The solution structure of SFTI-1 has been determined by H-1-NMR spectroscopy and compared with a synthetic acyclic permutant. The solution structures of both are remarkably similar. The lowest energy structures from each family of 20 structures of cyclic and acyclic SFTI-1 have an rmsd over the backbone and heavy atoms of 0.29 Angstrom and 0.66 Angstrom, respectively. The structures consist of two short antiparallel beta -strands joined by an extended loop containing the active site at one end. Cyclic SFTI-1 also has a hairpin turn completing the cycle. Both molecules contain particularly stable arrangements of cross-linking hydrogen bonds between the beta -strands and a single disulfide bridge, making them rigid and well defined in solution. These stable arrangements allow both the cyclic and acyclic variants of SFTI-1 to inhibit trypsin with very high potencies (0.5 nM and 12.1 nM, respectively). The cyclic nature of SFTI-1 appears to have evolved to provide higher trypsin inhibition as well as higher stability. The solution structures are similar to the crystal structure of the cyclic inhibitor in complex with trypsin. The lack of a major conformational change upon binding suggests that the structure of SFTI-1 is rigid and already pre-organized for maximal binding due to minimization of entropic losses compared to a more flexible ligand. These properties make SFTI-1 an ideal platform for the design of small peptidic pharmaceuticals or pesticides. (C) 2001 Academic Press.

A critical comparison of the external and internal boron requirements for contrasting species in boron-buffered solution culture

Despite reports that boron (B) requirements differ among plant species there is a shortage of critical evidence to demonstrate unequivocally whether species differ in internal or external B requirements or both. The present research was conducted to establish the external and internal B requirements of three contrasting species, a woody dicot (marri), an herbaceous dicot (sunflower) and a monocot (wheat) using B-buffered solution culture. Boron-buffered solution culture provided satisfactory control of external B concentrations ranging from 0.04 to 30 muM throughout the 20- (sunflower and wheat) or 40-day (marri) growth period. At low external B concentrations (less than or equal to 0.13 muM), the growth of marri and sunflower was severely depressed but by contrast the vegetative growth of wheat plants was satisfactory and free of B deficiency symptoms. Marri and sunflower plants achieved total maximum shoot growth at greater than or equal to1.2 muM B in solutions while wheat plants did so at greater than or equal to 0.6 muM B. The critical B concentrations (mg kg(-1) dry matter) in the youngest open leaf blades of marri, sunflower and wheat plants were 17.9, 19.7 and 1.2 on 20, 10 and 10 days after transplanting (DAT), respectively. Lower internal and external B requirements of wheat were matched by a lower uptake rate of B compared to marri and sunflower.

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 structure of BSTI: A new trypsin inhibitor from skin secretions of Bombina bombina

The three-dimensional solution structure of BSTI, a trypsin inhibitor from the European frog Bombina bombina, has been solved using H-1 NMR spectroscopy. The 60 amino acid protein contains five disulfide bonds, which were unambiguously determined to be Cvs (4-38), Cys (13-34), Cys (17-30), Cys (21-60), and Cys (40-54) by experimental restraints and subsequent structure calculations. The main elements of secondary structure are four beta -strands, arranged as two small antiparallel beta -sheets, The overall fold of BSTI is disk shaped and is characterized by the lack of a hydrophobic core. The presumed active site is located on a loop comprising residues 21-34, which is a relatively disordered region similar to that seen in many other protease inhibitors. However, the overall fold is different to other known protease inhibitors with the exception of a small family of inhibitors isolated from nematodes of the family Ascaris and recently also from the haemolymph of Apis mellifera. BSTI may thus be classified as a new member of this recently discovered family of protease inhibitors.

The solution structure of C1-T1, a two-domain proteinase inhibitor derived from a circular precursor protein from Nicotiana alata

A two-domain portion of the proteinase inhibitor precursor from Nicotiana alata (NaProPI) has been expressed and its structure determined by NMR spectroscopy. NaProPI contains six almost identical 53 amino acid repeats that fold into six highly similar domains; however, the sequence repeats do nut coincide with the structural domains. Five of the structural domains comprise the C-terminal portion of one repeat and the N-terminal portion of the next. The sixth domain contains the C-terminal portion of the sixth repeat and the N-terminal portion of the first repeat. Disulphide bonds link these C and N-terminal fragments to generate the clasped-bracelet fold of NaProPI. The three-dimensional structure of NaProPI is not known, but it is conceivable that adjacent domains in NaProPI interact to generate the circular bracelet with the N and C termini in close enough proximity to facilitate formation of the disulphide bonds that form the clasp The expressed protein, examined in the current study, comprises residues 25-135 of NaProPI and encompasses the first two contiguous structural domains, namely the chymotrypsin inhibitor C1 and the trypsin inhibitor T1, joined by a five-residue linker, and is referred to as C1-T1. The tertiary structure of each domain in C1-T1 is identical to that found in the isolated inhibitors. However, no nuclear Overhauser effect contacts are observed between the two domains and the five-residue linker adopts an extended conformation. The absence of interactions between the domains indicates that adjacent domains do not specifically interact to drive the circularisation of NaProPI. These results are in agreement with recent data which describe similar PI precursors from other members of the Solanaceae having two, three, or four repeats. The lack of strong interdomain association is likely to be important for the function of individual inhibitors by ensuring that there is no masking of reactive sites upon release from the precursor. (C) 2001 Academic Press.