A four-timescale algorithm for constrained stochastic optimization of RED

The overall performance of random early detection (RED) routers in the Internet is determined by the settings of their associated parameters. The non-availability of a functional relationship between the RED performance and its parameters makes it difficult to implement optimization techniques directly in order to optimize the RED parameters. In this paper, we formulate a generic optimization framework using a stochastically bounded delay metric to dynamically adapt the RED parameters. The constrained optimization problem thus formulated is solved using traditional nonlinear programming techniques. Here, we implement the barrier and penalty function approaches, respectively. We adopt a second-order nonlinear optimization framework and propose a novel four-timescale stochastic approximation algorithm to estimate the gradient and Hessian of the barrier and penalty objectives and update the RED parameters. A convergence analysis of the proposed algorithm is briefly sketched. We perform simulations to evaluate the performance of our algorithm with both barrier and penalty objectives and compare these with RED and a variant of it in the literature. We observe an improvement in performance using our proposed algorithm over RED, and the above variant of it.

SPSA algorithms with measurement reuse

Four algorithms, all variants of Simultaneous Perturbation Stochastic Approximation (SPSA), are proposed. The original one-measurement SPSA uses an estimate of the gradient of objective function L containing an additional bias term not seen in two-measurement SPSA. As a result, the asymptotic covariance matrix of the iterate convergence process has a bias term. We propose a one-measurement algorithm that eliminates this bias, and has asymptotic convergence properties making for easier comparison with the two-measurement SPSA. The algorithm, under certain conditions, outperforms both forms of SPSA with the only overhead being the storage of a single measurement. We also propose a similar algorithm that uses perturbations obtained from normalized Hadamard matrices. The convergence w.p. 1 of both algorithms is established. We extend measurement reuse to design two second-order SPSA algorithms and sketch the convergence analysis. Finally, we present simulation results on an illustrative minimization problem.

Dynamics and local structure of colossal magnetoresistance manganites

We report Extended X-ray Absorption Fine Structure and anelastic spectroscopy measurements on on hole doped manganese oxides La1-xCaxMnO3 which present the colossal magnetoresistance effect. EXAFS measurements were realized both in the absence and presence of an applied magnetic field of 1.1 Tesla, in a wide temperature range (between 330 and 77 K) and at various dopings (x = 0.25 and x = 0.33). The magnetic field orders the magnetic moments so favouring the electron mobility and the reduction of Mn-O octahedra distortions. We observe the presence of four short and two long Mn-O distances (1.93 and 2.05 Angstrom respectively) above and also below the metal-insulator phase transition. The overall distortion decreases but does not completely disappear in the metallic phase suggesting the possible coexistence of metallic and insulating regions at low temperatures. The magnetic field reduces the lattice distortions showing evidence of a microscopic counterpart of the macroscopic colossal magnetoresistance. We also present preliminary anelastic relaxation spectra in a wide temperature range from 900 K to 1 K on a sample with x = 0.40, in order to study the structural phase transitions and the lattice dynamics. A double peak has been observed at the metal-insulator transition in the imaginary part of Young's modulus. This double peak indicates that the metal-insulator transition could be a more complex phenomenon than a simple second order phase transition. In particular the peak at lower temperatures can be connected with the possible presence of inhomogeneous phase structures. Another intense dissipation peak has been observed corresponding to the structural orthorhombic-trigonal transition around 750 K.

Phase diagram of a hard-sphere system in a quenched random potential: A numerical study

We report numerical results for the phase diagram in the density-disorder plane of a hard-sphere system in the presence of quenched, random, pinning disorder. Local minima of a discretized version of the Ramakrishnan-Yussouff free energy functional are located numerically and their relative stability is studied as a function of the density and the strength of disorder. Regions in the phase diagram corresponding to liquid, glassy, and nearly crystalline states are mapped out, and the nature of the transitions is determined. The liquid to glass transition changes from first to second order as the strength of the disorder is increased. For weak disorder, the system undergoes a first-order crystallization transition as the density is increased. Beyond a critical value of the disorder strength, this transition is replaced by a continuous glass transition. Our numerical results are compared with those of analytical work on the same system. Implications of our results for the field-temperature phase diagram of type-II superconductors are discussed.

Extension of Darken's Quadratic Formalism to Dilute Multicomponent Solutions

Darken's quadratic formalism is extended to multicomponent solutions. Equations are developed for the representation of the integral and partial excess free energies, entropies and enthalpies in dilute multicomponent solutions. Quadratic formalism applied to multicomponent solutions is thermodynamically consistent. The formalism is compared with the conventional second order Maclaurin series or interaction parameter representation and the relations between them are derived. Advantages of the quadratic formalism are discussed.

A kinetic study of liquid-phase catalytic oxidation of styrene to benzaldehyde with Wilkinson complex

Liquid-phase homogeneous catalytic oxidation of styrene with Wilkinson complex by molecular oxygen in toluene medium gave selectively benzaldehyde and formaldehyde as the primary products. Higher temperatures and styrene conversions eventually led to acid formation due to co-oxidation of aldehyde.A reaction induction period and an initiation period, typical of free-radical reactions, characterized the oxidation process. The effects of temperature and catalyst and styrene concentrations on the conversion of styrene to benzaldehyde and acid formation have been studied. The optimum reaction parameters have been determined as a styrene-to-solvent mole ratio of 0.5, a catalyst-to-styrene mole ratio of 5.0 X lo4, and a reaction temperature of 75 "C. A reaction scheme based upon free-radical mechanism yielded a pseudo-first-order model which agreed well with the observed kinetic data in the absence of co-oxidation of aldehyde. A second-order model was found to fit the experimental data better in the case of aldehyde conversion to acid.

FTIR-spektroskopia kasviuutteiden analysoinnissa

FTIR-spektroskopia (Fourier-muunnosinfrapunaspektroskopia) on nopea analyysimenetelmä. Fourier-laitteissa interferometrin käyttäminen mahdollistaa koko infrapunataajuusalueen mittaamisen muutamassa sekunnissa. ATR-liitännäisellä varustetun FTIR-spektrometrin käyttö ei edellytä juuri näytteen valmistusta ja siksi menetelmä on käytössä myös helppo. ATR-liitännäinen mahdollistaa myös monien erilaisten näytteiden analysoinnin. Infrapunaspektrin mittaaminen onnistuu myös sellaisista näytteistä, joille perinteisiä näytteenvalmistusmenetelmiä ei voida käyttää. FTIR-spektroskopian avulla saatu tieto yhdistetään usein tilastollisiin monimuuttuja-analyyseihin. Klusterianalyysin avulla voidaan spektreistä saatu tieto ryhmitellä samanlaisuuteen perustuen. Hierarkkisessa klusterianalyysissa objektien välinen samanlaisuus määritetään laskemalla niiden välinen etäisyys. Pääkomponenttianalyysin avulla vähennetään datan ulotteisuutta ja luodaan uusia korreloimattomia pääkomponentteja. Pääkomponenttien tulee säilyttää mahdollisimman suuri määrä alkuperäisen datan variaatiosta. FTIR-spektroskopian ja monimuuttujamenetelmien sovellusmahdollisuuksia on tutkittu paljon. Elintarviketeollisuudessa sen soveltuvuutta esimerkiksi laadun valvontaan on tutkittu. Menetelmää on käytetty myös haihtuvien öljyjen kemiallisten koostumusten tunnistukseen sekä öljykasvien kemotyyppien havaitsemiseen. Tässä tutkimuksessa arvioitiin menetelmän käyttöä suoputken uutenäytteiden luokittelussa. Tutkimuksessa suoputken eri kasvinosien uutenäytteiden FTIR-spektrejä vertailtiin valikoiduista puhdasaineista mitattuihin FTIR-spektreihin. Puhdasaineiden FTIR-spektreistä tunnistettiin niiden tyypilliset absorptiovyöhykkeet. Furanokumariinien spektrien intensiivisten vyöhykkeiden aaltolukualueet valittiin monimuuttuja-analyyseihin. Monimuuttuja-analyysit tehtiin myös IR-spektrin sormenjälkialueelta aaltolukualueelta 1785-725 cm-1. Uutenäytteitä pyrittiin luokittelemaan niiden keräyspaikan ja kumariinipitoisuuden mukaan. Keräyspaikan mukaan ryhmittymistä oli havaittavissa, mikä selittyi vyöhykkeiden aaltolukualueiden mukaan tehdyissä analyyseissa pääosin kumariinipitoisuuksilla. Näissä analyyseissa uutenäytteet pääosin ryhmittyivät ja erottuivat kokonaiskumariinipitoisuuksien mukaan. Myös aaltolukualueen 1785-725 cm-1 analyyseissa havaittiin keräyspaikan mukaan ryhmittymistä, mitä kumariinipitoisuudet eivät kuitenkaan selittäneet. Näihin ryhmittymisiin vaikuttivat mahdollisesti muiden yhdisteiden samanlaiset pitoisuudet näytteissä. Analyyseissa käytettiin myös muita aaltolukualueita, mutta tulokset eivät juuri poikenneet aiemmista. 2. kertaluvun derivaattaspektrien monimuuttuja-analyysit sormenjälkialueelta eivät myöskään muuttaneet tuloksia havaittavasti. Jatkotutkimuksissa nyt käytettyä menetelmää on mahdollista edelleen kehittää esimerkiksi tutkimalla monimuuttuja-analyyseissa 2. kertaluvun derivaattaspektreistä suppeampia, tarkkaan valittuja aaltolukualueita.

Studies in nonlinear optical materials: structure of di-2-menthyl 2-(N,N'-dimethyl-2-imidazolidinylidene)malonate monohydrate

C28H48N2Oa.H2 O, Mr=494.7, orthorhombic,P2~2~2~, a = 7.634 (2), b = 11.370 (2), c=34. 167 (4) A, V = 2966 (2) A 3, Z = 4, D m = 1.095,D x -- 1. 108 g cm -3, Mo Kct, 2 -- 0.7107 ,/k, ~ =0.43 cm -~, F(000) = 1088.0, T= 293 K, R = 0.061 for 1578 significant reflections. The second-harmonicgeneration (SHG) efficiency of this compound is negligible (1/100th of the urea standard). The observed low second-order nonlinear response has been attributed to the unfavourable packing of the molecules in the crystal lattice.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Probability Distribution of the Eigenvalues of the Random String Equation

The probability distribution of the eigenvalues of a second-order stochastic boundary value problem is considered. The solution is characterized in terms of the zeros of an associated initial value problem. It is further shown that the probability distribution is related to the solution of a first-order nonlinear stochastic differential equation. Solutions of this equation based on the theory of Markov processes and also on the closure approximation are presented. A string with stochastic mass distribution is considered as an example for numerical work. The theoretical probability distribution functions are compared with digital simulation results. The comparison is found to be reasonably good.

Estimation of variables at ungaged locations by empirical orthogonal functions

A technique based on empirical orthogonal functions is used to estimate hydrologic time-series variables at ungaged locations. The technique is applied to estimate daily and monthly rainfall, temperature and runoff values. The accuracy of the method is tested by application to locations where data are available. The second-order characteristics of the estimated data are compared with those of the observed data. The results indicate that the method is quick and accurate.

Infinite dimensional slow modulations in a well known delayed model for cutting tool vibrations

We apply the method of multiple scales (MMS) to a well-known model of regenerative cutting vibrations in the large delay regime. By ``large'' we mean the delay is much larger than the timescale of typical cutting tool oscillations. The MMS up to second order, recently developed for such systems, is applied here to study tool dynamics in the large delay regime. The second order analysis is found to be much more accurate than the first order analysis. Numerical integration of the MMS slow flow is much faster than for the original equation, yet shows excellent accuracy in that plotted solutions of moderate amplitudes are visually near-indistinguishable. The advantages of the present analysis are that infinite dimensional dynamics is retained in the slow flow, while the more usual center manifold reduction gives a planar phase space; lower-dimensional dynamical features, such as Hopf bifurcations and families of periodic solutions, are also captured by the MMS; the strong sensitivity of the slow modulation dynamics to small changes in parameter values, peculiar to such systems with large delays, is seen clearly; and though certain parameters are treated as small (or, reciprocally, large), the analysis is not restricted to infinitesimal distances from the Hopf bifurcation.

Constructions and situations : a constructivist reading of Hegel's System

The object of this work is Hegel's Logic, which comprises the first third of his philosophical System that also includes the Philosophy of Nature and the Philosophy of Spirit. The work is divided into two parts, where the first part investigates Hegel s Logic in itself or without an explicit reference to rest of Hegel's System. It is argued in the first part that Hegel's Logic contains a methodology for constructing examples of basic ontological categories. The starting point on which this construction is based is a structure Hegel calls Nothing, which I argue to be identical with an empty situation, that is, a situation with no objects in it. Examples of further categories are constructed, firstly, by making previous structures objects of new situations. This rule makes it possible for Hegel to introduce examples of ontological structures that contain objects as constituents. Secondly, Hegel takes also the very constructions he uses as constituents of further structures: thus, he is able to exemplify ontological categories involving causal relations. The final result of Hegel's Logic should then be a model of Hegel s Logic itself, or at least of its basic methods. The second part of the work focuses on the relation of Hegel's Logic to the other parts of Hegel's System. My interpretation tries to avoid, firstly, the extreme of taking Hegel's System as a grand metaphysical attempt to deduce what exists through abstract thinking, and secondly, the extreme of seeing Hegel's System as mere diluted Kantianism or a second-order investigation of theories concerning objects instead of actual objects. I suggest a third manner of reading Hegel's System, based on extending the constructivism of Hegel's Logic to the whole of his philosophical System. According to this interpretation, transitions between parts of Hegel's System should not be understood as proofs of any sort, but as constructions of one structure or its model from another structure. Hence, these transitions involve at least, and especially within the Philosophy of Nature, modelling of one type of object or phenomenon through characteristics of an object or phenomenon of another type, and in the best case, and especially within the Philosophy of Spirit, transformations of an object or phenomenon of one type into an object or phenomenon of another type. Thus, the transitions and descriptions within Hegel's System concern actual objects and not mere theories, but they still involve no fallacious deductions.

Exact nonlinear travelling hydromagnetic wave solutions

Exact travelling wave solutions for hydromagnetic waves in an exponentially stratified incompressible medium are obtained. With the help of two integrals it becomes possible to reduce the system of seven nonlinear PDE's to a second order nonlinear ODE which describes an one dimensional harmonic oscillator with a nonlinear friction term. This equation is studied in detail in the phase plane. The travelling waves are periodic only when they propagate either horizontally or vertically. The reduced second order nonlinear differential equation describing the travelling waves in inhomogeneous conducting media has rather ubiquitous nature in that it also appears in other geophysical systems such as internal waves, Rossby waves and topographic Rossby waves in the ocean.

Attitude estimation using moment invariants

The problem of estimating the three-dimensional rotational parameters of a rigid body from its monocular image data has been considered using the method of moment invariants. Second- and third-order moment invariants are used to construct the feature vector for the scale and orientation independent identification of the camera view axis direction in the body-fixed reference frame. The camera rotation angle about the view axis is derived from second-order central moments. The relative attitude of the rigid body is then expressed in terms of quaternion parameters to model the outputs of a video sensor in attitude control simulations. Experimental results and simulation outputs are presented using the mathematical model of a spacecraft.