Measurement and structure of spiral wave response functions.

The rotating spiral waves that emerge in diverse natural and man-made systems typically exhibit a particle-like behaviour since their adjoint critical eigenmodes (response functions) are often seen to be localised around the spiral core. We present a simple method to numerically compute response functions for circular-core and meandering spirals by recording their drift response to many elementary perturbations. Although our method is computationally more expensive than solving the adjoint system, our technique is fully parallellisable, does not suffer from memory limitations and can be applied to experiments. For a cardiac tissue model with the linear spiral core, we find that the response functions are localised near the turning points of the trajectory.

Generalized minimal principle for rotor filaments.

To a reaction-diffusion medium with an inhomogeneous anisotropic diffusion tensor D, we add a fourth spatial dimension such that the determinant of the diffusion tensor is constant in four dimensions. We propose a generalized minimal principle for rotor filaments, stating that the scroll wave filament strives to minimize its surface area in the higher-dimensional space. As a consequence, stationary scroll wave filaments in the original 3D medium are geodesic curves with respect to the metric tensor G=det(D)D(-1). The theory is confirmed by numerical simulations for positive and negative filament tension and a model with a non-stationary spiral core. We conclude that filaments in cardiac tissue with positive tension preferentially reside or anchor in regions where cardiac cells are less interconnected, such as portions of the cardiac wall with a large number of cleavage planes.

Matrix product states for gauge field theories.

The matrix product state formalism is used to simulate Hamiltonian lattice gauge theories. To this end, we define matrix product state manifolds which are manifestly gauge invariant. As an application, we study (1+1)-dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model, and are able to determine very accurately the ground-state properties and elementary one-particle excitations in the continuum limit. In particular, a novel particle excitation in the form of a heavy vector boson is uncovered, compatible with the strong coupling expansion in the continuum. We also study full quantum nonequilibrium dynamics by simulating the real-time evolution of the system induced by a quench in the form of a uniform background electric field.

Entanglement renormalization for quantum fields in real space.

We show how to construct renormalization group (RG) flows of quantum field theories in real space, as opposed to the usual Wilsonian approach in momentum space. This is achieved by generalizing the multiscale entanglement renormalization ansatz to continuum theories. The variational class of wave functions arising from this RG flow are translation invariant and exhibits an entropy-area law. We illustrate the construction for a free nonrelativistic boson model, and argue that the full power of the construction should emerge in the case of interacting theories.

Buckling of scroll waves.

A scroll wave in a sufficiently thin layer of an excitable medium with negative filament tension can be stable nevertheless due to filament rigidity. Above a certain critical thickness of the medium, such a scroll wave will have a tendency to deform into a buckled, precessing state. Experimentally this will be seen as meandering of the spiral wave on the surface, the amplitude of which grows with the thickness of the layer, until a breakup to scroll wave turbulence happens. We present a simplified theory for this phenomenon and illustrate it with numerical examples.

Accurate eikonal-curvature relation for wave fronts in locally anisotropic reaction-diffusion systems.

The dependency of wave velocity in reaction-diffusion (RD) systems on the local front curvature determines not only the stability of wave propagation, but also the fundamental properties of other spatial configurations such as vortices. This Letter gives the first derivation of a covariant eikonal-curvature relation applicable to general RD systems with spatially varying anisotropic diffusion properties, such as cardiac tissue. The theoretical prediction that waves which seem planar can nevertheless possess a nonvanishing geometrical curvature induced by local anisotropy is confirmed by numerical simulations, which reveal deviations up to 20% from the nominal plane wave speed.

Time-dependent variational principle for quantum lattices.

We develop a new algorithm based on the time-dependent variational principle applied to matrix product states to efficiently simulate the real- and imaginary-time dynamics for infinite one-dimensional quantum lattices. This procedure (i) is argued to be optimal, (ii) does not rely on the Trotter decomposition and thus has no Trotter error, (iii) preserves all symmetries and conservation laws, and (iv) has low computational complexity. The algorithm is illustrated by using both an imaginary-time and a real-time example.

Applying the variational principle to (1+1)-dimensional quantum field theories.

We extend the recently introduced continuous matrix product state variational class to the setting of (1+1)-dimensional relativistic quantum field theories. This allows one to overcome the difficulties highlighted by Feynman concerning the application of the variational procedure to relativistic theories, and provides a new way to regularize quantum field theories. A fermionic version of the continuous matrix product state is introduced which is manifestly free of fermion doubling and sign problems. We illustrate the power of the formalism by studying the momentum occupation for free massive Dirac fermions and the chiral symmetry breaking in the Gross-Neveu model.

Variational principle for nonlinear wave propagation in dissipative systems.

The dynamics of many natural systems is dominated by nonlinear waves propagating through the medium. We show that in any extended system that supports nonlinear wave fronts with positive surface tension, the asymptotic wave-front dynamics can be formulated as a gradient system, even when the underlying evolution equations for the field variables cannot be written as a gradient system. The variational potential is simply given by a linear combination of the occupied volume and surface area of the wave front and changes monotonically over time.

Modified ionic models of cardiac tissue for efficient large scale computations.

Recirculation of excitation, or re-entry, is one of the most important mechanisms of life-threatening cardiac arrhythmias and fibrillation. Modelling these phenomena requires large scale computations in two and three-dimensional slabs of cardiac tissue. Because of computational constraints, most of the studies use simplified (non-ionic) models of cardiac tissue, which are electrophysiologically less accurate than the detailed ionic models. In this paper, we propose a method to modify ionic models of cardiac tissue into an intermediate class of models, which are almost as efficient for computations as simplified models, and retain most of the properties of the original ionic models, such as the shape of the action potential, the restitution of action potential duration and of the conduction velocity, as well as unchanged description of most of the ionic currents.

Effective dynamics of twisted and curved scroll waves using virtual filaments.

Scroll waves are three-dimensional excitation patterns that rotate around a central filament curve; they occur in many physical, biological, and chemical systems. We explicitly derive the equations of motion for scroll wave filaments in reaction-diffusion systems with isotropic diffusion up to third order in the filament's twist and curvature. The net drift components define at every instance of time a virtual filament which lies close to the instantaneous filament. Importantly, virtual filaments obey simpler, time-independent laws of motion which we analytically derive here and illustrate with numerical examples. Stability analysis of scroll waves is performed using virtual filaments, showing that filament curvature and twist add as quadratic terms to the nominal filament tension. Applications to oscillating chemical reactions and cardiac tissue are discussed.

Drift laws for spiral waves on curved anisotropic surfaces.

Rotating spiral waves organize spatial patterns in chemical, physical, and biological excitable systems. Factors affecting their dynamics, such as spatiotemporal drift, are of great interest for particular applications. Here, we propose a quantitative description for spiral wave dynamics on curved surfaces which shows that for a wide class of systems, including the Belousov-Zhabotinsky reaction and anisotropic cardiac tissue, the Ricci curvature scalar of the surface is the main determinant of spiral wave drift. The theory provides explicit equations for spiral wave drift direction, drift velocity, and the period of rotation. Depending on the parameters, the drift can be directed to the regions of either maximal or minimal Ricci scalar curvature, which was verified by direct numerical simulations.

Covariant stringlike dynamics of scroll wave filaments in anisotropic cardiac tissue.

It has been hypothesized that stationary scroll wave filaments in cardiac tissue describe a geodesic in a curved space whose metric is the inverse diffusion tensor. Several numerical studies support this hypothesis, but no analytical proof has been provided yet for general anisotropy. In this Letter, we derive dynamic equations for the filament in the case of general anisotropy. These equations are covariant under general spatial coordinate transformations and describe the motion of a stringlike object in a curved space whose metric tensor is the inverse diffusion tensor. Therefore the behavior of scroll wave filaments in excitable media with anisotropy is similar to the one of cosmic strings in a curved universe. Our dynamic equations are valid for thin filaments and for general anisotropy. We show that stationary filaments obey the geodesic equation.

Optimal teleportation with a mixed state of two qubits.

We consider a single copy of a mixed state of two qubits and derive the optimal trace-preserving local operations assisted by classical communication such as to maximize the fidelity of teleportation that can be achieved with this state. These optimal local operations turn out to be implementable by one-way communication and always yield a teleportation fidelity larger than 2/3 if the original state is entangled. This maximal achievable fidelity is an entanglement measure and turns out quantifying the minimal amount of mixing required to destroy the entanglement in a quantum state.