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.

Optimization of apodization functions in terahertz transient spectrometry

A quadratic programming optimization procedure for designing asymmetric apodization windows tailored to the shape of time-domain sample waveforms recorded using a terahertz transient spectrometer is proposed. By artificially degrading the waveforms, the performance of the designed window in both the time and the frequency domains is compared with that of conventional rectangular, triangular (Mertz), and Hamming windows. Examples of window optimization assuming Gaussian functions as the building elements of the apodization window are provided. The formulation is sufficiently general to accommodate other basis functions. (C) 2007 Optical Society of America

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.

Contribution of the double layer to transient faradaic processes: implications for hydrodynamic modulated voltammetry of nanostructures

The electrochemistry of Pt nanostructured electrodes is investigated using hydrodynamic modulated voltammetry (HMV). Here a liquid crystal templating process is used to produce platinum-modified electrodes with a range of surface areas (roughness factor 42.4-280.8). The electroreduction of molecular oxygen at these nanostructured platinum surfaces is used to demonstrate the ability of HMV to discriminate between faradaic and nonfaradaic electrode reactions. The HMV approach shows that the reduction of molecular oxygen experiences considerable signal loss within the high pseudocapacitive region of the voltammetry. Evidence for the contribution of the double layer to transient mass transfer events is presented. In addition, a model circuit and appropriate theoretical analysis are used to illustrate the transient responses of a time variant faradaic component. This in conjunction with the experimental evidence shows that, far from being a passive component in this system, the double layer can contribute to HMV faradaic reactions under certain conditions.

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.

Energetics of climate models: net energy balance and meridional enthalpy transport

We analyze the publicly released outputs of the simulations performed by climate models (CMs) in preindustrial (PI) and Special Report on Emissions Scenarios A1B (SRESA1B) conditions. In the PI simulations, most CMs feature biases of the order of 1 W m −2 for the net global and the net atmospheric, oceanic, and land energy balances. This does not result from transient effects but depends on the imperfect closure of the energy cycle in the fluid components and on inconsistencies over land. Thus, the planetary emission temperature is underestimated, which may explain the CMs' cold bias. In the PI scenario, CMs agree on the meridional atmospheric enthalpy transport's peak location (around 40°N/S), while discrepancies of ∼20% exist on the intensity. Disagreements on the oceanic transport peaks' location and intensity amount to ∼10° and ∼50%, respectively. In the SRESA1B runs, the atmospheric transport's peak shifts poleward, and its intensity increases up to ∼10% in both hemispheres. In most CMs, the Northern Hemispheric oceanic transport decreases, and the peaks shift equatorward in both hemispheres. The Bjerknes compensation mechanism is active both on climatological and interannual time scales. The total meridional transport peaks around 35° in both hemispheres and scenarios, whereas disagreements on the intensity reach ∼20%. With increased CO 2 concentration, the total transport increases up to ∼10%, thus contributing to polar amplification of global warming. Advances are needed for achieving a self-consistent representation of climate as a nonequilibrium thermodynamical system. This is crucial for improving the CMs' skill in representing past and future climate changes.

Reactions of NO3 with the man-made emissions 2-methylpent-2-ene, (Z)-3-methylpent-2-ene, ethyl vinyl ether, and the stress-induced plant emission ethyl vinyl ketone

Rate coefficients for reactions of nitrate radicals (NO3) with the anthropogenic emissions 2-methylpent-2-ene, (Z)-3-methylpent-2-ene.. ethyl vinyl ether, and the stress-induced plant emission ethyl vinyl ketone (pent-1-en-3-one) were determined to be (9.3 +/- 1.1) x 10(-12), (9.3 +/- 3.2) x 10(-12), (1.7 +/- 1.3) x 10(-12) and (9.4 + 2.7) x 10(-17) cm(3) molecule(-1) s(-1). We performed kinetic experiments at room temperature and atmospheric pressure using a relative-rate technique with GC-FID analysis. Experiments with ethyl vinyl ether required a modification of our established procedure that might introduce additional uncertainties, and the errors suggested reflect these difficulties. Rate coefficients are discussed in terms of electronic and steric influences. Atmospheric lifetimes with respect to important oxidants in the troposphere were calculated. NO3-initiated oxidation is found to be the strongly dominating degradation route for 2-methylpent-2-ene, (Z)-3-methylpent-2-ene and ethyl vinyl ether. Atmospheric concentrations of the alkenes and their relative contribution to the total NMHC emissions from trucks can be expected to increase if plans for the introduction of particle filters for diesel engines are implemented on a global scale. Thus more kinetic data are required to better evaluate the impact of these emissions.

The impact of traffic emissions on atmospheric ozone and OH: results from QUANTIFY

To estimate the impact of emissions by road, aircraft and ship traffic on ozone and OH in the present-day atmosphere six different atmospheric chemistry models have been used. Based on newly developed global emission inventories for road, ship and aircraft emission data sets each model performed sensitivity simulations reducing the emissions of each transport sector by 5%. The model results indicate that on global annual average lower tropospheric ozone responds most sensitive to ship emissions (50.6%±10.9% of the total traffic induced perturbation), followed by road (36.7%±9.3%) and aircraft exhausts (12.7%±2.9%), respectively. In the northern upper troposphere between 200–300 hPa at 30–60° N the maximum impact from road and ship are 93% and 73% of the maximum effect of aircraft, respectively. The latter is 0.185 ppbv for ozone (for the 5% case) or 3.69 ppbv when scaling to 100%. On the global average the impact of road even dominates in the UTLS-region. The sensitivity of ozone formation per NOx molecule emitted is highest for aircraft exhausts. The local maximum effect of the summed traffic emissions on the ozone column predicted by the models is 0.2 DU and occurs over the northern subtropical Atlantic extending to central Europe. Below 800 hPa both ozone and OH respond most sensitively to ship emissions in the marine lower troposphere over the Atlantic. Based on the 5% perturbation the effect on ozone can exceed 0.6% close to the marine surface (global zonal mean) which is 80% of the total traffic induced ozone perturbation. In the southern hemisphere ship emissions contribute relatively strongly to the total ozone perturbation by 60%–80% throughout the year. Methane lifetime changes against OH are affected strongest by ship emissions up to 0.21 (± 0.05)%, followed by road (0.08 (±0.01)%) and air traffic (0.05 (± 0.02)%). Based on the full scale ozone and methane perturbations positive radiative forcings were calculated for road emissions (7.3±6.2 mWm−2) and for aviation (2.9±2.3 mWm−2). Ship induced methane lifetime changes dominate over the ozone forcing and therefore lead to a net negative forcing (−25.5±13.2 mWm−2).

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.

Night-time tropospheric chemistry of the unsaturated alcohols (Z)-pent-2-en-1-ol and pent-1-en-3-ol: Kinetic studies of reactions of NO3 and N2O5 with stress-induced plant emissions

The night-time tropospheric chemistry of two stress-induced volatile organic compounds (VOCs), (Z)-pent-2-en-1-ol and pent-1-en-3-ol, has been studied at room temperature. Rate coefficients for reactions of the nitrate radical (NO3) with these pentenols were measured using the discharge-flow technique. Because of the relatively low volatility of these compounds, we employed off-axis continuous-wave cavity-enhanced absorption spectroscopy for detection of NO3 in order to be able to work in pseudo first-order conditions with the pentenols in large excess over NO3. The rate coefficients were determined to be (1.53 +/- 0.23) x 10(-13) and (1.39 +/- 0.19) x 10(-14) cm(3) molecule(-1) s(-1) for reactions of NO3 with (Z)-pent-2-en-1-ol and pent-1-en-3-ol. An attempt to study the kinetics of these reactions with a relative-rate technique, using N2O5 as source of NO3 resulted in significantly higher apparent rate coefficients. Performing relative-rate experiments in known excesses of NO2 allowed us to determine the rate coefficients for the N2O5 reactions to be (5.0 +/- 2.8) x 10(-19) cm(3) molecule(-1) s(-1) for (Z)-pent-2-en-1-ol, and (9.1 +/- 5.8) x 10(-19) cm(3) molecule(-1) s(-1) for pent-1-en-3-ol. We show that these relatively slow reactions can indeed interfere with rate determinations in conventional relative-rate experiments.

A potential lag between the open solar magnetic source flux and solar EUV and X-ray emissions as measured by the Earth's ionosphere during total solar eclipses

Measurements of the ionospheric E-region during total solar eclipses have been used to provide information about the evolution of the solar magnetic field and EUV and X-ray emissions from the solar corona and chromosphere. By measuring levels of ionisation during an eclipse and comparing these measurements with an estimate of the unperturbed ionisation levels (such as those made during a control day, where available) it is possible to estimate the percentage of ionising radiation being emitted by the solar corona and chromosphere. Previously unpublished data from the two eclipses presented here are particularly valuable as they provide information that supplements the data published to date. The eclipse of 23 October 1976 over Australia provides information in a data gap that would otherwise have spanned the years 1966 to 1991. The eclipse of 4 December 2002 over Southern Africa is important as it extends the published sequence of measurements. Comparing measurements from eclipses between 1932 and 2002 with the solar magnetic source flux reveals that changes in the solar EUV and X-ray flux lag the open source flux measurements by approximately 1.5 years. We suggest that this unexpected result comes about from changes to the relative size of the limb corona between eclipses, with the lag representing the time taken to populate the coronal field with plasma hot enough to emit the EUV and X-rays ionising our atmosphere.

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.