Stability of nonlinear difference inclusions

The stability of difference inclusions x(k+1) is an element of F(x(k)) is studied, where F(x) = {F(x, gimel) : is an element of Lambda} and the selections F(., gimel) : E -->E assume values in a Banach space E, partially ordered by a cone K. It is assumed that the operators F(.,gimel) are heterotone or pseudoconcave. The main results concern asymptotically stable absorbing sets, and include the case of a single equilibrium point. The results are applied to a number of practical problems.

Step size control for efficient discrete element simulation

The step size determines the accuracy of a discrete element simulation. The position and velocity updating calculation uses a pre-calculated table and hence the control of step size can not use the integration formulas for step size control. A step size control scheme for use with the table driven velocity and position calculation uses the difference between the calculation result from one big step and that from two small steps. This variable time step size method chooses the suitable time step size for each particle at each step automatically according to the conditions. Simulation using fixed time step method is compared with that of using variable time step method. The difference in computation time for the same accuracy using a variable step size (compared to the fixed step) depends on the particular problem. For a simple test case the times are roughly similar. However, the variable step size gives the required accuracy on the first run. A fixed step size may require several runs to check the simulation accuracy or a conservative step size that results in longer run times. (C) 2001 Elsevier Science Ltd. All rights reserved.

Elucidation of reaction scheme describing malondialdehyde-acetaldehyde-protein adduct formation

Malondialdehyde and acetaldehyde react together with proteins and form hybrid protein conjugates designated as MAA adducts, which have been detected in livers of ethanol-fed animals. Our previous studies have shown that MAA adducts are comprised of two distinct products. One adduct is composed of two molecules of malondialdehyde and one molecule of acetaldehyde and was identified as the 4-methpl-1,4-dihydropyridine-3,5-dicarbaldehyde derivative of an amino group (MHHDC adduct). The other adduct is a 1:1 adduct of malondialdehyde and acetaldehyde and was identified as the 2-formyl-3-(alkylamino)butanal derivative of an amino group (FAAB adduct). In this study, information on the mechanism of MAA adduct formation was obtained, focusing on whether the FAAB adduct serves as a precursor for the MDHDC adduct. Upon the basis of chemical analysis and NMR spectroscopy, two initial reaction steps appear to be a prerequisite for MDHDC formation. One step involves the reaction of one molecule of malondialdehyde and one of acetaldehyde with an amino group of a protein to form the FAAB product, while the other step involves the generation of a malondialdehyde-enamine. It appears that generation of the MDHDC adduct requires the FAAB moiety to be transferred to the nitrogen of the MDA-enamine. For efficient reaction of FAAB with the enamine to take place, additional experiments indicated that these two intermediates likely must be in positions on the protein of close proximity to each other. Further studies showed that the incubation of liver proteins from ethanol-fed rats with MDA resulted in a marked generation of MDHDC adducts, indicating the presence of a pool of FAAB adducts in the liver of ethanol-fed animals. Overall, these findings show that MDHDC-protein adduct formation occurs via the reaction of the FAAB moiety with a malondialdehyde-enamine, and further suggest that a similar mechanism may be operative in vivo in the liver during prolonged ethanol consumption.

Compact clinical MRI magnet design using a multilayer current density approach

In this work, a new method of optimization is successfully applied to the theoretical design of compact, actively shielded, clinical MRI magnets. The problem is formulated as a two-step process in which the desired current densities on multiple, cc-axial surface layers are first calculated by solving Fredholm equations of the first kind. Non-linear optimization methods with inequality constraints are then invoked to fit practical magnet coils to the desired current densities. The current density approach allows rapid prototyping of unusual magnet designs. The emphasis of this work is on the optimal design of short, actively-shielded MRI magnets for whole-body imaging. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric, and asymmetric MRI magnets. Magnet designs are presented for actively-shielded, symmetric magnets of coil length 1.0 m, which is considerably shorter than currently available designs of comparable dsv size. Novel, actively-shielded, asymmetric magnet designs are also presented in which the beginning of a 50-cm dsv is positioned just 11 cm from the end of the coil structure, allowing much improved access to the patient and reduced patient claustrophobia. Magn Reson Med 45:331540, 2001. (C) 2001 Wiley-Liss, Inc.

The stochastic design of force-minimized compact magnets for high-field magnetic resonance imaging applications

New designs for force-minimized compact high-field clinical MRI magnets are described. The design method is a modified simulated annealing (SA) procedure which includes Maxwell forces in the error function to be minimized. This permits an automated force reduction in the magnet designs while controlling the overall dimensions of the system. As SA optimization requires many iterations to achieve a final design, it is important that each iteration in the procedure is rapid. We have therefore developed a rapid force calculation algorithm. Novel designs for short 3- and 4-T clinical MRI systems are presented in which force reduction has been invoked. The final designs provide large homogeneous regions and reduced stray fields in remarkable short magnets. A shielded 4-T design that is approximately 30% shorter than current designs is presented. This novel magnet generates a full 50-cm diameter homogeneous region.

Towards a definition of "difference" in osteoarthritis

To assess existing information regarding detectable differences in osteoarthritis (OA), a systematic literature search was conducted up to December 1999. Thirty-three articles were considered methodologically relevant to the definition and categorization of detectable differences in OA. It was determined that the musculoskeletal literature contains a wealth of information that relates to observed changes, much of which is derived from the clinical trials literature, but there have been relatively few methodological studies that have systematically evaluated the nature, categorization, and relevance of the change. Furthermore, most of those that have been published take the perspective of an individual or groups of experts other than that of the patient. This summary of the current literature reveals that the diverse sources of information go part way towards developing an understanding of detectable differences and their importance in the area of OA research and clinical practice. Stakeholders' interests as well as factors that modulate perceptions of importance need to be taken under consideration. In particular, the patient's perspective of the importance of change at an individual level requires further evaluation. This area of clinical research is relatively underdeveloped, but there is considerable opportunity for progress.

Current density mapping approach for design of clinical magnetic resonance imaging magnets

Novel current density mapping (CDM) schemes are developed for the design of new actively shielded, clinical magnetic resonance imaging (MRI) magnets. This is an extended inverse method in which the entire potential solution space for the superconductors has been considered, rather than single current density layers. The solution provides an insight into the required superconducting coil pattern for a desired magnet configuration. This information is then used as an initial set of parameters for the magnet structure, and a previously developed hybrid numerical optimization technique is used to obtain the final geometry of the magnet. The CDM scheme is applied to the design of compact symmetric, asymmetric, and open architecture 1.0-1.5 T MRI magnet systems of novel geometry and utility. A new symmetric 1.0-T system that is just I m in length with a full 50-cm diameter of the active, or sensitive, volume (DSV) is detailed, as well as an asymmetric system in which a 50-cm DSV begins just 14 cm from the end of the coil structure. Finally a 1.0-T open magnet system with a full 50-cm DSV is presented. These new designs provide clinically useful homogeneous regions and have appropriately restricted stray fields but, in some of the designs, the DSV is much closer to the end of the magnet system than in conventional designs. These new designs have the potential to reduce patient claustrophobia and improve physician access to patients undergoing scans. (C) 2002 Wiley Periodicals, Inc.

Boundary value problems for systems of difference equations associated with systems of second-order ordinary differential equations

We study difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order ordinary differential equations. We formulate conditions which guarantee a priori bounds on first differences of solutions to the discretized problem. We establish existence results for solutions to the discretized boundary value problems subject to nonlinear boundary conditions. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. (C) 2002 Elsevier Science Ltd. All rights reserved.

Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations

We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) = (y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assume there exist strict lower and strict upper solutions and impose additional conditions on f and G which are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. We show that the discrete approximation also has solutions which approximate solutions of the continuous problem and converge to the solution of the continuous problem when it is unique, as the grid size goes to 0. Homotopy methods can be used to compute the solution of the discrete approximation. Our results were motivated by those of Gaines.

A new joint coding and modulation scheme for channels with clipping

A major limitation in any high-performance digital communication system is the linearity region of the transmitting amplifier. Nonlinearities typically lead to signal clipping. Efficient communication in such conditions requires maintaining a low peak-to-average power ratio (PAR) in the transmitted signal while achieving a high throughput of data. Excessive PAR leads either to frequent clipping or to inadequate resolution in the analog-to-digital or digital-to-analog converters. Currently proposed signaling schemes for future generation wireless communications suffer from a high PAR. This paper presents a new signaling scheme for channels with clipping which achieves a PAR as low as 3. For a given linear range in the transmitter's digital-to-analog converter, this scheme achieves a lower bit-error rate than existing multicarrier schemes, owing to increased separation between constellation points. We present the theoretical basis for this new scheme, approximations for the expected bit-error rate, and simulation results. (C) 2002 Elsevier Science (USA).

Practical scheme for quantum computation with any two-qubit entangling gate

Which gates are universal for quantum computation? Although it is well known that certain gates on two-level quantum systems (qubits), such as the controlled-NOT, are universal when assisted by arbitrary one-qubit gates, it has only recently become clear precisely what class of two-qubit gates is universal in this sense. We present an elementary proof that any entangling two-qubit gate is universal for quantum computation, when assisted by one-qubit gates. A proof of this result for systems of arbitrary finite dimension has been provided by Brylinski and Brylinski; however, their proof relies on a long argument using advanced mathematics. In contrast, our proof provides a simple constructive procedure which is close to optimal and experimentally practical.

Copper transfer from the Cu(I) chaperone, CopZ, to the repressor, Zn(II)CopY: Metal coordination environments and protein interactions

Extracellular copper regulates the DNA binding activity of the CopY repressor of Enterococcus hirae and thereby controls expression of the copper homeostatic genes encoded by the cop operon. CopY has a CxCxxxxCxC metal binding motif. CopZ, a copper chaperone belonging to a family of metallochaperones characterized by a MxCxxC metal binding motif, transfers copper to CopY. The copper binding stoichiometries of CopZ and CopY were determined by in vitro metal reconstitutions. The stoichiometries were found to be one copper(l) per CopZ and two copper(l) per CopY monomer. X-ray absorption studies suggested a mixture of two- and three-coordinate copper in Cu(1)CopZ, but a purely three-coordinate copper coordination with a Cu-Cu interaction for Cu(1)(2)CopY. The latter coordination is consistent with the formation of a compact binuclear Cu(l)-thiolate core in the CxCxxxxCxC binding motif of CopY. Displacement of zinc, by copper. from CopY was monitored with 2,4-pyridylazoresorcinol. Two copper(l) ions were required to release the single zinc(II) ion bound per CopY monomer. The specificity of copper transfer between CopZ and CopY was dependent on electrostatic interactions. Relative copper binding affinities of the proteins were investigated using the chelator, diethyldithiocarbamic acid (DDC). These data suggest that CopY has a higher affinity for copper than CopZ. However, this affinity difference is not the sole factor in the copper exchange: a charge-based interaction between the two proteins is required for the transfer reaction to proceed. Gain-of-function mutation of a CopZ homologue demonstrated the necessity of four lysine residues on the chaperone for the interaction with CopY. Taken together, these results suggest a mechanism for copper exchange between CopZ and CopY.