Multipole expansion of strongly focussed laser beams

Multipole expansion of an incident radiation field-that is, representation of the fields as sums of vector spherical wavefunctions-is essential for theoretical light scattering methods such as the T-matrix method and generalised Lorenz-Mie theory (GLMT). In general, it is theoretically straightforward to find a vector spherical wavefunction representation of an arbitrary radiation field. For example, a simple formula results in the useful case of an incident plane wave. Laser beams present some difficulties. These problems are not a result of any deficiency in the basic process of spherical wavefunction expansion, but are due to the fact that laser beams, in their standard representations, are not radiation fields, but only approximations of radiation fields. This results from the standard laser beam representations being solutions to the paraxial scalar wave equation. We present an efficient method for determining the multipole representation of an arbitrary focussed beam. (C) 2003 Elsevier Science Ltd. All rights reserved.

T-Q relation and exact solution for the XYZ chain with general non-diagonal boundary terms

We propose that the Baxter's Q-operator for the quantum XYZ spin chain with open boundary conditions is given by the j -> infinity limit of the corresponding transfer matrix with spin-j (i.e., (2j + I)-dimensional) auxiliary space. The associated T-Q relation is derived from the fusion hierarchy of the model. We use this relation to determine the Bethe Ansatz solution of the eigenvalues of the fundamental transfer matrix. The solution yields the complete spectrum of the Hamiltonian. (c) 2006 Elsevier B.V. All rights reserved.

Calculation of optical trapping landscapes

Manipulation of micrometer sized particles with optical tweezers can be precisely modeled with electrodynamic theory using Mie's solution for spherical particles or the T-matrix method for more complex objects. We model optical tweezers for a wide range of parameters including size, relative refractive index and objective numerical aperture. We present the resulting landscapes of the trap stiffness and maximum applicable trapping force in the parameter space. These landscapes give a detailed insight into the requirements and possibilities of optical trapping and provide detailed information on trapping of nanometer sized particles or trapping of high index particles like diamond.

Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

Effect of anode material and cavity design on the performance of microcavity OLEDs

We report on the effect of the replacement of the conventional ITO anode with the semitransparent metallic material on the performance of microcavity OLEDs. We performed comprehensive simulations of the emission from microcavity OLEDs consisting of widely used organic materials, N,N′-di(naphthalene-1- yl)-N,N′-diphenylbenzidine (NPB) as a hole transport layer and tris (8-hydroxyquinoline) (Alq3) as emitting and electron transporting layer. Silver and LiF/Al were considered as a cathode, while metallic (Au and Ag) anode was used and simulations were performed on devices with both the metallic and conventional ITO anode. The electroluminescence emission spectra, electric field distribution inside the device, carrier density, recombination rate and exciton density were calculated as a function of the position of the emission layer. The results show that the metallic anode enhances light output and that optimum emission from a microcavity OLED is achieved when the position of the recombination region is aligned with the antinode of the standing wave inside the cavity. The microcavity OLED devices with Ag/Ag and Ag/Au mirrors were fabricated and characterized. The experimental results have been compared to the simulations and the influence of the different anode, emission region width and position on the performance of microcavity OLEDs was discussed.

Electromagnetic interference from variable speed motor drives

A methodology is presented which can be used to produce the level of electromagnetic interference, in the form of conducted and radiated emissions, from variable speed drives, the drive that was modelled being a Eurotherm 583 drive. The conducted emissions are predicted using an accurate circuit model of the drive and its associated equipment. The circuit model was constructed from a number of different areas, these being: the power electronics of the drive, the line impedance stabilising network used during the experimental work to measure the conducted emissions, a model of an induction motor assuming near zero load, an accurate model of the shielded cable which connected the drive to the motor, and finally the parasitic capacitances that were present in the drive modelled. The conducted emissions were predicted with an error of +/-6dB over the frequency range 150kHz to 16MHz, which compares well with the limits set in the standards which specify a frequency range of 150kHz to 30MHz. The conducted emissions model was also used to predict the current and voltage sources which were used to predict the radiated emissions from the drive. Two methods for the prediction of the radiated emissions from the drive were investigated, the first being two-dimensional finite element analysis and the second three-dimensional transmission line matrix modelling. The finite element model took account of the features of the drive that were considered to produce the majority of the radiation, these features being the switching of the IGBT's in the inverter, the shielded cable which connected the drive to the motor as well as some of the cables that were present in the drive.The model also took account of the structure of the test rig used to measure the radiated emissions. It was found that the majority of the radiation produced came from the shielded cable and the common mode currents that were flowing in the shield, and that it was feasible to model the radiation from the drive by only modelling the shielded cable. The radiated emissions were correctly predicted in the frequency range 30MHz to 200MHz with an error of +10dB/-6dB. The transmission line matrix method modelled the shielded cable which connected the drive to the motor and also took account of the architecture of the test rig. Only limited simulations were performed using the transmission line matrix model as it was found to be a very slow method and not an ideal solution to the problem. However the limited results obtained were comparable, to within 5%, to the results obtained using the finite element model.

Study of Periodic and Quasi-Periodic Structures in Silicon-on-Insulator

In this thesis, a numerical design approach has been proposed and developed based on the transmission matrix method in order to characterize periodic and quasi-periodic photonic structures in silicon-on-insulator. The approach and its performance have been extensively tested with specific structures in 2D and its validity has been verified in 3D.

Estimation of sea surface temperature from sediment core PS1768-8

The quantitative diatom analysis of 218 surface sediment samples recovered in the Atlantic and western Indian sector of the Southern Ocean is used to define a base of reference data for paleotemperature estimations from diatom assemblages using the Imbrie and Kipp transfer function method. The criteria which justify the exclusion of samples and species out of the raw data set in order to define a reference database are outlined and discussed. Sensitivity tests with eight data sets were achieved evaluating the effects of overall dominance of single species, different methods of species abundance ranking, and no-analog conditions (e.g., Eucampia Antarctica) on the estimated paleotemperatures. The defined transfer functions were applied on a sediment core from the northern Antarctic zone. Overall dominance of Fragilariopsis kerguelensis in the diatom assemblages resulted in a close affinity between paleotemperature curve and relative abundance pattern of this species downcore. Logarithmic conversion of counting data applied with other ranking methods in order to compensate the dominance of F. kerguelensis revealed the best statistical results. A reliable diatom transfer function for future paleotemperature estimations is presented.

LGM, planktic foraminifera (EPILOG version)

[1] We used planktic foraminiferal assemblages in 70 sediment cores from the tropical and subtropical South Atlantic Ocean (10°N-37°S) to estimate annual mean sea surface temperatures (SSTs) and seasonality for the Last Glacial Maximum with a modified version of the Imbrie-Kipp transfer function method (IKTF) that takes into account the abundance of rare but temperature sensitive species. In contrast to CLIMAP Project Members [1981], the reconstructed SSTs indicate cooler glacial SSTs in the entire tropical/subtropical South Atlantic with strongest cooling in the upwelling region off Namibia (7-10°C) and smallest cooling (1-2°C) in the western subtropical gyre. In the western Atlantic, our data support recent temperature estimates from other proxies. In the upwelling regions in the eastern Atlantic, our data conflict with SST reconstructions from alkenones, which may be due to an environmental preference of the alkenone-producing algae or to an underestimation of foraminiferal SSTs due to anomalous high abundances of N. pachyderma (sinistral).

Pole placement method using the system matrix transfer function and sparsity

This paper presents a pole placement method using both the augmented Jacobian and the corresponding system transfer function matrices. From the manipulation of these matrices a straightforward approach results to get the coefficients of a non-linear system, whose solution gives the parameters of the stabilizers that can provide a pre-specified minimum damping to the system. (C) 2001 Elsevier B.V. Ltd. All rights reserved.

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.

Matrix Converter-Based Unified Power-Flow Controllers: Advanced Direct Power Control Method

This paper presents a direct power control (DPC) for three-phase matrix converters operating as unified power flow controllers (UPFCs). Matrix converters (MCs) allow the direct ac/ac power conversion without dc energy storage links; therefore, the MC-based UPFC (MC-UPFC) has reduced volume and cost, reduced capacitor power losses, together with higher reliability. Theoretical principles of direct power control (DPC) based on sliding mode control techniques are established for an MC-UPFC dynamic model including the input filter. As a result, line active and reactive power, together with ac supply reactive power, can be directly controlled by selecting an appropriate matrix converter switching state guaranteeing good steady-state and dynamic responses. Experimental results of DPC controllers for MC-UPFC show decoupled active and reactive power control, zero steady-state tracking error, and fast response times. Compared to an MC-UPFC using active and reactive power linear controllers based on a modified Venturini high-frequency PWM modulator, the experimental results of the advanced DPC-MC guarantee faster responses without overshoot and no steady-state error, presenting no cross-coupling in dynamic and steady-state responses.

Development of a new method for the characterisation of bioreactors concerning heat and mass transfer

Since the specific heat transfer coefficient (UA) and the volumetric mass transfer coefficient (kLa) play an important role for the design of biotechnological processes, different techniques were developed in the past for the determination of these parameters. However, these approaches often use imprecise dynamic methods for the description of stationary processes and are limited towards scale and geometry of the bioreactor. Therefore, the aim of this thesis was to develop a new method, which overcomes these restrictions. This new approach is based on a permanent production of heat and oxygen by the constant decomposition of hydrogen peroxide in continuous mode. Since the degradation of H2O2 at standard conditions only takes place by the support of a catalyst, different candidates were investigated for their potential (regarding safety issues and reaction kinetic). Manganese-(IV)-oxide was found to be suitable. To compensate the inactivation of MnO2, a continuous process with repeated feeds of fresh MnO2 was established. Subsequently, a scale-up was successfully carried out from 100 mL to a 5 litre glass bioreactor (UniVessel®)To show the applicability of this new method for the characterisation of bioreactors, it was compared with common approaches. With the newly established technique as well as with a conventional procedure, which is based on an electrical heat source, specific heat transfer coefficients were measured in the range of 17.1 – 24.8 W/K for power inputs of about 50 – 70 W/L. However, a first proof of concept regarding the mass transfer showed no constant kLa for different dilution rates up to 0.04 h-1.Based on this, consecutive studies concerning the mass transfer should be made with higher volume flows, due to more even inflow rates. In addition, further experiments are advisable, to analyse the heat transfer in single-use bioreactors and in larger common systems.

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack