title: The coolest 'stars' are free-floating planets

abstract: We show that the coolest known object that is probably formed in a star-like mode is a free-floating planet. We discovered recently that the free-floating planetary mass object OTS,44 (M9.5, ~12 Jupiter masses, age ~2 Myr) has significant accretion and a substantial disk. This demonstrates that the processes that characterize the canonical star-like mode of formation apply to isolated objects down to a few Jupiter masses. We detected in VLT/SINFONI spectra that OTS44 has strong, broad, and variable Paschen beta emission. This is the first evidence for active accretion of a free-floating planet. The object allows us to study accretion and disk physics at the extreme and can be seen as free-floating analog of accreting planets that orbit stars. Our analysis of OTS44 shows that the mass-accretion rate decreases continuously from stars of several solar masses down to free-floating planets. We determined, furthermore, the disk mass (10 Earth masses) and further disk properties of OTS44 through modeling its SED including Herschel far-IR data. We find that objects between 14 and 0.01 solar masses have the same ratio of the disk-to-central-mass of about 1%. Our results suggest that OTS44 is formed like a star and that the increasing number of young free-floating planets and ultra-cool T and Y field dwarfs are the low-mass extension of the stellar population.


title: Infrared Consistency and the Weak Gravity Conjecture

abstract: The weak gravity conjecture (WGC) asserts that an Abelian gauge theory coupled to gravity is inconsistent unless it contains a particle of charge $q$ and mass $m$ such that $q \geq m/m_{\rm Pl}$. This criterion is obeyed by all known ultraviolet completions and is needed to evade pathologies from stable black hole remnants. In this paper, we explore the WGC from the perspective of low-energy effective field theory. Below the charged particle threshold, the effective action describes a photon and graviton interacting via higher-dimension operators. We derive infrared consistency conditions on the parameters of the effective action using i) analyticity of light-by-light scattering, ii) unitarity of the dynamics of an arbitrary ultraviolet completion, and iii) absence of superluminality and causality violation in certain non-trivial backgrounds. For convenience, we begin our analysis in three spacetime dimensions, where gravity is non-dynamical but has a physical effect on photon-photon interactions. We then consider four dimensions, where propagating gravity substantially complicates all of our arguments, but bounds can still be derived. Operators in the effective action arise from two types of diagrams: those that involve electromagnetic interactions (parameterized by a charge-to-mass ratio $q/m$) and those that do not (parameterized by a coefficient $\gamma$). Infrared consistency implies that $q/m$ is bounded from below for small $\gamma$.


title: CLASH-VLT: Insights on the mass substructures in the Frontier Fields Cluster MACS J0416.1-2403 through accurate strong lens modeling

abstract: We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.


title: Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

abstract: We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.


title: Comparing Masses in Literature (CoMaLit)-I. Bias and scatter in weak lensing and X-ray mass estimates of clusters

abstract: The first building block to use galaxy clusters in astrophysics and cosmology is the accurate determination of their mass. Two of the most well regarded direct mass estimators are based on weak lensing (WL) determinations or X-ray analyses assuming hydrostatic equilibrium (HE). By comparing these two mass measurements in samples of rich clusters, we determined the intrinsic scatters, $\sigma_\mathrm{WL}\sim$15 per cent for WL masses and $\sigma_\mathrm{HE}\sim$25 per cent for HE masses. The certain assessment of the bias is hampered by differences as large as $\sim$40 per cent in either WL or HE mass estimates reported by different groups. If the intrinsic scatter in the mass estimate is not considered, the slope of any scaling relation `observable--mass' is biased towards shallower values, whereas the intrinsic scatter of the scaling is over-estimated.


title: CoMaLit - II. The scaling relation between mass and Sunyaev-Zel'dovich signal for Planck selected galaxy clusters

abstract: We discuss the scaling relation between mass and integrated Compton parameter of a sample of galaxy clusters from the all-sky {\it Planck} Sunyaev-Zel'dovich catalogue. Masses were measured with either weak lensing, caustics techniques, or assuming hydrostatic equilibrium. The retrieved $Y_{500}$-$M_{500}$ relation does not strongly depend on the calibration sample. We found a slope of 1.4-1.9, in agreement with self-similar predictions, with an intrinsic scatter of $20\pm10$ per cent. The absolute calibration of the relation can not be ascertained due to systematic differences of $\sim$20-40 per cent in mass estimates reported by distinct groups. Due to the scatter, the slope of the conditional scaling relation, to be used in cosmological studies of number counts, is shallower, $\sim$1.1-1.6. The regression methods employed account for intrinsic scatter in the mass measurements too. We found that Planck mass estimates suffer from a mass dependent bias.


title: Non-linear bond-operator theory and 1/d expansion for coupled-dimer magnets I: Paramagnetic phase

abstract: For coupled-dimer Heisenberg magnets, a paradigm of magnetic quantum phase transitions, we develop a systematic expansion in 1/d, the inverse number of space dimensions. The expansion employs a formulation of the bond-operator technique and is based on the observation that a suitably chosen product-state wavefunction yields exact zero-temperature expectation values of local observables in the d->infty limit, with corrections vanishing as 1/d. We demonstrate the approach for a model of dimers on a hypercubic lattice, which generalizes the square-lattice bilayer Heisenberg model to arbitrary d. In this paper, we use the 1/d expansion to calculate static and dynamic observables at zero temperature in the paramagnetic singlet phase, up to the quantum phase transition, and compare the results with numerical data available for d=2. Contact is also made with previously proposed refinements of bond-operator theory as well as with a perturbative expansion in the inter-dimer coupling. In a companion paper, the present 1/d expansion will be extended to the ordered phase, where it is shown to consistently describe the entire phase diagram including the quantum critical point.


title: Two Virasoro symmetries in stringy warped AdS$_3$

abstract: We study three-dimensional consistent truncations of type IIB supergravity which admit warped AdS$_3$ solutions. These theories contain subsectors that have no bulk dynamics. We show that the symplectic form for these theories, when restricted to the non-dynamical subsectors, equals the symplectic form for pure Einstein gravity in AdS$_3$. Consequently, for each consistent choice of boundary conditions in AdS$_3$, we can define a consistent phase space in warped AdS$_3$ with identical conserved charges. This way, we easily obtain a Virasoro $\times$ Virasoro asymptotic symmetry algebra in warped AdS$_3$; two different types of Virasoro $\times$ Ka\v{c}-Moody symmetries are also consistent alternatives. Next, we study the phase space of these theories when propagating modes are included. We show that, as long as one can define a conserved symplectic form without introducing instabilities, the Virasoro $\times$ Virasoro asymptotic symmetries can be extended to the entire (linearized) phase space. This implies that, at least at semi-classical level, consistent theories of gravity in warped AdS$_3$ are described by a two-dimensional conformal field theory, as long as stability is not an issue.


title: The Bishop-Phelps-Bollob\'as property for operators from $\mathcal{C}(K)$ to uniformly convex spaces

abstract: We show that the pair $(C(K),X)$ has the Bishop-Phelps-Bolloba\'as property for operators if $K$ is a compact Hausdorff space and $X$ is a uniformly convex space.


title: The Density Tur\'an problem

abstract: Let $H$ be a graph on $n$ vertices and let the blow-up graph $G[H]$ be defined as follows. We replace each vertex $v_i$ of $H$ by a cluster $A_i$ and connect some pairs of vertices of $A_i$ and $A_j$ if $(v_i,v_j)$ was an edge of the graph $H$. As usual, we define the edge density between $A_i$ and $A_j$ as $d(A_i,A_j)=\frac{e(A_i,A_j)}{|A_i||A_j|}.$ We study the following problem. Given densities $\gamma_{ij}$ for each edge $(i,j)\in E(H)$. Then one has to decide whether there exists a blow-up graph $G[H]$ with edge densities at least $\gamma_{ij}$ such that one cannot choose a vertex from each cluster so that the obtained graph is isomorphic to $H$, i.e, no $H$ appears as a transversal in $G[H]$. We call $d_{crit}(H)$ the maximal value for which there exists a blow-up graph $G[H]$ with edge densities $d(A_i,A_j)=d_{crit}(H)$ $((v_i,v_j)\in E(H))$ not containing $H$ in the above sense. Our main goal is to determine the critical edge density and to characterize the extremal graphs.


title: Bound-state formation for thermal relic dark matter and unitarity

abstract: We show that the relic abundance of thermal dark matter annihilating via a long-range interaction, is significantly affected by the formation and decay of dark matter bound states in the early universe, if the dark matter mass is above a few TeV. We determine the coupling required to obtain the observed dark matter density, taking into account both the direct 2-to-2 annihilations and the formation of bound states, and provide an analytical fit. We argue that the unitarity limit on the inelastic cross-section is realized only if dark matter annihilates via a long-range interaction, and we determine the upper bound on the mass of thermal-relic dark matter to be about 197 (139) TeV for (non)-self-conjugate dark matter.


title: The environment of low redshift quasar pairs

abstract: We investigate the properties of the galaxy environment of a sample of 14 low redshift (z $<$ 0.85) quasar physical pairs extracted from SDSS DR10 archives. The pairs have a systemic radial velocity difference $\Delta V_\parallel \leqslant$ 600 $km \ s^{-1}$ (based on [OIII]5007 \AA \ line) and projected distance $ R_\bot \leqslant$ 600 kpc. The physical association of the pairs is statistically confirmed at a level of $\sim$ 90 %. For most of the images of these quasars we are able to resolve their host galaxies that turn out to be on average similar to those of quasars not in pairs. We also found that quasars in a pair are on average in region of modest galaxy overdensity extending up 0.5 Mpc from the QSO. This galaxy overdensity is indistinguishable from that of a homogeneous sample of isolated quasars at the same redshift and with similar host galaxy luminosity. These results, albeit derived from a small (but homogeneous) sample of objects, suggest that the rare activation of two quasars with small physical separation does not require any extraordinary environment.


title: Scaling Laws for Dark Matter Halos in Late-Type and Dwarf Spheroidal Galaxies

abstract: Dark matter (DM) halos of Sc-Im galaxies satisfy scaling laws analogous to the fundamental plane relations for elliptical galaxies. Halos in less luminous galaxies have smaller core radii, higher central densities, and smaller central velocity dispersions. If dwarf spheroidal (dSph) and dwarf Magellanic irregular (dIm) galaxies lie on the extrapolations of these correlations, then we can estimate their baryon loss relative to that of brighter Sc-Im galaxies. We find that, if there had been no such enhanced baryon loss, then typical dSph and dIm galaxies would be brighter in absolute magnitude by 4 and 3.5 mag, respectively. Instead, these galaxies lost or retained as gas (in dIm galaxies) baryons that could have formed stars. Also, typical dSph and dIm galaxies have DM halos that are more massive than we thought, with velocity dispersions of about 30 km/s or circular-orbit rotation velocities of V_circ ~ 42 km/s. Comparison of DM and visible matter correlations confirms that, at V-band absolute magnitudes fainter than -18, dSph and dIm galaxies form a sequence of decreasing baryon-to-DM mass ratios in smaller dwarfs. We show explicitly that galaxy baryon content goes to (almost) zero at halo V_circ = 42 +- 4 km/s, in agreement with what we found from our estimate of baryon depletion. Our results suggest that there may be a large population of DM halos that are essentially dark and undiscovered. This helps to solve the problem that the fluctuation spectrum of cold DM predicts more dwarfs than we observe.


title: Noise Model of Relaxation Oscillators Due to Feedback Regeneration Based on Physical Phase Change

abstract: A new approach to investigate noise spikes due to regeneration in a relaxation oscillator is proposed. Noise spikes have not been satisfactorily accounted for in traditional phase noise models. This paper attempts to explain noise spikes/jump phenomenon by viewing it as phase change in the thermodynamic system(for example, from gas to liquid or magnetization of ferromagnet). Both are due to regeneration (positive feedback in oscillator as well as alignment of spin due to positive feedback in ferromagnet). The mathematical tool used is the partition function in thermodynamics, and the results mapped between thermodynamic system and relaxation oscillator. Theory is developed and formula derived to predict the magnitude of the jump, as a function of design parameter such as regeneration parameter or loop gain. Formulas show that noise increases sharply as regeneration parameter/loop gain approaches one, in much the same way when temperature approaches critical temperature in phase change. Simulations on circuits (Eldo) using CMOS as well as Monte Carlo simulations (Metropolis) on ferromagnet (Ising model) were performed and both show jump behaviour consistent with formula. Measurements on relaxation oscillators fabricated in 0.13um CMOS technology verify such behaviour, where the sharp increase in noise when regeneration parameter/loop gain is close to one, matches closely with the theoretical formula. Using the formula the designer can quantify the variation of noise spikes dependency on design parameters such as gm (device transconductance), R, I0, via their influence on regeneration parameter/loop gain.


title: Genera of non-algebraic leaves of polynomial foliations of $\mathbb C^2$

abstract: In this article, we prove two results. First, we construct a dense subset in the space of polynomial foliations of degree $n$ such that each foliation from this subset has a leaf with at least $\frac{(n+1)(n+2)}2-4$ handles. Next, we prove that for a generic foliation invariant under the map $(x, y)\mapsto (x, -y)$ all leaves have infinitely many handles.


title: Hot Jupiters and Cool Stars

abstract: Close-in planets are in jeopardy as their host stars evolve off the main sequence to the subgiant and red giant phases. In this paper, we explore the influences of the stellar mass (in the range 1.5--2\Mso ), mass-loss prescription, planet mass (from Neptune up to 10 Jupiter masses), and eccentricity, on the orbital evolution of planets as their parent stars evolve to become subgiants and Red Giants. We find that planet engulfment during the Red Giant Branch is not very sensitive to the stellar mass or mass-loss rates adopted in the calculations, but quite sensitive to the planetary mass. The range of initial separations for planet engulfment increases with decreasing mass-loss rates or stellar mass and increasing planetary masses. Regarding the planet's orbital eccentricity, we find that as the star evolves into the red giant phase, stellar tides start to dominate over planetary tides. As a consequence, a transient population of moderately eccentric close-in Jovian planets is created, that otherwise would have been expected to be absent from main sequence stars. We find that very eccentric and distant planets do not experience much eccentricity decay, and that planet engulfment is primarily determined by the pericenter distance and the maximum stellar radius.


title: Prediction of magnetic moment collapse in ZrFe$_2$ under hydrostatic pressure

abstract: Electronic structure and magnetic properties of ZrFe$_2$ with a cubic Laves phase are investigated by calculations based on the density functional theory. The total magnetic moment (m) of 3.14 $\mu_B$ per formula unit (\emph{f.u.}) is obtained at the experimental lattice constant (7.06 \AA), which is larger than 3.06 $\mu_B$/\emph{f.u.} obtained at the theoretical equilibrium lattice constant (6.85 \AA). The localized $3d$ magnetic moment is in negative diffusive sp background moment. We predict a two-step magnetic collapse under pressure: one is from 3.06 $\mu_B$/\emph{f.u.} to 1.26 $\mu_B$/\emph{f.u.} at about 3.6 GPa, and the other is from 0.5 $\mu_B$/\emph{f.u.} to nonmagnetic state at about 15 GPa. We understand this process by the changes of density of states. The magnetic moment decreases under the pressure in the vicinity of the experimental lattice constant with $d\ln m/dp=-0.038$ GPa$^{-1}$. The spontaneous volume magnetostriction is 0.015. We suggest that the Invar effect of this alloy may be understood when considering the magnetic moment variation according to the Weiss $2\gamma$-model.


title: Electronic Structure of Liquid Water and a Platinum Surface

abstract: Many-body perturbation theory within the G$_0$W$_0$ approximation is used to determine molecular orbital level alignment at a liquid water/Pt(111) interface generated through $ab~ initio$ molecular dynamics. Molecular orbital energy levels are shown to depend both on the position of H$_2$O molecules within the liquid relative to the surface and the details of their local bonding environment. Standard density functional theory calculations disagree qualitatively with level alignment predicted by many-body perturbation theory.


title: Distributed Maximum Matching in Bounded Degree Graphs

abstract: We present deterministic distributed algorithms for computing approximate maximum cardinality matchings and approximate maximum weight matchings. Our algorithm for the unweighted case computes a matching whose size is at least $(1-\eps)$ times the optimal in $\Delta^{O(1/\eps)} + O\left(\frac{1}{\eps^2}\right) \cdot\log^*(n)$ rounds where $n$ is the number of vertices in the graph and $\Delta$ is the maximum degree. Our algorithm for the edge-weighted case computes a matching whose weight is at least $(1-\eps)$ times the optimal in $\log(\min\{1/\wmin,n/\eps\})^{O(1/\eps)}\cdot(\Delta^{O(1/\eps)}+\log^*(n))$ rounds for edge-weights in $[\wmin,1]$. The best previous algorithms for both the unweighted case and the weighted case are by Lotker, Patt-Shamir, and Pettie~(SPAA 2008). For the unweighted case they give a randomized $(1-\eps)$-approximation algorithm that runs in $O((\log(n)) /\eps^3)$ rounds. For the weighted case they give a randomized $(1/2-\eps)$-approximation algorithm that runs in $O(\log(\eps^{-1}) \cdot \log(n))$ rounds. Hence, our results improve on the previous ones when the parameters $\Delta$, $\eps$ and $\wmin$ are constants (where we reduce the number of runs from $O(\log(n))$ to $O(\log^*(n))$), and more generally when $\Delta$, $1/\eps$ and $1/\wmin$ are sufficiently slowly increasing functions of $n$. Moreover, our algorithms are deterministic rather than randomized.


title: The one-loop worldsheet S-matrix for the AdS(n) x S(n) x T(10-2n) superstring

abstract: We compute the massive-sector worldsheet S-matrix for superstring theories in AdS(n) x S(n) x T(10-2n) (with n=2,3,5) in the near BMN expansion up to one-loop order in inverse string tension. We show that, after taking into account the wave function renormalization, the one-loop S-matrix is UV finite. In an appropriate regularization scheme the S-matrix is consistent with the underlying symmetries of the superstring theory, i.e. for the n=3,5 cases it coincides with the one implied by the light-cone gauge symmetries with the dressing phases determined from the crossing equations. For the n=2,3 cases we observe that the massless modes decouple from the one-loop calculation of massive mode scattering, i.e. the 2n-dimensional supercoset sigma model and the full 10-dimensional superstring happen to have the same massive one-loop S-matrix.


title: Response to Comment [arXiv:1407.6854] on Palmer et al., Nature, 510, 385, 2014

abstract: We respond to a Comment [arXiv:1407.6854 (2014)] on our recent Nature paper [Nature, 510, 385 (2014)]. We categorically disprove the arguments provided in arXiv:1407.6854 (2014) and thereby further substantiate the evidence we presented in our recent study, demonstrating the existence of a metastable liquid-liquid transition in a molecular model of water. We will make our code publicly available shortly along with proper user documentation that is currently under development.


title: Gravitational collapse of charged scalar fields

abstract: In order to study the gravitational collapse of charged matter we analyze the simple model of an self-gravitating massless scalar field coupled to the electromagnetic field in spherical symmetry. The evolution equations for the Maxwell-Klein-Gordon sector are derived in the 3+1 formalism, and coupled to gravity by means of the stress-energy tensor of these fields. To solve consistently the full system we employ a generalized Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of General Relativity that is adapted to spherical symmetry. We consider two sets of initial data that represent a time symmetric spherical thick shell of charged scalar field, and differ by the fact that one set has zero global electrical charge while the other has non-zero global charge. For compact enough initial shells we find that the configuration doesn't disperse and approaches a final state corresponding to a sub-extremal Reissner-N\"ordstrom black hole with $|Q|<M$. By increasing the fundamental charge of the scalar field $q$ we find that the final black hole tends to become more and more neutral. Our results support the cosmic censorship conjecture for the case of charged matter.


title: Refined Rotational Period, Pole Solution & Shape Model for (3200) Phaethon

abstract: (3200) Phaethon exhibits both comet- and asteroid-like properties, suggesting it could be a rare transitional object such as a dormant comet or previously volatile-rich asteroid. This justifies detailed study of (3200) Phaethon's physical properties, as a better understanding of asteroid-comet transition objects can provide insight into minor body evolution. We therefore acquired time-series photometry of (3200) Phaethon over 15 nights from 1994 to 2013, primarily using the Tektronix 2048x2048 pixel CCD on the University of Hawaii 2.2-m telescope. We utilized light curve inversion to: (1) refine (3200) Phaethon's rotational period to P=3.6032+/-0.0008 h; (2) estimate a rotational pole orientation of lambda=+85+/-13 degrees and beta=-20+/-10 degrees; and (3) derive a shape model. We also used our extensive light curve dataset to estimate the slope parameter of (3200) Phaethon's phase curve as G~0.06, consistent with C-type asteroids. We discuss how this highly oblique pole orientation with a negative ecliptic latitude supports previous evidence for (3200) Phaethon's origin in the inner main asteroid belt as well as the potential for deeply buried volatiles fueling impulsive yet rare cometary outbursts.


title: Going for Speed: Sublinear Algorithms for Dense r-CSPs

abstract: We give new sublinear and parallel algorithms for the extensively studied problem of approximating n-variable r-CSPs (constraint satisfaction problems with constraints of arity r up to an additive error. The running time of our algorithms is O(n/\epsilon^2) + 2^O(1/\epsilon^2) for Boolean r-CSPs and O(k^4 n / \epsilon^2) + 2^O(log k / \epsilon^2) for r-CSPs with constraints on variables over an alphabet of size k. For any constant k this gives optimal dependence on n in the running time unconditionally, while the exponent in the dependence on 1/\epsilon is polynomially close to the lower bound under the exponential-time hypothesis, which is 2^\Omega(\epsilon^(-1/2)). For Max-Cut this gives an exponential improvement in dependence on 1/\epsilon compared to the sublinear algorithms of Goldreich, Goldwasser and Ron (JACM'98) and a linear speedup in n compared to the algorithms of Mathieu and Schudy (SODA'08). For the maximization version of k-Correlation Clustering problem our running time is O(k^4 n / \epsilon^2) + k^O(1/\epsilon^2), improving the previously best n k^{O(1/\epsilon^3 log k/\epsilon) by Guruswami and Giotis (SODA'06).


title: Occupation times of long-range exclusion and connections to KPZ class exponents

abstract: With respect to a class of long-range exclusion processes on $\mathbb{Z}^d$, with single particle transition rates of order $|\cdot|^{-(d+\alpha)}$, starting under Bernoulli invariant measure $\nu_\rho$ with density $\rho$, we consider the fluctuation behavior of occupation times at a vertex and more general additive functionals. Part of our motivation is to investigate the dependence on $\alpha$, $d$ and $\rho$ with respect to the variance of these functionals and associated scaling limits. In the case the rates are symmetric, among other results, we find the scaling limits exhaust a range of fractional Brownian motions with Hurst parameter $H\in [1/2,3/4]$. However, in the asymmetric case, we study the asymptotics of the variances, which when $d=1$ and $\rho=1/2$ points to a curious dichotomy between long-range strength parameters $0<\alpha\leq 3/2$ and $\alpha>3/2$. In the former case, the order of the occupation time variance is the same as under the process with symmetrized transition rates, which are calculated exactly. In the latter situation, we provide consistent lower and upper bounds and other motivations that this variance order is the same as under the asymmetric short-range model, which is connected to KPZ class scalings of the space-time bulk mass density fluctuations.


title: Cooperation and Storage Tradeoffs in Power-Grids with Renewable Energy Resources

abstract: One of the most important challenges in smart grid systems is the integration of renewable energy resources into its design. In this work, two different techniques to mitigate the time varying and intermittent nature of renewable energy generation are considered. The first one is the use of storage, which smooths out the fluctuations in the renewable energy generation across time. The second technique is the concept of distributed generation combined with cooperation by exchanging energy among the distributed sources. This technique averages out the variation in energy production across space. This paper analyzes the trade-off between these two techniques. The problem is formulated as a stochastic optimization problem with the objective of minimizing the time average cost of energy exchange within the grid. First, an analytical model of the optimal cost is provided by investigating the steady state of the system for some specific scenarios. Then, an algorithm to solve the cost minimization problem using the technique of Lyapunov optimization is developed and results for the performance of the algorithm are provided. These results show that in the presence of limited storage devices, the grid can benefit greatly from cooperation, whereas in the presence of large storage capacity, cooperation does not yield much benefit. Further, it is observed that most of the gains from cooperation can be obtained by exchanging energy only among a few energy harvesting sources.


title: Variability Analysis of Complex Networks Measures based on Stochastic Distances

abstract: Complex networks can model the structure and dynamics of different types of systems. It has been shown that they are characterized by a set of measures. In this work, we evaluate the variability of complex networks measures face to perturbations and, for this purpose, we impose controlled perturbations and quantify their effect. We analyze theoretical models (random, small-world and scale-free) and real networks (a collaboration network and a metabolic networks) along with the shortest path length, vertex degree, local cluster coefficient and betweenness centrality measures. In such analysis, we propose the use of three stochastic quantifiers: the Kullback-Leibler divergence and the Jensen-Shannon and Hellinger distances. The sensitivity of these measures was analyzed with respect to the following perturbations: edge addition, edge removal, edge rewiring and node removal, all of them applied at different intensities. The results reveal that the evaluated measures are influenced by these perturbations. Additionally, hypotheses tests were performed to verify the behavior of the degree distribution to identify the intensity of the perturbations that leads to break this property.


title: Planck-scale soccer-ball problem: a case of mistaken identity

abstract: Over the last decade it has been found that nonlinear laws of composition of momenta are predicted by some alternative approaches to "real" 4D quantum gravity, and by all formulations of dimensionally-reduced (3D) quantum gravity coupled to matter. The possible relevance for rather different quantum-gravity models has motivated several studies, but this interest is being tempered by concerns that a nonlinear law of addition of momenta might inevitably produce a pathological description of the total momentum of a macroscopic body. I here show that such concerns are unjustified, finding that they are rooted in failure to appreciate the differences between two roles for laws composition of momentum in physics. Previous results relied exclusively on the role of a law of momentum composition in the description of spacetime locality. However, the notion of total momentum of a multi-particle system is not a manifestation of locality, but rather reflects translational invariance. By working within an illustrative example of quantum spacetime I show explicitly that spacetime locality is indeed reflected in a nonlinear law of composition of momenta, but translational invariance still results in an undeformed linear law of addition of momenta building up the total momentum of a multi-particle system.


abstract: Inspired by methods of N. P. Smart, we describe an algorithm to determine all Picard curves over Q with good reduction away from 3, up to Q-isomorphism. A correspondence between the isomorphism classes of such curves and certain quintic binary forms possessing a rational linear factor is established. An exhaustive list of integral models is determined, and an application to a question of Ihara is discussed.


title: On-sky vibration environment for the Gemini Planet Imager and mitigation effort

abstract: The Gemini Planet Imager (GPI) entered on-sky commissioning and had its first-light at the Gemini South (GS) telescope in November 2013. GPI is an extreme adaptive optics (XAO), high-contrast imager and integral-field spectrograph dedicated to the direct detection of hot exo-planets down to a Jupiter mass. The performance of the apodized pupil Lyot coronagraph depends critically upon the residual wavefront error (design goal of 60 nm RMS with 5 mas RMS tip/tilt), and therefore is most sensitive to vibration (internal or external) of Gemini's instrument suite. Excess vibration can be mitigated by a variety of methods such as passive or active dampening at the instrument or telescope structure or Kalman filtering of specific frequencies with the AO control loop. Understanding the sources, magnitudes and impact of vibration is key to mitigation. This paper gives an overview of related investigations based on instrument data (GPI AO module) as well as external data from accelerometer sensors placed at different locations on the GS telescope structure. We report the status of related mitigation efforts, and present corresponding results.


title: A primality criterion based on a Lucas' congruence

abstract: Let $p$ be a prime. In 1878 \'{E}. Lucas proved that the congruence $$ {p-1\choose k}\equiv (-1)^k\pmod{p}$$ holds for any nonnegative integer $k\in\{0,1,\ldots,p-1\}$. The converse statement was given in Problem 1494 of {\it Mathematics Magazine} proposed in 1997 by E. Deutsch and I.M. Gessel. In this note we generalize this converse assertion by the following result: If $n>1$ and $q>1$ are integers such that $$ {n-1\choose k}\equiv (-1)^k \pmod{q}$$ for every integer $k\in\{0,1,\ldots, n-1\}$, then $q$ is a prime and $n$ is a power of $q$.


title: On-sky low order non-common path correction of the GPI Calibration Unit

abstract: The Gemini Planet Imager (GPI) entered on-sky commissioning phase, and had its First Light at the Gemini South telescope in November 2013. Meanwhile, the fast loops for atmospheric correction of the Extreme Adaptive Optics (XAO) system have been closed on many dozen stars at different magnitudes (I=4-8), elevation angles and a variety of seeing conditions, and a stable loop performance was achieved from the beginning. Ultimate contrast performance requires a very low residual wavefront error (design goal 60 nm RMS), and optimization of the planet finding instrument on different ends has just begun to deepen and widen its dark hole region. Laboratory raw contrast benchmarks are in the order of 10^-6 or smaller. In the telescope environment and in standard operations new challenges are faced (changing gravity, temperature, vibrations) that are tackled by a variety of techniques such as Kalman filtering, open-loop models to keep alignment to within 5 mas, speckle nulling, and a calibration unit (CAL). The CAL unit was especially designed by the Jet Propulsion Laboratory to control slowly varying wavefront errors at the focal plane of the apodized Lyot coronagraph by the means of two wavefront sensors: 1) a 7x7 low order Shack-Hartmann SH wavefront sensor (LOWFS), and 2) a special Mach-Zehnder interferometer for mid-order spatial frequencies (HOWFS) - atypical in that the beam is split in the focal plane via a pinhole but recombined in the pupil plane with a beamsplitter. The original design goal aimed for sensing and correcting on a level of a few nm which is extremely challenging in a telescope environment. This paper focuses on non-common path low order wavefront correction as achieved through the CAL unit on sky. We will present the obtained results as well as explain challenges that we are facing.


title: Emergent gravitational dynamics in relativistic Bose--Einstein condensate

abstract: Analogue models of gravity have played a pivotal role in the past years by providing a test bench for many open issues in quantum field theory in curved spacetime such as the robustness of Hawking radiation and cosmological particle production. More recently, the same models have offered a valuable framework within which current ideas about the emergence of spacetime and its dynamics could be discussed via convenient toy models. In this context, we study here an analogue gravity system based on a relativistic Bose--Einstein condensate. We show that in a suitable limit this system provides not only an example of an emergent spacetime (with a massive and a massless relativistic fields propagating on it) but also that such spacetime is governed by an equation with geometric meaning that takes the familiar form of Nordstr{\"o}m theory of gravitation. In this equation the gravitational field is sourced by the expectation value of the trace of the effective stress energy tensor of the quasiparticles while the Newton and cosmological constants are functions of the fundamental scales of the microscopic system. This is the first example of analogue gravity in which a Lorentz invariant, geometric theory of semiclassical gravity emerges from an underlying quantum theory of matter in flat spacetime.


title: Measurement of D-meson production in pp, p-Pb, and Pb-Pb collisions at the LHC with the ALICE detector

abstract: Heavy quarks are a powerful probe for investigating the properties of the Quark-Gluon Plasma created in heavy-ion collisions, since they are produced in initial hard scattering processes and experience all the stages of the medium evolution. ALICE has measured the production of $\mathrm{D}^0$, $\mathrm{D}^+$, $\mathrm{D}^{*+}$, and $\mathrm{D}_{\mathrm{s}}^+$ mesons at central rapidity in their hadronic decay channels in various collision systems and energies. We present recent results for D-meson production measured by the ALICE Collaboration in pp collisions at $\sqrt{s}=7$ and 2.76 TeV, Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=2.76$ TeV and p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV.


title: Linear continuous functionals on FN-type spaces

abstract: By using the space of fuzzy numbers, in e.g. [5] have been considered several complete metric spaces (called here {\bf FN}-type spaces) endowed with addition and scalar multiplication, such that the metrics have nice properties but the spaces are not linear, i.e. are not groups with respect to addition and the scalar multiplication is not, in general, distributive with respect to usual scalar addition. This paper deals with the form of linear continuous functionals defined on these spaces.


title: Failure of steady state thermodynamics in lattice gases under nonuniform drive

abstract: To be useful, steady state thermodynamics (SST) must be self-consistent and have predictive value. Although consistency of SST was recently verified for driven lattice gases under global weak exchange, I show here that it does not predict the coexisting densities in athermal stochastic lattice gases under a nonuniform drive. I consider the lattice gas with nearest-neighbor exclusion on the square lattice, with nearest-neighbor hopping (NNE dynamics), and with hopping to both nearest and next-nearest neighbors (NNE2 dynamics). Part of the system is subject to a drive $D$ that favors hopping along one direction, while the other part is free of the drive. Thus the steady state represents coexistence between two subsystems, one far from equilibrium and the other in equilibrium, which exchange particles along the interfaces. The dimensionless chemical potential $\mu^*(\rho,D) = \mu/k_B T$ for the lattice gas with density $\rho$ is readily determined in studies of a uniform system. Under the nonuniform drive, however, equating the chemical potentials of the coexisting subsystems does not yield the coexisting densities. The steady state chemical potential is, moreover, different in the coexisting bulk regions, contrary to the basic principles of thermodynamics. These results cast serious doubt on the predictive value of SST.


title: Strict Self-Assembly of Fractals using Multiple Hands

abstract: In this paper we consider the problem of the strict self-assembly of infinite fractals within tile self-assembly. In particular, we provide tile assembly algorithms for the assembly of the discrete Sierpinski triangle and the discrete Sierpinski carpet within a class of models we term the \emph{$h$-handed assembly model} ($h$-HAM), which generalizes the 2-HAM to allow up to $h$ assemblies to combine in a single assembly step. Despite substantial consideration, no purely growth self-assembly model has yet been shown to strictly assemble an infinite fractal without significant modification to the fractal shape. In this paper we not only achieve this, but in the case of the Sierpinski carpet are able to achieve it within the 2-HAM, one of the most well studied tile assembly models in the literature. Our specific results are as follows. We provide a $6$-HAM construction for the Sierpinski triangle that works at scale factor 1, 30 tile types, and assembles the fractal in a \emph{near perfect} way in which all intermediate assemblies are finite-sized iterations of the recursive fractal. We further assemble the Sierpinski triangle within the $3$-HAM at scale factor 3 and 990 tile types. For the Sierpinski carpet, we present a $2$-HAM result that works at scale factor 3 and uses 1,216 tile types. We further include analysis showing that the aTAM is incapable of strictly assembling the non-tree Sierpinski triangle considered in this paper, and that multiple hands are needed for the near-perfect assembly of the Sierpinski triangle and the Sierpinski carpet.


title: User-friendly, theory-based web applet for rapidly predicting structure and thermodynamics of complex fluids

abstract: Based on a recently introduced analytical strategy [Hollingshead et al., J. Chem. Phys. 139, 161102 (2013)], we present a web applet that can quickly and semi-quantitatively estimate the equilibrium radial distribution function and related thermodynamic properties of a fluid from knowledge of its pair interaction. We present a detailed description of the applet's features and intended workflow, followed by a description of how the applet can be used to illustrate two (of many possible) concepts of interest for introductory statistical mechanics courses: the transition from ideal gas-like behavior to correlated-liquid behavior with increasing density and the tradeoff between dominant length scales with changing temperature in a system with ramp-shaped repulsions. The latter type of interaction qualitatively captures distinctive thermodynamic properties of liquid water because its energetic bias toward locally open structures mimics that of water's hydrogen-bond network.


title: Short effective intervals containing primes

abstract: We prove that if $x$ is large enough, namely $x\ge x_0$, then there exists a prime between $x(1- \Delta^{-1})$ and $x$, where $\Delta$ is an effective constant computed in terms of $x_0$. This improves some previous results of Ramar\'e and Saouter and is used to verify the Ternary Goldbach Conjecture for the integers up to $7.86\cdot 10^{27}$.


title: Poisson structure and stability analysis of a coupled system arising from the supersymmetric breaking of Super KdV

abstract: The Poisson structure of a coupled system arising from a supersymmetric breaking of N=1 Super KdV equations is obtained. The supersymmetric breaking is implemented by introducing a Clifford algebra instead of a Grassmann algebra. The Poisson structure follows from the Dirac brackets obtained by the constraint analysis of the hamiltonian of the system. The coupled system has multisolitonic solutions. We show that the one soliton solutions are Liapunov stable.


title: Anomalous Behavior of Dark States in Quantum Gases of 6Li

abstract: We create atom-molecule dark states in a degenerate Fermi gas of 6Li in both weakly and strongly interacting regimes using two-photon Raman scattering to couple fermion pairs to bound molecular states in the ground singlet and triplet potential. Near the unitarity point in the BEC-BCS crossover regime, the atom number revival height associated with the dark state abruptly and unexpectedly decreases and remains low for magnetic fields below the Feshbach resonance center at 832.2~G. With a weakly interacting Fermi gas at 0~G we perform precision dark-state spectroscopy of the least-bound vibrational levels of the lowest singlet and triplet potentials. From these spectra, we obtain binding energies of the $v''=9$, $N''=0$ level of the $a(1^3\Sigma_{u}^+)$ potential and the $v''=38$, $N''=0$ level of the $X(1^1\Sigma_{g}^+)$ potential with absolute uncertainty as low as $20$ kHz. For the triplet potential we resolve the molecular hyperfine structure.


title: The Spectrum and Morphology of the Fermi Bubbles

abstract: The Fermi bubbles are two large structures in the gamma-ray sky extending to $55^\circ$ above and below the Galactic center. We analyze 50 months of Fermi Large Area Telescope data between 100 MeV and 500 GeV above $10^\circ$ in Galactic latitude to derive the spectrum and morphology of the Fermi bubbles. We thoroughly explore the systematic uncertainties that arise when modeling the Galactic diffuse emission through two separate approaches. The gamma-ray spectrum is well described by either a log parabola or a power law with an exponential cutoff. We exclude a simple power law with more than 7$\sigma$ significance. The power law with an exponential cutoff has an index of $1.9 \pm 0.2$ and a cutoff energy of $110\pm 50$ GeV. We find that the gamma-ray luminosity of the bubbles is $4.4^{+2.4}_{-0.9} \times 10^{37}$ erg s$^{-1}$. We confirm a significant enhancement of gamma-ray emission in the south-eastern part of the bubbles, but we do not find significant evidence for a jet. No significant variation of the spectrum across the bubbles is detected. The width of the boundary of the bubbles is estimated to be $3.4^{+3.7}_{-2.6}$ deg. Both inverse Compton (IC) models and hadronic models including IC emission from secondary leptons fit the gamma-ray data well. In the IC scenario, the synchrotron emission from the same population of electrons can also explain the WMAP and Planck microwave haze with a magnetic field between 5 and 20 $\mu$G.


title: How Auto-Encoders Could Provide Credit Assignment in Deep Networks via Target Propagation

abstract: We propose to exploit {\em reconstruction} as a layer-local training signal for deep learning. Reconstructions can be propagated in a form of target propagation playing a role similar to back-propagation but helping to reduce the reliance on derivatives in order to perform credit assignment across many levels of possibly strong non-linearities (which is difficult for back-propagation). A regularized auto-encoder tends produce a reconstruction that is a more likely version of its input, i.e., a small move in the direction of higher likelihood. By generalizing gradients, target propagation may also allow to train deep networks with discrete hidden units. If the auto-encoder takes both a representation of input and target (or of any side information) in input, then its reconstruction of input representation provides a target towards a representation that is more likely, conditioned on all the side information. A deep auto-encoder decoding path generalizes gradient propagation in a learned way that can could thus handle not just infinitesimal changes but larger, discrete changes, hopefully allowing credit assignment through a long chain of non-linear operations. In addition to each layer being a good auto-encoder, the encoder also learns to please the upper layers by transforming the data into a space where it is easier to model by them, flattening manifolds and disentangling factors. The motivations and theoretical justifications for this approach are laid down in this paper, along with conjectures that will have to be verified either mathematically or experimentally, including a hypothesis stating that such auto-encoder mediated target propagation could play in brains the role of credit assignment through many non-linear, noisy and discrete transformations.


title: Hamiltonian structure of an operator valued extension of Super KdV equations

abstract: An extension of the super Korteweg-de Vries integrable system in terms of operator valued functions is obtained. In particular the extension contains the $N=1$ Super KdV and coupled systems with functions valued on a symplectic space. We introduce a Miura transformation for the extended system and obtain its hamiltonian structure. We also obtain an extended Gardner transformation which allows to find an infinite number of conserved quantities of the extended system.


title: Social Conformity Despite Individual Preferences for Distinctiveness

abstract: We demonstrate that individual behaviors directed at the attainment of distinctiveness can in fact produce complete social conformity. We thus offer an unexpected generative mechanism for this central social phenomenon. Specifically, we establish that agents who have fixed needs to be distinct and adapt their positions to achieve distinctiveness goals, can nevertheless self-organize to a limiting state of absolute conformity. This seemingly paradoxical result is deduced formally from a small number of natural assumptions, and is then explored at length computationally. Interesting departures from this conformity equilibrium are also possible, including divergence in positions. The effect of extremist minorities on these dynamics is discussed. A simple extension is then introduced, which allows the model to generate and maintain social diversity, including multimodal distinctiveness distributions. The paper contributes formal definitions, analytical deductions, and counterintuitive findings to the literature on individual distinctiveness and social conformity.


title: Using arbitrary parity-check matrices for quantum error correction assisted by less noisy qubits

abstract: Recently a framework for assisted quantum error correction was proposed in which a specific type of error is allowed to occur on auxiliary qubits, which is in contrast to standard entanglement assistance that requires noiseless auxiliary qubits. However, while the framework maintains the ability to import any binary or quaternary linear code without sacrificing active error correction power, it requires the code designer to turn a parity-check matrix of the underlying classical code into an equivalent one in standard form. This means that classical coding theoretic techniques that require parity-check matrices to be in specific form may not fully be exploitable. Another issue of the recently proposed scheme is that the error correction capabilities for bit errors and phase errors are generally equal, which is not ideal for asymmetric error models. This paper addresses these two problems. We generalize the framework in such a way that any parity-check matrix of any binary or quaternary linear code can be exploited. Our generalization also allows for importing a pair of distinct linear codes so that error correction capabilities become suitably asymmetric.


title: $\mathop{\rm PL}_+(I)$ is not a Polish group

abstract: The group $\mathop{\rm PL}_+(I)$ of increasing piecewise linear self-homeomorphisms of the interval $I=[0,1]$ may not be assigned a topology in such a way that it becomes a Polish group. The same statement holds for the groups $\mathop{\rm Homeo}_+^{Lip}(I)$ of bi-Lipschitz homeomorphisms of $I$, and $\mathop{\rm Diff}_+^{1+\epsilon}(I)$ of diffeomorphisms of $I$ whose derivatives are H\"older continuous with exponent $\epsilon$.


title: Synthesis and transport properties of ternary type-I Si clathrate K8Al7Si39

abstract: A ternary type-I Si clathrate, K8AlxSi46-x, which is a candidate functional material composed of abundant non-toxic elements, was synthesized and its transport properties were investigated at temperatures ranging from 10 to 320 K. The synthesized compound is confirmed to be the ternary type-I Si clathrate K8Al7Si39 with a lattice parameter of a = 10.442 A using neutron powder diffractometry and inductively coupled plasma optical emission spectrometry. Electrical resistivity and Hall coefficient measurements revealed that K8Al7Si39 is a metal with electrons as the dominant carriers at a density of approximately 1x10^27 /m3. The value of Seebeck coefficient for K8Al7Si39 is negative and its absolute value increases with the temperature. The temperature dependence of the thermal conductivity is similar to that for a crystalline solid. The dimensionless figure of merit is approximately 0.01 at 300 K, which is comparable to that for other ternary Si clathrates.


title: Stellar magnetic fields from four Stokes parameter observations

abstract: Magnetic Doppler imaging from four Stokes parameter observations has uncovered a new level of complexity of stellar magnetic fields previously not known. This new information is of interest to the evolution and structure of magnetic fields of intermediate and high-mass stars. Recently, we performed magnetic Doppler imaging study of the chemically peculiar star HD 24712 using phase-resolved spectropolarimetric observations of line-profiles in all four Stokes parameters. The observations were obtained with the HARPSpol instrument on the 3.6-m ESO telescope. Our analysis shows that the magnetic field of HD 24712 has very strong dipolar component with weak small-scale contributions. The finding of dipole-like field for HD 24712 gives evidence for the hypothesis that old Ap stars have predominantly dipolar magnetic fields. We plan on performing magnetic Doppler imaging from four Stokes parameter observations for other Ap/Bp stars as well, which will help us understand the nature of magnetic fields in these stars.


title: Spatiotemporal clustering, climate periodicity, and social-ecological risk factors for dengue during an outbreak in Machala, Ecuador, in 2010

abstract: The objective of this study was to characterize the spatiotemporal dynamics and climatic and social-ecological risk factors associated with the largest dengue epidemic to date in Machala, Ecuador, to inform the development of a dengue EWS. The following data were included in analyses: neighborhood-level georeferenced dengue cases, national census data, and entomological surveillance data from 2010; time series of weekly dengue cases (aggregated to the city-level) and meteorological data from 2003 to 2012. We applied LISA and Morans I to analyze the spatial distribution of the 2010 dengue cases, and developed multivariate logistic regression models through a multi-model selection process to identify census variables and entomological covariates associated with the presence of dengue at the neighborhood level. Using data aggregated at the city-level, we conducted a time-series (wavelet) analysis of weekly climate and dengue incidence (2003-2012) to identify significant time periods (e.g., annual, biannual) when climate co-varied with dengue, and to describe the climate conditions associated with the 2010 outbreak. We found significant hotspots of dengue transmission near the center of Machala. The best-fit model to predict the presence of dengue included older age and female gender of the head of the household, greater access to piped water in the home, poor housing condition, and less distance to the central hospital. Wavelet analyses revealed that dengue transmission co-varied with rainfall and minimum temperature at annual and biannual cycles, and we found that anomalously high rainfall and temperatures were associated with the 2010 outbreak. Our findings highlight the importance of geospatial information in dengue surveillance and the potential to develop a climate-driven spatiotemporal prediction models to inform disease prevention and control interventions.


title: Even and odd Kauffman bracket ideals for genus-1 tangles

abstract: This paper refines previous work by the first author. We study the question of which links in the 3-sphere can be obtained as closures of a given 1-manifold in an unknotted solid torus in the 3-sphere (or genus-1 tangle) by adjoining another 1-manifold in the complementary solid torus. We distinguish between even and odd closures, and define even and odd versions of the Kauffman bracket ideal. These even and odd Kauffman bracket ideals are used to obstruct even and odd tangle closures. Using a basis of Habiro's for the even Kauffman bracket skein module of the solid torus, we define bases for the even and odd skein module of the solid torus relative to two points. These even and odd bases allow us to compute a finite list of generators for the even and odd Kauffman bracket ideals of a genus-1 tangle. We do this explicitly for three examples. Furthermore, we use the even and odd Kauffman bracket ideals to conclude in some cases that the determinants of all even/odd closures of a genus-1 tangle possess a certain divisibility.


title: Metastability, excitations, fluctuations, and multiple-swallowtail structures of a superfluid in a Bose-Einstein condensate in the presence of a uniformly moving defect

abstract: We solve the Gross-Pitaevskii (GP) and Bogoliubov equations to investigate the metastability of superfluidity in a Bose-Einstein condensate in the presence of a uniformly moving defect potential in a two-dimensional torus. We calculate the total energy and momentum as functions of the driving velocity of the moving defect and find metastable states with negative effective-mass near the critical velocity. We also find that the first excited energy (energy gap) in the finite-sized torus closes at the critical velocity, that it obeys one-fourth power-law scaling, and that the dynamical fluctuation of the density (amplitude of the order parameter) is strongly enhanced near the critical velocity. We confirm the validity of our results near the critical velocity by calculating the quantum depletion. We find an unconventional swallowtail structure (multiple-swallowtail structure) through calculations of the unstable stationary solutions of the GP equation.


title: Some new applications of the Stanley-Macdonald Pieri Rules

abstract: In a seminal paper Richard Stanley derived Pieri rules for the Jack symmetric function basis. These rules were extended by Macdonald to his now famous symmetric function basis. The original form of these rules had a forbidding complexity that made them difficult to use in explicit calculations. In the early 90's it was discovered that, due to massive cancellations, the dual rule, which expresses skewing by $e_1$ the modified Macdonald polynomial ${\tilde H}_\mu[X;q,t]$, can be given a very simple combinatorial form in terms of corner weights of the Ferrers' diagram of $\mu$. A similar formula was later obtained by the last named author for the multiplication of ${\tilde H}_\mu[X;q,t]$ by $e_1$, but never published. In the years that followed we have seen some truly remarkable uses of these two Pieri rules in establishing highly non trivial combinatorial results in the Theory of Macdonald polynomials. This theory has recently been spectacularly enriched by various Algebraic Geometrical results in the works of Hikita, Schiffmann, Schiffmann-Vasserot, A. Negut and Gorsky-Negut. This development opens up the challenging task of deriving their results by purely Algebraic Combinatorial methods. In this paper we present the progress obtained by means of Pieri rules.


title: Top-Down Skiplists

abstract: We describe todolists (top-down skiplists), a variant of skiplists (Pugh 1990) that can execute searches using at most $\log_{2-\varepsilon} n + O(1)$ binary comparisons per search and that have amortized update time $O(\varepsilon^{-1}\log n)$. A variant of todolists, called working-todolists, can execute a search for any element $x$ using $\log_{2-\varepsilon} w(x) + o(\log w(x))$ binary comparisons and have amortized search time $O(\varepsilon^{-1}\log w(w))$. Here, $w(x)$ is the "working-set number" of $x$. No previous data structure is known to achieve a bound better than $4\log_2 w(x)$ comparisons. We show through experiments that, if implemented carefully, todolists are comparable to other common dictionary implementations in terms of insertion times and outperform them in terms of search times.


title: Exclusion process with slow boundary

abstract: We study the hydrodynamic and the hydrostatic behavior of the Simple Symmetric Exclusion Process with \emph{slow boundary}. The term \emph{slow boundary} means that particles can be born or die at the boundary sites, at a rate proportional to $N^{-\theta}$, where $\theta > 0$ and $N$ is the scaling parameter. In the bulk, the particles exchange rate is equal to $1$. In the hydrostatic scenario, we obtain three different linear profiles, depending on the value of the parameter $\theta$; in the hydrodynamic scenario, we obtain that the time evolution of the spatial density of particles, in the diffusive scaling, is given by the weak solution of the heat equation, with boundary conditions that depend on $ \theta $. If $\theta \in(0,1)$, we get Dirichlet boundary conditions, (which is the same behavior if $\theta=0$, see \cite{f}); if $\theta=1$, we get Robin boundary conditions; and, if $\theta\in(1,\infty)$, we get Neumann boundary conditions.


title: Particle Motion in Monopoles and Geodesics on Cones

abstract: The equations of motion of a charged particle in the field of Yang's $\mathrm{SU}(2)$ monopole in 5-dimensional Euclidean space are derived by applying the Kaluza-Klein formalism to the principal bundle $\mathbb{R}^8\setminus\{0\}\to\mathbb{R}^5\setminus\{0\}$ obtained by radially extending the Hopf fibration $S^7\to S^4$, and solved by elementary methods. The main result is that for every particle trajectory $\mathbf{r}:I\to\mathbb{R}^5\setminus\{0\}$, there is a 4-dimensional cone with vertex at the origin on which $\mathbf{r}$ is a geodesic. We give an explicit expression of the cone for any initial conditions.


title: Lebesgue and Hardy spaces for symmetric norms I

abstract: In this paper, we define and study a class $\mathcal{R}_{c}$ of norms on $L^{\infty}\left( \mathbb{T}\right) $, called $continuous\ rotationally\ symmetric \ norms$, which properly contains the class $\left \{ \left \Vert \cdot \right \Vert _{p}:1\leq p<\infty \right \} .$ For $\alpha \in \mathcal{R}% _{c}$ we define $L^{\alpha}\left( \mathbb{T}\right) $ and the Hardy space $H^{\alpha}\left( \mathbb{T}\right) $, and we extend many of the classical results, including the dominated convergence theorem, convolution theorems, dual spaces, Beurling-type invariant spaces, inner-outer factorizations, characterizing the multipliers and the closed densely-defined operators commuting with multiplication by $z$. We also prove a duality theorem for a version of $L^{\alpha}$ in the setting of von Neumann algebras.


title: Distributed event-triggered coordination for average consensus on weight-balanced digraphs

abstract: This paper proposes a novel distributed event-triggered algorithmic solution to the multi-agent average consensus problem for networks whose communication topology is described by weight-balanced, strongly connected digraphs. The proposed event-triggered communication and control strategy does not rely on individual agents having continuous or periodic access to information about the state of their neighbors. In addition, it does not require the agents to have a priori knowledge of any global parameter. We show that, under the proposed law, events cannot be triggered an infinite number of times in any finite period (i.e., no Zeno behavior), and that the resulting network executions provably converge to the average of the initial agents' states exponentially fast. We also provide weaker conditions on connectivity under which convergence is guaranteed when the communication topology is switching. Simulations illustrate our results.


title: On the canonical volume of 3-folds of general type with $P_{12}\geq 2$

abstract: Let $V$ be a nonsingular projective 3-fold of general type. When the pluricanonical section index $\delta(V)>12$, Chen-Chen \cite{Chen3} has a complete list of the possibility for the weighted basket ${\mathbb B}(V)$. However the possibility of ${\mathbb B}(V)$ could be infinite in the situation $\delta(V)\leq 12$, which is the main challenge to the classification. In this paper we mainly study the case with $\delta(V)=12$ and show that $Vol(V)\geq \frac{31}{48048}$, which improves the corresponding result of Chen--Chen \cite[Prop. 4.9(5)]{Chen3}.


title: Divisibility of Weil Sums of Binomials

abstract: Consider the Weil sum $W_{F,d}(u)=\sum_{x \in F} \psi(x^d+u x)$, where $F$ is a finite field of characteristic $p$, $\psi$ is the canonical additive character of $F$, $d$ is coprime to $|F^*|$, and $u \in F^*$. We say that $W_{F,d}(u)$ is three-valued when it assumes precisely three distinct values as $u$ runs through $F^*$: this is the minimum number of distinct values in the nondegenerate case, and three-valued $W_{F,d}$ are rare and desirable. When $W_{F,d}$ is three-valued, we give a lower bound on the $p$-adic valuation of the values. This enables us to prove the characteristic $3$ case of a 1976 conjecture of Helleseth: when $p=3$ and $[F:{\mathbb F}_3]$ is a power of $2$, we show that $W_{F,d}$ cannot be three-valued.


title: Chance Constrained Optimization for Targeted Internet Advertising

abstract: We introduce a chance constrained optimization model for the fulfillment of guaranteed display Internet advertising campaigns. The proposed formulation for the allocation of display inventory takes into account the uncertainty of the supply of Internet viewers. We discuss and present theoretical and computational features of the model via Monte Carlo sampling and convex approximations. Theoretical upper and lower bounds are presented along with a numerical substantiation.


title: The quark mass gap in strong magnetic fields

abstract: Quarks in strong magnetic fields (|eB|>>Lambda_QCD^2 ~ 0.04 GeV^2) acquire enhanced infrared phase space proportional to |eB|. Accordingly they provide larger chiral condensates and stronger backreactions to the gluon dynamics. Confronting theories with lattice data at various values of |eB|, one can test theoretical ideas as well as validity of various approximations, domain of applicability of the effective models, and so on. The particularly interesting findings on the lattice are inverse magnetic catalysis and linear growth of the chiral condensate as a function of |eB|, which pose theoretical challenges. In this talk we propose a scenario to explain both phenomena, claiming that the quark mass gap should stay at around ~ Lambda_QCD, instead of ~|eB|^{1/2} which has been supposed from dimensional arguments and/or effective model calculations. The contrast between infrared and ultraviolet behaviors of the interaction is a key ingredient to obtain the mass gap of ~Lambda_QCD.


title: OGLE-2013-BLG-0102LA,B: Microlensing binary with components at star/brown-dwarf and brown-dwarf/planet boundaries

abstract: We present the analysis of the gravitational microlensing event OGLE-2013-BLG-0102. The light curve of the event is characterized by a strong short-term anomaly superposed on a smoothly varying lensing curve with a moderate magnification $A_{\rm max}\sim 1.5$. It is found that the event was produced by a binary lens with a mass ratio between the components of $q = 0.13$ and the anomaly was caused by the passage of the source trajectory over a caustic located away from the barycenter of the binary. From the analysis of the effects on the light curve due to the finite size of the source and the parallactic motion of the Earth, the physical parameters of the lens system are determined. The measured masses of the lens components are $M_{1} = 0.096 \pm 0.013~M_{\odot}$ and $M_{2} = 0.012 \pm 0.002~M_{\odot}$, which correspond to near the hydrogen-burning and deuterium-burning mass limits, respectively. The distance to the lens is $3.04 \pm 0.31~{\rm kpc}$ and the projected separation between the lens components is $0.80 \pm 0.08~{\rm AU}$.


title: Nonuniform Dichotomy Spectrum and Normal Forms for Nonautonomous Differential Systems

abstract: The aim of this paper is to study the normal forms of nonautonomous differential systems. For doing so, we first investigate the nonuniform dichotomy spectrum of the linear evolution operators that admit a nonuniform exponential dichotomy, where the linear evolution operators are defined by nonautonomous differential equations $\dot x=A(t)x$ in $\R^n$. Using the nonuniform dichotomy spectrum we obtain the normal forms of the nonautonomous linear differential equations. Finally we establish the finite jet normal forms of the nonlinear differential systems $\dot x=A(t)x+f(t,x)$ in $\R^n$, which is based on the nonuniform dichotomy spectrum and the normal forms of the nonautonomous linear systems.


title: Generating and Solving Symbolic Parity Games

abstract: We present a new tool for verification of modal mu-calculus formulae for process specifications, based on symbolic parity games. It enhances an existing method, that first encodes the problem to a Parameterised Boolean Equation System (PBES) and then instantiates the PBES to a parity game. We improved the translation from specification to PBES to preserve the structure of the specification in the PBES, we extended LTSmin to instantiate PBESs to symbolic parity games, and implemented the recursive parity game solving algorithm by Zielonka for symbolic parity games. We use Multi-valued Decision Diagrams (MDDs) to represent sets and relations, thus enabling the tools to deal with very large systems. The transition relation is partitioned based on the structure of the specification, which allows for efficient manipulation of the MDDs. We performed two case studies on modular specifications, that demonstrate that the new method has better time and memory performance than existing PBES based tools and can be faster (but slightly less memory efficient) than the symbolic model checker NuSMV.


title: Strategic Port Graph Rewriting: An Interactive Modelling and Analysis Framework

abstract: We present strategic portgraph rewriting as a basis for the implementation of visual modelling and analysis tools. The goal is to facilitate the specification, analysis and simulation of complex systems, using port graphs. A system is represented by an initial graph and a collection of graph rewriting rules, together with a user-defined strategy to control the application of rules. The strategy language includes constructs to deal with graph traversal and management of rewriting positions in the graph. We give a small-step operational semantics for the language, and describe its implementation in the graph transformation and visualisation tool PORGY.


title: Entity-Linking via Graph-Distance Minimization

abstract: Entity-linking is a natural-language-processing task that consists in identifying the entities mentioned in a piece of text, linking each to an appropriate item in some knowledge base; when the knowledge base is Wikipedia, the problem comes to be known as wikification (in this case, items are wikipedia articles). One instance of entity-linking can be formalized as an optimization problem on the underlying concept graph, where the quantity to be optimized is the average distance between chosen items. Inspired by this application, we define a new graph problem which is a natural variant of the Maximum Capacity Representative Set. We prove that our problem is NP-hard for general graphs; nonetheless, under some restrictive assumptions, it turns out to be solvable in linear time. For the general case, we propose two heuristics: one tries to enforce the above assumptions and another one is based on the notion of hitting distance; we show experimentally how these approaches perform with respect to some baselines on a real-world dataset.


title: Graph- versus Vector-Based Analysis of a Consensus Protocol

abstract: The Paxos distributed consensus algorithm is a challenging case-study for standard, vector-based model checking techniques. Due to asynchronous communication, exhaustive analysis may generate very large state spaces already for small model instances. In this paper, we show the advantages of graph transformation as an alternative modelling technique. We model Paxos in a rich declarative transformation language, featuring (among other things) nested quantifiers, and we validate our model using the GROOVE model checker, a graph-based tool that exploits isomorphism as a natural way to prune the state space via symmetry reductions. We compare the results with those obtained by the standard model checker Spin on the basis of a vector-based encoding of the algorithm.


title: Specifying and Executing Optimizations for Parallel Programs

abstract: Compiler optimizations, usually expressed as rewrites on program graphs, are a core part of all modern compilers. However, even production compilers have bugs, and these bugs are difficult to detect and resolve. The problem only becomes more complex when compiling parallel programs; from the choice of graph representation to the possibility of race conditions, optimization designers have a range of factors to consider that do not appear when dealing with single-threaded programs. In this paper we present PTRANS, a domain-specific language for formal specification of compiler transformations, and describe its executable semantics. The fundamental approach of PTRANS is to describe program transformations as rewrites on control flow graphs with temporal logic side conditions. The syntax of PTRANS allows cleaner, more comprehensible specification of program optimizations; its executable semantics allows these specifications to act as prototypes for the optimizations themselves, so that candidate optimizations can be tested and refined before going on to include them in a compiler. We demonstrate the use of PTRANS to state, test, and refine the specification of a redundant store elimination optimization on parallel programs.


title: Graph Transformation Planning via Abstraction

abstract: Modern software systems increasingly incorporate self-* behavior to adapt to changes in the environment at runtime. Such adaptations often involve reconfiguring the software architecture of the system. Many systems also need to manage their architecture themselves, i.e., they need a planning component to autonomously decide which reconfigurations to execute to reach a desired target configuration. For the specification of reconfigurations, we employ graph transformations systems (GTS) due to the close relation of graphs and UML object diagrams. We solve the resulting planning problems with a planning system that works directly on a GTS. It features a domain-independent heuristic that uses the solution length of an abstraction of the original problem as an estimate. Finally, we provide experimental results on two different domains, which confirm that our heuristic performs better than another domain-independent heuristic which resembles heuristics employed in related work.


title: Backwards State-space Reduction for Planning in Dynamic Knowledge Bases

abstract: In this paper we address the problem of planning in rich domains, where knowledge representation is a key aspect for managing the complexity and size of the planning domain. We follow the approach of Description Logic (DL) based Dynamic Knowledge Bases, where a state of the world is represented concisely by a (possibly changing) ABox and a (fixed) TBox containing the axioms, and actions that allow to change the content of the ABox. The plan goal is given in terms of satisfaction of a DL query. In this paper we start from a traditional forward planning algorithm and we propose a much more efficient variant by combining backward and forward search. In particular, we propose a Backward State-space Reduction technique that consists in two phases: first, an Abstract Planning Graph P is created by using the Abstract Backward Planning Algorithm (ABP), then the abstract planning graph P is instantiated into a corresponding planning graph P by using the Forward Plan Instantiation Algorithm (FPI). The advantage is that in the preliminary ABP phase we produce a symbolic plan that is a pattern to direct the search of the concrete plan. This can be seen as a kind of informed search where the preliminary backward phase is useful to discover properties of the state-space that can be used to direct the subsequent forward phase. We evaluate the effectiveness of our ABP+FPI algorithm in the reduction of the explored planning domain by comparing it to a standard forward planning algorithm and applying both of them to a concrete business case study.


title: Existence of strong-pairing quantum Hall phase in bilayer cold atom systems with dipolar interactions

abstract: We study bilayer fermionic cold atom systems with dipolar interactions, as well as a two-component tunable pseudopotential (TCTP) model which keeps only the zeroth and first Haldane pseudopotentials, at total Landau level filling factor 1/2. Our numerical results on the TCTP model indicates that Haldane-Rezayi state describes the critical point between strong and weak d-wave pairing quantum Hall phases. Further increasing the attractive zeroth pseudopotentials, the system transits from the strong-pairing phase to a stripe phase, and then to a cluster phase (or phase separation). The dipolar interaction can be mapped onto the TCTP model in the strong-pairing phase, if high order pseudopotentials are ignored. Our numerical results show that this is indeed the case, so the strong-pairing phase can be realized in the cold atom system.


title: Kinetic Simulations of Rayleigh-Taylor Instabilities

abstract: We report on an ongoing project to develop a large scale Direct Simulation Monte Carlo code. The code is primarily aimed towards applications in astrophysics such as simulations of core-collapse supernovae. It has been tested on shock wave phenomena in the continuum limit and for matter out of equilibrium. In the current work we focus on the study of fluid instabilities. Like shock waves these are routinely used as test-cases for hydrodynamic codes and are discussed to play an important role in the explosion mechanism of core-collapse supernovae. As a first test we study the evolution of a single-mode Rayleigh-Taylor instability at the interface of a light and a heavy fluid in the presence of a gravitational acceleration. To suppress small-wavelength instabilities caused by the irregularity in the separation layer we use a large particle mean free path. The latter leads to the development of a diffusion layer as particles propagate from one fluid into the other. For small amplitudes, when the instability is in the linear regime, we compare its position and shape to the analytic prediction. Despite the broadening of the fluid interface we see a good agreement with the analytic solution. At later times we observe the development of a mushroom like shape caused by secondary Kelvin-Helmholtz instabilities as seen in hydrodynamic simulations and consistent with experimental observations.


title: Learning Economic Parameters from Revealed Preferences

abstract: A recent line of work, starting with Beigman and Vohra (2006) and Zadimoghaddam and Roth (2012), has addressed the problem of {\em learning} a utility function from revealed preference data. The goal here is to make use of past data describing the purchases of a utility maximizing agent when faced with certain prices and budget constraints in order to produce a hypothesis function that can accurately forecast the {\em future} behavior of the agent. In this work we advance this line of work by providing sample complexity guarantees and efficient algorithms for a number of important classes. By drawing a connection to recent advances in multi-class learning, we provide a computationally efficient algorithm with tight sample complexity guarantees ($\Theta(d/\epsilon)$ for the case of $d$ goods) for learning linear utility functions under a linear price model. This solves an open question in Zadimoghaddam and Roth (2012). Our technique yields numerous generalizations including the ability to learn other well-studied classes of utility functions, to deal with a misspecified model, and with non-linear prices.


title: HD314884: A Slowly Pulsating B star in a Close Binary

abstract: We present the results of a spectroscopic and photometric analysis of HD314884, a slowly pulsating B star (SPB) in a binary system with detected soft X-ray emission. We spectrally classify the B star as a B5V-B6V star with T$_{eff}$ = 15,490 $\pm$ 310 K, log $g$ = 3.75 $\pm$ 0.25 dex, and a photometric period of P$_{0}$ = 0.889521(12) days. A spectroscopic period search reveals an orbital period for the system of P$_{orb}$ = 1.3654(11) days. The discrepancy in the two periods and the identification of a second and third distinct frequency in the photometric fourier transform at P$_1$ = 3.1347(56) and P$_2$ = 1.517(28) days provides evidence that HD314884 is a slowly pulsating B star (SPB) with at least three oscillation frequencies. These frequencies appear to originate from higher-order, non-linear tidal pulsations. Using the dynamical parameters obtained from the radial velocity curve, we find the most probable companion mass to be M$_1$ = $\sim$0.8 M$_{\odot}$ assuming a typical mass for the B star and most probable inclination. We conclude that the X-ray source companion to HD314884 is most likely a coronally active G-type star or a white dwarf (WD), with no apparent emission lines in the optical spectrum. The mass probability distribution of the companion star mass spans 0.6-2.3 M$_{\odot}$ at 99$\%$ confidence which allows the possibility of a neutron star companion. The X-ray source is unlikely to be a black hole unless it is of a very low mass or low binary inclination.


title: Quasi-particle energy spectra in local reduced density matrix functional theory

abstract: Recently, we introduced (e-print arXiv:1407.7128) {\em local reduced density matrix functional theory} (local RDMFT), a theoretical scheme capable of incorporating static correlation effects in Kohn-Sham equations. Here, we apply local RDMFT to molecular systems of relatively large size, as a demonstration of its computational efficiency and its accuracy in predicting single-electron properties from the eigenvalue spectrum of the single-particle Hamiltonian with a local effective potential. We present encouraging results on the photoelectron spectrum of molecular systems and the relative stability of C$_{20}$ isotopes. In addition, we propose a modelling of the fractional occupancies as functions of the orbital energies that further improves the efficiency of the method useful in applications to large systems and solids.


title: Room Temperature Spin Pumping in Topological Insulator Bi2Se3

abstract: Three-dimensional (3D) topological insulators are known for their strong spin-orbit coupling and the existence of spin-textured topological surface states which could be potentially exploited for spintronics. Here, we investigate spin pumping from a metallic ferromagnet (CoFeB) into a 3D topological insulator (Bi2Se3) and demonstrate successful spin injection from CoFeB into Bi2Se3 and the direct detection of the electromotive force generated by the inverse spin Hal effect (ISHE) at room temperature. The spin pumping, driven by the magnetization dynamics of the metallic ferromagnet, introduces a spin current into the topological insulator layer, resulting in a broadening of the ferromagnetic resonance (FMR) linewidth. We find that the FMR linewidth more than quintuples, the spin mixing conductance can be as large as 3.4*10^20m^-2 and the spin Hall angle can be as large as 0.23 in the Bi2Se3 layer.


title: Global structure of quaternion polynomial differential equations

abstract: In this paper we mainly study the global structure of the quaternion Bernoulli equations $\dot q=aq+bq^n$ for $q\in \mathbb H$ the quaternion field and also some other form of cubic quaternion differential equations. By using the Liouvillian theorem of integrability and the topological characterization of $2$--dimensional torus: orientable compact connected surface of genus one, we prove that the quaternion Bernoulli equations may have invariant tori, which possesses a full Lebesgue measure subset of $\mathbb H$. Moreover, if $n=2$ all the invariant tori are full of periodic orbits; if $n=3$ there are nfiinitely many invariant tori fulfilling periodic orbits and also infinitely many invariant ones fulfilling dense orbits.


title: Inverse Jacobian multipliers and Hopf bifurcation on center manifolds

abstract: In this paper we consider a class of higher dimensional differential systems in $\mathbb R^n$ which have a two dimensional center manifold at the origin with a pair of pure imaginary eigenvalues. First we characterize the existence of either analytic or $C^\infty$ inverse Jacobian multipliers of the systems around the origin, which is either a center or a focus on the center manifold. Later we study the cyclicity of the system at the origin through Hopf bifurcation by using the vanishing multiplicity of the inverse Jacobian multiplier.


title: Second random-phase approximation, Thouless' theorem and the stability condition reexamined and clarified

abstract: It has been revealed through numerical calculations that the Second Random Phase Approximation (SRPA) with the Hartree-Fock solution as its reference state results in 1) spurious states at genuinely finite energy, contrary to common expectation, and 2) unstable solutions, which within the first-order Random Phase Approximation correspond to real low-energy collective vibrations. In the present work, these shortcomings of SRPA are shown to not contradict Thouless' theorem about the energy-weighted sum rule, and their origin is traced to the violation of the stability condition. A more general theorem is proven. Formal arguments are elucidated through numerical examples. Implications for the validity of SRPA are discussed.


title: Analytic normalization of analytically integrable differential systems near a periodic orbit

abstract: For an analytic differential system in $\mathbb R^n$ with a periodic orbit, we will prove that if the system is analytically integrable around the periodic orbit, i.e. it has $n-1$ functionally independent analytic first integrals defined in a neighborhood of the periodic orbit, then the system is analytically equivalent to its Poincar\'e--Dulac type normal form. This result is an extension for analytic integrable differential systems around a singularity to the ones around a periodic orbit.


title: Analytic integrable systems: Analytic normalization and embedding flows

abstract: In this paper we mainly study the existence of analytic normalization and the normal form of finite dimensional complete analytic integrable dynamical systems. More details, we will prove that any complete analytic integrable diffeomorphism $F(x)=Bx+f(x)$ in $(\mathbb C^n,0)$ with $B$ having eigenvalues not modulus $1$ and $f(x)=O(|x|^2)$ is locally analytically conjugate to its normal form. Meanwhile, we also prove that any complete analytic integrable differential system $\dot x=Ax+f(x)$ in $(\mathbb C^n,0)$ with $A$ having nonzero eigenvalues and $f(x)=O(|x|^2)$ is locally analytically conjugate to its normal form. Furthermore we will prove that any complete analytic integrable diffeomorphism defined on an analytic manifold can be embedded in a complete analytic integrable flow. We note that parts of our results are the improvement of Moser's one in {\it Comm. Pure Appl. Math.} 9$($1956$)$, 673--692 and of Poincar\'e's one in {\it Rendiconti del circolo matematico di Palermo} 5$($1897$)$, 193--239. These results also improve the ones in {\it J. Diff. Eqns.} 244$($2008$)$, 1080--1092 in the sense that the linear part of the systems can be nonhyperbolic, and the one in {\it Math. Res. Lett.} 9$($2002$)$, 217--228 in the way that our paper presents the concrete expression of the normal form in a restricted case.


title: The $16$th Hilbert problem on algebraic limit cycles

abstract: For real planar polynomial differential systems there appeared a simple version of the $16$th Hilbert problem on algebraic limit cycles: {\it Is there an upper bound on the number of algebraic limit cycles of all polynomial vector fields of degree $m$?} In [J. Differential Equations, 248(2010), 1401--1409] Llibre, Ram\'irez and Sadovskia solved the problem, providing an exact upper bound, in the case of invariant algebraic curves generic for the vector fields, and they posed the following conjecture: {\it Is $1+(m-1)(m-2)/2$ the maximal number of algebraic limit cycles that a polynomial vector field of degree $m$ can have?} In this paper we will prove this conjecture for planar polynomial vector fields having only nodal invariant algebraic curves. This result includes the Llibre {\it et al}\,'s as a special one. For the polynomial vector fields having only non--dicritical invariant algebraic curves we answer the simple version of the 16th Hilbert problem.


title: A comparison of possibilities to measure the coordinates of a moving hot body using an infrared telescope or Doppler radar

abstract: A comparison of the sensitivities of methods which allow us to determine the coordinates of a moving hot body is made.


title: Generalized rational first integrals of analytic differential systems

abstract: In this paper we mainly study the necessary conditions for the existence of functionally independent generalized rational first integrals of ordinary differential systems via the resonances. The main results extend some of the previous related ones, for instance the classical Poincar\'e's one \cite{Po}, the Furta's one, part of Chen's ones, and the Shi's one. The key point in the proof of our main results is that functionally independence of generalized rational functions implies the functionally independence of their lowest order rational homogeneous terms.


title: The embedding flows of $C^\infty$ hyperbolic diffeomorphisms

abstract: In [{\it American J. Mathematics}, 124(2002), 107--127] we proved that for a germ of $C^\infty$ hyperbolic diffeomorphisms $F(x)=Ax+f(x)$ in $(\mathbb R^n,0)$, if $A$ has a real logarithm with its eigenvalues weakly nonresonant, then $F(x)$ can be embedded in a $C^\infty$ autonomous differential system. Its proof was very complicated, which involved the existence of embedding periodic vector field of $F(x)$ and the extension of the Floquet's theory to nonlinear $C^\infty$ periodic differential systems. In this paper we shall provide a simple and direct proof to this last result. Next we shall show that the weakly nonresonant condition in the last result on the real logarithm of $A$ is necessary for some $C^\infty$ diffeomorphisms $F(x)=Ax+f(x)$ to have $C^\infty$ embedding flows. Finally we shall prove that a germ of $C^\infty$ hyperbolic diffeomorphisms $F(x)=Ax+f(x)$ with $f(x)=O(|x|^2)$ in $(\mathbb R^2,0)$ has a $C^\infty$ embedding flow if and only if either $A$ has no negative eigenvalues or $A$ has two equal negative eigenvalues and it can be diagonalizable.


title: Quantum phases of attractive bosons on a Bose-Hubbard ladder with three-body constraint

abstract: We obtain the complete quantum phase diagram of bosons on a two-leg ladder in the presence of attractive onsite and repulsive interchain nearest-neighbor interactions by imposing the onsite three-body constraint. We find three distinct phases; namely, the atomic superfluid (ASF), dimer superfluid (DSF), and the dimer rung insulator (DRI). In the absence of the interchain nearest-neighbor repulsion, the system exhibits a transition from the ASF to the DSF phase with increasing onsite attraction. However, the presence of the interchain nearest-neighbor repulsion stabilizes a gapped DRI phase, which is flanked by the DSF phase. We also obtain the phase diagram of the system for different values of the interchain nearest-neighbor interaction. By evaluating different order parameters, we obtain the complete phase diagram and the properties of the phase transitions using the self-consistent cluster mean-field theory.


title: Stability Analysis of Ecomorphodynamic Equations

abstract: Although riparian vegetation is present in or along many water courses of the world, its active role resulting from the interaction with flow and sediment processes has only recently become an active field of research. Especially, the role of vegetation in the process of river pattern formation has been explored and demonstrated mostly experimentally and numerically until now. In the present work, we shed light on this subject by performing a linear stability analysis on a simple model for riverbed vegetation dynamics coupled with the set of classical river morphodynamic equations. The vegetation model only accounts for logistic growth, local positive feedback through seeding and resprouting, and mortality by means of uprooting through flow shear stress. Due to the simplicity of the model, we can transform the set of equations into an eigenvalue problem and assess the stability of the linearized equations when slightly perturbated away from a spatially homogeneous solution. If we couple vegetation dynamics with a 1D morphodynamic framework, we observe that instability towards long sediment waves is possible due to competitive interaction between vegetation growth and mortality. Moreover, the domain in the parameter space where perturbations are amplified was found to be simply connected. Subsequently, we proceed to the analysis of vegetation dynamics coupled with a 2D morphodynamic framework, which can be used to evaluate instability towards alternate and multiple bars. It is found that two kinds of instabilities, which are discriminated mainly by the Froude number, occur in a connected domain in the parameter space. At lower Froude number, instability is mainly governed by sediment dynamics and leads to the formation of alternate and multiple bars while at higher Froude number instability is driven by vegetation dynamics, which only allows for alternate bars.


title: Time-dependent modeling of extended thin decretion disks of critically rotating stars

abstract: During their evolution massive stars can reach the phase of critical rotation when a further increase in rotational speed is no longer possible. Direct centrifugal ejection from a critically or near-critically rotating surface forms a gaseous equatorial decretion disk. Anomalous viscosity provides the efficient mechanism for transporting the angular momentum outwards. The outer part of the disk can extend up to a very large distance from the parent star. We study the evolution of density, radial and azimuthal velocity, and angular momentum loss rate of equatorial decretion disks out to very distant regions. We investigate how the physical characteristics of the disk depend on the distribution of temperature and viscosity. We calculated stationary models using the Newton-Raphson method. For time-dependent hydrodynamic modeling we developed the numerical code based on an explicit finite difference scheme on an Eulerian grid including full Navier-Stokes shear viscosity. The sonic point distance and the maximum angular momentum loss rate strongly depend on the temperature profile and are almost independent of viscosity. The rotational velocity at large radii rapidly drops accordingly to temperature and viscosity distribution. The total amount of disk mass and the disk angular momentum increase with decreasing temperature and viscosity. The time-dependent one-dimensional models basically confirm the results obtained in the stationary models as well as the assumptions of the analytical approximations. Including full Navier-Stokes viscosity we systematically avoid the rotational velocity sign change at large radii. The unphysical drop of the rotational velocity and angular momentum loss at large radii (present in some models) can be avoided in the models with decreasing temperature and viscosity.