title: Gr\"obner Bases over Algebraic Number Fields

abstract: Although Buchberger's algorithm, in theory, allows us to compute Gr\"obner bases over any field, in practice, however, the computational efficiency depends on the arithmetic of the ground field. Consider a field $K = \mathbb{Q}(\alpha)$, a simple extension of $\mathbb{Q}$, where $\alpha$ is an algebraic number, and let $f \in \mathbb{Q}[t]$ be the minimal polynomial of $\alpha$. In this paper we present a new efficient method to compute Gr\"obner bases in polynomial rings over the algebraic number field $K$. Starting from the ideas of Noro [Noro, 2006], we proceed by joining $f$ to the ideal to be considered, adding $t$ as an extra variable. But instead of avoiding superfluous S-pair reductions by inverting algebraic numbers, we achieve the same goal by applying modular methods as in [Arnold, 2003; B\"ohm et al., 2015; Idrees et al., 2011], that is, by inferring information in characteristic zero from information in characteristic $p > 0$. For suitable primes $p$, the minimal polynomial $f$ is reducible over $\mathbb{F}_p$. This allows us to apply modular methods once again, on a second level, with respect to the factors of $f$. The algorithm thus resembles a divide and conquer strategy and is in particular easily parallelizable. At current state, the algorithm is probabilistic in the sense that, as for other modular Gr\"obner basis computations, an effective final verification test is only known for homogeneous ideals or for local monomial orderings. The presented timings show that for most examples, our algorithm, which has been implemented in SINGULAR, outperforms other known methods by far.


title: Real-time correction of panoramic images using hyperbolic M\"obius transformations

abstract: Wide-angle images gained a huge popularity in the last years due to the development of computational photography and imaging technological advances. They present the information of a scene in a way which is more natural for the human eye but, on the other hand, they introduce artifacts such as bent lines. These artifacts become more and more unnatural as the field of view increases. In this work, we present a technique aimed to improve the perceptual quality of panorama visualization. The main ingredients of our approach are, on one hand, considering the viewing sphere as a Riemann sphere, what makes natural the application of M\"obius (complex) transformations to the input image, and, on the other hand, a projection scheme which changes in function of the field of view used. We also introduce an implementation of our method, compare it against images produced with other methods and show that the transformations can be done in real-time, which makes our technique very appealing for new settings, as well as for existing interactive panorama applications.


title: Connecting the latent multinomial

abstract: Link et al. (2010) define a general framework for analyzing capture-recapture data with potential misidentifications. In this framework, the observed vector of counts, $y$, is considered as a linear function of a vector of latent counts, $x$, such that $y = A x$, with $x$ assumed to follow a multinomial distribution conditional on the model parameters, $\theta$. Bayesian methods are then applied by sampling from the joint posterior distribution of both $x$ and $\theta$. In particular, Link et al. (2010) propose a Metropolis-Hastings algorithm to sample from the full conditional distribution of $x$, where new proposals are generated by sequentially adding elements from a basis of the null space (kernel) of $A$. We consider this algorithm and show that using elements from a simple basis for the kernel of $A$ may not produce an irreducible Markov chain. Instead, we require a Markov basis, as defined by Diaconis and Sturmfels (1998). We illustrate the importance of Markov bases with three capture-recapture examples. We prove that a specific lattice basis is a Markov basis for a class of models including the original model considered by Link et al. (2010) and confirm that the specific basis used by Link et al. (2010) for their example with two sampling occasions is a Markov basis. The constructive nature of our proof provides an immediate method to obtain a Markov basis for any model in this class.


title: Energy level properties of 4p$^6$4d$^3$, 4p$^6$4d$^2$4f and 4p$^5$4d$^4$ configurations of W$^{35+}$ ion

abstract: The ab initio quasirelativistic Hartree-Fock method developed specifically for the calculation of spectroscopic parameters of heavy atoms and highly charged ions was used to derive spectral data for the multicharged tungsten ion W$^{35+}$. The configuration interaction method was applied to include the electron-correlation effects. The relativistic effects were taken into account in the Breit-Pauli approximation for quasirelativistic Hartree-Fock radial orbitals. The energy level spectra, radiative lifetimes $\tau$, Lande $g$-factors are calculated for the $\mathrm{4p^64d^3}$, $\mathrm{4p^64d^24f}$ and $\mathrm{4p^54d^4}$ configurations of the W$^{35+}$ ion.


title: Thickness-shear Vibration Frequencies of an Infinite Plate with a Generalized Material Property Grading along the Thickness

abstract: For quartz crystal resonators of thickness-shear type, the vibration frequency and mode shapes, which are key features of resonators in circuit applications, reflect the basic material and structural properties of the quartz plate and its variation with time under various factors such as erosive gases and liquids that can cause surface and internal damages and degradation of crystal blanks. The accumulated effects eventually will change the surface conditions in terms of elastic constants and stiffness and more importantly, the gradient of such properties along the thickness. This is a typical functionally graded materials (FGM) structure and has been studied extensively for structural applications under multiple loadings such as thermal and electromagnetic fields in recent years. For acoustic wave resonators, such studies are equally important and the wave propagation in FGM structures can be used in the evaluation and assessment of performance, reliability, and life of sensors based on acoustic waves such as the quartz crystal microbalances (QCM). Now we studied the thickness-shear vibrations of FGM plates with properties of AT-cut quartz crystal varying along the thickness in a general pattern represented by a trigonometric function with both sine and cosine functions of the thickness coordinate. The solutions are obtained by using Fourier expansion of the plate deformation. We also obtained the frequency changes of the fundamental and overtone modes which are strongly coupled for the evaluation of resonator structures with property variation or design to take advantages of FGM in novel applications.


title: A formula for the numerical range of an elementary operator on C*-algebra

abstract: Let $\A$ be a $C^*-$algebra with unit element $1$ and unitary group $U.$ Let $a=(a_1, ..., a_k)$ and $b=(b_1, ..., b_k)$ two $k-$tuples of elements in $\A.$ The elementary operator associated to $a$ and $b$ is defined by $ R_{a, b} (x)= \sum_{i=1}^ {k} a_ixb_i.$ In this paper we prove the following formula for the numerical range of $R_{a,b}$: $$V(R_{a, b}, B(\A))= [\cup \{ V(\sum_{i=1}^ {k}u^*a_iub_i, \A): u\in U\} ]^-.$$ This formulat solves the problem 4.5 of \cite{Fi}.


title: Nematic resonance in the Raman response of iron-based superconductors

abstract: In a fully-gapped superconductor the electronic Raman response has a pair-breaking peak at twice the superconducting gap $\Delta$, if the Bogoliubov excitations are uncorrelated. Motivated by the iron based superconductors, we study how this peak is modified if the superconducting phase hosts a nematic-structural quantum critical point. We show that, upon approaching this point by tuning, e.g., doping, the growth of nematic correlations between the quasiparticles transforms the pair-breaking peak into a nematic resonance. The mode energy is below 2$\Delta$, and stays finite at the quantum critical point, where its spectral weight is sharply enhanced. The latter is consistent with recent experiments on electron-doped iron based superconductors and provides direct evidence of nematic correlations in their superconducting phases.


title: Ensemble Forecasting of Major Solar Flares -- First Results

abstract: We present the results from the first ensemble prediction model for major solar flares (M and X classes). The primary aim of this investigation is to explore the construction of an ensemble for an initial prototyping of this new concept. Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemble forecast by linearly combining the flaring probabilities from all four methods. Performance-based combination weights were calculated using a Monte-Carlo-type algorithm that applies a decision threshold $P_{th}$ to the combined probabilities and maximizing the Heidke Skill Score (HSS). Using the data for 13 recent solar active regions between years 2012 - 2014, we found that linear combination methods can improve the overall probabilistic prediction and improve the categorical prediction for certain values of decision thresholds. Combination weights vary with the applied threshold and none of the tested individual forecasting models seem to provide more accurate predictions than the others for all values of $P_{th}$. According to the maximum values of HSS, a performance-based weights calculated by averaging over the sample, performed similarly to a equally weighted model. The values $P_{th}$ for which the ensemble forecast performs the best are 25 \% for M-class flares and 15 \% for X-class flares. When the human-adjusted probabilities from NOAA are excluded from the ensemble, the ensemble performance in terms of the Heidke score, is reduced.


title: Integrable light-cone lattice discretizations from the universal R-matrix

abstract: Our goal is to develop a more general scheme for constructing integrable lattice regularisations of integrable quantum field theories. Considering the affine Toda theories as examples, we show how to construct such lattice regularisations using the representation theory of quantum affine algebras. This requires us to clarify in particular the relations between the light-cone approach to integrable lattice models and the representation theory of quantum affine algebras. Both are found to be related in a very natural way, suggesting a general scheme for the construction of generalised Baxter Q-operators. One of the main difficulties we need to deal with is coming from the infinite-dimensionality of the relevant families of representations. It is handled by means of suitable renormalisation prescriptions defining what may be called the modular double of quantum affine algebras. This framework allows us to give a representation-theoretic proof of finite-difference equations generalising the Baxter equation.


title: Representations of the Kauffamn skein algebra of small surfaces

abstract: We prove a uniqueness result for finite-dimensional representations of the Kauffman skein algebra $\mathcal{S}_A(S)$ of a surface $S$, when $A$ is a root of unity and when the surface $S$ is a sphere with at most four punctures or a torus with at most one puncture. We show that, if two irreducible representations of $\mathcal{S}_A(S)$ have the same classical shadow and the same puncture invariants, and if this classical shadow is sufficiently generic in the character variety $\mathcal{X}_{SL_2(\mathbb{C})} (S)$, then the two representations are isomorphic.


title: Tidal Love numbers of membrane worlds: Europa, Titan, and Co

abstract: Under tidal forcing, icy satellites with subsurface oceans deform as if the surface were a membrane stretched around a fluid layer. `Membrane worlds' is thus a fitting name for these bodies and membrane theory provides the perfect toolbox to predict tidal effects. I describe here a new membrane approach to tidal perturbations based on the general theory of viscoelastic-gravitational deformations of spherically symmetric bodies. The massive membrane approach leads to explicit formulas for viscoelastic tidal Love numbers which are exact in the limit of zero crust thickness. Formulas for load Love numbers come as a bonus. The accuracy on $k_2$ and $h_2$ is better than one percent if the crust thickness is less than five percents of the surface radius, which is probably the case for Europa and Titan. The new approach allows for density differences between crust and ocean and correctly includes crust compressibility. This last feature makes it more accurate than the incompressible propagator matrix method. Membrane formulas factorize shallow and deep interior contributions, the latter affecting Love numbers mainly through density stratification. I show that a screening effect explains why ocean stratification typically increases Love numbers instead of reducing them. For Titan, a thin and dense liquid layer at the bottom of a light ocean can raise $k_2$ by more than ten percents. The membrane approach can also deal with dynamical tides in a non-rotating body. I show that a dynamical resonance significantly decreases the tilt factor and may thus lead to underestimating Europa's crust thickness. Finally, the dynamical resonance increases tidal deformations and tidal heating in the crust if the ocean thickness is of the order of a few hundred meters.


title: Witnessing entanglement by proxy

abstract: Entanglement is a ubiquitous feature of low temperature systems and believed to be highly relevant for the dynamics of condensed matter properties and quantum computation even at higher temperatures. The experimental certification of this paradigmatic quantum effect in macroscopic high temperature systems is constrained by the limited access to the quantum state of the system. In this paper we show how macroscopic observables beyond the energy of the system can be exploited as proxy witnesses for entanglement detection. Using linear and semi-definite relaxations we show that all previous approaches to this problem can be outperformed by our proxies, i.e. entanglement can be certified at higher temperatures without access to any local observable. For an efficient computation of proxy witnesses one can resort to a generalized grand canonical ensemble, enabling entanglement certification even in complex systems with macroscopic particle numbers.


title: Decomposability of nonnegative r-potent operators on L2(X)

abstract: We investigate the decomposability of nonnegative compact r-potent operators on a separable Hilbert space L2(X). We provide a constructive algorithm to prove that basis functions of range spaces of nonnegative r-potent operators can be chosen to be all nonnegative and mutually orthogonal. We use this orthogonality to establish that nonnegative compact r-potent operators with range spaces of dimension strictly greater than r-1 are decomposable.


title: Topological Nodal Line Semimetal and Dirac Semimetal State in Antiperovskite Cu$_3$PdN

abstract: Based on first-principles calculation and effective model analysis, we propose that the cubic antiperovskite material Cu$_3$PdN can host a three-dimensional (3D) topological nodal line semimetal state when spin-orbit coupling (SOC) is ignored, which is protected by coexistence of time-reversal and inversion symmetry. There are three nodal line circles in total due to the cubic symmetry. "Drumhead"-like surface flat bands are also derived. When SOC is included, each nodal line evolves into a pair of stable 3D Dirac points as protected by C$_4$ crystal symmetry. This is remarkably distinguished from the Dirac semimetals known so far, such as Na$_3$Bi and Cd$_3$As$_2$, both having only one pair of Dirac points. Once C$_4$ symmetry is broken, the Dirac points are gapped and the system becomes a strong topological insulator with (1;111) Z$_2$ indices.


title: Ultracold Dipolar Molecules Composed of Strongly Magnetic Atoms

abstract: In a combined experimental and theoretical effort, we demonstrate a novel type of dipolar system made of ultracold bosonic dipolar molecules with large magnetic dipole moments. Our dipolar molecules are formed in weakly bound Feshbach molecular states from a sample of strongly magnetic bosonic erbium atoms. We show that the ultracold magnetic molecules can carry very large dipole moments and we demonstrate how to create and characterize them, and how to change their orientation. Finally, we confirm that the relaxation rates of molecules in a quasi-two dimensional geometry can be reduced by using the anisotropy of the dipole-dipole interaction and that this reduction follows a universal dipolar behavior.


title: An Intermediate Language and Estimator for Automated Design Space Exploration on FPGAs

abstract: We present the TyTra-IR, a new intermediate language intended as a compilation target for high-level language compilers and a front-end for HDL code generators. We develop the requirements of this new language based on the design-space of FPGAs that it should be able to express and the estimation-space in which each configuration from the design-space should be mappable in an automated design flow. We use a simple kernel to illustrate multiple configurations using the semantics of TyTra-IR. The key novelty of this work is the cost model for resource-costs and throughput for different configurations of interest for a particular kernel. Through the realistic example of a Successive Over-Relaxation kernel implemented both in TyTra-IR and HDL, we demonstrate both the expressiveness of the IR and the accuracy of our cost model.


title: Robust estimation of U-statistics

abstract: An important part of the legacy of Evarist Gin\'e is his fundamental contributions to our understanding of $U$-statistics and $U$-processes. In this paper we discuss the estimation of the mean of multivariate functions in case of possibly heavy-tailed distributions. In such situations, reliable estimates of the mean cannot be obtained by usual $U$-statistics. We introduce a new estimator, based on the so-called median-of-means technique. We develop performance bounds for this new estimator that generalizes an estimate of Arcones and Gin\'e (1993), showing that the new estimator performs, under minimal moment conditions, as well as classical $U$-statistics for bounded random variables. We discuss an application of this estimator to clustering.


title: Dirac Processes and Default Risk

abstract: We introduce Dirac processes, using Dirac delta functions, for short-rate-type pricing of financial derivatives. Dirac processes add spikes to the existing building blocks of diffusions and jumps. Dirac processes are Generalized Processes, which have not been used directly before because the dollar value of non-Real numbers is meaningless. However, short-rate pricing is based on integrals so Dirac processes are natural. This integration directly implies that jumps are redundant whilst Dirac processes expand expressivity of short-rate approaches. Practically, we demonstrate that Dirac processes enable high implied volatility for CDS swaptions that has been otherwise problematic in hazard rate setups.


title: Dynamical evolution of quintessence cosmology in a physical phase space

abstract: The phase space analysis of cosmological parameters $\Omega_{\phi}$ and $\gamma_{\phi}$ is given. Based on this, the well-known quintessence cosmology is studied with an exponential potential $V(\phi)=V_{0}\exp(-\lambda\phi)$. Given observational data, the current state of universe could be pinpointed in the phase diagrams, thus making the diagrams more informative. The scaling solution of quintessence usually is not supposed to give the cosmic accelerating expansion, but we prove it could educe the transient acceleration. We also find that the differential equations of system used widely in study of scalar field are incomplete, and then a numerical method is used to figure out the range of application.


title: Planck Intermediate Results. XXXVI. Optical identification and redshifts of Planck SZ sources with telescopes in the Canary Islands Observatories

abstract: We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).


title: Diophantine Equations and Congruent Number Equation Solutions

abstract: By using pairs of nontrivial rational solutions of congruent number equation $$ C_N:\;\;y^2=x^3-N^2x, $$ constructed are pairs of rational right (Pythagorean) triangles with one common side and the other sides equal to the sum and difference of the squares of the same rational numbers. The parametrizations are found for following Diophantine systems: \begin{align*} (p^2\pm q^2)^2-a^2 & =\square_{1,2}\,, \\[0.2cm] c^2-(p^2\pm q^2)^2 & =\square_{1,2}\,, \\[0.2cm] a^2+(p^2\pm q^2)^2 & =\square_{1,2}\,, \\[0.2cm] (p^2\pm q^2)^2-a^2 & =(r^2\pm s^2)^2. \end{align*}


title: Decomposability of Nonnegative r-Potent Matrices

abstract: We consider nonnegative r-potent matrices with finite dimensions and study their decomposability. We derive the precise conditions under which an r-potent matrix is decomposable. We further determine a general structure for the r-potent matrices based on their decomposability. Finally, we establish that semigroups of r-potent matrices are also decomposable.


title: A Reconfigurable Vector Instruction Processor for Accelerating a Convection Parametrization Model on FPGAs

abstract: High Performance Computing (HPC) platforms allow scientists to model computationally intensive algorithms. HPC clusters increasingly use General-Purpose Graphics Processing Units (GPGPUs) as accelerators; FPGAs provide an attractive alternative to GPGPUs for use as co-processors, but they are still far from being mainstream due to a number of challenges faced when using FPGA-based platforms. Our research aims to make FPGA-based high performance computing more accessible to the scientific community. In this work we present the results of investigating the acceleration of a particular atmospheric model, Flexpart, on FPGAs. We focus on accelerating the most computationally intensive kernel from this model. The key contribution of our work is the architectural exploration we undertook to arrive at a solution that best exploits the parallelism available in the legacy code, and is also convenient to program, so that eventually the compilation of high-level legacy code to our architecture can be fully automated. We present the three different types of architecture, comparing their resource utilization and performance, and propose that an architecture where there are a number of computational cores, each built along the lines of a vector instruction processor, works best in this particular scenario, and is a promising candidate for a generic FPGA-based platform for scientific computation. We also present the results of experiments done with various configuration parameters of the proposed architecture, to show its utility in adapting to a range of scientific applications.


title: Isomorphism classes of $k$-involutions of algebraic groups of type $E_6$

abstract: Automorphisms of order $2$ are studied in order to understand generalized symmetric spaces. The groups of type $E_6$ we consider here can be realized as both the group of linear maps that leave a certain determinant invariant, and also as the identity component of the automorphism group of a class of structurable algebras known as Brown algebras. We will classify the $k$-involutions of these groups of type $E_6$ using aspects of both descriptions, and give detailed descriptions of representatives over certain fields including algebraically closed fields, $\mathbb{R}$, $\mathbb{F}_p$, and $\mathbb{Q}_p$.


title: The Nataf-Beta Random Field Classifier: An Extension of the Beta Conjugate Prior to Classification Problems

abstract: This paper presents the Nataf-Beta Random Field Classifier, a discriminative approach that extends the applicability of the Beta conjugate prior to classification problems. The approach's key feature is to model the probability of a class conditional on attribute values as a random field whose marginals are Beta distributed, and where the parameters of marginals are themselves described by random fields. Although the classification accuracy of the approach proposed does not statistically outperform the best accuracies reported in the literature, it ranks among the top tier for the six benchmark datasets tested. The Nataf-Beta Random Field Classifier is suited as a general purpose classification approach for real-continuous and real-integer attribute value problems.


title: A Two-Level Approach to Large Mixed-Integer Programs with Application to Cogeneration in Energy-Efficient Buildings

abstract: We study a two-stage mixed-integer linear program (MILP) with more than 1 million binary variables in the second stage. We develop a two-level approach by constructing a semi-coarse model (coarsened with respect to variables) and a coarse model (coarsened with respect to both variables and constraints). We coarsen binary variables by selecting a small number of pre-specified daily on/off profiles. We aggregate constraints by partitioning them into groups and summing over each group. With an appropriate choice of coarsened profiles, the semi-coarse model is guaranteed to find a feasible solution of the original problem and hence provides an upper bound on the optimal solution. We show that solving a sequence of coarse models converges to the same upper bound with proven finite steps. This is achieved by adding violated constraints to coarse models until all constraints in the semi-coarse model are satisfied. We demonstrate the effectiveness of our approach in cogeneration for buildings. The coarsened models allow us to obtain good approximate solutions at a fraction of the time required by solving the original problem. Extensive numerical experiments show that the two-level approach scales to large problems that are beyond the capacity of state-of-the-art commercial MILP solvers.


title: Finite-temperature scaling of spin correlations in an experimental realization of the one-dimensional Ising quantum critical point

abstract: Inelastic neutron scattering is used to study the finite-temperature scaling behavior of the local dynamic structure factor in the quasi-one-dimensional quantum antiferromagnet NTENP ($\text{Ni(N,N'-bis(3-aminopropyl)propane-1,3-diamine)(}\mu\text{-NO}_2\text{)ClO}_4$), at its field-induced Ising quantum critical point. The validity and the limitations of the theoretically predicted scaling relations are tested.


title: Persistence and permanence for a class of functional differential equations with infinite delay

abstract: The paper deals with a class of cooperative functional differential equations (FDEs) with infinite delay, for which sufficient conditions for persistence and permanence are established. Here, the persistence refers to all solutions with initial conditions that are positive, continuous and bounded. The present method applies to a very broad class of abstract systems of FDEs with infinite delay, both autonomous and non-autonomous, which include many important models used in mathematical biology. Moreover, the hypotheses imposed are in general very easy to check. The results are illustrated with some selected examples.


title: Non-existence of time-periodic vacuum spacetimes

abstract: We prove that smooth asymptotically flat solutions to the Einstein vacuum equations which are assumed to be periodic in time, are in fact stationary in a neighborhood of infinity. Our result applies under physically relevant regularity assumptions purely at the level of the initial data. In particular, our work removes the assumption of analyticity up to null infinity in [Bicak, Scholtz, and Tod; 2010]. The proof relies on extending a suitably constructed "candidate" Killing vector field from null infinity, via Carleman-type estimates obtained in [Alexakis, Schlue, Shao; 2013].


title: Involutive automorphisms of $N_\circ^\circ$ groups of finite Morley rank

abstract: We classify a large class of small groups of finite Morley rank: $N_\circ^\circ$-groups which are the infinite analogues of Thompson's $N$-groups. More precisely, we constrain the $2$-structure of groups of finite Morley rank containing a definable, normal, non-soluble, $N_\circ^\circ$-subgroup.


title: A Posteriori Error Estimator for a Front-Fixing Finite Difference Scheme for American Options

abstract: For the numerical solution of the American option valuation problem, we provide a script written in MATLAB implementing an explicit finite difference scheme. Our main contribute is the definition of a posteriori error estimator for the American options pricing which is based on Richardson's extrapolation theory. This error estimator allows us to find a suitable grid where the computed solution, both the option price field variable and the free boundary position, verify a prefixed error tolerance.


title: Random-projection ensemble classification

abstract: We introduce a very general method for high-dimensional classification, based on careful combination of the results of applying an arbitrary base classifier to random projections of the feature vectors into a lower-dimensional space. In one special case that we study in detail, the random projections are divided into disjoint groups, and within each group we select the projection yielding the smallest estimate of the test error. Our random projection ensemble classifier then aggregates the results of applying the base classifier on the selected projections, with a data-driven voting threshold to determine the final assignment. Our theoretical results elucidate the effect on performance of increasing the number of projections. Moreover, under a boundary condition implied by the sufficient dimension reduction assumption, we show that the test excess risk of the random projection ensemble classifier can be controlled by terms that do not depend on the original data dimension and a term that becomes negligible as the number of projections increases. The classifier is also compared empirically with several other popular high-dimensional classifiers via an extensive simulation study, which reveals its excellent finite-sample performance.


title: Structural Learning of Diverse Ranking

abstract: Relevance and diversity are both crucial criteria for an effective search system. In this paper, we propose a unified learning framework for simultaneously optimizing both relevance and diversity. Specifically, the problem is formalized as a structural learning framework optimizing Diversity-Correlated Evaluation Measures (DCEM), such as ERR-IA, a-NDCG and NRBP. Within this framework, the discriminant function is defined to be a bi-criteria objective maximizing the sum of the relevance scores and dissimilarities (or diversity) among the documents. Relevance and diversity features are utilized to define the relevance scores and dissimilarities, respectively. Compared with traditional methods, the advantages of our approach lie in that: (1) Directly optimizing DCEM as the loss function is more fundamental for the task; (2) Our framework does not rely on explicit diversity information such as subtopics, thus is more adaptive to real application; (3) The representation of diversity as the feature-based scoring function is more flexible to incorporate rich diversity-based features into the learning framework. Extensive experiments on the public TREC datasets show that our approach significantly outperforms state-of-the-art diversification approaches, which validate the above advantages.


title: Fire in the sky: The southern lights in Indigenous oral traditions

abstract: Parts of Australia have been privileged to see dazzling lights in the night sky as the Aurora Australis (known as the southern lights) puts on a show this year. Aurorae are significant in Australian Indigenous astronomical traditions. Aboriginal people associate aurorae with fire, death, blood, and omens, sharing many similarities with Native American communities.


title: Ultrametric fixed points in reduced axiomatic systems

abstract: The Brezis-Browder ordering principle [Advances Math., 21 (1976), 355-364] is used to get a proof, in the reduced axiomatic system (ZF-AC+DC), of a fixed point result [in the complete axiomatic system (ZF)] over Cantor complete ultrametric spaces due to Petalas and Vidalis [Proc. Amer. Math. Soc., 118 (1993), 819-821].


title: Testing Closeness With Unequal Sized Samples

abstract: We consider the problem of closeness testing for two discrete distributions in the practically relevant setting of \emph{unequal} sized samples drawn from each of them. Specifically, given a target error parameter $\varepsilon > 0$, $m_1$ independent draws from an unknown distribution $p,$ and $m_2$ draws from an unknown distribution $q$, we describe a test for distinguishing the case that $p=q$ from the case that $||p-q||_1 \geq \varepsilon$. If $p$ and $q$ are supported on at most $n$ elements, then our test is successful with high probability provided $m_1\geq n^{2/3}/\varepsilon^{4/3}$ and $m_2 = \Omega(\max\{\frac{n}{\sqrt m_1\varepsilon^2}, \frac{\sqrt n}{\varepsilon^2}\});$ we show that this tradeoff is optimal throughout this range, to constant factors. These results extend the recent work of Chan et al. who established the sample complexity when the two samples have equal sizes, and tightens the results of Acharya et al. by polynomials factors in both $n$ and $\varepsilon$. As a consequence, we obtain an algorithm for estimating the mixing time of a Markov chain on $n$ states up to a $\log n$ factor that uses $\tilde{O}(n^{3/2} \tau_{mix})$ queries to a "next node" oracle, improving upon the $\tilde{O}(n^{5/3}\tau_{mix})$ query algorithm of Batu et al. Finally, we note that the core of our testing algorithm is a relatively simple statistic that seems to perform well in practice, both on synthetic data and on natural language data.

categories: cs.LG cs.IT math.IT math.ST stat.ML stat.TH


title: The Decay of the Standard Model Higgs after Inflation

abstract: We study the nonperturbative dynamics of the Standard Model (SM) after inflation, in the regime where the SM is decoupled from (or weakly coupled to) the inflationary sector. We use classical lattice simulations in an expanding box in (3+1) dimensions, modeling the SM gauge interactions with both global and Abelian-Higgs analogue scenarios. We consider different post-inflationary expansion rates. During inflation, the Higgs forms a condensate, which starts oscillating soon after inflation ends. Via nonperturbative effects, the oscillations lead to a fast decay of the Higgs into the SM species, transferring most of the energy into $Z$ and $W^{\pm}$ bosons. All species are initially excited far away from equilibrium, but their interactions lead them into a stationary stage, with exact equipartition among the different energy components. From there on the system eventually reaches equilibrium. We have characterized in detail, in the different expansion histories considered, the evolution of the Higgs and of its dominant decay products, until equipartition is established. We provide a useful mapping between simulations with different parameters, from where we derive a master formula for the Higgs decay time, as a function of the coupling constants, Higgs initial amplitude and postinflationary expansion rate.


title: Electrically Engineered Band Gap in Two-Dimensional Ge, Sn, and Pb: A First-Principles and Tight-Binding Approach

abstract: First-principles calculations were performed to investigate the electronic structure of two-dimensional (2-D) Ge, Sn, and Pb without and with the presence of an external electric field in combination with spin-orbit coupling. Tight-binding calculations based on four orbitals per atom and an effective single orbital are presented to match with the results obtained from first-principles calculations. In particular, the electronic band structure and the band splitting are investigated with both models. Moreover, the simple $k\cdot p$ model is also considered in order to understand the band splitting in the presence of an external electric field and spin-orbit coupling. A large splitting is obtained, which is expected to be useful for spintronic devices. The fair agreement between the first-principle, $k\cdot p$ model, and tight-binding approaches leads to a table of parameters for future tight-binding studies on hexagonal 2-D nanostructures. By using the tight binding parameters, the transport properties of typical 0-D triangular quantum dots between two semi-infinite electrodes in the presence of spin-orbit coupling are addressed.


title: Magnetism of Ta Dichalcogenide Monolayers Tuned by Strain and Hydrogenation

abstract: The effects of strain and hydrogenation on the electronic and magnetic properties of monolayers of Ta based dichalcogenides (TaX2; X = S, Se, Te) are investigated using density-functional theo-ry. We predict a complex scenario of strain-dependent magnetic phase transitions involving par-amagnetic, ferromagnetic, and modulated antiferromagnetic states. Covering one of the two chalcogenide surfaces with hydrogen switches the antiferromagnetic/nonmagnetic TaX2 mono-layers to a semiconductor. Our research opens new pathways towards the manipulation of mag-netic properties for future optoelectronics and spintronics applications.


title: Soliton models for thick branes

abstract: In this work, we present new soliton solutions for thick branes in $4+1$ dimensions. In particular, we consider brane models based on the sine-Gordon ($SG$), $\varphi^{4}$ and $\varphi^{6}$ scalar fields, which have broken $Z_{2}$ symmetry in some cases, and are responsible for supporting and stabilizing the thick branes. The origin of the symmetry breaking in these models resides in the fact that the modified scalar field potential may have non-degenerate vacuua. These vacuua determine the cosmological constant on both sides of the brane. We also study the geodesic equations along the fifth dimension, in order to explore the particle motion in the neighbourhood of the brane. Furthermore, we examine the stability of the thick branes, by determining the sign of the $w^2$ term in the expansion of the potential for the resulting Schrodinger-like equation, where $w$ is the 5-dimensional coordinate. It turns out that the $\varphi^4$ brane is stable, while there are unstable modes for certain ranges of the model parameters in the SG and $\varphi^6$ branes.


title: Upper limits on DM annihilation cross sections from the first AMS-02 antiproton data

abstract: The first measurement on the antiproton to proton ratio made by the AMS-02 collaboration agrees with the expection from conventional cosmic-ray secondaries in the kinetic energy range $\sim 10-100$ GeV, which can be turned into stringent upper limits on the dark matter (DM) annihilation cross sections above $\sim 300$ GeV. Using the GALPROP code, we derive the upper limits in various propagation models and DM profiles. We show that in the "conventional" propagation model, for the $q\bar q$, $b\bar b$, and $WW$ final states, the constraints can be more stringent than that derived from the recent Ferm-LAT gamma-ray data on the dwarf spheroidal satellite galaxies. Making use of the typical minimal, median and maximal models obtained from a previous global fit, we show that the variation of the upper limits is around a factor of five. The possibility of DM contribution to the high energy $\bar p/p$ data is discussed.


title: Analysis of events with $b$-jets and a pair of leptons of the same charge in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector

abstract: An analysis is presented of events containing jets including at least one $b$-tagged jet, sizeable missing transverse momentum, and at least two leptons including a pair of the same electric charge, with the scalar sum of the jet and lepton transverse momenta being large. A data sample with an integrated luminosity of 20.3 fb$^{-1}$ of $pp$ collisions at $\sqrt{s} = 8$ TeV recorded by the ATLAS detector at the Large Hadron Collider is used. Standard Model processes rarely produce these final states, but there are several models of physics beyond the Standard Model that predict an enhanced rate of production of such events; the ones considered here are production of vector-like quarks, enhanced four-top-quark production, pair production of chiral $b^\prime$-quarks, and production of two positively charged top quarks. Eleven signal regions are defined; subsets of these regions are combined when searching for each class of models. In the three signal regions primarily sensitive to positively charged top quark pair production, the data yield is consistent with the background expectation. There are more data events than expected from background in the set of eight signal regions defined for searching for vector-like quarks and chiral $b^\prime$-quarks, but the significance of the discrepancy is less than two standard deviations. The discrepancy reaches 2.5 standard deviations in the set of five signal regions defined for searching for four-top-quark production. The results are used to set 95% CL limits on various models.


title: Level Lines of Gaussian Free Field II: Whole-Plane GFF

abstract: We study the level lines of GFF starting from interior points. We show that the level line of GFF starting from an interior point turns out to be a sequence of level loops. The sequence of level loops satisfies "target-independent" property. All sequences of level loops starting from interior points give a tree-structure of the plane. We also introduce the continuum exploration process of GFF starting from interior. The continuum exploration process of whole-plane GFF satisfies "reversibility".


title: A perturbative theory for Brownian vortexes

abstract: Brownian vortexes are stochastic machines that use static non-conservative force fields to bias random thermal fluctuations into steadily circulating currents. The archetype for this class of systems is a colloidal sphere in an optical tweezer. Trapped near the focus of a strongly converging beam of light, the particle is displaced by random thermal kicks into the nonconservative part of the optical force field arising from radiation pressure, which then biases its diffusion. Assuming the particle remains localized within the trap, its time-averaged trajectory traces out a toroidal vortex. Unlike trivial Brownian vortexes, such as the biased Brownian pendulum, which circulate preferentially in the direction of the bias, the general Brownian vortex can change direction and even topology in response to temperature changes. Here we introduce a theory based on a perturbative expansion of the Fokker-Planck equation for weak non-conservative driving. The first-order solution takes the form of a modified Boltzmann relation and accounts for the rich phenomenology observed in experiments on micrometer-scale colloidal spheres in optical tweezers.


title: Influence of disordered porous media in the anomalous properties of a simple water model

abstract: The thermodynamic, dynamic and structural behavior of a water-like system confined in a matrix is analyzed for increasing confining geometries. The liquid is modeled by a two dimensional associating lattice gas model that exhibits density and diffusion anomalies, in similarity to the anomalies present in liquid water. The matrix is a triangular lattice in which fixed obstacles impose restrictions to the occupation of the particles. We show that obstacules shortens all lines, including the phase coexistence, the critical and the anomalous lines. The inclusion of a very dense matrix not only suppress the anomalies but also the liquid-liquid critical point.


title: The small and large lags of the elastic and anelastic tides. The virtual identity of two rheophysical theories

abstract: The aim of this letter is to discuss the virtual identity of two recent tidal theories: the creep tide theory of Ferraz-Mello (Cel. Mech. Dyn. Astron. 116, 109, 2013) and the Maxwell model developed by Correia et al. (Astron. Astrophys. 571, A50, 2014). It includes the discussion of the basic equations of the theories, which, in both cases, include an elastic and an anelastic component, and shows that the basic equations of the two theories are equivalent and differ by only a numerical factor in the anelastic tide. It also includes a discussion of the lags: the lag of the full tide (geodetic), dominated by the elastic component, and the phase of the anelastic tide. In rotating rocky bodies not trapped in a spin-orbit resonance (e.g., the Earth) the geodetic lag is close to zero and the phase of the semi-diurnal argument in the anelastic tide is close to 90 degrees. The results obtained from combining tidal solutions from satellite tracking data and from Topex/Poseidon satellite altimeter data, by Ray et al., are extended to determine the phase of the semi-diurnal argument in the Earth's anelastic tide as sigma_0=89.80 \pm 0.05 degrees.


title: Perspectives on gravity-induced radiative processes in astrophysics

abstract: Single-vertex Feynman diagrams represent the dominant contribution to physical processes, but are frequently forbidden kinematically. This is changed when the particles involved propagate in a gravitational background and acquire an effective mass. Procedures are introduced that allow the calculation of lowest order diagrams, their corresponding transition probabilities, emission powers and spectra to all orders in the metric deviation, for particles of any spin propagating in gravitational fields described by any metric. Physical properties of the "space-time medium" are also discussed. It is shown in particular that a small dissipation term in the particle wave equations can trigger a strong back-reaction that introduces resonances in the radiative process and affects the resulting gravitational background.


title: The Araucaria Project: The First-Overtone Classical Cepheid in the Eclipsing System OGLE-LMC-CEP-2532

abstract: We present here the first spectroscopic and photometric analysis of the double-lined eclipsing binary containing the classical, first-overtone Cepheid OGLE-LMC-CEP-2532 (MACHO 81.8997.87). The system has an orbital period of 800 days and the Cepheid is pulsating with a period of 2.035 days. Using spectroscopic data from three high-class telescopes and photometry from three surveys spanning 7500 days we are able to derive the dynamical masses for both stars with an accuracy better than 3%. This makes the Cepheid in this system one of a few classical Cepheids with an accurate dynamical mass determination (M_1=3.90 +/- 0.10 M_sun). The companion is probably slightly less massive (3.82 +/- 0.10 M_sun), but may have the same mass within errors (M_2/M_1= 0.981 +/- 0.015). The system has an age of about 185 million years and the Cepheid is in a more advanced evolutionary stage. For the first time precise parameters are derived for both stars in this system. Due to the lack of the secondary eclipse for many years not much was known about the Cepheid's companion. In our analysis we used extra information from the pulsations and the orbital solution from the radial velocity curve. The best model predicts a grazing secondary eclipse shallower than 1 mmag, hence undetectable in the data, about 370 days after the primary eclipse. The dynamical mass obtained here is the most accurate known for a first-overtone Cepheid and will contribute to the solution of the Cepheid mass discrepancy problem.


title: Quaternionic Second-Order Freeness and the Fluctuations of Large Symplectically Invariant Random Matrices

abstract: We present a definition for second-order freeness in the quaternionic case. We demonstrate that this definition on a second-order probability space is asymptotically satisfied by independent symplectically invariant quaternionic matrices. This definition is different from the natural definition for complex and real second-order probability spaces, those motivated by the asymptotic behaviour of unitarily invariant and orthogonally invariant random matrices respectively. Most notably, because the quaternionic trace does not have the cyclic property of a trace over a commutative field, the asymmetries which appear in the multi-matrix context result in an asymmetric contribution from the terms which appear symmetrically in the complex and real cases.


title: Characterization of the Inner Knot of the Crab: The Site of the Gamma-ray Flares?

abstract: One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.


title: 6d String Chains

abstract: We consider bound states of strings which arise in 6d (1,0) SCFTs that are realized in F-theory in terms of linear chains of spheres with negative self-intersections 1,2, and 4. These include the strings associated to N small E8 instantons, as well as the ones associated to M5 branes probing A and D type singularities in M-theory or D5 branes probing ADE singularities in Type IIB string theory. We find that these bound states of strings admit (0,4) supersymmetric quiver descriptions and show how one can compute their elliptic genera.


title: MISO Broadcast Channel with Hybrid CSIT: Beyond Two Users

abstract: We study the impact of heterogeneity of channel-state-information available at the transmitters (CSIT) on the capacity of broadcast channels with a multiple-antenna transmitter and $k$ single-antenna receivers (MISO BC). In particular, we consider the $k$-user MISO BC, where the CSIT with respect to each receiver can be either instantaneous/perfect, delayed, or not available; and we study the impact of this heterogeneity of CSIT on the degrees-of-freedom (DoF) of such network. We first focus on the $3$-user MISO BC; and we completely characterize the DoF region for all possible heterogeneous CSIT configurations, assuming linear encoding strategies at the transmitters. The result shows that the state-of-the-art achievable schemes in the literature are indeed sum-DoF optimal, when restricted to linear encoding schemes. To prove the result, we develop a novel bound, called Interference Decomposition Bound, which provides a lower bound on the interference dimension at a receiver which supplies delayed CSIT based on the average dimension of constituents of that interference, thereby decomposing the interference into its individual components. Furthermore, we extend our outer bound on the DoF region to the general $k$-user MISO BC, and demonstrate that it leads to an approximate characterization of linear sum-DoF to within an additive gap of $0.5$ for a broad range of CSIT configurations. Moreover, for the special case where only one receiver supplies delayed CSIT, we completely characterize the linear sum-DoF.


title: On the readability of overlap digraphs

abstract: We introduce the graph parameter readability and study it as a function of the number of vertices in a graph. Given a digraph D, an injective overlap labeling assigns a unique string to each vertex such that there is an arc from x to y if and only if x properly overlaps y. The readability of D is the minimum string length for which an injective overlap labeling exists. In applications that utilize overlap digraphs (e.g., in bioinformatics), readability reflects the length of the strings from which the overlap digraph is constructed. We study the asymptotic behaviour of readability by casting it in purely graph theoretic terms (without any reference to strings). We prove upper and lower bounds on readability for certain graph families and general graphs


title: Quantum Coding with Finite Resources

abstract: The quantum capacity of a memoryless channel is often used as a single figure of merit to characterize its ability to transmit quantum information coherently. The capacity determines the maximal rate at which we can code reliably over asymptotically many uses of the channel. We argue that this asymptotic treatment is insufficient to the point of being irrelevant in the quantum setting where decoherence severely limits our ability to manipulate large quantum systems in the encoder and decoder. For all practical purposes we should instead focus on the trade-off between three parameters: the rate of the code, the number of coherent uses of the channel, and the fidelity of the transmission. The aim is then to specify the region determined by allowed combinations of these parameters. Towards this goal, we find approximate and exact characterizations of the region of allowed triplets for the qubit dephasing channel and for the erasure channel with classical post-processing assistance. In each case the region is parametrized by a second channel parameter, the quantum channel dispersion. In the process we also develop several general inner (achievable) and outer (converse) bounds on the coding region that are valid for all finite-dimensional quantum channels and can be computed efficiently. Applied to the depolarizing channel, this allows us to determine a lower bound on the number of coherent uses of the channel necessary to witness super-additivity of the coherent information.


title: Critical analysis and remedy of switching failures in straintronic logic using Bennett clocking in the presence of thermal fluctuations

abstract: Straintronic logic is a promising platform for beyond Moore's law computing. Using Bennett clocking mechanism, information can propagate through an array of strain-mediated multiferroic nanomagnets exploiting the dipolar coupling between the magnets without having to physically interconnect them. Here we perform a critical analysis of switching failures, i.e., error in information propagation due to thermal fluctuations through a chain of such straintronic devices. We solved stochastic Landau-Lifshitz-Gilbert equation considering room-temperature thermal perturbations and show that magnetization switching may fail due to inherent magnetization dynamics accompanied by thermally broadened switching delay distribution. Avenues available to circumvent such issue are proposed.


title: Optical aperture synthesis with electronically connected telescopes

abstract: Highest resolution imaging in astronomy is achieved by interferometry, connecting telescopes over increasingly longer distances, and at successively shorter wavelengths. Here, we present the first diffraction-limited images in visual light, produced by an array of independent optical telescopes, connected electronically only, with no optical links between them. With an array of small telescopes, second-order optical coherence of the sources is measured through intensity interferometry over 180 baselines between pairs of telescopes, and two-dimensional images reconstructed. The technique aims at diffraction-limited optical aperture synthesis over kilometre-long baselines to reach resolutions showing details on stellar surfaces and perhaps even the silhouettes of transiting exoplanets. Intensity interferometry circumvents problems of atmospheric turbulence that constrain ordinary interferometry. Since the electronic signal can be copied, many baselines can be built up between dispersed telescopes, and over long distances. Using arrays of air Cherenkov telescopes, this should enable the optical equivalent of interferometric arrays currently operating at radio wavelengths.


title: Diffuse coronae in cosmological simulations of Milky Way-sized galaxies

abstract: We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.


title: Discovering Higgs boson pair production through rare final states at a 100 TeV collider

abstract: We consider Higgs boson pair production at a future proton collider with centre-of-mass energy of 100 TeV, focusing on rare final states that include a bottom-anti-bottom quark pair and multiple isolated leptons: $hh \rightarrow (b\bar{b}) + n \ell + X$, $n = \{2,4\}$, $X = \{ E_T^\mathrm{miss}, \gamma, -\}$. We construct experimental search strategies for observing the process through these channels and make suggestions on the desired requirements for the detector design of the future collider.


title: Static Self-Forces in a Five-Dimensional Black Hole Spacetime

abstract: We obtain the electric field and scalar field for a static point charge in closed form in the 5D Schwarzschild-Tangherlini black hole spacetime. We then compute the static self-force in each of these cases by assuming that the appropriate singular field is a 4D Hadamard Green's function on the constant time Riemannian slice. It is well known that the Hadamard Green's function involves an arbitrary regular biscalar $W_{0}(x,x')$, whose coincidence limit $w(x)$ appears in the expression for the self-force. We develop an axiomatic approach to reduce this arbitrary function to a single arbitrary dimensionless coefficient. We show that in the context of this approach to regularization, the self-force does not depend on any undetermined length-scale and need not depend on the internal structure of the charge.


title: Cosmological Tests of Modified Gravity

abstract: We review recent progress in the construction of modified gravity models as alternatives to dark energy as well as the development of cosmological tests of gravity. Einstein's theory of General Relativity (GR) has been tested accurately within the local universe i.e. the Solar System, but this leaves the possibility open that it is not a good description of gravity at the largest scales in the Universe. This being said, the standard model of cosmology assumes GR on all scales. In 1998, astronomers made the surprising discovery that the expansion of the Universe is accelerating, not slowing down. This late-time acceleration of the Universe has become the most challenging problem in theoretical physics. Within the framework of GR, the acceleration would originate from an unknown dark energy. Alternatively, it could be that there is no dark energy and GR itself is in error on cosmological scales. In this review, we first give an overview of recent developments in modified gravity theories including $f(R)$ gravity, braneworld gravity, Horndeski theory and massive/bigravity theory. We then focus on common properties these models share, such as screening mechanisms they use to evade the stringent Solar System tests. Once armed with a theoretical knowledge of modified gravity models, we move on to discuss how we can test modifications of gravity on cosmological scales. We present tests of gravity using linear cosmological perturbations and review the latest constraints on deviations from the standard $\Lambda$CDM model. Since screening mechanisms leave distinct signatures in the non-linear structure formation, we also review novel astrophysical tests of gravity using clusters, dwarf galaxies and stars.


title: Signatures of planets and protoplanets in the Galactic center: a clue to understand the G2 cloud?

abstract: Several hundred young stars lie in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might capture planets orbiting these stars, and bring them onto nearly radial orbits. The same fate might occur to planetary embryos (PEs), i.e. protoplanets born from gravitational instabilities in protoplanetary disks. In this paper, we investigate the emission properties of rogue planets and PEs in the Galactic center. In particular, we study the effects of photoevaporation, caused by the ultraviolet background. Rogue planets can hardly be detected by current or forthcoming facilities, unless they are tidally disrupted and accrete onto the SMBH. In contrast, photoevaporation of PEs (especially if the PE is being tidally stripped) might lead to a recombination rate as high as ~10^45 s^-1, corresponding to a Brackett-gamma luminosity ~10^31 erg s^-1, very similar to the observed luminosity of the dusty object G2. We critically discuss the possibility that G2 is a rogue PE, and the major uncertainties of this model.


title: Kepler-539: a young extrasolar system with two giant planets on wide orbits and in gravitational interaction

abstract: We confirm the planetary nature of Kepler-539b (aka Kepler object of interest K00372.01), a giant transiting exoplanet orbiting a solar-analogue G2 V star. The mass of Kepler-539b was accurately derived thanks to a series of precise radial velocity measurements obtained with the CAFE spectrograph mounted on the CAHA 2.2m telescope. A simultaneous fit of the radial-velocity data and Kepler photometry revealed that Kepler-539b is a dense Jupiter-like planet with a mass of Mp = 0.97 Mjup and a radius of Rp = 0.747 Rjup, making a complete circular revolution around its parent star in 125.6 days. The semi-major axis of the orbit is roughly 0.5 au, implying that the planet is at roughly 0.45 au from the habitable zone. By analysing the mid-transit times of the 12 transit events of Kepler-539b recorded by the Kepler spacecraft, we found a clear modulated transit time variation (TTV), which is attributable to the presence of a planet c in a wider orbit. The few timings available do not allow us to precisely estimate the properties of Kepler-539c and our analysis suggests that it has a mass between 1.2 and 3.6 Mjup, revolving on a very eccentric orbit (0.4<e<0.6) with a period larger than 1000 days. The high eccentricity of planet c is the probable cause of the TTV modulation of planet b. The analysis of the CAFE spectra revealed a relatively high photospheric lithium content, A(Li)=2.48 dex, which, together with both a gyrochronological and isochronal analysis, suggests that the parent star is relatively young.


title: Observable fractions of core-collapse supernova light curves brightened by binary companions

abstract: Many core-collapse supernova progenitors are presumed to be in binary systems. If a star explodes in a binary system, the early supernova light curve can be brightened by the collision of the supernova ejecta with the companion star. The early brightening can be observed when the observer is in the direction of the hole created by the collision. Based on a population synthesis model, we estimate the fractions of core-collapse supernovae in which the light-curve brightening by the collision can be observed. We find that 0.19% of core-collapse supernova light curves can be observed with the collisional brightening. Type Ibc supernova light curves are more likely to be brightened by the collision (0.53%) because of the high fraction of the progenitors being in binary systems and their proximity to the companion stars. Type II and IIb supernova light curves are less affected (~1e-3% and ~1e-2%, respectively). Although the early, slow light-curve declines of some Type IIb and Ibc supernovae are argued to be caused by the collision with the companion star (e.g. SN 2008D), the small expected fraction, as well as the unrealistically small separation required, disfavour the argument. The future transient survey by the Large Synoptic Survey Telescope is expected to detect ~10 Type Ibc supernovae with the early collisional brightening per year, and they will be able to provide information on supernova progenitors in binary systems.


title: Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field

abstract: Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We compare the predicted peculiar velocities from 2M++ to Tully-Fisher and SNe peculiar velocities. We find a value of $\beta^* \equiv \Omega_{\rm{m}}^{0.55}/b^* = 0.431 \pm 0.021$, suggesting $\Omega_{\rm{m}}^{0.55}\sigma_{\rm{8,lin}} = 0.401 \pm 0.024$, in good agreement with other probes. The predicted peculiar velocity of the Local Group arising from the 2M++ volume alone is $540 \pm 40$ km/s, towards $l = 268 \pm 4, b= 38 \pm 6$, only $10^\circ$ out of alignment with the Cosmic Microwave Background dipole. To account for velocity contributions arising from sources outside the 2M++ volume, we fit simultaneously for $\beta^*$ and an external bulk flow in our analysis. We find that an external bulk flow is preferred at the 5.1$\sigma$ level, and the best fit has a velocity of $159\pm23$ km/s towards $l=304 \pm 11, b = 6 \pm 13$. Finally, the predicted bulk flow of a 50 Mpc/h Gaussian-weighted volume centred on the Local Group is $230 \pm 30$ km/s, in the direction $l=293\pm 8, b = 14 \pm 10$, in agreement with predictions from $\Lambda$CDM.


title: Nucleon electromagnetic form factors in two-flavour QCD

abstract: We present results for the nucleon electromagnetic form factors, including the momentum transfer dependence and derived quantities (charge radii and magnetic moment). The analysis is performed using O(a) improved Wilson fermions in Nf=2 QCD measured on the CLS ensembles. Particular focus is placed on a systematic evaluation of the influence of excited states in three-point correlation functions, which lead to a biased evaluation, if not accounted for correctly. We argue that the use of summed operator insertions and fit ans\"atze including excited states allow us to suppress and control this effect. We employ a novel method to perform joint chiral and continuum extrapolations, by fitting the form factors directly to the expressions of covariant baryonic chiral effective field theory. The final results for the charge radii and magnetic moment from our lattice calculations include, for the first time, a full error budget. We find that our estimates are compatible with experimental results within their overall uncertainties.


title: Transient Sulfate Aerosols as a Signature of Exoplanet Volcanism

abstract: Geological activity is thought to be important for the origin of life and for maintaining planetary habitability. We show that transient sulfate aerosols could be a signature of exoplanet volcanism, and therefore a geologically active world. A detection of transient aerosols, if linked to volcanism, could thus aid in habitability evaluations of the exoplanet. On Earth, subduction-induced explosive eruptions inject SO2 directly into the stratosphere, leading to the formation of sulfate aerosols with lifetimes of months to years. We demonstrate that the rapid increase and gradual decrease in sulfate aerosol loading associated with these eruptions may be detectable in transit transmission spectra with future large-aperture telescopes, such as the James Webb Space Telescope (JWST) and European Extremely-Large Telescope (E-ELT) for a planetary system at a distance of 10 pc, assuming an Earth-like atmosphere, bulk composition, and size. Specifically, we find that a S/N of 12.1 and 7.1 could be achieved with E-ELT (assuming photon-limited noise) for an Earth-analog orbiting a Sun-like star and M5V star, respectively, even without multiple transits binned together. We propose that the detection of this transient signal would strongly suggest an exoplanet volcanic eruption, if potential false positives such as dust storms or bolide impacts can be ruled out. Furthermore, because scenarios exist in which O2 can form abiotically in the absence of volcanic activity, a detection of transient aerosols that can be linked to volcanism, along with a detection of O2, would be a more robust biosignature than O2 alone.


title: The Hunt for the Rest of the Higgs Bosons

abstract: We assess the current state of searches at the LHC for additional Higgs bosons in light of both direct limits and indirect bounds coming from coupling measurements of the Standard Model-like Higgs boson. Given current constraints, we identify and study three LHC searches that are critical components of a comprehensive program to investigate extended electroweak symmetry breaking sectors: production of a heavy scalar or pseudoscalar with decay to $t \bar t$; $b \bar b$ and $t \bar t$ associated production of a heavy scalar or pseudoscalar with decay to invisible final states; and $t \bar b$ associated production of a charged Higgs with decay to $\bar t b$. Systematic experimental searches in these channels would contribute to robust coverage of the possible single production modes of additional heavy Higgs bosons.


title: A regularity result for the nonlocal Fokker-Planck equation with Ornstein-Uhlenbeck drift

abstract: Despite there are numerous theoretical studies of stochastic differential equations with a symmetric $\alpha$-stable L\'evy noise, very few regularity results exist in the case of $0<\alpha\leq1$. In this paper, we study the fractional Fokker-Planck equation with Ornstein-Uhlenbeck drift, and prove that there exists a unique solution, which is $C^\infty$ in space for $t>0$ when $\alpha\in (0, 2]$.


title: Likelihood-Ratio Ranking Statistic for Compact Binary Coalescence Candidates with Rate Estimation

abstract: We present a new likelihood-ratio ranking statistic for use in searches for gravitational waves from the inspiral and merger of compact object binaries. Expanding on previous work, the ranking statistic incorporates a model for the correlations in the signal-to-noise ratios with which signals will be seen in a network of ground-based antennas while retaining an algebraic procedure for mapping ranking statistic values to false-alarm probability. Additionally, the ranking statistic enables the implementation of a rigorous signal rate estimation technique. We implement the ranking statistic and demonstrate its use including signal rate estimation in an analysis of a simulated signal population in simulated noise.


title: Comparison of DAMA/LIBRA and DM ICE Results using Information Theory to Rule out Dark Matter Claims

abstract: We study the details of the DAMA/LIBRA results and compare those with the recent published DM Ice results of ICE Cube. In various recent papers, it was shown that the 40K peak on DAMA/LIBRA data leaves no room for a Dark Matter signal in the bulk of the data. Using Information Theory for the different types of detection environments, we show that annual variation calculations and the DM Ice data reinforce the claims that the DAMA/LIBRA detector is not observing Dark Matter WIMPs.


title: Magnetocaloric Effect of Gadolinium at Adiabatic and Quasi-Isothermal Conditions in High Magnetic Fields

abstract: High cooling power of magnetocaloric refrigeration can be achieved only at large amounts of heat, which can be transferred in one cycle from cold end to hot end at quasi-isothermal conditions. The simple experimental method for direct measurement of the transferred heat from material with magnetocaloric effect (MCE) to massive nonmagnetic block at quasi-isothermal conditions was proposed. The vacuum calorimeter was designed for the simultaneous measurements of MCE both at adiabatic conditions (deltaT) and quasi-isothermal conditions (deltaQ) in the magnetic fields of Bitter coil magnet. This calorimeter was tested on samples of pure polycrystalline Gd with direct MCE. The maximal obtained values were deltaT = 17.7 K and deltaQ = 5900 J/kg at initial temperature 20 0C in magnetic field 140 kOe.


title: The Khavinson-Shapiro Conjecture for the Bergman Projection in One and Several Complex Variables

abstract: We reveal a complex analogue to a result about polynomial solutions to the Dirichlet Problem on ellipsoids in $\mathbb{R}^n$ by showing that the Bergman projection on any ellipsoid in $\mathbb{C}^n$ is such that the projection of any polynomial function of degree at most $N$ is a holomorphic polynomial function of degree at most $N$. The discussion is motivated by a connection between the Bergman projection and the Khavinson-Shapiro conjecture in $\mathbb{C}$. We also relate the Khavinson-Shapiro conjecture to polyharmonic Bergman projections in $\mathbb{R}^n$ by showing that these projections take polynomials to polynomials on ellipsoids.


title: Consistent Learning by Composite Proximal Thresholding

abstract: We investigate the modeling and the numerical solution of machine learning problems with prediction functions which are linear combinations of elements of a possibly infinite-dimensional dictionary. We propose a novel flexible composite regularization model, which makes it possible to incorporate various priors on the coefficients of the prediction function, including sparsity and hard constraints. We show that the estimators obtained by minimizing the regularized empirical risk are consistent in a statistical sense, and we design an error-tolerant composite proximal thresholding algorithm for computing such estimators. New results on the asymptotic behavior of the proximal forward-backward splitting method are derived and exploited to establish the convergence properties of the proposed algorithm. In particular, our method features a $o(1/m)$ convergence rate in objective values.


title: Height variation of the vector magnetic field in solar spicules

abstract: Proving the magnetic configuration of solar spicules has hitherto been difficult due to the lack of spatial resolution and image stability during off-limb ground-based observations. We report spectropolarimetric observations of spicules taken in the He I 1083 nm spectral region with the Tenerife Infrared Polarimeter II at the German Vacuum Tower Telescope of the Observatorio del Teide (Tenerife; Canary Islands; Spain). The data provide the variation with geometrical height of the Stokes I, Q, U, and V profiles whose encoded information allows the determination of the magnetic field vector by means of the HAZEL inversion code. The inferred results show that the average magnetic field strength at the base of solar spicules is about 80 gauss and then it decreases rapidly with height to about 30 gauss at a height of 3000 km above the visible solar surface. Moreover, the magnetic field vector is close to vertical at the base of the chromosphere and has mid inclinations (about 50 degree) above 2 Mm height.


title: Emergent hyperuniformity in periodically-driven emulsions

abstract: We report the emergence of large-scale hyperuniformity in microfluidic emulsions. Upon periodic driving confined emulsions undergo a first-order transition from a reversible to an irreversible dynamics. We evidence that this dynamical transition is accompanied by structural changes at all scales yielding macroscopic yet finite hyperuniform structures. Numerical simulations are performed to single out the very ingredients responsible for the suppression of density fluctuations. We show that as opposed to equilibrium systems the long-range nature of the hydrodynamic interactions are not required for the formation of hyperuniform patterns, thereby suggesting a robust relation between reversibility and hyperuniformity which should hold in a broad class of periodically driven materials.


title: Pairwise-like models for non-Markovian epidemics on networks

abstract: In this letter, a generalization of pairwise models to non-Markovian epidemics on networks is presented. For the case of infectious periods of fixed length, the resulting pairwise model is a system of delay differential equations (DDEs), which shows excellent agreement with results based on explicit stochastic simulations of non-Markovian epidemics on networks. Furthermore, we analytically compute a new R0-like threshold quantity and an implicit analytical relation between this and the final epidemic size. In addition we show that the pairwise model and the analytic calculations can be generalized in terms of integro-differential equations to any distribution of the infectious period, and we illustrate this by presenting a closed form expression for the final epidemic size. By showing the rigorous mathematical link between non-Markovian network epidemics and pairwise DDEs, we provide the framework for a deeper and more rigorous understanding of the impact of non-Markovian dynamics with explicit results for final epidemic size and threshold quantities.


title: The Influence of Malloc Placement on TSX Hardware Transactional Memory

abstract: The hardware transactional memory (HTM) implementation in Intel's i7-4770 "Haswell" processor tracks the transactional read-set in the L1 (level-1), L2 (level-2) and L3 (level-3) caches and the write-set in the L1 cache. Displacement or eviction of read-set entries from the cache hierarchy or write-set entries from the L1 results in abort. We show that the placement policies of dynamic storage allocators -- such as those found in common "malloc" implementations -- can influence the L1 conflict miss rate in the L1. Conflict misses -- sometimes called mapping misses -- arise because of less than ideal associativity and represent imbalanced distribution of active memory blocks over the set of available L1 indices. Under transactional execution conflict misses may manifest as aborts, representing wasted or futile effort instead of a simple stall as would occur in normal execution mode. Furthermore, when HTM is used for transactional lock elision (TLE), persistent aborts arising from conflict misses can force the offending thread through the so-called "slow path". The slow path is undesirable as the thread must acquire the lock and run the critical section in normal execution mode, precluding the concurrent execution of threads in the "fast path" that monitor that same lock and run their critical sections in transactional mode. For a given lock, multiple threads can concurrently use the transactional fast path, but at most one thread can use the non-transactional slow path at any given time. Threads in the slow path preclude safe concurrent fast path execution. Aborts rising from placement policies and L1 index imbalance can thus result in loss of concurrency and reduced aggregate throughput.


title: Cool White Dwarfs Selection with Pan-STARRS Proper Motions

abstract: The use of Reduced Proper Motion in identifying isolated white dwarfs has long been used as a proxy for the absolute magnitude in a population with known kinematics. This, however, introduces a proper motion detection limit on top of the existing photometric limit. How the survey volume is hampered by this extra parameter is discussed in Hambly et al. 2012. In this work, we discuss some robust outlier rejection methods in order to minimise the proper motion limit and hence maximise the survey volume. The generalised volume, corrected for the distance of the Sun from the Galactic Plane, is integrated explicitly.


title: Commensurability Classes of Fake Quadrics

abstract: A fake quadric is a smooth projective surface that has the same rational cohomology as a smooth quadric surface but is not biholomorphic to one. We provide an explicit classification of all irreducible fake quadrics according to the commensurability class of their fundamental group. To accomplish this task, we develop a number of new techniques that explicitly bound the arithmetic invariants of a fake quadric and more generally of an arithmetic manifold of bounded volume arising from a form of SL_2 over a number field.


title: A maximum volume density estimator generalized over a proper motion-limited sample

abstract: The traditional Schmidt density estimator has been proven to be unbiased and effective in a magnitude-limited sample. Previously, efforts have been made to generalize it for populations with non-uniform density and proper motion-limited cases. This work shows that the then-good assumptions for a proper motion-limited sample are no longer sufficient to cope with modern data. Populations with larger differences in the kinematics as compared to the local standard of rest are most severely affected. We show that this systematic bias can be removed by treating the discovery fraction inseparable from the generalized maximum volume integrand. The treatment can be applied to any proper motion-limited sample with good knowledge of the kinematics. This work demonstrates the method through application to a mock catalogue of a white dwarf-only solar neighbourhood for various scenarios and compared against the traditional treatment using a survey with Pan-STARRS-like characteristics.


title: Gender Issues & Information Communication Technology for Development (ICT4D): Prospects and Challenges for Women in Nigeria

abstract: Information and Communication Technology is a compendium of interrelated applications, products and services that can either deepen or alter gender equality. It has successfully transformed education, businesses, healthcare, entertainment, politics and good governance within the Global North; providing equitable access to developmental framework driven by ICTs. However, in drafting the core developmental objectives of sustainable development, there has been considerable gender digital divide limiting women access to resources based on their gender, ethnicity, socio-cultural bias and the rights to utilize such resources for development. In realization of the United Nation Millennium Development Goals within the Global South specifically Nigeria, women are often marginalized or excluded from ICT policy drafting and imbalance associated with ICT. This paper identifies and evaluates gender issues and information communication technology, with focus on the challenges and prospects for women empowerment in Nigeria. The study critically examined research literatures and conducted research survey on the prospects and challenges of promoting gender equality and women empowerment through ICTs; and identifies policy implication for Nigeria. The research survey used a random sampling technique with a target sample size of eighty respondents. Data gathered from the questionnaire was analyzed using Statistical Package for Social Science version 19, and the result was presented using ANOVA, and descriptive analysis. The study reveal gender inclusiveness in policy drafting as a key driver for socio-economic development, improved healthcare and women empowerment in Nigeria. We recommend a deliberate ICT policy that attract and encourages women participation in ICT developmental framework.


title: Piecewise equidistant meshes for quasilinear turning point problems: Technical report

abstract: A class of quasilinear singularly perturbed boundary value problems with a turning point of attractive type is considered. The problems are solved numerically by a finite-difference scheme on a special discretization mesh which is dense near the turning point. The scheme is a combination of the standard central and midpoint schemes and is practically second-order accurate. Pointwise accuracy is uniform in the perturbation parameter and, moreover, L1 errors decrease when the perturbation parameter tends to 0. This is achieved by the use of meshes which generalize the piecewise equidistant Shishkin mesh. Two particular types of meshes are considered and compared.


title: Performance Evaluation of Machine Learning Algorithms in Post-operative Life Expectancy in the Lung Cancer Patients

abstract: The nature of clinical data makes it difficult to quickly select, tune and apply machine learning algorithms to clinical prognosis. As a result, a lot of time is spent searching for the most appropriate machine learning algorithms applicable in clinical prognosis that contains either binary-valued or multi-valued attributes. The study set out to identify and evaluate the performance of machine learning classification schemes applied in clinical prognosis of post-operative life expectancy in the lung cancer patients. Multilayer Perceptron, J48, and the Naive Bayes algorithms were used to train and test models on Thoracic Surgery datasets obtained from the University of California Irvine machine learning repository. Stratified 10-fold cross-validation was used to evaluate baseline performance accuracy of the classifiers. The comparative analysis shows that multilayer perceptron performed best with classification accuracy of 82.3%, J48 came out second with classification accuracy of 81.8%, and Naive Bayes came out the worst with classification accuracy of 74.4%. The quality and outcome of the chosen machine learning algorithms depends on the ingenuity of the clinical miner.


title: A uniform nilsequence Wiener-Wintner theorem for bilinear ergodic averages

abstract: We show that a $k$-linear pointwise ergodic theorem on an ergodic measure-preserving system implies a uniform $k$-linear nilsequence Wiener-Wintner theorem on that system. The assumption is known to hold for arbitrary systems and $k=2$ (due to Bourgain) and for distal systems and arbitrary $k$ (due to Huang, Shao, and Ye).


title: On equivariant asymptotic dimension

abstract: The work discusses equivariant asymptotic dimension (also known as "wide equivariant covers", "$N$-$\mathcal F$-amenability" or "amenability dimension", and "$d$-BLR condition") and its generalisation, transfer reducibility, which are versions of asymptotic dimension invented for the proofs of the Farrell--Jones and Borel conjectures. We prove that groups of null equivariant asymptotic dimension are exactly virtually cyclic groups. We show that a covering of the boundary always extends to a covering of the whole compactification. We provide a number of characterisations of equivariant asymptotic dimension in the general setting of homotopy actions, including equivariant counterparts of classical characterisations of asymptotic dimension. Finally, we strengthen the result of Mole and Rueping (arXiv:1308.2799) about equivariant refinements from finite groups to infinite groups.


title: Van der Waals forces from first principles for periodic systems: Application to graphene-water interactions

abstract: We extend the method of Silvestrelli [P. L. Silvestrelli, J. Chem. Phys. 139, 054106 (2013)] to approximate long-range van der Waals interactions at the density functional theory level based on maximally localized Wannier functions combined with the quantum harmonic oscillator model, to periodic systems. Applying this scheme to study London dispersion forces between graphene and water layers, we demonstrate that collective many-body effects beyond simple additive pair-wise interactions are essential to accurately describe van der Waals forces.

categories: cond-mat.mtrl-sci cond-mat.dis-nn physics.chem-ph physics.comp-ph quant-ph


title: A Faster FPTAS for the Unbounded Knapsack Problem

abstract: The Unbounded Knapsack Problem (UKP) is a well-known variant of the famous 0-1 Knapsack Problem (0-1 KP). In contrast to 0-1 KP, an arbitrary number of copies of every item can be taken in UKP. Since UKP is NP-hard, fully polynomial time approximation schemes (FPTAS) are of great interest. Such algorithms find a solution arbitrarily close to the optimum $\mathrm{OPT}(I)$, i.e. of value at least $(1-\varepsilon) \mathrm{OPT}(I)$ for $\varepsilon > 0$, and have a running time polynomial in the input length and $\frac{1}{\varepsilon}$. For over thirty years, the best FPTAS was due to Lawler with a running time in $O(n + \frac{1}{\varepsilon^3})$ and a space complexity in $O(n + \frac{1}{\varepsilon^2})$, where $n$ is the number of knapsack items. We present an improved FPTAS with a running time in $O(n + \frac{1}{\varepsilon^2} \log^3 \frac{1}{\varepsilon})$ and a space bound in $O(n + \frac{1}{\varepsilon} \log^2 \frac{1}{\varepsilon})$. This directly improves the running time of the fastest known approximation schemes for Bin Packing and Strip Packing, which have to approximately solve UKP instances as subproblems.


title: Biometrics for Child Vaccination and Welfare: Persistence of Fingerprint Recognition for Infants and Toddlers

abstract: With a number of emerging applications requiring biometric recognition of children (e.g., tracking child vaccination schedules, identifying missing children and preventing newborn baby swaps in hospitals), investigating the temporal stability of biometric recognition accuracy for children is important. The persistence of recognition accuracy of three of the most commonly used biometric traits (fingerprints, face and iris) has been investigated for adults. However, persistence of biometric recognition accuracy has not been studied systematically for children in the age group of 0-4 years. Given that very young children are often uncooperative and do not comprehend or follow instructions, in our opinion, among all biometric modalities, fingerprints are the most viable for recognizing children. This is primarily because it is easier to capture fingerprints of young children compared to other biometric traits, e.g., iris, where a child needs to stare directly towards the camera to initiate iris capture. In this report, we detail our initiative to investigate the persistence of fingerprint recognition for children in the age group of 0-4 years. Based on preliminary results obtained for the data collected in the first phase of our study, use of fingerprints for recognition of 0-4 year-old children appears promising.


title: Defect-induced phase transition in the asymmetric simple exclusion process

abstract: We reconsider the long-standing question of the critical defect hopping rate $r_c$ in the one-dimensional totally asymmetric exclusion process (TASEP) with a slow bond (defect). For $r< r_c$ a phase separated state is observed due to queuing at the defect site whereas for $r\geq r_c$ the defect site has only local effects on the stationary state of the homogeneous system. Mean-field theory predicts $r_c=1$ (when hopping rates outside the defect bond are equal to 1) but numerical investigations seem to indicate $r_c \approx 0.80(2)$. Here we improve the numerics to show that $r_c > 0.99$ and give strong evidence that indeed $r_c=1$ as predicted by mean-field theory, and anticipated by recent theoretical findings.


title: SymPix: A spherical grid for efficient sampling of rotationally invariant operators

abstract: We present SymPix, a special-purpose spherical grid optimized for efficient sampling of rotationally invariant linear operators. This grid is conceptually similar to the Gauss-Legendre (GL) grid, aligning sample points with iso-latitude rings located on Legendre polynomial zeros. Unlike the GL grid, however, the number of grid points per ring varies as a function of latitude, avoiding expensive over-sampling near the poles and ensuring nearly equal sky area per grid point. The ratio between the number of grid points in two neighbouring rings is required to be a low-order rational number (3, 2, 1, 4/3, 5/4 or 6/5) to maintain a high degree of symmetries. Our main motivation for this grid is to solve linear systems using multi-grid methods, and to construct efficient preconditioners through pixel-space sampling of the linear operator in question. The GL grid is not suitable for these purposes due to its massive over-sampling near the poles, leading to nearly degenerate linear systems, while HEALPix, another commonly used spherical grid, exhibits few symmetries, and is therefore computationally inefficient for these purposes. As a benchmark and representative example, we compute a preconditioner for a linear system with both HEALPix and SymPix that involves the operator $D + B^T N^{-1} B$, where $B$ and $D$ may be described as both local and rotationally invariant operators, and $N$ is diagonal in pixel domain. For a bandwidth limit of $\ell_\text{max}=3000$, we find that SymPix, due to its higher number of internal symmetries, yields average speed-ups of 360 and 23 for $B^T N^{-1} B$ and $D$, respectively, relative to HEALPix.


title: Information without rolling dice

abstract: The deterministic notions of capacity and entropy are studied in the context of communication and storage of information using square-integrable, bandlimited signals subject to perturbation. The $(\epsilon,\delta)$-capacity, that extends the Kolmogorov $\epsilon$-capacity to packing sets of overlap at most $\delta$, is introduced and compared to the Shannon capacity. The functional form of the results indicates that in both Kolmogorov and Shannon's settings, capacity and entropy grow linearly with the number of degrees of freedom, but only logarithmically with the signal to noise ratio. This basic insight transcends the details of the stochastic or deterministic description of the information-theoretic model. For $\delta=0$, the analysis leads to new bounds on the Kolmogorov $\epsilon$-capacity, and to a tight asymptotic expression of the Kolmogorov $\epsilon$-entropy of bandlimited signals. A deterministic notion of error exponent is introduced. Applications of the theory are briefly discussed.


title: Ground States for a nonlinear Schr\"odinger system with sublinear coupling terms

abstract: We study the existence of ground states for the coupled Schr\"odinger system \begin{equation} \left\{\begin{array}{lll} \displaystyle -\Delta u_i+\lambda_i u_i= \mu_i |u_i|^{2q-2}u_i+\sum_{j\neq i}b_{ij} |u_j|^q|u_i|^{q-2}u_i \\ u_i\in H^1(\mathbb{R}^n), \quad i=1,\ldots, d, \end{array}\right. \end{equation} $n\geq 1$, for $\lambda_i,\mu_i >0$, $b_{ij}=b_{ji}>0$ (the so-called "symmetric attractive case") and $1<q<n/(n-2)^+$. We prove the existence of a nonnegative ground state $(u_1^*,\ldots,u_d^*)$ with $u_i^*$ radially decreasing. Moreover we show that, for $1<q<2$, such ground states are positive in all dimensions and for all values of the parameters.


abstract: A detailed Faddeev-Jackiw quantization of an Abelian and non-Abelian exotic action for gravity in three dimensions is performed. We obtain for the theories under study the constraints, the gauge transformations, the generalized Faddeev-Jackiw brackets and we perform the counting of physical degrees of freedom. In addition, we compare our results with those found in the literature where the canonical analysis is developed, in particular, we show that both the generalized Faddeev-Jackiw brackets and Dirac's brackets coincide to each other. Finally we discuss some remarks and prospects.


title: Kra\'skiewicz-Pragacz modules and Ringel duality

abstract: Kra\'skiewicz and Pragacz introduced representations of the upper-triangular Lie algebras whose characters are Schubert polynomials. In a previous work the author studied the structure of Kra\'skiewicz-Pragacz modules using the theory of highest weight categories. From the results there, in particular we obtain a certain highest weight category whose standard modules are KP modules. In this paper we show that this highest weight category is self Ringel-dual: this leads to an interesting symmetry relation on Ext groups between KP modules. We also show that the tensor product operation on b-modules is compatible with Ringel duality functor.