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.

Fractional order and non-linear system identification algorithms for biomedical applications

We discuss the modelling of dielectric responses of amorphous biological samples. Such samples are commonly encountered in impedance spectroscopy studies as well as in UV, IR, optical and THz transient spectroscopy experiments and in pump-probe studies. In many occasions, the samples may display quenched absorption bands. A systems identification framework may be developed to provide parsimonious representations of such responses. To achieve this, it is appropriate to augment the standard models found in the identification literature to incorporate fractional order dynamics. Extensions of models using the forward shift operator, state space models as well as their non-linear Hammerstein-Wiener counterpart models are highlighted. We also discuss the need to extend the theory of electromagnetically excited networks which can account for fractional order behaviour in the non-linear regime by incorporating nonlinear elements to account for the observed non-linearities. The proposed approach leads to the development of a range of new chemometrics tools for biomedical data analysis and classification.

Ultrafast charge separation in a photoreactive rhenium-appended porphyrin assembly monitored by picosecond transient infrared spectroscopy

Rhenium(bipyridine)(tricarbonyl)(picoline) units have been linked covalently to tetraphenylmetalloporphyrins of magnesium and zinc via an amide bond between the bipyridine and one phenyl substituent of the porphyrin. The resulting complexes, abbreviated as [Re(CO)(3)(Pic)Bpy-MgTPP][OTf] and [Re(CO)(3)(Pic)Bpy-ZnTPP][OTf], exhibit no signs of electronic interaction between the Re(CO)(3)(bpy) units and the metalloporphyrin units in their ground states. However, emission spectroscopy reveals solvent-dependent quenching of porphyrin emission on irradiation into the long-wavelength absorption bands localized on the porphyrin. The characteristics of the excited states have been probed by picosecond time-resolved absorption (TRVIS) spectroscopy and time-resolved infrared (TRIR) spectroscopy in nitrile solvents. The presence of the charge-separated state involving electron transfer from MgTPP or ZnTPP to Re(bpy) is signaled in the TRIR spectra by a low-frequency shift in the nu(CO) bands of the Re(CO)(3) moiety similar to that observed by spectroelectrochemical reduction. Long-wavelength excitation of [Re(CO)(3)(Pic)Bpy-MTPP][OTf] results in characteristic TRVIS spectra of the S-1 state of the porphyrin that decay with a time constant of 17 ps (M = Mg) or 24 ps (M = Zn). The IR bands of the CS state appear on a time scale of less than 1 ps (Mg) or ca. 5 ps (Zn) and decay giving way to a vibrationally excited (i.e., hot) ground state via back electron transfer. The IR bands of the precursors recover with a time constant of 35 ps (Mg) or 55 ps (Zn). The short lifetimes of the charge-transfer states carry implications for the mechanism of reaction in the presence of triethylamine.

A narrow band pulsed eddy current instrument for materials testing

The principles of operation of an experimental prototype instrument known as J-SCAN are described along with the derivation of formulae for the rapid calculation of normalized impedances; the structure of the instrument; relevant probe design parameters; digital quantization errors; and approaches for the optimization of single frequency operation. An eddy current probe is used As the inductance element of a passive tuned-circuit which is repeatedly excited with short impulses. Each impulse excites an oscillation which is subject to decay dependent upon the values of the tuned-circuit components: resistance, inductance and capacitance. Changing conditions under the probe that affect the resistance and inductance of this circuit will thus be detected through changes in the transient response. These changes in transient response, oscillation frequency and rate of decay, are digitized, and then normalized values for probe resistance and inductance changes are calculated immediately in a micro processor. This approach coupled with a minimum analogue processing and maximum of digital processing has advantages compared with the conventional approaches to eddy current instruments. In particular there are: the absence of an out of balance condition and the flexibility and stability of digital data processing.

A multispacecraft analysis of a small-scale transient entrained by solar wind streams

The images taken by the Heliospheric Imagers (HIs), part of the SECCHI imaging package onboard the pair of STEREO spacecraft, provide information on the radial and latitudinal evolution of the plasma compressed inside corotating interaction regions (CIRs). A plasma density wave imaged by the HI instrument onboard STEREO-B was found to propagate towards STEREO-A, enabling a comparison between simultaneous remotesensing and in situ observations of its structure to be performed. In situ measurements made by STEREO-A show that the plasma density wave is associated with the passage of a CIR. The magnetic field compressed after the CIR stream interface (SI) is found to have a planar distribution. Minimum variance analysis of the magnetic field vectors shows that the SI is inclined at 54° to the orbital plane of the STEREO-A spacecraft. This inclination of the CIR SI is comparable to the inclination of the associated plasma density wave observed by HI. A small-scale magnetic cloud with a flux rope topology and radial extent of 0.08 AU is also embedded prior to the SI. The pitch-angle distribution of suprathermal electrons measured by the STEREO-A SWEA instrument shows that an open magnetic field topology in the cloud replaced the heliospheric current sheet locally. These observations confirm that HI observes CIRs in difference images when a small-scale transient is caught up in the compression region.

Comparison of subspace and ARX models of a waveguide's terahertz transient response after optimal wavelet filtering

A quasi-optical deembedding technique for characterizing waveguides is demonstrated using wide-band time-resolved terahertz spectroscopy. A transfer function representation is adopted for the description of the signal in the input and output port of the waveguides. The time-domain responses were discretized and the waveguide transfer function was obtained through a parametric approach in the z-domain after describing the system with an AutoRegressive with eXogenous input (ARX), as well as with a state-space model. Prior to the identification procedure, filtering was performed in the wavelet domain to minimize both signal distortion, as well as the noise propagating in the ARX and subspace models. The optimal filtering procedure used in the wavelet domain for the recorded time-domain signatures is described in detail. The effect of filtering prior to the identification procedures is elucidated with the aid of pole-zero diagrams. Models derived from measurements of terahertz transients in a precision WR-8 waveguide adjustable short are presented.

High signal to noise ratio THz spectroscopy with ASOPS and signal processing schemes for mapping and controlling molecular and bulk relaxation processes

Asynchronous Optical Sampling has the potential to improve signal to noise ratio in THz transient sperctrometry. The design of an inexpensive control scheme for synchronising two femtosecond pulse frequency comb generators at an offset frequency of 20 kHz is discussed. The suitability of a range of signal processing schemes adopted from the Systems Identification and Control Theory community for further processing recorded THz transients in the time and frequency domain are outlined. Finally, possibilities for femtosecond pulse shaping using genetic algorithms are mentioned.

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

Application of near and mid-infrared spectroscopy to determine cheese quality and authenticity

This paper reviews the current state of development of both near-infrared (NIR) and mid-infrared (MIR) spectroscopic techniques for process monitoring, quality control, and authenticity determination in cheese processing. Infrared spectroscopy has been identified as an ideal process analytical technology tool, and recent publications have demonstrated the potential of both NIR and MIR spectroscopy, coupled with chemometric techniques, for monitoring coagulation, syneresis, and ripening as well as determination of authenticity, composition, sensory, and rheological parameters. Recent research is reviewed and compared on the basis of experimental design, spectroscopic and chemometric methods employed to assess the potential of infrared spectroscopy as a technology for improving process control and quality in cheese manufacture. Emerging research areas for these technologies, such as cheese authenticity and food chain traceability, are also discussed.

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120kyr have been carried out using FAMOUS, a fast atmosphere-ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth's orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic Meridional Overturning Circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth's orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings.

Protease-activated receptor 2 sensitizes the transient receptor potential vanilloid 4 ion channel to cause mechanical hyperalgesia in mice

Exacerbated sensitivity to mechanical stimuli that are normally innocuous or mildly painful (mechanical allodynia and hyperalgesia) occurs during inflammation and underlies painful diseases. Proteases that are generated during inflammation and disease cleave protease-activated receptor 2 (PAR2) on afferent nerves to cause mechanical hyperalgesia in the skin and intestine by unknown mechanisms. We hypothesized that PAR2-mediated mechanical hyperalgesia requires sensitization of the ion channel transient receptor potential vanilloid 4 (TRPV4). Immunoreactive TRPV4 was coexpressed by rat dorsal root ganglia (DRG) neurons with PAR2, substance P (SP) and calcitonin gene-related peptide (CGRP), mediators of pain transmission. In PAR2-expressing cell lines that either naturally expressed TRPV4 (bronchial epithelial cells) or that were transfected to express TRPV4 (HEK cells), pretreatment with a PAR2 agonist enhanced Ca2+ and current responses to the TRPV4 agonists phorbol ester 4alpha-phorbol 12,13-didecanoate (4alphaPDD) and hypotonic solutions. PAR2-agonist similarly sensitized TRPV4 Ca2+ signals and currents in DRG neurons. Antagonists of phospholipase Cbeta and protein kinases A, C and D inhibited PAR2-induced sensitization of TRPV4 Ca2+ signals and currents. 4alphaPDD and hypotonic solutions stimulated SP and CGRP release from dorsal horn of rat spinal cord, and pretreatment with PAR2 agonist sensitized TRPV4-dependent peptide release. Intraplantar injection of PAR2 agonist caused mechanical hyperalgesia in mice and sensitized pain responses to the TRPV4 agonists 4alphaPDD and hypotonic solutions. Deletion of TRPV4 prevented PAR2 agonist-induced mechanical hyperalgesia and sensitization. This novel mechanism, by which PAR2 activates a second messenger to sensitize TRPV4-dependent release of nociceptive peptides and induce mechanical hyperalgesia, may underlie inflammatory hyperalgesia in diseases where proteases are activated and released.

Coherent EISCAT Svalbard Radar spectra from the dayside cusp/cleft and their implications for transient field-aligned currents

Naturally enhanced incoherent scatter spectra from the vicinity of the dayside cusp/cleft, interpreted as being due to plasma turbulence driven by short bursts of intense field-aligned current, are compared with high-resolution narrow-angle auroral images and meridian scanning photometer data. Enhanced spectra have been observed on many occasions in association with nightside aurora, but there has been only one report of such spectra seen in the cusp/cleft region. Narrow-angle images show considerable change in the aurora on timescales shorter than the 10-s radar integration period, which could explain spectra observed with both ion lines simultaneously enhanced. Enhanced radar spectra are generally seen inside or beside regions of 630-nm auroral emission, indicative of sharp F region conductivity gradients, but there appears also to be a correlation with dynamic, small-scale auroral forms of order 100 m and less in width.

Implications of the altitude of transient 630-nm dayside auroral emissions

The altitude from which transient 630-nm (“red line”) light is emitted in transient dayside auroral breakup events is discussed. Theoretically, the emissions should normally originate from approximately 250 to 550 km. Because the luminosity in dayside breakup events moves in a way that is consistent with newly opened field lines, they have been interpreted as the ionospheric signatures of transient reconnection at the dayside magnetopause. For this model the importance of these events for convection can be assessed from the rate of change of their area. The area derived from analysis of images from an all-sky camera and meridian scans from a photometer, however, depends on the square of the assumed emission altitude. From field line mapping, it is shown for both a westward and an eastward moving event, that the main 557.7-nm emission comes from the edge of the 630 nm transient, where a flux transfer event model would place the upward field-aligned current (on the poleward and equatorward edge, respectively). The observing geometry for the two cases presented is such that this is true, irrespective of the 630-nm emission altitude. From comparisons with the European incoherent scatter radar data for the westward (interplanetary magnetic field By > 0) event on January 12, 1988, the 630-nm emission appears to emanate from an altitude of 250 km, and to be accompanied by some 557.7-nm “green-line” emission. However, for a large, eastward moving event observed on January 9, 1989, there is evidence that the emission altitude is considerably greater and, in this case, the only 557.7-nm emission is that on the equatorward edge of the event, consistent with a higher altitude 630-nm excitation source. Assuming an emission altitude of 250 km for this event yields a reconnection voltage of >50 kV during the reconnection burst but a contribution to the convection voltage of >15 kV. However, from the motion of the event we infer that the luminosity peaks at an altitude in the range of 400 and 500 km, and for the top of this range the reconnection and average convection voltages would be increased to >200 kV and >60 kV, respectively. (These are all minimum estimates because the event extends in longitude beyond the field-of-view of the camera). Hence the higher-emission altitude has a highly significant implication, namely that the reconnection bursts which cause the dayside breakup events could explain most of the voltage placed across the magnetosphere and polar cap by the solar wind flow. Analysis of the plasma density and temperatures during the event on January 9, 1989, predicts the required thermal excitation of significant 630-nm intensities at altitudes of 400-500 km.

The water vapour continuum in near-infrared windows – current understanding and prospects for its inclusion in spectroscopic databases

Spectroscopic catalogues, such as GEISA and HITRAN, do not yet include information on the water vapour continuum that pervades visible, infrared and microwave spectral regions. This is partly because, in some spectral regions, there are rather few laboratory measurements in conditions close to those in the Earth’s atmosphere; hence understanding of the characteristics of the continuum absorption is still emerging. This is particularly so in the near-infrared and visible, where there has been renewed interest and activity in recent years. In this paper we present a critical review focusing on recent laboratory measurements in two near-infrared window regions (centred on 4700 and 6300 cm−1) and include reference to the window centred on 2600 cm−1 where more measurements have been reported. The rather few available measurements, have used Fourier transform spectroscopy (FTS), cavity ring down spectroscopy, optical-feedback – cavity enhanced laser spectroscopy and, in very narrow regions, calorimetric interferometry. These systems have different advantages and disadvantages. Fourier Transform Spectroscopy can measure the continuum across both these and neighbouring windows; by contrast, the cavity laser techniques are limited to fewer wavenumbers, but have a much higher inherent sensitivity. The available results present a diverse view of the characteristics of continuum absorption, with differences in continuum strength exceeding a factor of 10 in the cores of these windows. In individual windows, the temperature dependence of the water vapour self-continuum differs significantly in the few sets of measurements that allow an analysis. The available data also indicate that the temperature dependence differs significantly between different near-infrared windows. These pioneering measurements provide an impetus for further measurements. Improvements and/or extensions in existing techniques would aid progress to a full characterisation of the continuum – as an example, we report pilot measurements of the water vapour self-continuum using a supercontinuum laser source coupled to an FTS. Such improvements, as well as additional measurements and analyses in other laboratories, would enable the inclusion of the water vapour continuum in future spectroscopic databases, and therefore allow for a more reliable forward modelling of the radiative properties of the atmosphere. It would also allow a more confident assessment of different theoretical descriptions of the underlying cause or causes of continuum absorption.