The Dynamics of Quadropoly: League Position in English Football between 1888 and 2010

The paper explores competitive balance in top tier English league football from its inception in 1888. It examines the extent to which finishing in the top four positions in successive seasons is the preserve of a small number of clubs. Using a range of statistical measures, the analysis shows that the current high levels of competitive imbalance are not new phenomena. The overall pattern approximates a ‘U curve’: current patterns parallel those in the 1890s. In the early years of English league football, differences in resources between clubs soon became apparent. Clubs from the larger conurbations generated consistently larger revenues than their counterparts in the smaller industrial towns. This was primarily the result of the larger crowds that they could attract to their home games. This enabled them to entice the best players to their clubs away from their smaller rivals. The introduction of the maximum wage in 1901 and the transfer system helped to stem these increasing inequalities between clubs. This coincided with a massive wave of new stadia construction which enabled all the clubs to compete on an increasingly level playing field. These conjunctural changes to English football before 1915 produced the era of relatively competitive football during the inter-war years. This continued until the abolition of the maximum wage in 1961. Since that time, competitive balance has reversed and become increasingly restricted. English top-tier football has re-entered an era of extreme competitive imbalance.

Reactions of atomic and molecular ions with acetone, 1,1,1- trifluoroacetone, and hexafluoroacetone: An investigation of the effects of molecular structure on the dynamics and kinetics of ion-molecule reactions

A selected ion flow tube study of the reactions of a series of gas-phase atomic cations (S⁺, Xe⁺, O⁺, Kr⁺, N⁺, Ar⁺ and Ne⁺) and molecular ions (SF _n⁺ (n = 1-5), CF_n⁺ (n = 1-3), CF₂Cl⁺, H₃O⁺, NO⁺, N ₂O⁺, CO₂⁺, CO⁺, and N₂⁺) spanning a large range of recombination energies (6.3-21.6 eV), with acetone, 1,1,1-trifluoroacetone, and hexafluoroacetone has been undertaken with the objective of exploring the nature of the reaction ion chemistry as the methyl groups in acetone are substituted for CF₃. The reaction rate coefficients and product ion branching ratios for all 66 reactions, measured at 298 K, are reported. The experimental reaction rate coefficients are compared to theoretically calculated collisional values. Several distinct reaction processes were observed among the large number of reactions studied, including charge transfer (non-dissociative and dissociative), abstraction, ion-molecule associations and, in the case of the reactions involving the reagent ion H₃O⁺, proton transfer. 

Direct Field Measurement of the Dynamic Amplification in a Bridge

In this paper, the level of dynamics, as described by the Assessment Dynamic Ratio (ADR), is measured directly through a field test on a bridge in the United Kingdom. The bridge was instrumented using fiber optic strain sensors and piezo-polymer weigh-in-motion sensors were installed in the pavement on the approach road. Field measurements of static and static-plus-dynamic strains were taken over 45 days. The results show that, while dynamic amplification is large for many loading events, these tend not to be the critical events. ADR, the allowance that should be made for dynamics in an assessment of safety, is small.

The Role of Corals on Tsunami Dynamics in an Island Setting: A Case Study of Tutuila Island

Thesis (Ph.D.)--University of Washington, 2016-07

Local Field Potentials in the Gustatory Cortex Carry Taste Information

It has been recently shownthat localfield potentials (LFPs)fromthe auditory and visual cortices carry information about sensory stimuli, but whether this is a universal property of sensory cortices remains to be determined. Moreover, little is known about the temporal dynamics of sensory information contained in LFPs following stimulus onset. Here we investigated the time course of the amount of stimulus information in LFPs and spikes from the gustatory cortex of awake rats subjected to tastants and water delivery on the tongue. We found that the phase and amplitude of multiple LFP frequencies carry information about stimuli, which have specific time courses after stimulus delivery. The information carried by LFP phase and amplitude was independent within frequency bands, since the joint information exhibited neither synergy nor redundancy. Tastant information in LFPs was also independent and had a different time course from the information carried by spikes. These findings support the hypothesis that the brain uses different frequency channels to dynamically code for multiple features of a stimulus.

Dynamics of coupled modes in two-section semiconductor lasers with saturable absorption

Semiconductor lasers have the potential to address a number of critical applications in advanced telecommunications and signal processing. These include applications that require pulsed output that can be obtained from self-pulsing and mode-locked states of two-section devices with saturable absorption. Many modern applications place stringent performance requirements on the laser source, and a thorough understanding of the physical mechanisms underlying these pulsed modes of operation is therefore highly desirable. In this thesis, we present experimental measurements and numerical simulations of a variety of self-pulsation phenomena in two-section semiconductor lasers with saturable absorption. Our theoretical and numerical results will be based on rate equations for the field intensities and the carrier densities in the two sections of the device, and we establish typical parameter ranges and assess the level of agreement with experiment that can be expected from our models. For each of the physical examples that we consider, our model parameters are consistent with the physical net gain and absorption of the studied devices. Following our introductory chapter, the first system that we consider is a two-section Fabry-Pérot laser. This example serves to introduce our method for obtaining model parameters from the measured material dispersion, and it also allows us to present a detailed discussion of the bifurcation structure that governs the appearance of selfpulsations in two-section devices. In the following two chapters, we present two distinct examples of experimental measurements from dual-mode two-section devices. In each case we have found that single mode self-pulsations evolve into complex coupled dualmode states following a characteristic series of bifurcations. We present optical and mode resolved power spectra as well as a series of characteristic intensity time traces illustrating this progression for each example. Using the results from our study of a twosection Fabry-Pérot device as a guide, we find physically appropriate model parameters that provide qualitative agreement with our experimental results. We highlight the role played by material dispersion and the underlying single mode self-pulsing orbits in determining the observed dynamics, and we use numerical continuation methods to provide a global picture of the governing bifurcation structure. In our concluding chapter we summarise our work, and we discuss how the presented results can inform the development of optimised mode-locked lasers for performance applications in integrated optics.

(Quantum-) critical dynamics of magnetoelectric materials

In this thesis the critical dynamics of several magnetoelectric compounds at their phase transition were examined. Mostly measurements of the dielectric properties in the frequency range of below 1 Hz up to 5 GHz were employed to evaluate the critical exponents for both magnetic field and temperature-dependent measurements. Most of the materials that are part of this work show anomalous behavior, especially at very low temperatures where quantum fluctuations are of the order of or even dominate those induced thermally. This anomalous behavior manifests in different forms. In Dy2Ti2O7 we demonstrate the existence of electric dipoles on magnetic monopoles. Here the dynamics at the critical endpoint located at 0.36K and in a magnetic field of 1T parallel to the [111] direction are of special interest. At this critical endpoint the expected critical slowing down of the dynamics could not only not be observed but instead the opposite, critical speeding-up by several orders of magnitude, could be demonstrated. Furthermore, we show that the phase diagram of Dy2Ti2O7 in this field direction can be reproduced solely from the dynamical properties, for example the resonance frequency of the observed relaxation that is connected to the monopole movement. Away from this point of the phase diagram the dynamics are slowing-down with reduction of temperature as one would expect. Additional measurements on Y2Ti2O7, a structurally identical but non-magnetic material, show only slowing down with reduction of temperature and no additional features. A possible explanation for the observed critical speeding-up is a coherent movement of magnetic monopoles close to the critical field that increases the resonance frequency by reducing the damping of the process. LiCuVO4 on the other hand behaves normally at its phase transition as long as the temperature is higher than 0.4 K. In this temperature regime the dynamics show critical slowing-down analogous to classical ferroelectric materials. This analogy extends also towards higher frequencies where the permittivity displays a ‘dispersion’ minimum that is temperature-dependent but of the order of 2 GHz. Below 0.4K the observed behavior changes drastically. Here we found no longer relaxational behavior but instead an excitation with very low energy. This low energy excitation was predicted by theory and is caused by nearly gapless soliton excitations within the 1D Cu2+ chains of LiCuVO4. Finally, in TbMnO3 the dynamics of the phase transition into the multiferroic phase was observed at roughly 27 K, a much higher temperature compared to the other materials. Here the expected critical slowing-down was observed, even though in low-frequency measurements this transition into the ferroelectric phase is overshadowed by the so-called c-axis relaxation. Therefore, only frequencies above 1MHz could be used to determine the critical exponents for both temperatureand magnetic-field-dependent measurements. This was done for both the peak frequency as well as the relaxation strength. In TbMnO3 an electromagnetic soft-mode with small optical weight causes the observed fluctuations, similar to the case of multiferroic MnWO4.

Enhancing the numerical performance of fire field models (CMS Paper)

This paper describes two new techniques designed to enhance the performance of fire field modelling software. The two techniques are "group solvers" and automated dynamic control of the solution process, both of which are currently under development within the SMARTFIRE Computational Fluid Dynamics environment. The "group solver" is a derivation of common solver techniques used to obtain numerical solutions to the algebraic equations associated with fire field modelling. The purpose of "group solvers" is to reduce the computational overheads associated with traditional numerical solvers typically used in fire field modelling applications. In an example, discussed in this paper, the group solver is shown to provide a 37% saving in computational time compared with a traditional solver. The second technique is the automated dynamic control of the solution process, which is achieved through the use of artificial intelligence techniques. This is designed to improve the convergence capabilities of the software while further decreasing the computational overheads. The technique automatically controls solver relaxation using an integrated production rule engine with a blackboard to monitor and implement the required control changes during solution processing. Initial results for a two-dimensional fire simulation are presented that demonstrate the potential for considerable savings in simulation run-times when compared with control sets from various sources. Furthermore, the results demonstrate the potential for enhanced solution reliability due to obtaining acceptable convergence within each time step, unlike some of the comparison simulations.

Waves, bumps, and patterns in neural field theories

Neural field models of firing rate activity have had a major impact in helping to develop an understanding of the dynamics seen in brain slice preparations. These models typically take the form of integro-differential equations. Their non-local nature has led to the development of a set of analytical and numerical tools for the study of waves, bumps and patterns, based around natural extensions of those used for local differential equation models. In this paper we present a review of such techniques and show how recent advances have opened the way for future studies of neural fields in both one and two dimensions that can incorporate realistic forms of axo-dendritic interactions and the slow intrinsic currents that underlie bursting behaviour in single neurons.

Continuum tensor network field states, path integral representations and spatial symmetries

A natural way to generalize tensor network variational classes to quantum field systems is via a continuous tensor contraction. This approach is first illustrated for the class of quantum field states known as continuous matrix-product states (cMPS). As a simple example of the path-integral representation we show that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. A completeness argument is also provided that shows that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity. Beyond this, we obtain a well-behaved field limit of projected entangled-pair states (PEPS) in two dimensions that provide an abstract class of quantum field states with natural symmetries. We demonstrate how symmetries of the physical field state are encoded within the dynamics of an auxiliary field system of one dimension less. In particular, the imposition of Euclidean symmetries on the physical system requires that the auxiliary system involved in the class' definition must be Lorentz-invariant. The physical field states automatically inherit entropy area laws from the PEPS class, and are fully described by the dissipative dynamics of a lower dimensional virtual field system. Our results lie at the intersection many-body physics, quantum field theory and quantum information theory, and facilitate future exchanges of ideas and insights between these disciplines.

Dynamics of a large submersed plant bed in upper Chesapeake Bay

A large SAV bed in upper Chesapeake Bay has experienced several abrupt shifts over the past half-century, beginning with near-complete loss after a record-breaking flood in 1972, followed by an unexpected, rapid resurgence in the early 2000’s, then partial decline in 2011 following another major flood event. Together, these trends and events provide a unique opportunity to study a recovering SAV ecosystem from several different perspectives. First, I analyzed and synthesized existing time series datasets to make inferences about what factors prompted the recovery. Next, I analyzed existing datasets, together with field samples and a simple hydrodynamic model to investigate mechanisms of SAV bed loss and resilience to storm events. Finally, I conducted field deployments and experiments to explore how the bed affects internal physical and biogeochemical processes and what implications those effects have for the dynamics of the system. I found that modest reductions in nutrient loading, coupled with several consecutive dry years likely facilitated the SAV resurgence. Furthermore, positive feedback processes may have played a role in the sudden nature of the recovery because they could have reinforced the state of the bed before and after the abrupt shift. I also found that scour and poor water clarity associated with sediment deposition during the 2011 flood event were mechanisms of plant loss. However, interactions between the bed, water flow, and waves served as mechanisms of resilience because these processes created favorable growing conditions (i.e., clear water, low flow velocities) in the inner core of the bed. Finally, I found that that interactions between physical and biogeochemical processes led to low nutrient concentrations inside the bed relative to outside the bed, which created conditions that precluded algal growth and reinforced vascular plant dominance. This work demonstrates that positive feedbacks play a central role in SAV resilience to both chronic eutrophication as well as acute storm events. Furthermore, I show that analysis of long-term ecological monitoring data, together with field measurements and experiments, can be an effective approach for understanding the mechanisms underlying ecosystem dynamics.

Variational ensemble Kalman filtering applied to data assimilation problems in computational fluid dynamics

One challenge on data assimilation (DA) methods is how the error covariance for the model state is computed. Ensemble methods have been proposed for producing error covariance estimates, as error is propagated in time using the non-linear model. Variational methods, on the other hand, use the concepts of control theory, whereby the state estimate is optimized from both the background and the measurements. Numerical optimization schemes are applied which solve the problem of memory storage and huge matrix inversion needed by classical Kalman filter methods. Variational Ensemble Kalman filter (VEnKF), as a method inspired the Variational Kalman Filter (VKF), enjoys the benefits from both ensemble methods and variational methods. It avoids filter inbreeding problems which emerge when the ensemble spread underestimates the true error covariance. In VEnKF this is tackled by resampling the ensemble every time measurements are available. One advantage of VEnKF over VKF is that it needs neither tangent linear code nor adjoint code. In this thesis, VEnKF has been applied to a two-dimensional shallow water model simulating a dam-break experiment. The model is a public code with water height measurements recorded in seven stations along the 21:2 m long 1:4 m wide flume’s mid-line. Because the data were too sparse to assimilate the 30 171 model state vector, we chose to interpolate the data both in time and in space. The results of the assimilation were compared with that of a pure simulation. We have found that the results revealed by the VEnKF were more realistic, without numerical artifacts present in the pure simulation. Creating a wrapper code for a model and DA scheme might be challenging, especially when the two were designed independently or are poorly documented. In this thesis we have presented a non-intrusive approach of coupling the model and a DA scheme. An external program is used to send and receive information between the model and DA procedure using files. The advantage of this method is that the model code changes needed are minimal, only a few lines which facilitate input and output. Apart from being simple to coupling, the approach can be employed even if the two were written in different programming languages, because the communication is not through code. The non-intrusive approach is made to accommodate parallel computing by just telling the control program to wait until all the processes have ended before the DA procedure is invoked. It is worth mentioning the overhead increase caused by the approach, as at every assimilation cycle both the model and the DA procedure have to be initialized. Nonetheless, the method can be an ideal approach for a benchmark platform in testing DA methods. The non-intrusive VEnKF has been applied to a multi-purpose hydrodynamic model COHERENS to assimilate Total Suspended Matter (TSM) in lake Säkylän Pyhäjärvi. The lake has an area of 154 km2 with an average depth of 5:4 m. Turbidity and chlorophyll-a concentrations from MERIS satellite images for 7 days between May 16 and July 6 2009 were available. The effect of the organic matter has been computationally eliminated to obtain TSM data. Because of computational demands from both COHERENS and VEnKF, we have chosen to use 1 km grid resolution. The results of the VEnKF have been compared with the measurements recorded at an automatic station located at the North-Western part of the lake. However, due to TSM data sparsity in both time and space, it could not be well matched. The use of multiple automatic stations with real time data is important to elude the time sparsity problem. With DA, this will help in better understanding the environmental hazard variables for instance. We have found that using a very high ensemble size does not necessarily improve the results, because there is a limit whereby additional ensemble members add very little to the performance. Successful implementation of the non-intrusive VEnKF and the ensemble size limit for performance leads to an emerging area of Reduced Order Modeling (ROM). To save computational resources, running full-blown model in ROM is avoided. When the ROM is applied with the non-intrusive DA approach, it might result in a cheaper algorithm that will relax computation challenges existing in the field of modelling and DA.

Shape and chemically anisotropic colloidal particles: A theoretical study of their structure, phase behavior, and slow dynamics

A detailed non-equilibrium state diagram of shape-anisotropic particle fluids is constructed. The effects of particle shape are explored using Naive Mode Coupling Theory (NMCT), and a single particle Non-linear Langevin Equation (NLE) theory. The dynamical behavior of non-ergodic fluids are discussed. We employ a rotationally frozen approach to NMCT in order to determine a transition to center of mass (translational) localization. Both ideal and kinetic glass transitions are found to be highly shape dependent, and uniformly increase with particle dimensionality. The glass transition volume fraction of quasi 1- and 2- dimensional particles fall monotonically with the number of sites (aspect ratio), while 3-dimensional particles display a non-monotonic dependence of glassy vitrification on the number of sites. Introducing interparticle attractions results in a far more complex state diagram. The ideal non-ergodic boundary shows a glass-fluid-gel re-entrance previously predicted for spherical particle fluids. The non-ergodic region of the state diagram presents qualitatively different dynamics in different regimes. They are qualified by the different behaviors of the NLE dynamic free energy. The caging dominated, repulsive glass regime is characterized by long localization lengths and barrier locations, dictated by repulsive hard core interactions, while the bonding dominated gel region has short localization lengths (commensurate with the attraction range), and barrier locations. There exists a small region of the state diagram which is qualified by both glassy and gel localization lengths in the dynamic free energy. A much larger (high volume fraction, and high attraction strength) region of phase space is characterized by short gel-like localization lengths, and long barrier locations. The region is called the attractive glass and represents a 2-step relaxation process whereby a particle first breaks attractive physical bonds, and then escapes its topological cage. The dynamic fragility of fluids are highly particle shape dependent. It increases with particle dimensionality and falls with aspect ratio for quasi 1- and 2- dimentional particles. An ultralocal limit analysis of the NLE theory predicts universalities in the behavior of relaxation times, and elastic moduli. The equlibrium phase diagram of chemically anisotropic Janus spheres and Janus rods are calculated employing a mean field Random Phase Approximation. The calculations for Janus rods are corroborated by the full liquid state Reference Interaction Site Model theory. The Janus particles consist of attractive and repulsive regions. Both rods and spheres display rich phase behavior. The phase diagrams of these systems display fluid, macrophase separated, attraction driven microphase separated, repulsion driven microphase separated and crystalline regimes. Macrophase separation is predicted in highly attractive low volume fraction systems. Attraction driven microphase separation is charaterized by long length scale divergences, where the ordering length scale determines the microphase ordered structures. The ordering length scale of repulsion driven microphase separation is determined by the repulsive range. At the high volume fractions, particles forgo the enthalpic considerations of attractions and repulsions to satisfy hard core constraints and maximize vibrational entropy. This results in site length scale ordering in rods, and the sphere length scale ordering in Janus spheres, i.e., crystallization. A change in the Janus balance of both rods and spheres results in quantitative changes in spinodal temperatures and the position of phase boundaries. However, a change in the block sequence of Janus rods causes qualitative changes in the type of microphase ordered state, and induces prominent features (such as the Lifshitz point) in the phase diagrams of these systems. A detailed study of the number of nearest neighbors in Janus rod systems reflect a deep connection between this local measure of structure, and the structure factor which represents the most global measure of order.

Exotic dynamics in a firing rate model of neural tissue with threshold accommodation

Many of the equations describing the dynamics of neural systems are written in terms of firing rate functions, which themselves are often taken to be threshold functions of synaptic activity. Dating back to work by Hill in 1936 it has been recognized that more realistic models of neural tissue can be obtained with the introduction of state-dependent dynamic thresholds. In this paper we treat a specific phenomenological model of threshold accommodation that mimics many of the properties originally described by Hill. Importantly we explore the consequences of this dynamic threshold at the tissue level, by modifying a standard neural field model of Wilson-Cowan type. As in the case without threshold accommodation classical Mexican-Hat connectivity is shown to allow for the existence of spatially localized states (bumps) in both one and two dimensions. Importantly an analysis of bump stability in one dimension, using recent Evans function techniques, shows that bumps may undergo instabilities leading to the emergence of both breathers and traveling waves. Moreover, a similar analysis for traveling pulses leads to the conditions necessary to observe a stable traveling breather. In the regime where a bump solution does not exist direct numerical simulations show the possibility of self-replicating bumps via a form of bump splitting. Simulations in two space dimensions show analogous localized and traveling solutions to those seen in one dimension. Indeed dynamical behavior in this neural model appears reminiscent of that seen in other dissipative systems that support localized structures, and in particular those of coupled cubic complex Ginzburg-Landau equations. Further numerical explorations illustrate that the traveling pulses in this model exhibit particle like properties, similar to those of dispersive solitons observed in some three component reaction-diffusion systems. A preliminary account of this work first appeared in S Coombes and M R Owen, Bumps, breathers, and waves in a neural network with spike frequency adaptation, Physical Review Letters 94 (2005), 148102(1-4).