Maze solvers demystified and some other thoughts

There is a growing interest towards implementation of maze solving in spatially-extended physical, chemical and living systems. Several reports of prototypes attracted great publicity, e.g. maze solving with slime mould and epithelial cells, maze navigating droplets. We show that most prototypes utilise one of two phenomena: a shortest path in a maze is a path of the least resistance for fluid and current flow, and a shortest path is a path of the steepest gradient of chemoattractants. We discuss that substrates with so-called maze-solving capabilities simply trace flow currents or chemical diffusion gradients. We illustrate our thoughts with a model of flow and experiments with slime mould. The chapter ends with a discussion of experiments on maze solving with plant roots and leeches which show limitations of the chemical diffusion maze-solving approach.Comment: This is a preliminary version of the chapter to be published in Adamatzky A. (Ed.) Shortest path solvers. From software to wetware. Springer, 201

The study of the functionality of cardiomyocytes obtained from induced pluripotent stem cells for the modeling of cardiac arrhythmias based on long QT syndrome

There are risk factors that lead the normal conduction of excitation in the heart into a chaotic one. These factors include hereditary and acquired channelopathies. Many dangerous changes in the work of the heart can be identified using the patient’s electrocardiogram. Such relatively easily detectable changes include the long QT interval syndrome (LQTS). Despite a relatively high prevalence of hereditary LQTS, to which it is necessary to add both hereditary and induced LQTS as well as the ease of detection on the ECG, the mechanism of reentry formation in this syndrome is still un­known. What should be noted is a high variability of the hereditary syndrome and the fact of the connection between the increase in the heart rate and the risk of cardiac arrest. After an electrophysiological study on individual cardiac cells from patients with the LQT syndrome, it became apparent that the search for a mechanism for the transition of the normal heart rhythm to chaotic and fibrillation cannot be limited to recording ion currents in single cells. To solve this problem, we need a model of the behavior of cardiac tissue which reflects the relationship of various factors and the risk of reentry. In order to create an experimental model of LQTS in our work, the iPSC of a pati­ent-specific line from a healthy patient was differentiated into a monolayer of cardiac cells and the parameters of the excitation propagation were studied depending on the stage of differentiation. It was shown that a stable value of the propagation velocity and the response to periodic stimulation in the range of physiological values, are reached after the 30th day of dif­ferentiation

Nuclear spin driven quantum relaxation in LiY_0.998Ho_0.002F_4

Staircase hysteresis loops of the magnetization of a LiY_0.998Ho_0.002F_4 single crystal are observed at subkelvin temperatures and low field sweep rates. This behavior results from quantum dynamics at avoided level crossings of the energy spectrum of single Ho^{3+} ions in the presence of hyperfine interactions. Enhanced quantum relaxation in constant transverse fields allows the study of the relative magnitude of tunnel splittings. At faster sweep rates, non-equilibrated spin-phonon and spin-spin transitions, mediated by weak dipolar interactions, lead to magnetization oscillations and additional steps.Comment: 5 pages, 5 eps figures, using RevTe

Theory of spiral wave dynamics in weakly excitable media: asymptotic reduction to a kinematic model and applications

In a weakly excitable medium, characterized by a large threshold stimulus, the free end of an isolated broken plane wave (wave tip) can either rotate (steadily or unsteadily) around a large excitable core, thereby producing a spiral pattern, or retract causing the wave to vanish at boundaries. An asymptotic analysis of spiral motion and retraction is carried out in this weakly excitable large core regime starting from the free-boundary limit of the reaction-diffusion models, valid when the excited region is delimited by a thin interface. The wave description is shown to naturally split between the tip region and a far region that are smoothly matched on an intermediate scale. This separation allows us to rigorously derive an equation of motion for the wave tip, with the large scale motion of the spiral wavefront slaved to the tip. This kinematic description provides both a physical picture and exact predictions for a wide range of wave behavior, including: (i) steady rotation (frequency and core radius), (ii) exact treatment of the meandering instability in the free-boundary limit with the prediction that the frequency of unstable motion is half the primary steady frequency (iii) drift under external actions (external field with application to axisymmetric scroll ring motion in three-dimensions, and spatial or/and time-dependent variation of excitability), and (iv) the dynamics of multi-armed spiral waves with the new prediction that steadily rotating waves with two or more arms are linearly unstable. Numerical simulations of FitzHug-Nagumo kinetics are used to test several aspects of our results. In addition, we discuss the semi-quantitative extension of this theory to finite cores and pinpoint mathematical subtleties related to the thin interface limit of singly diffusive reaction-diffusion models

Scroll waves in isotropic excitable media : linear instabilities, bifurcations and restabilized states

Scroll waves are three-dimensional analogs of spiral waves. The linear stability spectrum of untwisted and twisted scroll waves is computed for a two-variable reaction-diffusion model of an excitable medium. Different bands of modes are seen to be unstable in different regions of parameter space. The corresponding bifurcations and bifurcated states are characterized by performing direct numerical simulations. In addition, computations of the adjoint linear stability operator eigenmodes are also performed and serve to obtain a number of matrix elements characterizing the long-wavelength deformations of scroll waves.Comment: 30 pages 16 figures, submitted to Phys. Rev.

Sub-millimeter to centimeter excess emission from the Magellanic Clouds. I. Global spectral energy distribution

In order to reconstruct the global SEDs of the Magellanic Clouds over eight decades in spectral range, we combined literature flux densities representing the entire LMC and SMC respectively, and complemented these with maps extracted from the WMAP and COBE databases covering the missing the 23--90 GHz (13--3.2 mm) and the poorly sampled 1.25--250 THz (240--1.25 micron). We have discovered a pronounced excess of emission from both Magellanic Clouds, but especially the SMC, at millimeter and sub-millimeter wavelengths. We also determined accurate thermal radio fluxes and very low global extinctions for both LMC and SMC. Possible explanations are briefly considered but as long as the nature of the excess emission is unknown, the total dust masses and gas-to-dust ratios of the Magellanic Clouds cannot reliably be determined.Comment: Accepted for publication by A&

Order Parameter Equations for Front Transitions: Planar and Circular Fronts

Near a parity breaking front bifurcation, small perturbations may reverse the propagation direction of fronts. Often this results in nonsteady asymptotic motion such as breathing and domain breakup. Exploiting the time scale differences of an activator-inhibitor model and the proximity to the front bifurcation, we derive equations of motion for planar and circular fronts. The equations involve a translational degree of freedom and an order parameter describing transitions between left and right propagating fronts. Perturbations, such as a space dependent advective field or uniform curvature (axisymmetric spots), couple these two degrees of freedom. In both cases this leads to a transition from stationary to oscillating fronts as the parity breaking bifurcation is approached. For axisymmetric spots, two additional dynamic behaviors are found: rebound and collapse.Comment: 9 pages. Aric Hagberg: http://t7.lanl.gov/People/Aric/; Ehud Meron: http://www.bgu.ac.il/BIDR/research/staff/meron.htm

Reappraising foreground contamination in the COBE-DMR data

With the advent of all-sky H-Alpha surveys it is possible to determine a reliable free-free template of the diffuse interstellar medium (Dickinson, Davies & Davis 2003) which can be used in conjunction with the synchrotron and dust templates to correct CMB observations for diffuse Galactic foregrounds. We have used the COBE-DMR data at 31.5, 53 and 90 GHz and employed cross- correlation techniques to re-evaluate the foreground contributions, particularly that due to dust which is known to be correlated with H-Alpha (and free-free) emission.The DMR microwave maps are found to contain, as well as the expected synchrotron and free-free components, a component tightly correlated to the 140 micron dust maps of DIRBE. At 31.5, 53 and 90 GHz this emission is 6.3 +/- 0.6, 2.4 +/- 0.4 and 2.2 +/- 0.4 microK/(MJy/sr) at 140 microns respectively. When corrected for the contribution from thermal dust following model 7 of Finkbeiner, Davis & Schlegel (1999), a strong anomalous dust emission component remains, which is well-fitted by a frequency spectrum of the form $\nu^{-\beta}$ where beta ~ 2.5 in the DMR frequency range; this is the dominant foreground at 31.5 GHz. The result implies the presence of an emission component with a dust-like morphology but a synchrotron-like spectrum. We discuss the possible origins of this component and compare it with the recent WMAP interpretation. The better knowledge of the individual foregrounds provided by the present study enables a larger area of the sky (|b|>15) to be used to re-appraise the CMB quadrupole normalisation, Qrms-PS, and the scalar perturbations spectral index, n. We find Qrms-PS=15.2^{+2.8}_{-2.3} with a power-law spectral index of n=1.2 +/- 0.2. These values are consistent with previous COBE-DMR analyses and the WMAP 1-year analysis.Comment: 17 pages, 6 figures, replaced with version accepted for publication in MNRAS, contains new section comparing results to the WMAP foreground studie

A submillimeter exponential disk in M~51: evidence for an extended cold dust disk

A 850 micron map of the interacting spiral galaxy M51 shows well-defined spiral arms, closely resembling the structures seen in CO and HI emission. However, most of the 850 micron emission originates in an underlying exponential disk, a component that has not been observed before in a face-on galaxy at these wavelengths. The scale-length of this disk is 5.45 kpc, which is somewhat larger than the scale-length of the stellar disk, but somewhat smaller than that of atomic hydrogen. Its profile can not be explained solely by a radial disk temperature gradient but requires the underlying dust to have an exponential distribution as well. This reinforces the view that the submm emission from spiral galaxy disks traces total hydrogen column density, i.e.the sum of H2 and HI. A canonical gas-to-dust ratio of 100+/-26 is obtained for kappa(850)=1.2 g**-1 cm**2, where kappa(850) is the dust opacity at 850 micron.Comment: accepted by A&A, 8 page

Fitting the integrated Spectral Energy Distributions of Galaxies

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details ofdust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (sedfitting.org), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.Comment: 54 pages, 26 figures, Accepted for publication in Astrophysics & Space Scienc