From error bounds to the complexity of first-order descent methods for convex functions

This paper shows that error bounds can be used as effective tools for deriving complexity results for first-order descent methods in convex minimization. In a first stage, this objective led us to revisit the interplay between error bounds and the Kurdyka-\L ojasiewicz (KL) inequality. One can show the equivalence between the two concepts for convex functions having a moderately flat profile near the set of minimizers (as those of functions with H\"olderian growth). A counterexample shows that the equivalence is no longer true for extremely flat functions. This fact reveals the relevance of an approach based on KL inequality. In a second stage, we show how KL inequalities can in turn be employed to compute new complexity bounds for a wealth of descent methods for convex problems. Our approach is completely original and makes use of a one-dimensional worst-case proximal sequence in the spirit of the famous majorant method of Kantorovich. Our result applies to a very simple abstract scheme that covers a wide class of descent methods. As a byproduct of our study, we also provide new results for the globalization of KL inequalities in the convex framework. Our main results inaugurate a simple methodology: derive an error bound, compute the desingularizing function whenever possible, identify essential constants in the descent method and finally compute the complexity using the one-dimensional worst case proximal sequence. Our method is illustrated through projection methods for feasibility problems, and through the famous iterative shrinkage thresholding algorithm (ISTA), for which we show that the complexity bound is of the form $O(q^{k})$ where the constituents of the bound only depend on error bound constants obtained for an arbitrary least squares objective with $\ell^1$ regularization

Bayesian and Markov chain Monte Carlo methods for identifying nonlinear systems in the presence of uncertainty

In this paper, the authors outline the general principles behind an approach to Bayesian system identification and highlight the benefits of adopting a Bayesian framework when attempting to identify models of nonlinear dynamical systems in the presence of uncertainty. It is then described how, through a summary of some key algorithms, many of the potential difficulties associated with a Bayesian approach can be overcome through the use of Markov chain Monte Carlo (MCMC) methods. The paper concludes with a case study, where an MCMC algorithm is used to facilitate the Bayesian system identification of a nonlinear dynamical system from experimentally observed acceleration time histories

Benzocaine-Induced Methemoglobinemia: A Condition of which all Endoscopists Should Be Aware

Methemoglobinemia is a rare complication that can occur with the use of benzocaine-containing compounds. Two cases of methemoglobinemia are reported, and the pathophysiology and treatment of methemoglobinemia are reviewed. Both patients received topical 20% benzocaine spray before endoscopy. Immediately following the procedure, there was a reduction in O2 saturation assessed by pulse oximetry that was refractory to O2 therapy. Dramatic peripheral and central cyanosis developed. O2 saturation measured by pulse oximetry ranged from 83% to 87% on O2 by nasal prongs and 100% O2 by a nonrebreathing mask. Both patients were mildly confused and one patient complained of a significant headache. The diagnosis of methemoglobinemia was considered and arterial blood gas sampling was performed. In both patients, the arterial blood had a chocolate brown colour. A methemoglobin level of 48% and 18% was noted in patient 1 and patient 2, respectively. Both patients were treated with methylene blue, resulting in a significant improvement with gradual normalization of their O2 saturation within 10 min to 30 min. The use of benzocaine spray may not markedly alter the patient's perception of endoscopy and thus, the routine use of these agents should be questioned. If such agents are used, the physician must be aware of this association to prevent a delay in the diagnosis and management of this rare, but potentially lethal, condition

Invasion speeds for structured populations in fluctuating environments

We live in a time where climate models predict future increases in environmental variability and biological invasions are becoming increasingly frequent. A key to developing effective responses to biological invasions in increasingly variable environments will be estimates of their rates of spatial spread and the associated uncertainty of these estimates. Using stochastic, stage-structured, integro-difference equation models, we show analytically that invasion speeds are asymptotically normally distributed with a variance that decreases in time. We apply our methods to a simple juvenile-adult model with stochastic variation in reproduction and an illustrative example with published data for the perennial herb, \emph{Calathea ovandensis}. These examples buttressed by additional analysis reveal that increased variability in vital rates simultaneously slow down invasions yet generate greater uncertainty about rates of spatial spread. Moreover, while temporal autocorrelations in vital rates inflate variability in invasion speeds, the effect of these autocorrelations on the average invasion speed can be positive or negative depending on life history traits and how well vital rates ``remember'' the past

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.Comment: 4 pages, 1 figur

Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at s√=8 TeV with the ATLAS detector

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of 20.1fb−1 of proton–proton collision data at √s = 8 TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via t˜→tχ˜01 or t˜→ bχ˜±1 →bW(∗)χ˜01 , where χ˜01 (χ˜±1 ) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of t˜ → tχ˜01 . For a branching fraction of 100%, top squark masses in the range 270–645 GeV are excluded for χ˜01 masses below 30 GeV. For a branching fraction of 50% to either t˜ → tχ˜01 or t˜ → bχ˜±1 , and assuming the χ˜±1 mass to be twice the χ˜01 mass, top squark masses in the range 250–550 GeV are excluded for χ˜01 masses below 60 GeV

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT

Creatine ingestion augments dietary carbohydrate mediated muscle glycogen supercomposition during the initial 24 hrs of recovery following prolonged exhaustive exercise in humans

Muscle glycogen availability can limit endurance exercise performance. We previously demonstrated 5 days of creatine (Cr) and carbohydrate (CHO) ingestion augmented post-exercise muscle glycogen storage compared to CHO feeding alone in healthy volunteers. Here we aimed to characterise the time-course of this Cr-induced response under more stringent and controlled experimental conditions and identify potential mechanisms underpinning this phenomenon. Fourteen healthy, male volunteers cycled to exhaustion at 70% VO2peak. Muscle biopsies were obtained at rest immediately post-exercise and after 1, 3 and 6 days of recovery, during which Cr or placebo supplements (20g.day-1) were ingested along with a prescribed high CHO diet (37.5 kcal.kg body mass-1.day-1, >80% calories CHO). Oral-glucose tolerance tests (oral-GTT) were performed pre-exercise and after 1, 3 and 6 days of Cr and placebo supplementation. Exercise depleted muscle glycogen content to the same extent in both treatment groups. Creatine supplementation increased muscle total-Cr, free-Cr and phosphocreatine (PCr) content above placebo following 1, 3 and 6 days of supplementation (all P<0.05). Creatine supplementation also increased muscle glycogen content noticeably above placebo after 1 day of supplementation (P<0.05), which was sustained thereafter. This study confirmed dietary Cr augments post-exercise muscle glycogen super-compensation, and demonstrates this occurred during the initial 24 h of post-exercise recovery (when muscle total-Cr had increased by <10%). This marked response ensued without apparent treatment differences in muscle insulin sensitivity (oral-GTT, muscle GLUT4 mRNA), osmotic stress (muscle c-fos and HSP72 mRNA) or muscle cell volume (muscle water content) responses, such that another mechanism must be causative