Numerical evaluation of the fields induced by body motion in or near high-field MRI scanners

In modern magnetic resonance imaging, both patients and health care workers are exposed to strong. non-uniform static magnetic fields inside and outside of the scanner. In which body movement may be able to induce electric currents in tissues which could be potentially harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in a tissue-equivalent human model when moving at various positions around the magnet. The numerical calculations are based on an efficient. quasi-static, finite-difference scheme. Three-dimensional field profiles from an actively shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The methodology presented herein can be extrapolated to very high field strengths for the evaluation of the effects of motion at a variety of field strengths and velocities. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical evaluation of e-fields induced by body motion near high-field MRI scanner

In modern magnetic resonance imaging (MRI), both patients and radiologists are exposed to strong, nonuniform static magnetic fields inside or outside of the scanner, in which the body movement may be able to induce electric currents in tissues which could be possibly harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in the human model when moving at various positions around the magnet. The numerical calculations are based on an efficient, quasistatic, finite-difference scheme and an anatomically realistic, full-body, male model. 3D field profiles from an actively-shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The results are easy to extrapolate to very high field strengths for the safety evaluation at a variety of field strengths and motion velocities.

Modelling the input-output behaviour of piezoelectric structural health monitoring systems for composite plates

This paper investigates the input-output characteristics of structural health monitoring systems for composite plates based on permanently attached piezoelectric transmitter and sensor elements. Using dynamic piezoelectricity theory and a multiple integral transform method to describe the propagating and scattered flexural waves an electro-mechanical model for simulating the voltage input-output transfer function for circular piezoelectric transmitters and sensors adhesively attached to an orthotropic composite plate is developed. The method enables the characterization of all three physical processes, i.e. wave generation, wave propagation and wave reception. The influence of transducer, plate and attached electrical circuit characteristics on the voltage output behaviour of the system is examined through numerical calculations, both in frequency and the time domain. The results show that the input-output behaviour of the system is not properly predicted by the transducers' properties alone. Coupling effects between the transducers and the tested structure have to be taken into account, and adding backing materials to the piezoelectric elements can significantly improve the sensitivity of the system. It is shown that in order to achieve maximum sensitivity, particular piezoelectric transmitters and sensors need to be designed according to the structure to be monitored and the specific frequency regime of interest.

Exact solution, scaling behaviour and quantum dynamics of a model of an atom-molecule Bose-Einstein condensate

0We study the exact solution for a two-mode model describing coherent coupling between atomic and molecular Bose-Einstein condensates (BEC), in the context of the Bethe ansatz. By combining an asymptotic and numerical analysis, we identify the scaling behaviour of the model and determine the zero temperature expectation value for the coherence and average atomic occupation. The threshold coupling for production of the molecular BEC is identified as the point at which the energy gap is minimum. Our numerical results indicate a parity effect for the energy gap between ground and first excited state depending on whether the total atomic number is odd or even. The numerical calculations for the quantum dynamics reveals a smooth transition from the atomic to the molecular BEC.

Effects of hydrogen-air non-equilibrium chemistry on the performance of a model scramjet thrust nozzle

Two aspects of hydrogen-air non-equilibrium chemistry related to scramjets are nozzle freezing and a process called 'kinetic afterburning' which involves continuation of combustion after expansion in the nozzle. These effects were investigated numerically and experimentally with a model scramjet combustion chamber and thrust nozzle combination. The overall model length was 0.5m, while precombustion Mach numbers of 3.1 +/- 0.3 and precombustion temperatures ranging from 740K to 1,400K were involved. Nozzle freezing was investigated at precombustion pressures of 190kPa and higher, and it was found that the nozzle thrusts were within 6% of values obtained from finite rate numerical calculations, which were within 7% of equilibrium calculations. When precombustion pressures of 70kPa or less were used, kinetic afterburning was found to be partly responsible for thrust production, in both the numerical calculations and the experiments. Kinetic afterburning offers a means of extending the operating Mach number range of a fixed geometry scramjet.

Entanglement and bifurcations in Jahn-Teller models

We compare and contrast the entanglement in the ground state of two Jahn-Teller models. The Exbeta system models the coupling of a two-level electronic system, or qubit, to a single-oscillator mode, while the Exepsilon models the qubit coupled to two independent, degenerate oscillator modes. In the absence of a transverse magnetic field applied to the qubit, both systems exhibit a degenerate ground state. Whereas there always exists a completely separable ground state in the Exbeta system, the ground states of the Exepsilon model always exhibit entanglement. For the Exbeta case we aim to clarify results from previous work, alluding to a link between the ground-state entanglement characteristics and a bifurcation of a fixed point in the classical analog. In the Exepsilon case we make use of an ansatz for the ground state. We compare this ansatz to exact numerical calculations and use it to investigate how the entanglement is shared between the three system degrees of freedom.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

Orbital currents and charge density waves in a generalized Hubbard ladder

We study a generalized Hubbard model on the two-leg ladder at zero temperature, focusing on a parameter region with staggered flux (SF)/d-density wave (DDW) order. To guide our numerical calculations, we first investigate the location of a SF/DDW phase in the phase diagram of the half-filled weakly interacting ladder using a perturbative renormalization group (RG) and bosonization approach. For hole doping 6 away from half-filling, finite-system density-matrix renormalizationgroup (DMRG) calculations are used to study ladders with up to 200 rungs for intermediate-strength interactions. In the doped SF/DDW phase, the staggered rung current and the rung electron density both show periodic spatial oscillations, with characteristic wavelengths 2/delta and 1/delta, respectively, corresponding to ordering wavevectors 2k(F) and 4k(F) for the currents and densities, where 2k(F) = pi(1 - delta). The density minima are located at the anti-phase domain walls of the staggered current. For sufficiently large dopings, SF/DDW order is suppressed. The rung density modulation also exists in neighboring phases where currents decay exponentially. We show that most of the DMRG results can be qualitatively understood from weak-coupling RG/bosonization arguments. However, while these arguments seem to suggest a crossover from non-decaying correlations to power-law decay at a length scale of order 1/delta, the DMRG results are consistent with a true long-range order scenario for the currents and densities. (c) 2005 Elsevier Inc. All rights reserved.

Distillability of Werner states using entanglement witnesses and robust semidefinite programs

We use robust semidefinite programs and entanglement witnesses to study the distillability of Werner states. We perform exact numerical calculations that show two-undistillability in a region of the state space, which was previously conjectured to be undistillable. We also introduce bases that yield interesting expressions for the distillability witnesses and for a tensor product of Werner states with an arbitrary number of copies.

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f-shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self-interaction correction and LIDA + U, the f-(non-f)-mixing interaction is decreased by spectral weights of the many-body states of the f-ion. (c) 2005 Wiley Periodicals, Inc.