Classification of lactose and mandelic acid THz spectra using subspace and wavelet-packet algorithms

This work compares classification results of lactose, mandelic acid and dl-mandelic acid, obtained on the basis of their respective THz transients. The performance of three different pre-processing algorithms applied to the time-domain signatures obtained using a THz-transient spectrometer are contrasted by evaluating the classifier performance. A range of amplitudes of zero-mean white Gaussian noise are used to artificially degrade the signal-to-noise ratio of the time-domain signatures to generate the data sets that are presented to the classifier for both learning and validation purposes. This gradual degradation of interferograms by increasing the noise level is equivalent to performing measurements assuming a reduced integration time. Three signal processing algorithms were adopted for the evaluation of the complex insertion loss function of the samples under study; a) standard evaluation by ratioing the sample with the background spectra, b) a subspace identification algorithm and c) a novel wavelet-packet identification procedure. Within class and between class dispersion metrics are adopted for the three data sets. A discrimination metric evaluates how well the three classes can be distinguished within the frequency range 0. 1 - 1.0 THz using the above algorithms.

Optimal discrimination and classification of THz spectra in the wavelet domain

In rapid scan Fourier transform spectrometry, we show that the noise in the wavelet coefficients resulting from the filter bank decomposition of the complex insertion loss function is linearly related to the noise power in the sample interferogram by a noise amplification factor. By maximizing an objective function composed of the power of the wavelet coefficients divided by the noise amplification factor, optimal feature extraction in the wavelet domain is performed. The performance of a classifier based on the output of a filter bank is shown to be considerably better than that of an Euclidean distance classifier in the original spectral domain. An optimization procedure results in a further improvement of the wavelet classifier. The procedure is suitable for enhancing the contrast or classifying spectra acquired by either continuous wave or THz transient spectrometers as well as for increasing the dynamic range of THz imaging systems. (C) 2003 Optical Society of America.

On the spectra of certain integro-differential-delay problems with applications in neurodynamics

We investigate the spectrum of certain integro-differential-delay equations (IDDEs) which arise naturally within spatially distributed, nonlocal, pattern formation problems. Our approach is based on the reformulation of the relevant dispersion relations with the use of the Lambert function. As a particular application of this approach, we consider the case of the Amari delay neural field equation which describes the local activity of a population of neurons taking into consideration the finite propagation speed of the electric signal. We show that if the kernel appearing in this equation is symmetric around some point a= 0 or consists of a sum of such terms, then the relevant dispersion relation yields spectra with an infinite number of branches, as opposed to finite sets of eigenvalues considered in previous works. Also, in earlier works the focus has been on the most rightward part of the spectrum and the possibility of an instability driven pattern formation. Here, we numerically survey the structure of the entire spectra and argue that a detailed knowledge of this structure is important within neurodynamical applications. Indeed, the Amari IDDE acts as a filter with the ability to recognise and respond whenever it is excited in such a way so as to resonate with one of its rightward modes, thereby amplifying such inputs and dampening others. Finally, we discuss how these results can be generalised to the case of systems of IDDEs.

Ferric and ferrous iron in nitroso-myoglobin: computer simulations of stable and metastable states and their infrared spectra

The binding of NO to iron is involved in the biological function of many heme proteins. Contrary to ligands like CO and O-2, which only bind to ferrous (Fe-II) iron, NO binds to both ferrous and ferric (Fe-II) iron. In a particular protein, the natural oxidation state can therefore be expected to be tailored to the required function. Herein, we present an ob initio potential-energy surface for ferric iron interacting with NO. This potential-energy surface exhibits three minima corresponding to eta'-NO coordination (the global minimum), eta(1)-ON coordination and eta(2) coordination. This contrasts with the potential-energy surface for Fe-II-NO, which ex- hibits only two minima (the eta(2) coordination mode for Fe-II is a transition state, not a minimum). In addition, the binding energies of NO are substantially larger for Fe-III than for Fe-II. We have performed molecular dynamics simulations for NO bound to ferric myoglobin (Mb(III)) and compare these with results obtained for Mb(II). Over the duration of our simulations (1.5 ns), all three binding modes are found to be stable at 200 K and transiently stable at 300 K, with eventual transformation to the eta(1)-NO global-minimum conformation. We discuss the implication of these results related to studies of rebinding processes in myoglobin.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

A new interpretation of the bonding properties and UV–vis spectra of [M3(CO)12] clusters (M = Ru, Os): a TD-DFT study

DFT and TD-DFT calculations (ADF program) were performed in order to analyze the electronic structure of the [M-3(CO)(12)] clusters (M = Ru, Os) and interpret their electronic spectra. The highest occupied molecular orbitals are M-M bonding (sigma) involving different M-M bonds, both for Ru and Os. They participate in low-energy excitation processes and their depopulation should weaken M-M bonds in general. While the LUMO is M-NI and M-CO anti-bonding (sigma*), the next, higher-lying empty orbitals have a main contribution from CO (pi*) and either a small (Ru) or an almost negligible one (Os) from the metal atoms. The main difference between the two clusters comes from the different nature of these low-energy unoccupied orbitals that have a larger metal contribution in the case of ruthenium. The photochemical reactivity of the two clusters is reexamined and compared to earlier interpretations.

Infrared absorption spectra, radiative efficiencies, and global warming potentials of perfluorocarbons: Comparison between experiment and theory

Experimentally and theoretically determined infrared spectra are reported for a series of straight-chain perfluorocarbons: C2F6, C3F8, C4F10, C5F12, C6F14, and C8F18. Theoretical spectra were determined using both density functional (DFT) and ab initio methods. Radiative efficiencies (REs) were determined using the method of Pinnock et al. (1995) and combined with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). Theoretically determined absorption cross sections were within 10% of experimentally determined values. Despite being much less computationally expensive, DFT calculations were generally found to perform better than ab initio methods. There is a strong wavenumber dependence of radiative forcing in the region of the fundamental C-F vibration, and small differences in wavelength between band positions determined by theory and experiment have a significant impact on the REs. We apply an empirical correction to the theoretical spectra and then test this correction on a number of branched chain and cyclic perfluoroalkanes. We then compute absorption cross sections, REs, and GWPs for an additional set of perfluoroalkenes.

Origin of electronic absorption spectra of MLCT-excited and one-electron reduced 2,2′-bipyridine and 1,10-phenanthroline complexes

UV–Vis absorption spectra of one-electron reduction products and 3MLCT excited states of [ReICl(CO)3- (N,N)] (N,N = 2,20-bipyridine, bpy; 1,10-phenanthroline, phen) have been measured by low-temperature spectroelectrochemistry and UV–Vis transient absorption spectroscopy, respectively, and assigned by open-shell TD-DFT calculations. The characters of the electronic transitions are visualized and analyzed using electron density redistribution maps. It follows that reduced and excited states can be approximately formulated as [ReICl(CO)3(N,Nÿ)]ÿ and ⁄[ReIICl(CO)3(N,Nÿ)], respectively. UV–Vis spectra of the reduced complexes are dominated by IL transitions, plus weaker MLCT contributions. Excited-state spectra show an intense band in the UV region of 50% IL origin mixed with LMCT (bpy, 373 nm) or MLCT (phen, 307 nm) excitations. Because of the significant IL contribution, this spectral feature is akin to the principal IL band of the anions. In contrast, the excited-state visible spectral pattern arises from predominantly LMCT transitions, any resemblance with the reduced-state visible spectra being coincidental. The Re complexes studied herein are representatives of a broad class of metal a-diimines, for which similar spectroscopic behavior can be expected.

Structure and spectra of dichloro(hydroxy-methoxy-di(2-pyridyl)methane)copper(II)

The non-electrolyte dichloro(hydroxy-methoxy-di(2-pyridylmethane)copper(II), resulting from the reaction of di(2-pyridyl)ketone and copper(II) chloride in methanol solution, was isolated and characterized and its structure was determined by X-ray diffraction. The pyridyl nitrogens and the chloride anions virtually from a basal plane in which lies the copper atom, while the oxygen of the methoxy group is in an apical position at a distance of 2.497 (3)Å. The nitrogenous base adopts the boat conformation with the pyridyl rings forming a dihedral angle of 108.72 (14)°. The nearest interatomic copper distance of 3.940(3)Å precludes copper-copper interactions, while the proximity of copper to the out-of-plane chlorine atoms [3.109(3)Å] suggests weakly bound chloro-bridged dimers. Spectral changes indicate that protic molecules displace the methoxy group and water affords the corresponding 1,1-diol.

Influence of hydrogen bonding in the structure of mixed-ligand copper(II) chelates. Spectra and structure of [aquo((2,2'-dipyridylamine)-(3-chloro-2,4-pentanedionato)copper(II)] nitrate

The IR, the ligand field spectra and the crystal structure of the mixed-ligand compound [(aquo)2,2P1 , a = 8.718(5), b = 9.407(5), c = 13.484 (7) Å, = 94.17(4)°, = 105.12(5)°, = 119.75(5)°, Z = 2, R = 0.0332, R W = 0.0869).

Synthesis, characterisation, Mössbauer spectra, and structures of some trinuclear iron-tin clusters

Reactions of [Fe3(CO)12] with diaryltin species SnR2(R1= 2,4,6-triisopropylphenyl, R2= 2,6-diethylphenyl, R3= pentamethylphenyl) and with Sn[CH(PPh2)2]2 have been investigated. The tin reagents SnR2(R = R1 or R2) reacted under mild conditions to give in moderate yields the trinuclear species [Fe2(CO)8(µ-SnR12)]1 or [Fe2(CO)8(µ-SnR22)]2, as orange-red crystalline solids, which decompose in air on prolonged exposure. The compound [Fe2(CO)8(µ-SnR42)]3(R4= 2,4,6-triphenylphenyl) can be similarly obtained. Prolonged treatment of the carbonyl with the novel tin reagent SnR32, by contrast, afforded the known compound spiro-[(OC)8Fe2SnFe2(CO)8]4 for which data are briefly reported. Reactions with tin or lead reagents M[CH(PPh2)2]2(M = Sn or Pb) afforded [Fe2(CO)6(µ-CO)(µ-dppm)][dppm = 1,2-bis(diphenylphosphino)methane] rapidly and almost quantitatively. Full crystal and molecular structural data are reported for [Fe2(CO)8(µ-SnR12)] and [Fe2(CO)8(µ-SnR22)]. Mössbauer data are also presented for compounds 1–3, and interpreted in terms of the structural data for these and other systems.

Raytracing simulation for predicting light pipe transmittance

Daylighting systems can offer energy savings primarily by reducing electric lighting usage. Accurate predictive models of daylighting system performances are crucial for effective design and implementation of this renewable energy technology. A comparative study of predictive methods was performed and the use of a commercial raytracing software program was validated as a method of predicting light pipe performance. Raytracing simulation was shown to more accurately predict transmission effi ciency than existing analytical methods.

The troposphere-to-stratosphere transition in kinetic energy spectra and nonlinear spectral fluxes as seen in ECMWF analyses

Global horizontal wavenumber kinetic energy spectra and spectral fluxes of rotational kinetic energy and enstrophy are computed for a range of vertical levels using a T799 ECMWF operational analysis. Above 250 hPa, the kinetic energy spectra exhibit a distinct break between steep and shallow spectral ranges, reminiscent of dual power-law spectra seen in aircraft data and high-resolution general circulation models. The break separates a large-scale ‘‘balanced’’ regime in which rotational flow strongly dominates divergent flow and a mesoscale ‘‘unbalanced’’ regime where divergent energy is comparable to or larger than rotational energy. Between 230 and 100 hPa, the spectral break shifts to larger scales (from n 5 60 to n 5 20, where n is spherical harmonic index) as the balanced component of the flow preferentially decays. The location of the break remains fairly stable throughout the stratosphere. The spectral break in the analysis occurs at somewhat larger scales than the break seen in aircraft data. Nonlinear spectral fluxes defined for the rotational component of the flow maximize between about 300 and 200 hPa. Large-scale turbulence thus centers on the extratropical tropopause region, within which there are two distinct mechanisms of upscale energy transfer: eddy–eddy interactions sourcing the transient energy peak in synoptic scales, and zonal mean–eddy interactions forcing the zonal flow. A well-defined downscale enstrophy flux is clearly evident at these altitudes. In the stratosphere, the transient energy peak moves to planetary scales and zonal mean–eddy interactions become dominant.