Multistep variable methods for exact integration of perturbed stiff linear systems

A family of real and analytical functions with values within the ring of M(m, R) is introduced. The solution for linear systems of differential equations is expressed as a series of Φ-functions. This new multistep method is defined for variable-step and variable-order, maintains the good properties of the Φ-function series method. It incorporates to compute the coefficients of the algorithm a recurrent algebraic procedure, based in the existing relation between the divided differences and the elemental and complete symmetrical functions. In addition, under certain hypotheses, the new multistep method calculates the exact solution of the perturbed problem. The new method is implemented in a computational algorithm which enables us to resolve in a general manner some physics and engineering IVP’s modeled by means systems of differential equations. The good behaviour and precision of the method is evidenced by contrasting the results with other-reputed algorithms and even with methods based on Scheifele’s G-functions.This work has been supported by GRE09-13 project of the University of Alicante and the project of the Generalitat Valenciana GV/2011/032

Dynamics of Structured Systems

[no abstract available

Non-reciprocal phase transitions

Out of equilibrium, the lack of reciprocity is the rule rather than the exception. Non-reciprocal interactions occur, for instance, in networks of neurons, directional growth of interfaces, and synthetic active materials. While wave propagation in non-reciprocal media has recently been under intense study, less is known about the consequences of non-reciprocity on the collective behavior of many-body systems. Here, we show that non-reciprocity leads to time-dependent phases where spontaneously broken symmetries are dynamically restored. The resulting phase transitions are controlled by spectral singularities called exceptional points. We describe the emergence of these phases using insights from bifurcation theory and non-Hermitian quantum mechanics. Our approach captures non-reciprocal generalizations of three archetypal classes of self-organization out of equilibrium: synchronization, flocking and pattern formation. Collective phenomena in these non-reciprocal systems range from active time-(quasi)crystals to exceptional-point enforced pattern-formation and hysteresis. Our work paves the way towards a general theory of critical phenomena in non-reciprocal matter.Comment: Supplementary movies at https://home.uchicago.edu/~vitelli/videos.htm

Impact of Adaptation Currents on Synchronization of Coupled Exponential Integrate-and-Fire Neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies

A Trigonometrically Fitted Block Method for Solving Oscillatory Second-Order Initial Value Problems and Hamiltonian Systems

In this paper, we present a block hybrid trigonometrically fitted Runge-Kutta-Nyström method (BHTRKNM), whose coefficients are functions of the frequency and the step-size for directly solving general second-order initial value problems (IVPs), including Hamiltonian systems such as the energy conserving equations and systems arising from the semidiscretization of partial differential equations (PDEs). Four discrete hybrid formulas used to formulate the BHTRKNM are provided by a continuous one-step hybrid trigonometrically fitted method with an off-grid point. We implement BHTRKNM in a block-by-block fashion; in this way, the method does not suffer from the disadvantages of requiring starting values and predictors which are inherent in predictor-corrector methods. The stability property of the BHTRKNM is discussed and the performance of the method is demonstrated on some numerical examples to show accuracy and efficiency advantages

Theoretical study and experimental implementation of an ion-nanowire hybrid system

Numerous efforts are currently aiming at the implementation of hybrid quantum systems in the quest for developing new devices for quantum technologies. The motivation for a solid state-atomic interface emerges from the possibility to couple otherwise incompatible platforms in order to combine the advantages of both in a complementary fashion. Under this scope, we report in this thesis the implementation of a new quantum device for the interface of an ultracold trapped ion and a conductive nanowire. The experimental measurements obtained here seem promising for the realization of the hybrid dynamics of the ion-nanowire hybrid system. Additionally, we performed numerical calculations in order to characterize the perturbation of the trapping potential for the ion by the nanowire. From our classical and quantum dynamics calculations we explored possibilities of generating Gaussian and non-Gaussian quantum states of the ion from the nanowire drive. Our modelling indicated, also, that sympathetic cooling and quantum entanglement can be realized when both subsystems operate in the quantum regime. The present ion-nanowire hybrid system might prove promising as a new quantum device for quantum information and quantum sensing experiments, for spectroscopy and for mass spectrometry

Activities of the Institute for Computer Applications in Science and Engineering

Research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, numerical analysis, and computer science during the period April 1, 1985 through October 2, 1985 is summarized