Stochastic point processes with refractoriness appear frequently in the quantitative analysis of physical and biological systems, such as the generation of action potentials by nerve cells, the release and reuptake of vesicles at a synapse, and the counting of particles by detector devices. Here we present an extension of renewal theory to describe ensembles of point processes with time varying input. This is made possible by a representation in terms of occupation numbers of two states: Active and refractory. The dynamics of these occupation numbers follows a distributed delay differential equation. In particular, our theory enables us to uncover the effect of refractoriness on the time-dependent rate of an ensemble of encoding point processes in response to modulation of the input. We present exact solutions that demonstrate generic features, such as stochastic transients and oscillations in the step response as well as resonances, phase jumps and frequency doubling in the transfer of periodic signals. We show that a large class of renewal processes can indeed be regarded as special cases of the model we analyze. Hence our approach represents a widely applicable framework to define and analyze non-stationary renewal processes.

Moritz Deger Moritz Helias Stefano Cardanobile Fatihcan M. Atay Stefan Rotter

The continuous integration of experimental data into coherent models of the brain is an increasing challenge of modern neuroscience. Such models provide a bridge between structure and activity, and identify the mechanisms giving rise to experimental observations. Nevertheless, structurally realistic network models of spiking neurons are necessarily underconstrained even if experimental data on brain connectivity are incorporated to the best of our knowledge. Guided by physiological observations, any model must therefore explore the parameter ranges within the uncertainty of the data. Based on simulation results alone, however, the mechanisms underlying stable and physiologically realistic activity often remain obscure. We here employ a mean-field reduction of the dynamics, which allows us to include activity constraints into the process of model construction. We shape the phase space of a multi-scale network model of the vision-related areas of macaque cortex by systematically refining its connectivity. Fundamental constraints on the activity, i.e., prohibiting quiescence and requiring global stability, prove sufficient to obtain realistic layer- and area-specific activity. Only small adaptations of the structure are required, showing that the network operates close to an instability. The procedure identifies components of the network critical to its collective dynamics and creates hypotheses for structural data and future experiments. The method can be applied to networks involving any neuron model with a known gain function.

Jannis Schuecker Maximilian Schmidt Sacha J. van Albada Markus Diesmann Moritz Helias

The ongoing development of electromagnets based on High Temperature Superconductors has led to the conceptual exploration of high-magnetic-field fusion reactors of the tokamak type, operating at on-axis fields above 10 T. In this work we explore the consequences of the potential future availability of high-field three-dimensional electromagnets on the physics design point of a stellarator reactor. We find that, when an increase in the magnetic field strength $B$ is used to maximally reduce the device linear size $R\sim B^{-4/3}$ (with otherwise fixed magnetic geometry), the physics design point is largely independent of the chosen field strength/device size. A similar degree of optimization is to be imposed on the magnetohydrodynamic, transport and fast ion confinement properties of the magnetic configuration of that family of reactor design points. Additionally, we show that the family shares an invariant operation map of fusion power output as a function of the auxiliary power and relative density variation. The effects of magnetic field over-engineering and the $R(B)$ scaling of design points with constant neutron wall loading are also inspected. In this study we use geometric parameters characteristic of the helias reactor, but most results apply to other stellarator configurations

J. A. Alonso I. Calvo D. Carralero J. L. Velasco J. M. García-Regaña I. Palermo D. Rapisarda

Advanced stellarators require convoluted modular coils to produce a plasma with satisfactory performance. Moreover, the number of coils is sometimes high to decrease the modular ripple created by the coils. For reactor stellarators, these requirements imply relatively small ports for in-vessel access and maintenance, i.e. in comparison with tokamaks. The blankets and divertor modules will have to be replaced periodically (about each 1-4 years depending on the design) due to neutron damage, and also erosion of divertor targets. Blanket modules are activated, thus, all the maintenance operations have to be produced remotely. In order to reduce the shutdown time and cost during component replacement, and to reduce the number, speed and other specifications of the remote maintenance equipment, the number of blanket modules in the reactor should be low and thus, the blanket modules should be large (in relation to the minor and major radius). Nevertheless, the size of the openings between coils limits the maximum size of the blanket and divertor modules, though several potential enhancements have been proposed in the past for stellarators, like straightening the outboard segments of the coils and the movement and/or expansion of certain coils to have wider access. The present work reports on a coil geometry for the 'Helias Stellarator Reactor' (HSR) of three periods (HSR3) with coils located far from the plasma at the outboard region of the straight-like sector. This feature creates natural wide openings at such regions of the coils, which may be utilized to allow access to large blanket and divertor modules.

V. Queral V. Tribaldos J. M. Reynolds I. Fernandez

After dimensional reduction the stationary spherically symmetric sector of Einstein's gravity is identified with an SL(2,R)/SO(2) Sigma model coupled to a one dimensional gravitational remnant. The space of classical solutions consists of a one parameter family interpolating between the Schwarzschild and the Taub-NUT solution. A Dirac Quantization of this system is performed and the observables -- the Schwarzschild mass and the Taub-NUT charge operator -- are shown to be self-adjoint operators with a continuous spectrum ranging from $-\infty$ to $\infty$. The Hilbert space is constructed explicitely using a harmonic space approach.

Helia Hollmann

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, turbulent transport poses a serious limitation to the global performance of the machine. Several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of $T_{i} \simeq 1.7$ keV, except in the few scenarios in which turbulence can be suppressed. In order to understand turbulent mechanisms at play, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. In this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ($k_\perp\rho_i\simeq 1$) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients $a/L_{ne}$ and $a/L_{Ti}$. Furthermore, we discuss the importance of other potentially relevant parameters such as $T_e/T_i$, $E_r$ or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio $\eta_i=L_{ne}/L_{Ti}$, as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which $T_{i,core} > 1.7$ keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.

D. Carralero T. Estrada E. Maragkoudakis T. Windisch J. A. Alonso M. Beurskens S. Bozhenkov I. Calvo H. Damm O. Ford G. Fuchert J. M. García-Regaña N. Pablant E. Sánchez E. Pasch J. L. Velasco the Wendelstein 7-X team

Neural dynamics is often investigated with tools from bifurcation theory. However, many neuron models are stochastic, mimicking fluctuations in the input from unknown parts of the brain or the spiking nature of signals. Noise changes the dynamics with respect to the deterministic model; in particular bifurcation theory cannot be applied. We formulate stochastic neuronal dynamics in the Martin-Siggia-Rose de Dominicis-Janssen (MSRDJ) formalism and present the fluctuation expansion of the effective action and the functional renormalization group (fRG) as two systematic ways to incorporate corrections to the mean dynamics and time-dependent statistics due to fluctuations in the presence of nonlinear neuronal gain. To formulate self-consistency equations, we derive a fundamental link between the effective action in the Onsager-Machlup(OM) formalism, which allows the study of phase transitions, and the MSRDJ effective action, which is computationally advantageous. These results in particular allow the derivation of an OM effective action for systems with non-Gaussian noise. This approach naturally leads to effective deterministic equations for the first moment of the stochastic system; they explain how nonlinearities and noise cooperate to produce memory effects. Moreover, the MSRDJ formulation yields an effective linear system that has identical power spectra and linear response. Starting from the better known loopwise approximation, we then discuss the use of the fRG as a method to obtain self-consistency beyond the mean. We present a new efficient truncation scheme for the hierarchy of flow equations for the vertex functions by adapting the Blaizot, M\'endez and Wschebor approximation from the derivative expansion to the vertex expansion. The methods are presented by means of the simplest possible example of a stochastic differential equation that has generic features of neuronal dynamics.

Jonas Stapmanns Tobias Kühn David Dahmen Thomas Luu Carsten Honerkamp Moritz Helias

The algorithms-with-predictions framework has been used extensively to develop online algorithms with improved beyond-worst-case competitive ratios. Recently, there is growing interest in leveraging predictions for designing data structures with improved beyond-worst-case running times. In this paper, we study the fundamental data structure problem of maintaining approximate shortest paths in incremental graphs in the algorithms-with-predictions model. Given a sequence $\sigma$ of edges that are inserted one at a time, the goal is to maintain approximate shortest paths from the source to each vertex in the graph at each time step. Before any edges arrive, the data structure is given a prediction of the online edge sequence $\hat{\sigma}$ which is used to ``warm start'' its state. As our main result, we design a learned algorithm that maintains $(1+\epsilon)$-approximate single-source shortest paths, which runs in $\tilde{O}(m \eta \log W/\epsilon)$ time, where $W$ is the weight of the heaviest edge and $\eta$ is the prediction error. We show these techniques immediately extend to the all-pairs shortest-path setting as well. Our algorithms are consistent (performing nearly as fast as the offline algorithm) when predictions are nearly perfect, have a smooth degradation in performance with respect to the prediction error and, in the worst case, match the best offline algorithm up to logarithmic factors. As a building block, we study the offline incremental approximate single-source shortest-paths problem. In this problem, the edge sequence $\sigma$ is known a priori and the goal is to efficiently return the length of the shortest paths in the intermediate graph $G_t$ consisting of the first $t$ edges, for all $t$. Note that the offline incremental problem is defined in the worst-case setting (without predictions) and is of independent interest.

Samuel McCauley Benjamin Moseley Aidin Niaparast Helia Niaparast Shikha Singh

Special Issue: New Stellarators New stellarators are being constructed throughout the world as part of a well-coordinated international program. This issue concentrates on the engineering and construction details of Large Helical Device (LHD), Toki, Japan Wendelstein 7-X (W7-X), Greifswald, Germany Helias Stellarator Experiment (HSX), Madison Wisconsin, USA TJ-II, Madrid, Spain

James A. Rome

The elementary low energy excitations in the fractional quantum Hall (FQH) regime are known to be fractionally charged quasiparticles and quasiholes. This work focusses on quasiholes in a finite system of a few electrons treated by exact numerical diagonalization. Recent experimental [Chung et al., PRB 67, 201104 (2003)] and theoretical [Kane, Fisher, PRB 67, 045307 (2003)] work on quantum point contacts in the FQH regime showed strong evidence for single quasiparticle tunneling to be forbidden at low temperature. We want to gain insight into such constricted FQH systems by a numerical approach. Initially the basics of numerically treating a system in rectangular geometry and procedures to insert a quasihole excitation are presented. The excitations' stability and fractional charge is confirmed for different interactions. The question of quasihole tunneling is approached in two different systems: 1. A system posessing bound, localized states of a quasihole shows evidence for tunnel coupling between these states to exist. 2. Different potential forms to model a quantum point contact are investigated and the injection of a single quasihole near the constriction is surveyed in the system's time evolution.

Moritz Helias Daniela Pfannkuche