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.

Argon predissociation spectroscopy of the OH-center dot H2O and Cl-center dot H2O complexes in the 1000-1900 cm(-1) region: Intramolecular bending transitions and the search for the shared-proton fundamental in the hydroxide monohydrate

We present argon predissociation vibrational spectra of the OH-.H2O and Cl-.H2O complexes in the 1000-1900 cm(-1) energy range, far below the OH stretching region reported in previous studies. This extension allows us to explore the fundamental transitions of the intramolecular bending vibrations associated with the water molecule, as well as that of the shared proton inferred from previous assignments of overtones in the higher energy region. Although the water bending fundamental in the Cl-.H2O spectrum is in very good agreement with expectations, the OH-.H2O spectrum is quite different than anticipated, being dominated by a strong feature at 1090 cm(-1). New full-diniensionality calculations of the OH-.H2O vibrational level structure using diffusion Monte Carlo and the VSCF/CI methods indicate this band arises from excitation of the shared proton.

The gas-phase reactions Of SiCl4 and Si2Cl6 with CH3OH and C2H5OH: an investigation by mass spectrometry and matrix-isolation infrared spectroscopy

The gas phase reactions Of SiCl4 and Si2Cl6 With CH3OH and C2H5OH have been investigated using both mass spectrometry and matrix isolation techniques. SiCl4 reacts with both CH3OH and C2H5OH upon mixing of the vapours for times in excess of 3 h to generate the HCl-elimination products SiCl3OR (R = CH3 or C2H5). The identity of these products is confirmed by deuteration experiments and by ab initio calculations at the HF/6-31G(d) level. Further products are generated when the mixture is passed through a tube heated to 750degreesC. Si2Cl6 reacts with CH3OH and C2H5OH via a different mechanism in which the Si-Si bond is cleaved to yield SiCl3OR and HCl. Other products of the type SiCl4-n(OCH3)(n) are tentatively identified by a combination of mass spectrometric and matrix isolation measurements. These latter products indicate further replacement of Cl atoms by OR groups as a result of reaction of CH3OH or C2H5OH with the initial product.

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.

Protein kinase C down-regulation, and not transient activation, correlates with melanocyte growth

The nontumorigenic, immortal line of murine melanocytes, Mel-ab, requires the continual presence of biologically active phorbol esters for growth (R. E. Wilson et al., Cancer Res., 49: 711–716, 1989). Comparable treatments of B16 murine melanoma cells result in partial inhibition of cell proliferation. The role of protein kinase C (PKC) in the modulation of growth of cells from these two melanocytic cell lines has been investigated. Significant levels of PKC were present in quiescent Mel-ab cells as determined by Western blotting, whereas no immunoreactive protein was detected in cell extracts from either proliferating Mel-ab or B16.F1 cells. Phosphorylation of a Mr 80,000 protein, which by one- and two-dimensional gel analysis comigrated with the known Mr 80,000 protein substrate of PKC in fibroblasts, was induced in 12-O-tetradecanoylphorbol-13-acetate-stimulated quiescent Mel-ab cells but not in proliferating Mel-ab cells or B16.F1 melanoma cells. Direct measurement of PKC activity in these cells demonstrated a 10-fold greater level of activity in quiescent Mel-ab cells (262 ± 50 pmol/min/mg SD) compared with growing cells (22.8 ± 11.8 pmol/min/mg SD). An intermediate level of activity was detected in proliferating B16.F1 melanoma cells (148.5 ± 20.4 pmol/min/mg SD). The subcellular distribution of PKC was dependent upon the growth state of the cells such that quiescent Mel-ab cells displayed a higher level of activity in the cytosol, whereas growing Melab cells displayed greater activity in the particulate fraction. Like many other transformed lines, B16.F1 melanoma cells constitutively expressed the majority of enzyme activity in the particulate fraction. Measurement of [3H]phorbol ester binding in intact cells paralleled the PKC activation data such that quiescent Mel-ab cells displayed binding of 1612 ± 147 cpm/106 cells, whereas proliferating Mel-ab and B16.F1 melanoma cells displayed binding of 652 ± 28 and 947 ± 81 cpm/106 cells, respectively. Membrane-permeant diacylglycerol analogues, which activated but did not down-regulate PKC, were devoid of growth-stimulating effects on melanocytes, even in the presence of the specific diacylglycerol kinase inhibitor, R59022. Together, these data show that PKC down-regulation, and not activation, correlates with the growth of melanocytes in culture.

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.

Deposition of charged inhaled aerosols with transient airflow in sequential lung airway model

Transport and deposition of charged inhaled aerosols in double planar bifurcation representing generation three to five of human respiratory system has been studied under a light activity breathing condition. Both steady and oscillatory laminar inhalation airflow is considered. Particle trajectories are calculated using a Lagrangian reference frame, which is dominated by the fluid force driven by airflow, gravity force and electrostatic forces (both of space and image charge forces). The particle-mesh method is selected to calculate the space charge force. This numerical study investigates the deposition efficiency in the three-dimensional model under various particle sizes, charge values, and inlet particle distribution. Numerical results indicate that particles carrying an adequate level of charge can improve deposition efficiency in the airway model.

Ocean heat uptake and its consequences for the magnitude of sea level rise and climate change

Under increasing greenhouse gas concentrations, ocean heat uptake moderates the rate of climate change, and thermal expansion makes a substantial contribution to sea level rise. In this paper we quantify the differences in projections among atmosphere-ocean general circulation models of the Coupled Model Intercomparison Project in terms of transient climate response, ocean heat uptake efficiency and expansion efficiency of heat. The CMIP3 and CMIP5 ensembles have statistically indistinguishable distributions in these parameters. The ocean heat uptake efficiency varies by a factor of two across the models, explaining about 50% of the spread in ocean heat uptake in CMIP5 models with CO2 increasing at 1%/year. It correlates with the ocean global-mean vertical profiles both of temperature and of temperature change, and comparison with observations suggests the models may overestimate ocean heat uptake and underestimate surface warming, because their stratification is too weak. The models agree on the location of maxima of shallow ocean heat uptake (above 700 m) in the Southern Ocean and the North Atlantic, and on deep ocean heat uptake (below 2000 m) in areas of the Southern Ocean, in some places amounting to 40% of the top-to-bottom integral in the CMIP3 SRES A1B scenario. The Southern Ocean dominates global ocean heat uptake; consequently the eddy-induced thickness diffusivity parameter, which is particularly influential in the Southern Ocean, correlates with the ocean heat uptake efficiency. The thermal expansion produced by ocean heat uptake is 0.12 m YJ−1, with an uncertainty of about 10% (1 YJ = 1024 J).

Assessment of winter cyclone activity in a transient ECHAM4-OPYC3 GHG experiment

Winter cyclone activity over the Northern Hemisphere is investigated in an ECHAM4/OPYC3 greenhouse gas scenario simulation. The goal of this investigation is to identify changes in cyclone activity associated with increasing concentrations. To this aim, two 50-year time periods are analysed, one representing present day climate conditions and the other a perturbed climate when CO2 concentrations exceed twice the present concentrations. Cyclone activity is assessed using an automatic algorithm, which identifies and tracks cyclones based on sea level pressure fields. The algorithm detects not only large and long living cyclones over the main ocean basins, but also their smaller counterparts in secondary storm track regions like the Mediterranean Basin. For the present climate, results show a good agreement with NCEP-reanalysis, provided that the spectral and time resolutions of the reanalysis are reduced to those available for the model. Several prominent changes in cyclone activity are observed for the scenario period in comparison to the present day climate, especially over the main ocean basins. A significant decrease of overall cyclone track density is found between 35 and 55 degrees North, together with a small increase polewards. These changes result from two different signals for deep and medium cyclones: for deep cyclones (core pressure below 990 hPa) there is a poleward shift in the greenhouse gas scenario, while for medium cyclones (core pressure between 990 and 1010 hPa) a general decrease in cyclone counts is found. The same kind of changes (a shift for intense cyclones and an overall decrease for the weaker ones) are detected when distinguishing cyclones from their intensity, quantified in terms of ∇2p. Thus, the simulated changes can not solely be attributed to alterations in mean sea level pressure. Instead, corresponding increases in upper-tropospheric baroclinicity suggest more favourable conditions for the development of stronger systems at higher latitudes, especially at the delta regions of the North Atlantic and the North Pacific storm tracks.

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers. Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling

Understanding neurovascular coupling is a prerequisite for the interpretation of results obtained from modern neuroimaging techniques. This study investigated the hemodynamic and neural responses in rat somatosensory cortex elicited by 16 seconds electrical whisker stimuli. Hemodynamics were measured by optical imaging spectroscopy and neural activity by multichannel electrophysiology. Previous studies have suggested that the whisker-evoked hemodynamic response contains two mechanisms, a transient ‘backwards’ dilation of the middle cerebral artery, followed by an increase in blood volume localized to the site of neural activity. To distinguish between the mechanisms responsible for these aspects of the response, we presented whisker stimuli during normocapnia (‘control’), and during a high level of hypercapnia. Hypercapnia was used to ‘predilate’ arteries and thus possibly ‘inhibit’ aspects of the response related to the ‘early’ mechanism. Indeed, hemodynamic data suggested that the transient stimulus-evoked response was absent under hypercapnia. However, evoked neural responses were also altered during hypercapnia and convolution of the neural responses from both the normocapnic and hypercapnic conditions with a canonical impulse response function, suggested that neurovascular coupling was similar in both conditions. Although data did not clearly dissociate early and late vascular responses, they suggest that the neurovascular coupling relationship is neurogenic in origin.