title: Meta-food-chains as a many-layer epidemic process on networks

abstract: Notable recent works have focused on the multi-layer properties of coevolving diseases. We point out that very similar systems play an important role in population ecology. Specifically we study a meta food-web model that was recently proposed by Pillai et al. This model describes a network of species connected by feeding interactions, which spread over a network of spatial patches. Focusing on the essential case, where the network of feeding interactions is a chain, we develop an analytical approach for the computation of the degree distributions of colonized spatial patches for the different species in the chain. This framework allows us to address ecologically relevant questions. Considering configuration model ensembles of spatial networks, we find that there is an upper bound for the fraction of patches that a given species can occupy, which depends only on the networks mean degree. For a given mean degree there is then an optimal degree distribution that comes closest to the upper bound. Notably scale-free degree distributions perform worse than more homogeneous degree distributions if the mean degree is sufficiently high. Because species experience the underlying network differently the optimal degree distribution for one particular species is generally not the optimal distribution for the other species in the same food web. These results are of interest for conservation ecology, where, for instance, the task of selecting areas of old-growth forest to preserve in an agricultural landscape, amounts to the design of a patch network.


title: Minimal cosmography

abstract: The minimal requirement for cosmography - a nondynamical description of the universe - is a prescription for calculating null geodesics, and timelike geodesics as a function of their proper time. In this paper, we consider the most general linear connection compatible with homogeneity and isotropy, but not necessarily with a metric. A light-cone structure is assigned by choosing a set of geodesics representing light rays. This defines a "scale factor" and a local notion of distance, as that travelled by light in a given proper time interval. We find that the velocities and relativistic energies of free-falling bodies decrease in time as a consequence of cosmic expansion, but at a rate that can be different than that dictated by the usual metric framework. By extrapolating this behavior to photons redshift, we find that the latter is in principle independent of the "scale factor". Interestingly, redshift-distance relations and other standard geometric observables are modified in this extended framework, in a way that could be experimentally tested. An extremely tight constraint on the model, however, is represented by the blackbody-ness of the Cosmic Microwave Background. Finally, as a check, we also consider the effects of a non-metric connection in a different set-up, namely, that of a static, spherically symmetric spacetime.


title: The relative wavelength independence of IR lags in active galactic nuclei: implications for the distribution of the hot dust

abstract: We show that, contrary to simple predictions, most AGNs show at best only a small increase of lags in the J, H, K, and L bands with increasing wavelength. We suggest that a possible cause of this near simultaneity of the variability from the near-IR to the mid-IR is that the hot dust is in a hollow bi-conical outflow of which we only see the near side. Although most AGNs show near simultaneity of IR variability, there was at least one epoch when NGC 4151 showed the sharply increasing IR lag with the increase of the wavelength. This behaviour might also be present in GQ Comae. We discuss these results briefly. The relative wavelength independence of IR lags simplifies the use of IR lags for estimating cosmological parameters.


title: Multi-Threaded Dense Linear Algebra Libraries for Low-Power Asymmetric Multicore Processors

abstract: Dense linear algebra libraries, such as BLAS and LAPACK, provide a relevant collection of numerical tools for many scientific and engineering applications. While there exist high performance implementations of the BLAS (and LAPACK) functionality for many current multi-threaded architectures,the adaption of these libraries for asymmetric multicore processors (AMPs)is still pending. In this paper we address this challenge by developing an asymmetry-aware implementation of the BLAS, based on the BLIS framework, and tailored for AMPs equipped with two types of cores: fast/power hungry versus slow/energy efficient. For this purpose, we integrate coarse-grain and fine-grain parallelization strategies into the library routines which, respectively, dynamically distribute the workload between the two core types and statically repartition this work among the cores of the same type. Our results on an ARM big.LITTLE processor embedded in the Exynos 5422 SoC, using the asymmetry-aware version of the BLAS and a plain migration of the legacy version of LAPACK, experimentally assess the benefits, limitations, and potential of this approach.


title: Theoretical and Numerical Analyses of a Slit-Masked Chicane for Modulated Bunch Generation

abstract: Density modulations on electron beams can improve machine performance of beam-driven accelerators and FELs with resonance beam-wave coupling. The beam modulation is studied with a masked chicane by the analytic model and simulations with the beam parameters of the Fermilab Accelerator Science and Technology (FAST) facility. With the chicane design parameters (bending angle of 18{\deg}, bending radius of 0.95 m and R56 ~ - 0.19 m) and a nominal beam of 3 ps bunch length, the analytic model showed that a slit-mask with slit period 900 {\mu}m and aperture width 300 {\mu}m induces a modulation of bunch-to-bunch spacing ~100 {\mu}m to the bunch with 2.4% correlated energy spread. With the designed slit mask and a 3 ps bunch, particle-in-cell (PIC) simulations, including nonlinear energy distributions, space charge force, and coherent synchrotron radiation (CSR) effect, also result in beam modulation with bunch-to-bunch distance around 100 {\mu}m and a corresponding modulation frequency of 3 THz. The beam modulation has been extensively examined with three different beam conditions, 2.25 ps (0.25 nC), 3.25 ps (1 nC), and 4.75 ps (3.2 nC) by tracking code Elegant. The simulation analysis indicates that the sliced beam by the slit-mask with 3 ~ 6% correlated energy spread has modulation lengths about 187 {\mu}m (0.25 nC), 270 {\mu}m (1 nC) and 325 {\mu}m (3.2 nC). The theoretical and numerical data proved the capability of the designed masked chicane in producing modulated bunch train with micro-bunch length around 100 fs.


title: Minimal surfaces in the soliton surface approach

abstract: The main objective of this paper is to derive the Enneper-Weierstrass representation of minimal surfaces in $\mathbb{E}^3$ using the soliton surface approach. We exploit the Bryant-type representation of conformally parametrized surfaces in the hyperbolic space $H^3(\lambda)$ of curvature $-\lambda^2$, which can be interpreted as a 2 by 2 linear problem involving the spectral parameter $\lambda$. In the particular case of constant mean curvature-$\lambda$ surfaces a special limiting procedure $(\lambda\rightarrow 0)$, different from that of Umehara and Yamada [33], allows us to recover the Enneper-Weierstrass representation. Applying such a limiting procedure to the previously known cases, we obtain Sym-type formulas. Finally we exploit the relation between the Bryant representation of constant mean curvature-$\lambda$ surfaces and second-order linear ordinary differential equations. We illustrate this approach by the example of the error function equation.


title: Evidence for a Gapped Spin-Liquid Ground State in a Kagome Heisenberg Antiferromagnet

abstract: The kagome Heisenberg antiferromagnet is a leading candidate in the search for a spin system with a quantum spin-liquid ground state. The nature of its ground state remains a matter of great debate. We conducted 17-O single crystal NMR measurements of the S=1/2 kagome lattice in herbertsmithite ZnCu$_3$(OH)$_6$Cl$_2$, which is known to exhibit a spinon continuum in the spin excitation spectrum. We demonstrate that the intrinsic local spin susceptibility $\chi_{kagome}$ deduced from the 17-O NMR frequency shift asymptotes to zero below temperature T ~ 0.03 J, where J ~ 200 K is the Cu-Cu super-exchange interaction. Combined with the magnetic field dependence of $\chi_{kagome}$ we observed at low temperatures, these results imply that the kagome Heisenberg antiferromagnet has a spin-liquid ground state with a finite gap.


title: Methods of Class Field Theory to Separate Logics over Finite Residue Classes and Circuit Complexity

abstract: Separations among the first order logic ${\cal R}ing(0,+,*)$ of finite residue class rings, its extensions with generalized quantifiers, and in the presence of a built-in order are shown, using algebraic methods from class field theory. These methods include classification of spectra of sentences over finite residue classes as systems of congruences, and the study of their $h$-densities over the set of all prime numbers, for various functions $h$ on the natural numbers. Over ordered structures the logic of finite residue class rings and extensions are known to capture DLOGTIME-uniform circuit complexity classes ranging from $AC^0$ to $TC^0$. Separating these circuit complexity classes is directly related to classifying the $h$-density of spectra of sentences in the corresponding logics of finite residue classes. We further give general conditions under which a logic over the finite residue class rings has a sentence whose spectrum has no $h$-density. One application of this result is that in ${\cal R}ing(0,+,*,<) + M$, the logic of finite residue class rings with built-in order and extended with the majority quantifier $M$, there are sentences whose spectrum have no exponential density.


title: Optimal Non-Asymptotic Lower Bound on the Minimax Regret of Learning with Expert Advice

abstract: We prove non-asymptotic lower bounds on the expectation of the maximum of $d$ independent Gaussian variables and the expectation of the maximum of $d$ independent symmetric random walks. Both lower bounds recover the optimal leading constant in the limit. A simple application of the lower bound for random walks is an (asymptotically optimal) non-asymptotic lower bound on the minimax regret of online learning with expert advice.


title: The $\mathbb{Z}_2^n$ Dirac-Dunkl operator and a higher rank Bannai-Ito algebra

abstract: The kernel of the $\mathbb{Z}_2^{n}$ Dirac-Dunkl operator is examined. The symmetry algebra $\mathcal{A}_{n}$ of the associated Dirac-Dunkl equation on $\mathbb{S}^{n-1}$ is determined and is seen to correspond to a higher rank generalization of the Bannai-Ito algebra. A basis for the polynomial null-solutions of the Dirac-Dunkl operator is constructed. The basis elements are joint eigenfunctions of a maximal commutative subalgebra of $\mathcal{A}_{n}$ and are given explicitly in terms of Jacobi polynomials. The symmetry algebra is shown to act irreducibly on this basis via raising/lowering operators. A scalar realization of $\mathcal{A}_{n}$ is proposed.


title: A Bayesian Consistent Dual Ensemble Kalman Filter for State-Parameter Estimation in Subsurface Hydrology

abstract: Ensemble Kalman filtering (EnKF) is an efficient approach to addressing uncertainties in subsurface groundwater models. The EnKF sequentially integrates field data into simulation models to obtain a better characterization of the model's state and parameters. These are generally estimated following joint and dual filtering strategies, in which, at each assimilation cycle, a forecast step by the model is followed by an update step with incoming observations. The Joint-EnKF directly updates the augmented state-parameter vector while the Dual-EnKF employs two separate filters, first estimating the parameters and then estimating the state based on the updated parameters. In this paper, we reverse the order of the forecast-update steps following the one-step-ahead (OSA) smoothing formulation of the Bayesian filtering problem, based on which we propose a new dual EnKF scheme, the Dual-EnKF$_{\rm OSA}$. Compared to the Dual-EnKF, this introduces a new update step to the state in a fully consistent Bayesian framework, which is shown to enhance the performance of the dual filtering approach without any significant increase in the computational cost. Numerical experiments are conducted with a two-dimensional synthetic groundwater aquifer model to assess the performance and robustness of the proposed Dual-EnKF$_{\rm OSA}$, and to evaluate its results against those of the Joint- and Dual-EnKFs. The proposed scheme is able to successfully recover both the hydraulic head and the aquifer conductivity, further providing reliable estimates of their uncertainties. Compared with the standard Joint- and Dual-EnKFs, the proposed scheme is found more robust to different assimilation settings, such as the spatial and temporal distribution of the observations, and the level of noise in the data. Based on our experimental setups, it yields up to 25% more accurate state and parameters estimates.


title: Irrational Base Counting

abstract: We will provide algorithmic implementation with proofs of existence and uniqueness for the Absolute and Alternating Ostrowski Numeration Systems.


title: Stabilization of a Rigid Body Payload with Multiple Cooperative Quadrotors

abstract: This paper presents the full dynamics and control of arbitrary number of quadrotor unmanned aerial vehicles (UAV) transporting a rigid body. The rigid body is connected to the quadrotors via flexible cables where each flexible cable is modeled as a system of arbitrary number of serially connected links. It is shown that a coordinate-free form of equations of motion can be derived for the complete model without any simplicity assumptions that commonly appear in other literature, according to Lagrangian mechanics on a manifold. A geometric nonlinear controller is presented to transport the rigid body to a fixed desired position while aligning all of the links along the vertical direction. A rigorous mathematical stability proof is given and the desirable features of the proposed controller are illustrated by numerical examples and experimental results.


title: Production of a forward J/psi and a backward jet at the LHC

abstract: We study the production of a forward $J/\psi$ meson and a backward jet with a large rapidity separation at the LHC using the BFKL formalism. We compare the predictions given by the Non Relativistic QCD (NRQCD) approach to charmonium prediction and by the Color Evaporation Model. In NRQCD, we find that the $^3 S_1^{\, 8}$ part of the onium wavefunction is completely dominating the process. NRQCD and the color evaporation model give similar results, although a discrepancy seems to appear as the value of the transverse momenta of the charmonium and of the jet decrease.


title: Optimizing Damper Connectors for Adjacent Buildings

abstract: Many theoretical and experimental studies have used heuristic methods to investigate the dynamic behaviour of the passive coupling of adjacent structures. However, few papers have used optimization techniques with guaranteed convergence in order to increase the efficiency of the passive coupling of adjacent structures. In this paper, the combined problem of optimal arrangement and mechanical properties of dampers placed between two adjacent buildings is considered. A new bi-level optimization approach is presented. The outer-loop of the approach optimizes damper configuration and is solved using the ``inserting dampers'' method, which was recently shown to be a very effective heuristic method. Under the assumption that the dampers have varying damper coefficients, the inner-loop finds the optimal damper coefficients by solving an $n$-dimensional optimization problem, where derivative information of the objective function is not available. Three different non-gradient methods are compared for solving the inner loop: a genetic algorithm (GA), the mesh adaptive direct search (MADS) algorithm, and the robust approximate gradient sampling (RAGS) algorithm. It is shown that by exploiting this new bi-level problem formulation, modern derivative free optimization techniques with guaranteed convergence (such as MADS and RAGS) can be used. The results indicate a great increase in the efficiency of the retrofitting system, as well as the existence of a threshold on the number of dampers inserted with respect to the efficiency of the retrofitting system.


title: Transcranial electrical stimulation: An introduction

abstract: The main objective of the electrical stimulation of the brain is to generate action potentials from the application of electromagnetic fields. Among the available techniques, transcranial electrical stimulation (TES) represents a popular method of administration that has the advantage of being non-invasive and economically more affordable. This article aims to briefly introduce the reader into the understanding of TES in terms of the physics involved as well as for some of the relevant results of studies applying this technique.


title: Moduli of Continuity for Viscosity Solutions

abstract: In this paper, we investigate the moduli of continuity for viscosity solutions of a wide class of nonsingular quasilinear evolution equations and also for the level set mean curvature flow, which is an example of singular degenerate equations. We prove that the modulus of continuity is a viscosity subsolution of some one dimensional equation. This work extends B. Andrews' recent result on moduli of continuity for smooth spatially periodic solutions.


title: Moduli of Continuity for Viscosity Solutions on Manifolds

abstract: We establish the estimates of modulus of continuity for viscosity solutions of nonlinear evolution equations on manifolds, extending previous work of B. Andrews and J. Clutterbuck for regular solutions on manifolds \cite{AC3} and the first author's recent work for viscosity solutions in Euclidean spaces \cite{me1}.


title: Accelerating Adaptive IDW Interpolation Algorithm on a Single GPU

abstract: This paper focuses on the design and implementing of GPU-accelerated Adaptive Inverse Distance Weighting (AIDW) interpolation algorithm. The AIDW is an improved version of the standard IDW, which can adaptively determine the power parameter according to the spatial points distribution pattern and achieve more accurate predictions than those by IDW. In this paper, we first present two versions of the GPU accelerated AIDW, the naive version without profiting from shared memory and the tiled version taking advantage of shared memory. We also implement the naive version and the tiled version using the data layouts, Structure of Arrays (AoS) and Array of aligned Structures (AoaS), on single and double precision. We then evaluate the performance of the GPU-accelerated AIDW by comparing it with its original CPU version. Experimental results show that: on single precision the naive version and the tiled version can achieve the speedups of approximately 270 and 400, respectively. In addition, on single precision the implementations using the layout SoA are always slightly faster than those using layout AoaS. However, on double precision, the speedup is only about 8; and we have also observed that: (1) there are no performance gains obtained from the tiled version against the naive version; and (2) the use of SoA and AoaS does not lead to significant differences in computational efficiency.


title: Streaming regularization parameter selection via stochastic gradient descent

abstract: We propose a framework to perform streaming covariance selection. Our approach employs regularization constraints where a time-varying sparsity parameter is iteratively estimated via stochastic gradient descent. This allows for the regularization parameter to be efficiently learnt in an online manner. The proposed framework is developed for linear regression models and extended to graphical models via neighbourhood selection. Under mild assumptions, we are able to obtain convergence results in a non-stochastic setting. The capabilities of such an approach are demonstrated using both synthetic data as well as neuroimaging data.


title: Efficient Functional-Based Adaptation for CFD Applications

abstract: Adjoint methods have gained popularity in recent years for driving adaptation procedures which aim to reduce error in solution functionals. While adjoint methods have been proven effective for functional-based adaptation, the practical implementation of an adjoint method can be quite burdensome since code developers constantly need to ensure and maintain a dual consistent discretization as updates are made. Also, since most engineering problems consider multiple functionals, an adjoint solution must be obtained for each functional of interest which can increase the overall computational cost significantly. In this paper, an alternative to adjoints is presented which uses a sparse approximate inverse of the Jacobian of the residual to obtain approximate adjoint sensitivities for functional-based adaptation indicators. Since the approximate inverse need only be computed once, it can be recycled for any number of functionals making the new approach more efficient than a conventional adjoint method. This new method for functional-based adaptation will be tested using the quasi-1D nozzle problem, and results are presented for functionals of integrated pressure and entropy.


title: Stop or Continue Data Collection: A Nonignorable Missing Data Approach for Continuous Variables

abstract: We present an approach to inform decisions about nonresponse follow-up sampling. The basic idea is (i) to create completed samples by imputing nonrespondents' data under various assumptions about the nonresponse mechanisms, (ii) take hypothetical samples of varying sizes from the completed samples, and (iii) compute and compare measures of accuracy and cost for different proposed sample sizes. As part of the methodology, we present a new approach for generating imputations for multivariate continuous data with nonignorable unit nonresponse. We fit mixtures of multivariate normal distributions to the respondents' data, and adjust the probabilities of the mixture components to generate nonrespondents' distributions with desired features. We illustrate the approaches using data from the 2007 U. S. Census of Manufactures.


title: Parallel projection methods for variational inequalities involving common fixed point problems

abstract: In this paper, we introduce two novel parallel projection methods for finding a solution of a system of variational inequalities which is also a common fixed point of a family of (asymptotically) $\kappa$ - strict pseudocontractive mappings. A technical extension in the proposed algorithms helps in computing practical numerical experiments when the number of subproblems is large. Some numerical examples are implemented to demonstrate the efficiency of parallel computations.


title: The warm ionized gas in CALIFA early-type galaxies: 2D emission-line patterns and kinematics for 32 galaxies

abstract: The morphological, spectroscopic and kinematical properties of the warm interstellar medium (wim) in early-type galaxies (ETGs) hold key observational constraints to nuclear activity and the buildup history of these massive, quiescent systems. High-quality integral field spectroscopy (IFS) data with a wide spectral and spatial coverage, such as those from the CALIFA survey, offer an unprecedented opportunity for advancing our understanding of the wim in ETGs. This article centers on a 2D investigation of the wim component in 32 nearby (<~150Mpc) ETGs from CALIFA, complementing a previous 1D analysis of the same sample (Papaderos et al. 2013; P13). We include here H\alpha\ intensity and equivalent width (EW) maps and radial profiles, diagnostic emission-line ratios, besides ionized-gas and stellar kinematics. This study is supplemented by \tau-ratio maps as an efficient means to quantify the role of photoionization by pAGB stars, as compared to other mechanisms (e.g., AGN, low-level star formation). Additionally, we extend the tentative classification proposed in P13 by the type i+, which is assigned to a subset of type i ETGs exhibiting ongoing low-level star-formation (SF) in their periphery. This finding along with faint traces of localized SF in the extranuclear component of several of our sample ETGs points to a non-negligible contribution by OB stars to the total ionizing budget. We also demonstrate that, at the typical emission-line detection threshold of ~2\AA in previous studies, most of the extranuclear wim emission in an ETG may evade detection, which could in turn prompt its classification as an entirely gas-devoid system. This study adds further observational evidence for a considerable heterogeneity among ETGs with regard to the physical properties and 2D kinematics of the wim component, and underscores the importance of IFS studies over their entire optical extent.


title: Quantum Memristors

abstract: Technology based on memristors, resistors with memory whose resistance depends on the history of the crossing charges, has lately enhanced the classical paradigm of computation with neuromorphic architectures. However, in contrast to the known quantized models of passive circuit elements, such as inductors, capacitors or resistors, the design and realization of a quantum memristor is still missing. Here, we introduce the concept of a quantum memristor as a quantum dissipative device, whose decoherence mechanism is controlled by a continuous-measurement feedback scheme, which accounts for the memory. Indeed, we provide numerical simulations showing that memory effects actually persist in the quantum regime. Our quantization method, specifically designed for superconducting circuits, may be extended to other quantum platforms, allowing for memristor-type constructions in different quantum technologies. The proposed quantum memristor is then a building block for neuromorphic quantum computation and quantum simulations of non-Markovian systems.


title: Curvature and torsion effects in the spin-current driven domain wall motion

abstract: The domain wall motion along a helix-shaped nanowire is studied for the case of spin-current driving via Bazaliy-Zhang-Li mechanism. The analysis is based on collective variable approach. Two new effects are ascertained: (i) the curvature results in appearance of the Walker limit for a uniaxial wire, (ii) the torsion results in effective shift of the nonadiabatic spin torque parameter $\beta$. The latter effect changes considerably the domain wall velocity and can result in negative domain wall mobility. This effect can be also used for an experimental determination of the nonadiabatic parameter $\beta$ and damping coefficient $\alpha$.


title: Exact conditions on the temperature dependence of density functionals

abstract: Universal exact conditions guided the construction of most ground-state density functional approximations in use today. We derive the relation between the entropy and Mermin free energy density functionals for thermal density functional theory. Both the entropy and sum of kinetic and electron-electron repulsion functionals are shown to be monotonically increasing with temperature, while the Mermin functional is concave downwards. Analogous relations are found for both exchange and correlation. The importance of these conditions is illustrated in two extremes: the Hubbard dimer and the uniform gas.


title: Implications Of The Higgs Discovery On Physics Beyond The Standard Model

abstract: In this thesis, we investigate the implications of the LHC Higgs data on different BSM scenarios. Since the data seem to agree with the SM expectations, any nonstandard couplings will be strongly constrained. First we investigate, in a model independent way, the constraints on the nonstandard Higgs couplings with the fermions and the vector bosons in view of high energy unitarity and the measured value of the Higgs to diphoton signal strength. Then we concentrate on a particular BSM scenario, namely, the two Higgs-doublet models (2HDMs). Consistency of the Higgs data with the corresponding SM predictions strongly motivates us to work in the {\em alignment limit}. In this limit, including the informations of Higgs mass and its SM-like nature, we find many new constraints on the nonstandard masses and $\tan\beta$. We also study the constraints on the charged scalar mass arising from the $h\to \gamma\gamma$ signal strength measurements and observe that the charged scalar does not necessarily decouple from the diphoton decay width. We then move on to some particular variants of 2HDMs, known as BGL models, and study the flavor constraints on these models. Here we find that lighter than conventionally allowed nonstandard scalars can successfully negotiate the stringent bounds coming from flavor physics data and can leave unconventional decay signatures that can be used as distinctive features of these models. We also analyze the stability and unitarity constraints in a three Higgs-doublet model (3HDM) with $S_3$ symmetry and find that there must be many more nonstandard particles below 1~TeV. We also observe that the nondecoupling feature of the charged scalar in the context of $h\to \gamma\gamma$ is not unique to the 2HDMs only, instead it is a general property of the multi doublet extensions of the SM with an exact discrete symmetry.


title: Evaluating Protein-protein Interaction Predictors with a Novel 3-Dimensional Metric

abstract: In order for the predicted interactions to be directly adopted by biologists, the ma- chine learning predictions have to be of high precision, regardless of recall. This aspect cannot be evaluated or numerically represented well by traditional metrics like accuracy, ROC, or precision-recall curve. In this work, we start from the alignment in sensitivity of ROC and recall of precision-recall curve, and propose an evaluation metric focusing on the ability of a model to be adopted by biologists. This metric evaluates the ability of a machine learning algorithm to predict only new interactions, meanwhile, it eliminates the influence of test dataset. In the experiment of evaluating different classifiers with a same data set and evaluating the same predictor with different datasets, our new metric fulfills the evaluation task of our interest while two widely recognized metrics, ROC and precision-recall curve fail the tasks for different reasons.


title: A Novel Visualization of the Geometry of Special Relativity

abstract: The mathematical treatment and graphical representation of Special Relativity (SR) are well established, yet carry deep implications that remain hard to visualize. This paper presents a new graphical interpretation of the geometry of SR that may, by complementing the standard works, aid the understanding of SR and its fundamental principles in a more intuitive way. From the axiom that the velocity of light remains constant to any inertial observer, the geodesic is presented as a line of constant angle on the complex plane across a set of diverging reference frames. The resultant curve is a logarithmic spiral, and this view of the geodesic is extended to illustrate the relativistic Doppler effect, time dilation, length contraction, the twin paradox, and relativistic radar distance in an original way, whilst retaining the essential mathematical relationships of SR. Using a computer-generated graphical representation of photon trajectories allows a visual comparison between the relativistic relationships and their classical counterparts, to visualize the consequences of SR as velocities become relativistic. The model may readily be extended to other situations, and may be found useful in presenting a fresh understanding of SR through geometric visualization.


title: Theory of packaged entangled states

abstract: The entangled states that include every physical properties of particles would be important for both theoretical and applied physics. However, the existence and properties of such entangled states are unclear at present. Here we theoretically show that a particle-antiparticle pair can form the so-called packaged entangled states which encapsulate all the necessary physical quantities for completely identifying the particles. The particles in the packaged entangled states are indeterminate and exhibit unusual properties. Thereafter, we discussed the possible applications of these new entangled states, i.e., the protocol for teleporting the entire quantum state of a particle (or an antiparticle) to an arbitrarily large distance without a classical channel, transfer of the new entangled states from a particle pair to another particle pair, and new interpretation to the matter-antimatter asymmetry of the observable universe.


title: The Poisson Gamma Belief Network

abstract: To infer a multilayer representation of high-dimensional count vectors, we propose the Poisson gamma belief network (PGBN) that factorizes each of its layers into the product of a connection weight matrix and the nonnegative real hidden units of the next layer. The PGBN's hidden layers are jointly trained with an upward-downward Gibbs sampler, each iteration of which upward samples Dirichlet distributed connection weight vectors starting from the first layer (bottom data layer), and then downward samples gamma distributed hidden units starting from the top hidden layer. The gamma-negative binomial process combined with a layer-wise training strategy allows the PGBN to infer the width of each layer given a fixed budget on the width of the first layer. The PGBN with a single hidden layer reduces to Poisson factor analysis. Example results on text analysis illustrate interesting relationships between the width of the first layer and the inferred network structure, and demonstrate that the PGBN, whose hidden units are imposed with correlated gamma priors, can add more layers to increase its performance gains over Poisson factor analysis, given the same limit on the width of the first layer.


title: Proposal to Search for Magnetically Charged Particles with Magnetic Charge 1e

abstract: A model for composite elementary Standard Model (SM) particles based upon magnetically bound vorton pairs, briefly introduced here, predicts the existence of a complete family of magnetically charged particles, as well as their neutral isotopic partners (all counterparts to the SM elementary particles), in which the lowest mass (charged) particle would be an electrically neutral stable lepton, but which carries a magnetic charge equivalent to 1e. This new particle, which we call a magneticon (a counterpart to the electron) would be pair produced at all e+e- colliders at an Ecm above twice its mass. In addition, PP and PPbar colliders should also be able to produce these new particles through the Drell-Yan process. To our knowledge, no monopole search experiment has been sensitive to such a low-charged magnetic monopole above a particle mass of about 5 GeV/c2. Hence, we propose that a search for such a stable particle of magnetic charge 1e should be undertaken. We have taken the ATLAS detector at the LHC as an example in which this search might be done. To this end, we have modeled the magnetic fields and muon trigger chambers of this detector. We show results from a simple Monte Carlo simulation program to indicate how these particles might look in the detector and describe how one might search for these new particles in the ATLAS data stream.


title: Cooperation and Competition among Energy Storages

abstract: We study competition and cooperation among a group of storage units. We show that as the number of storages increases, the profit of storages approaches zero under Nash competition. We propose two ways in which storages can achieve non-zero profit and show that they are optimal in the sense that storages achieve the maximum possible profit. The first is a decentralized approach in which storages are exposed to artificial cost functions that incentivize them to behavior as a coalition. No private information needs to be exchanged between the storages to calculate the artificial function. The second is a centralized approach in which an aggregator coordinates and splits profits with storages in order to achieve maximum profit. We use Nash's axiomatic bargaining problem to model and predict the profit split between aggregator and storages.


title: 1-Grothendieck $C(K)$ spaces

abstract: A Banach space is said to be Grothendieck if weak and weak$^*$ convergent sequences in the dual space coincide. This notion has been quantificated by H. Bendov\'{a}. She has proved that $\ell_\infty$ has the quantitative Grothendieck property, namely, it is 1-Grothendieck. Our aim is to show that Banach spaces from a certain wider class are 1-Grothendieck, precisely, $C(K)$ is 1-Grothendieck provided $K$ is a totally disconnected compact space such that its algebra of clopen subsets has the so called Subsequential completeness property.


title: Spherical Tropicalization

abstract: We extend tropicalization and tropical compactification of subvarieties of algebraic tori to subvarieties of spherical homogeneous spaces. Given a tropical compactification of a subvariety, we show that the support of the colored fan of the ambient spherical variety agrees with the tropicalization of the subvariety. The proof is based on our equivariant version of the attening by blow-up theorem. We provide many examples.


title: An Extended Frank-Wolfe Method with "In-Face" Directions, and its Application to Low-Rank Matrix Completion

abstract: Motivated principally by the low-rank matrix completion problem, we present an extension of the Frank-Wolfe method that is designed to induce near-optimal solutions on low-dimensional faces of the feasible region. This is accomplished by a new approach to generating ``in-face" directions at each iteration, as well as through new choice rules for selecting between in-face and ``regular" Frank-Wolfe steps. Our framework for generating in-face directions generalizes the notion of away-steps introduced by Wolfe. In particular, the in-face directions always keep the next iterate within the minimal face containing the current iterate. We present computational guarantees for the new method that trade off efficiency in computing near-optimal solutions with upper bounds on the dimension of minimal faces of iterates. We apply the new method to the matrix completion problem, where low-dimensional faces correspond to low-rank matrices. We present computational results that demonstrate the effectiveness of our methodological approach at producing nearly-optimal solutions of very low rank. On both artificial and real datasets, we demonstrate significant speed-ups in computing very low-rank nearly-optimal solutions as compared to either the Frank-Wolfe method or its traditional away-step variant.


title: Entanglement scaling of excited states in large one-dimensional many-body localized systems

abstract: We study the properties of excited states in one-dimensional many-body localized (MBL) systems using a matrix product state algorithm. First, the method is tested for a large disordered non-interacting system, where for comparison we compute a quasi-exact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of L~100 sites and large bond dimensions chi~1700, which allows us to quantitatively analyze the scaling behavior of the entanglement S in the system. The MBL phase is characterized by a logarithmic growth (L)~log(L) over a large scale separating the regimes where volume and area laws hold. We check the validity of the eigenstate thermalization hypothesis. Our results are consistent with the existence of a mobility edge.


title: Lower Bounds for Enumerative Counts of Positive-Genus Real Curves

abstract: We transform the positive-genus real Gromov-Witten invariants of many real-orientable symplectic threefolds into signed counts of curves. These integer invariants provide lower bounds for counts of real curves of a given genus that pass through conjugate pairs of constraints. We conclude with some implications and related conjectures for one- and two-partition Hodge integrals.


title: Norm-inflation results for the BBM equation

abstract: Considered here is the periodic initial-value probem for the regularized long-wave (BBM) equation \[u_t+u_x+uu_x-u_{xxt}=0.\] Adding to previous work in the literature, it is shown here that for any $s < 0$, there is smooth initial data that is small in the $L_2$-based Sobolev spaces $H^s$, but the solution emanating from it becomes arbitrarily large in arbitrarily small time. This so called {\it norm inflation} result has as a consequence the previously determined conclusion that this problem is ill-posed in these negative-norm spaces.


title: The HI supershell GS 118+01-44 and its role in the interstellar medium

abstract: We carry out a multiwavelength study to characterize the HI supershell designated GS 118+01-44, and to analyse its possible origin. A multiwavlength study has been carried out to study the supershell and its environs. We performed an analysis of the HI, CO, radio continuum, and infrared emission distributions. The Canadian Galactic Plane Survey (CGPS) HI data reveals that GS 118+01-44 is centred at (l, b) = (117.7, 1.4) with a systemic velocity of -44.3 km/s. According to Galactic rotation models this structure is located at 3.0 +- 0.6 kpc from the Sun. There are several HII regions and three supernova remnants (SNRs) catalogued in the region. On the other hand, the analysis of the temperature spectral index distribution shows that in the region there is a predominance of non-thermal emission. Infrared emission shows that cool temperatures dominate the area of the supershell. Concerning the origin of the structure, we found that even though several OB stars belonging to Cas OB5 are located in the interior of GS 118+01-44, an analysis of the energy injected by these stars through their stellar winds indicates that they do not have sufficient energy to create GS 118+01-44. Therefore, an additional energy source is needed to explain the genesis of GS 118+01-44. On the other hand, the presence of several HII regions and young stellar object candidates in the edges of GS 118+01-44 shows that the region is still active in forming new stars.


title: Farrell-Jones Conjecture for fundamental groups of graphs of virtually cyclic groups

abstract: In this note, we prove the K- and L-theoretic Farrell-Jones Conjecture with coefficients in an additive category for fundamental groups of graphs of virtually cyclic groups.


title: Learning Optimized Or's of And's

abstract: Or's of And's (OA) models are comprised of a small number of disjunctions of conjunctions, also called disjunctive normal form. An example of an OA model is as follows: If ($x_1 = $ `blue' AND $x_2=$ `middle') OR ($x_1 = $ `yellow'), then predict $Y=1$, else predict $Y=0$. Or's of And's models have the advantage of being interpretable to human experts, since they are a set of conditions that concisely capture the characteristics of a specific subset of data. We present two optimization-based machine learning frameworks for constructing OA models, Optimized OA (OOA) and its faster version, Optimized OA with Approximations (OOAx). We prove theoretical bounds on the properties of patterns in an OA model. We build OA models as a diagnostic screening tool for obstructive sleep apnea, that achieves high accuracy with a substantial gain in interpretability over other methods.


title: A sharp lower bound for choosing the maximum of an independent sequence

abstract: This paper considers a variation of the full-information secretary problem where the random variables to be observed are independent but not necessary identically distributed. The main result is a sharp lower bound for the optimal win probability. Precisely, if $X_1,\dots,X_n$ are independent random variables with known continuous distributions and $V_n(X_1,\dots,X_n):=\sup_\tau P(X_\tau=M_n)$, where $M_n:=\max\{X_1,\dots,X_n\}$ and the supremum is over all stopping times adapted to $X_1,\dots,X_n$, then $$V_n(X_1,\dots,X_n)\geq \left(1-\frac{1}{n}\right)^{n-1},$$ and this bound is attained. The method of proof consists in reducing the problem to that of a sequence of two-valued random variables, and then applying Bruss' sum-the-odds theorem (2000). In order to obtain a sharp bound for each $n$, we improve Bruss' lower bound (2003) for the sum-the-odds problem.


title: How Large is $A_g(\mathbb{F}_q)$?

abstract: Let $B(g,p)$ denote the number of isomorphism classes of $g$-dimensional abelian varieties over the finite field of size $p.$ Let $A(g,p)$ denote the number of isomorphism classes of principally polarized $g$ dimensional abelian varieties over the finite field of size $p.$ We derive upper bounds for $B(g,p)$ and lower bounds for $A(g,p)$ for $p$ fixed and $g$ increasing. The extremely large gap between the lower bound for $A(g,p)$ and the upper bound $B(g,p)$ implies some statistically counterintuitive behavior for abelian varieties of large dimension over a fixed finite field.


title: Fate of the Wiedemann-Franz law near quantum critical points of electron systems in solids

abstract: We introduce and analyze two different scenarios for violation of the Wiedemann-Franz law in strongly correlated electron systems of solids, close to a topological quantum critical point (TQCP) where the density of states $N(0)$ diverges. The first, applicable to the Fermi-liquid (FL) side of the TQCP, involves a transverse zero-sound collective mode that opens a new channel for the thermal conductivity, thereby enhancing the Lorenz number $L(0)$ relative to the value $L_0=\pi^2k^2_B/3e^2$ dictated by conventional FL theory. The second mechanism for violation of the WF law, relevant to the non-Fermi-liquid (NFL) side of the TQCP, involves the formation of a flat band and leads instead to a reduction of the Lorenz number.


title: The linear response function of an idealized atmosphere. Part 1: Construction using Green's functions and applications

abstract: A linear response function (LRF) determines the mean-response of a nonlinear climate system to weak imposed forcings, and an eddy flux matrix (EFM) determines the eddy momentum and heat flux responses to mean-flow changes. Neither LRF nor EFM can be calculated from first principles due the lack of a complete theory for turbulent eddies. Here the LRF and EFM for an idealized dry atmosphere are computed by applying numerous localized weak forcings, one at a time, to a GCM with Held-Suarez physics and calculating the mean-responses. The LRF and EFM for zonally-averaged responses are then constructed using these forcings and responses through matrix inversion. Tests demonstrate that LRF and EFM are fairly accurate. Spectral analysis of the LRF shows that the most excitable dynamical mode, the neutral vector, strongly resembles the model's Annular Mode. The framework described here can be employed to compute the LRF/EFM for zonally-asymmetric responses and more complex GCMs. The potential applications of the LRF/EFM constructed here are i) forcing a specified mean-flow for hypothesis-testing, ii) isolating/quantifying the eddy-feedbacks in complex eddy-mean flow interaction problems, and iii) evaluating/improving more generally-applicable methods currently used to construct LRFs or diagnose eddy-feedbacks in comprehensive GCMs or observations. As an example for iii, in Part 2, the LRF is also computed using the fluctuation-dissipation theorem (FDT), and the previously-calculated LRF is exploited to investigate why FDT performs poorly in some cases. It is shown that dimension-reduction using leading EOFs, which is commonly used to construct LRFs from the FDT, can significantly degrade the accuracy due to the non-normality of the operator.


title: Operator analysis of effective spin-flavor interactions for L=1 excited baryons

abstract: We match the non-relativistic quark model, with both flavor dependent and flavor independent effective quark-quark interactions, to the spin-flavor operator basis of the 1/Nc expansion for the L=1 non-strange baryons. We obtain analytic expressions for the coefficients of the 1/Nc operators in terms of radial integrals that depend on the shape and relative strength of the spin-spin, spin-orbit and tensor interactions of the model, which are left unspecified. We obtain several new, parameter-free relations between the seven masses and the two mixing angles that can discriminate between different spin-flavor structures of the effective quark-quark interaction. We discuss in detail how a general parametrization of the mass matrix depends on the mixing angles and is constrained by the assumptions on the effective quark-quark interaction. We find that, within the present experimental uncertainties, consistency with the best values of the mixing angles as determined by a recent global fit to masses and decays does not exclude any of the two most extreme possibilities of flavor dependent (independent) quark-quark interactions, as generated by meson (gluon) exchange interactions.


title: Atomic Scale Visualization of Quantum Interference on a Weyl Semimetal Surface by Scanning Tunneling Microscopy/Spectroscopy

abstract: Weyl semimetals may open a new era in condensed matter physics, materials science and nanotech after graphene and topological insulators. We report the first atomic scale view of the surface states of a Weyl semimetal (NbP) using scanning tunneling microscopy/spectroscopy. We observe coherent quantum interference patterns that arise from the scattering of quasiparticles near point defects on the surface. The measurements reveal the surface electronic structure both below and above the chemical potential in both real and reciprocal spaces. Moreover, the interference maps uncover the scattering processes of NbP's exotic surface states. Through comparison between experimental data and theoretical calculations, we further discover that the scattering channels are largely restricted by the orbital and/or spin texture of the surface band. The visualization of the scattering processes can help design novel transport effects and electronics on the topological surface of a Weyl semimetal.


title: Exact Bremsstrahlung and Effective Couplings

abstract: We calculate supersymmetric Wilson loops on the ellipsoid for a large class of $\mathcal{N}=2$ SCFT using the localization formula of Hama and Hosomichi. From them we extract the radiation emitted by an accelerating heavy probe quark as well as the entanglement entropy following the recent works of Lewkowycz-Maldacena and Fiol-Gerchkovitz-Komargodski. Comparing our results with the $\mathcal{N}=4$ SYM ones, we obtain interpolating functions $f(g^2)$ such that a given $\mathcal{N}=2$ SCFT observable is obtained by replacing in the corresponding $\mathcal{N}=4$ SYM result the coupling constant by $f(g^2)$. These "exact effective couplings" encode the finite, relative renormalization between the $\mathcal{N}=2$ and the $\mathcal{N}=4$ gluon propagator, they interpolate between the weak and the strong coupling. We discuss the range of their applicability.


title: Comparison of nanomechanical properties of in vivo and in vitro keratin 8/18 networks

abstract: In our work we compare the mechanical properties of the extracted keratin cytoskeleton of pancreatic carcinoma cells with the mechanical properties of in vitro assembled keratin 8/18 networks. For this purpose we use microrheology measurements with embedded tracer beads. This method is a suitable tool, because the size of the beads compared to the meshsize of the network allows us to treat the network as a continuum. Observing the beads motion with a CCD-High-Speed-Camera then leads to the dynamic shear modulus. Our measurements show lower storage moduli with increasing distance between the rim of the nucleus and the bead, but no clear tendency for the loss modulus. The same measurement method applied to in vitro assembled keratin 8/18 networks shows different characteristics of storage and loss moduli. The storage modulus is one order of magnitude lower than that of the extracted cytoskeleton and the loss modulus is higher. We draw conclusions on the network topology of both keratin network types based on the mechanical behaviour.


title: Self-organized magnetic particles to tune the mechanical behaviour of a granular system

abstract: Above a certain density a granular material jams. This property can be controlled by either tuning a global property, such as the packing fraction or by applying shear strain, or at the micro-scale by tuning grain shape, inter-particle friction or externally controlled organization. Here, we introduce a novel way to change a local granular property by adding a weak anisotropic magnetic interaction between particles. We measure the evolution of the pressure, $P$, and coordination number, $Z$, for a packing of 2D photo-elastic disks, subject to uniaxial compression. Some of the particles have embedded cuboidal magnets. The strength of the magnetic interactions between particles are too weak to have a strong direct effect on $P$ or $Z$ when the system is jammed. However, the magnetic interactions play an important role in the evolution of latent force networks when systems containing a large enough fraction of the particles with magnets are driven through unjammed states. In this case, a statistically stable network of magnetic chains self-organizes and overlaps with force chains, strengthening the granular medium. We believe this property can be used to reversibly control mechanical properties of granular materials.


title: Second-order cosmological perturbation theory and initial conditions for $N$-body simulations

abstract: We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.


title: On correlations of certain multiplicative functions

abstract: In this paper, we study sums of shifted products $\sum\limits_{n \leq x} F(n) G(n-h)$ for any $|h| \leq x/2$ and arithmetic functions $F=f*1$ and $G=g*1$, with $f$ and $g$ small. We obtain asymptotic formula for different orders of magnitude of $f$ and $g$. We also provide asymptotic formula for sums of the type $\sum\limits_{n \leq x} \mu^2(n) G(n-h)$, where $G=g*1$ and $g$ is small. For small order of magnitudes of $f$ and $g$, we improve the error terms and make them independent of $h$.


title: Deep Kernel Learning

abstract: We introduce scalable deep kernels, which combine the structural properties of deep learning architectures with the non-parametric flexibility of kernel methods. Specifically, we transform the inputs of a spectral mixture base kernel with a deep architecture, using local kernel interpolation, inducing points, and structure exploiting (Kronecker and Toeplitz) algebra for a scalable kernel representation. These closed-form kernels can be used as drop-in replacements for standard kernels, with benefits in expressive power and scalability. We jointly learn the properties of these kernels through the marginal likelihood of a Gaussian process. Inference and learning cost $O(n)$ for $n$ training points, and predictions cost $O(1)$ per test point. On a large and diverse collection of applications, including a dataset with 2 million examples, we show improved performance over scalable Gaussian processes with flexible kernel learning models, and stand-alone deep architectures.


title: Non-classicality indicators for entangled and squeezed number states

abstract: Many non-classicality indicators are used to measure quantum effects of different systems. Kenfack's and Sadeghi's non-classicality indicators are introduced in terms of the amount of Wigner function's negativities and interferences in phase space quantum mechanics, respectively. They are, effectively, applied for some real distribution functions. In this paper, we compared these non-classicality indicators for the entangled state of photonic number states in the Wigner, Husimi and Rivier representations. It is shown that for a two-level entangled state, Sadeghi's indicator has more benefits with respect to the Kenfack's indicator. For the two-level entangled state, we show a correspondence between the Sadeghi's non-classicality indicator and the Von Neumann entropy. It is also shown that for the superposition of squeezed number state the Sadeghi's (and no Kenfack's) non-classicality indicators is sensible with respect to the squeezing parameter for a superposition of squeezed number states.


title: Search for a dark photon at BABAR

abstract: We present the results of a search for a dark photon decaying to $e^+e^-$ or $\mu^+\mu^-$ in the \babar\ detector. We find no evidence for such a dark photon and set upper limits on the dark mixing as a function mass, tightening the constraints on the parameter space of several dark sector theories.


title: Doublon dynamics and polar molecule production in an optical lattice

abstract: Ultracold polar molecules provide an excellent platform to study quantum many-body spin dynamics, which has become accessible in the recently realized low entropy quantum gas of polar molecules in an optical lattice. To obtain a detailed understanding for the molecular formation process in the lattice, we prepare a density distribution where lattice sites are either empty or occupied by a doublon composed of a bosonic atom interacting with a fermionic atom. By letting this disordered, out-of-equilibrium system evolve from a well-defined initial condition, we observe clear effects on pairing that arise from inter-species interactions, a higher partial wave Feshbach resonance, and excited Bloch-band population. When only the lighter fermions are allowed to tunnel in the three-dimensional (3D) lattice, the system dynamics can be well described by theory. However, in a regime where both fermions and bosons can tunnel, we encounter correlated dynamics that is beyond the current capability of numerical simulations. Furthermore, we show that we can probe the microscopic distribution of the atomic gases in the lattice by measuring the inelastic loss of doublons. These techniques realize tools that are generically applicable to heteronuclear diatomic systems in optical lattices and can shed light on molecule production as well as dynamics of a Bose-Fermi mixture.


title: Classification and detection of symmetry fractionalization in chiral spin liquids

abstract: In this work we study the crystal symmetry fractionalization in chiral spin liquids with the chiral-semion topological order. We show that if such a chiral spin liquid is realized in a two-dimensional lattice model with odd number of spin-$\frac12$ per unit cell and the state preserves spin rotation symmetry and translation symmetries, the semion excitation must carry both half-integer spin and fractional crystal quantum numbers. As a result, only a unique symmetry enriched topological phase can be realized in chiral spin liquids in a spin-$\frac12$ kagome lattice model. These fractional symmetry quantum numbers are confirmed numerically using ground state wave functions obtained with the density matrix renormalization group method.


title: Semiclassical electron transport at the edge of a 2D topological insulator: Interplay of protected and unprotected modes

abstract: We study electron transport at the edge of a generic disordered two-dimensional topological insulator, where some channels are topologically protected from backscattering. Assuming the total number of channels is large, we consider the edge as a quasi-one-dimensional quantum wire and describe it in terms of a non-linear sigma model with a topological term. Neglecting localization effects, we calculate the average distribution function of transmission probabilities as a function of the sample length. We mainly focus on the two experimentally relevant cases: a junction between two quantum Hall (QH) states with different filling factors (unitary class) and a relatively thick quantum well exhibiting quantum spin Hall (QSH) effect (symplectic class). In a QH sample, the presence of topologically protected modes leads to a strong suppression of diffusion in the other channels already at the scales much shorter than the localization length. On the semiclassical level, this is accompanied by the formation of a gap in the spectrum of transmission probabilities close to unit transmission, thereby suppressing shot noise and conductance fluctuations. In the case of a QSH system, there is at most one topologically protected edge channel leading to weaker transport effects. In order to describe `topological' suppression of nearly perfect transparencies, we develop an exact mapping of the semiclassical limit of the one-dimensional sigma model onto a zero-dimensional sigma model of a different symmetry class, allowing us to identify the distribution of transmission probabilities with the average spectral density of a certain random-matrix ensemble. We extend our results to other symmetry classes with topologically protected edges in two dimensions.


title: Seven ways to improve example-based single image super resolution

abstract: In this paper we present seven techniques that everybody should know to improve example-based single image super resolution (SR): 1) augmentation of data, 2) use of large dictionaries with efficient search structures, 3) cascading, 4) image self-similarities, 5) back projection refinement, 6) enhanced prediction by consistency check, and 7) context reasoning. We validate our seven techniques on standard SR benchmarks (i.e. Set5, Set14, B100) and methods (i.e. A+, SRCNN, ANR, Zeyde, Yang) and achieve substantial improvements.The techniques are widely applicable and require no changes or only minor adjustments of the SR methods. Moreover, our Improved A+ (IA) method sets new state-of-the-art results outperforming A+ by up to 0.9dB on average PSNR whilst maintaining a low time complexity.


title: Wage gap between men and women in Tunisia

abstract: This paper focuses on estimating wage differences between males and females in Tunisia by using the Oaxaca-Blinder decomposition, a technical that isolates wage gap due to characteristics, from wage gap due to discrimination against women. The data used in the analysis is obtained from the Tunisian Population and Employment Survey 2005. It is estimated that, the gender wage gap is about 19% and the results ascertain that the gender wage gap is mostly attributed to discrimination, especially to underestimation of females'caracteristics on the labor market.


title: Doppler-imaging of the planetary debris disc at the white dwarf SDSS J122859.93+104032.9

abstract: Debris discs which orbit white dwarfs are signatures of remnant planetary systems. We present twelve years of optical spectroscopy of the metal-polluted white dwarf SDSS J1228+1040, which shows a steady variation in the morphology of the 8600 {\AA} Ca II triplet line profiles from the gaseous component of its debris disc. We identify additional emission lines of O I, Mg I, Mg II, Fe II and Ca II in the deep co-added spectra. These emission features (including Ca H & K) exhibit a wide range in strength and morphology with respect to each other and to the Ca II triplet, indicating different intensity distributions of these ionic species within the disc. Using Doppler tomography we show that the evolution of the Ca II triplet profile can be interpreted as the precession of a fixed emission pattern with a period in the range 24-30 years. The Ca II line profiles vary on time-scales that are broadly consistent with general relativistic precession of the debris disc.


title: Towards a self-consistent halo model for the nonlinear large-scale structure

abstract: The halo model is a theoretically and empirically well-motivated framework for predicting the statistics of the nonlinear matter distribution in the Universe. However, current incarnations of the halo model suffer from two major deficiencies: $(i)$ they do not enforce the stress-energy conservation of matter; $(ii)$ they are not guaranteed to recover exact perturbation theory results on large scales. Here, we provide a formulation of the halo model ("EHM") that remedies both drawbacks in a consistent way, while attempting to maintain the predictivity of the approach. In the formulation presented here, mass and momentum conservation are guaranteed on large scales, and results of perturbation theory and the effective field theory can in principle be matched to any desired order on large scales. We find that a key ingredient in the halo model power spectrum is the halo stochasticity covariance, which has been studied to a much lesser extent than other ingredients such as mass function, bias, and profiles of halos. As written here, this approach still does not describe the transition regime between perturbation theory and halo scales realistically, which is left as an open problem. We also show explicitly that, when implemented consistently, halo model predictions do not depend on any properties of low-mass halos that are smaller than the scales of interest.


title: Flaring of tidally compressed dark-matter clumps

abstract: We explore the physics and observational consequences of tidal compression events (TCEs) of dark-matter clumps (DMCs) by supermassive black holes (SMBHs). Our analytic calculations show that a DMC approaching a SMBH much closer than the tidal radius undergoes significant compression along the axis perpendicular to the orbital plane, shortly after pericenter passage. For DMCs composed of self-annihilating dark-matter particles, we find that the boosted DMC density and velocity dispersion lead to a flaring of the annihilation rate, most pronounced for a velocity- dependent annihilation cross section. If the end products of the annihilation are photons, this results in a gamma-ray flare, detectable (and possibly already detected) by the Fermi telescope for a range of model parameters. If the end products of dark-matter annihilation are relativistic electrons and positrons and the local magnetic field is large enough, TCEs of DMCs can lead to flares of synchrotron radiation. Finally, TCEs of DMCs lead to a burst of gravitational waves, in addition to the ones radiated by the orbital motion alone, and with a different frequency spectrum. These transient phenomena provide interesting new avenues to explore the properties of dark matter.


title: The Spatial Distribution of the Young Stellar Clusters in the Star Forming Galaxy NGC 628

abstract: We present a study of the spatial distribution of the stellar cluster populations in the star forming galaxy NGC 628. Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey), we have identified 1392 potential young (<100 Myr) stellar clusters within the galaxy, identified from a combination of visual inspection and automatic selection. We investigate the clustering of these young stellar clusters and quantify the strength and change of clustering strength with scale using the two-point correlation function. We also investigate how image boundary conditions and dust lanes affect the observed clustering. The distribution of the clusters is well fit by a broken power law with negative exponent $\alpha$. We recover a weighted mean index of $\alpha$ ~ -0.8 for all spatial scales below the break at 3".3 (158 pc at a distance of 9.9 Mpc) and an index of $\alpha$ ~ -0.18 above 158 pc for the accumulation of all cluster types. The strength of the clustering increases with decreasing age and clusters older than 40 Myr lose their clustered structure very rapidly and tend to be randomly distributed in this galaxy whereas the mass of the star cluster has little effect on the clustering strength. This is consistent with results from other studies that the morphological hierarchy in stellar clustering resembles the same hierarchy as the turbulent interstellar medium.


title: The Low-Mass Astrometric Binary LSR1610-0040

abstract: Even though it was discovered more than a decade ago, LSR1610-0040 remains an enigma. This object has a peculiar spectrum that exhibits some features typically found in L subdwarfs, and others common in the spectra of more massive M dwarf stars. It is also a binary system with a known astrometric orbital solution. Given the available data, it remains a challenge to reconcile the observed properties of the combined light of LSR1610-0040AB with current theoretical models of low-mass stars and brown dwarfs. We present the results of a joint fit to both astrometric and radial velocity measurements of this unresolved, low-mass binary. We find that the photocentric orbit has a period $P = 633.0 \pm 1.7$ days, somewhat longer than previous results, with eccentricity of $e = 0.42 \pm 0.03$, and we estimate that the semi-major axis of the orbit of the primary is $a_1 \approx 0.32$ AU, consistent with previous results. While a complete characterization of the system is limited by our small number of radial velocity measurements, we establish a likely primary mass range of $0.09 - 0.10 M_\odot$ from photometric and color-magnitude data. For a primary mass in this range, the secondary is constrained to be $0.06 - 0.075 M_\odot$, making a negligible contribution to the total I-band luminosity. This effectively rules out the possibility of the secondary being a compact object such as an old, low-mass white dwarf. Based on our analysis, we predict a likely angular separation at apoapsis comparable to the resolution limits of current high-resolution imaging systems. Measuring the angular separation of the A & B components would finally enable a full, unambiguous solution for the masses of the components of this system.


title: NAND-Trees, Average Choice Complexity, and Effective Resistance

abstract: We show that the quantum query complexity of evaluating NAND-tree instances with average choice complexity at most $W$ is $O(W)$, where average choice complexity is a measure of the difficulty of winning the associated two-player game. This generalizes a superpolynomial speedup over classical query complexity due to Zhan et al. [Zhan et al., ITCS 2012, 249-265]. We further show that the player with a winning strategy for the two-player game associated with the NAND-tree can win the game with an expected $\widetilde{O}(N^{1/4}\sqrt{{\cal C}(x)})$ quantum queries against a random opponent, where ${\cal C }(x)$ is the average choice complexity of the instance. This gives an improvement over the query complexity of the naive strategy, which costs $\widetilde{O}(\sqrt{N})$ queries. The results rely on a connection between NAND-tree evaluation and $st$-connectivity problems on certain graphs, and span programs for $st$-connectivity problems. Our results follow from relating average choice complexity to the effective resistance of these graphs, which itself corresponds to the span program witness size.


title: Velocity anisotropy in tidally limited star clusters

abstract: We explore the long-term evolution of the anisotropy in the velocity space of star clusters starting with different structural and kinematical properties. We show that the evolution of the radial anisotropy strength and its radial variation within a cluster contain distinct imprints of the cluster initial structural properties, dynamical history, and of the external tidal field of its host galaxy. Initially isotropic and compact clusters with small initial values of the ratio of the half-mass to Jacobi radius, $r_h/r_J$, develop a strong radial anisotropy during their long-term dynamical evolution. Many clusters, if formed with small values of $r_h/r_J$, should now be characterized by a significant radial anisotropy increasing with the distance from the cluster centre, reaching its maximum at a distance between 0.2 $r_J$ and 0.4 $r_J$, and then becoming more isotropic or mildly tangentially anisotropic in the outermost regions. A similar radial variation of the anisotropy can also result from an early violent relaxation phase. In both cases, as a cluster continues its evolution and loses mass, the anisotropy eventually starts to decrease and the system evolves toward an isotropic velocity distribution. However, in order to completely erase the strong anisotropy developed by these compact systems during their evolution, they must be in the advanced stages of their evolution and lose a large fraction of their initial mass. Clusters that are initially isotropic and characterized by larger initial values of $r_h/r_J$, on the other hand, never develop a significant radial anisotropy.


title: Instanton-dyon Ensemble with two Dynamical Quarks: the Chiral Symmetry Breaking

abstract: This is the second paper of the series aimed at understanding the ensemble of instanton-dyons, now with two flavors of light dynamical quarks. The partition function is appended by the fermionic factor, $(det T)^{N_f}$ and Dirac eigenvalue spectra at small values are derived from the numerical simulation of 64 and 128 dyons. Those spectra show clear chiral symmetry breaking pattern at high dyon density.


title: Twisting algebraically special solutions in five dimensions

abstract: We determine the general form of the solutions of the five-dimensional vacuum Einstein equations with cosmological constant for which (i) the Weyl tensor is everywhere type II or more special in the null alignment classification of Coley et al., and (ii) the $3 \times 3$ matrix encoding the expansion, shear and twist of the aligned null direction has rank 2. The dependence of the solution on 2 coordinates is determined explicitly, so the Einstein equation reduces to PDEs in the 3 remaining coordinates, just as for four-dimensional algebraically special solutions. The solutions fall into several families. One of these consists of warped products of four-dimensional algebraically special solutions. The others are new.


title: Zoomed high-resolution simulations of Multi-coupled Dark Energy: cored galaxy density profiles at high redshift

abstract: We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct dark matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard Lambda-CDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ~ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.


title: Climbing Mont Blanc - A Training Site for Energy Efficient Programming on Heterogeneous Multicore Processors

abstract: Climbing Mont Blanc (CMB) is an open online judge used for training in energy efficient programming of state-of-the-art heterogeneous multicores. It uses an Odroid-XU3 board from Hardkernel with an Exynos Octa processor and integrated power sensors. This processor is three-way heterogeneous containing 14 different cores of three different types. The board currently accepts C and C++ programs, with support for OpenCL v1.1, OpenMP 4.0 and Pthreads. Programs submitted using the graphical user interface are evaluated with respect to time, energy used, and energy-efficiency (EDP). A small and varied set of problems are available, and the system is currently in use in a medium sized course on parallel computing at NTNU. Other online programming judges exist, but we are not aware of any similar system that also reports energy-efficiency.


title: General Relativity in Post Independence India

abstract: The most outstanding contribution to general relativity in this era came in 1953 (published in 1955 \cite{akr}) in the form of the Raychaudhri equation. It is in 1960s that the observations began to confront the eupherial theory and thus began exploration of GR as a legitimate physical theory in right earnest. The remarkable discoveries of cosmic microwave background radiation, quasars, rotating Kerr black hole and the powerful singularity theorems heralded a new canvas of relativistic astrophysics and cosmology. I would attempt to give a brief account of Indian participation in these exciting times.

categories: physics.hist-ph astro-ph.HE gr-qc hep-th


title: New asymptotic heat transfer model in thin liquid films

abstract: In this article, we present a model of heat transfer occurring through a li\-quid film flowing down a vertical wall. This new model is formally derived using the method of asymptotic expansions by introducing appropriately chosen dimensionless variables. In our study the small parameter, known as the film parameter, is chosen as the ratio of the flow depth to the characteristic wavelength. A new Nusselt solution should be explained, taking into account the hydrodynamic free surface variations and the contributions of the higher order terms coming from temperature variation effects. Comparisons are made with numerical solutions of the full Fourier equations in a steady state frame. The flow and heat transfer are coupled through Marangoni and temperature dependent viscosity effects. Even if these effects have been considered separately before, here a fully coupled model is proposed. Another novelty consists in the asymptotic approach in contrast to the weighted residual approach which have been formerly applied to these problems.

categories: physics.flu-dyn math.NA physics.class-ph physics.comp-ph


title: Nonintegrability of NS5-like Interface in $\mathcal{N}=4$ Supersymmetric Yang-Mills

abstract: Four-dimensional $\mathcal{N}=4$ super Yang-Mills theory admits interfaces that preserve integrability of the theory. It was shown that addition of fundamental hypermultiplets living on a co-dimension one defect is an example of such an interface. We consider NS5 interface as a different example of half-BPS defect in this theory and prove that integrability is generically broken already at one loop. We show that one-loop dilatation operator acting on boundary local operators does not lead to an integrable spin chain. The same is true also for D5-like and NS5-like boundaries.


title: Understanding diboson anomalies

abstract: We conduct a model-independent effective theory analysis of hypercharged fields with various spin structures towards understanding the diboson excess found in LHC run I, as well as possible future anomalies involving $WZ$ and $WH$ modes. Within the assumption of no additional physics beyond the standard model up to the scale of the possible diboson resonance, we show that a hypercharged scalar and a spin 2 particle do not have tree-level $WZ$ and $WH$ decay channels up to dimension 5 operators, and cannot therefore account for the anomaly, whereas a hypercharged vector is a viable candidate provided we also introduce a $Z'$ in order to satisfy electroweak precision constraints. We calculate bounds on the $Z'$ mass consistent with the Atlas/CMS diboson signals as well as electroweak precision data, taking into account both LHC run I and II data.


title: Galaxy And Mass Assembly (GAMA): growing up in a bad neighbourhood - how do low-mass galaxies become passive?

abstract: Both theoretical predictions and observations of the very nearby Universe suggest that low-mass galaxies (log$_{10}$[M$_{*}$/M$_{\odot}$]<9.5) are likely to remain star-forming unless they are affected by their local environment. To test this premise, we compare and contrast the local environment of both passive and star-forming galaxies as a function of stellar mass, using the Galaxy and Mass Assembly survey. We find that passive fractions are higher in both interacting pair and group galaxies than the field at all stellar masses, and that this effect is most apparent in the lowest mass galaxies. We also find that essentially all passive log$_{10}$[M$_{*}$/M$_{\odot}$]<8.5 galaxies are found in pair/group environments, suggesting that local interactions with a more massive neighbour cause them to cease forming new stars. We find that the effects of immediate environment (local galaxy-galaxy interactions) in forming passive systems increases with decreasing stellar mass, and highlight that this is potentially due to increasing interaction timescales giving sufficient time for the galaxy to become passive via starvation. We then present a simplistic model to test this premise, and show that given our speculative assumptions, it is consistent with our observed results.


title: Lectures on the classical moment problem and its noncommutative generalization

abstract: These notes contain a presentation of the noncommutative generalization of the classical moment problem introduced in [10] and [12]. They also contain a short summary of the classical moment problem in infinite dimension.


title: Optimal control of complex atomic quantum systems

abstract: Quantum technologies will ultimately require manipulating many-body quantum systems with high precision. Cold atom experiments represent a stepping stone in that direction: a high degree of control has been achieved on systems of increasing complexity, however, this control is still sub-optimal. Optimal control theory is the ideal candidate to bridge the gap between early stage and optimal experimental protocols, particularly since it was extended to encompass many-body quantum dynamics. Here, we experimentally demonstrate optimal control applied to two dynamical processes subject to interactions: the coherent manipulation of motional states of an atomic Bose-Einstein condensate and the crossing of a quantum phase transition in small systems of cold atoms in optical lattices. We show theoretically that these transformations can be made fast and robust with respect to perturbations, including temperature and atom number fluctuations, resulting in a good agreement between theoretical predictions and experimental results.


title: Self-replication with magnetic dipolar colloids

abstract: Colloidal self-replication represents an exciting research frontier in soft matter physics. Currently, all reported self-replication schemes involve coating colloidal particles with stimuli-responsive molecules to allow switchable interactions. In this paper, we introduce a scheme using ferromagnetic dipolar colloids and pre-programmed external magnetic fields to create an autonomous self-replication system. Interparticle dipole-dipole forces and periodically varying weak-strong magnetic fields cooperate to drive colloid monomers from the solute onto templates, bind them into replicas, and dissolve template complexes. We present three general design principles for autonomous linear replicators, derived from a focused study of a minimalist sphere-dimer magnetic system in which single binding sites allow formation of dimeric templates. We show via statistical models and computer simulations that our system exhibits nonlinear growth of templates and produces nearly exponential growth (low error rate) upon adding an optimized competing electrostatic potential. We devise experimental strategies for constructing the required magnetic colloids based on documented laboratory techniques. We also present qualitative ideas about building more complex self-replicating structures utilizing magnetic colloids.


title: Tricomplex dynamical systems generated by polynomials of odd degree

abstract: In this article, we give the exact interval of the cross section of the Multibrot sets generated by the polynomial $z^p+c$ where $z$ and $c$ are complex numbers and $p > 2$ is an odd integer. Furthermore, we show that the same Multibrots defined on the hyperbolic numbers are always squares. Moreover, we give a generalized 3D version of the hyperbolic Multibrot set and prove that our generalization is an octahedron for a specific 3D slice of the dynamical system generated by the tricomplex polynomial $\eta^p+c$ where $p > 2$ is an odd integer.


title: Applications of Physics to Archery

abstract: Archery lends itself to scientific analysis. In this paper we discuss physics laws that relate to the mechanics of bow and arrow, to the shooting process and to the flight of the arrow. In parallel, we describe experiments that address these laws. The detailed results of these measurements, performed with a specific bow and arrow, provide insight into many aspects of archery and illustrate the importance of quantitative information in the scientific process. Most of the proposed experiments use only modest tools and can be carried out by archers with their own equipment.


title: Learning Visual Features from Large Weakly Supervised Data

abstract: Convolutional networks trained on large supervised dataset produce visual features which form the basis for the state-of-the-art in many computer-vision problems. Further improvements of these visual features will likely require even larger manually labeled data sets, which severely limits the pace at which progress can be made. In this paper, we explore the potential of leveraging massive, weakly-labeled image collections for learning good visual features. We train convolutional networks on a dataset of 100 million Flickr photos and captions, and show that these networks produce features that perform well in a range of vision problems. We also show that the networks appropriately capture word similarity, and learn correspondences between different languages.


title: The weirdest martensite: Smectic liquid crystal microstructure and Weyl-Poincar\'e invariance

abstract: Smectic liquid crystals are remarkable, beautiful examples of materials microstructure, with ordered patterns of geometrically perfect ellipses and hyperbolas. The solution of the complex problem of filling three-dimensional space with domains of focal conics under constraining boundary conditions yields a set of strict rules, which are similar to the compatibility conditions in a martensitic crystal. Here we present the rules giving compatible conditions for the concentric circle domains found at two-dimensional smectic interfaces with planar boundary conditions. Using configurations generated by numerical simulations, we develop a clustering algorithm to decompose the planar boundaries into domains. The interfaces between different domains agree well with the smectic compatibility conditions. We also discuss generalizations of our approach to describe the full three-dimensional smectic domains, where the variant symmetry group is the Weyl-Poincar\'e group of Lorentz boosts, translations, rotations, and dilatations.


title: The Variationally Fitted Electron-Electron Potential

abstract: Perhaps the simplest first-principles approach to electronic structure is to fit the charge distribution of each orbital pair and use those fits wherever they appear in the entire electron-electron (EE) interaction energy. The charge distributions in quantum chemistry are typically represented as a sums over products of Gaussian orbital basis functions. If fitted, they are also represented as a sum over single-center Gaussian fitting basis functions. With two representations of the charge distributions, the proper definition of energy is ambiguous. To remedy this, we require that the variation of the energy with respect to a product of orbitals generates a fitted potential. This makes the quantum-mechanical energy robust, i.e. corrected to first order for the error made using an incomplete fitting basis. The coupled orbital and fitting equations are then the result of making the energy stationary with respect to two independent sets of variables. We define the potentials and unique energies for methods based on the Hartree Fock model and variationally fit the full EE interaction in DFT. We compare implementations of variational fitting in DFT at six different levels for three different functionals. Our calculations are performed on transition metal atoms, for which first-order Coulomb errors, due to an incomplete fitting basis sets, are significant. Variational first-order exchange and correlation errors have similar magnitude in all cases. Robust energy differences are much smaller, particularly in the local density approximation.


title: Active Perceptual Similarity Modeling with Auxiliary Information

abstract: Learning a model of perceptual similarity from a collection of objects is a fundamental task in machine learning underlying numerous applications. A common way to learn such a model is from relative comparisons in the form of triplets: responses to queries of the form "Is object a more similar to b than it is to c?". If no consideration is made in the determination of which queries to ask, existing similarity learning methods can require a prohibitively large number of responses. In this work, we consider the problem of actively learning from triplets -finding which queries are most useful for learning. Different from previous active triplet learning approaches, we incorporate auxiliary information into our similarity model and introduce an active learning scheme to find queries that are informative for quickly learning both the relevant aspects of auxiliary data and the directly-learned similarity components. Compared to prior approaches, we show that we can learn just as effectively with much fewer queries. For evaluation, we introduce a new dataset of exhaustive triplet comparisons obtained from humans and demonstrate improved performance for different types of auxiliary information.


title: Iso-oriented monolayer {\alpha}-MoO3(010) films epitaxially grown on SrTiO3(001)

abstract: The ability to synthesis well-ordered two-dimensional materials under ultra-high vacuum and directly characterize them by other techniques in-situ can greatly advance our current understanding on their physical and chemical properties. In this paper, we demonstrate that iso-oriented {\alpha}-MoO3 films with as low as single monolayer thickness can be reproducibly grown on SrTiO3(001) (STO) substrates by molecular beam epitaxy ( (010)MoO3 || (001)STO, [100]MoO3 || [100]STO or [010]STO) through a self-limiting process. While one in-plane lattice parameter of the MoO3 is very close to that of the SrTiO3 (aMoO3 = 3.96 {\AA}, aSTO = 3.905 {\AA}), the lattice mismatch along other direction is large (~5%, cMoO3 = 3.70 {\AA}), which leads to relaxation as clearly observed from the splitting of streaks in reflection high-energy electron diffraction (RHEED) patterns. A narrow range in the growth temperature is found to be optimal for the growth of monolayer {\alpha}-MoO3 films. Increasing deposition time will not lead to further increase in thickness, which is explained by a balance between deposition and thermal desorption due to the weak van der Waals force between {\alpha}-MoO3 layers. Lowering growth temperature after the initial iso-oriented {\alpha}-MoO3 monolayer leads to thicker {\alpha}-MoO3(010) films with excellent crystallinity.


title: On the Optimality of Uncoded Cache Placement

abstract: Caching is an efficient way to reduce peak-hour network traffic congestion by storing some contents at user's local cache without knowledge of later demands. Maddah-Ali and Niesen initiated a fundamental study of caching systems; they proposed a scheme (with uncoded cache placement and linear network coding delivery) that is provably optimal to within a factor 12. In this paper, by noticing that when the cache contents and the demands are fixed, the caching problem can be seen as an index coding problem, we show the optimality of Maddah-Ali and Niesen's scheme assuming that cache placement is restricted to be uncoded and the number of users is not less than the number of files. Furthermore, this result states that further improvement to the Maddah-Ali and Niesen's scheme in this regimes can be obtained only by coded cache placement.


title: Measurement of scintillation and ionization yield with high-pressure gaseous mixtures of Xe and TMA for improved neutrinoless double beta decay and dark matter searches

abstract: The gaseous Xenon(Xe) time projection chamber (TPC) is an attractive detector technique for neutrinoless double beta decay and WIMP dark matter searches. While it is less dense compared to Liquid Xe detectors, it has intrinsic advantages in tracking capability and better energy resolution. The performance of gaseous Xe can be further improved by molecular additives such as trimethylamine(TMA), which is expected to (1) cool down the ionization electrons, (2) convert Xe excitation energy to TMA ionizations through Penning transfer, and (3) produce scintillation and electroluminescence light in a more easily detectable wavelength (300 nm). In order to test the feasibility of the performance improvements with TMA, we made the first direct measurement of Penning and fluorescence transfer efficiency with gaseous mixtures of Xe and TMA. While we observed a Penning transfer efficiency up to ~35%, we found strong suppression of primary scintillation light with TMA. We also found that the primary scintillation light with Xe and TMA mixture can be well characterized by ~3% fluorescence transfer from Xe to TMA, with further suppression due to TMA self-quenching. No evidence of the scintillation light produced by recombination of TMA ions was found. This strong suppression of scintillation light makes dark matter searches quite challenging, while the possibility of improved neutrinoless double beta decay searches remains open. This work has been carried out within the context of the NEXT collaboration.


title: Efficient Multiscale Gaussian Process Regression using Hierarchical Clustering

abstract: Standard Gaussian Process (GP) regression, a powerful machine learning tool, is computationally expensive when it is applied to large datasets, and potentially inaccurate when data points are sparsely distributed in a high-dimensional feature space. To address these challenges, a new multiscale, sparsified GP algorithm is formulated, with the goal of application to large scientific computing datasets. In this approach, the data is partitioned into clusters and the cluster centers are used to define a reduced training set, resulting in an improvement over standard GPs in terms of training and evaluation costs. Further, a hierarchical technique is used to adaptively map the local covariance representation to the underlying sparsity of the feature space, leading to improved prediction accuracy when the data distribution is highly non-uniform. A theoretical investigation of the computational complexity of the algorithm is presented. The efficacy of this method is then demonstrated on smooth and discontinuous analytical functions and on data from a direct numerical simulation of turbulent combustion.


title: Deterministic Annealing Based Optimization for Zero-Delay Source-Channel Coding in Networks

abstract: This paper studies the problem of global optimization of zero-delay source-channel codes that map between the source space and the channel space, under a given transmission power constraint and for the mean square error distortion. Particularly, we focus on two well known network settings: the Wyner-Ziv setting where only a decoder has access to side information and the distributed setting where independent encoders transmit over independent channels to a central decoder. Prior work derived the necessary conditions for optimality of the encoder and decoder mappings, along with a greedy optimization algorithm that imposes these conditions iteratively, in conjunction with the heuristic noisy channel relaxation method to mitigate poor local minima. While noisy channel relaxation is arguably effective in simple settings, it fails to provide accurate global optimization in more complicated settings considered in this paper. We propose a powerful non-convex optimization method based on the concept of deterministic annealing -- which is derived from information theoretic principles and was successfully employed in several problems including vector quantization, classification and regression. We present comparative numerical results that show strict superiority of the proposed method over greedy optimization methods as well as prior approaches in literature.