title: EIT and TRACE responses to flare plasma

abstract: Aims: To understand the contribution of active region and flare plasmas to the $\lambda$195 channels of SOHO/EIT (Extreme-ultraviolet Imaging Telescope) and TRACE (Transition Region and Coronal Explorer). Methods: We have analysed an M8 flare simultaneously observed by the Coronal Diagnostic Spectrometer (CDS), EIT, TRACE and RHESSI. We obtained synthetic spectra for the flaring region and an outer region using the differential emission measures (DEM) of emitting plasma based on CDS and RHESSI observations and the CHIANTI atomic database. We then predicted the EIT and TRACE count rates. Results: For the flaring region, both EIT and TRACE images taken through the $\lambda$195 filter are dominated by Fe ${\rm XXIV}$ (formed at about 20 MK). However, in the outer region, the emission was primarily due to the Fe${\rm XII}$, with substantial contributions from other lines. The average count rate for the outer region was within 25% the observed value for EIT, while for TRACE it was a factor of two higher. For the flare region, the predicted count rate was a factor of two (in case of EIT) and a factor of three (in case of TRACE) higher than the actual count rate. Conclusions: During a solar flare, both TRACE and EIT $\lambda$195 channels are found to be dominated by Fe ${\rm XXIV}$ emission. Reasonable agreement between predictions and observations is found, however some discrepancies need to be further investigated.


title: On the interpretative essence of the term "interaction-free measurement": The role of entanglement

abstract: The polemical term "interaction-free measurement" (IFM) is analyzed in its interpretative nature. Two seminal works proposing the term are revisited and their underlying interpretations are assessed. The role played by nonlocal quantum correlations (entanglement) is formally discussed and some controversial conceptions in the original treatments are identified. As a result the term IFM is shown to be consistent neither with the standard interpretation of quantum mechanics nor with the lessons provided by the EPR debate.


title: Vortices in brain waves

abstract: Interactions by mutual excitation in neural populations in human and animal brains create a mesoscopic order parameter that is recorded in brain waves (electroencephalogram, EEG). Spatially and spectrally distributed oscillations are imposed on the background activity by inhibitory feedback in the gamma range (30-80 Hz). Beats recur at theta rates (3-7 Hz), at which the order parameter transiently approaches zero and microscopic activity becomes disordered. After these null spikes, the order parameter resurges and initiates a frame bearing a mesoscopic spatial pattern of gamma amplitude modulation that governs the microscopic activity, and that is correlated with behavior. The brain waves also reveal a spatial pattern of phase modulation in the form of a cone. Using the formalism of the dissipative many-body model of brain, we describe the null spikes and the accompanying phase cones as vortices.


title: The Discrete Hilbert Transform for Non-Periodic Signals

abstract: This note investigates the size of the guard band for non-periodic discrete Hilbert transform, which has recently been proposed for data hiding and security applications. It is shown that a guard band equal to the duration of the message is sufficient for a variety of analog signals and is, therefore, likely to be adequate for discrete or digital data.


title: Oscillations of tori in the pseudo-Newtonian potential

abstract: Context. The high-frequency quasi-periodic oscillations (HF QPOs) in neutron star and stellar-mass black hole X-ray binaries may be the result of a resonance between the radial and vertical epicyclic oscillations in strong gravity. Aims. In this paper we investigate the resonant coupling between the epicyclic modes in a torus in a strong gravitational field. Methods. We perform numerical simulations of axisymmetric constant angular momentum tori in the pseudo-Newtonian potential. The epicyclic motion is excited by adding a constant radial velocity to the torus. Results. We verify that slender tori perform epicyclic motions at the frequencies of free particles, but the epicyclic frequencies decrease as the tori grow thicker. More importantly, and in contrast to previous numerical studies, we do not find a coupling between the radial and vertical epicyclic motions. The appearance of other modes than the radial epicyclic motion in our simulations is rather due to small numerical deviations from exact equilibrium in the initial state of our torus. Conclusions. We find that there is no pressure coupling between the two axisymmetric epicyclic modes as long as the torus is symmetric with respect to the equatorial plane. However we also find that there are other modes in the disc that may be more attractive for explaining the HF QPOs.


title: Three lectures on classical integrable systems and gauge field theories

abstract: In these lectures I consider the Hitchin integrable systems and their relations with the self-duality equations and the twisted super-symmetric Yang-Mills theory in four dimension follow Hitchin and Kapustin-Witten. I define the Symplectic Hecke correspondence between different integrable systems. As an example I consider Elliptic Calogero-Moser system and integrable Euler-Arnold top on coadjoint orbits of the group GL(N,C) and explain the Symplectic Hecke correspondence for these systems.


title: Constraints on the rare tau decays from mu --> e gamma in the supersymmetric see-saw model

abstract: It is now a firmly established fact that all family lepton numbers are violated in Nature. In this paper we discuss the implications of this observation for future searches for rare tau decays in the supersymmetric see-saw model. Using the two loop renormalization group evolution of the soft terms and the Yukawa couplings we show that there exists a lower bound on the rate of the rare process mu --> e gamma of the form BR(mu --> e gamma) > C BR(tau --> mu gamma) BR(tau --> e gamma), where C is a constant that depends on supersymmetric parameters. Our only assumption is the absence of cancellations among the high-energy see-saw parameters. We also discuss the implications of this bound for future searches for rare tau decays. In particular, for large regions of the mSUGRA parameter space, we show that present B-factories could discover either tau --> mu gamma or tau --> e gamma, but not both.


title: A homotopical algebra of graphs related to zeta series

abstract: The purpose of this paper is to develop a homotopical algebra for graphs, relevant to zeta series and spectra of finite graphs. More precisely, we define a Quillen model structure in a category of graphs (directed and possibly infinite, with loops and multiple arcs allowed). The weak equivalences for this model structure are the Acyclics (graph morphisms which preserve cycles). The cofibrations and fibrations for the model are determined from the class of Whiskerings (graph morphisms produced by grafting trees). Our model structure seems to fit well with the importance of acyclic directed graphs in many applications. In addition to the weak factorization systems which form this model structure, we also describe two Freyd-Kelly factorization systems based on Folding, Injecting, and Covering graph morphisms.


title: Separating NOF communication complexity classes RP and NP

abstract: We provide a non-explicit separation of the number-on-forehead communication complexity classes RP and NP when the number of players is up to \delta log(n) for any \delta<1. Recent lower bounds on Set-Disjointness [LS08,CA08] provide an explicit separation between these classes when the number of players is only up to o(loglog(n)).


title: The open mapping theorem for regular quaternionic functions

abstract: The basic results of a new theory of regular functions of a quaternionic variable have been recently stated, following an idea of Cullen. In this paper we prove the minimum modulus principle and the open mapping theorem for regular functions. The proofs involve some peculiar geometric properties of such functions which are of independent interest.


title: Process POVM: A mathematical framework for the description of process tomography experiments

abstract: In this paper we shall introduce the mathematical framework for the description of measurements of quantum processes. Using this framework the process estimation problems can be treated in the similar way as the state estimation problems, only replacing the concept of positive operator valued measure (POVM) by the concept of process POVM (PPOVM). In particular, we will show that any measurement of qudit channels can be described by a collection of effects (positive operators) defined on two-qudit system. However, the effects forming a PPOVM are not normalized in the usual sense. We will demonstrate the usage of this formalism in discrimination problems by showing that perfect channel discrimination is equivalent to a specific unambiguous state discrimination.


title: A novel technique to make Ohmic contact to a buried two-dimensional electron gas in a molecular-beam-epitaxy grown $GaAs/Al_{0.3}Ga_{0.7}As$ heterostructure with Mn $\delta$-doping

abstract: We report on the growth and characterization of a new Diluted Magnetic Semiconductor (DMS) heterostructure that presents a Two-Dimensional Electron Gas (2DEG) with a carrier density $n \sim 1.08 \times 10^{12} cm^{-2}$ and a mobility $\mu \sim 600 cm^{2} / (Vs)$ at T $\sim$ 4.2K. As far as we know this is the highest mobility value reported in the literature for GaMnAs systems. A novel technique was developed to make Ohmic contact to the buried 2DEG without destroying the magnetic properties of our crystal.


title: HST/STIS Optical Transit Transmission Spectra of the hot-Jupiter HD209458b

abstract: We present the transmission spectra of the hot-Jupiter HD209458b taken with the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. Our analysis combines data at two resolutions and applies a complete pixel-by-pixel limb-darkening correction to fully reveal the spectral line shapes of atmospheric absorption features. Terrestrial-based Na I and H I contamination are identified which mask the strong exoplanetary absorption signature in the Na core, which we find reaches total absorption levels of ~0.11% in a 4.4 Ang band. The Na spectral line profile is characterized by a wide absorption profile at the lowest absorption depths, and a sharp transition to a narrow absorption profile at higher absorption values. The transmission spectra also shows the presence of an additional absorber at ~6,250 Ang, observed at both medium and low resolutions. We performed various limb-darkening tests, including using high precision limb-darkening measurements of the sun to characterize a general trend of Atlas models to slightly overestimate the amount of limb-darkening at all wavelengths, likely due to the limitations of the model's one-dimensional nature. We conclude that, despite these limitations, Atlas models can still successfully model limb-darkening in high signal-to-noise transits of solar-type stars, like HD209458, to a high level of precision over the entire optical regime (3,000-10,000 Ang) at transit phases between 2nd and 3rd contact.


title: Generalization of One of Lie's Theorems

abstract: We prove a generalization of one of Lie's Theorems in the context of Lie-like algebras$^{2-nd}$.


title: Influence of Non-Uniform Distribution of Acoustic Wavefield Strength on Time-Distance Helioseismology Measurements

abstract: By analyzing numerically simulated solar oscillation data, we study the influence of non-uniform distribution of acoustic wave amplitude, acoustic source strength, and perturbations of the sound speed on the shifts of acoustic travel times measured by the time-distance helioseismology method. It is found that for short distances, the contribution to the mean travel time shift caused by non-uniform distribution of acoustic sources in sunspots may be comparable to (but smaller than) the contribution from the sound speed perturbation in sunspots, and that it has the opposite sign to the sound-speed effect. This effect may cause some underestimation of the negative sound-speed perturbations in sunspots just below the surface, that was found in previous time-distance helioseismology inferences. This effect cannot be corrected by artificially increasing the amplitude of oscillations in sunspots. For large time-distance annuli, the non-uniform distribution of wavefields does not have significant effects on the mean travel times, and thus the sound-speed inversion results. The measured travel time differences, which are used to determine the mass flows beneath sunspots, can also be systematically shifted by this effect, but only in an insignificant magnitude.


title: Non-uniform liquid--crystalline phases of parallel hard rod-shaped particles: From ellipsoids to cylinders

abstract: In this article we consider systems of parallel hard {\it superellipsoids}, which can be viewed as a possible interpolation between ellipsoids of revolution and cylinders. Superellipsoids are characterized by an aspect ratio and an exponent $\alpha$ (shape parameter) which takes care of the geometry, with $\alpha=1$ corresponding to ellipsoids of revolution, while $\alpha=\infty$ is the limit of cylinders. It is well known that, while hard parallel cylinders exhibit nematic, smectic, and solid phases, hard parallel ellipsoids do not stabilize the smectic phase, the nematic phase transforming directly into a solid as density is increased. We use computer simulation to find evidence that for $\alpha\ge\alpha_c$, where $\alpha_c$ is a critical value which the simulations estimate to be in the interval 1.2--1.3, the smectic phase is stabilized. This is surprisingly close to the ellipsoidal case. In addition, we use a density-functional approach, based on the Parsons--Lee approximation, to describe smectic and columnar ordering. In combination with a free--volume theory for the crystalline phase, a theoretical phase diagram is predicted. While some qualitative features, such as the enhancement of smectic stability for increasing $\alpha$, and the probable absence of a stable columnar phase, are correct, the precise location of coexistence densities are quantitatively incorrect.


title: Early-A stars from IPHAS, and their distribution in and around the Cyg OB2 association

abstract: Stellar photometry derived from the INT/WFC Photometric H$\alpha$ Survey of the Northern Galactic Plane (IPHAS) can be used to identify large, reliable samples of A0-A5 dwarfs. For every A star, so identified, it is also possible to derive individual reddening and distance estimates, under the assumption that most selected objects are on or near the main sequence, at a mean absolute r' magnitude of 1.5 -- 1.6. This study presents the method for obtaining such samples and shows that the known reddenings and distances to the open clusters NGC 7510 and NGC 7790 are successfully recovered. A sample of over 1000 A stars is then obtained from IPHAS data in the magnitude range 13.5 < r' < 20 from the region of sky including the massive northern OB association Cyg OB2. Analysis of these data reveals a concentration of ~200 A stars over an area about a degree across, offset mainly to the south of the known 1--3 Myr old OB stars in Cyg OB2: their dereddened r' magnitudes fall in the range 11.8 to 12.5. These are consistent with a ~7 Myr old stellar population at DM = 10.8, or with an age of ~5 Myr at DM = 11.2. The number of A stars found in this clustering alone is consistent with a lower limit to the cluster mass of ~10000 M-sun.


title: Superspace formulations of ten-dimensional supergravity

abstract: We present a new formulation for N=1, D=10 supergravity in superspace, in presence of a Lorentz Chern-Simons-form. This formulation entails the following properties: it furnishes a solution of the Bianchi identities that is algebraically consistent to all orders in alpha'; at first order it is the simplest formulation proposed so far, and it is therefore most suitable for an explicit higher order analysis; it allows a well defined perturbative expansion in alpha', in which no poltergeist fields appear; it reconciles the two different classes of first order solutions available in the literature, that until now appeared physically inequivalent.


title: Gradient NLW on curved background in 4+1 dimensions

abstract: We obtain a sharp local well-posedness result for the Gradient Nonlinear Wave Equation on a nonsmooth curved background. In the process we introduce variable coefficient versions of Bourgain's $X^{s,b}$ spaces, and use a trilinear multiscale wave packet decomposition in order to prove a key trilinear estimate.


title: Electron-mediated ferromagnetism and small spin-orbit interaction in a molecular-beam-epitaxy grown n-type $GaAs/Al_{0.3}Ga_{0.7}As$ heterostructure with Mn $\delta$-doping

abstract: We report the first evidence of electron-mediated ferromagnetism in a molecular-beam-epitaxy (MBE) grown $GaAs/Al_{0.3}Ga_{0.7}As$ heterostructure with Mn $\delta$-doping. The interaction between the magnetic dopants (Mn) and the Two-Dimensional Electron Gas (2DEG) realizes magnetic ordering when the temperature is below the Curie temperature ($T_{C} \sim 1.7K$) and the 2DEG is brought in close proximity to the Mn layer by gating. The Anomalous Hall Effect (AHE) contribution to the total Hall resistance is shown to be about three to four orders of magnitude smaller than in the case of hole-mediated ferromagnetism indicating the presence of small spin-orbit interaction.


title: On the endomorphisms of Weyl modules over affine Kac-Moody algebras at the critical level

abstract: We present an independent short proof of the main result of arXiv:0706.3725 that the algebra of endomorphisms of a Weyl module of critical level is isomorphic to the algebra of functions on the space of monodromy-free opers on the disc with regular singularity and residue determined by the highest weight of the Weyl module. We derive this from the results of arXiv:0712.1183 about the shift of argument subalgebras.


title: Powerful flares from recoiling black holes in quasars

abstract: Mergers of spinning black holes can give recoil velocities from gravitational radiation up to several thousand km/s. A recoiling supermassive black hole in an AGN retains the inner part of its accretion disk. Marginally bound material rejoining the disk around the moving black hole releases a large amount of energy in shocks in a short time, leading to a flare in thermal soft X-rays with a luminosity approaching the Eddington limit. Reprocessing of the X-rays by the infalling material gives strong optical and ultraviolet emission lines with a distinctive spectrum. Despite the short lifetime of the flare (~10^4 yr), as many as 100 flares may be in play at the present time in QSOs at redshifts ~ 1 to 3. These flares provide a means to identify high velocity recoils.


title: On low rank perturbation of matrices

abstract: The article is devoted to different aspects of the question "What can be done with a matrix by low rank perturbation?" It is proved that one can change a geometrically simple spectrum drastically by a rank 1 permutation, but the situation is quite different if one restricts oneself to normal matrices. Also, the Jordan normal form of a perturbed matrix is considered. It is proved that with respect to rank as a distance all almost unitary matrices are near unitary.


title: Are complex systems hard to evolve?

abstract: Evolutionary complexity is here measured by the number of trials/evaluations needed for evolving a logical gate in a non-linear medium. Behavioural complexity of the gates evolved is characterised in terms of cellular automata behaviour. We speculate that hierarchies of behavioural and evolutionary complexities are isomorphic up to some degree, subject to substrate specificity of evolution and the spectrum of evolution parameters.


title: The evolution of microRNA-regulation in duplicated genes facilitates expression divergence

abstract: Background: The evolution of microRNA regulation in metazoans is a mysterious process: MicroRNA sequences are highly conserved among distal organisms, but on the other hand, there is no evident conservation of their targets. Results: We study this extensive rewiring of the microRNA regulatory network by analyzing the evolutionary trajectories of duplicated genes in D. melanogatser. We find that in general microRNA-targeted genes tend to avoid gene duplication. However, in cases where gene duplication is evident, we find that the gene that displays high divergence from the ancestral gene at the sequence level is also likely to be associated in an opposing manner with the microRNA regulatory system - if the ancestral gene is a miRNA target then the divergent gene tends not to be, and vice versa. Conclusions: This suggests that miRNAs not only have a role in conferring expression robustness, as was suggested by previous works, but are also an accessible tool in evolving expression divergence.


title: Rigorous Derivation of the Gross-Pitaevskii Equation with a Large Interaction Potential

abstract: Consider a system of $N$ bosons in three dimensions interacting via a repulsive short range pair potential $N^2V(N(x_i-x_j))$, where $\bx=(x_1, >..., x_N)$ denotes the positions of the particles. Let $H_N$ denote the Hamiltonian of the system and let $\psi_{N,t}$ be the solution to the Schr\"odinger equation. Suppose that the initial data $\psi_{N,0}$ satisfies the energy condition \[ < \psi_{N,0}, H_N \psi_{N,0} > \leq C N >. \] and that the one-particle density matrix converges to a projection as $N \to \infty$. Then, we prove that the $k$-particle density matrices of $\psi_{N,t}$ factorize in the limit $N \to \infty$. Moreover, the one particle orbital wave function solves the time-dependent Gross-Pitaevskii equation, a cubic non-linear Schr\"odinger equation with the coupling constant proportional to the scattering length of the potential $V$. In \cite{ESY}, we proved the same statement under the condition that the interaction potential $V$ is sufficiently small; in the present work we develop a new approach that requires no restriction on the size of the potential.


title: A Survey of Hydroxyl Toward Supernova Remnants: Evidence for Extended 1720 MHz Maser Emission

abstract: We present the results of GBT observations of all four ground-state hydroxyl (OH) transitions toward 15 supernova remnants (SNRs) which show OH(1720 MHz) maser emission. This species of maser is well established as an excellent tracer of an ongoing interaction between the SNR and dense molecular material. For the majority of these objects we detect significantly higher flux densities with a single dish than has been reported with interferometric observations. We infer that spatially extended, low level maser emission is a common phenomenon that traces the large-scale interaction in maser-emitting SNRs. Additionally we use a collisional pumping model to fit the physical conditions under which OH is excited behind the SNR shock front. We find the observed OH gas associated with the SNR interaction having columns less than approximately 10^17 per square cm, temperatures of 20 to 125 K, and densities 10^5 per cubic cm.


title: Intrinsic and extrinsic noise effects on the phase transition of swarming systems and related network models

abstract: We analyze order-disorder phase transitions driven by noise that occur in two kinds of network models closely related to the self-propelled model proposed by Vicsek et. al. to describe the collective motion of groups of organisms [\emph{Phys. Rev. Lett.} {\bf 75}:1226 (1995)]. Two different types of noise, which we call intrinsic and extrinsic, are considered. The intrinsic noise, the one used by Vicsek et. al. in their original work, is related to the decision mechanism through which the particles update their positions. In contrast, the extrinsic noise, later introduced by Gr\'egoire and Chat\'e [\emph{Phys. Rev. Lett.} {\bf 92}:025702 (2004)], affects the signal that the particles receive from the environment. The network models presented here can be considered as the mean-field representation of the self-propelled model. We show analytically and numerically that, for these two network models, the phase transitions driven by the intrinsic noise are continuous, whereas the extrinsic noise produces discontinuous phase transitions. This is true even for the small-world topology, which induces strong spatial correlations between the network elements. We also analyze the case where both types of noise are present simultaneously. In this situation, the phase transition can be continuous or discontinuous depending upon the amplitude of each type of noise.


title: Doped AB_2 Hubbard Chain: Spiral, Nagaoka and RVB States, Phase Separation and Luttinger Liquid Behavior

abstract: We present an extensive numerical study of the Hubbard model on the doped AB$_2$ chain, both in the weak coupling and the infinite-U limit. Due to the special unit cell topology, this system displays a rich variety of phases as function of hole doping ($\delta$) away from half-filling. Near half-filling, spiral states develop in the weak coupling regime, while Nagaoka itinerant ferromagnetism is observed in the infinite-U limit. For higher doping the system phase-separates before reaching a Mott insulating phase of short-range RVB states at $\delta=1/3$. Moreover, for $\delta>1/3$ we observe a crossover, which anticipates the Luttinger liquid behavior for $\delta > 2/3$.


title: Deriving Sorting Algorithms

abstract: This paper proposes new derivations of three well-known sorting algorithms, in their functional formulation. The approach we use is based on three main ingredients: first, the algorithms are derived from a simpler algorithm, i.e. the specification is already a solution to the problem (in this sense our derivations are program transformations). Secondly, a mixture of inductive and coinductive arguments are used in a uniform, algebraic style in our reasoning. Finally, the approach uses structural invariants so as to strengthen the equational reasoning with logical arguments that cannot be captured in the algebraic framework.


title: Initial Measurements of Black Hole Spin in GX 339-4 from Suzaku Spectroscopy

abstract: We report on a deep Suzaku observation of the stellar-mass black hole GX 339-4 in outburst. A clear, strong, relativistically-shaped iron emission line from the inner accretion disk is observed. The broad-band disk reflection spectrum revealed is one of the most sensitive yet obtained from an accreting black hole. We fit the Suzaku spectra with a physically-motivated disk reflection model, blurred by a new relativistic line function in which the black hole spin parameter is a variable. This procedure yielded a black hole spin parameter of a = 0.89 +/- 0.04. Joint modeling of these Suzaku spectra and prior XMM-Newton spectra obtained in two different outburst phases yields a spin parameter of a = 0.93 +/- 0.01. The degree of consistency between these results suggests that disk reflection models allow for spin measurements that are not strongly biased by scattering effects. We suggest that the best value of the black hole spin parameter is a = 0.93 +/- 0.01 (statistical) +/- 0.04 (systematic). Although preliminary, these results represent the first direct measurement of non-zero spin in a stellar-mass black hole using relativistic line modeling.


title: Genetic selection of neutron star structure matching the X-ray observations

abstract: Assuming a resonant origin of the quasiperiodic oscillations observed in the X-ray neutron star binary systems, we apply a genetic algorithm method for selection of neutron star models. It was suggested that pairs of kilo-Hertz peaks in the X-ray Fourier power density spectra of some neutron stars reflect a non-linear resonance between two modes of accretion disk oscillations. In several specific models, the two modes are related to physically plausible combinations of Keplerian, vertical and radial frequencies of geodesic orbital motion. We investigate this concept for a specific neutron star source, a fixed pair of modes and various neutron star equations of state. Each neutron star model is characterized by the equation of state (EOS), rotation frequency ($\Omega$) and central energy density ($\rho_\mathrm c$). These determine the spacetime structure governing geodesic motion and position dependent radial and vertical epicyclic oscillations related to the stable circular geodesics. When the parameters of neutron star model are fixed, the two considered modes imply a frequency-frequency relation which can be compared to the observation in order to eliminate the unsatisfactory sets (KR,$\rho_\mathrm c, \Omega$, EOS). For the elimination we use the advanced genetic algorithm. Genetic algorithm comes out from the method of natural selection when subjects with the best adaptation to assigned conditions have best chances to survive. The chosen genetic algorithm with sexual reproduction contains one chromosome with restricted lifetime, uniform crossing and genes of type 3/3/5. For encryption of physical description (KR,$\rho_\mathrm c, \Omega$, EOS) into chromosome we use the Gray code. As a fitness function we use correspondence between the observed and calculated pairs of eigenfrequencies.


title: Are Aryabhata's and Galilean Relativity Equivalent?

abstract: This paper considers the question of whether Aryabhata's direct and indirect references to relativity of motion are equivalent to Galilean Relativity. Examining verses fom different sections of the text Aryabhatiya, we show that there is explicit mention of relativity of space and motion, although there is no explicit reference to the motion being uniform.


title: Rich, Sturmian, and trapezoidal words

abstract: In this paper we explore various interconnections between rich words, Sturmian words, and trapezoidal words. Rich words, first introduced in arXiv:0801.1656 by the second and third authors together with J. Justin and S. Widmer, constitute a new class of finite and infinite words characterized by having the maximal number of palindromic factors. Every finite Sturmian word is rich, but not conversely. Trapezoidal words were first introduced by the first author in studying the behavior of the subword complexity of finite Sturmian words. Unfortunately this property does not characterize finite Sturmian words. In this note we show that the only trapezoidal palindromes are Sturmian. More generally we show that Sturmian palindromes can be characterized either in terms of their subword complexity (the trapezoidal property) or in terms of their palindromic complexity. We also obtain a similar characterization of rich palindromes in terms of a relation between palindromic complexity and subword complexity.


title: Galaxy Evolution in Dense Environments: Properties of Interacting Galaxies in the A901/02 Supercluster

abstract: We present a study of galaxies in the STAGES survey of the Abell 901/902 supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps. We characterize galaxies with strong, externally-triggered morphological distortions and normal, relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We report here our preliminary results.


title: Search for Pair Production of Scalar Top Quarks Decaying to a tau Lepton and a b Quark in ppbar Collisions at sqrt{s}=1.96 TeV

abstract: We search for pair production of supersymmetric top quarks (~t_1), followed by R-parity violating decay ~t_1 -> tau b with a branching ratio beta, using 322 pb^-1 of ppbar collisions at sqrt{s}=1.96 TeV collected by the CDF II detector at Fermilab. Two candidate events pass our final selection criteria, consistent with the standard model expectation. We set upper limits on the cross section sigma(~t_1 ~tbar_1)*beta^2 as a function of the stop mass m(~t_1). Assuming beta=1, we set a 95% confidence level limit m(~t_1)>153 GeV/c^2. The limits are also applicable to the case of a third generation scalar leptoquark (LQ_3) decaying LQ_3 -> tau b.


title: Directive words of episturmian words: equivalences and normalization

abstract: Episturmian morphisms constitute a powerful tool to study episturmian words. Indeed, any episturmian word can be infinitely decomposed over the set of pure episturmian morphisms. Thus, an episturmian word can be defined by one of its morphic decompositions or, equivalently, by a certain directive word. Here we characterize pairs of words directing a common episturmian word. We also propose a way to uniquely define any episturmian word through a normalization of its directive words. As a consequence of these results, we characterize episturmian words having a unique directive word.


title: Polygones de Newton de certaines sommes de caract\`eres et s\'eries de Poincar\'e

abstract: In this paper, we shall precise the asymptotic behaviour of Newton polygons of $L$ functions associated to character sums, coming from some $n$ variable Laurent polynomials. In order to do this, we use the free sum on convex polytopes. This operation allows the determination of the limit of generic Newton polygons for the sum $\Delta=\Delta_1\oplus \Delta_2$ when we know the limit of generic Newton polygons for each factor. To our knowledge, these are the first results concerning the asymptotic behaviour of Newton polygons for multivariable polynomials.


title: PGA Tour Scores as a Gaussian Random Variable

abstract: In this paper it is demonstrated that the scoring at each PGA Tour stroke play event can be reasonably modeled as a Gaussian random variable. All 46 stroke play events in the 2007 season are analyzed. The distributions of scores are favorably compared with a Gaussian distribution using the Kolmogorov-Smirnov test. This observation suggests performance tracking on the PGA tour should be done in terms of the z-score, calculated by subtracting the mean from the raw score and dividing by the standard deviation. This methodology measures performance relative to the field of competitors, independent of the venue, and in terms of a statistic that has quantitative meaning. Several examples of the use of this scoring methodology are provided, including a calculation of the probability that Tiger Woods will break Byron Nelson's record of eleven consecutive PGA Tour victories.


title: On logical gates in precipitating medium: cellular automaton model

abstract: We study a two-dimensional semi-totalistic binary cell-state cellular automaton, which imitates a reversible precipitation in an abstract chemical medium. The systems exhibits a non-trivial growth and nucleation. We demonstrate how basic computational operation can be realized in the system when the propagation of the growing patterns is self-restricted by stationary localizations. We show that precipitating patterns of different morphology compete between each other and thus implement serial and non-serial logical gates.


title: Incomplete quantum process tomography and principle of maximal entropy

abstract: The main goal of this paper is to extend and apply the principle of maximum entropy (MaxEnt) to incomplete quantum process estimation tasks. We will define a so-called process entropy function being the von Neumann entropy of the state associated with the quantum process via Choi-Jamiolkowski isomorphism. It will be shown that an arbitrary process estimation experiment can be reformulated in a unified framework and MaxEnt principle can be consistently exploited. We will argue that the suggested choice for the process entropy satisfies natural list of properties and it reduces to the state MaxEnt principle, if applied to preparator devices.


title: The cross-correlation of the CMB polarisation and the 21-cm line fluctuations from cosmic reionisation

abstract: The cosmic microwave background (CMB) polarisation and the 21 cm line fluctuations are powerful probes of cosmological reionisation. We study how the cross-correlation between the CMB polarisation (E-modes) and the 21 cm line fluctuations can be used to gain further understanding of the reionisation history, within the framework of inhomogeneous reionisation. Since the E-mode polarisation reflects the amplitude of the quadrupole component of the CMB temperature fluctuations, the angular power spectrum of the cross-correlation exhibits oscillations at all multipoles. The first peak of the power spectrum appears at the scale corresponding to the quadrupole at the redshift that is probed by the 21 cm line fluctuations. The peak reaches its maximum value in redshift when the average ionisation fraction of the universe is about half. On the other hand, on small scales, there is a damping that depends on the duration of reionisation. Thus, the cross-correlation between the CMB polarisation and the 21 cm line fluctuations has the potential to constrain accurately the epoch and the duration of reionisation.


title: Pattern formation in mixtures of ultracold atoms in optical lattices

abstract: Regular pattern formation is ubiquitous in nature; it occurs in biological, physical, and materials science systems. Here we propose a set of experiments with ultracold atoms that show how to examine different types of pattern formation. In particular, we show how one can see the analog of labyrinthine patterns (so-called quantum emulsions) in mixtures of light and heavy atoms (that tend to phase separate) by tuning the trap potential and we show how complex geometrically ordered patterns emerge (when the mixtures do not phase separate), which could be employed for low-temperature thermometry. The complex physical mechanisms for the pattern formation at zero temperature are understood within a theoretical analysis called the local density approximation.


title: Complexity Metrics for Spreadsheet Models

abstract: Several complexity metrics are described which are related to logic structure, data structure and size of spreadsheet models. They primarily concentrate on the dispersion of cell references and cell paths. Most metrics are newly defined, while some are adapted from traditional software engineering. Their purpose is the identification of cells which are liable to errors. In addition, they can be used to estimate the values of dependent process metrics, such as the development duration and effort, and especially to adjust the cell error rate in accordance with the contents of each individual cell, in order to accurately asses the reliability of a model. Finally, two conceptual constructs - the reference branching condition cell and the condition block - are discussed, aiming at improving the reliability, modifiability, auditability and comprehensibility of logical tests.


title: Optimal tracking for pairs of qubit states

abstract: In classical control theory, tracking refers to the ability to perform measurements and feedback on a classical system in order to enforce some desired dynamics. In this paper we investigate a simple version of quantum tracking, namely, we look at how to optimally transform the state of a single qubit into a given target state, when the system can be prepared in two different ways, and the target state depends on the choice of preparation. We propose a tracking strategy that is proved to be optimal for any input and target states. Applications in the context of state discrimination, state purification, state stabilization and state-dependent quantum cloning are presented, where existing optimality results are recovered and extended.


title: Consistency test of general relativity from large scale structure of the Universe

abstract: We construct a consistency test of General Relativity (GR) on cosmological scales. This test enables us to distinguish between the two alternatives to explain the late-time accelerated expansion of the universe, that is, dark energy models based on GR and modified gravity models without dark energy. We derive the consistency relation in GR which is written only in terms of observables - the Hubble parameter, the density perturbations, the peculiar velocities and the lensing potential. The breakdown of this consistency relation implies that the Newton constant which governs large-scale structure is different from that in the background cosmology, which is a typical feature in modified gravity models. We propose a method to perform this test by reconstructing the weak lensing spectrum from measured density perturbations and peculiar velocities. This reconstruction relies on Poisson's equation in GR to convert the density perturbations to the lensing potential. Hence any inconsistency between the reconstructed lensing spectrum and the measured lensing spectrum indicates the failure of GR on cosmological scales. The difficulties in performing this test using actual observations are discussed.


title: Constraining Inflation

abstract: Slow roll reconstruction is derived from the Hamilton-Jacobi formulation of inflationary dynamics. It automatically includes information from sub-leading terms in slow roll, and facilitatesthe inclusion of priors based on the duration on inflation. We show that at low inflationary scales the Hamilton-Jacobi equations simplify considerably. We provide a new classification scheme for inflationary models, based solely on the number of parameters needed to specify the potential, and provide forecasts for likely bounds on the slow roll parameters from future datasets. A minimal running of the spectral index, induced solely by the first two slow roll parameters (\epsilon and \eta) appears to be effectively undetectable by realistic Cosmic Microwave Background experiments. However, we show that the ability to detect this signal increases with the lever arm in comoving wavenumber, and we conjecture that high redshift 21 cm data may allow tests of second order consistency conditions on inflation. Finally, we point out that the second order corrections to the spectral index are correlated with the inflationary scale, and thus the amplitude of the CMB B-mode.


title: The orbital poles of Milky Way satellite galaxies: a rotationally supported disc-of-satellites

abstract: Available proper motion measurements of Milky Way (MW) satellite galaxies are used to calculate their orbital poles and projected uncertainties. These are compared to a set of recent cold dark-matter (CDM) simulations, tailored specifically to solve the MW satellite problem. We show that the CDM satellite orbital poles are fully consistent with being drawn from a random distribution, while the MW satellite orbital poles indicate that the disc-of-satellites of the Milky Way is rotationally supported. Furthermore, the bootstrapping analysis of the spatial distribution of theoretical CDM satellites also shows that they are consistent with being randomly drawn. The theoretical CDM satellite population thus shows a significantly different orbital and spatial distribution than the MW satellites, most probably indicating that the majority of the latter are of tidal origin rather than being DM dominated sub-structures. A statistic is presented that can be used to test a possible correlation of satellite galaxy orbits with their spatial distribution.


title: Infinite randomness fixed point of the superconductor-metal quantum phase transition

abstract: We examine the influence of quenched disorder on the superconductor-metal transition, as described by a theory of overdamped Cooper pairs which repel each other. The self-consistent pairing eigenmodes of a quasi-one dimensional wire are determined numerically. Our results support the recent proposal by Hoyos et al. (arXiv:0705.1865) that the transition is characterized by the same strong disorder fixed point describing the onset of ferromagnetism in the random quantum Ising chain in a transverse field.


title: Frequency and Impact of Galaxy Mergers and Interactions over the last 7 Gyr

abstract: We explore the history and impact of galaxy mergers and interactions over z~0.24 to 0.80, based on HST ACS, Combo-17, and Spitzer 24 mu data of ~4500 galaxies in the GEMS survey. Using visual and quantitative parameters,we identify galaxies with strong distortions indicative of recent strong interactions and mergers versus normal galaxies (E/S0, Sa, Sb-Sc, Sd/Irr). Our results are: (1) The observed fraction F of strongly disturbed systems among high mass (M>=2.5E10 Msun) galaxies is ~9% to 12% in every Gyr bin over z~0.24 to 0.80. The corresponding merger rate is a few times 10^-4 galaxies Gyr-1 Mpc-3. The fraction F shows fair agreement with the merger fraction of mass ratio >=1:10 predicted by several LCDM-based simulations. (2) For M>=1E9 Msun systems, the average star formation rate (SFR) of strongly disturbed systems is only modestly enhanced with respect to normal galaxies, in agreement with recent simulations. In fact, over z~0.24 to 0.80, strongly disturbed systems only account for a small fraction (<~30%) of the total SFR density. This suggests that the behaviour of the cosmic SFR density over the last seven billion years is predominantly shaped by normal galaxies.


title: Electronic and optical properties of electromigrated molecular junctions

abstract: Electromigrated nanoscale junctions have proven very useful for studying electronic transport at the single-molecule scale. However, confirming that conduction is through precisely the molecule of interest and not some contaminant or metal nanoparticle has remained a persistent challenge, typically requiring a statistical analysis of many devices. We review how transport mechanisms in both purely electronic and optical measurements can be used to infer information about the nanoscale junction configuration. The electronic response to optical excitation is particularly revealing. We briefly discuss surface-enhanced Raman spectroscopy on such junctions, and present new results showing that currents due to optical rectification can provide a means of estimating the local electric field at the junction due to illumination.


title: Constraining Galaxy Evolution With Bulge-Disk-Bar Decomposition

abstract: Structural decomposition of galaxies into bulge, disk, and bar components is important to address a number of scientific problems. Measuring bulge, disk, and bar structural parameters will set constraints on the violent and secular processes of galaxy assembly and recurrent bar formation and dissolution models. It can also help to quantify the fraction and properties of bulgeless galaxies (those systems having no bulge or only a relatively insignificant disky-pseudobulges), which defy galaxy formation paradigms requiring almost every disk galaxy to have a classical bulge at its core. We demonstrate a proof of concept and show early results of our ongoing three-component bulge-disk-bar decomposition of NIR images for a sample of three complementary samples spanning different epochs and different environments (field and cluster). In contrast to most early studies, which only attempt two-component bulge-disk decomposition, we fit three components using GALFIT: a bulge, a disk, and a bar. We show that it is important to include the bar component, as this can significantly lower the bulge-to-total luminosity ratio (B/T), in many cases by a factor of two or more, thus effectively changing the Hubble type of a galaxy from early to late.


title: The mass-metallicity relation of interacting galaxies

abstract: We study the mass-metallicity relation of galaxies in pairs and in isolation taken from the SDSS-DR4 using the stellar masses and oxygen abundances derived by Tremonti et al. (2004). Close galaxy pairs, defined by projected separation r_p < 25kpc/h and radial velocity Delta_V < 350 km/s, are morphologically classified according to the strength of the interaction signs. We find that only for pairs showing signs of strong interactions, the mass-metallicity relation differs significantly from that of galaxies in isolation. In such pairs, the mean gas-phase oxygen abundances of galaxies with low stellar masses (Mstar ~< 10^9 Msun/h) exhibit an excess of 0.2 dex. Conversely, at larger masses (Mstar >~ 10^10 Msun/h) galaxies have a systematically lower metallicity, although with a smaller difference (-0.05 dex). Similar trends are obtained if g-band magnitudes are used instead of stellar masses. In minor interactions, we find that the less massive member is systematically enriched, while a galaxy in interaction with a comparable stellar mass companion shows a metallicity decrement with respect to galaxies in isolation. We argue that metal-rich starbursts triggered by a more massive component, and inflows of low metallicity gas induced by comparable or less massive companion galaxies, provide a natural scenario to explain our findings.


title: A note on Gribov copies in 3D Chern-Simons theory

abstract: Using powerful tools of harmonic maps and integrable systems, all the Gribov copies in the Coulomb gauge in 3D Chern-Simons theory are constructed. Some issues about the Gribov and the modular re- gions are shortly discussed. The Gribov copies of the vacuum in 3D QCD in the Coulomb gauge are described. An interesting implication of the presence of Gribov copies is briefy pointed out.


title: Flavor Physics in SO(10) GUTs with Suppressed Proton decay Due to Gauged Discrete Symmetry

abstract: Generic SO(10) GUT models suffer from the problem that Planck scale induced non-renormalizable proton decay operators require extreme suppression of their couplings to be compatible with present experimental upper limits. One way to resolve this problem is to supplement SO(10) by simple gauged discrete symmetries which can also simultaneously suppress the renormalizable R-parity violating ones when they occur and make the theory "more natural". Here we discuss the phenomenological viability of such models. We first show that for both classes of models, e.g the ones that use ${\bf 16}_H$ or ${\bf 126}_H$ to break B-L symmetry, the minimal Higgs content which is sufficient for proton decay suppression is inadequate for explaining fermion masses despite the presence of all apparently needed couplings. We then present an extended ${\bf 16}_H$ model, with three {\bf 10} and three {\bf 45}-Higgs, where is free of this problem. We propose this as a realistic and "natural" model for fermion unification and discuss the phenomenology of this model e.g. its predictions for neutrino mixings and lepton flavor violation.


title: Fast beam stacking using RF barriers

abstract: Two barrier RF systems were fabricated, tested and installed in the Fermilab Main Injector. Each can provide 8 kV rectangular pulses (the RF barriers) at 90 kHz. When a stationary barrier is combined with a moving barrier, injected beams from the Booster can be continuously deflected, folded and stacked in the Main Injector, which leads to doubling of the beam intensity. This paper gives a report on the beam experiment using this novel technology.


title: Morphological Transformations of Galaxies in the A901/02 Supercluster from STAGES

abstract: We present a study of galaxies in the Abell 901/902 Supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps, as part of the STAGES survey. We characterize galaxies with strong externally-triggered morphological distortions and normal relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We revisit the morphology density relation, which postulates a higher fraction of early type galaxies in dense environments, by considering separately galaxies with a low bulge-to-disk (B/D) ratio and a low gas content as these two parameters may not be correlated in clusters. We report here on our preliminary analysis.


title: Update of eta-eta' mixing from J/psi->VP decays

abstract: The $\eta$-$\eta^\prime$ mixing angle is deduced from an updated phenomenological analysis of $J/\psi$ decays into a vector and a pseudoscalar meson. Corrections due to non-ideal $\omega$-$\phi$ mixing are confirmed to be crucial to find $\theta_P=(-16.9\pm 1.7)^\circ$, in agreement with most recent analyses. The newly reported values of $B(J/\psi\to\rho\pi)$ by the BABAR and BES Collaborations make more difficult a reasonable description of data.


title: Characterizing Barred Galaxies in the Abell 901/902 Supercluster

abstract: In dense clusters, higher densities at early epochs as well as physical processes, such as ram pressure stripping and tidal interactions become important, and can have direct consequences for the evolution of bars and their host disks. To study bars and disks as a function of environment, we are using the STAGES ACS HST survey of the Abell 901/902 supercluster (z~0.165), along with earlier field studies based the SDSS and the Ohio State University Bright Spiral Galaxy Survey (OSUBSGS). We explore the limitations of traditional methods for characterizing the bar fraction, and in particular highlight uncertainties in disk galaxy selection in cluster environments. We present an alternative approach for exploring the proportion of bars, and investigate the properties of bars as a function of host galaxy color, Sersic index, stellar mass, star formation rate (SFR), specific SFR, and morphology.


title: Operational aspects of the Main Injector large aperture quadrupole (WQB)

abstract: A two-year Large Aperture Quadrupole (WQB) Project was completed in the summer of 2006 at Fermilab. Nine WQBs were designed, fabricated and bench-tested by the Technical Division. Seven of them were installed in the Main Injector and the other two for spares. They perform well. The aperture increase meets the design goal and the perturbation to the lattice is minimal. The machine acceptance in the injection and extraction regions is increased from 40pi to 60pi mm-mrad. This paper gives a brief report of the operation and performance of these magnets.


title: Alignment Timescale of the Microquasar GRO J1655-40

abstract: The microquasar GRO J1655-40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady state disc warped by Lense-Thirring precession and find the timescale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in the latter case but we find that, in both cases, the lifetime of the mass transfer state is at most a few times the alignment timescale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter-aligned after its supernova kick but that it is most likely to be close to alignment rather than counteralignment now.


title: The Scattered Disk as the source of the Jupiter Family comets

abstract: The short period Jupiter family comets (JFCs) are thought to originate in the Kuiper Belt; specifically, a dynamical subclass of the Kuiper Belt known as the `scattered disk' is argued to be the dominant source of JFCs. However, the best estimates from observational surveys indicate that this source may fall short by more than two orders of magnitude the estimates obtained from theoretical models of the dynamical evolution of Kuiper belt objects into JFCs. We re-examine the scattered disk as a source of the JFCs and make a rigorous estimate of the discrepancy. We find that the uncertainties in the dynamical models combined with a change in the size distribution function of the scattered disk at faint magnitudes (small sizes) beyond the current observational limit offer a possible but problematic resolution to the discrepancy. We discuss several other possibilities: that the present population of JFCs is a large fluctuation above their long term average, that larger scattered disk objects tidally break-up into multiple fragments during close planetary encounters as their orbits evolve from the trans-Neptune zone to near Jupiter, or that there are alternative source populations that contribute significantly to the JFCs. Well-characterized observational investigations of the Centaurs, objects that are transitioning between the trans-Neptune Kuiper belt region and the inner solar system, can test the predictions of the non-steady state and the tidal break-up hypotheses. The classical and resonant classes of the Kuiper belt are worth re-consideration as significant additional or alternate sources of the JFCs.


title: Bimodality in Damped Lyman alpha Systems

abstract: We report evidence for a bimodality in damped Ly systems (DLAs). Using [C II] 158 mu cooling rates, lc, we find a distribution with peaks at lc=10^-27.4 and 10^-26.6 ergs s^-1 H^-1 separated by a trough at lc^crit ~= lc < 10^-27.0 ergs s^-1 H^-1. We divide the sample into low cool DLAs with lc < lc^crit and high cool DLAs with lc > lc^crit and find the Kolmogorv-Smirnov probabilities that velocity width, metallicity, dust-to-gas ratio, and Si II equivalent width in the two subsamples are drawn from the same parent population are small. All these quantities are significantly larger in the high cool population, while the H I column densities are indistinguishable in the two populations. We find that heating by X-ray and FUV background radiation is insufficient to balance the cooling rates of either population. Rather, the DLA gas is heated by local radiation fields. The rare appearance of faint, extended objects in the Hubble Ultra Deep Field rules out in situ star formation as the dominant star-formation mode for the high cool population, but is compatible with in situ star formation as the dominant mode for the low cool population. Star formation in the high cool DLAs likely arises in Lyman Break galaxies. We investigate whether these properties of DLAs are analogous to the bimodal properties of nearby galaxies. Using Si II equivalent width as a mass indicator, we construct bivariate distributions of metallicity, lc, and areal SFR versus the mass indicators. Tentative evidence is found for correlations and parallel sequences, which suggest similarities between DLAs and nearby galaxies. We suggest that the transition-mass model provides a plausible scenario for the bimodality we have found. As a result, the bimodality in current galaxies may have originated in DLAs.


title: GoArrays: highly dynamic and efficient microarray probe design

abstract: MOTIVATION: The use of oligonucleotide microarray technology requires a very detailed attention to the design of specific probes spotted on the solid phase. These problems are far from being commonplace since they refer to complex physicochemical constraints. Whereas there are more and more publicly available programs for microarray oligonucleotide design, most of them use the same algorithm or criteria to design oligos, with only little variation. RESULTS: We show that classical approaches used in oligo design software may be inefficient under certain experimental conditions, especially when dealing with complex target mixtures. Indeed, our biological model is a human obligate parasite, the microsporidia Encephalitozoon cuniculi. Targets that are extracted from biological samples are composed of a mixture of pathogen transcripts and host cell transcripts. We propose a new approach to design oligonucleotides which combines good specificity with a potentially high sensitivity. This approach is original in the biological point of view as well as in the algorithmic point of view. We also present an experimental validation of this new strategy by comparing results obtained with standard oligos and with our composite oligos. A specific E.cuniculi microarray will overcome the difficulty to discriminate the parasite mRNAs from the host cell mRNAs demonstrating the power of the microarray approach to elucidate the lifestyle of an intracellular pathogen using mix mRNAs.


title: A fast Bayesian approach to discrete object detection in astronomical datasets - PowellSnakes I

abstract: A new fast Bayesian approach is introduced for the detection of discrete objects immersed in a diffuse background. This new method, called PowellSnakes, speeds up traditional Bayesian techniques by: i) replacing the standard form of the likelihood for the parameters characterizing the discrete objects by an alternative exact form that is much quicker to evaluate; ii) using a simultaneous multiple minimization code based on Powell's direction set algorithm to locate rapidly the local maxima in the posterior; and iii) deciding whether each located posterior peak corresponds to a real object by performing a Bayesian model selection using an approximate evidence value based on a local Gaussian approximation to the peak. The construction of this Gaussian approximation also provides the covariance matrix of the uncertainties in the derived parameter values for the object in question. This new approach provides a speed up in performance by a factor of `hundreds' as compared to existing Bayesian source extraction methods that use MCMC to explore the parameter space, such as that presented by Hobson & McLachlan. We illustrate the capabilities of the method by applying to some simplified toy models. Furthermore PowellSnakes has the advantage of consistently defining the threshold for acceptance/rejection based on priors which cannot be said of the frequentist methods. We present here the first implementation of this technique (Version-I). Further improvements to this implementation are currently under investigation and will be published shortly. The application of the method to realistic simulated Planck observations will be presented in a forthcoming publication.


title: Exploring the Impact of Galaxy Interactions over Seven Billion Years with CAS

abstract: We explore galaxy assembly over the last seven billion years by characterizing "normal" galaxies along the Hubble sequence, against strongly disturbed merging/interacting galaxies with the widely used CAS system of concentration (C), asymmetry (A), and 'clumpiness' (S) parameters, as well as visual classification. We analyze Hubble Space Telescope (HST) ACS images of ~4000 intermediate and high mass (> 10^9 solar masses) galaxies from the GEMS survey, one of the largest HST surveys conducted to date in two filters. We explore the effectiveness of the CAS criteria [A>S and A>~0.35] in separating normal and strongly disturbed galaxies at different redshifts, and quantify the recovery and contamination rate. We also compare the average star formation rate and the cosmic star formation rate density as a function of redshift between normal and interacting systems identified by CAS.


title: The Impact of HD Cooling on the Formation of the First Stars

abstract: We use numerical simulations to investigate the importance of HD formation and cooling on the first generation of metal-free stars in a LCDM cosmology. We have implemented and tested non-equilibrium HD chemistry in an adaptive mesh refinement simulation code and applied it to two situations. (1) It is first applied to the formation of 10^5 - 10^6 Msun halos which form in the absence of any ionizing source. We show, in agreement with previous work, that HD cooling is of only marginal importance for most halos; however, we find that for the lowest mass halos, HD cooling can equal or surpass the H2 cooling rate. This leads to a population of stars formed in halos with effective HD cooling that are less massive by a factor of ~6 compared to halos dominated by H2 cooling. (2) In the second part of the paper, we ionize the halos in order to explore the impact of HD cooling in the presence of an ample population of free electrons. This leads to cooler temperatures (due to the electron- catalyzed production of H2) and somewhat lower resulting proto-stellar mass. Adding HD chemistry lowers the temperature further, to the level of the CMB. We find that HD cooling dominates over H2 cooling in the density range 10^2 cm^-3 to 10^6 cm^-3, but above this density, the temperature rises and H2 cooling dominates again. Because of this, the accretion rate on to the protostar is almost the same as in the H2 case (at least for accreted masses below 50-100 Msun), therefore we argue that HD cooling in ionized halos will probably not result in a population of significantly lower mass stars.


title: A Paradigm for Spreadsheet Engineering Methodologies

abstract: Spreadsheet engineering methodologies are diverse and sometimes contradictory. It is difficult for spreadsheet developers to identify a spreadsheet engineering methodology that is appropriate for their class of spreadsheet, with its unique combination of goals, type of problem, and available time and resources. There is a lack of well-organized, proven methodologies with known costs and benefits for well-defined spreadsheet classes. It is difficult to compare and critically evaluate methodologies. We present a paradigm for organizing and interpreting spreadsheet engineering recommendations. It systematically addresses the myriad choices made when developing a spreadsheet, and explicitly considers resource constraints and other development parameters. This paradigm provides a framework for evaluation, comparison, and selection of methodologies, and a list of essential elements for developers or codifiers of new methodologies. This paradigm identifies gaps in our knowledge that merit further research.


title: Acceleration radiation and the Planck scale

abstract: A uniformly accelerating observer perceives the Minkowski vacuum state as a thermal bath of radiation. We point out that this field-theory effect can be derived, for any dimension higher than two, without actually invoking very high energy physics. This supports the view that this phenomenon is robust against Planck-scale physics and, therefore, should be compatible with any underlying microscopic theory.


title: The commutants of certain Toeplitz operators on weighted Bergman spaces

abstract: For $\alpha>-1$, let $A^2_{\alpha}$ be the corresponding weighted Bergman space of the unit ball in $\mathbb{C}^n$. For a bounded measurable function $f$, let $T_f$ be the Toeplitz operator with symbol $f$ on $A^2_{\alpha}$. This paper describes all the functions $f$ for which $T_f$ commutes with a given $T_g$, where $g(z)=z_{1}^{L_1}... z_{n}^{L_n}$ for strictly positive integers $L_1,..., L_n$, or $g(z)=|z_1|^{s_1}... |z_n|^{s_n}h(|z|)$ for non-negative real numbers $s_1,..., s_n$ and a bounded measurable function $h$ on $[0,1)$.


title: Constrained Consensus

abstract: We present distributed algorithms that can be used by multiple agents to align their estimates with a particular value over a network with time-varying connectivity. Our framework is general in that this value can represent a consensus value among multiple agents or an optimal solution of an optimization problem, where the global objective function is a combination of local agent objective functions. Our main focus is on constrained problems where the estimate of each agent is restricted to lie in a different constraint set. To highlight the effects of constraints, we first consider a constrained consensus problem and present a distributed ``projected consensus algorithm'' in which agents combine their local averaging operation with projection on their individual constraint sets. This algorithm can be viewed as a version of an alternating projection method with weights that are varying over time and across agents. We establish convergence and convergence rate results for the projected consensus algorithm. We next study a constrained optimization problem for optimizing the sum of local objective functions of the agents subject to the intersection of their local constraint sets. We present a distributed ``projected subgradient algorithm'' which involves each agent performing a local averaging operation, taking a subgradient step to minimize its own objective function, and projecting on its constraint set. We show that, with an appropriately selected stepsize rule, the agent estimates generated by this algorithm converge to the same optimal solution for the cases when the weights are constant and equal, and when the weights are time-varying but all agents have the same constraint set.


title: On Super-Planckian Fields at Sub-Planckian Energies

abstract: For a light scalar coupled to gravity, I study the gravitational backreaction associated with large field variations. I show a generic obstruction in sourcing a super-Planckian scalar profile without making the whole experiment collapse into a black hole. In empty space the scalar variation obeys an absolute bound, of order of the Planck scale. A Newtonian analysis suggests that inside its sources the scalar can undergo arbitrarily large variations without causing large gravitational backreactions. However the maximum attainable \Delta \phi increases only logarithmically with the size of the source. The bound straightforwardly generalizes to any number of dimensions, and to moduli space-like cases, where it applies to the invariant length in field space as measured by the kinetic metric.


title: A Toolkit for Scalable Spreadsheet Visualization

abstract: This paper presents a toolkit for spreadsheet visualization based on logical areas, semantic classes and data modules. Logical areas, semantic classes and data modules are abstract representations of spreadsheet programs that are meant to reduce the auditing and comprehension effort, especially for large and regular spreadsheets. The toolkit is integrated as a plug-in in the Gnumeric spreadsheet system for Linux. It can process large, industry scale spreadsheet programs in reasonable time and is tightly integrated with its host spreadsheet system. Users can generate hierarchical and graph-based representations of their spreadsheets. This allows them to spot conceptual similarities in different regions of the spreadsheet, that would otherwise not fit on a screen. As it is assumed that the learning effort for effective use of such a tool should be kept low, we aim for intuitive handling of most of the tool's functions.


title: A refined Luecking's theorem and finite-rank products of Toeplitz operators

abstract: For any function $f$ in $L^{\infty}(\mathbb{D})$, let $T_f$ denote the corresponding Toeplitz operator the Bergman space $A^2(\mathbb{D})$. A recent result of D. Luecking shows that if $T_f$ has finite rank then $f$ must be the zero function. Using a refined version of this result, we show that if all except possibly one of the functions $f_1,..., f_{m}$ are radial and $T_{f_1}... T_{f_m}$ has finite rank, then one of these functions must be zero.


title: Stochastic Network Model of Receptor Cross-Talk Predicts Anti-Angiogenic Effects

abstract: Cancer invasion and metastasis depend on angiogenesis. The cellular processes (growth, migration, and apoptosis) that occur during angiogenesis are tightly regulated by signaling molecules. Thus, understanding how cells synthesize multiple biochemical signals initiated by key external stimuli can lead to the development of novel therapeutic strategies to combat cancer. In the face of large amounts of disjoint experimental data generated from multitudes of laboratories using various assays, theoretical signal transduction models provide a framework to distill this vast amount of data. Such models offer an opportunity to formulate and test new hypotheses, and can be used to make experimentally verifiable predictions. This study is the first to propose a network model that highlights the cross-talk between the key receptors involved in angiogenesis, namely growth factor, integrin, and cadherin receptors. From available experimental data, we construct a stochastic Boolean network model of receptor cross-talk, and systematically analyze the dynamical stability of the network under continuous-time Boolean dynamics with a noisy production function. We find that the signal transduction network exhibits a robust and fast response to external signals, independent of the internal cell state. We derive an input-output table that maps external stimuli to cell phenotypes, which is extraordinarily stable against molecular noise with one important exception: an oscillatory feedback loop between the key signaling molecules RhoA and Rac1 is unstable under arbitrarily low noise, leading to erratic, dysfunctional cell motion. Finally, we show that the network exhibits an apoptotic response rate that increases with noise, suggesting that the probability of programmed cell death depends on cell health.

categories: q-bio.MN cond-mat.dis-nn nlin.CG q-bio.QM


title: Kinematics in Matrix Gravity

abstract: We develop the kinematics in Matrix Gravity, which is a modified theory of gravity obtained by a non-commutative deformation of General Relativity. In this model the usual interpretation of gravity as Riemannian geometry is replaced by a new kind of geometry, which is equivalent to a collection of Finsler geometries with several Finsler metrics depending both on the position and on the velocity. As a result the Riemannian geodesic flow is replaced by a collection of Finsler flows. This naturally leads to a model in which a particle is described by several mass parameters. If these mass parameters are different then the equivalence principle is violated. In the non-relativistic limit this also leads to corrections to the Newton's gravitational potential. We find the first and second order corrections to the usual Riemannian geodesic flow and evaluate the anomalous nongeodesic acceleration in a particular case of static spherically symmetric background.


title: Non-Markovian Dynamics of Charge Carriers in Quantum Dots

abstract: We have investigated the dynamics of bound particles in multilevel current-carrying quantum dots. We look specifically in the regime of resonant tunnelling transport, where several channels are available for transport. Through a non-Markovian formalism under the Born approximation, we investigate the real-time evolution of the confined particles including transport-induced decoherence and relaxation. In the case of a coherent superposition between states with different particle number, we find that a Fock-space coherence may be preserved even in the presence of tunneling into and out of the dot. Real-time results are presented for various asymmetries of tunneling rates into different orbitals.


title: Quark-nova remnants III:Application to RRATs

abstract: This is the third paper of a series of papers where we explore the evolution of iron-rich ejecta from quark-novae. In the first paper, we explored the case where a quark-nova ejecta forms a degenerate shell, supported by the star's magnetic field, with applications to SGRs. In the second paper we considered the case where the ejecta would have sufficient angular momentum to form a degenerate Keplerian torus and applied such a system to two AXPs, namely 1E2259 + 586 and 4U0142 + 615. Here we explore the late evolution of the degenerate torus and find that it can remain unchanged for $\sim 10^{6}$ years before it becomes non-degenerate. This transition from a degenerate torus (accretion dominated) to a non-degenerate disk (no accretion), occurs about $10^{6}$ years following the quark-nova, and exhibits features that are reminiscent of observed properties of RRATs. Using this model, we can account for the duration of both the radio bursts and the quiet phase, as well as the observed radio flux from RRATs. The unique on and off activity of the radio pulsar PSR B1931 + 24 is similar to that of "old RRATs" in our model. For old RRATs, in our model, the spin-down rate during the radio loud phase is about a factor 1.6 larger than when it is quiet, remarkably similar to what has been measured for PSR B1931 + 24. We discuss a connection between XDINs and RRATs and argue that some XDINs may be ``dead RRATs'' that have already consumed their non-degenerate disk.


title: Growth Gap vs. smoothness for diffeomorphisms of the interval

abstract: Given a diffeomorphism of the interval, consider the uniform norm of the derivative of its n-th iteration. We get a sequence of real numbers called the growth sequence. Its asymptotic behavior is an invariant which naturally appears both in smooth dynamics and in geometry of the diffeomorphisms groups. We find sharp estimates for the growth sequence of a given diffeomorphism in terms of the modulus of continuity of its derivative. These estimates extend previous results of Polterovich and Sodin, and Borichev.


title: The Energetics of a Flaring Solar Active Region, and Observed Flare Statistics

abstract: A stochastic model for the energy of a flaring solar active region is presented, generalising and extending the approach of Wheatland & Glukhov (1998). The probability distribution for the free energy of an active region is described by the solution to a master equation involving deterministic energy input and random jump transitions downwards in energy (solar flares). It is shown how two observable distributions, the flare frequency-energy distribution and the flare waiting-time distribution, may be derived from the steady-state solution to the master equation, for given choices for the energy input and for the rates of flare transitions. An efficient method of numerical solution of the steady-state master equation is presented. Solutions appropriate for flaring, involving a constant rate of energy input and power-law distributed jump transition rates, are numerically investigated. The flare-like solutions exhibit power-law flare frequency-energy distributions below a high energy rollover, set by the largest energy the active region is likely to have. The solutions also exhibit approximately exponential (i.e. Poisson) waiting-time distributions, despite the rate of flaring depending on the free energy of the system.


title: The lipid bilayer at the mesoscale: a physical continuum model

abstract: We study a continuum model of the lipid bilayer based on minimizing the free energy of a mixture of water and lipid molecules. This paper extends previous work by Blom & Peletier (2004) in the following ways. (a) It formulates a more physical model of the hydrophobic effect to facilitate connections with microscale simulations. (b) It clarifies the meaning of the model parameters. (c) It outlines a method for determining parameter values so that physically-realistic bilayer density profiles can be obtained, for example for use in macroscale simulations. Points (a)-(c) suggest that the model has potential to robustly connect some micro- and macroscale levels of multiscale blood flow simulations. The mathematical modelling in point (a) is based upon a consideration of the underlying physics of inter-molecular forces. The governing equations thus obtained are minimized by gradient flows via a novel numerical approach; this enables point (b). The numerical results are shown to behave physically in terms of the effect of background concentration, in contrast to the earlier model which is shown here to not display the expected behaviour. A "short-tail" approximation of the lipid molecules also gives an analytical tool which yields critical values of some parameters under certain conditions. Point (c) involves the first quantitative comparison of the numerical data with physical experimental results.


title: Dynamical systems gradient method for solving ill-conditioned linear algebraic systems

abstract: A version of the Dynamical Systems Method (DSM) for solving ill-conditioned linear algebraic systems is studied in this paper. An {\it a priori} and {\it a posteriori} stopping rules are justified. An algorithm for computing the solution using a spectral decomposition of the left-hand side matrix is proposed. Numerical results show that when a spectral decompositon of the left-hand side matrix is available or not computationally expensive to obtain the new method can be considered as an alternative to the Variational Regularization.


title: Stochastic Tamed 3D Navier-Stokes Equations: Existence, Uniqueness and Ergodicity

abstract: In this paper, we prove the existence of a unique strong solution to a stochastic tamed 3D Navier-Stokes equation in the whole space as well as in the periodic boundary case. Then, we also study the Feller property of solutions, and prove the existence of invariant measures for the corresponding Feller semigroup in the case of periodic conditions. Moreover, in the case of periodic boundary and degenerated additive noise, using the notion of asymptotic strong Feller property proposed by Hairer and Mattingly \cite{Ha-Ma}, we prove the uniqueness of invariant measures for the corresponding transition semigroup.


title: New phase transition in $Pr_{1-x}Ca_xMnO_3$ system: evidence for electrical polarization in charge ordered manganites

abstract: A detailed study of the electric field gradient (EFG) across the Pr1-xCaxMnO3 phase diagram and its temperature dependence is given. Clearly, distinct EFG behaviours for samples outside/inside the charge order (CO) region are observed. The EFG temperature dependence evidences a new phase transition occurring over the broad CO region of the phase diagram. This transition is discontinuous and occurs at temperatures between the charge ordering and the Neel temperatures. The features observed in the EFG are associated with polar atomic vibrations which eventually lead to a spontaneous local electric polarization below CO transition.


title: Magnetic field evolution of the quasiparticle interference in a d-wave superconductor

abstract: Quasiparticle interference in a d-wave superconductor with weak disorder produces distinctive peaks in the Fourier-transformed local density of states measured by scanning tunneling spectroscopy. We predict that amplitudes of these peaks can be enhanced or suppressed by applied magnetic field according to a very specific pattern governed by the symmetry of the superconducting order parameter. This calculated pattern agrees with the recent experimental measurement and suggests that the technique could be useful for probing the underlying normal state at high fields.


title: Methods for measuring energy dissipation rate in anisotropic turbulence

abstract: Energy dissipation rate is an important parameter for nearly every experiment on turbulent flow. Mathematically precise relationships between energy dissipation rate and other measurable statistics for the case of anisotropic turbulence are useful to experimentalists. Such relationships are obtained for which the measurable statistics are the 3rd-order and 2nd-order velocity structure functions as well as the acceleration-velocity structure function. The relationships are derived using the Navier-Stokes equation without approximation. Approximate versions are obtained on the basis of local stationarity and local homogeneity. The latter are valid for arbitrary Reynolds numbers for the case of stationary, homogeneous turbulence. Precise use of the mathematics requires care noted in the Discussion section.


title: Localization theorems in topological Hochschild homology and topological cyclic homology

abstract: We construct localization cofiber sequences for the topological Hochschild homology (THH) and topological cyclic homology (TC) of spectral categories. Using a global construction of the THH and TC of a scheme in terms of the perfect complexes in a spectrally enriched version of the category of unbounded complexes, the sequences specialize to localization cofiber sequences associated to the inclusion of an open subscheme. These are the targets of the cyclotomic trace from the localization sequence of Thomason-Trobaugh in K-theory. We also deduce versions of Thomason's blow-up formula and the projective bundle formula for THH and TC.


title: Using Layout Information for Spreadsheet Visualization

abstract: This paper extends a spreadsheet visualization technique by using layout information. The original approach identifies logically or semantically related cells by relying exclusively on the content of cells for identifying semantic classes. A disadvantage of semantic classes is that users have to supply parameters which describe the possible shapes of these blocks. The correct parametrization requires a certain degree of experience and is thus not suitable for untrained users. To avoid this constraint, the approach reported in this paper uses row/column-labels as well as common format information for locating areas with common, recurring semantics. Heuristics are provided to distinguish between cell groups with intended common semantics and cell groups related in an ad-hoc manner.


title: Spreadsheet Structure Discovery with Logic Programming

abstract: Our term "structure discovery" denotes the recovery of structure, such as the grouping of cells, that was intended by a spreadsheet's author but is not explicit in the spreadsheet. We are implementing structure discovery tools in the logic-programming language Prolog for our spreadsheet analysis program Model Master, by writing grammars for spreadsheet structures. The objective is an "intelligent structure monitor" to run beside Excel, allowing users to reconfigure spreadsheets to the representational needs of the task at hand. This could revolutionise spreadsheet "best practice". We also describe a formulation of spreadsheet reverse-engineering based on "arrows".


title: Generalized entanglement in static and dynamic quantum phase transitions

abstract: We investigate a class of one-dimensional, exactly solvable anisotropic XY spin-1/2 models in an alternating transverse magnetic field from an entanglement perspective. We find that a physically motivated Lie-algebraic generalized entanglement measure faithfully portraits the static phase diagram -- including second- and fourth-order quantum phase transitions belonging to distinct universality classes. In the simplest time-dependent scenario of a slow quench across a quantum critical point, we identify parameter regimes where entanglement exhibits universal dynamical scaling relative to the static limit.


title: Strichartz estimates on Schwarzschild black hole backgrounds

abstract: We study dispersive properties for the wave equation in the Schwarzschild space-time. The first result we obtain is a local energy estimate. This is then used, following the spirit of earlier work of Metcalfe-Tataru, in order to establish global-in-time Strichartz estimates. A considerable part of the paper is devoted to a precise analysis of solutions near the trapping region, namely the photon sphere.


title: Hidden particle production at the ILC

abstract: In a class of new physics models, new physics sector is completely or partly hidden, namely, singlet under the Standard Model (SM) gauge group. Hidden fields included in such new physics models communicate with the Standard Model sector through higher dimensional operators. If a cutoff lies in the TeV range, such hidden fields can be produced at future colliders. We consider a scalar filed as an example of the hidden fields. Collider phenomenology on this hidden scalar is similar to that of the SM Higgs boson, but there are several features quite different from those of the Higgs boson. We investigate productions of the hidden scalar at the International Linear Collider (ILC) and study the feasibility of its measurements, in particular, how well the ILC distinguishes the scalar from the Higgs boson, through realistic Monte Carlo simulations.


title: Simulating the ISM in Global Disk Galaxies

abstract: Until recently, simulations that modeled entire galaxies were restricted to an isothermal or fixed 2- or 3-phase interstellar medium (ISM). This obscured the full role of the ISM in shaping the observed galactic-scale star formation relations. In particular the Kennicutt relation suggests that star formation rates depend in a simple way on global galactic quantities, such as mean gas mass surface density and dynamical time. Contrary to this, observations of nearby star-forming regions, including images from the Spitzer telescope, show that all the way down to ``local-scales'' ($\sim$ parsecs) star formation is a highly clustered process with the gas existing at a wide range of densities and multiple phases, and seemingly decoupled from the larger-scale Galaxy. Many different physical processes appear to be influencing the star formation rate, including heating and cooling of the gas, ``turbulent'' energy injection from outflows, winds and supernova blastwaves, ionization, and magnetic field support. Here we present results from global simulations of disk galaxy star formation that include a fully multiphase ISM and compare to simpler isothermal models. We also investigate the effect of background heating. We discuss the process of giant molecular cloud formation as one of the possible links between the galactic and local scales of star formation.


title: Population imbalanced fermions in harmonically trapped optical lattices

abstract: The attractive Fermi-Hubbard Hamiltonian is solved via the Bogoliubov-de Gennes formalism to analyze the ground state phases of population imbalanced fermion mixtures in harmonically trapped two-dimensional optical lattices. In the low density limit the superfluid order parameter modulates in the radial direction towards the trap edges to accommodate the unpaired fermions that are pushed away from the trap center with a single peak in their density. However in the high density limit while the order parameter modulates in the radial direction towards the trap center for low imbalance, it also modulates towards the trap edges with increasing imbalance until the superfluid to normal phase transition occurs beyond a critical imbalance. This leads to a single peak in the density of unpaired fermions for low and high imbalance but leads to double peaks for intermediate imbalance.


title: A Link between Binomial Parameters and Means of Bounded Random Variables

abstract: In this paper, we establish a fundamental connection between binomial parameters and means of bounded random variables. Such connection finds applications in statistical inference of means of bounded variables.


title: El Roque de Los Muchachos Site Characteristics. III. Analysis of Atmospheric Dust and Aerosol Extinction

abstract: Canary Islands are normally interested by dominant North-East winds that, in some meteorological conditions, can transport sand at high altitude from the Sahara desert. The dust may affect the efficiency of the telescopes and decreases the transparency of the sky. In order to maximize the scientific return of the telescopes located at the ORM, we present an analysis of the atmospheric dust content and its effects on astronomical observations. B, V and I dust aerosol astronomical extinction are derived. Using a 5 years series database of data taken from the four channel TNG dust monitor, we compute a mean hourly and daily values of the dust content. We have detected particles having size 0.3, 0.5, 1.0 and 5.0 um. Using a power law we have derived the content of 10.0 um particles. We found a typical local dust concentration ranging from 3x10^6 particles per cubic meter at 0.3 um, to 10^3 at 5.0 um and 10 at 10.0 um, increasing up to 3 order of magnitudes during the dust storms, with a relative higher increase of 1.0, 5.0 and 10.0 um particles. The number of local dust storm events is the same in winter- and summertime, but, the average background and storm-related increases in the dust concentration in summer are significantly higher than in winter. In a uniform approximation, during the dust storms, an average height of the dust layer of 2.5 km above the telescope is inferred. During the sand storms La Palma Island is affected by an almost uniform layer extending up to 5 km above the sea level, down, at least the height of the telescope. The visible extinction is dominated by particles at 1.0, 5.0 and 10.0 um. In agreement with the results from Carlsberg Automatic Meridian Circle (CAMC) we find a typical extinction during dust storms of about 0.2 mag/airmass.


title: Generating functions for colored 3D Young diagrams and the Donaldson-Thomas invariants of orbifolds

abstract: We derive two multivariate generating functions for three-dimensional Young diagrams (also called plane partitions). The variables correspond to a colouring of the boxes according to a finite Abelian subgroup G of SO(3). We use the vertex operator methods of Okounkov--Reshetikhin--Vafa for the easy case G = Z/n; to handle the considerably more difficult case G=Z/2 x Z/2, we will also use a refinement of the author's recent q--enumeration of pyramid partitions. In the appendix, we relate the diagram generating functions to the Donaldson-Thomas partition functions of the orbifold C^3/G. We find a relationship between the Donaldson-Thomas partition functions of the orbifold and its G-Hilbert scheme resolution. We formulate a crepant resolution conjecture for the Donaldson-Thomas theory of local orbifolds satisfying the Hard Lefschetz condition.


title: Adaptation dynamics of the quasispecies model

abstract: We study the adaptation dynamics of an initially maladapted population evolving via the elementary processes of mutation and selection. The evolution occurs on rugged fitness landscapes which are defined on the multi-dimensional genotypic space and have many local peaks separated by low fitness valleys. We mainly focus on the Eigen's model that describes the deterministic dynamics of an infinite number of self-replicating molecules. In the stationary state, for small mutation rates such a population forms a {\it quasispecies} which consists of the fittest genotype and its closely related mutants. The quasispecies dynamics on rugged fitness landscape follow a punctuated (or step-like) pattern in which a population jumps from a low fitness peak to a higher one, stays there for a considerable time before shifting the peak again and eventually reaches the global maximum of the fitness landscape. We calculate exactly several properties of this dynamical process within a simplified version of the quasispecies model.


title: Belief Propagation and Loop Series on Planar Graphs

abstract: We discuss a generic model of Bayesian inference with binary variables defined on edges of a planar graph. The Loop Calculus approach of [1, 2] is used to evaluate the resulting series expansion for the partition function. We show that, for planar graphs, truncating the series at single-connected loops reduces, via a map reminiscent of the Fisher transformation [3], to evaluating the partition function of the dimer matching model on an auxiliary planar graph. Thus, the truncated series can be easily re-summed, using the Pfaffian formula of Kasteleyn [4]. This allows to identify a big class of computationally tractable planar models reducible to a dimer model via the Belief Propagation (gauge) transformation. The Pfaffian representation can also be extended to the full Loop Series, in which case the expansion becomes a sum of Pfaffian contributions, each associated with dimer matchings on an extension to a subgraph of the original graph. Algorithmic consequences of the Pfaffian representation, as well as relations to quantum and non-planar models, are discussed.


title: On bi-hamiltonian geometry of the Lagrange top

abstract: We consider three different incompatible bi-Hamiltonian structures for the Lagrange top, which have the same foliation by symplectic leaves. These bivectors may be associated with the different 2-coboundaries in the Poisson-Lichnerowicz cohomology defined by canonical bivector on $e^*(3)$.


title: Delta Expansion on the Lattice and Dilated Scaling Region

abstract: A new kind of delta expansion is applied on the lattice to the d=2 non-linear sigma model at N=infinity and N=1 which corresponds to the Ising model. We introduce the parameter delta for the dilation of the scaling region of the model with the replacement of the lattice spacing a to (1-delta)^{1/2}a. Then, we demonstrate that the expansion in delta admits an approximation of the scaling behavior of the model at both limits of N from the information at a large lattice spacing a.


title: Super-Weyl invariance in 5D supergravity

abstract: We propose a superspace formulation for the Weyl multiplet of N=1 conformal supergravity in five dimensions. The corresponding superspace constraints are invariant under super-Weyl transformations generated by a real scalar parameter. The minimal supergravity multiplet, which was introduced by Howe in 1981, emerges if one couples the Weyl multiplet to an Abelian vector multiplet and then breaks the super-Weyl invariance by imposing the gauge condition W=1, with W the field strength of the vector multiplet. The geometry of superspace is shown to allow the existence of a large family of off-shell supermultiplets that possess uniquely determined super-Weyl transformation laws and can be used to describe supersymmetric matter. Many of these supermultiplets have not appeared within the superconformal tensor calculus. We formulate a manifestly locally supersymmetric and super-Weyl invariant action principle. In the super-Weyl gauge W=1, this action reduces to that constructed in arXiv:0712.3102. We also present a superspace formulation for the dilaton Weyl multiplet.


title: Dirac particles tunneling from BTZ black hole

abstract: We calculated the Dirac Particles' Hawking radiation from the outer horizon of BTZ black hole via tunneling formalism. Applying WKB approximation to the Dirac equation in (2+1) dimensional BTZ spacetime background, we obtain the radiation spectrum for fermions and Hawking temperature of BTZ black hole. The results obtained by taking the fermion tunneling into account is consistent with the previous literatures.