Subspace system identification framework for the analysis of multimoded propagation of THz-transient signals

We provide a system identification framework for the analysis of THz-transient data. The subspace identification algorithm for both deterministic and stochastic systems is used to model the time-domain responses of structures under broadband excitation. Structures with additional time delays can be modelled within the state-space framework using additional state variables. We compare the numerical stability of the commonly used least-squares ARX models to that of the subspace N4SID algorithm by using examples of fourth-order and eighth-order systems under pulse and chirp excitation conditions. These models correspond to structures having two and four modes simultaneously propagating respectively. We show that chirp excitation combined with the subspace identification algorithm can provide a better identification of the underlying mode dynamics than the ARX model does as the complexity of the system increases. The use of an identified state-space model for mode demixing, upon transformation to a decoupled realization form is illustrated. Applications of state-space models and the N4SID algorithm to THz transient spectroscopy as well as to optical systems are highlighted.

Wigner Transport models of the electron-phonon kinetics in quantum wires

Two quantum-kinetic models of ultrafast electron transport in quantum wires are derived from the generalized electron-phonon Wigner equation. The various assumptions and approximations allowing one to find closed equations for the reduced electron Wigner function are discussed with an emphasis on their physical relevance. The models correspond to the Levinson and Barker-Ferry equations, now generalized to account for a space-dependent evolution. They are applied to study the quantum effects in the dynamics of an initial packet of highly nonequilibrium carriers, locally generated in the wire. The properties of the two model equations are compared and analyzed.

Global asymptotic stability in a rational recursive sequence

In this paper, we study the global stability of the difference equation x(n) = a + bx(n-1) + cx(n-1)(2)/d - x(n-2), n = 1,2,....., where a, b greater than or equal to 0 and c, d > 0. We show that one nonnegative equilibrium point of the equation is a global attractor with a basin that is determined by the parameters, and every positive Solution of the equation in the basin exponentially converges to the attractor. (C) 2003 Elsevier Inc. All rights reserved.

The demand for a healthy diet: estimating the almost ideal demand system with infrequency of purchase

A Bayesian method of estimating multivariate sample selection models is introduced and applied to the estimation of a demand system for food in the UK to account for censoring arising from infrequency of purchase. We show how it is possible to impose identifying restrictions on the sample selection equations and that, unlike a maximum likelihood framework, the imposition of adding up at both latent and observed levels is straightforward. Our results emphasise the role played by low incomes and socio-economic circumstances in leading to poor diets and also indicate that the presence of children in a household has a negative impact on dietary quality.

Neurofuzzy control using Kalman filtering state feedback with coloured noise for unknown non-linear processes

This paper presents a controller design scheme for a priori unknown non-linear dynamical processes that are identified via an operating point neurofuzzy system from process data. Based on a neurofuzzy design and model construction algorithm (NeuDec) for a non-linear dynamical process, a neurofuzzy state-space model of controllable form is initially constructed. The control scheme based on closed-loop pole assignment is then utilized to ensure the time invariance and linearization of the state equations so that the system stability can be guaranteed under some mild assumptions, even in the presence of modelling error. The proposed approach requires a known state vector for the application of pole assignment state feedback. For this purpose, a generalized Kalman filtering algorithm with coloured noise is developed on the basis of the neurofuzzy state-space model to obtain an optimal state vector estimation. The derived controller is applied in typical output tracking problems by minimizing the tracking error. Simulation examples are included to demonstrate the operation and effectiveness of the new approach.

Changes in the intestinal microbiota after a short period of dietary over-indulgence, representative of a holiday or festival season

The effects on the intestinal microbiota of a short period of marginal over-eating, characteristic of holiday or festival periods, were investigated in a pilot study. Fourteen healthy male subjects consumed a diet rich in animal protein and fat for seven days. During this period, the subjects significantly increased their dietary energy, protein, carbohydrate and fat intakes by 56, 59, 53 and 58%, respectively (all P < 0.05). The mean weight gain of 0.27 kg was less than the expected 1 kg, but this was consistent with a degree of under-reporting on the baseline diet. Fluorescence in situ hybridisation analysis confirmed the relative stability of each individual’s faecal microbiota but showed considerable variations between them. The diet was associated with a significant increase in numbers of total faecal bacteria and the bacteroides group, as detected by the universal bacterial probe (DAPI) and Bacteroides probe (Bac 303), respectively. Overall, there was a decrease in numbers of the Lactobacillus/Enterococcus group (Lab 158 probe; 2.8 ± 3.0% to 1.8 ± 1.8%) and the Bifidobacterium group (Bif 164 probe; 3.0 ± 3.7% to 1.7 ± 1.2%), although there was considerable inter-individual variation. Analysis of the relative proportions of each bacterial group as a percentage of the subject’s total bacteria showed a trend for a change in the intestinal microbiota that might be considered potentially unhealthy.

Anisotropic refinement of the structure of Thermoascus aurantiacusxylanase I

The isotropic crystallographic model of the structure of xylanase I from Thermoascus aurantiacus (TAXI) has now been refined anisotropically at 1.14 Å resolution to a standard residual of R = 11.1% for all data. TAXI is amongst the five largest proteins deposited in the Protein Data Bank to have been refined with anisotropic displacement parameters (ADPs) at this level of resolution. The anisotropy analysis revealed a more isotropic distribution of anisotropy than usually observed previously. Adding ADPs resulted in high-quality electron-density maps which revealed discrepancies from the previously suggested primary sequences for this enzyme. Side-chain conformational disorder was modelled for 16 residues, including Trp275, a bulky residue at the active site. An unrestrained refinement was consistent with the protonation of the catalytic acid/base glutamate and the deprotonation of the nucleophile glutamate, as required for catalysis. The thermal stability of TAXI is reinterpreted in the light of the new refined model.

Efficiency of pseudo-spectral algorithms with Anderson mixing for the SCFT of periodic block-copolymer phases

This study examines the numerical accuracy, computational cost, and memory requirements of self-consistent field theory (SCFT) calculations when the diffusion equations are solved with various pseudo-spectral methods and the mean field equations are iterated with Anderson mixing. The different methods are tested on the triply-periodic gyroid and spherical phases of a diblock-copolymer melt over a range of intermediate segregations. Anderson mixing is found to be somewhat less effective than when combined with the full-spectral method, but it nevertheless functions admirably well provided that a large number of histories is used. Of the different pseudo-spectral algorithms, the 4th-order one of Ranjan, Qin and Morse performs best, although not quite as efficiently as the full-spectral method.

Comment on “A modified leapfrog scheme for shallow water equations” by Wen-Yih Sun and Oliver M.T. Sun

A recent paper published in this journal considers the numerical integration of the shallow-water equations using the leapfrog time-stepping scheme [Sun Wen-Yih, Sun Oliver MT. A modified leapfrog scheme for shallow water equations. Comput Fluids 2011;52:69–72]. The authors of that paper propose using the time-averaged height in the numerical calculation of the pressure-gradient force, instead of the instantaneous height at the middle time step. The authors show that this modification doubles the maximum Courant number (and hence the maximum time step) at which the integrations are stable, doubling the computational efficiency. Unfortunately, the pressure-averaging technique proposed by the authors is not original. It was devised and published by Shuman [5] and has been widely used in the atmosphere and ocean modelling community for over 40 years.

A global view of the extratropical tropopause transition layer from Atmospheric Chemistry Experiment Fourier Transform Spectrometer O3, H2O, and CO

The global behavior of the extratropical tropopause transition layer (ExTL) is investigated using O3, H2O, and CO measurements from the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) on Canada’s SCISAT-1 satellite obtained between February 2004 and May 2007. The ExTL depth is derived using H2O-O3 and CO-O3 correlations. The ExTL top derived from H2O-O3 shows an increase from roughly 1–1.5 km above the thermal tropopause in the subtropics to 3–4 km (2.5–3.5 km) in the north (south) polar region, implying somewhat weaker tropospherestratosphere- transport in the Southern Hemisphere. The ExTL bottom extends ~1 km below the thermal tropopause, indicating a persistent stratospheric influence on the troposphere at all latitudes. The ExTL top derived from the CO-O3 correlation is lower, at 2 km or ~345 K (1.5 km or ~335 K) in the Northern (Southern) Hemisphere. Its annual mean coincides with the relative temperature maximum just above the thermal tropopause. The vertical CO gradient maximizes at the thermal tropopause, indicating a local minimum in mixing within the tropopause region. The seasonal changes in and the scales of the vertical H2O gradients show a similar pattern as the static stability structure of the tropopause inversion layer (TIL), which provides observational support for the hypothesis that H2O plays a radiative role in forcing and maintaining the structure of the TIL.

Tuning the solubilities of bis-triazinylphenanthroline ligands (BTPhens) and their complexes

A series of bis-triazinylphenanthroline ligands (BTPhens) was synthesized by modifying the triazine substituents. It was found that varying these substituents altered the solubilities of the ligands in a number of non-polar solvents. Thus C5-BTPhen showed significantly higher solubility in octanol than C1-BTPhen. The high solubility of C5-BTPhen and its complexes was exploited to facilitate the NMR titration experiments. These experiments shown that the dominant species in solution were the 1:2 complexes [Ln(III)(BTPhen)2], even at high Ln concentrations, and that the relative stability of the 2:1 to 1:1 BTPhen-Ln complexes varied with different lanthanides. C5-BTPhen therefore shows considerable promise for a once-through selective actinide separation process.

Correlations between two sets of angular relation equations

It is shown here that the angular relation equations between direct and reciprocal vectors are very similar to the angular relation equations in Euler's theorem. These two sets of equations are usually treated separately as unrelated equations in different fields. In this careful study, the connection between the two sets of angular equations is revealed by considering the cosine rule for the spherical triangle. It is found that understanding of the correlation is hindered by the facts that the same variables are defined differently and different symbols are used to represent them in the two fields. Understanding the connection between different concepts is not only stimulating and beneficial, but also a fundamental tool in innovation and research, and has historical significance. The background of the work presented here contains elements of many scientific disciplines. This work illustrates the common ground of two theories usually considered separately and is therefore of benefit not only for its own sake but also to illustrate a general principle that a theory relevant to one discipline can often be used in another. The paper works with chemistry related concepts using mathematical methodologies unfamiliar to the usual audience of mainstream experimental and theoretical chemists.

Commodity futures prices: more evidence on forecast power, risk premia and the theory of storage

In this paper, we examine the temporal stability of the evidence for two commodity futures pricing theories. We investigate whether the forecast power of commodity futures can be attributed to the extent to which they exhibit seasonality and we also consider whether there are time varying parameters or structural breaks in these pricing relationships. Compared to previous studies, we find stronger evidence of seasonality in the basis, which supports the theory of storage. The power of the basis to forecast subsequent price changes is also strengthened, while results on the presence of a risk premium are inconclusive. In addition, we show that the forecasting power of commodity futures cannot be attributed to the extent to which they exhibit seasonality. We find that in most cases where structural breaks occur, only changes in the intercepts and not the slopes are detected, illustrating that the forecast power of the basis is stable over different economic environments.

A unified theory of balance in the extratropics

Many physical systems exhibit dynamics with vastly different time scales. Often the different motions interact only weakly and the slow dynamics is naturally constrained to a subspace of phase space, in the vicinity of a slow manifold. In geophysical fluid dynamics this reduction in phase space is called balance. Classically, balance is understood by way of the Rossby number R or the Froude number F; either R ≪ 1 or F ≪ 1. We examined the shallow-water equations and Boussinesq equations on an f -plane and determined a dimensionless parameter _, small values of which imply a time-scale separation. In terms of R and F, ∈= RF/√(R^2+R^2 ) We then developed a unified theory of (extratropical) balance based on _ that includes all cases of small R and/or small F. The leading-order systems are ensured to be Hamiltonian and turn out to be governed by the quasi-geostrophic potential-vorticity equation. However, the height field is not necessarily in geostrophic balance, so the leading-order dynamics are more general than in quasi-geostrophy. Thus the quasi-geostrophic potential-vorticity equation (as distinct from the quasi-geostrophic dynamics) is valid more generally than its traditional derivation would suggest. In the case of the Boussinesq equations, we have found that balanced dynamics generally implies hydrostatic balance without any assumption on the aspect ratio; only when the Froude number is not small and it is the Rossby number that guarantees a timescale separation must we impose the requirement of a small aspect ratio to ensure hydrostatic balance.

On wave action and phase in the non-canonical Hamiltonian formulation

The long time–evolution of disturbances to slowly–varying solutions of partial differential equations is subject to the adiabatic invariance of the wave action. Generally, this approximate conservation law is obtained under the assumption that the partial differential equations are derived from a variational principle or have a canonical Hamiltonian structure. Here, the wave action conservation is examined for equations that possess a non–canonical (Poisson) Hamiltonian structure. The linear evolution of disturbances in the form of slowly varying wavetrains is studied using a WKB expansion. The properties of the original Hamiltonian system strongly constrain the linear equations that are derived, and this is shown to lead to the adiabatic invariance of a wave action. The connection between this (approximate) invariance and the (exact) conservation laws of pseudo–energy and pseudomomentum that exist when the basic solution is exactly time and space independent is discussed. An evolution equation for the slowly varying phase of the wavetrain is also derived and related to Berry's phase.