## You are here

# Applied Dynamics Seminar Series

# Applied Dynamics Seminar Series

## Thursdays, 12:30 p.m.

## IREAP Large Conference Room (ERF 1207)

Subscribe to our mailing list for announcements sending an email to __listserv@listserv.umd.edu__ with no subject, and **SUBSCRIBE APPLIED-DYNAMICS [Your full name] **or **SUBSCRIBE APPLIED-DYNAMICS [Anonymous]** (no square brackets) in the body of the email. An email signature might prevent your subsription command from working.

### September 7, 2017

#### Ultrafast Large-scale Neural Network Processor on a Chip

Daniel Lathrop and Itamar Shani

University of Maryland | Department of Physics / IREAP

*Abstract: Neural networks allow machines to imitate the way in which human intelligence solves problems by inferring from past experience. These networks are composed of large arrays of communicating neurons, each one performing a simple non-linear operation. When combined and trained by variation of connection weights, the network can perform complex perceptive computational tasks such as image and voice recognition and complex pattern predictions. When implementing neural networks on conventional digital processing hardware such as those at the core of our PCs, an immense inefficiency stands out: neural network computations are inherently parallel, while computers were designed to perform computations serially. This leads to slow computation times and a high toll of energy consumption. Here we report a way to overcome this challenge. We implement a silicon chip with thousands, and potentially millions, of processing interconnected ‘neurons’, each one operating at 200ps rates. The chip is, thus, capable of performing fully parallel, highly efficient computation. For the design of our network, we follow the well-established machine learning algorithms in which interconnections are described by a random sparse directed graph. We show preliminary laboratory measurements of the network dynamics on a chip and discuss its software variations. *

### September 14, 2017

#### Synchronization of Quantum Dipoles

University of Colorado Boulder | Department of Applied Mathematics

*Abstract: In this talk we discuss the emergence of synchronization in arrays of quantum radiating dipoles coupled only via anisotropic and long-range dipolar interactions. It is found that in the presence of an incoherent energy source, dipolar interactions can lead to a resilient synchronized steady-state. A classical mean-field description of the model results in equations similar to the classical Kuramoto model for synchronization of phase oscillators. Using the mean-field formulation for the all-to-all coupled case, the synchronized state can be studied, and it is found that it exists only for a finite range of the external energy source rates. Results obtained from the mean-field model are compared with numerical simulations of the quantum system and it is found that synchronization is robust to quantum fluctuations and spatially decaying coupling. Additional nonstationary synchronization patterns and bistability are discussed.*

### September 21, 2017

#### Basin Entropy: A Measure of the Final State Unpredictability and Applications to Some Physical Systems

King Juan Carlos University | Department of Physics

*Abstract: In nonlinear dynamics, basins of attraction link a given set of initial conditions to its corresponding final states. This notion appears in a broad range of applications where several outcomes are possible, which is a common situation in neuroscience, economy, astronomy, ecology and many other disciplines. Depending on the nature of the basins, prediction can be difficult even in systems that evolve under deterministic rules. From this respect, a proper classification of this unpredictability is clearly required. To address this issue, we introduce the basin entropy, a measure to quantify this uncertainty. Its application is illustrated with several paradigmatic examples that allow us to identify the ingredients that hinder the prediction of the final state. The basin entropy provides an efficient method to probe the behavior of a system when different parameters are varied.Additionally, we provide a sufficient condition for the existence of fractal basin boundaries: when the basin entropy of the boundaries is larger than log 2, the basin is fractal. These ideas have been applied to some physical systems such as experiments of chaotic scattering of cold atoms, models of shadows of binary black holes, and classical and relativistic chaotic scattering associated to the Hénon-Heiles Hamiltonian system in astrophysics.*

### September 28, 2017

#### A Model-Free Machine Learning Technique for Studying High Dimensional Spatiotemporal Chaos

University of Maryland | IREAP

*Abstract: Networks of nonlinearly interacting neuron-like units have the capacity to approximately reproduce the dynamical behavior of a wide variety of dynamical systems. We demonstrate the use of such neural networks for reconstruction of chaotic attractors from limited time series data using a machine learning technique known as reservoir computing. The orbits of the reconstructed attractor can be used to obtain approximate estimates of the ergodic properties of the original system. As a specific example, we focus on the task of determining the Lyapunov exponents of a system from limited time series data. Using the example of the Kuramoto-Sivashinsky system, we show that this technique offers a robust estimate of a large number of Lyapunov exponents of a high dimensional spatiotemporal chaotic system. We further develop an effective, computationally parallelizable technique for model-free prediction of spatiotemporal chaotic systems of arbitrarily large spatial extent and dimension purely from observations of the system's past evolution.*

### October 5, 2017

#### Complexity and Self Organization in Superconducting Metamaterials

University of Crete | Department of Physics

*Abstract: Macroscopic quantum devices are becoming reality not only for computational purposes but also as sensors and for other general applications In this talk we will focus on superconducting technology and analyze the emergence of coherence in coupled networks of meta-atoms made of units such as SQUIDS and Josephson junctions. These networks may operate classically in a negative permeability regime[1], induce intrinsic nonlinear localized modes and partial coherence in the form of chimeras[2], tame disorder through hysteretic loops or transmit through nonlinear frequency bands. In the quantum regime, on the other hand, meta-atoms may interact through injected electromagnetic fields and form propagating “quantum breathers”, i.e. compound semi-classical propagating modes induced by the nonlinearity of the qubit-field interaction [3]. These coherent modes generate self-induced transparency in the medium and in certain cases may also induce super-radiance. [1] N. Lazarides and G. P. Tsironis, rf SQUID metamaterials, Appl. Phys. Lett. 90, 163501 (2007). [2] N. Lazarides, G. Neofotistos, and G. P. Tsironis, Chimeras in SQUID metamaterials, Physical Review B 91, 054303 (2015). [3] Z. Ivic, N. Lazarides, and G. P. Tsironis, Qubit lattice coherence induced by electromagnetic pulses in superconducting metamaterials, Scientific Reports 6, 29374(2016).*

### October 12, 2017

#### Granular Dynamics in Low Gravity

University of Maryland | Department of Astronomy

*Abstract: Small solar system bodies are generally covered in layers of particulate "regolith" with largely unknown properties. The effective gravities on these bodies can sometimes tend to zero at the equator, and have been measured to be negative in a few cases. Space agencies and commercial enterprises show increasing interest in visiting, landing on, and sampling from such bodies, so it is important to understand how the regolith will respond to intrusion. Due to the difficulty and expense of carrying out experiments in low gravity, we turn to computer simulations of granular dynamics to provide insight into the conditions that missions to other worlds may encounter and to help interpret observations of these bodies. As high-end computing resources become more readily available, granular dynamics simulations have become more sophisticated, treating particle collisions as finite-duration, multi-contact events with explicit friction forces between irregularly shaped grains subject to external forces, boundary conditions, and cohesion. It is imperative that such simulation methods be calibrated and validated at laboratory scales to give confidence in their application to other environments. Here I review our group's approach to simulating granular dynamics in low gravity, with examples that include impacts into granular materials, vibration-induced segregation, landslides, models of samplers and landers, and, on larger scales, simulations of entire granular bodies (rubble piles) and planetary rings.*

### Special Seminar: FRIDAY October 13, 2017 at 3:00pm

#### Statistical Description of Mixed Systems (Chaotic and Regular), Correlations and "Thermalization"

Technion | Department of Physics

*Abstract: We discuss a statistical theory for Hamiltonian dynamics with a mixed phase space, where in some parts of phase space the dynamics is chaotic while in other parts it is regular. Transport in phase space is dominated by sticking to complicated structures and its distribution is universal. The survival probability in the vicinity of the initial point is a power law in time with a universal exponent. We calculate this exponent in the framework of the Markov Tree model proposed by Meiss and Ott in 1986. It turns out that, inspite of many approximations, it predicts important results quantitatively. The calculations are extended to the quantum regime where correlation functions and observables are studied. The seminar will be very informal and some work still in progress will be reported. The work reported is in collaboration with Or Alus, James Meiss and Mark Srednicki.*

### October 19, 2017

#### Ultra-High Intensity Laser Physics and Applications

University of Maryland | Department of Electrical and Computer Engineering / Department of Physics

*Abstract: This talk will cover ultra-high field physics phenomena and interactions associated with high intensity short pulse lasers. The operating parameter regime of these lasers covers a wide range, e.g., peak powers of ~ 1012 - 1015 W, pulse lengths of ~ 10-12 – 10-14 sec, and intensities of ~ 1014 - 1023 W/cm2. These lasers are used in high-field physics research and have a number of unique applications. The physical processes associated with USPL interactions and propagation include: photo ionization, Kerr and Raman effects, self-phase modulation, optical and plasma filamentation, optical shocks, frequency chirping, propagation through atmospheric turbulence, etc. This talk will discuss these interrelated physical processes and some unique applications, such as: laser driven acceleration, UV/X-ray generation, underwater acoustic sources, atmospheric spark formation, detection of radioactive materials and atmospheric lasing.*

### October 26, 2017

#### Neural Representation of Speech in Human Auditory Cortex

University of Maryland | Department of Electrical and Computer Engineering

*Abstract: We investigate how continuous speech is represented in human auditory cortex. We use magnetoencephalography (MEG) to record the neural responses of listeners to natural, continuous speech, in a variety of auditory scenes. Systems analysis, which quantitatively compares a speech signal to its evoked cortical responses, allows us to determine the cortical representations of the speech. Interestingly, the cortical representation allows the time-varying envelope of the speech to be reconstructed from the observed neural response to the speech. We find that cortical representations of continuous speech are very robust to interference from competing speakers, and many other kinds of noise, consistent with our ability to understand speech even in a noisy room (the "Cocktail Party" problem). Indeed, individual neural representations of the speech of both the foreground and background speaker are observed, with each being selectively time-locked to the rhythm of the corresponding speech, but the with the foreground speech represented more faithfully than the background.*

### November 2, 2017

#### Modeling the Network Dynamics of Pulse-Coupled Neurons

University of Maryland | IREAP

*Abstract: Computer modeling of neural dynamics is an important component of the long-term goal of understanding the brain. A barrier to such modeling is the practical limit on computer resources given the enormous number of neurons in the human brain (about 10^11.) Our work addresses this problem by developing a method for obtaining low dimensional macroscopic descriptions for functional groups consisting of many neurons. Specifically, we formulate a mean-field approximation to investigate macroscopic network effects on the dynamics of large systems of pulse-coupled neurons and derive a reduced system of ordinary differential equations describing the dynamics. We find that solutions of the reduced system agree with those of the full network. This dimensional reduction allows for more efficient characterization of system phase transitions and attractors. Our results show the utility of these dimensional reduction techniques for analyzing the effects of network topology on macroscopic behavior in neuronal networks.*

* *

### November 9, 2017

#### Optimal Bounds and Extremal Trajectories for Time Averages in Dynamical Systems

University of Michigan | Mathematics / Center for the Study of Complex Systems

*Abstract: For any quantity of interest in a dynamical system governed by ordinary differential equations it is natural to seek the largest (or smallest) long-time average among solution trajectories. Upper bounds can be proved a priori using auxiliary functions, the optimal choice of which is a convex optimization. The problems of finding maximal trajectories and minimal auxiliary functions are in fact strongly dual so auxiliary functions can produce arbitrarily sharp upper bounds on maximal time averages. They also define volumes in phase space where maximal trajectories must lie. For polynomial equations of motion auxiliary functions can be constructed by semidefinite programming, which we illustrate using the Lorenz system. This is joint work with Ian Tobasco and David Goluskin.*

* *

### Special Seminar: TUESDAY November 14, 2017 at 11:00am

#### Sturm-Liouville Framework for Dynamical Reconstruction by Delay Embedding

Kyoto University | Department of Mathematics

*Abstract: Delay embedding is well-known for non-linear time-series analysis, and it is used in several research fields such as physics, informatics, neuroscience and so on. The celebrated theorem of Takens ensures validity of the delay embedding analysis: embedded data preserves topological properties, which the original dynamics possesses, if one embeds into some phase space with sufficiently high dimension. This means that, for example, an attractor can be reconstructed by the delay coordinate system topologically. However, configuration of an embedded dataset may easily vary with the delay width and the delay dimension, namely, ``the way of embedding". In a practical sense, this sensitivity may cause degradation of reliability of the method, therefore it is natural to require robustness of the result obtained by the embedding method in certain sense. In this study, we investigate the mathematical structure of the framework of delay-embedding analysis to provide Ansatz to choose the appropriate way of embedding, in order to utilize for time-series prediction. In short, mathematical theories of the Hilbert-Schmidt integral operator and the corresponding Sturm-Liouville eigenvalue problem underlie the framework. Using these mathematical theories, one can derive error estimates of mode decomposition obtained by the present method and can obtain the phase-space reconstruction by using the leading modes of the decomposition. In this talk, we will show some results for some numerical and experimental datasets to validate the present method.*

### November 16, 2017

#### Adaptive Coding for Sensory Inference in Dynamic Environments

HHMI Janelia Research Campus

*Abstract: Making reliable inferences about the environment is crucial for survival. In order to escape a hawk, for example, a mouse might need to infer the hawk’s position and velocity from patterns of light that fall on its retina. Such inferences require large ensembles of sensory neurons whose activity is metabolically expensive. A growing body of evidence suggests that sensory systems reduce metabolic costs by limiting the fidelity with which some stimuli are encoded in neural responses. Limited coding fidelity, however, can lead to inaccuracies in inference. Here, we derive a framework for dynamically balancing the cost of encoding with the error that encoding can induce in inference. We model a system that must use minimal metabolic resources to maintain an accurate estimate of a nonstationary environment, and we show that the optimal system should adapt the fidelity with which stimuli are encoded in neural responses based on a changing estimate of the environment. We use this framework to illustrate how a range of neuronal and behavioral phenomena can be understood as signatures of adaptive encoding for accurate inference.*

### November 23, 2017

#### Thanksgiving Break - No Seminar

### November 30, 2017

#### No Equations, No Variables, No Parameters, No Space, No Time: Data and the Modeling of Complex Systems

Johns Hopkins Univeristy | Department Chemical and Biomolecular Engineering / Department of Applied Mathematics and Statistics

*Abstract: Obtaining predictive dynamical equations from data lies at the heart of science and engineering modeling, and is the linchpin of our technology. In mathematical modeling one typically progresses from observations of the world (and some serious thinking!) first to equations for a model, and then to the analysis of the model to make predictions. Good mathematical models give good predictions (and inaccurate ones do not) - but the computational tools for analyzing them are the same: algorithms that are typically based on closed form equations. While the skeleton of the process remains the same, today we witness the development of mathematical techniques that operate directly on observations -data-, and appear to circumvent the serious thinking that goes into selecting variables and parameters and deriving accurate equations. The process then may appear to the user a little like making predictions by "looking in a crystal ball". Yet the "serious thinking" is still there and uses the same -and some new- mathematics: it goes into building algorithms that "jump directly" from data to the analysis of the model (which is now not available in closed form) so as to make predictions. Our work here presents a couple of efforts that illustrate this ``new” path from data to predictions. It really is the same old path, but it is travelled by new means.*

### December 7, 2017

#### Optical Control of Excitation Waves in A Biological Excitable Medium

McGill University | Department of Physiology

*Abstract: Macroscopic excitation waves are found in a diverse range of settings including chemical reactions and the heart and brain. In the case of living biological tissue, the spatiotemporal patterns formed by these excitation waves are different in healthy and diseased states. Current electrical and pharmacological methods for wave modulation lack the spatiotemporal precision needed to control these patterns. Optical methods have the potential to overcome these limitations, but until recently have only been demonstrated in simple systems. Here I discuss results using a new dye-free optical imaging modality with optogenetic actuation to achieve dynamic control of cardiac excitation waves. Illumination with patterned light is demonstrated to optically control the direction, speed, and spiral chirality of such waves in cardiac tissue. This all-optical approach offers a new experimental platform for the study and control of pattern formation in complex biological excitable systems.*

* *