Normal forms, stability and splitting of invariant manifolds I. Gevrey Hamiltonians

In this paper, we give a new construction of resonant normal forms with a small remainder for near-integrable Hamiltonians at a quasi-periodic frequency. The construction is based on the special case of a periodic frequency, a Diophantine result concerning the approximation of a vector by independent periodic vectors and a technique of composition of periodic averaging. It enables us to deal with non-analytic Hamiltonians, and in this first part we will focus on Gevrey Hamiltonians and derive normal forms with an exponentially small remainder. This extends a result which was known for analytic Hamiltonians, and only in the periodic case for Gevrey Hamiltonians. As applications, we obtain an exponentially large upper bound on the stability time for the evolution of the action variables and an exponentially small upper bound on the splitting of invariant manifolds for hyperbolic tori, generalizing corresponding results for analytic Hamiltonians.

Normal forms, stability and splitting of invariant manifolds II. Finitely differentiable Hamiltonians

This paper is a sequel to ``Normal forms, stability and splitting of invariant manifolds I. Gevrey Hamiltonians", in which we gave a new construction of resonant normal forms with an exponentially small remainder for near-integrable Gevrey Hamiltonians at a quasi-periodic frequency, using a method of periodic approximations. In this second part we focus on finitely differentiable Hamiltonians, and we derive normal forms with a polynomially small remainder. As applications, we obtain a polynomially large upper bound on the stability time for the evolution of the action variables and a polynomially small upper bound on the splitting of invariant manifolds for hyperbolic tori.

Finite field treatment of vibrational polarizabilities and hyperpolarizabilities: on the role of the Eckart conditions, their implementation, and their use in characterizing key vibrations

In the finite field (FF) treatment of vibrational polarizabilities and hyperpolarizabilities, the field-free Eckart conditions must be enforced in order to prevent molecular reorientation during geometry optimization. These conditions are implemented for the first time. Our procedure facilities identification of field-induced internal coordinates that make the major contribution to the vibrational properties. Using only two of these coordinates, quantitative accuracy for nuclear relaxation polarizabilities and hyperpolarizabilities is achieved in π-conjugated systems. From these two coordinates a single most efficient natural conjugation coordinate (NCC) can be extracted. The limitations of this one coordinate approach are discussed. It is shown that the Eckart conditions can lead to an isotope effect that is comparable to the isotope effect on zero-point vibrational averaging, but with a different mass-dependence

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Methodology for studying the numerical speed of sound in finite differences schemes

The space and time discretization inherent to all FDTD schemesintroduce non-physical dispersion errors, i.e. deviations ofthe speed of sound from the theoretical value predicted bythe governing Euler differential equations. A generalmethodologyfor computing this dispersion error via straightforwardnumerical simulations of the FDTD schemes is presented.The method is shown to provide remarkable accuraciesof the order of 1/1000 in a wide variety of twodimensionalfinite difference schemes.

Impedance boundary conditions for pseudo-spectral time-domain methods in room acoustics

The pseudo-spectral time-domain (PSTD) method is an alternative time-marching method to classicalleapfrog finite difference schemes in the simulation of wave-like propagating phenomena. It is basedon the fundamentals of the Fourier transform to compute the spatial derivatives of hyperbolic differential equations. Therefore, it results in an isotropic operator that can be implemented in an efficient way for room acoustics simulations. However, one of the first issues to be solved consists on modeling wallabsorption. Unfortunately, there are no references in the technical literature concerning to that problem. In this paper, assuming real and constant locally reacting impedances, several proposals to overcome this problem are presented, validated and compared to analytical solutions in different scenarios.

Impedance boundary conditions for pseudo-spectral time-domain methods in room acoustics

The Pseudo-Spectral Time Domain (PSTD) method is an alternative time-marching method to classical leapfrog finite difference schemes inthe simulation of wave-like propagating phenomena. It is based on the fundamentals of the Fourier transform to compute the spatial derivativesof hyperbolic differential equations. Therefore, it results in an isotropic operator that can be implemented in an efficient way for room acousticssimulations. However, one of the first issues to be solved consists on modeling wall absorption. Unfortunately, there are no references in thetechnical literature concerning to that problem. In this paper, assuming real and constant locally reacting impedances, several proposals toovercome this problem are presented, validated and compared to analytical solutions in different scenarios.

The control of microtubule stability in vitro and in transfected cells by MAP1B and SCG10.

In neurons, the regulation of microtubules plays an important role for neurite outgrowth, axonal elongation, and growth cone steering. SCG10 family proteins are the only known neuronal proteins that have a strong destabilizing effect, are highly enriched in growth cones and are thought to play an important role during axonal elongation. MAP1B, a microtubule-stabilizing protein, is found in growth cones as well, therefore it was important to test their effect on microtubules in the presence of both proteins. We used recombinant proteins in microtubule assembly assays and in transfected COS-7 cells to analyze their combined effects in vitro and in living cells, respectively. Individually, both proteins showed their expected activities in microtubule stabilization and destruction respectively. In MAP1B/SCG10 double-transfected cells, MAP1B could not protect microtubules from SCG10-induced disassembly in most cells, in particular not in cells that contained high levels of SCG10. This suggests that SCG10 is more potent to destabilize microtubules than MAP1B to rescue them. In microtubule assembly assays, MAP1B promoted microtubule formation at a ratio of 1 MAP1B per 70 tubulin dimers while a ratio of 1 SCG10 per two tubulin dimers was needed to destroy microtubules. In addition to its known binding to tubulin dimers, SCG10 binds also to purified microtubules in growth cones of dorsal root ganglion neurons in culture. In conclusion, neuronal microtubules are regulated by antagonistic effects of MAP1B and SCG10 and a fine tuning of the balance of these proteins may be critical for the regulation of microtubule dynamics in growth cones.

Review of Stability Berm Alternatives for Environmentally Senstitive Areas, 2005

Stability berms are commonly constructed where roadway embankments cross soft or unstable ground conditions. Under certain circumstances, the construction of stability berms cause unfavorable environmental impacts, either directly or indirectly, through their effect on wetlands, endangered species habitat, stream channelization, longer culvert lengths, larger right-of-way purchases, and construction access limits. Due to an ever more restrictive regulatory environment, these impacts are problematic. The result is the loss of valuable natural resources to the public, lengthy permitting review processes for the department of transportation and permitting agencies, and the additional expenditures of time and money for all parties. The purpose of this project was to review existing stability berm alternatives for potential use in environmentally sensitive areas. The project also evaluates how stabilization technologies are made feasible, desirable, and cost-effective for transportation projects and determines which alternatives afford practical solutions for avoiding and minimizing impacts to environmentally sensitive areas. An online survey of engineers at state departments of transportation was also conducted to assess the frequency and cost effectiveness of the various stabilization technologies. Geotechnical engineers that responded to the survey overwhelmingly use geosynthetic reinforcement as a suitable and cost-effective solution for stabilizing embankments and cut slopes. Alternatively, chemical stabilization and installation of lime/cement columns is rarely a remediation measure employed by state departments of transportation.

Sequential estimation of time-varying multipath channel for MIMO-OFDM systems

In this paper, we introduce a pilot-aided multipath channel estimator for Multiple-Input Multiple-Output (MIMO) Orthogonal Frequency Division Multiplexing (OFDM) systems. Typical estimation algorithms assume the number of multipath components and delays to be known and constant, while theiramplitudes may vary in time. In this work, we focus on the more realistic assumption that also the number of channel taps is unknown and time-varying. The estimation problem arising from this assumption is solved using Random Set Theory (RST), which is a probability theory of finite sets. Due to the lack of a closed form of the optimal filter, a Rao-Blackwellized Particle Filter (RBPF) implementation of the channel estimator is derived. Simulation results demonstrate the estimator effectiveness.

Determination of Vasopressin and Desmopressin in urine by means of liquid chromatography coupled to quadrupole time-of-flight mass spectrometry for doping control purposes.

The anti-diuretic neurohypophysial hormone Vasopressin (Vp) and its synthetic analogue Desmopressin (Dp, 1-desamino-vasopressin) have received considerable attention from doping control authorities due to their impact on physiological blood parameters. Accordingly, the illicit use of Desmopressin in elite sport is sanctioned by the World Anti-Doping Agency (WADA) and the drug is classified as masking agent. Vp and Dp are small (8-9 amino acids) peptides administered orally as well as intranasally. Within the present study a method to determine Dp and Vp in urinary doping control samples by means of liquid chromatography coupled to quadrupole high resolution time-of-flight mass spectrometry was developed. After addition of Lys-Vasopressin as internal standard and efficient sample clean up with a mixed mode solid phase extraction (weak cation exchange), the samples were directly injected into the LC-MS system. The method was validated considering the parameters specificity, linearity, recovery (80-100%), accuracy, robustness, limit of detection/quantification (20/50 pg mL(-1)), precision (inter/intra-day<10%), ion suppression and stability. The analysis of administration study urine samples collected after a single intranasal or oral application of Dp yielded in detection windows for the unchanged target analyte for up to 20 h at concentrations between 50 and 600 pg mL(-1). Endogenous Vp was detected in concentrations of approximately 20-200 pg mL(-1) in spontaneous urine samples obtained from healthy volunteers. The general requirements of the developed method provide the characteristics for an easy transfer to other anti-doping laboratories and support closing another potential gap for cheating athletes.

Non-Quasi-Competitiveness and Stability in Cournot's Model

This paper retakes previous work of the authors, about the relationship between non-quasi-competitiveness (the increase in price caused by an increase in the number of oligopolists) and stability of the equilibrium in the classical Cournot oligopoly model. Though it has been widely accepted in the literature that the loss of quasi-competitiveness is linked, in the long run as new firms entered the market, to instability of the model, the authors in their previous work put forward a model in which a situation of monopoly changed to duopoly losing quasi-competitiveness but maintaining the stability of the equilibrium. That model could not, at the time, be extended to any number of oligopolists. The present paper exhibits such an extension. An oligopoly model is shown in which the loss of quasi-competitiveness resists the presence in the market of as many firms as one wishes and where the successive Cournot's equilibrium points are unique and asymptotically stable. In this way, for the first time, the conjecture that non-quasi- competitiveness and instability were equivalent in the long run, is proved false.

Do orthopaedic shoes improve local dynamic stability of gait? An observational study in patients with chronic foot and ankle injuries.

BACKGROUND: Complex foot and ankle fractures, such as calcaneum fractures or Lisfranc dislocations, are often associated with a poor outcome, especially in terms of gait capacity. Indeed, degenerative changes often lead to chronic pain and chronic functional limitations. Prescription footwear represents an important therapeutic tool during the rehabilitation process. Local Dynamic Stability (LDS) is the ability of locomotor system to maintain continuous walking by accommodating small perturbations that occur naturally during walking. Because it reflects the degree of control over the gait, LDS has been advocated as a relevant indicator for evaluating different conditions and pathologies. The aim of this study was to analyze changes in LDS induced by orthopaedic shoes in patients with persistent foot and ankle injuries. We hypothesised that footwear adaptation might help patients to improve gait control, which could lead to higher LDS: METHODS: Twenty-five middle-aged inpatients (5 females, 20 males) participated in the study. They were treated for chronic post-traumatic disabilities following ankle and/or foot fractures in a Swiss rehabilitation clinic. During their stay, included inpatients received orthopaedic shoes with custom-made orthoses (insoles). They performed two 30s walking trials with standard shoes and two 30s trials with orthopaedic shoes. A triaxial motion sensor recorded 3D accelerations at the lower back level. LDS was assessed by computing divergence exponents in the acceleration signals (maximal Lyapunov exponents). Pain was evaluated with Visual Analogue Scale (VAS). LDS and pain differences between the trials with standard shoes and the trials with orthopaedic shoes were assessed. RESULTS: Orthopaedic shoes significantly improved LDS in the three axes (medio-lateral: 10% relative change, paired t-test p < 0.001; vertical: 9%, p = 0.03; antero-posterior: 7%, p = 0.04). A significant decrease in pain level (VAS score -29%) was observed. CONCLUSIONS: Footwear adaptation led to pain relief and to improved foot & ankle proprioception. It is likely that that enhancement allows patients to better control foot placement. As a result, higher dynamic stability has been observed. LDS seems therefore a valuable index that could be used in early evaluation of footwear outcome in clinical settings.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search.

BACKGROUND: The reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used, highly sensitive laboratory technique to rapidly and easily detect, identify and quantify gene expression. Reliable RT-qPCR data necessitates accurate normalization with validated control genes (reference genes) whose expression is constant in all studied conditions. This stability has to be demonstrated.We performed a literature search for studies using quantitative or semi-quantitative PCR in the rat spared nerve injury (SNI) model of neuropathic pain to verify whether any reference genes had previously been validated. We then analyzed the stability over time of 7 commonly used reference genes in the nervous system - specifically in the spinal cord dorsal horn and the dorsal root ganglion (DRG). These were: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) and L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and hydroxymethylbilane synthase (HMBS). We compared the candidate genes and established a stability ranking using the geNorm algorithm. Finally, we assessed the number of reference genes necessary for accurate normalization in this neuropathic pain model. RESULTS: We found GAPDH, HMBS, Actb, HPRT1 and 18S cited as reference genes in literature on studies using the SNI model. Only HPRT1 and 18S had been once previously demonstrated as stable in RT-qPCR arrays. All the genes tested in this study, using the geNorm algorithm, presented gene stability values (M-value) acceptable enough for them to qualify as potential reference genes in both DRG and spinal cord. Using the coefficient of variation, 18S failed the 50% cut-off with a value of 61% in the DRG. The two most stable genes in the dorsal horn were RPL29 and RPL13a; in the DRG they were HPRT1 and Actb. Using a 0.15 cut-off for pairwise variations we found that any pair of stable reference gene was sufficient for the normalization process. CONCLUSIONS: In the rat SNI model, we validated and ranked Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 and 18S as good reference genes in the spinal cord. In the DRG, 18S did not fulfill stability criteria. The combination of any two stable reference genes was sufficient to provide an accurate normalization.

Hot vs. cold: Sequential responses and preference stability in experimental games

In experiments with two-person sequential games we analyzewhether responses to favorable and unfavorable actions dependon the elicitation procedure. In our hot treatment thesecond player responds to the first player s observed actionwhile in our cold treatment we follow the strategy method and have the second player decide on a contingent action foreach and every possible first player move, without firstobserving this move. Our analysis centers on the degree towhich subjects deviate from the maximization of their pecuniaryrewards, as a response to others actions. Our results show nodifference in behavior between the two treatments. We also findevidence of the stability of subjects preferences with respectto their behavior over time and to the consistency of theirchoices as first and second mover.