title: A Polynomial Delay Algorithm for Enumerating Minimal Dominating Sets in Chordal Graphs

abstract: An output-polynomial algorithm for the listing of minimal dominating sets in graphs is a challenging open problem and is known to be equivalent to the well-known Transversal problem which asks for an output-polynomial algorithm for listing the set of minimal hitting sets in hypergraphs. We give a polynomial delay algorithm to list the set of minimal dominating sets in chordal graphs, an important and well-studied graph class where such an algorithm was open for a while.


title: Which of the world's institutions employ the most highly cited researchers? An analysis of the data from highlycited.com

abstract: A few weeks ago, Thomson Reuters published a list of the highly cited researchers worldwide (highlycited.com). Since the data is freely available for downloading and includes the names of the researchers' institutions, we produced a ranking of the institutions on the basis of the number of highly cited researchers per institution. This ranking is intended to be a helpful amendment of other available institutional rankings.


title: Non thermal and purely electronic resistive transition in narrow gap Mott insulators

abstract: Mott insulator to metal transitions under electric field are currently the subject of numerous fundamental and applied studies. This puzzling effect, which involves non-trivial out-of-equilibrium effects in correlated systems, is indeed at play in the operation of a new class of electronic memories, the Mott memories. However the combined electronic and thermal effects are difficult to disentangle in Mott insulators undergoing such transitions. We report here a comparison between the properties under electric field of a canonical Mott insulator and a model built on a realistic 2D resistor network able to capture both thermal effects and electronic transitions. This comparison made specifically on the family of narrow gap Mott insulators AM4Q8, (A = Ga or Ge; M=V, Nb or Ta, and Q = S or Se) unambiguously establishes that the resistive transition experimentally observed under electric field arises from a purely electronic mechanism.


title: Amplitude modulated phase in a Bose-Einstein condensate: the role of non-local interactions

abstract: We consider a Gross-Pitaevskii model of BEC with non-local interactions of range of the order of the s-wave scattering length. With this model, we study the density modulated phase in 1D and 2D, which are solutions of this modified model along with the usual uniform density state. We find an exact free energy functional for our model and show that the 1D density modulated state can have lower energy than the uniform density state. Although, the density modulated state can be made to be energetically favourable, we show also that, this state is inherently dynamically unstable due to the coupling of instabilities to the spatial order.


title: The density matrix renormalization group for ab initio quantum chemistry

abstract: During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. Its underlying wavefunction ansatz, the matrix product state (MPS), is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS, the rank of the decomposition, controls the size of the corner of the many-body Hilbert space that can be reached with the ansatz. This parameter can be systematically increased until numerical convergence is reached. The MPS ansatz naturally captures exponentially decaying correlation functions. Therefore DMRG works extremely well for noncritical one-dimensional systems. The active orbital spaces in quantum chemistry are however often far from one-dimensional, and relatively large virtual dimensions are required to use DMRG for ab initio quantum chemistry (QC-DMRG). The QC-DMRG algorithm, its computational cost, and its properties are discussed. Two important aspects to reduce the computational cost are given special attention: the orbital choice and ordering, and the exploitation of the symmetry group of the Hamiltonian. With these considerations, the QC-DMRG algorithm allows to find numerically exact solutions in active spaces of up to 40 electrons in 40 orbitals.


title: ImpNet: Programming Software-Defied Networks Using Imperative Techniques

abstract: Software and hardware components are basic parts of modern networks. However the software compo- nent is typical sealed and function-oriented. Therefore it is very difficult to modify these components. This badly affected networking innovations. Moreover, this resulted in network policies having complex interfaces that are not user-friendly and hence resulted in huge and complicated flow tables on physical switches of networks. This greatly degrades the network performance in many cases. Software-Defined Networks (SDNs) is a modern architecture of networks to overcome issues mentioned above. The idea of SDN is to add to the network a controller device that manages all the other devices on the network including physical switches of the network. One of the main tasks of the managing process is switch learning; achieved via programming physical switches of the network by adding or removing rules for packet-processing to/from switches, more specifically to/from their flow tables. A high-level imperative network programming language, called ImpNet, is presented in this paper. ImpNet enables writing efficient, yet simple, and powerful programs to run on the controller to control all other network devices including switches. ImpNet is compositional, simply-structured, expressive, and more importantly imperative. The syntax of ImpNet together two types of operational semantics to contracts of ImpNet are presented in the paper. The proposed semantics are of the static and dynamic types. Two modern application programmed using ImpNet are shown in the paper as well. The semantics of the applications are shown in the paper also.


title: Time integration of tensor trains

abstract: A robust and efficient time integrator for dynamical tensor approximation in the tensor train or matrix product state format is presented. The method is based on splitting the projector onto the tangent space of the tensor manifold. The algorithm can be used for updating time-dependent tensors in the given data-sparse tensor train / matrix product state format and for computing an approximate solution to high-dimensional tensor differential equations within this data-sparse format. The formulation, implementation and theoretical properties of the proposed integrator are studied, and numerical experiments with problems from quantum molecular dynamics and with iterative processes in the tensor train format are included.


title: Off-diagonal matrix elements of local operators in many-body quantum systems

abstract: In the time evolution of isolated quantum systems out of equilibrium, local observables generally relax to a long-time asymptotic value, governed by the expectation values (diagonal matrix elements) of the corresponding operator in the eigenstates of the system. The temporal fluctuations around this value, response to further perturbations, and the relaxation toward this asymptotic value, are all determined by the off-diagonal matrix elements. Motivated by this non-equilibrium role, we present generic statistical properties of off-diagonal matrix elements of local observables in two families of interacting many-body systems with local interactions. Since integrability (or lack thereof) is an important ingredient in the relaxation process, we analyze models that can be continuously tuned to integrability. We show that, for generic non-integrable systems, the distribution of off-diagonal matrix elements is a gaussian centered at zero. As one approaches integrability, the peak around zero becomes sharper, so that the distribution is approximately a combination of two gaussians. We characterize the proximity to integrability through the deviation of this distribution from a gaussian shape. We also determine the scaling dependence on system size of the average magnitude of off-diagonal matrix elements.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el quant-ph


title: Tracking Individual Targets in High Density Crowd Scenes Analysis of a Video Recording in Hajj 2009

abstract: In this paper we present a number of methods (manual, semi-automatic and automatic) for tracking individual targets in high density crowd scenes where thousand of people are gathered. The necessary data about the motion of individuals and a lot of other physical information can be extracted from consecutive image sequences in different ways, including optical flow and block motion estimation. One of the famous methods for tracking moving objects is the block matching method. This way to estimate subject motion requires the specification of a comparison window which determines the scale of the estimate. In this work we present a real-time method for pedestrian recognition and tracking in sequences of high resolution images obtained by a stationary (high definition) camera located in different places on the Haram mosque in Mecca. The objective is to estimate pedestrian velocities as a function of the local density.The resulting data of tracking moving pedestrians based on video sequences are presented in the following section. Through the evaluated system the spatio-temporal coordinates of each pedestrian during the Tawaf ritual are established. The pilgrim velocities as function of the local densities in the Mataf area (Haram Mosque Mecca) are illustrated and very precisely documented.


title: Cassini ISS Mutual Event Astrometry of the Mid-sized Saturnian Satellites 2005-2012

abstract: We present astrometric observations of the Saturnian satellites Mimas, Enceladus, Tethys, Dione and Rhea from Cassini Imaging Science Subsystem (ISS) narrow-angle camera (NAC) images. Image sequences were designed to observe mutual occultations between these satellites. The positions of satellite centres were estimated by fitting ellipsoidal shape models to the measured limbs of the imaged satellites. Spacecraft pointing corrections were computed using the UCAC2 star catalogue. We provide a total of 2303 astrometric observations, resulting in 976 pairs, the remainder consisting of observations of a single satellite. Mean residuals for the individual satellite positions relative to the SAT360 ephemeris were 4.3 km in the line direction and -2.4 km in the sample direction, with standard deviations of 5.6 and 7.0 km respectively, an order of magnitude improvement in precision compared to published HST observations. By considering inter-satellite separations, uncertainties in camera pointing and spacecraft positioning along with possible biases in the individual positions of the satellites can be largely eliminated, resulting in an order-of-magnitude increase in accuracy compared to that achievable using the individual satellite positions. We show how factors relating to the viewing geometry cause small biases in the individual positions of order 0.28 pixel to become systematic across the dataset as a whole and discuss options for reducing their effects . The reduced astrometric data are provided in the form of individual positions for each satellite, together with the measured positions of reference stars, in order to allow more flexibility in the processing of the observations, taking into account possible future advances in limb-fitting techniques as well as the future availability of more accurate star catalogues, such as those from the GAIA mission.


title: Comparison between thermal and current driven spin-transfer torque in nanopillar metallic spin valves

abstract: We investigate the relation between thermal spin-transfer torque (TSTT) and the spin-dependent Seebeck effect (SDSE), which produces a spin current when a temperature gradient is applied across a metallic ferromagnet, in nanopillar metallic spin valves. Comparing its angular dependence (aSDSE) with the angle dependent magnetoresistance (aMR) measurements on the same device, we are able to verify that a small spin heat accumulation builds up in our devices. From the SDSE measurement and the observed current driven STT switching current of 0.8 mA in our spin valve devices, it was estimated that a temperature difference of 230 K is needed to produce an equal amount of TSTT. Experiments specifically focused on investigating TSTT show a response that is dominated by overall heating of the magnetic layer. Comparing it to the current driven STT experiments we estimate that only ~10% of the response is due to TSTT. This leads us to conclude that switching dominated by TSTT requires a direct coupling to a perfect heat sink to minimize the effect of overall heating. Nevertheless the combined effect of heating, STT and TSTT could prove useful for inducing magnetization switching when further investigated and optimized.


title: Thermal, chemical and spectral equilibration in heavy-ion collisions

abstract: We have considered the equilibration in a relativistic heavy ion collision using our transport model. We applied periodic boundary conditions to close the system in a box. We found that the thermal equilibration takes place in the first 20-40 fm/c which time is comparable to the duration of a heavy ion collision. The chemical equilibration is a much slower process and the system does not equilibrate in a heavy ion collision. We have also shown that the mass spectra of broad resonances immediately follows their in-medium spectral functions.


title: The supergiant shell with triggered star formation in Irr galaxy IC 2574: neutral and ionized gas kinematics

abstract: We analyse the ionized gas kinematics in the star formation regions of the supergiant shell (SGS) of the IC 2574 galaxy using observations with the Fabry-Perot interferometer at the 6-m telescope of SAO RAS; the data of the THINGS survey are used to analyze the neutral gas kinematics in the area. We perform the 'derotation' of the H-alpha and HI data cubes and show its efficiency in kinematics analysis. We confirm the SGS expansion velocity 25 km/s obtained by Walter & Brinks (1999) and conclude that the SGS is located at the far side of the galactic disc plane. We determine the expansion velocities, kinematic ages, and the required mechanical energy input rates for four star formation complexes in the walls of the SGS; for the remaining ones we give the limiting values of the above parameters. A comparison with the age and energy input of the complexes' stellar population shows that sufficient energy is fed to all HII regions except one. We discuss in detail the possible nature of this region and that of another one, which was believed to be an SNR according to radio observations. We measured the expansion velocity of the latter and confirm its identification as an old SNR. Our observations allowed us to identify a faint diffuse H-alpha emission inside the SGS which was never observed before.


title: Effect of divalent Ba cation substitution with Sr on coupled multiglass state in the magnetoelectric multiferroic compound Ba3NbFe3Si2O14

abstract: (Ba/Sr)3NbFe3Si2O14 is a magnetoelectric multiferroic with an incommensurate antiferromagnetic spiral magnetic structure which induces electric polarization at 26 K. The structure, as revealed by x-rays and neutrons, as well as static and dynamic magnetic and dielectric properties of these compounds have been studied down to 6 K under different conditions. Both the compounds have similar crystal structure but with different lattice constants down to 6 K. The Ba-and Sr-compounds exhibit a transition at 26 K and 25 K respectively, as indicated by the specific heat capacity and dc specific magnetization, into an antiferromagnetic state. Although Ba and Sr are isovalent, they exhibit very different static and dynamic magnetic behavior. The Ba-compound exhibits both thermal and magnetic field hysteresis with the thermal hysteresis decreasing with increasing magnetic field, a behavior typical of glasses. The glassy behavior is also clearly seen in the ac susceptibility studies which show a dispersive peak in the range 40 K to 90 K in the frequency range 10^1 Hz to 10^4 Hz. The dispersive behavior follows a cluster glass critical slowing dynamics with a freezing temperature of 35 K and a critical exponent of 3.9, a value close to the 3-D Ising model. The Sr-compound however does not exhibit any dispersive behavior except for the invariant transition at 25 K in ac susceptibility with no magnetic field hysteresis at all temperatures. The dielectric constant studied in the frequency range 10^1 Hz to 10^6 Hz also reflects the magnetic behavior of the two compounds. The Ba-compound has two distinct dispersive peaks near TN and in the range 40 K to 125 K while the Sr-compound has a single dispersive peak in the range 40 K to 80 K. The activation energy of the high temperature dispersive peak in both compounds however is found to be similar, 71 meV and 65 meV respectively for Ba- and Sr-compounds.


title: On the particle-hole symmetry of the fermionic spinless Hubbard model in $D=1$

abstract: We revisit the particle-hole symmetry of the one-dimensional ($D=1$) fermionic spinless Hubbard model, associating that symmetry to the invariance of the Helmholtz free energy of the one-dimensional spin-1/2 $XXZ$ Heisenberg model, under reversal of the longitudinal magnetic field and at any finite temperature. Upon comparing two regimes of that chain model so that the number of particles in one regime equals the number of holes in the other, one finds that, in general, their thermodynamics is similar, but not identical: both models share the specific heat and entropy functions, but not the internal energy per site, the first-neighbor correlation functions, and the number of particles per site. Due to that symmetry, the difference between the first-neighbor correlation functions is proportional to the $z$-component of magnetization of the $XXZ$ Heisenberg model. The results presented in this paper are valid for any value of the interaction strength parameter $V$, which describes the attractive/null/repulsive interaction of neighboring fermions.


title: Completely positive maps within the framework of direct-sum decomposition of state space

abstract: We investigate completely positive maps for an open system interacting with its environment. The families of the initial states for which the reduced dynamics can be described by a completely positive map are identified within the framework of direct-sum decomposition of state space. They includes not only separable states with vanishing or nonvanishing quantum discord but also entangled states. A general expression of the families as well as the Kraus operators for the completely positive maps are explicitly given. It significantly extends the previous results.


title: MPI-AMRVAC for Solar and Astrophysics

abstract: In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure.


title: On the Enumeration of Minimal Dominating Sets and Related Notions

abstract: A dominating set $D$ in a graph is a subset of its vertex set such that each vertex is either in $D$ or has a neighbour in $D$. In this paper, we are interested in the enumeration of (inclusion-wise) minimal dominating sets in graphs, called the Dom-Enum problem. It is well known that this problem can be polynomially reduced to the Trans-Enum problem in hypergraphs, i.e., the problem of enumerating all minimal transversals in a hypergraph. Firstly we show that the Trans-Enum problem can be polynomially reduced to the Dom-Enum problem. As a consequence there exists an output-polynomial time algorithm for the Trans-Enum problem if and only if there exists one for the Dom-Enum problem. Secondly, we study the Dom-Enum problem in some graph classes. We give an output-polynomial time algorithm for the Dom-Enum problem in split graphs, and introduce the completion of a graph to obtain an output-polynomial time algorithm for the Dom-Enum problem in $P_6$-free chordal graphs, a proper superclass of split graphs. Finally, we investigate the complexity of the enumeration of (inclusion-wise) minimal connected dominating sets and minimal total dominating sets of graphs. We show that there exists an output-polynomial time algorithm for the Dom-Enum problem (or equivalently Trans-Enum problem) if and only if there exists one for the following enumeration problems: minimal total dominating sets, minimal total dominating sets in split graphs, minimal connected dominating sets in split graphs, minimal dominating sets in co-bipartite graphs.


title: Electronic inhomogeneities in the superconducting phase of CaFe1.96Ni0.04As2 single crystals

abstract: Superconductivity in Iron-Arsenic based pnictides emerges in close proximity to an antiferromagnetic (AFM)ordered parent state and the AFM phase overlaps with superconducting (SC) phase in some pnictides for certain range of doping. CaFe2-xNixAs2 belongs to this category, where both the phases overlap. Here we use scanning tunneling microscopy and spectroscopy to investigate the local electronic properties of uderdoped CaFe1.96Ni0.04As2 single crystals in the vicinity of the boundary of the two phases. Both resistivity and magnetic measurements show that a tiny portion (~ 1.2 %) of this compound becomes superconductor below the SC onset temperature TC ~ 15 K. Topographic images show reasonably flat surface with signatures of atomic resolution. High temperature spectra are spatially homogeneous and show signatures of spin density wave (SDW) gap with a finite density of states near the Fermi energy. Below TC, spectra show significant spatial inhomogeneity with a SDW gap everywhere but at some locations we also see an asymmetric or symmetric depression in ~ +/- 5 meV energy range together with the SDW gap. Inhomogeneity reduces significantly as the temperature goes above TC and disappears completely far above TC. These observations are discussed in terms of an inhomogeneous electronic phase that may exist due to the vicinity of this composition to the SC dome boundary on the underdoped side of the phase diagram.


title: Molecular gas heating in Arp 299

abstract: Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest molecules in the Herschel wavelength range. CO transitions are observed with Herschel, and together, these lines trace the excitation of CO. We study Arp 299, a colliding galaxy group, with one component harboring an AGN and two more undergoing intense star formation. For Arp 299 A, we present PACS spectrometer observations of high-J CO lines up to J=20-19 and JCMT observations of $^{13}$CO and HCN to discern between UV heating and alternative heating mechanisms. There is an immediately noticeable difference in the spectra of Arp 299 A and Arp 299 B+C, with source A having brighter high-J CO transitions. This is reflected in their respective spectral energy line distributions. We find that photon-dominated regions (PDRs) are unlikely to heat all the gas since a very extreme PDR is necessary to fit the high-J CO lines. In addition, this extreme PDR does not fit the HCN observations, and the dust spectral energy distribution shows that there is not enough hot dust to match the amount expected from such an extreme PDR. Therefore, we determine that the high-J CO and HCN transitions are heated by an additional mechanism, namely cosmic ray heating, mechanical heating, or X-ray heating. We find that mechanical heating, in combination with UV heating, is the only mechanism that fits all molecular transitions. We also constrain the molecular gas mass of Arp 299 A to 3e9 Msun and find that we need 4% of the total heating to be mechanical heating, with the rest UV heating.


title: Separating Para and Ortho Water

abstract: Water exists as two nuclear-spin isomers, para and ortho, determined by the overall spin of its two hydrogen nuclei. For isolated water molecules the conversion between these isomers is forbidden and they act as different molecular species. Yet, these species are not readily separable, and little is known about their specific physical and chemical properties, conversion mechanisms, or interactions. Here we demonstrate the production of isolated samples of both spin isomers in pure beams of para and ortho water, with both species in their respective absolute ground state. These single-quantum-state samples are ideal targets for unraveling spin-conversion mechanisms, for precision spectroscopy and fundamental-symmetry-breaking studies, and for spin-enhanced applications, e. g., laboratory astrophysics and -chemistry or hypersensitized NMR experiments.


title: Duality between the trigonometric BC(n) Sutherland system and a completed rational Ruijsenaars-Schneider-van Diejen system

abstract: We present a new case of duality between integrable many-body systems, where two systems live on the action-angle phase spaces of each other in such a way that the action variables of each system serve as the particle positions of the other one. Our investigation utilizes an idea that was exploited previously to provide group-theoretic interpretation for several dualities discovered originally by Ruijsenaars. In the group-theoretic framework one applies Hamiltonian reduction to two Abelian Poisson algebras of invariants on a higher dimensional phase space and identifies their reductions as action and position variables of two integrable systems living on two different models of the single reduced phase space. Taking the cotangent bundle of U(2n) as the upstairs space, we demonstrate how this mechanism leads to a new dual pair involving the BC(n) trigonometric Sutherland system. Thereby we generalize earlier results pertaining to the A(n) trigonometric Sutherland system as well as a recent work by Pusztai on the hyperbolic BC(n) Sutherland system.


title: Large dynamics of Yang--Mills theory: mean dimension formula

abstract: This paper studies the Yang--Mills ASD equation over the cylinder as a non-linear evolution equation. We consider a dynamical system consisting of bounded orbits of this evolution equation. This system contains many chaotic orbits, and moreover it becomes an infinite dimensional and infinite entropy system. We study the mean dimension of this huge dynamical system. Mean dimension is a topological invariant of dynamical systems introduced by Gromov. We prove the exact formula of the mean dimension by developing a new technique based on the metric mean dimension theory of Lindenstrauss--Weiss.


title: High frequency homogenisation for elastic lattices

abstract: A complete methodology, based on a two-scale asymptotic approach, that enables the homogenisation of elastic lattices at non-zero frequencies is developed. Elastic lattices are distinguished from scalar lattices in that two or more types of coupled waves exist, even at low frequencies. Such a theory enables the determination of effective material properties at both low and high frequencies. The theoretical framework is developed for the propagation of waves through lattices of arbitrary geometry and dimension. The asymptotic approach provides a method through which the dispersive properties of lattices at frequencies near standing waves can be described; the theory accurately describes both the dispersion curves and the response of the lattice near the edges of the Brillouin zone. The leading order solution is expressed as a product between the standing wave solution and long-scale envelope functions that are eigensolutions of the homogenised partial differential equation. The general theory is supplemented by a pair of illustrative examples for two archetypal classes of two-dimensional elastic lattices. The efficiency of the asymptotic approach in accurately describing several interesting phenomena is demonstrated, including dynamic anisotropy and Dirac cones.


title: On the superconductivity of graphite interfaces

abstract: We propose an explanation for the appearance of superconductivity at the interfaces of graphite with Bernal stacking order. A network of line defects with flat bands appears at the interfaces between two slightly twisted graphite structures. Due to the flat band the probability to find high temperature superconductivity at these quasi one-dimensional corridors is strongly enhanced. When the network of superconducting lines is dense it becomes effectively two-dimensional. The model provides an explanation for several reports on the observation of superconductivity up to room temperature in different oriented graphite samples, graphite powders as well as graphite-composite samples published in the past.


title: Luttinger-field approach to thermoelectric transport in nanoscale conductors

abstract: Thermoelectric transport in nanoscale conductors is analyzed in terms of the response of the system to a thermo-mechanical field, first introduced by Luttinger, which couples to the electronic energy density. While in this approach the temperature remains spatially uniform, we show that a spatially varying thermo-mechanical field effectively simulates a temperature gradient across the system and allows us to calculate the electric and thermal currents that flow due to the thermo-mechanical field. In particular, we show that, in the long-time limit, the currents thus calculated reduce to those that one obtains from the Landauer-B\"uttiker formula, suitably generalized to allow for different temperatures in the reservoirs, if the thermo-mechanical field is applied to prepare the system, and subsequently turned off at ${t=0}$. Alternatively, we can drive the system out of equilibrium by switching the thermo-mechanical field after the initial preparation. We compare these two scenarios, employing a model noninteracting Hamiltonian, in the linear regime, in which they coincide, and in the nonlinear regime in which they show marked differences. We also show how an operationally defined local effective temperature can be computed within this formalism.


title: Surface bundles over surfaces with arbitrarily many fiberings

abstract: In this paper we give the first example of a surface bundle over a surface with at least three fiberings. In fact, for each $n \ge 3$ we construct $4$-manifolds $E$ admitting at least $n$ distinct fiberings $p_i: E \to \Sigma_{g_i}$ as a surface bundle over a surface with base and fiber both closed surfaces of negative Euler characteristic. We give examples of surface bundles admitting multiple fiberings for which the monodromy representation has image in the Torelli group, showing the necessity of all of the assumptions made in the main theorem of our recent paper [arXiv:1404.0066]. Our examples show that the number of surface bundle structures that can be realized on a $4$-manifold $E$ with Euler characteristic $d$ grows exponentially with $d$.


title: Approximation and Streaming Algorithms for Projective Clustering via Random Projections

abstract: Let $P$ be a set of $n$ points in $\mathbb{R}^d$. In the projective clustering problem, given $k, q$ and norm $\rho \in [1,\infty]$, we have to compute a set $\mathcal{F}$ of $k$ $q$-dimensional flats such that $(\sum_{p\in P}d(p, \mathcal{F})^\rho)^{1/\rho}$ is minimized; here $d(p, \mathcal{F})$ represents the (Euclidean) distance of $p$ to the closest flat in $\mathcal{F}$. We let $f_k^q(P,\rho)$ denote the minimal value and interpret $f_k^q(P,\infty)$ to be $\max_{r\in P}d(r, \mathcal{F})$. When $\rho=1,2$ and $\infty$ and $q=0$, the problem corresponds to the $k$-median, $k$-mean and the $k$-center clustering problems respectively. For every $0 < \epsilon < 1$, $S\subset P$ and $\rho \ge 1$, we show that the orthogonal projection of $P$ onto a randomly chosen flat of dimension $O(((q+1)^2\log(1/\epsilon)/\epsilon^3) \log n)$ will $\epsilon$-approximate $f_1^q(S,\rho)$. This result combines the concepts of geometric coresets and subspace embeddings based on the Johnson-Lindenstrauss Lemma. As a consequence, an orthogonal projection of $P$ to an $O(((q+1)^2 \log ((q+1)/\epsilon)/\epsilon^3) \log n)$ dimensional randomly chosen subspace $\epsilon$-approximates projective clusterings for every $k$ and $\rho$ simultaneously. Note that the dimension of this subspace is independent of the number of clusters~$k$. Using this dimension reduction result, we obtain new approximation and streaming algorithms for projective clustering problems. For example, given a stream of $n$ points, we show how to compute an $\epsilon$-approximate projective clustering for every $k$ and $\rho$ simultaneously using only $O((n+d)((q+1)^2\log ((q+1)/\epsilon))/\epsilon^3 \log n)$ space. Compared to standard streaming algorithms with $\Omega(kd)$ space requirement, our approach is a significant improvement when the number of input points and their dimensions are of the same order of magnitude.


title: Robustness against Disorder of Relativistic Spectral Properties in Chalcogenide Alloys

abstract: In order to carefully address the interplay between substitutional disorder and spin-orbit-coupling in IV-VI alloys, we propose a novel theoretical approach that integrates the reliability of plane-wave based density-functional theory beyond the local-density approximation with the Coherent Potential Approximation. By applying the proposed method to ternary chalcogenide alloys, we predict a substantial robustness of spectral features close to the Fermi energy against substitutional disorder. Supplementing our first-principles calculations with the analysis of the $k \cdot p$ model for rock-salt chalcogenides, we show that the disorder self-energy is vanishingly small close to the band gap, thus allowing for bulk Rashba-like spin splitting to be observed in ferroelectric alloys, such as PbS$_x$Te$_{1-x}$, and protecting the band-character inversion related to the topological transition in the recently discovered Topological Crystalline Insulator Pb$_{1-x}$Sn$_x$Te.


title: Fluctuating Currents in Stochastic Thermodynamics I. Gauge Invariance of Asymptotic Statistics

abstract: Stochastic Thermodynamics uses Markovian jump processes to model random transitions between observable mesoscopic states. Physical currents are obtained from anti-symmetric jump observables defined on the edges of the graph representing the network of states. The asymptotic statistics of such currents are characterized by scaled cumulants. In the present work, we use the algebraic and topological structure of Markovian models to prove a gauge invariance of the scaled cumulant-generating function. Exploiting this invariance yields an efficient algorithm for practical calculations of asymptotic averages and correlation integrals. We discuss how our approach generalizes the Schnakenberg decomposition of the average entropy-production rate, and how it unifies previous work. The application of our results to concrete models is presented in an accompanying publication.


title: Impact of occultations of stellar active regions on transmission spectra: Can occultation of a plage mimic the signature of a blue sky?

abstract: Transmission spectroscopy during planetary transits, which is based on the measurements of the variations of planet-to-star radius ratio as a function of wavelength, is a powerful technique to study the atmospheric properties of transiting planets. One of the main limitation of this technique is the effects of stellar activity, which up until now, have been taken into account only by assessing the effect of non-occulted stellar spots on the estimates of planet-to-star radius ratio. In this paper, we study, for the first time, the impact of the occultation of a stellar spot and plage on the transmission spectra of transiting exoplanets. We simulated this effect by generating a large number of transit light curves for different transiting planets, stellar spectral types, and for different wavelengths. Results of our simulations indicate that the anomalies inside the transit light curve can lead to a significant underestimation or overestimation of the planet-to-star radius ratio as a function of wavelength. At short wavelengths, the effect can reach to a difference of up to 10% in the planet-to-star radius ratio, mimicking the signature of light scattering in the planetary atmosphere. Atmospheric scattering has been proposed to interpret the increasing slopes of transmission spectra toward blue for exoplanets HD 189733b and GJ 3470b. Here we show that these signatures can be alternatively interpreted by the occultation of stellar plages. Results also suggest that the best strategy to identify and quantify the effects of stellar activities on the transmission spectrum of a planet is to perform several observations during the transit epoch at the same wavelength. This will allow for identifying the possible variations in transit depth as a function of time due to stellar activity variability.


title: Nonminimally coupled scalar field in teleparallel gravity: boson stars

abstract: We study the nonminimally coupled complex scalar field within the framework of teleparallel gravity. Coupling of the field nonminimally to the torsion scalar destroys the Lorentz invariance of the theory in the sense that the resulting equations of motion depend on the choice of a tetrad. For the assumed static spherically symmetric spacetime, we find a tetrad which leads to a self-consistent set of equations, and we construct the self-gravitating configurations of the scalar field---boson stars. The resulting configurations develop anisotropic principal pressures and satisfy the dominant energy condition. An interesting property of the configurations obtained with sufficiently large field-to-torsion coupling constant is the outwardly increasing energy density, followed by an abrupt drop towards the usual asymptotic tail. This feature is not present in the boson stars with the field minimally or nonminimally coupled to the curvature scalar, and therefore appears to be a torsion--only effect.


title: Data-driven controller design for nonlinear systems: a two degrees of freedom architecture

abstract: In this paper, the D2-IBC (Data-Driven Inversion Based Control) approach for nonlinear control is introduced and analyzed. The method does not require any a-priori knowledge of the system dynamics and relies on a two degrees of freedom scheme, with a nonlinear controller and a linear controller running in parallel. In particular, the former is devoted to stabilize the system around a trajectory of interest, whereas the latter is used to boost the closed-loop performance. The paper also presents a thorough stability and performance analysis of the closed-loop system.


title: Occupation numbers in strongly polarized Fermi gases and the Luttinger theorem

abstract: We study a two-component Fermi gas that is so strongly polarized that it remains normal fluid at zero temperature. We calculate the occupation numbers within the particle-particle random-phase approximation, which is similar to the Nozieres-Schmitt-Rink approach. We show that the Luttinger theorem is fulfilled in this approach. We also study the change of the chemical potentials which allows us to extract, in the limit of extreme polarization, the polaron energy.


title: Non-destructive electron microscopy imaging and analysis of biological samples with graphene coating

abstract: In electron microscopy, charging of non-conductive biological samples by focused electron beams hinders their high-resolution imaging. Gold or platinum coatings have been commonly used to prevent such sample charging, but it disables further quantitative and qualitative chemical analyses by energy dispersive spectroscopy (EDS). Here we report that graphene-coating on biological samples enables non-destructive high-resolution imaging by scanning electron microscopy (SEM) as well as chemical analysis by EDS, utilizing graphene's transparency to electron beams, high conductivity, outstanding mechanical strength, and flexibility. We believe that the graphene-coated imaging and analysis would provide us a new opportunity to explore various biological phenomena unseen before due to the limitation in sample preparation and image resolution, which will broaden our understanding on the life mechanism of various living organisms.


title: Generalized classical dynamical Yang-Baxter equations and moduli spaces of flat connections on surfaces

abstract: In this paper, we explain how generalized dynamical r-matrices can be obtained by (quasi-)Poisson reduction. New examples of Poisson structures and Poisson groupoid actions naturally appear in this setting. As an application, we use a generalized dynamical r-matrix induced by the gauge fixing procedure to give a new finite dimensional description of the Atiyah-Bott symplectic structure on the moduli space of flat connections on a surface.


title: Engineering the microstructure and magnetism of La$_2$CoMnO$_6$ thin films by tailoring oxygen stoichiometry

abstract: We report on the magnetic and structural properties of ferromagnetic-insulating La$_2$CoMnO$_6$ thin films grown on top of (001) STO substrates by means of RF sputtering technique. Careful structural analysis, by using synchrotron X-ray diffraction, allows identifying two different crystallographic orientations that are closely related to oxygen stoichiometry and to the features (coercive fields and remanence) of the hysteresis loops. Both Curie temperature and magnetic hysteresis turn out to be dependent on the oxygen stoichiometry. In situ annealing conditions allow tailoring the oxygen content of the films, therefore controlling their microstructure and magnetic properties.


title: MIMTool: A Tool for Drawing Molecular Interaction Maps

abstract: Background: To understand protein function, it is important to study protein- protein interaction networks. These networks can be represented in network diagrams called protein interaction maps that can lead to better understanding by visualization. We address the problem of drawing of protein interactions in Kohn's Molecular Interaction Map (MIM) notation. Even though there are some existing tools for graphical visualization of protein interactions in general, there is no tool that can draw protein interactions with MIM notation with full support. Results: MIMTool was developed for drawing protein interaction maps in Kohn's MIM notation. MIMTool was developed using the Qt toolkit libraries and introduces several unique features such as full interactivity, object dragging, ability to export files in MIMML, SBML and line drawing with automatic bending and crossover minimization, which are not available in other diagram editors. MIMTool also has a unique orthogonal edge drawing method that is both easy and more flexible than other orthogonal drawing methods present in other interaction drawing tools. Conclusions: MIMTool facilitates faster drawing, updating and exchanging of MIMs. Among its several features, it also includes a semi-automatic drawing algorithm that makes use of shortest path algorithm for constructing lines with small number of bends and crossings. MIMTool contributes a much needed software tool that was missing and will facilitate wider adoption of Kohn's MIM notation.


title: Visualization of Large Volumetric Multi-Channel Microscopy Data Streams on Standard PCs

abstract: Background: Visualization of multi-channel microscopy data plays a vital role in biological research. With the ever-increasing resolution of modern microscopes the data set size of the scanned specimen grows steadily. On commodity hardware this size easily exceeds the available main memory and the even more limited GPU memory. Common volume rendering techniques require the entire data set to be present in the GPU memory. Existing out-of-core rendering approaches for large volume data sets either are limited to single-channel volumes, or require a computer cluster, or have long preprocessing times. Results: We introduce a ray-casting technique for rendering large volumetric multi-channel microscopy data streams on commodity hardware. The volumetric data is managed at different levels of detail by an octree structure. In contrast to previous octree-based techniques, the octree is built incrementally and therefore supports streamed microscopy data as well as data set sizes exceeding the available main memory. Furthermore, our approach allows the user to interact with the partially rendered data set at all stages of the octree construction. After a detailed description of our method, we present performance results for different multi-channel data sets with a size of up to 24 GB on a standard desktop PC. Conclusions: Our rendering technique allows biologists to visualize their scanned specimen on their standard desktop computers without high-end hardware requirements. Furthermore, the user can interact with the data set during the initial loading to explore the already loaded parts, change rendering parameters like color maps or adjust clipping planes. Thus, the time of biologists being idle is reduced. Also, streamed data can be visualized to detect and stop flawed scans early during the scan process.


title: An ansatz to the quantum phase transition in a dissipative two-qubit system

abstract: By means of a unitary transformation, we propose an ansatz to study quantum phase transitions in the ground state of a two-qubit system interacting with a dissipative reservoir. First, the ground state phase diagram is analyzed in the presence of the Ohmic and sub-Ohmic bath using an analytic ground state wave function which takes into account the competition between intrasite tunneling and intersite correlation. The quantum critical point is determined as the transition point from non-degenerate to degenerate ground state and our calculated critical coupling strength $\alpha_c$ agrees with that from the numerical renormalization group method. Moreover, by computing the entanglement entropy between the qubits and the bath as well as the qubit-qubit correlation function in the ground state, we explore the nature of the quantum phase transition between the delocalized and localized states.


title: Shifted convolution and the Titchmarsh divisor problem over F_q[t]

abstract: In this paper we solve a function field analogue of classical problems in analytic number theory, concerning the auto-correlations of divisor functions, in the limit of a large finite field.


title: A Cyber-Physical System-based Approach for Industrial Automation Systems

abstract: Industrial automation systems (IASs) are commonly developed using the languages defined by the IEC 61131 standard and are executed on PLCs. In this paper, a system-based approach for the development of IASs is adopted. A framework is described to refine the UML model of the software part, which is extracted from the SysML system model, and get the implementation code. Two implementation alternatives are considered to exploit PLCs but also the recent deluge of embedded boards in the market. For PLC targets, the new version of IEC 61131 that supports Object-Orientation is adopted, while Java is used for embedded boards. The case study was developed as a lab exercise for teaching the various technologies that address challenges in the domain of cyber-physical systems where Internet of Things (IoT) would be the glue regarding their cyber interfaces.


title: Shear accelerated crystallization in a supercooled atomic liquid

abstract: A bulk metallic glass forming alloy is subjected to shear flow in its supercooled state by compression of a short rod to produce a flat disc. The resulting material exhibits enhanced crystallization kinetics during isothermal annealing as reflected in the decrease of the crystallization time relative to the non-deformed case. The transition from quiescent to shear-accelerated crystallization is linked to strain accumulated during shear flow above a critical shear rate $\dot\gamma_c\approx 0.3$ s$^{-1}$ which corresponds to P\'{e}clet number, $Pe\sim\mathcal{O}(1)$. The observation of shear accelerated crystallization in an atomic system at modest shear rates is uncommon. It is made possible here by the substantial viscosity of the supercooled liquid which increases strongly with temperature in the approach to the glass transition. We may therefore anticipate the encounter of non-trivial shear-related effects during thermoplastic deformation of similar systems.


title: Enhancement of the lifetime of metastable states in Er-doped Si nanocrystals by external coloured noise

abstract: The changes in the lifetime of a metastable energy level in Er-doped Si nanocrystals in the presence of an external source of colored noise are analyzed for different values of noise intensity and correlation time. Exciton dynamics is simulated by a set of phenomenological rate equations which take into account all the possible phenomena inherent to the energy states of Si nanocrystals and Er$^{3+}$ ions in the host material of Si oxide. The electronic deexcitation is studied by examining the decay of the initial population of the Er atoms in the first excitation level $^4$I$_{13/2}$ through the fluorescence and the cooperative upconversion by energy transfer. Our results show that the deexcitation process of the level $^4$I$_{13/2}$ is slowed down within wide ranges of noise intensity and correlation time. Moreover, a nonmonotonic behavior of the lifetime with the amplitude of the fluctuations is found, characterized by a maximum variation for values of the noise correlation time comparable to the deexcitation time. The indirect influence of the colored noise on the efficiency of the energy transfer upconversion activated from the level $^4$I$_{13/2}$ is also discussed.


title: Solvable and/or integrable many-body models on a circle

abstract: Various many-body models are treated, which describe $N$ points confined to move on a plane circle. Their Newtonian equations of motion ("accelerations equal forces") are integrable, i. e. they allow the explicit exhibition of $N$ constants of motion in terms of the dependent variables and their time-derivatives. Some of these models are moreover solvable by purely algebraic operations, by (explicitly performable) quadratures and, finally, by functional inversions. The techniques to manufacture these models are not new; some of these models are themselves new; others are reinterpretations of known models.


title: The inner boundary of random walk range

abstract: In this paper, we deal with the inner boundary of random walk range, that is, the set of those points in a random walk range which have at least one neighbor site outside the range. If $L_n$ be the number of the inner boundary points of random walk range in the $n$ steps, we prove $\lim_{n\to \infty}\frac{L_n}{n}$ exists with probability one. Also, we obtain some large deviation result for transient walk. We find that the expectation of the number of the inner boundary points of simple random walk on two dimensionnal square lattice is of the same order as $\frac{n}{(\log n)^2}$.


title: FPGA Based Efficient Multiplier for Image Processing Applications Using Recursive Error Free Mitchell Log Multiplier and KOM Architecture

abstract: The Digital Image processing applications like medical imaging, satellite imaging, Biometric trait images etc., rely on multipliers to improve the quality of image. However, existing multiplication techniques introduce errors in the output with consumption of more time, hence error free high speed multipliers has to be designed. In this paper we propose FPGA based Recursive Error Free Mitchell Log Multiplier (REFMLM) for image Filters. The 2x2 error free Mitchell log multiplier is designed with zero error by introducing error correction term is used in higher order Karastuba-Ofman Multiplier (KOM) Architectures. The higher order KOM multipliers is decomposed into number of lower order multipliers using radix 2 till basic multiplier block of order 2x2 which is designed by error free Mitchell log multiplier. The 8x8 REFMLM is tested for Gaussian filter to remove noise in fingerprint image. The Multiplier is synthesized using Spartan 3 FPGA family device XC3S1500-5fg320. It is observed that the performance parameters such as area utilization, speed, error and PSNR are better in the case of proposed architecture compared to existing architectures


title: Last zero time or Maximum time of the winding number of Brownian motions

abstract: In this paper we consider the winding number, $\theta(s)$, of planar Brownian motion and study asymptotic behavior of the process of the maximum time, the time when $\theta(s)$ attains the maximum in the interval $0\le s \le t$. We find the limit law of its logarithm with a suitable normalization factor and the upper growth rate of the maximum time process itself. We also show that the process of the last zero time of $\theta(s)$ in $[0,t]$ has the same law as the maximum time process.


title: Analyzing Chromatin Using Tiled Binned Scatterplot Matrices

abstract: Background: Over the last years, more and more biological data became available. Besides the pure amount of new data, also its dimensionality - the number of different attributes per data point - increased. Recently, especially the amount of data on chromatin and its modifications increased considerably. In the field of epigenetics, appropriate visualization tools designed for highlighting the different aspects of epigenetic data are currently not available. Results: We present a tool called TiBi-Scatter enabling correlation analysis in 2D. This approach allows for analyzing multidimensional data while keeping the use of resources such as memory small. Thus, it is in particular applicable to large data sets. Conclusions: TiBi-Scatter is a resource-friendly and easy to use tool that allows for the hypothesis-free analysis of large multidimensional biological data sets.


title: The scientific potential of space-based gravitational wave detectors

abstract: The millihertz gravitational wave band can only be accessed with a space-based interferometer, but it is one of the richest in potential sources. Observations in this band have amazing scientific potential. The mergers between massive black holes with mass in the range 10 thousand to 10 million solar masses, which are expected to occur following the mergers of their host galaxies, produce strong millihertz gravitational radiation. Observations of these systems will trace the hierarchical assembly of structure in the Universe in a mass range that is very difficult to probe electromagnetically. Stellar mass compact objects falling into such black holes in the centres of galaxies generate detectable gravitational radiation for several years prior to the final plunge and merger with the central black hole. Measurements of these systems offer an unprecedented opportunity to probe the predictions of general relativity in the strong-field and dynamical regime. Millihertz gravitational waves are also generated by millions of ultra-compact binaries in the Milky Way, providing a new way to probe galactic stellar populations. ESA has recognised this great scientific potential by selecting The Gravitational Universe as its theme for the L3 large satellite mission, scheduled for launch in ~2034. In this article we will review the likely sources for millihertz gravitational wave detectors and describe the wide applications that observations of these sources could have for astrophysics, cosmology and fundamental physics.


title: Quantum Gravity and the Cosmological Constant Problem

abstract: A finite and unitary nonlocal formulation of quantum gravity is applied to the cosmological constant problem. The entire functions in momentum space at the graviton-standard model particle loop vertices generate an exponential suppression of the vacuum density and the cosmological constant to produce agreement with their observational bounds.


title: Riemannian Center of Mass and so called karcher mean

abstract: The Riemannian center of mass was constructed in [GrKa] (1973). In [GKR1, GKR2, Gr, Ka, BuKa] (1974-1981) it was successfully applied with more refined estimates. Probably in 1990 someone renamed it without justification into karcher mean and references to the older papers were omitted by those using the new name. As a consequence newcomers started to reprove results from the above papers. - Here I explain the older history.


title: Fractal Structure of Shortest Interaction Paths in Native Proteins and Determination of Residues on a Given Shortest Path

abstract: Fractal structure of shortest paths depends strongly on interresidue interaction cutoff distance. The dimensionality of shortest paths is calculated as a function of interaction cutoff distance. Shortest paths are self similar with a fractal dimension of 1.12 when calculated with step lengths larger than 6.8 {\AA}. Paths are multifractal below 6.8 {\AA}. The number of steps to traverse a shortest path is a discontinuous function of cutoff size at short cutoff values, showing abrupt decreases to smaller values as cutoff distance increases. As information progresses along the direction of a shortest path a large set of residues are affected because they are interacting neighbors to the residues of the shortest path. Thus, several residues are involved diffusively in information transport which may be identified with the present model. An algorithm is introduced to determine the residues of a given shortest path. The shortest path residues are the highly visited residues during information transport. These paths are shown to lie on the high entropy landscape of the protein where entropy is taken to increase with abundance of visits to nodes during signal transport.


title: Visualization and Correction of Automated Segmentation, Tracking and Lineaging from 5-D Stem Cell Image Sequences

abstract: Results: We present an application that enables the quantitative analysis of multichannel 5-D (x, y, z, t, channel) and large montage confocal fluorescence microscopy images. The image sequences show stem cells together with blood vessels, enabling quantification of the dynamic behaviors of stem cells in relation to their vascular niche, with applications in developmental and cancer biology. Our application automatically segments, tracks, and lineages the image sequence data and then allows the user to view and edit the results of automated algorithms in a stereoscopic 3-D window while simultaneously viewing the stem cell lineage tree in a 2-D window. Using the GPU to store and render the image sequence data enables a hybrid computational approach. An inference-based approach utilizing user-provided edits to automatically correct related mistakes executes interactively on the system CPU while the GPU handles 3-D visualization tasks. Conclusions: By exploiting commodity computer gaming hardware, we have developed an application that can be run in the laboratory to facilitate rapid iteration through biological experiments. There is a pressing need for visualization and analysis tools for 5-D live cell image data. We combine accurate unsupervised processes with an intuitive visualization of the results. Our validation interface allows for each data set to be corrected to 100% accuracy, ensuring that downstream data analysis is accurate and verifiable. Our tool is the first to combine all of these aspects, leveraging the synergies obtained by utilizing validation information from stereo visualization to improve the low level image processing tasks.


title: Some remarks on Lefschetz thimbles and complex Langevin dynamics

abstract: Lefschetz thimbles and complex Langevin dynamics both provide a means to tackle the numerical sign problem prevalent in theories with a complex weight in the partition function, e.g. due to nonzero chemical potential. Here we collect some findings for the quartic model, and for U(1) and SU(2) models in the presence of a determinant, which have some features not discussed before, due to a singular drift. We find evidence for a relation between classical runaways and stable thimbles, and give an example of a degenerate fixed point. We typically find that the distributions sampled in complex Langevin dynamics are related to the thimble(s), but with some important caveats, for instance due to the presence of unstable fixed points in the Langevin dynamics.


title: The plastic number and its generalized polynomial

abstract: The polynomial $X^{3}-X-1$ has a unique positive root known as plastic number, which is denoted by $\rho$ and is approximately equal to $1.32471795$. In this note we study the zeroes of the generalized polynomial $X^{k}-\sum_{j=0}^{k-2}X^{j}$ for $k\geq 3$ and prove that its unique positive root $\lambda_{k}$ tends to the golden ratio $\phi=\frac{1+\sqrt{5}}{2}$ as $k \to \infty$. We also derive bounds on $\lambda_{k}$ in terms of Fibonacci numbers.


title: On Hilbert's Tenth Problem

abstract: Using an iterated Horner schema for evaluation of diophantine polynomials, we define a partial $\mu$-recursive "decision" algorithm decis as a "race" for a first nullstelle versus a first (internal) proof of non-nullity for such a polynomial -- within a given theory T extending Peano Arithmetique PA. If T is diophantine sound, i.e., if (internal) provability implies truth -- for diophantine formulae --, then the T-map decis gives correct results when applied to the codes of polynomial inequalities $D(x_1,...,x_m) \neq 0$. The additional hypothesis that T be diophantine complete (in the syntactical sense) would guarantee in addition termination of decis on these formula, i.e., decis would constitute a decision algorithm for diophantine formulae in the sense of Hilbert's 10th problem. From Matiyasevich's impossibility for such a decision it follows, that a consistent theory T extending PA cannot be both diophantine sound and diophantine complete. We infer from this the existence of a diophantine formulae which is undecidable by T. Diophantine correctness is inherited by the diophantine completion T~ of T, and within this extension decis terminates on all externally given diophantine polynomials, correctly. Matiyasevich's theorem -- for the strengthening T~ of T -- then shows that T~, and hence T, cannot be diophantine sound. But since the internal consistency formula Con_T for T implies -- within PA -- diophantine soundness of T, we get that PA derives \neg Con_T, in particular PA must derive its own internal inconsistency formula.


title: Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

abstract: This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.


title: On the star discrepancy of sequences in the unit interval

abstract: It is known that there is a constant $c > 0$ such that for every sequence $x_1, x_2, \ldots$ in $[0,1)$ we have for the star discrepancy $D_N^*$ of the first $N$ elements of the sequence that $N D_N^* \ge c \cdot \log N$ holds for infinitely many $N$. Let $c^*$ be the supremum of all such $c$ with this property. We show $c^* > 0.0646363$, thereby improving the until now known estimates.


title: On Flexible Prismatic Polyhedra

abstract: We demonstrate the existence of four types of flexible prismatic polyhedra that can be derived or inferred from a consideration of Bricard octahedra and generalizations of Bricard octahedra. These flexible polyhedra are of genus 0 and 1, have dihedral angles all of which are non-constant under flexion, have faces that exhibit self-intersections, have an indefinite number of faces and have vertex motion that can treated as a function of a single variable.


title: Effect of parameter mismatch on the synchronization of strongly coupled self sustained oscillators

abstract: In this paper we present an experimental setup and an associated mathematical model to study the synchronization of two self sustained strongly coupled mechanical oscillators (metronomes). The effects of a small detuning in the internal parameters, namely damping and frequency, have been studied. Our experimental system is a pair of spring wound mechanical metronomes, coupled by placing them on a common base, free to move along a horizontal direction. In our system the mass of the oscillating pendula form a significant fraction of the total mass of the system, leading to strong coupling of the oscillators. We modified the internal mechanism of the spring-wound "clockwork" slightly, such that the natural frequency and the internal damping could be independently tuned. Stable synchronized and anti-synchronized states were observed as the difference in the parameters was varied. We designed a photodiode array based non-contact, non-magnetic position detection system driven by a microcontroller to record the instantaneous angular displacement of each oscillator and the small linear displacement of the base coupling the two. Our results indicate that such a system can be made to stabilize in both in-phase anti-phase synchronized state by tuning the parameter mismatch. Results from both numerical simulations and experimental observations are in qualitative agreement and are both reported in the present work.


title: Multi-level interference resonances in strongly-driven three-level systems

abstract: We study multi-photon resonances in a strongly-driven three-level quantum system, where one level is periodically swept through a pair of levels with constant energy separation $E$. Near the multi-photon resonance condition $n\hbar\omega = E$, where $n$ is an integer, we find qualitatively different behavior for $n$ even or odd. We explain this phenomenon in terms of families of interfering trajectories of the multi-level system. Remarkably, the behavior is insensitive to fluctuations of the energy of the driven level, and survives deep into the strong dephasing regime. The setup can be relevant for a variety of solid state and atomic or molecular systems. In particular, it provides a clear mechanism to explain recent puzzling experimental observations in strongly-driven double quantum dots.


title: inPHAP: Interactive visualization of genotype and phased haplotype data

abstract: Background: To understand individual genomes it is necessary to look at the variations that lead to changes in phenotype and possibly to disease. However, genotype information alone is often not sufficient and additional knowledge regarding the phase of the variation is needed to make correct interpretations. Interactive visualizations, that allow the user to explore the data in various ways, can be of great assistance in the process of making well informed decisions. But, currently there is a lack for visualizations that are able to deal with phased haplotype data. Results: We present inPHAP, an interactive visualization tool for genotype and phased haplotype data. inPHAP features a variety of interaction possibilities such as zooming, sorting, filtering and aggregation of rows in order to explore patterns hidden in large genetic data sets. As a proof of concept, we apply inPHAP to the phased haplotype data set of Phase 1 of the 1000 Genomes Project. Thereby, inPHAP's ability to show genetic variations on the population as well as on the individuals level is demonstrated for several disease related loci. Conclusions: As of today, inPHAP is the only visual analytical tool that allows the user to explore unphased and phased haplotype data interactively. Due to its highly scalable design, inPHAP can be applied to large datasets with up to 100 GB of data, enabling users to visualize even large scale input data. inPHAP closes the gap between common visualization tools for unphased genotype data and introduces several new features, such as the visualization of phased data.


title: Non-local kinetic and macroscopic models for self-organised animal aggregations

abstract: The last two decades have seen a surge in kinetic and macroscopic models derived to investigate the multi-scale aspects of self-organised biological aggregations. Because the individual-level details incorporated into the kinetic models (e.g., individual speeds and turning rates) make them somewhat difficult to investigate, one is interested in transforming these models into simpler macroscopic models, by using various scaling techniques that are imposed by the biological assumptions of the models. Here, we consider three scaling approaches (parabolic, hydrodynamic and grazing collision limits) that can be used to reduce a class of non-local 1D and 2D models for biological aggregations to simpler models existent in the literature. Next, we investigate how some of the spatio-temporal patterns exhibited by the original kinetic models are preserved via these scalings. To this end, we focus on the parabolic scaling for non-local 1D models and apply asymptotic preserving numerical methods, which allow us to analyse changes in the patterns as the scaling coefficient $\epsilon$ is varied from $\epsilon=1$ (for 1D transport models) to $\epsilon=0$ (for 1D parabolic models). We show that some patterns (describing stationary aggregations) are preserved in the limit $\epsilon\to 0$, while other patterns (describing moving aggregations) are lost in this limit. To understand the loss of these patterns, we construct bifurcation diagrams.


title: Large-order aspects of the delta-expansion in low-dimensional Ising models

abstract: We investigate the large order aspects of the delta-expansion under the estimation procession of the critical quantities. As illustrative examples, we revisit one-dimensional Ising model for the analytic study and two-dimensional square Ising model in the high temperature phase for the numerical experiment to large orders. In both models, proposed fundamental base on which the estimation protocol should be constructed is investigated in details and confirmed to be valid. In the square lattice model, we present a new protocol for the estimation of critical exponents and temperature.


title: eXamine: Exploring annotated modules in networks

abstract: Background: Biological networks have a growing importance for the interpretation of high-throughput omics data. Integrative network analysis makes use of statistical and combinatorial methods to extract smaller subnetwork modules, and performs enrichment analysis to annotate the modules with ontology terms or other available knowledge. This process results in an annotated module, which retains the original network structure and includes enrichment information as a set system. A major bottleneck is a lack of tools that allow exploring both network structure of extracted modules and its annotations. Results: Thispaperpresentsavisualanalysisapproachthattargetssmallmoduleswithmanyset-based annotations, and which displays the annotations as contours on top of a node-link diagram. We introduce an extension of self-organizing maps to lay out nodes, links, and contours in a unified way. An implementation of this approach is freely available as the Cytoscape app eXamine. Conclusions: eXamine accurately conveys small and annotated modules consisting of several dozens of proteins and annotations. We demonstrate that eXamine facilitates the interpretation of integrative network analysis results in a guided case study. This study has resulted in a novel biological insight regarding the virally-encoded G-protein coupled receptor US28.


title: A new approach to determine optically thick H2 cooling and its effect on primordial star formation

abstract: We present a new method for estimating the H2 cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities, to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability, because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of 40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disc has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above 10^{10}/ccm at a distance of 200AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation.


title: On a class of biorthogonal polynomials on the unit circle

abstract: A system of biorthogonal polynomials with respect to a complex valued measure supported on the unit circle is considered and all the terms with bounds are explicitly given for the remainder of an asymptotic formula given by R. Askey for this system. An electrostatic interpretation for the zeros of a class of para-orthogonal polynomials associated with the biorthogonal system is also considered.


title: Decomposition with respect to outputs for Boolean control networks

abstract: This paper investigates the problem of decomposition with respect to outputs for Boolean control networks (BCNs). First, with the linear expression of BCNs and the matrix semi-tensor product, some algebraic equivalent conditions for the decomposition are obtained. Second, a necessary and sufficient graphical condition for the decomposition with respect to outputs is given. Third, an effective method is proposed to reduce the computational burden in the realization of the decomposition. Finally, some examples are addressed to validate the effectiveness of the proposed method.


title: Bound states and expansion dynamics of interacting bosons on a one-dimensional lattice

abstract: The expansion dynamics of bosonic gases in optical lattices has recently been the focus of increasing attention, both experimental and theoretical. We consider, by means of numerical Bethe ansatz, the expansion dynamics of initially confined wave packets of two interacting bosons on a lattice. We show that a correspondence between the asymptotic expansion velocities and the projection of the evolved wave function over the bound states of the system exists, clarifying the existing picture for such situations. Moreover, we investigate the role of the lattice in this kind of evolution.


title: Checking Termination of Bottom-Up Evaluation of Logic Programs with Function Symbols

abstract: Recently, there has been an increasing interest in the bottom-up evaluation of the semantics of logic programs with complex terms. The presence of function symbols in the program may render the ground instantiation infinite, and finiteness of models and termination of the evaluation procedure, in the general case, are not guaranteed anymore. Since the program termination problem is undecidable in the general case, several decidable criteria (called program termination criteria) have been recently proposed. However, current conditions are not able to identify even simple programs, whose bottom-up execution always terminates. The paper introduces new decidable criteria for checking termination of logic programs with function symbols under bottom-up evaluation, by deeply analyzing the program structure. First, we analyze the propagation of complex terms among arguments by means of the extended version of the argument graph called propagation graph. The resulting criterion, called Gamma-acyclicity, generalizes most of the decidable criteria proposed so far. Next, we study how rules may activate each other and define a more powerful criterion, called safety. This criterion uses the so-called safety function able to analyze how rules may activate each other and how the presence of some arguments in a rule limits its activation. We also study the application of the proposed criteria to bound queries and show that the safety criterion is well-suited to identify relevant classes of programs and bound queries. Finally, we propose a hierarchy of classes of terminating programs, called k-safety, where the k-safe class strictly includes the (k-1)-safe class. Note: To appear in Theory and Practice of Logic Programming (TPLP).


title: iGPSe: A Visual Analytic System for Integrative Genomic Based Cancer Patient Stratification

abstract: Background: Cancers are highly heterogeneous with different subtypes. These subtypes often possess different genetic variants, present different pathological phenotypes, and most importantly, show various clinical outcomes such as varied prognosis and response to treatment and likelihood for recurrence and metastasis. Recently, integrative genomics (or panomics) approaches are often adopted with the goal of combining multiple types of omics data to identify integrative biomarkers for stratification of patients into groups with different clinical outcomes. Results: In this paper we present a visual analytic system called Interactive Genomics Patient Stratification explorer (iGPSe) which significantly reduces the computing burden for biomedical researchers in the process of exploring complicated integrative genomics data. Our system integrates unsupervised clustering with graph and parallel sets visualization and allows direct comparison of clinical outcomes via survival analysis. Using a breast cancer dataset obtained from the The Cancer Genome Atlas (TCGA) project, we are able to quickly explore different combinations of gene expression (mRNA) and microRNA features and identify potential combined markers for survival prediction. Conclusions: Visualization plays an important role in the process of stratifying given population patients. Visual tools allowed for the selection of possibly features across various datasets for the given patient population. We essentially made a case for visualization for a very important problem in translational informatics.


title: An error analysis for polynomial optimization over the simplex based on the multivariate hypergeometric distribution

abstract: We study the minimization of fixed-degree polynomials over the simplex. This problem is well-known to be NP-hard, as it contains the maximum stable set problem in graph theory as a special case. In this paper, we consider a rational approximation by taking the minimum over the regular grid, which consists of rational points with denominator $r$ (for given $r$). We show that the associated convergence rate is $O(1/r^2)$ for quadratic polynomials. For general polynomials, if there exists a rational global minimizer over the simplex, we show that the convergence rate is also of the order $O(1/r^2)$. Our results answer a question posed by De Klerk et al. (2013) and improves on previously known $O(1/r)$ bounds in the quadratic case.


title: Planar Graphs: Random Walks and Bipartiteness Testing

abstract: We initiate the study of property testing in arbitrary planar graphs. We prove that bipartiteness can be tested in constant time, improving on the previous bound of $\tilde{O}(\sqrt{n})$ for graphs on $n$ vertices. The constant-time testability was only known for planar graphs with bounded degree. Our algorithm is based on random walks. Since planar graphs have good separators, i.e., bad expansion, our analysis diverges from standard techniques that involve the fast convergence of random walks on expanders. We reduce the problem to the task of detecting an odd-parity cycle in a multigraph induced by constant-length cycles. We iteratively reduce the length of cycles while preserving the detection probability, until the multigraph collapses to a collection of easily discoverable self-loops. Our approach extends to arbitrary minor-free graphs. We also believe that our techniques will find applications to testing other properties in arbitrary minor-free graphs.


title: Addressing the unmet need for visualizing Conditional Random Fields in Biological Data

abstract: Background: The biological world is replete with phenomena that appear to be ideally modeled and analyzed by one archetypal statistical framework - the Graphical Probabilistic Model (GPM). The structure of GPMs is a uniquely good match for biological problems that range from aligning sequences to modeling the genome-to-phenome relationship. The fundamental questions that GPMs address involve making decisions based on a complex web of interacting factors. Unfortunately, while GPMs ideally fit many questions in biology, they are not an easy solution to apply. Building a GPM is not a simple task for an end user. Moreover, applying GPMs is also impeded by the insidious fact that the complex web of interacting factors inherent to a problem might be easy to define and also intractable to compute upon. Discussion: We propose that the visualization sciences can contribute to many domains of the bio-sciences, by developing tools to address archetypal representation and user interaction issues in GPMs, and in particular a variety of GPM called a Conditional Random Field(CRF). CRFs bring additional power, and additional complexity, because the CRF dependency network can be conditioned on the query data. Conclusions: In this manuscript we examine the shared features of several biological problems that are amenable to modeling with CRFs, highlight the challenges that existing visualization and visual analytics paradigms induce for these data, and document an experimental solution called StickWRLD which, while leaving room for improvement, has been successfully applied in several biological research projects.


title: Tube filling technique influence on superconductivity in MgB2-based conductors prepared using laboratory made boron

abstract: Here we report a systematic study of superconductive properties of mono-filamentary MgB2-based wires, manufactured with different filling technique. A detailed comparison of the influence of filling technique and final heat treatment on superconducting properties has been given. Boron used was synthesized in laboratory following magnesiothermic reduction of boron oxide, purified thanks to several acid leaching and heat treated at high temperature, to enhance crystalline degree. Critical current density and critical temperature have been investigated at different synthesis/sintering temperature in order to evaluate their dependence to applied final heat treatment. Critical current density has been evaluated on short wire pieces by magnetic measurement in a MPMS 5.5T Quantum Design SQUID, while critical temperature has been measured with a four probe system, by the drop of resistivity during the cooling of the sample in a liquid helium dewar.


title: MCA: Multiresolution Correlation Analysis, a graphical tool for subpopulation identification in single-cell gene expression data

abstract: Background: Biological data often originate from samples containing mixtures of subpopulations, corresponding e.g. to distinct cellular phenotypes. However, identification of distinct subpopulations may be difficult if biological measurements yield distributions that are not easily separable. Results: We present Multiresolution Correlation Analysis (MCA), a method for visually identifying subpopulations based on the local pairwise correlation between covariates, without needing to define an a priori interaction scale. We demonstrate that MCA facilitates the identification of differentially regulated subpopulations in simulated data from a small gene regulatory network, followed by application to previously published single-cell qPCR data from mouse embryonic stem cells. We show that MCA recovers previously identified subpopulations, provides additional insight into the underlying correlation structure, reveals potentially spurious compartmentalizations, and provides insight into novel subpopulations. Conclusions: MCA is a useful method for the identification of subpopulations in low-dimensional expression data, as emerging from qPCR or FACS measurements. With MCA it is possible to investigate the robustness of covariate correlations with respect subpopulations, graphically identify outliers, and identify factors contributing to differential regulation between pairs of covariates. MCA thus provides a framework for investigation of expression correlations for genes of interests and biological hypothesis generation.


title: A Strained Silicon Cold Electron Bolometer using Schottky Contacts

abstract: We describe optical characterisation of a Strained Silicon Cold Electron Bolometer (CEB), operating on a $350~\mathrm{mK}$ stage, designed for absorption of millimetre-wave radiation. The silicon Cold Electron Bolometer utilises Schottky contacts between a superconductor and an n++ doped silicon island to detect changes in the temperature of the charge carriers in the silicon, due to variations in absorbed radiation. By using strained silicon as the absorber, we decrease the electron-phonon coupling in the device and increase the responsivity to incoming power. The strained silicon absorber is coupled to a planar aluminium twin-slot antenna designed to couple to $160~\mathrm{GHz}$ and that serves as the superconducting contacts. From the measured optical responsivity and spectral response, we calculate a maximum optical efficiency of $50~\%$ for radiation coupled into the device by the planar antenna and an overall noise equivalent power (NEP), referred to absorbed optical power, of $1.1 \times 10^{-16}~\mathrm{\mbox{W Hz}^{-1/2}}$ when the detector is observing a $300~\mathrm{K}$ source through a $4~\mathrm{K}$ throughput limiting aperture. Even though this optical system is not optimised we measure a system noise equivalent temperature difference (NETD) of $6~\mathrm{\mbox{mK Hz}^{-1/2}}$. We measure the noise of the device using a cross-correlation of time stream data measured simultaneously with two junction field-effect transistor (JFET) amplifiers, with a base correlated noise level of $300~\mathrm{\mbox{pV Hz}^{-1/2}}$ and find that the total noise is consistent with a combination of photon noise, current shot noise and electron-phonon thermal noise.


title: Liquid-induced damping of mechanical feedback effects in single electron tunneling through a suspended carbon nanotube

abstract: In single electron tunneling through clean, suspended carbon nanotube devices at low temperature, distinct switching phenomena have regularly been observed. These can be explained via strong interaction of single electron tunneling and vibrational motion of the nanotube. We present measurements on a highly stable nanotube device, subsequently recorded in the vacuum chamber of a dilution refrigerator and immersed in the 3He/4He mixture of a second dilution refrigerator. The switching phenomena are absent when the sample is kept in the viscous liquid, additionally supporting the interpretation of dc-driven vibration. Transport measurements in liquid helium can thus be used for finite bias spectroscopy where otherwise the mechanical effects would dominate the current.


title: Measurement of high-k density fluctuation wavenumber spectrum in MAST and Doppler backscattering for spherical tokamaks

abstract: The high-k ($7 \lesssim k_{\bot} \rho_i \lesssim 11$) wavenumber spectrum of density fluctuations has been measured for the first time in MAST [B. Lloyd et al, Nucl. Fusion 43, 1665 (2003)]. This was accomplished with the first implementation of Doppler backscattering (DBS) for core measurements in a spherical tokamak. DBS has become a well-established and versatile diagnostic technique for the measurement of intermediate-k ($k_{\bot} \rho_i \sim 1$, and higher) density fluctuations and flows in magnetically confined fusion experiments. A novel implementation with 2D steering was necessary to enable DBS measurements in MAST, where the large magnetic field pitch angle presents a challenge. We report on the scattering considerations and ray tracing calculations used to optimize the design and present data demonstrating measurement capabilities. Initial results confirm the applicability of the design and implementation approaches, showing the strong dependence of scattering alignment on toroidal launch angle. We also present comparisons of DBS plasma velocity measurements with charge exchange recombination and beam emission spectroscopy measurements, which show reasonable agreement over most of the minor radius, but imply large poloidal flows approaching the magnetic axis in a discharge with an internal transport barrier. The 2D steering is shown to enable high-k measurements with DBS, at $k_{\bot}>20\ \mathrm{cm}^{-1}$ ($k_{\bot} \rho_i>10$) for launch frequencies less than 75 GHz; this capability is used to measure the wavenumber spectrum of turbulence and we find $|n(k_{\bot})|^2 \propto k_{\bot}^{- 4.7 \pm 0.2}$ for $k_{\bot} \rho_i \approx 7-11$, which is similar to the expectation for the turbulent kinetic cascade of $|n(k_{\bot})|^2 \propto k_{\bot}^{- 13/3}$.


title: On the discretization of nonholonomic dynamics in $\mathbb{R}^n$

abstract: In this paper we explore the nonholonomic Lagrangian setting of mechanical systems in local coordinates on finite-dimensional configuration manifolds. We prove existence and uniqueness of solutions by reducing the basic equations of motion to a set of ordinary differential equations on the underlying distribution manifold $D$. Moreover, we show that any $D-$preserving discretization may be understood as beeing generated by the exact evolution map of a time-periodic non-autonomous perturbation of the original continuous-time nonholonomic system. By means of discretizing the corresponding Lagrange-d'Alembert principle, we construct geometric integrators for the original nonholonomic system. We give precise conditions under which these integrators generate a discrete flow preserving the distribution $D$. Also, we derive corresponding consistency estimates. Finally, we carefully treat the example of the nonholonomic particle, showing how to discretize the equations of motion in a reasonable way, particularly regarding the nonholonomic constraints. The exploration in this paper lays the ground to analyze the dynamics of appropriate discretizations of nonholonomic mechnical systems in the Lagrangian framework and to relate that dynamics to the dynamics of the original nonholonomic systems. We postpone this analysis to a series of forthcoming papers.


title: Visualization of gene expression information within the context of the mouse anatomy

abstract: Background: The eMouse Atlas of Gene Expression (EMAGE) is an online resource that publishes the results of in situ gene expression experiments on the developmental mouse. The resource provides comprehensive search facilities, but few analytical tools or visual mechanisms for navigating the data set. To deal with the missing visual navigation, this paper explores the application of sunburst and icicle visualizations within EMAGE. Results: A prototype solution delivered a simple user interface that helps the user query EMAGE and generate a sunburst/icicle diagram. An evaluation featuring test subjects from the EMAGE staff studied the visualizations and provided a range of suggested improvements. Moreover the evaluation discovered that in addition to providing a visual means of walking through the data, when grouped, the sunburst delivers an interactive overview that assists with analysing sets of related genes. Conclusions: The sunburst and icicle visualizations have been shown to be effective tools for summarising gene expression data. The sunburst with its space saving radial layout was found especially useful for providing an overview of gene families or pathways. Work is ongoing to integrate these visualizations into EMAGE.


title: Mid-Infrared Spectroscopy of Uranus from the Spitzer Infrared Spectrometer: 2. Determination of the Mean Composition of the Upper Troposphere and Stratosphere

abstract: Mid-infrared spectral observations Uranus acquired with the Infrared Spectrometer (IRS) on the Spitzer Space Telescope are used to determine the abundances of C2H2, C2H6, CH3C2H, C4H2, CO2, and tentatively CH3 on Uranus at the time of the 2007 equinox. For vertically uniform eddy diffusion coefficients in the range 2200-2600 cm2 s-1, photochemical models that reproduce the observed methane emission also predict C2H6 profiles that compare well with emission in the 11.6-12.5 micron wavelength region, where the nu9 band of C2H6 is prominent. Our nominal model with a uniform eddy diffusion coefficient Kzz = 2430 cm2 sec-1 and a CH4 tropopause mole fraction of 1.6x10-5 provides a good fit to other hydrocarbon emission features, such as those of C2H2 and C4H2, but the model profile for CH3C2H must be scaled by a factor of 0.43, suggesting that improvements are needed in the chemical reaction mechanism for C3Hx species. The nominal model is consistent with a CH3D/CH4 ratio of 3.0+-0.2x10-4. From the best-fit scaling of these photochemical-model profiles, we derive column abundances above the 10-mbar level of 4.5+01.1/-0.8 x 10+19 molecule-cm-2 for CH4, 6.2 +- 1.0 x 10+16 molecule-cm-2 for C2H2 (with a value 24% higher from a different longitudinal sampling), 3.1 +- 0.3 x 10+16 molecule-cm-2 for C2H6, 8.6 +- 2.6 x 10+13 molecule-cm-2 for CH3C2H, 1.8 +- 0.3 x 10+13 molecule-cm-2 for C4H2, and 1.7 +- 0.4 x 10+13 molecule-cm-2 for CO2 on Uranus. Our results have implications with respect to the influx rate of exogenic oxygen species and the production rate of stratospheric hazes on Uranus, as well as the C4H2 vapor pressure over C4H2 ice at low temperatures.