Spectral induced polarization signatures of abiotic FeS precipitation

In recent years, geophysical methods have been shown to be sensitive to microbial-induced mineralization processes. The spectral induced-polarization (SIP) method appears to be very promising for monitoring mineralization and microbial processes. With this work, we study the links of mineralization and SIP signals, in the absence of microbial activity. We recorded the SIP response during abiotic FeS precipitation. We show that the SIP signals are diagnostic of FeS mineralization and can be differentiated from SIP signals from biomineralization processes. More specifically, the imaginary conductivity shows almost linear dependence on the amount of FeS precipitating out of solution, above the threshold value 0.006 gr under our experimental conditions. This research has direct implications for the use of the SIP method as a monitoring and decision-making tool for sustainable remediation of metals in contaminated soils and groundwater.

A Phase Grating Approach to Modeling Surface Diffusion in FDTD Room Acoustics Simulations

In this paper, a method for modeling diffusive boundaries in finite difference time domain (FDTD) room acoustics simulations with the use of impedance filters is presented. The proposed technique is based on the concept of phase grating diffusers, and realized by designing boundary impedance filters from normal-incidence reflection filters with added delay. These added delays, that correspond to the diffuser well depths, are varied across the boundary surface, and implemented using Thiran allpass filters. The proposed method for simulating sound scattering is suitable for modeling high frequency diffusion caused by small variations in surface roughness and, more generally, diffusers characterized by narrow wells with infinitely thin separators. This concept is also applicable to other wave-based modeling techniques. The approach is validated by comparing numerical results for Schroeder diffusers to measured data. In addition, it is proposed that irregular surfaces are modeled by shaping them with Brownian noise, giving good control over the sound scattering properties of the simulated boundary through two parameters, namely the spectral density exponent and the maximum well depth.

On the limits of spectral methods for frequency estimation

An algorithm is presented which generates pairs of oscillatory random time series which have identical periodograms but differ in the number of oscillations. This result indicates the intrinsic limitations of spectral methods when it comes to the task of measuring frequencies. Other examples, one from medicine and one from bifurcation theory, are given, which also exhibit these limitations of spectral methods. For two methods of spectral estimation it is verified that the particular way end points are treated, which is specific to each method, is, for long enough time series, not relevant for the main result.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

A 330-360 GHz spectral survey of G34.3+0.15 - III. The outer halo

We have undertaken a 330-360 GHz molecular line survey of the halo gas surrounding the hot core associated with G34.26+0.15. In contrast to our molecular line survey of the hot core itself, where 338 lines from at least 38 species were detected, only 18 lines from 9 species were detected in the halo. The lines are mainly single transitions of simple di atomic and triatomic molecules. Lower limits to their column densities have been evaluated by an LTE method. In the case of methanol, where four transitions were detected, the rotation temperature and column density have been evaluated by the rotation diagram technique. We have modified the previous depth-dependent chemical model developed in Paper II to calculate the column densities observed along a general line of sight drawn through the model cloud. The model is also extended to produce beam-averaged column densities for better comparison with those observed. We compare the model column densities with those observed and make recommendations for future depth-dependent chemical modelling of hot cores.

Design and evaluation of tinnitus synthesis methods: From spectral to spatial matching

Purpose
This study was designed to investigate methods to help patients suffering from unilateral tinnitus synthesizing an auditory replica of their tinnitus.

Materials and methods
Two semi-automatic methods (A and B) derived from the auditory threshold of the patient and a method (C) combining a pure tone and a narrow band-pass noise centred on an adjustable frequency were devised and rated on their likeness over two test sessions. A third test evaluated the stability over time of the synthesized tinnitus replica built with method C, and its proneness to merge with the patient's tinnitus. Patients were then asked to try and control the lateralisation of this single percept through the adjustment of the tinnitus replica level.

Results
The first two tests showed that seven out of ten patients chose the tinnitus replica built with method C as their preferred one. The third test, performed on twelve patients, revealed pitch tuning was rather stable over a week interval. It showed that eight patients were able to consistently match the central frequency of the synthesized tinnitus (presented to the contralateral ear) to their own tinnitus, which leaded to a unique tinnitus percept. The lateralisation displacement was consistent across patients and revealed an average range of 29dB to obtain a full lateral shift from the ipsilateral to the contralateral side.

Conclusions
Although spectrally simpler than the semi-automatic methods, method C could replicate patients' tinnitus, to some extent. When a unique percept between synthesized tinnitus and patients' tinnitus arose, lateralisation of this percept was achieved.

Dosimetry and spectral analysis of a radiobiological experiment using laser-driven proton beams

Laser-driven proton and ion acceleration is an area of increasing research interest given the recent development of short pulse-high intensity lasers. Several groups have reported experiments to understand whether a laser-driven beam can be applied for radiobiological purposes and in each of these, the method to obtain dose and spectral analysis was slightly different. The difficulty with these studies is that the very large instantaneous dose rate is a challenge for commonly used dosimetry techniques, so that other more sophisticated procedures need to be explored. This paper aims to explain a method for obtaining the energetic spectrum and the dose of a laser-driven proton beam irradiating a cell dish used for radiobiology studies. The procedure includes the use of a magnet to have charge and energy separation of the laser-driven beam, Gafchromic films to have information on dose and partially on energy, and a Monte Carlo code to expand the measured data in order to obtain specific details of the proton spectrum on the cells. Two specific correction factors have to be calculated: one to take into account the variation of the dose response of the films as a function of the proton energy and the other to obtain the dose to the cell layer starting from the dose measured on the films. This method, particularly suited to irradiation delivered in a single laser shot, can be applied in any other radiobiological experiment performed with laser-driven proton beams, with the only condition that the initial proton spectrum has to be at least roughly known. The method was tested in an experiment conducted at Queen's University of Belfast using the TARANIS laser, where the mean energy of the protons crossing the cells was between 0.9 and 5 MeV, the instantaneous dose rate was estimated to be close to 10(9) Gy s(-1) and doses between 0.8 and 5 Gy were delivered to the cells in a single laser shot. The combination of the applied corrections modified the initial estimate of dose by up to 40%.

The VLT-FLAMES survey of massive stars: wind properties and evolution of hot massive stars in the Large Magellanic Cloud

We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star forming region N11. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar parameters are determined using an automated fitting method ( Mokiem et al. 2005), which combines the stellar atmosphere code fastwind ( Puls et al. 2005) with the genetic algorithm based optimisation routine PIKAIA ( Charbonneau 1995). We derive an age of 7.0 +/- 1.0 and 3.0 +/- 1.0 Myr for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation scenario in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2. From helium and hydrogen line fitting we find the hottest three of these stars to be similar to 49- 54 kK ( compared to similar to 45- 46 kK for O3 stars). Detailed determination of the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. ( 2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non-rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z similar to 1/2 Z(circle dot) environment sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D-mom as a function of luminosity with predictions for LMC metallicities by Vink et al. ( 2001) yields good agreement in the entire luminosity range that was investigated, i.e. 5.0

An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

Physical and numerical constraints in source modeling for finite difference simulation of room acoustics

In finite difference time domain simulation of room acoustics, source functions are subject to various constraints. These depend on the way sources are injected into the grid and on the chosen parameters of the numerical scheme being used. This paper addresses the issue of selecting and designing sources for finite difference simulation, by first reviewing associated aims and constraints, and evaluating existing source models against these criteria. The process of exciting a model is generalized by introducing a system of three cascaded filters, respectively, characterizing the driving pulse, the source mechanics, and the injection of the resulting source function into the grid. It is shown that hard, soft, and transparent sources can be seen as special cases within this unified approach. Starting from the mechanics of a small pulsating sphere, a parametric source model is formulated by specifying suitable filters. This physically constrained source model is numerically consistent, does not scatter incoming waves, and is free from zero- and low-frequency artifacts. Simulation results are employed for comparison with existing source formulations in terms of meeting the spectral and temporal requirements on the outward propagating wave.

A proposed model membrane and test method for microneedle insertion studies

A commercial polymeric film (Parafilm M (R), a blend of a hydrocarbon wax and a polyolefin) was evaluated as a model membrane for microneedle (MN) insertion studies. Polymeric MN arrays were inserted into Parafilm M (R) (PF) and also into excised neonatal porcine skin. Parafilm M (R) was folded before the insertions to closely approximate thickness of the excised skin. Insertion depths were evaluated using optical coherence tomography (OCT) using either a force applied by a Texture Analyser or by a group of human volunteers. The obtained insertion depths were, in general, slightly lower, especially for higher forces, for PF than for skin. However, this difference was not a large, being less than the 10% of the needle length. Therefore, all these data indicate that this model membrane could be a good alternative to biological tissue for MN insertion studies. As an alternative method to OCT, light microscopy was used to evaluate the insertion depths of MN in the model membrane. This provided a rapid, simple method to compare different MN formulations. The use of Parafilm M (R), in conjunction with a standardised force/time profile applied by a Texture Analyser, could provide the basis for a rapid MN quality control test suitable for in-process use. It could also be used as a comparative test of insertion efficiency between candidate MN formulations. 

CLUSTER LEED - A VERY FAST METHOD FOR SURFACE-STRUCTURE DETERMINATION

A very fast method, cluster low-energy electron diffraction (LEED) is proposed for LEED I-V spectral analysis, in which three appproximations are introduced: the small-atom approximation, omission of the structure factors, and truncation of higher order ( > 2) scattering events. The method has been tested using a total of four sets of I-V spectra calculated by fully dynamic LEED for (i) the simple overlayer system, O on Ni{100}, and (ii) the reconstructed system, Cu on W{100}, and also one set of experimental data from W{100}-c(2 X 2)-Cu. In each case the correct structural parameters are recovered. It is suggested that for complex systems cluster LEED provides an efficient fast route to trial structures, which could be refined by automated tenser LEED.

Low complexity greedy detection method with generalized multicarrier index keying OFDM

Multicarrier Index Keying (MCIK) is a recently developed technique that modulates subcarriers but also indices of the subcarriers. In this paper a novel low-complexity detection scheme of subcarrier indices is proposed for an MCIK system and addresses a substantial reduction in complexity over the optimalmaximum likelihood (ML) detection. For the performance evaluation, a closed-form expression for the pairwise error probability (PEP) of an active subcarrier index, and a tight approximation of the average PEP of multiple subcarrier indices are derived in closed-form. The theoretical outcomes are validated usingsimulations, at a difference of less than 0.1dB. Compared to the optimal ML, the proposed detection achieves a substantial reduction in complexity with small loss in error performance (<= 0.6dB).