Control theoretic smoothing splines

Some of the relationships between optimal control and statistics are examined. We produce generalized, smoothing splines by solving an optimal control problem for linear control systems, minimizing the L²-norm of the control signal, while driving the scalar output of the control system close to given, prespecified interpolation points. We then prove a convergence result for the smoothing splines, using results from the theory of numerical quadrature. Finally, we show, in simulations, that our approach works in practice as well as in theory     

Local Dynamic Stability Assessment of Motion Impaired Elderly Using Electronic Textile Pants

A clear association has been demonstrated between gait stability and falls in the elderly. Integration of wearable computing and human dynamic stability measures into home automation systems may help differentiate fall-prone individuals in a residential environment. The objective of the current study was to evaluate the capability of a pair of electronic textile (e-textile) pants system to assess local dynamic stability and to differentiate motion-impaired elderly from their healthy counterparts. A pair of e-textile pants comprised of numerous e-TAGs at locations corresponding to lower extremity joints was developed to collect acceleration, angular velocity and piezoelectric data. Four motion-impaired elderly together with nine healthy individuals (both young and old) participated in treadmill walking with a motion capture system simultaneously collecting kinematic data. Local dynamic stability, characterized by maximum Lyapunov exponent, was computed based on vertical acceleration and angular velocity at lower extremity joints for the measurements from both e-textile and motion capture systems. Results indicated that the motion-impaired elderly had significantly higher maximum Lyapunov exponents (computed from vertical acceleration data) than healthy individuals at the right ankle and hip joints. In addition, maximum Lyapunov exponents assessed by the motion capture system were found to be significantly higher than those assessed by the e-textile system. Despite the difference between these measurement techniques, attaching accelerometers at the ankle and hip joints was shown to be an effective sensor configuration. It was concluded that the e-textile pants system, via dynamic stability assessment, has the potential to identify motion-impaired elderly.     

Fabrication of planar cobalt electrodes separated by a sub-10nm gap using high resolution electron beam lithography with negative PMMA

We present a fabrication process of cobalt nanoelectrodes compatible with spin-dependent transport measurements through a few or a single nano-object. It consists in etching a cobalt thin layer into pairs of planar nanoelectrodes separated by a nanometric gap using a negative Poly-MethylMethAcrylate (PMMA) mask patterned by high resolution electron beam lithography (HREBL). The irradiation parameters of 200keV HREBL on PMMA have been investigated using atomic force microscopy (AFM) to define accurately the PMMA transformation from positive to negative tone. The influence of the electron dose and the designed gap on the final gap between electrodes is presented. This complete study proves that PMMA can be used as a HREBL negative resist to fabricate nanoelectrodes separated by a controlled and reproducible gap ranging from 5nm to several tens of nanometers.     

tweezercalib 2.0: Faster version of MatLab package for precise calibration of optical tweezers

We present a vectorized version of the MatLab (MathWorks Inc.) package tweezercalib for calibration of optical tweezers with precision. The calibration is based on the power spectrum of the Brownian motion of a dielectric bead trapped in the tweezers. Precision is achieved by accounting for a number of factors that affect this power spectrum, as described in vs. 1 of the package [I.M. Toli?-Nørrelykke, K. Berg-Sørensen, H. Flyvbjerg, Matlab program for precision calibration of optical tweezers, Comput. Phys. Comm. 159 (2004) 225-240]. The graphical user interface allows the user to include or leave out each of these factors. Several “health tests” are applied to the experimental data during calibration, and test results are displayed graphically. Thus, the user can easily see whether the data comply with the theory used for their interpretation. Final calibration results are given with statistical errors and covariance matrix.

New version program summary

Title of program: tweezercalibCatalogue identifier: ADTV_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADTV_v2_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandReference in CPC to previous version: I.M. Toli?-Nørrelykke, K. Berg-Sørensen, H. Flyvbjerg, Comput. Phys. Comm. 159 (2004) 225Catalogue identifier of previous version: ADTVDoes the new version supersede the original program: YesComputer for which the program is designed and others on which it has been tested: General computer running MatLab (Mathworks Inc.)Operating systems under with the program has been tested: Windows2000, Windows-XP, LinuxProgramming language used: MatLab (Mathworks Inc.), standard licenseMemory required to execute with typical data: Of order four times the size of the data fileHigh speed storage required: noneNo. of lines in distributed program, including test data, etc.: 135 989No. of bytes in distributed program, including test data, etc.: 1 527 611Distribution format: tar. gzNature of physical problem: Calibrate optical tweezers with precision by fitting theory to experimental power spectrum of position of bead doing Brownian motion in incompressible fluid, possibly near microscope cover slip, while trapped in optical tweezers. Thereby determine spring constant of optical trap and conversion factor for arbitrary-units-to-nanometers for detection system.Method of solution: Elimination of cross-talk between quadrant photo-diode's output channels for positions (optional). Check that distribution of recorded positions agrees with Boltzmann distribution of bead in harmonic trap. Data compression and noise reduction by blocking method applied to power spectrum. Full accounting for hydrodynamic effects: Frequency-dependent drag force and interaction with nearby cover slip (optional). Full accounting for electronic filters (optional), for “virtual filtering” caused by detection system (optional). Full accounting for aliasing caused by finite sampling rate (optional). Standard non-linear least-squares fitting. Statistical support for fit is given, with several plots facilitating inspection of consistency and quality of data and fit.Summary of revisions: A faster fitting routine, adapted from [J. Nocedal, Y.x. Yuan, Combining trust region and line search techniques, Technical Report OTC 98/04, Optimization Technology Center, 1998; W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, Numerical Recipes. The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986], is applied. It uses fewer function evaluations, and the remaining function evaluations have been vectorized. Calls to routines in Toolboxes not included with a standard MatLab license have been replaced by calls to routines that are included in the present package. Fitting parameters are rescaled to ensure that they are all of roughly the same size (of order 1) while being fitted. Generally, the program package has been updated to comply with MatLab, vs. 7.0, and optimized for speed.Restrictions on the complexity of the problem: Data should be positions of bead doing Brownian motion while held by optical tweezers. For high precision in final results, data should be time series measured over a long time, with sufficiently high experimental sampling rate: The sampling rate should be well above the characteristic frequency of the trap, the so-called corner frequency. Thus, the sampling frequency should typically be larger than 10 kHz. The Fast Fourier Transform used works optimally when the time series contain n² data points, and long measurement time is obtained with n>12-15. Finally, the optics should be set to ensure a harmonic trapping potential in the range of positions visited by the bead. The fitting procedure checks for harmonic potential.Typical running time: SecondsUnusual features of the program: NoneReferences: The theoretical underpinnings for the procedure are found in [K. Berg-Sørensen, H. Flyvbjerg, Power spectrum analysis for optical tweezers, Rev. Sci. Ins. 75 (2004) 594-612].     

Hydration of tricalcium germanate

The hydration characteristics of tricalcium germanate were examined. Kinetics of tricalcium germanate hydration, and hydration product morphologies and compositions were determined at 5 °C intervals between 5 and 70 °C. Complete hydration was rapidly achieved, X-ray diffraction indicated that tricalcium germanate reacted completely within 2 h at all the temperatures investigated. However, a constant value for an Arrhenius activation energy could not be determined. At hydration temperatures above 45 °C a small heat peak, which preceded the main calorimetric peak, appeared in the calorimetric curves. Calcium hydroxide was rarely observed by SEM for hydration at temperatures below 45 °C. Dense regions of calcium hydroxide were readily observed for hydration above this temperature. In contrast to the hydration of tricalcium silicate, an induction period was not observed nor was an extended period of diffusionally controlled reaction. Calcium germanate hydrate was fibrous with the fibre thickness exhibiting a temperature dependence. Hydration product compositions were determined by thermal gravimetry. Expressing the composition of calcium germanate hydrate as (CaO)_3–x GeO₂·nH₂O, the value of 3–x decreases from 1.68 to 1.59 with increasing temperature from 5–70 °C. The values obtained forn varied inconsistently between 2.4 and 3.3. The unit cell of the calcium germanate hydrate was determined to be monoclinic. Cell parameters werea=1.851,b=1.147,c=0.531 nm and =98.10 °.     

Producing persistent, high-current, high-duty-factor H- beams for routine 1 MW operation of Spallation Neutron Source (invited)

Since 2009, the Spallation Neutron Source (SNS) has been producing neutrons with ion beam powers near 1 MW, which requires the extraction of ～50 mA H(-) ions from the ion source with a ～5% duty factor. The 50 mA are achieved after an initial dose of ～3 mg of Cs and heating the Cs collar to ～170 °C. The 50 mA normally persist for the entire 4-week source service cycles. Fundamental processes are reviewed to elucidate the persistence of the SNS H(-) beams without a steady feed of Cs and why the Cs collar temperature may have to be kept near 170 °C.     

Total Travel Cost in Stochastic Assignment Models

Deterministic assignment models are sometimes used to approximate properties of more complex stochastic models. One property that is of particular interest from a system optimization viewpoint is total travel cost. This paper looks at the approximation of mean total travel cost. It is shown that deterministic models will underestimate this quantity in many common situations. Furthermore, discrepancies between total travel cost under the different modelling frameworks can lead to situations in which network modifications which are detrimental according to a stochastic model appear beneficial when using the natural deterministic approximation. We conclude that estimation of mean travel cost in stochastic assignment is often best done using simulation. Some suggestions are made regarding the implementation of traffic assignment simulation.     

Mobile phones as computing devices: the viruses are coming!

Understanding existing threats against mobile phones helps us better protect our information and prepare for future dangers. Mobile phones have evolved from their roots as analog walkie-talkies to full-scale Internet-enabled computers. Today, mobile phone handsets are arguably the dominant computer form factor consumers purchase. But having such powerful networked computers leads to a new class of malware: viruses, worms, and trojans specifically designed for the mobile environment. Merely enumerating known attacks doesn't help us understand all possible threats. Instead, we present a taxonomy of attacks against mobile phones that shows known as well as potential attacks.     

Microscopy under pressure--an optical chamber system for fluorescence microscopic analysis of living cells under high hydrostatic pressure

High hydrostatic pressure (HHP) becomes more and more interesting for life science research, since it can be employed to inactivate various cells. To directly monitor "cells under pressure," the development of an optical high-pressure chamber is required. Therefore, an optical pressure chamber that can be used for up to 300 MPa was constructed. This chamber has already been described as a tool for in situ observation of dynamic changes of microscopic structures in bright field as well as phase contrast. In combination with an inverted microscope, we obtained brilliant microscopic color pictures with an optical resolution more than 0.56 microm. Here, we demonstrate the capabilities of the HHP cell, in combination with epifluorescence microscopy. Using a nonadherent human B-cell line (Raji, ATCC CCL 86), stained with the fluorescent dyes propidium iodide, Hoechst 33342, or dihexyloxacarbocyanine iodide, we were able to show that the system is suitable to perform fluorescence microscopic analyses, with pressures up to 300 MPa, with viable mammalian cells.