title: On the structure of finite groups isospectral to finite simple groups

abstract: Finite groups are said to be isospectral if they have the same sets of element orders. A finite nonabelian simple group $L$ is said to be almost recognizable by spectrum if every finite group isospectral to $L$ is an almost simple group with socle isomorphic to $L$. It is known that all finite simple sporadic, alternating and exceptional groups of Lie type, except $J_2$, $A_6$, $A_{10}$ and $^3D_4(2)$, are almost recognizable by spectrum. The present paper is the final step in the proof of the following conjecture due to V.D. Mazurov: there exists a positive integer $d_0$ such that every finite simple classical group of dimension larger than $d_0$ is almost recognizable by spectrum. Namely, we prove that a nonabelian composition factor of a~finite group isospectral to a finite simple symplectic or orthogonal group $L$ of dimension at least 10, is either isomorphic to $L$ or not a group of Lie type in the same characteristic as $L$, and combining this result with earlier work, we deduce that Mazurov's conjecture holds with $d_0=60$.


title: Scattering amplitudes in gauge theories with and without supersymmetry

abstract: This thesis aims at providing better understanding of the perturbative expansion of gauge theories with and without supersymmetry. At tree level, the BCFW recursion relations are analyzed with respect to their validity for general off-shell objects in Yang-Mills theory, which is a significant step away from their established zone of applicability. Unphysical poles constitute a new potential problem in addition to the boundary behavior issue, common to the on-shell case as well. For an infinite family of massive fermion currents, both obstacles are shown to be avoided under the certain conditions, which provides a natural recursion relation. At one loop, scattering amplitudes can be calculated from unitarity cuts through their expansion into known scalar integrals with free coefficients. A powerful method to obtain these coefficients, namely spinor integration, is discussed and rederived in a somewhat novel form. It is then used to compute analytically the infinite series of one-loop gluon amplitudes in N = 1 super-Yang-Mills theory with exactly three negative helicities. The final part of this thesis concerns the intriguing relationship between gluon and graviton scattering amplitudes, which involves a beautiful duality between the color and kinematic content of gauge theories. This BCJ duality is extended to include particles in the fundamental representation of the gauge group, which is shown to relieve the restriction of the BCJ construction to factorizable gravities and thus give access to amplitudes in generic (super-)gravity theories.


title: Nonstandard estimates for a class of 1D dispersive equations and applications to linearized water waves

abstract: In this work, we obtain decay bounds for a class of ID dispersive equations that includes the linearized water wave. These decay bounds display a surprising growth factor, which we show is sharp, The proofs rely on careful analysis of certain oscillatory integral operators. In addition, these results have applications to the linearized water wave operator for low regularity data.


title: Localization and magnetism of the resonant impurity states in Ti doped PbTe

abstract: The problem of localization of the resonant impurity states is discussed for an illustrative example of Titanium doped Lead Telluride. Electronic structure of PbTe:Ti is studied using first principles methods, densities of states and Bloch spectral functions are analyzed. We show that Ti creates resonant states in the conduction band of PbTe, however, spectral functions of the system strongly suggest localization of these states and show poor hybridization with PbTe electronic structure. The contrast between results presented here and previously reported spectral functions for PbTe:Tl correlate very well with the different effect of those impurities on thermopower (S) of PbTe, which is large increase is S for PbTe:Tl and almost no effect on S for PbTe:Ti. Moreover, magnetic properties of the system are studied and formation of magnetic moments on Ti atoms is found, both for ordered (ferromagnetic) and disordered (paramagnetic-like) phases, showing that PbTe:Ti can be a magnetic semiconductor.


title: Absence of a transport signature of spin-orbit coupling in graphene with indium adatoms

abstract: Enhancement of the spin-orbit coupling in graphene may lead to various topological phenomena and also find applications in spintronics. Adatom absorption has been proposed as an effective way to achieve the goal. In particular, great hope has been held for indium in strengthening the spin-orbit coupling and realizing the quantum spin Hall effect. To search for evidence of the spin-orbit coupling in graphene absorbed with indium adatoms, we carry out extensive transport measurements, i.e., weak localization magnetoresistance, quantum Hall effect and non-local spin Hall effect. No signature of the spin-orbit coupling is found. Possible explanations are discussed.


title: On Central Skew Armendariz rings

abstract: For a ring endomirphism, we introduce the central skew Armendariz rings, which are a generalization of skew Armendariz rings and central Armendariz rings, and investigate their properties.


title: Unconventional magnetism in the layered oxide LaSrRhO$_4$

abstract: We have prepared polycrystalline samples of LaSrRh$_{1-x}$Ga$_x$O$_4$ and LaSr$_{1-x}$Ca$_x$RhO$_4$,and have measured the x-ray diffraction, resistivity, Seebeck coefficient, magnetization and electron spin resonance in order to evaluate their electronic states. The energy gap evaluated from the resistivity and the Seebeck coefficient systematically changes with the Ga concentration, and suggests that the system changes from a small polaron insulator to a band insulator. We find that all the samples show Curie-Weiss-like susceptibility with a small Weiss temperature of the order of 1 K, which is seriously incompatible with the collective wisdom that a trivalent rhodium ion is nonmagnetic. We have determined the $g$ factor to be $g$=2.3 from the electron spin resonance, and the spin number to be $S$=1 from the magnetization-field curves by fitting with a modified Brillouin function. The fraction of the $S$=1 spins is 2--5%, which depends on the degree of disorder in the La/Sr/Ca-site, which implies that disorder near the apical oxygen is related to the magnetism of this system. A possible origin for the magnetic Rh$^{3+}$ ions is discussed.


title: The sorting index on colored permutations and even-signed permutations

abstract: We define a new statistic $\mathsf{sor}$ on the set of colored permutations $\mathsf{G}_{r,n}$ and prove that it has the same distribution as the length function. For the set of restricted colored permutations corresponding to the arrangements of $n$ non-attacking rooks on a fixed Ferrers shape we show that the following two sequences of set-valued statistics are joint equidistributed: $(\ell,\mathsf{Rmil}^0,\mathsf{Rmil}^1,...,\mathsf{Rmil}^{r-1}$, $\mathsf{Lmil}^0,\mathsf{Lmil}^1,...,\mathsf{Lmil}^{r-1}$, $\mathsf{Lmal}^0,\mathsf{Lmal}^1,...,\mathsf{Lmal}^{r-1}$, $\mathsf{Lmap}^0,\mathsf{Lmap}^1,...,\mathsf{Lmap}^{r-1})$ and $(\mathsf{sor},\mathsf{Cyc}^0,\mathsf{Cyc}^{r-1},...,\mathsf{Cyc}^{1}$, $\mathsf{Lmic}^0,\mathsf{Lmic}^{r-1},...,\mathsf{Lmic}^{1}$, $\mathsf{Lmal}^0,\mathsf{Lmal}^1,...,\mathsf{Lmal}^{r-1}$, $\mathsf{Lmap}^0,\mathsf{Lmap}^1,...,\mathsf{Lmap}^{r-1})$. Analogous results are also obtained for Coxeter group of type $D$. Our results extend recent results of Petersen, Chen-Gong-Guo and Poznanovi\'{c}.


title: On quasi-ergodic distribution for one-dimensional diffusions

abstract: In this paper, we study quasi-ergodicity for one-dimensional diffusion $X$ killed at 0, when 0 is an exit boundary and $+\infty$ is an entrance boundary. Using the spectral theory tool, we show that if the killed semigroup is intrinsically ultracontractive, then there exists a unique quasi-ergodic distribution for $X$. An example is given to illustrate the result. Moreover, the ultracontractivity of the killed semigroup is also studied.


title: The combinatorics of real double Hurwitz numbers with real positive branch points

abstract: We investigate the combinatorics of real double Hurwitz numbers with real positive branch points using the symmetric group. Our main focus is twofold. First, we prove correspondence theorems relating these numbers to counts of tropical real covers and study the structure of real double Hurwitz numbers with the help of the tropical count. Second, we express the numbers as counts of paths in a subgraph of the Cayley graph of the symmetric group. By restricting to real double Hurwitz numbers with real positive branch points, we obtain a concise translation of the counting problem in terms of tuples of elements of the symmetric group that enables us to uncover the beautiful combinatorics of these numbers both in tropical geometry and in the Cayley graph.


title: Robustness of Controlled Quantum Dynamics

abstract: Control of multi-level quantum systems is sensitive to implementation errors in the control field and uncertainties associated with system Hamiltonian parameters. A small variation in the control field spectrum or the system Hamiltonian can cause an otherwise optimal field to deviate from controlling desired quantum state transitions and reaching a particular objective. An accurate analysis of robustness is thus essential in understanding and achieving model-based quantum control, such as in control of chemical reactions based on ab initio or experimental estimates of the molecular Hamiltonian. In this paper, theoretical foundations for quantum control robustness analysis are presented from both a distributional perspective - in terms of moments of the transition amplitude, interferences, and transition probability - and a worst-case perspective. Based on this theory, analytical expressions and a computationally efficient method for determining the robustness of coherently controlled quantum dynamics are derived. The robustness analysis reveals that there generally exists a set of control pathways that are more resistant to destructive interferences in the presence of control field and system parameter uncertainty. These robust pathways interfere and combine to yield a relatively accurate transition amplitude and high transition probability when uncertainty is present.


title: Hydrogen assisted growth of high quality epitaxial graphene on the C-face of 4H-SiC

abstract: We demonstrate hydrogen assisted growth of high quality epitaxial graphene on the C-face of 4H-SiC. Compared with the conventional thermal decomposition technique, the size of the growth domain by this method is substantially increased and the thickness variation is reduced. Based on the morphology of epitaxial graphene, the role of hydrogen is revealed. It is found that hydrogen acts as a carbon etchant. It suppresses the defect formation and nucleation of graphene. It also improves the kinetics of carbon atoms via hydrocarbon species. These effects lead to increase of the domain size and the structure quality. The consequent capping effect results in smooth surface morphology and suppression of multilayer growth. Our method provides a viable route to fine tune the growth kinetics of epitaxial graphene on SiC.


title: Workflow Partitioning and Deployment on the Cloud using Orchestra

abstract: Orchestrating service-oriented workflows is typically based on a design model that routes both data and control through a single point - the centralised workflow engine. This causes scalability problems that include the unnecessary consumption of the network bandwidth, high latency in transmitting data between the services, and performance bottlenecks. These problems are highly prominent when orchestrating workflows that are composed from services dispersed across distant geographical locations. This paper presents a novel workflow partitioning approach, which attempts to improve the scalability of orchestrating large-scale workflows. It permits the workflow computation to be moved towards the services providing the data in order to garner optimal performance results. This is achieved by decomposing the workflow into smaller sub workflows for parallel execution, and determining the most appropriate network locations to which these sub workflows are transmitted and subsequently executed. This paper demonstrates the efficiency of our approach using a set of experimental workflows that are orchestrated over Amazon EC2 and across several geographic network regions.


title: Ideal structure of the C*-algebra of Thompson group T

abstract: In a recent paper Uffe Haagerup and Kristian Knudsen Olesen show that for Richard Thompson's group $T$, if there exists a finite set $H$ which can be decomposed as disjoint union of sets $H_1$ and $H_2$ with $\sum_{g\in H_1}\pi(g)=\sum_{h\in H_2}\pi(h)$ and such that the closed ideal generated by $\sum_{g\in H_1}\lambda(g)-\sum_{h\in H_2}\lambda(h)$ coincides with $C^*_\lambda(T)$, then the Richard Thompson group $F$ is not amenable. In particular, if $C_{\lambda}^*(T)$ is simple then $F$ is not amenable. Here we prove the converse, namely, if $F$ is not amenable then we can find two sets $H_1$ and $H_2$ with the above properties. The only currently available tool for proving simplicity of group $C^*$-algebra is Power's condition. We show that it fails for $C_{\lambda}^*(T)$ and present an apparent weakening of that condition which could potentially be used for various new groups $H$ to show the simplicity of $C_{\lambda}^*(H)$. While we use our weakening in the proof of the first result, we also show that the new condition is still too strong to be used to show the simplicity of $C_{\lambda}^*(T)$. Along the way, we give a new application of the Ping-Pong Lemma to find free groups as subgroups in groups of homeomorphisms of the circle generated by elements with rational rotation number.


title: Anisotropy in the Hubble constant as modeled by density gradients

abstract: The all-sky maps of the observed variation of the Hubble constant can be reproduced from a theoretical point of view by introducing an intergalactic plasma with a variable number density of electrons. The observed averaged value and variance of the Hubble constant are reproduced by adopting a rim model, an auto-gravitating model, and a Voronoi diagrams model as the backbone for an auto-gravitating medium. We also analyze an astronomer's model based on the 3D spatial distribution of galaxies as given by the 2MASS Redshift Survey and an auto-gravitating Lane--Emden ($n=5$) profile of the electrons. The simulation which involves the Voronoi diagrams is done in a cubic box with sides of 100 Mpc. The simulation which involves the 2MASS covers the range of redshift smaller than 0.05.


title: Radio-Gamma-ray connection and spectral evolution in 4C +49.22 (S4 1150+49): the Fermi, Swift and Planck view

abstract: The Large Area Telescope on board the Fermi Gamma-ray Space Telescope detected a strong gamma-ray flare on 2011 May 15 from a source identified as 4C 49.22, a flat spectrum radio quasar also known as S4 1150+49. This blazar, characterised by a prominent radio-optical-X-ray jet, was in a low gamma-ray activity state during the first years of Fermi observations. Simultaneous observations during the quiescent, outburst and post-flare gamma-ray states were obtained by Swift, Planck and optical-IR-radio telescopes (INAOE, Catalina CSS, VLBA, Metsahovi). The flare is observed from microwave to X-ray bands with correlated variability and the Fermi, Swift and Planck data for this FSRQ show some features more typical of BL Lac objects, like the synchrotron peak in the optical band that outshines the thermal blue-bump emission, and the X-ray spectral softening. Multi-epoch VLBA observations show the ejection of a new component close in time with the GeV gamma-ray flare. The radio-to-gamma-ray spectral energy distribution is modeled and fitted successfully for the outburst and the post-flare epochs using either a single flaring blob with two emission processes (synchrotron self Compton, and external-radiation Compton), and a two-zone model with SSC-only mechanism.


title: Boundedness of large-time solutions to a chemotaxis model with nonlocal and semilinear flux

abstract: A semilinear version of parabolic-elliptic Keller-Segel system with the \emph{critical} nonlocal diffusion is considered in one space dimension. We show boundedness of weak solutions under very general conditions on our semilinearity. It can degenerate, but has to provide a stronger dissipation for large values of a solution than in the critical linear case or we need to assume certain (explicit) data smallness. Moreover, when one considers a logistic term with a parameter $r$, we obtain our results even for diffusions slightly weaker than the critical linear one and for arbitrarily large initial datum, provided $r >1$. For a mild logistic dampening, we can improve the smallness condition on the initial datum up to $\sim \frac{1}{1-r}$.


title: Domain of attraction of quasi-stationary distribution for one-dimensional diffusions

abstract: We study quasi-stationarity for one-dimensional diffusions killed at 0, when 0 is a regular boundary and $+\infty$ is an entrance boundary. We give a necessary and sufficient condition for the existence of exactly one quasi-stationary distribution, and we also show that this distribution attracts all initial distributions.


title: Capacity Region of MISO Broadcast Channel for Simultaneous Wireless Information and Power Transfer

abstract: This paper studies a multiple-input single-output (MISO) broadcast channel (BC) featuring simultaneous wireless information and power transfer (SWIPT), where a multi-antenna access point (AP) delivers both information and energy via radio signals to multiple single-antenna receivers simultaneously, and each receiver implements either information decoding (ID) or energy harvesting (EH). In particular, pseudo-random sequences that are {\it a priori} known and therefore can be cancelled at each ID receiver is used as the energy signals, and the information-theoretically optimal dirty paper coding (DPC) is employed for the information transmission. We characterize the capacity region for ID receivers under given energy requirements for EH receivers, by solving a sequence of weighted sum-rate (WSR) maximization (WSRMax) problems subject to a maximum sum-power constraint for the AP, and a set of minimum harvested power constraints for individual EH receivers. The problem corresponds to a new form of WSRMax problem in MISO-BC with combined maximum and minimum linear transmit covariance constraints (MaxLTCCs and MinLTCCs), which differs from the celebrated capacity region characterization problem for MISO-BC under a set of MaxLTCCs only and is challenging to solve. By extending the general BC-multiple access channel (MAC) duality, which is only applicable to WSRMax problems with MaxLTCCs, and applying the ellipsoid method, we propose an efficient algorithm to solve this problem globally optimally. Furthermore, we also propose two suboptimal algorithms with lower complexity by assuming that the information and energy signals are designed separately. Finally, numerical results are provided to validate our proposed algorithms.


title: The volume of random polytopes circumscribed around a convex body

abstract: Let $K$ be a convex body in $\mathbb{R}^d$ which slides freely in a ball. Let $K^{(n)}$ denote the intersection of $n$ closed half-spaces containing $K$ whose bounding hyperplanes are independent and identically distributed according to a certain prescribed probability distribution. We prove an asymptotic formula for the expectation of the difference of the volumes of $K^{(n)}$ and $K$, and an asymptotic upper bound on the variance of the volume of $K^{(n)}$. We achieve these results by first proving similar statements for weighted mean width approximations of convex bodies that admit a rolling ball by inscribed random polytopes and then by polarizing these results.


title: Comment on the numerical measurements of the magnetohydrodynamic turbulence spectrum by A. Beresnyak (Phys. Rev. Lett. 106 (2011) 075001; MNRAS 422 (2012) 3495; ApJ 784 (2014) L20)

abstract: The inertial-interval energy spectrum of strong magnetohydrodynamic (MHD) turbulence with a uniform background magnetic field was observed numerically to be close to $k^{-3/2}$ by a number of independent groups. A dissenting opinion has been voiced by Beresnyak, A. 2011, PRL, 106, 075001-. 2012, MNRAS, 422, 3495-. 2014, ApJ, 784, L20 that the spectral scaling is close to $k^{-5/3}$. The conclusions of these papers are however incorrect as they are based on numerical simulations that are drastically unresolved, so that the discrete numerical scheme does not approximate the physical solution at the scales where the measurements are performed. These results have been rebutted in our more detailed papers Perez, J. C., Mason, J., Boldyrev, S., & Cattaneo, F. 2012, PRX, 2, 041005-. 2014, ApJL, 793, L13; here, by popular demand, we present a brief and simple explanation of our major criticism of Beresnyak's work.


title: Hot spot manifestation in eclipsing dwarf nova HT Cassiopeiae

abstract: We report the detection of the hot spot in light curves of the eclipsing dwarf nova HT Cassiopeiae during its superoutburst in 2010 November. Analysis of eight reconstructed light curves of the hot spot eclipses showed directly that the brightness of the hot spot was changing significantly during the superoutburst. Thereby, detected hot spot manifestation in HT Cas is the newest observational evidence for the EMT model for dwarf novae.


title: $p$-vanishing conjugacy classes of symmetric groups

abstract: For a prime $p$, we say that a conjugacy class of a finite group $G$ is $p$-vanishing if every irreducible character of $G$ of degree divisible by $p$ takes value 0 on that conjugacy class. In this paper we completely classify 2-vanishing and 3-vanishing conjugacy classes for the symmetric group and do some work in the classification of $p$-vanishing conjugacy classes of the symmetric group for $p\geq 5$. This answers a question by Navarro for $p=2$ and $p=3$ and partly answers it for $p\geq 5$.


title: Sequential Monte Carlo samplers for semilinear inverse problems and application to magnetoencephalography

abstract: We discuss the use of a recent class of sequential Monte Carlo methods for solving inverse problems characterized by a semi-linear structure, i.e. where the data depend linearly on a subset of variables and nonlinearly on the remaining ones. In this type of problems, under proper Gaussian assumptions one can marginalize the linear variables. This means that the Monte Carlo procedure needs only to be applied to the nonlinear variables, while the linear ones can be treated analytically; as a result, the Monte Carlo variance and/or the computational cost decrease. We use this approach to solve the inverse problem of magnetoencephalography, with a multi-dipole model for the sources. Here, data depend nonlinearly on the number of sources and their locations, and depend linearly on their current vectors. The semi-analytic approach enables us to estimate the number of dipoles and their location from a whole time-series, rather than a single time point, while keeping a low computational cost.


title: The principal series of $p$-adic groups with disconnected centre

abstract: Let G be a split connected reductive group over a local non-archimedean field. We classify all irreducible complex G-representations in the principal series, irrespective of the (dis)connectedness of the centre of G. This leads to a local Langlands correspondence for principal series representations, which satisfies all expected properties. We also prove that the ABPS conjecture about the geometric structure of Bernstein components is valid throughout the principal series of G.


title: Resolving the shocked gas in HH54 with Herschel: CO line mapping at high spatial and spectral resolution

abstract: The HH54 shock is a Herbig-Haro object, located in the nearby Chamaeleon II cloud. Observed CO line profiles are due to a complex distribution in density, temperature, velocity, and geometry. Resolving the HH54 shock wave in the far-infrared cooling lines of CO constrain the kinematics, morphology, and physical conditions of the shocked region. We used the PACS and SPIRE instruments on board the Herschel space observatory to map the full FIR spectrum in a region covering the HH54 shock wave. Complementary Herschel-HIFI, APEX, and Spitzer data are used in the analysis as well. The observed features in the line profiles are reproduced using a 3D radiative transfer model of a bow-shock, constructed with the Line Modeling Engine code (LIME). The FIR emission is confined to the HH54 region and a coherent displacement of the location of the emission maximum of CO with increasing J is observed. The peak positions of the high-J CO lines are shifted upstream from the lower J CO lines and coincide with the position of the spectral feature identified previously in CO(10-9) profiles with HIFI. This indicates a hotter molecular component in the upstream gas with distinct dynamics. The coherent displacement with increasing J for CO is consistent with a scenario where IRAS12500-7658 is the exciting source of the flow, and the 180 K bow-shock is accompanied by a hot (800 K) molecular component located upstream from the apex of the shock and blueshifted by -7 km s$^{-1}$. The spatial proximity of this knot to the peaks of the atomic fine-structure emission lines observed with Spitzer and PACS ([OI]63, 145 $\mu$m) suggests that it may be associated with the dissociative shock as the jet impacts slower moving gas in the HH54 bow-shock.


title: Efficient high-quality motion planning by fast all-pairs r-nearest-neighbors

abstract: Sampling-based motion-planning algorithms typically rely on nearest-neighbor (NN) queries when constructing a roadmap. Recent results suggest that in various settings NN queries may be the computational bottleneck of such algorithms. Moreover, in several asymptotically-optimal algorithms these NN queries are of a specific form: Given a set of points and a radius r report all pairs of points whose distance is at most r. This calls for an application-specific NN data structure tailored to efficiently answering this type of queries. Randomly transformed grids (RTG) were recently proposed by Aiger et al. as a tool to answer such queries and have been shown to outperform common implementations of NN data structures in this context. In this work we employ RTG for sampling-based motion-planning algorithms and describe an efficient implementation of the approach. We show that for motion-planning, RTG allow for faster convergence to high-quality solutions when compared with existing NN data structures. Additionally, RTG enable significantly shorter construction times for batched-PRM variants; specifically, we demonstrate a speedup by a factor of two to three for some scenarios.


title: Results and Future Plans on Central Exclusive Production with the LHCb Detector

abstract: The LHCb detector and LHC running conditions are ideally suited to measuring central exclusive production. Several recent measurements of exclusive dimuon, single and double charmonia are reviewed. The potential for future measurements across a broad range of physics channels is discussed.


title: Graph properties of graph associahedra

abstract: A graph associahedron is a simple polytope whose face lattice encodes the nested structure of the connected subgraphs of a given graph. In this paper, we study certain graph properties of the 1-skeleta of graph associahedra, such as their diameter and their Hamiltonicity. Our results extend known results for the classical associahedra (path associahedra) and permutahedra (complete graph associahedra). We also discuss partial extensions to the family of nestohedra.


title: Muon Anomalous Magnetic Moment in a $SU(4) \otimes U(1)_N$ Model

abstract: We study the muon anomalous magnetic moment in an electroweak model based on the gauge symmetry $SU(4)_L\otimes U(1)_N$, which has right-handed neutrinos in its spectrum, and no flavor changing neutral currents at tree level. We discuss relevant collider and electroweak constraints on the model, and derive the most stringent upper bounds on the scale of symmetry breaking based on the corrections to the muon magnetic moment. We conclude that a scale of symmetry breaking of around 2TeV might explain the muon magnetic moment anomaly. In case the anomaly is otherwise resolved, using the current and projected sensitive of g-2 experiments, we rule out scales of symmetry breaking smaller than 3.5TeV (5.5TeV) at $1\sigma$ level.


title: PoisFFT - A Free Parallel Fast Poisson Solver

abstract: A fast Poisson solver software package PoisFFT is presented. It is available as a free software licensed under the GNU GPL license version 3. The package uses the fast Fourier transform to directly solve the Poisson equation on a uniform orthogonal grid. It can solve the pseudo-spectral approximation and the second order finite difference approximation of the continuous solution. The paper reviews the mathematical methods for the fast Poisson solver and discusses the software implementation and parallelization. The use of PoisFFT in an incompressible flow solver is also demonstrated.


title: Validation of the magnetic energy vs. helicity scaling in solar magnetic structures

abstract: We assess the validity of the free magnetic energy - relative magnetic helicity diagram for solar magnetic structures. We used two different methods of calculating the free magnetic energy and the relative magnetic helicity budgets: a classical, volume-calculation nonlinear force-free (NLFF) method applied to finite coronal magnetic structures and a surface-calculation NLFF derivation that relies on a single photospheric or chromospheric vector magnetogram. Both methods were applied to two different data sets, namely synthetic active-region cases obtained by three-dimensional magneto-hydrodynamic (MHD) simulations and observed active-region cases, which include both eruptive and noneruptive magnetic structures. The derived energy--helicity diagram shows a consistent monotonic scaling between relative helicity and free energy with a scaling index 0.84$\pm$0.05 for both data sets and calculation methods. It also confirms the segregation between noneruptive and eruptive active regions and the existence of thresholds in both free energy and relative helicity for active regions to enter eruptive territory. We consider the previously reported energy-helicity diagram of solar magnetic structures as adequately validated and envision a significant role of the uncovered scaling in future studies of solar magnetism.


title: Some identities involving Prouhet-Thue-Morse sequence and its relatives

abstract: Let $s_{k}(n)$ denote the sum of digits of an integer $n$ in base $k$. Motivated by certain identities of Nieto, and Bateman and Bradley involving sums of the form $\sum_{i=0}^{2^{n}-1}(-1)^{s_{2}(i)}(x+i)^{m}$ for $m=n$ and $m=n+1$, we consider the sequence of polynomials \begin{equation*} f_{m,n}^{\mathbf u}(x)=\sum_{i=0}^{k^{n}-1}\zeta_{k}^{s_{k}(i)}(x+{\mathbf u}(i))^{m}. \end{equation*} defined for sequences ${\bf u}(i)$ satisfying a certain recurrence relation. We prove that computing these polynomials is essentially equivalent with computing their constant term and we find an explicit formula for this number. This allows us to prove several interesting identities involving sums of binary digits. We also prove some related results which are of independent interests and can be seen as further generalizations of certain sums involving Prouhet-Thue-Morse sequence.


title: Scaling analysis of time series of daily prices from stock markets of transitional economies in the Western Balkans

abstract: In this paper we have analyzed scaling properties of time series of stock market indices (SMIs) of developing economies of Western Balkans, and have compared the results we have obtained with the results from more developed economies. We have used three different techniques of data analysis to obtain and verify our findings: Detrended Fluctuation Analysis (DFA) method, Detrended Moving Average (DMA) method, and Wavelet Transformation (WT) analysis. We have found scaling behavior in all SMI data sets that we have analyzed. The scaling of our SMI series changes from long-range correlated to slightly anti-correlated behavior with the change in growth or maturity of the economy the stock market is embedded in. We also report the presence of effects of potential periodic-like influences on the SMI data that we have analyzed. One such influence is visible in all our SMI series, and appears at a period $T_{p}\approx 90$ days. We propose that the existence of various periodic-like influences on SMI data may partially explain the observed difference in types of correlated behavior of corresponding scaling functions.


title: Extracting femtoscopic radii in the presence of significant additional correlation sources

abstract: The Large Hadron Collider has provided large amounts of data on collisions of small systems, such as proton-proton and proton-lead at unprecedented collision energies. Their space-time size and structure can be inferred from the measurement of the femtoscopic correlations for pairs of identical particles. The analysis is complicated by the presence of significant additional sources of two-particle correlations, which influence the correlation function in the region of the femtoscopic effect. In this work we use p-Pb events generated in a model that includes such additional correlation sources to characterize them and propose a robust method of taking them into account in the extraction of the femtoscopic information.


title: Properties of optically selected BL Lac candidates from the SDSS

abstract: \textbf{Context.} Deep optical surveys open the avenue for find large numbers of BL Lac objects that are hard to identify because they lack the unique properties classifying them as such. While radio or X-ray surveys typically reveal dozens of sources, recent compilations based on optical criteria alone have increased the number of BL Lac candidates considerably. However, these compilations are subject to biases and may contain a substantial number of contaminating sources. \textbf{Aims.} In this paper we extend our analysis of 182 optically selected BL Lac object candidates from the SDSS with respect to an earlier study. The main goal is to determine the number of bona fide BL Lac objects in this sample. \textbf{Methods.} We examine their variability characteristics, determine their broad-band radio-UV SEDs, and search for the presence of a host galaxy. In addition we present new optical spectra for 27 targets with improved S/N with respect to the SDSS spectra. \textbf{Results.} At least 59% of our targets have shown variability between SDSS DR2 and our observations by more than 0.1-0.27 mag de- pending on the telescope used. A host galaxy was detected in 36% of our targets. The host galaxy type and luminosities are consistent with earlier studies of BL Lac host galaxies. Simple fits to broad-band SEDS for 104 targets of our sample derived synchrotron peak frequencies between $13.5 \leq \mathrm{log}_{10}(\nu_{\mathrm{peak}}) \leq 16$ with a peak at $\mathrm{log}_{10} \sim 14.5$. Our new optical spectra do not reveal any new redshift for any of our objects. Thus the sample contains a large number of bona fide BL Lac objects and seems to contain a substantial fraction of intermediate-frequency peaked BL Lacs.


title: Borderline gradient continuity of minima

abstract: The gradient of any local minimiser of functionals of the type $$ w \mapsto \int_\Omega f(x,w,Dw)\,dx+\int_\Omega w\mu\,dx, $$ where $f$ has $p$-growth, $p>1$, and $\Omega \subset \mathbb R^n$, is continuous provided the optimal Lorentz space condition $\mu \in L(n,1)$ is satisfied and $x\to f(x, \cdot)$ is suitably Dini-continuous.


title: The number of the maximal triangle-free graphs

abstract: Paul Erd\H{o}s suggested the following problem: Determine or estimate the number of maximal triangle-free graphs on $n$ vertices. Here we show that the number of maximal triangle-free graphs is at most $2^{n^2/8+o(n^2)}$, which matches the previously known lower bound. Our proof uses among others the Ruzsa-Szemer\'{e}di triangle removal lemma, and recent results on characterizing of the structure of independent sets in hypergraphs.


title: Self-Gravitating Relativistic Models of Fermions with Anisotropy and Cutoff Energy in their Distribution Function

abstract: In this paper we study the equilibrium configurations of anisotropic self-gravitating fermions, by extending to general relativity the solutions obtained in a previous paper. This treatment also generalizes to anisotropic systems the relativistic self-gravitating Fermi gas model, by considering different degrees of anisotropy. We discuss some important characteristics of the models and the obtained density profiles, and generalize the relation between the anisotropy and the mass of particles in the relativistic regime. These relativistic models may also be applied to the study of superdense neutron stars with anisotropic pressure or super-Chandrasekhar white dwarfs generated by the presence of a magnetic field.


title: Random Access Protocols with Collision Resolution in a Noncoherent Setting

abstract: Wireless systems are increasingly used for Machine-Type Communication (MTC), where the users sporadically send very short messages. In such a setting, the overhead imposed by channel estimation is substantial, thereby demanding noncoherent communication. In this paper we consider a noncoherent setup in which users randomly access the medium to send short messages to a common receiver. We propose a transmission scheme based on Gabor frames, where each user has a dedicated codebook of M possible codewords, while the codebook simultaneously serves as an ID for the user. The scheme is used as a basis for a simple protocol for collision resolution.


title: Measurability of open orbits in flag supermanifolds

abstract: Let $G_\mathbb{R}$ be a Lie group and $G$ its complexification. An open $G_\mathbb{R}$-orbit in a $G$-flag manifold is measurable whenever it carries a $G_\mathbb{R}$-invariant volume element. In this paper the notion of measurability is generalized to the supersymmetric setting and a classification of measurable open orbits in flag supermanifolds is given.


title: Efimov like phase of a three stranded DNA (Efimov-DNA) and the renormalization group limit cycle

abstract: A three-stranded DNA with short range base pairings only is known to exhibit a classical analog of the quantum Efimov effect, viz., a three chain bound state at the two chain melting point where no two are bound. By using a non-perturbative renormalization group method for a rigid duplex DNA and a flexible third strand, with base pairings and strand exchange, we show that the Efimov-DNA is associated with a limit cycle type behavior of the flow of an effective three chain interaction. The analysis also shows that thermally generated bubbles play an essential role in producing the effect. A toy model for the flow equations shows the limit cycle in an extended three dimensional parameter space of the two-chain coupling and a complex three chain interaction.


title: Electron thermal conduction as a possible mechanism to make the inner heliosheath thnner

abstract: We show that the electron thermal conductivity may strongly affect the heliosheath plasma flow and the global pattern of the solar wind (SW) interaction with the local interstellar medium (LISM). In particular, it leads to strong reduction of the inner heliosheath thickness that makes possible to explain (qualitatively) why Voyager 1 (V1) has crossed the heliopause at unexpectedly small heliocentric distance of 122 AU. To estimate the effect of thermal conductivity we consider a limiting case when thermal conduction is very effective. To do that we assume the plasma flow in the entire heliosphere is nearly isothermal. Due to this effect, the heliospheric distance of the termination shock has increased by about 15 AU in V1 direction compared to the adiabatic case with gamma = 5/3. The heliospheric distance of the heliopause has decreased by about 27 AU. As a result, the thickness of the inner heliosheath in the model has decreased by about 42 AU and become equal to 32 AU.


title: Collaborative sparse regression using spatially correlated supports - Application to hyperspectral unmixing

abstract: This paper presents a new Bayesian collaborative sparse regression method for linear unmixing of hyperspectral images. Our contribution is twofold; first, we propose a new Bayesian model for structured sparse regression in which the supports of the sparse abundance vectors are a priori spatially correlated across pixels (i.e., materials are spatially organised rather than randomly distributed at a pixel level). This prior information is encoded in the model through a truncated multivariate Ising Markov random field, which also takes into consideration the facts that pixels cannot be empty (i.e, there is at least one material present in each pixel), and that different materials may exhibit different degrees of spatial regularity. Secondly, we propose an advanced Markov chain Monte Carlo algorithm to estimate the posterior probabilities that materials are present or absent in each pixel, and, conditionally to the maximum marginal a posteriori configuration of the support, compute the MMSE estimates of the abundance vectors. A remarkable property of this algorithm is that it self-adjusts the values of the parameters of the Markov random field, thus relieving practitioners from setting regularisation parameters by cross-validation. The performance of the proposed methodology is finally demonstrated through a series of experiments with synthetic and real data and comparisons with other algorithms from the literature.


title: A primal-dual mimetic finite element scheme for the rotating shallow water equations on polygonal spherical meshes

abstract: A new numerical method is presented for solving the rotating shallow water equations on a rotating sphere using quasi-uniform polygonal meshes. The method uses special families of finite element function spaces to mimic key mathematical properties of the continuous equations and thereby capture several desirable physical properties related to balance and conservation. The method relies on two novel features. The first is the use of compound finite elements to provide suitable finite element spaces on general polygonal meshes. The second is the use of dual finite element spaces on the dual of the original mesh, along with suitably defined discrete Hodge star operators to map between the primal and dual meshes, enabling the use of a finite volume scheme on the dual mesh to compute potential vorticity fluxes. The resulting method has the same mimetic properties as a finite volume method presented previously, but is more accurate on a number of standard test cases.


title: VLT/UVES observations of peculiar alpha abundances in a sub-DLA at z ~ 1.8 towards the quasar B1101-26

abstract: We present a detailed analysis of chemical abundances in a sub-damped Lyman alpha absorber at z=1.839 towards the quasar B1101-26, based on a very-high-resolution (R ~ 75,000) and high-signal-to-noise (S/N >100) spectrum observed with the UV Visual Echelle spectrograph (UVES) installed on the ESO Very Large Telescope (VLT). The absorption line profiles are resolved into a maximum of eleven velocity components spanning a rest-frame velocity range of 200 km/s. Detected ions include CII, CIV, NII, OI, MgI, MgII, AlII, AlIII, SiII, SiIII, SiIV, FeII, and possibly SII. The total neutral hydrogen column density is log N(HI) = 19.48 +/- 0.01. From measurements of column densities and Doppler parameters we estimate element abundances of the above-given elements. The overall metallicity, as traced by [OI/HI], is -1.56 +/- 0.01. For the nitrogen-to-oxygen ratio we derive an upper limit of [NI/OI] < -0.65, which suggests a chemically young absorption line system. This is supported by a supersolar alpha/Fe ratio of [SiII/FeII] ~ 0.5. The most striking feature in the observed abundance pattern is an unusually high sulphur-to-oxygen ratio of 0.69 < [SII/OI] < 1.26. We calculate detailed photoionisation models for two subcomponents with Cloudy, and can rule out that ionisation effects alone are responsible for the high S/O ratio. We instead speculate that the high S/O ratio is caused by the combination of several effects, such as specific ionisation conditions in multi-phase gas, unusual relative abundances of heavy elements, and/or dust depletion in a local gas environment that is not well mixed and/or that might be related to star-formation activity in the host galaxy. We discuss the implications of our findings for the interpretation of alpha-element abundances in metal absorbers at high redshift.


title: Some experiments with Ramanujan-Nagell type Diophantine equations

abstract: Stiller proved that the Diophantine equation $x^2+119=15\cdot 2^{n}$ has exactly six solutions in positive integers. Motivated by this result we are interested in constructions of Diophantine equations of Ramanujan-Nagell type $x^2=Ak^{n}+B$ with many solutions. Here, $A,B\in\Z$ (thus $A, B$ are not necessarily positive) and $k\in\Z_{\geq 2}$ are given integers. In particular, we prove that for each $k$ there exists an infinite set $\cal{S}$ containing pairs of integers $(A, B)$ such that for each $(A,B)\in \cal{S}$ we have $\gcd(A,B)$ is square-free and the Diophantine equation $x^2=Ak^n+B$ has at least four solutions in positive integers. Moreover, we construct several Diophantine equations of the form $x^2=Ak^n+B$ with $k>2$, each containing five solutions in non-negative integers. %For example the equation $y^2=130\cdot 3^{n}+5550606$ has exactly five solutions with $n=0, 6, 11, 15, 16$. We also find new examples of equations $x^2=A2^{n}+B$ having six solutions in positive integers, e.g. the following Diophantine equations has exactly six solutions: \begin{equation*} \begin{array}{ll} x^2= 57\cdot 2^{n}+117440512 & n=0, 14, 16, 20, 24, 25, x^2= 165\cdot 2^{n}+26404 & n=0, 5, 7, 8, 10, 12. \end{array} \end{equation*} Moreover, based on an extensive numerical calculations we state several conjectures on the number of solutions of certain parametric families of the Diophantine equations of Ramanujan-Nagell type.


title: The Maximum k-Differential Coloring Problem

abstract: Given an $n$-vertex graph $G$ and two positive integers $d,k \in \mathbb{N}$, the ($d,kn$)-differential coloring problem asks for a coloring of the vertices of $G$ (if one exists) with distinct numbers from 1 to $kn$ (treated as \emph{colors}), such that the minimum difference between the two colors of any adjacent vertices is at least $d$. While it was known that the problem of determining whether a general graph is ($2,n$)-differential colorable is NP-complete, our main contribution is a complete characterization of bipartite, planar and outerplanar graphs that admit ($2,kn$)-differential colorings. For practical reasons, we consider also color ranges larger than $n$, i.e., $k > 1$. We show that it is NP-complete to determine whether a graph admits a ($3,2n$)-differential coloring. The same negative result holds for the ($\lfloor 2n/3 \rfloor, 2n$-differential coloring problem, even in the case where the input graph is planar.


title: Limit theorems of a two-phase quantum walk with one defect

abstract: We treat a position dependent quantum walk (QW) on the line which we assign two different time-evolution operators to positive and negative parts respectively. We call the model "the two-phase QW" here, which has been expected to be a mathematical model of the topological insulator. We obtain the stationary and time-averaged limit measures related to localization for the two-phase QW with one defect. This is the first result on localization for the two-phase QW. The analytical methods are mainly based on the splitted generating function of the solution for the eigenvalue problem, and the generating function of the weight of the passages of the model. In this paper, we call the methods "the splitted generating function method" and "the generating function method", respectively. The explicit expression of the stationary measure is asymmetric for the origin, and depends on the initial state and the choice of the parameters of the model. On the other hand, the time-averaged limit measure has a starting point symmetry and localization effect heavily depends on the initial state and the parameters of the model. Regardless of the strong effect of the initial state and the parameters, the time-averaged limit measure also suggests that localization can be always observed for our two-phase QW. Furthermore, our results imply that there is an interesting relation between the stationary and time-averaged limit measures when the parameters of the model have specific periodicities, which suggests that there is a possibility that we can analyze localization of the two-phase QW with one defect from the stationary measure.


title: A note on the Minimum Norm Point algorithm

abstract: We present a provably more efficient implementation of the Minimum Norm Point Algorithm conceived by Fujishige than the one presented in \cite{FUJI06}. The algorithm solves the minimization problem for a class of functions known as submodular. Many important functions, such as minimum cut in the graph, have the so called submodular property \cite{FUJI82}. It is known that the problem can also be efficiently solved in strongly polynomial time \cite{IWAT01}, however known theoretical bounds are far from being practical. We present an improved implementation of the algorithm, for which unfortunately no worst case bounds are know, but which performs very well in practice. With the modifications presented, the algorithm performs an order of magnitude faster for certain submodular functions.


title: Which spacetimes admit conformal compactifications?

abstract: We consider the future causal boundary as a tool to find obstructions to conformal extensions, the latter being a slight generalization to conformal compactifications.


title: A statistical noise model for a class of Physically Unclonable Functions

abstract: The interest in "Physically Unclonable Function"-devices has increased rapidly over the last few years, as they have several interesting properties for system security related applications like, for example, the management of cryptographic keys. Unfortunately, the output provided by these devices is noisy and needs to be corrected for these applications. Related error correcting mechanisms are typically constructed on the basis of an equal error probability for each output bit. This assumption does not hold for Physically Unclonable Functions, where varying error probabilities can be observed. This results in a generalized binomial distribution for the number of errors in the output. The intention of this paper is to discuss a novel Bayesian statistical model for the noise of an especially wide-spread class of Physically Unclonable Functions, which properly handles the varying output stability and also reflects the different noise behaviors observed in a collection of such devices. Furthermore, we compare several different methods for estimating the model parameters and apply the proposed model to concrete measurements obtained within the CODES research project in order to evaluate typical correction and stabilization approaches.


title: Baryons made of heavy quarks at the LHC

abstract: Baryons made of heavy quarks are extremely interesting and could be seen at the LHC.


title: Holographic Estimate of Electromagnetic Mass

abstract: Using the gauge/gravity duality, we calculate the electromagnetic contributions to hadron masses, where mass generates dynamically by strong QCD interactions. Based on the Sakai-Sugimoto model of holographic QCD we find that the electromagnetic mass of proton is $0.48~{\rm MeV}$ larger than that of neutron, which is in agreement with recent lattice results. Similarly for pions we obtain $m_{\pi^{\pm}}-m_{\pi^0}=1.8~{\rm MeV}$, roughly half of the experimental value. The electromagnetic mass of pions is found to be independent of $N_c$ and 't Hooft coupling and its scale is set only by the Kaluza-Klein scale of the model, $M_{\rm KK}=949~{\rm MeV}$.


title: Where do rivers grow? Path selection and growth in a harmonic field

abstract: River networks exhibit a complex ramified structure that has inspired decades of studies. Yet, an understanding of the propagation of a single stream remains elusive. Here we invoke a criterion for path selection from fracture mechanics and apply it to the growth of streams in a diffusion field. We show that a stream will follow local symmetry in order to maximize the water flux and that its trajectory is defined by the local field in its vicinity. We also study the growth of a real network. We use this principle to construct the history of a network and to find a growth law associated with it. The results show that the deterministic growth of a single channel based on its local environment can be used to characterize the structure of river networks.


title: An upper bound on the extremal version of Hajnal's triangle-free game

abstract: A game starts with the empty graph on $n$ vertices, and two player alternate adding edges to the graph. Only moves which do not create a triangle are valid. The game ends when a maximal triangle-free graph is reached. The goal of one player is to end the game as soon as possible, while the other player is trying to prolong the game. With optimal play, the length of the game (number of edges played) is called the $K_3$ game saturation number. In this paper we prove an upper bound for this number.


title: Fe I Oscillator Strengths for Transitions from High-lying Even-Parity Levels

abstract: New radiative lifetimes, measured to $\pm$ 5 % accuracy, are reported for 31 even-parity levels of Fe I ranging from 45061 cm$^{-1}$ to 56842 cm$^{-1}$. These lifetimes have been measured using single-step and two-step time-resolved laser-induced fluorescence on a slow atomic beam of iron atoms. Branching fractions have been attempted for all of these levels, and completed for 20 levels. This set of levels represents an extension of the collaborative work reported in Ruffoni et al. (2014). The radiative lifetimes combined with the branching fractions yields new oscillator strengths for 203 lines of Fe I. Utilizing a 1D-LTE model of the solar photosphere, spectral syntheses for a subset of these lines which are unblended in the solar spectrum yields a mean iron abundance of <log[{\epsilon}(Fe)]> = 7.45 $\pm$ 0.06.


title: Modulated traveling fronts for a nonlocal Fisher-KPP equation: a dynamical systems approach

abstract: We consider a nonlocal generalization of the Fisher-KPP equation in one spatial dimension. As a parameter is varied the system undergoes a Turing bifurcation. We study the dynamics near this Turing bifurcation. Our results are two-fold. First, we prove the existence of a two-parameter family of bifurcating stationary periodic solutions and derive a rigorous asymptotic approximation of these solutions. We also study the spectral stability of the bifurcating stationary periodic solutions with respect to almost co-periodic perturbations. Secondly, we restrict to a specific class of exponential kernels for which the nonlocal problem is transformed into a higher order partial differential equation. In this context, we prove the existence of modulated traveling fronts near the Turing bifurcation that describe the invasion of the Turing unstable homogeneous state by the periodic pattern established in the first part. Both results rely on a center manifold reduction to a finite dimensional ordinary differential equation.


title: Endeavours in Discrete Lorentzian Geometry; A thesis in five papers

abstract: To solve the path integral for quantum gravity, one needs to regularise the space-times that are summed over. This regularisation usually is a discretisation, which makes it necessary to give up some paradigms or symmetries of continuum physics. Causal dynamical triangulations regularises the path integral through a simplicial discretisation that introduces a preferred time foliation. The first part of this thesis presents three articles on causal dynamical triangulations. The first article shows how to obtain a multicritical 2d model by coupling the theory to hard dimers. The second explores the connection to Ho\v{r}ava-Lifshitz gravity that is suggested by the time foliation and establishes that in 2d the theories are equivalent. The last article does not directly concern causal dynamical triangulations but Euclidian dynamical triangulations with an additional measure term, which are examined to understand whether they contain an extended phase without the need for a preferred time foliation. Causal set theory uses an explicitly Lorentz invariant discretisation, which introduces non-local effects. The second part of this thesis presents two articles in causal set theory. The first explicitly calculates closed form expressions for the d'Alembertian operator in any dimension, which can be implemented in computer simulations. The second develops a ruler to examine the manifoldlikeness of small regions in a causal set, and can be used to recover locality.


title: Some properties of plurisubharmonic functions

abstract: Two properties of plurisubharmonic functions are proven. The first result is a Skoda type integrability theorem with respect to a Monge-Amp\`ere mass with H\"older continuous potential. The second one says that locally, a p.s.h. function is $k$-Lipschitz outside a set of Lebesgue measure smaller that $c/k^2$.


title: From High-Level Modeling Towards Efficient and Trustworthy Circuits

abstract: Behavior-Interaction-Priority (BIP) is a layered embedded system design and verification framework that provides separation of functionality, synchronization, and priority concerns to simplify system design and to establish correctness by construction. The framework comes with a runtime engine and a suite of verification tools that uses D-Finder and NuSMV as model checkers. In this paper we provide a method and a supporting tool that takes a BIP system and a set of invariants and computes a reduced sequential circuit with a system-specific scheduler and with a designated output that is true when the invariants hold. Our method uses ABC, a sequential circuit synthesis and verification framework to (1) generate an efficient FPGA implementation of the system, and to (2) verify the system and debug it in case a counterexample was found. Moreover we generate a concurrent C implementation of the circuit that can be directly used as a simulator. We evaluated our method with two large systems and our results outperform those possible with existing techniques.


title: Majorization-minimization procedures and convergence of SQP methods for semi-algebraic and tame programs

abstract: In view of solving nonsmooth and nonconvex problems involving complex constraints (like standard NLP problems), we study general maximization-minimization procedures produced by families of strongly convex sub-problems. Using techniques from semi-algebraic geometry and variational analysis -in particular Lojasiewicz inequality- we establish the convergence of sequences generated by this type of schemes to critical points. The broad applicability of this process is illustrated in the context of NLP. In that case critical points coincide with KKT points. When the data are semi-algebraic or real analytic our method applies (for instance) to the study of various SQP methods: the moving balls method, Sl1QP, ESQP. Under standard qualification conditions, this provides -to the best of our knowledge- the first general convergence results for general nonlinear programming problems. We emphasize the fact that, unlike most works on this subject, no second-order assumption and/or convexity assumptions whatsoever are made. Rate of convergence are shown to be of the same form as those commonly encountered with first order methods.


title: Mini-review on QCD spectral sum rules

abstract: Taking the example of the most popular and well-established Borel / Laplace / Exponential sum rule (LSR), I shortly review some of its recent applications in hadron physics namely the estimates of non-perturbative condensates, the determination of the light and heavy quark masses, the extraction of the heavy-light decay constants, the estimates of charmonium and bottomium molecule masses and the properties of scalar gluonium. In addition to the standard \tau (sum rule variable) and t_c (QCD continuum threshold) stablity criteria, I introduce a new stability criterion versus the arbitrary QCD subtraction point \mu for extracting the optimal results when radiative QCD corrections are included. Future improvements on further uses of QCD spectral sum rules (QSSR) are commented.


title: Prospects for joint gravitational wave and short gamma-ray burst observations

abstract: We present a detailed evaluation of the expected rate of joint gravitational-wave and short gamma-ray burst (GRB) observations over the coming years. We begin by evaluating the improvement in distance sensitivity of the gravitational wave search that arises from using the GRB observation to restrict the time and sky location of the source. We argue that this gives a 25% increase in sensitivity when compared to an all-sky, all-time search, corresponding to more than doubling the number of detectable gravitational wave signals associated with GRBs. Using this, we present the expected rate of joint observations with the advanced LIGO and Virgo instruments, taking into account the expected evolution of the gravitational wave detector network. We show that in the early advanced gravitational wave detector observing runs, from 2015-2017, there is only a small chance of a joint observation. However, as the detectors approach their design sensitivities, there is a good chance of joint observations provided wide field GRB satellites, such as Fermi and the Interplanetary Network, continue operation. The rate will also depend critically upon the nature of the progenitor, with neutron star--black hole systems observable to greater distances than double neutron star systems. The relative rate of binary mergers and GRBs will depend upon the jet opening angle of GRBs. Consequently, joint observations, as well as accurate measurement of both the GRB rate and binary merger rates, will allow for an improved estimation of the opening angle of GRBs.


title: Near-optimal estimation of jump activity in semimartingales

abstract: In quantitative finance, we often model asset prices as semimartingales, with drift, diffusion and jump components. The jump activity index measures the strength of the jumps at high frequencies, and is of interest both in model selection and fitting, and in volatility estimation. In this paper, we give a novel estimate of the jump activity, together with corresponding confidence intervals. Our estimate improves upon previous work, achieving near-optimal rates of convergence, and good finite-sample performance in Monte-Carlo experiments.


title: Running of scalar spectral index in multi-field inflation

abstract: We compute the running of the scalar spectral index in general multi-field slow-roll inflation. By incorporating explicit momentum dependence at the moment of horizon crossing, we can find the running straightforwardly. At the same time, we can distinguish the contributions from the quasi de Sitter background and the super-horizon evolution of the field fluctuations.


title: Controversy and Sentiment in Online News

abstract: How do news sources tackle controversial issues? In this work, we take a data-driven approach to understand how controversy interplays with emotional expression and biased language in the news. We begin by introducing a new dataset of controversial and non-controversial terms collected using crowdsourcing. Then, focusing on 15 major U.S. news outlets, we compare millions of articles discussing controversial and non-controversial issues over a span of 7 months. We find that in general, when it comes to controversial issues, the use of negative affect and biased language is prevalent, while the use of strong emotion is tempered. We also observe many differences across news sources. Using these findings, we show that we can indicate to what extent an issue is controversial, by comparing it with other issues in terms of how they are portrayed across different media.


title: Looking for good Hofmeister and Braunschadel bases

abstract: A_k = (1, a_2, ... a_k} is an h-basis for n if every positive integer not exceeding n can be expressed as the sum of no more than h values a_i. An "extremal" h-basis A_k is one for which n is as large as possible. Computing extremal bases has become known as the "global" Postage Stamp Problem. This paper describes the author's early attempts to identify extremal Hofmeister and Braunschadel bases for large t (where t = 12h + r for some 0 <= r <= 11), and also includes a transcription of Prof. Selmer's correspondence with the author about this work which has not been published before. See (Selmer, E.S., Private Communication, [8], 1992/1993).


title: A Schur-Weyl Duality Approach to Walking on Cubes

abstract: Walks on the representation graph $\mathcal R_{\mathsf{V}}(\mathsf{G})$ determined by a group $\mathsf{G}$ and a $\mathsf{G}$-module $\mathsf{V}$ are related to the centralizer algebras of the action of $\mathsf{G}$ on the tensor powers $\mathsf{V}^{\otimes k}$ via Schur-Weyl duality. This paper explores that connection when the group is $\mathbb{Z}_2^n$ and the module $\mathsf{V}$ is chosen so the representation graph is the $n$-cube. We describe a basis for the centralizer algebras in terms of labeled partition diagrams. We obtain an expression for the number of walks by counting certain partitions and determine the exponential generating functions for the number of walks


title: The p-adic Shintani modular symbol and evil Eisenstein series

abstract: We compute the p-adic L-function (in the sense of Pollack-Stevens and Bellaiche) of the "evil twin" Eisenstein series. We use Solomon's Shintani cocycle to construct a p-adic Shintani modular symbol, and show how to specialize this modular symbol to overconvergent modular symbols with the same p-adic L-function as evil/critical-slope Eisenstein series. This gives a new proof of recent results of Bellaiche and Dasgupta.


title: Kepler-424 b: A "Lonely" Hot Jupiter That Found A Companion

abstract: Hot Jupiter systems provide unique observational constraints for migration models in multiple systems and binaries. We report on the discovery of the Kepler-424 (KOI-214) two-planet system, which consists of a transiting hot Jupiter (Kepler-424b) in a 3.31-d orbit accompanied by a more massive outer companion in an eccentric (e=0.3) 223-d orbit. The outer giant planet, Kepler-424c, is not detected to transit the host star. The masses of both planets and the orbital parameters for the second planet were determined using precise radial velocity (RV) measurements from the Hobby-Eberly Telescope (HET) and its High Resolution Spectrograph (HRS). In stark contrast to smaller planets, hot Jupiters are predominantly found to be lacking any nearby additional planets, the appear to be "lonely" (e.g. Steffen et al.~2012). This might be a consequence of a highly dynamical past of these systems. The Kepler-424 planetary system is a system with a hot Jupiter in a multiple system, similar to upsilon Andromedae. We also present our results for Kepler-422 (KOI-22), Kepler-77 (KOI-127; Gandolfi et al.~2013), Kepler-43 (KOI-135; Bonomo et al.~2012), and Kepler-423 (KOI-183). These results are based on spectroscopic data collected with the Nordic Optical Telescope (NOT), the Keck 1 telescope and HET. For all systems we rule out false positives based on various follow-up observations, confirming the planetary nature of these companions. We performed a comparison with planetary evolutionary models which indicate that these five hot Jupiters have a heavy elements content between 20 and 120 M_Earth.


title: Modelling CO emission from hydrodynamic simulations of nearby spirals, starbursting mergers, and high-redshift galaxies

abstract: We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.


title: CO excitation of normal star forming galaxies out to z=1.5 as regulated by the properties of their interstellar medium

abstract: We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.


title: Pursuit on a Graph Using Partial Information

abstract: The optimal control of a "blind" pursuer searching for an evader moving on a road network and heading at a known speed toward a set of goal vertices is considered. To aid the "blind" pursuer, certain roads in the network have been instrumented with Unattended Ground Sensors (UGSs) that detect the evader's passage. When the pursuer arrives at an instrumented node, the UGS therein informs the pursuer if and when the evader visited the node. The pursuer's motion is not restricted to the road network. In addition, the pursuer can choose to wait/loiter for an arbitrary time at any UGS location/node. At time 0, the evader passes by an entry node on his way towards one of the exit nodes. The pursuer also arrives at this entry node after some delay and is thus informed about the presence of the intruder/evader in the network, whereupon the chase is on - the pursuer is tasked with capturing the evader. Because the pursuer is "blind", capture entails the pursuer and evader being collocated at an UGS location. If this happens, the UGS is triggered and this information is instantaneously relayed to the pursuer, thereby enabling capture. On the other hand, if the evader reaches one of the exit nodes without being captured, he is deemed to have escaped. We provide an algorithm that computes the maximum initial delay at the entry node for which capture is guaranteed. The algorithm also returns the corresponding optimal pursuit policy.


title: Characterization of a multimode coplanar waveguide parametric amplifier

abstract: We characterize a novel Josephson parametric amplifier based on a flux-tunable quarter-wavelength resonator. The fundamental resonance frequency is ~1GHz, but we use higher modes of the resonator for our measurements. An on-chip tuning line allows for magnetic flux pumping of the amplifier. We investigate and compare degenerate parametric amplification, involving a single mode, and nondegenerate parametric amplification, using a pair of modes. We show that we reach quantum-limited noise performance in both cases, and we show that the added noise can be less than 0.5 added photons in the case of low gain.


title: Long-term BVRI light curves of 5 pre-main sequence stars in the field of "Gulf of Mexico"

abstract: We present new data from BVRI photometric observations of five PMS stars during the period from April 2013 to July 2014. The stars are located in the field of NGC 7000/IC 5070 ("Gulf of Mexico") - a region with active star formation. The presented paper is a continuation of our long-term photometric investigations of the young stellar objects in this region. The long-term multicolor photometric observations of PMS stars are very important for their exact classification. Our results show that the studied stars exhibit different types of photometric variability in all bands. We tried to classify them using our data from the long-term photometry and data published by other authors.


title: Ultra-High-Energy Cosmic Rays from Low-Luminosity Active Galactic Nuclei

abstract: We investigate the production of ultra-high-energy cosmic ray (UHECR) in relativistic jets from low-luminosity active galactic nuclei (LLAGN). We start by proposing a model for the UHECR contribution from the black holes (BHs) in LLAGN, which present a jet power $P_{\mathrm{j}} \leqslant 10^{46}$ erg s$^{-1}$. This is in contrast to the opinion that only high-luminosity AGN can accelerate particles to energies $ \geqslant 50$ EeV. We rewrite the equations which describe the synchrotron self-absorbed emission of a non-thermal particle distribution to obtain the observed radio flux density from sources with a flat-spectrum core and its relationship to the jet power. We find that the UHECR flux is dependent on the {\it observed radio flux density, the distance to the AGN, and the BH mass}, where the particle acceleration regions can be sustained by the magnetic energy extraction from the BH at the center of the AGN. We use a complete sample of 29 radio sources with a total flux density at 5 GHz greater than 0.5 Jy to make predictions for the maximum particle energy, luminosity, and flux of the UHECRs from nearby AGN. These predictions are then used in a semi-analytical code developed in Mathematica (SAM code) as inputs for the Monte-Carlo simulations to obtain the distribution of the arrival direction at the Earth and the energy spectrum of the UHECRs, taking into account their deflection in the intergalactic magnetic fields. For comparison, we also use the CRPropa code with the same initial conditions as for the SAM code. Importantly, to calculate the energy spectrum we also include the weighting of the UHECR flux per each UHECR source. Next, we compare the energy spectrum of the UHECRs with that obtained by the Pierre Auger Observatory.


title: Development of the Method of Averaging in Clifford Geometric Algebras

abstract: We develop the method of averaging in Clifford (geometric) algebras suggested by the author in previous papers. We consider operators constructed using two different sets of anticommuting elements of real or complexified Clifford algebras. These operators generalize Reynolds operators from the representation theory of finite groups. We prove a number of new properties of these operators. Using the generalized Reynolds operators, we give a complete proof of the generalization of Pauli's theorem to the case of Clifford algebras of arbitrary dimension. The results can be used in geometry, physics, engineering, computer science, and other applications.


title: The iEBE-VISHNU code package for relativistic heavy-ion collisions

abstract: The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion.


title: Theory for Baryon Number and Dark Matter at the LHC

abstract: We investigate the possibility to test the simplest theory for spontaneous baryon number violation at the Large Hadron Collider. In this context the baryon number is a local gauge symmetry spontaneously broken at the low scale through the Brout-Englert-Higgs mechanism. This theory predicts the existence of a leptophobic neutral gauge boson and a fermionic dark matter candidate with baryon number. We study the gauge boson and Higgs decays, and explore the connection between collider signatures and constraints coming from dark matter experiments. We point out an upper bound on the symmetry breaking scale using the relic density constraints which tells us that this model can be tested or ruled out at current or future collider experiments.


title: Stream Processing for Solar Physics: Applications and Implications for Big Solar Data

abstract: Modern advances in space technology have enabled the capture and recording of unprecedented volumes of data. In the field of solar physics this is most readily apparent with the advent of the Solar Dynamics Observatory (SDO), which returns in excess of 1 terabyte of data daily. While we now have sufficient capability to capture, transmit and store this information, the solar physics community now faces the new challenge of analysis and mining of high-volume and potentially boundless data sets such as this: a task known to the computer science community as stream mining. In this paper, we survey existing and established stream mining methods in the context of solar physics, with a goal of providing an introductory overview of stream mining algorithms employed by the computer science fields. We consider key concepts surrounding stream mining that are applicable to solar physics, outlining existing algorithms developed to address this problem in other fields of study, and discuss their applicability to massive solar data sets. We also discuss the considerations and trade-offs that may need to be made when applying stream mining methods to solar data. We find that while no one single solution is readily available, many of the methods now employed in other data streaming applications could successfully be modified to apply to solar data and prove invaluable for successful analysis and mining of this new source.


title: On Hoelder-continuity of Oseledets subspaces

abstract: For Hoelder cocycles over a Lipschitz base transformation, possibly non-invertible, we show that the subbundles given by the Oseledets Theorem are Hoelder-continuous on compact sets of measure arbitrarily close to 1. The results extend to vector bundle automorphisms, as well as to the Kontsevich-Zorich cocycle over the Teichmueller flow on the moduli space of abelian differentials. Following a recent result of Chaika-Eskin, our results also extend to any given Teichmueller disk.


title: Quaternion Derivatives: The GHR Calculus

abstract: Quaternion derivatives in the mathematical literature are typically defined only for analytic (regular) functions. However, in engineering problems, functions of interest are often real-valued and thus not analytic, such as the standard cost function. The HR calculus is a convenient way to calculate formal derivatives of both analytic and non-analytic functions of quaternion variables, however, both the HR and other functional calculus in quaternion analysis have encountered an essential technical obstacle, that is, the traditional product rule is invalid due to the non- commutativity of the quaternion algebra. To address this issue, a generalized form of the HR derivative is proposed based on a general orthogonal system. The so introduced generalization, called the generalized HR (GHR) calculus, encompasses not just the left- and right-hand versions of quaternion derivative, but also enables solutions to some long standing problems, such as the novel product rule, the chain rule, the mean-valued theorem and Taylor's theorem. At the core of the proposed approach is the quaternion rotation, which can naturally be applied to other functional calculi in non-commutative settings. Examples on using the GHR calculus in adaptive signal processing support the analysis.


title: Holderian weak invariance principle for stationary mixing sequences

abstract: We provide some sufficient mixing conditions on a strictly stationary sequence in order to guarantee the weak invariance principle in H\"older spaces. Strong mixing and $\rho$-mixing conditions are investigated as well as $\tau$-dependent sequences. The main tools are Fuk-Nagaev type inequalities for mixing sequences and a truncation argument.


title: Effective dynamics of strongly dissipative Rydberg gases

abstract: We investigate the evolution of interacting Rydberg gases in the limit of strong noise and dissipation. Starting from a description in terms of a Markovian quantum master equation we derive effective equations of motion that govern the dynamics on a "coarse-grained" timescale where fast dissipative degrees of freedom have been adiabatically eliminated. Specifically, we consider two scenarios which are of relevance for current theoretical and experimental studies --- Rydberg atoms in a two-level (spin) approximation subject to strong dephasing noise as well as Rydberg atoms under so-called electromagnetically induced transparency (EIT) conditions and fast radiative decay. In the former case we find that the effective dynamics is described by classical rate equations up to second order in an appropriate perturbative expansion. This drastically reduces the computational complexity of numerical simulations in comparison to the full quantum master equation. When accounting for the fourth order correction in this expansion, however, we find that the resulting equation breaks the preservation of positivity and thus cannot be interpreted as a proper classical master rate equation. In the EIT system we find that the expansion up to second order retains information not only on the "classical" observables, but also on some quantum coherences. Nevertheless, this perturbative treatment still achieves a non-trivial reduction of complexity with respect to the original problem.


title: An Analysis of BitTorrent Cross-Swarm Peer Participation and Geolocational Distribution

abstract: Peer-to-Peer (P2P) file-sharing is becoming increasingly popular in recent years. In 2012, it was reported that P2P traffic consumed over 5,374 petabytes per month, which accounted for approximately 20.5% of consumer internet traffic. TV is the popular content type on The Pirate Bay (the world's largest BitTorrent indexing website). In this paper, an analysis of the swarms of the most popular pirated TV shows is conducted. The purpose of this data gathering exercise is to enumerate the peer distribution at different geolocational levels, to measure the temporal trend of the swarm and to discover the amount of cross-swarm peer participation. Snapshots containing peer related information involved in the unauthorised distribution of this content were collected at a high frequency resulting in a more accurate landscape of the total involvement. The volume of data collected throughout the monitoring of the network exceeded 2 terabytes. The presented analysis and the results presented can aid in network usage prediction, bandwidth provisioning and future network design.


title: Some Banach spaces added by a Cohen real

abstract: We study certain Banach spaces that are added in the extension by one Cohen real. Specifically, we show that adding just one Cohen real to any model adds a Banach space of density $\aleph_1$ which does not embed into any such space in the ground model such a Banach space can be chosen to be UG This has consequences on the the isomorphic universality number for Banach spaces of density $\aleph_1$, which is hence equal to $\aleph_2$ in the standard Cohen model and the same is true for UG spaces. Analogous universality results for Banach spaces are true for other cardinals, by a different proof.


title: Triple product formula and mass equidistribution on modular curves of level N

abstract: It was shown in previous works that the measure associated to holomorphic newforms of weight $k$ and level $q$ will tend weakly to the Haar measure on modular curve of level 1, as $qk\rightarrow \infty$. In this paper we proved that this phenomenon is also true on modular curves of general level N.


title: BitTorrent Sync: First Impressions and Digital Forensic Implications

abstract: With professional and home Internet users becoming increasingly concerned with data protection and privacy, the privacy afforded by popular cloud file synchronisation services, such as Dropbox, OneDrive and Google Drive, is coming under scrutiny in the press. A number of these services have recently been reported as sharing information with governmental security agencies without warrants. BitTorrent Sync is seen as an alternative by many and has gathered over two million users by December 2013 (doubling since the previous month). The service is completely decentralised, offers much of the same synchronisation functionality of cloud powered services and utilises encryption for data transmission (and optionally for remote storage). The importance of understanding BitTorrent Sync and its resulting digital investigative implications for law enforcement and forensic investigators will be paramount to future investigations. This paper outlines the client application, its detected network traffic and identifies artefacts that may be of value as evidence for future digital investigations.


title: Kinetic Model of Mass Exchange with Dynamic Arrhenius Transition Rates

abstract: We study a nonlinear kinetic model of mass exchange between interacting grains. The transition rates follow the Arrhenius equation with an activation energy that depends on the grain mass. We show that the activation parameter can be absorbed in the initial conditions for the grain masses, and that the total mass is conserved. We obtain numerical solutions of the coupled, nonlinear, ordinary differential equations of mass exchange for the two-grain system, and we compare them with approximate theoretical solutions in specific neighborhoods of the phase space. Using phase plane methods, we determine that the system exhibits regimes of diffusive and growth-decay (reverse diffusion) kinetics. The equilibrium states are determined by the mass equipartition and separation nullcline curves. If the transfer rates are perturbed by white noise, numerical simulations show that the system still exhibits diffusive and growth-decay regimes, although the noise can reverse the sign of equilibrium mass difference. Finally, we present theoretical analysis and numerical simulations of a system with many interacting grains. Diffusive and growth-decay regimes are established as well, but the approach to equilibrium is considerably slower. Potential applications of the mass exchange model involve coarse-graining during sintering and wealth exchange in econophysics.


title: Aloof Baby Skyrmions

abstract: We show that a suitable choice for the potential term in the two-dimensional baby Skyrme model yields solitons that have a short-range repulsion and a long-range attraction. The solitons are therefore aloof, in the sense that static multi-soliton bound states have constituents that preserve their individual identities and are sufficiently far apart that tail interactions yield small binding energies. The static multi-soliton solutions are found to have a cluster structure that is reproduced by a simple binary species particle model. In the standard three-dimensional Skyrme model of nuclei, solitons are too tightly bound and are often too symmetric, due to symmetry enhancement as solitons coalesce to form bound states. The aloof baby Skyrmion results endorse a way to resolve these issues and provides motivation for a detailed study of the related three-dimensional version of the Skyrme model.


title: A Combinatorial Formula for Certain Elements of Upper Cluster Algebras

abstract: We develop an elementary formula for certain non-trivial elements of upper cluster algebras. These elements have positive coefficients. We show that when the cluster algebra is acyclic these elements form a basis. Using this formula, we show that each non-acyclic skew-symmetric cluster algebra of rank 3 is properly contained in its upper cluster algebra.


title: Accounting for Classical Hardware in the Control of Quantum Devices

abstract: High fidelity coherent control of quantum systems is critical to building quantum devices and quantum computers. We provide a general optimal control framework for designing control sequences that account for hardware control distortions while maintaining robustness to environmental noise. We demonstrate the utility of our algorithm by presenting examples of robust quantum gates optimized in the presence of nonlinear distortions. We show that nonlinear classical controllers do not necessarily incur additional computational cost to pulse optimization, enabling more powerful quantum devices.


title: Nonconvex Sorted $\ell_1$ Minimization for Sparse Approximation

abstract: The $\ell_1$ norm is the tight convex relaxation for the $\ell_0$ "norm" and has been successfully applied for recovering sparse signals. For problems with fewer samplings, one needs to enhance the sparsity by nonconvex penalties such as $\ell_p$ "norm". As one method for solving $\ell_p$ minimization problems, iteratively reweighted $\ell_1$ minimization updates the weight for each component based on the value of the same component at the previous iteration. It assigns large weights on small components in magnitude and small weights on large components in magnitude. In this paper, we consider a weighted $\ell_1$ penalty with the set of the weights fixed and the weights are assigned based on the sort of all the components in magnitude. The smallest weight is assigned to the largest component in magnitude. This new penalty is called nonconvex sorted $\ell_1$. Then we propose two methods for solving nonconvex sorted $\ell_1$ minimization problems: iteratively reweighted $\ell_1$ minimization and iterative sorted thresholding, and prove that both methods will converge to a local optimum. We also show that both methods are generalizations of iterative support detection and iterative hard thresholding respectively. The numerical experiments demonstrate the better performance of assigning weights by sort compared to $\ell_p$ minimization.


title: Localized shocks

abstract: We study products of precursors of spatially local operators, $W_{x_{n}}(t_{n}) ... W_{x_1}(t_1)$, where $W_x(t) = e^{-iHt} W_x e^{iHt}$. Using chaotic spin-chain numerics and gauge/gravity duality, we show that a single precursor fills a spatial region that grows linearly in $t$. In a lattice system, products of such operators can be represented using tensor networks. In gauge/gravity duality, they are related to Einstein-Rosen bridges supported by localized shock waves. We find a geometrical correspondence between these two descriptions, generalizing earlier work in the spatially homogeneous case.


title: Modified energy-momentum conservation laws and vacuum Cherenkov radiation

abstract: We present a general parametrization for the leading order terms in a momentum power expansion of a non-universal Lorentz-violating, but rotational invariant, kinematics and its implications for two-body decay thresholds. The considered framework includes not only modified dispersion relations for particles, but also modified energy-momentum conservation laws, something which goes beyond effective field theory. As a particular and relevant example, bounds on the departures from special relativistic kinematics from the non-observation of vacuum Cherenkov radiation are discussed and compared with those obtained within the effective field theory scenario.


title: Graphene as Gain Medium for Broadband Lasers

abstract: In contrast to conventional structures, efficient non-radiative carrier recombination counteracts the appearance of optical gain in graphene. Based on a microscopic and fully quantum-mechanical study of the coupled carrier, phonon, and photon dynamics in graphene, we present a strategy to obtain a long-lived gain: Integrating graphene into a photonic crystal nanocavity and applying a high-dielectric substrate gives rise to pronounced coherent light emission suggesting the design of graphene-based laser devices covering a broad spectral range.


title: An Improved Constraint-Tightening Approach for Stochastic MPC

abstract: The problem of achieving a good trade-off in Stochastic Model Predictive Control between the competing goals of improving the average performance and reducing conservativeness, while still guaranteeing recursive feasibility and low computational complexity, is addressed. We propose a novel, less restrictive scheme which is based on considering stability and recursive feasibility separately. Through an explicit first step constraint we guarantee recursive feasibility. In particular we guarantee the existence of a feasible input trajectory at each time instant, but we only require that the input sequence computed at time $k$ remains feasible at time $k+1$ for most disturbances but not necessarily for all, which suffices for stability. To overcome the computational complexity of probabilistic constraints, we propose an offline constraint-tightening procedure, which can be efficiently solved via a sampling approach to the desired accuracy. The online computational complexity of the resulting Model Predictive Control (MPC) algorithm is similar to that of a nominal MPC with terminal region. A numerical example, which provides a comparison with classical, recursively feasible Stochastic MPC and Robust MPC, shows the efficacy of the proposed approach.


title: Totally disconnected locally compact groups with a linear open subgroup

abstract: We describe the global structure of totally disconnected locally compact groups having a linear open compact subgroup. Among the applications, we show that if a non-discrete, compactly generated, topologically simple, totally disconnected locally compact group is locally linear, then it is a simple algebraic group over a local field.


title: Adaptive Low-Complexity Sequential Inference for Dirichlet Process Mixture Models

abstract: We develop a sequential low-complexity inference procedure for Dirichlet process mixtures of Gaussians for online clustering and parameter estimation when the number of clusters are unknown a-priori. We present an easily computable, closed form parametric expression for the conditional likelihood, in which hyperparameters are recursively updated as a function of the streaming data assuming conjugate priors. Motivated by large-sample asymptotics, we propose a novel adaptive low-complexity design for the Dirichlet process concentration parameter and show that the number of classes grow at most at a logarithmic rate. We further prove that in the large-sample limit, the conditional likelihood and data predictive distribution become asymptotically Gaussian. We demonstrate through experiments on synthetic and real data sets that our approach is superior to other online state-of-the-art methods.