title: Data and uncertainty in extreme risks - a nonlinear expectations approach

abstract: Estimation of tail quantities, such as expected shortfall or Value at Risk, is a difficult problem. We show how the theory of nonlinear expectations, in particular the Data-robust expectation introduced in [5], can assist in the quantification of statistical uncertainty for these problems. However, when we are in a heavy-tailed context (in particular when our data are described by a Pareto distribution, as is common in much of extreme value theory), the theory of [5] is insufficient, and requires an additional regularization step which we introduce. By asking whether this regularization is possible, we obtain a qualitative requirement for reliable estimation of tail quantities and risk measures, in a Pareto setting.


title: Anatomically Constrained Neural Networks (ACNN): Application to Cardiac Image Enhancement and Segmentation

abstract: Incorporation of prior knowledge about organ shape and location is key to improve performance of image analysis approaches. In particular, priors can be useful in cases where images are corrupted and contain artefacts due to limitations in image acquisition. The highly constrained nature of anatomical objects can be well captured with learning based techniques. However, in most recent and promising techniques such as CNN based segmentation it is not obvious how to incorporate such prior knowledge. State-of-the-art methods operate as pixel-wise classifiers where the training objectives do not incorporate the structure and inter-dependencies of the output. To overcome this limitation, we propose a generic training strategy that incorporates anatomical prior knowledge into CNNs through a new regularisation model, which is trained end-to-end. The new framework encourages models to follow the global anatomical properties of the underlying anatomy (e.g. shape, label structure) via learned non-linear representations of the shape. We show that the proposed approach can be easily adapted to different analysis tasks (e.g. image enhancement, segmentation) and improve the prediction accuracy of the state-of-the-art models. The applicability of our approach is shown on multi-modal cardiac datasets and public benchmarks. Additionally, we demonstrate how the learned deep models of 3D shapes can be interpreted and used as biomarkers for classification of cardiac pathologies.


title: Modeling variational inpainting methods with splines

abstract: Mathematical methods of image inpainting involve the discretization of given continuous models. We present a method that avoids the standard pointwise discretization by modeling known variational approaches, in particular total variation (TV), using a finite dimensional spline space. Besides the analysis of the resulting model, we present a numerical implementation based on the alternating method of multipliers. We compare the results numerically with classical TV inpainting and give examples of applications.


title: Learning Optimal Routing for the Uplink in LPWANs Using Similarity-enhanced epsilon-greedy

abstract: Despite being a relatively new communication technology, Low-Power Wide Area Networks (LPWANs) have shown their suitability to empower a major part of Internet of Things applications. Nonetheless, most LPWAN solutions are built on star topology (or single-hop) networks, often causing lifetime shortening in stations located far from the gateway. In this respect, recent studies show that multi-hop routing for uplink communications can reduce LPWANs' energy consumption significantly. However, it is a troublesome task to identify such energetically optimal routings through trial-and-error brute-force approaches because of time and, especially, energy consumption constraints. In this work we show the benefits of facing this exploration/exploitation problem by running centralized variations of the multi-arm bandit's epsilon-greedy, a well-known online decision-making method that combines best known action selection and knowledge expansion. Important energy savings are achieved when proper randomness parameters are set, which are often improved when conveniently applying similarity, a concept introduced in this work that allows harnessing the gathered knowledge by sporadically selecting unexplored routing combinations akin to the best known one.


title: A note on Fibonacci number of even index

abstract: We prove that a positive integer $n$ is a Fibonacci number of even index if and only if $\langle n\varphi\rangle+\frac{1}{n}>1$.


title: Compact representation of the full Broyden class of quasi-Newton updates

abstract: In this paper, we present the compact representation for matrices belonging to the the Broyden class of quasi-Newton updates, where each update may be either rank-one or rank-two. This work extends previous results solely for the restricted Broyden class of rank-two updates. In this article, it is not assumed the same Broyden update is used each iteration; rather, different members of the Broyden class may be used each iteration. Numerical experiments suggest that a practical implementation of the compact representation is able to accurately represent matrices belonging to the Broyden class of updates. Furthermore, we demonstrate how to compute the compact representation for the inverse of these matrices, as well as a practical algorithm for solving linear systems with members of the Broyden class of updates. We demonstrate through numerical experiments that the proposed linear solver is able to efficiently solve linear systems with members of the Broyden class of matrices to high accuracy. As an immediate consequence of this work, it is now possible to efficiently compute the eigenvalues of any limited-memory member of the Broyden class of matrices, allowing for the computation of condition numbers and the ability perform sensitivity analysis.


title: Tunnelling mechanism in non-commutative space with generalized uncertainty principle and Bohr-like black hole

abstract: The paper deals with non-thermal radiation spectrum by tunnelling mechanism with correction due to the generalized uncertainty principle (GUP) in the background of non-commutative geometry. Considering the reformulation of the tunnelling mechanism by Banerjee and Majhi, the Hawking radiation spectrum is evaluated through the density matrix for the outgoing modes. The GUP corrected effective temperature and the corresponding GUP corrected effective metric in non-commutative geometry are determined using Hawking's periodicity arguments. Thus, we obtain further corrections to the non-strictly thermal black hole (BH) radiation spectrum which give new final distributions. Then, we show that the GUP and the non-commutative geometry modify the Bohr-like BH recently discussed in a series of papers in the literature. In particular, we find the intriguing result that the famous law of Bekenstein on the area quantization is affected neither by non-commutative geometry nor by the GUP. This is a clear indication of the universality of Bekenstein's result. In addition, we find that both the Bekentsein-Hawking entropy and the total BH entropy to third order approximation are still functions of the BH quantum level.


title: Tropical moduli spaces of stable maps to a curve

abstract: We construct moduli spaces of rational covers of an arbitrary smooth tropical curve in R^r as tropical varieties. They are contained in the balanced fan parametrizing tropical stable maps of the appropriate degree to R^r. The weights of the top-dimensional polyhedra are given in terms of certain lattice indices and local Hurwitz numbers.


title: Radiative neutrino mass in an alternative $U(1)_{B-L}$ gauge symmetry

abstract: We propose a neutrino model in which neutrino masses are generated at one loop level and three right-handed fermions have non-trivial charges under $U(1)_{B-L}$ gauge symmetry in no conflict with anomaly cancellation. After the spontaneously symmetry breaking, a remnant $Z_2$ symmetry is induced and plays an role in assuring the stability of dark matter candidate.


title: D-vine quantile regression with discrete variables

abstract: Quantile regression, the prediction of conditional quantiles, finds applications in various fields. Often, some or all of the variables are discrete. The authors propose two new quantile regression approaches to handle such mixed discrete-continuous data. Both of them generalize the continuous D-vine quantile regression, where the dependence between the response and the covariates is modeled by a parametric D-vine. D-vine quantile regression provides very flexible models, that enable accurate and fast predictions. Moreover, it automatically takes care of major issues of classical quantile regression, such as quantile crossing and interactions between the covariates. The first approach keeps the parametric estimation of the D-vines, but modifies the formulas to account for the discreteness. The second approach estimates the D-vine using continuous convolution to make the discrete variables continuous and then estimates the D-vine nonparametrically. A simulation study is presented examining for which scenarios the discrete-continuous D-vine quantile regression can provide superior prediction abilities. Lastly, the functionality of the two introduced methods is demonstrated by a real-world example predicting the number of bike rentals.


title: Nuclear magnetic resonance line shapes of Wigner crystals in $^{13}$C-enriched graphene

abstract: Assuming that the nuclear magnetic resonance (NMR) signal from a $^{13}$C isotope enriched layer of graphene can be made sufficiently intense to be measured, we compute the NMR\ lineshape of the different crystals ground states that are expected to occur in graphene in a strong magnetic field. We first show that in nonuniform states, there is, in addition to the frequency shift due to the spin hyperfine interaction, a second contribution of equal importance from the coupling between the orbital motion of the electrons and the nuclei. We then show that, if the linewidth of the bare signal can be made sufficiently small, the Wigner and bubble crystals have line shapes that differ qualitatively from that of the uniform state at the same density while crystal states that have spin or valley pseudospin textures do not. Finally, we find that a relatively small value of the bare linewidth is sufficient to wash out the distinctive signature of the crystal states in the NMR line shape.


title: A functional perspective on emergent supersymmetry

abstract: We investigate the emergence of ${\cal N}=1$ supersymmetry in the long-range behavior of three-dimensional parity-symmetric Yukawa systems. We discuss a renormalization approach that manifestly preserves supersymmetry whenever such symmetry is realized, and use it to prove that supersymmetry-breaking operators are irrelevant, thus proving that such operators are suppressed in the infrared. All our findings are illustrated with the aid of the $\epsilon$-expansion and a functional variant of perturbation theory, but we provide numerical estimates of critical exponents that are based on the non-perturbative functional renormalization group.


title: Comparison of Gaia and asteroseismic distances

abstract: Asteroseismology provides fundamental properties (mass, radius and effective temperature) of solar-like oscillating stars using so-called scaling relations. These properties allow the computation of the asteroseismic distance of stars. We compare the asteroseismic distances with the recently released Gaia distances for 74 stars studied in Y{\i}ld{\i}z et al. There is a very good agreement between these two distances; for 64 of these stars, the difference is less than 10 per cent. However, a systematic difference is seen if we use the effective temperature obtained by spectroscopic methods; the Gaia distances are about 5 per cent greater than the asteroseismic distances.


title: Fusion of Head and Full-Body Detectors for Multi-Object Tracking

abstract: In order to track all persons in a scene, the tracking-by-detection paradigm has proven to be a very effective approach. Yet, relying solely on a single detector is also a major limitation, as useful image information might be ignored. Consequently, this work demonstrates how to fuse two detectors into a tracking system. To obtain the trajectories, we propose to formulate tracking as a weighted graph labeling problem, resulting in a binary quadratic program. As such problems are NP-hard, the solution can only be approximated. Based on the Frank-Wolfe algorithm, we present a new solver that is crucial to handle such difficult problems. Evaluation on pedestrian tracking is provided for multiple scenarios, showing superior results over single detector tracking and standard QP-solvers. Finally, our tracker ranks 2nd on the MOT16 benchmark and 1st on the new MOT17 benchmark, outperforming over 90 trackers.


title: Origin of time before inflation from a topological phase transition

abstract: We study the origin of the universe (or pre-inflation) by suggesting that the primordial space-time in the universe suffered a global topological phase transition, from a 4D Euclidean manifold to an asymptotic 4D hyperbolic one. We introduce a complex time, $\tau$, such that its real part becomes dominant after started the topological phase transition. Before the big bang, $\tau$ is a space-like coordinate, so that can be considered as a reversal variable. After the phase transition is converted in a causal variable. The formalism solves in a natural manner the quantum to classical transition of the geometrical relativistic quantum fluctuations: $\sigma$, which has a geometric origin.


title: Nordhaus-Gaddum-type theorem for conflict-free connection number of graphs

abstract: An edge-colored graph $G$ is \emph{conflict-free connected} if, between each pair of distinct vertices, there exists a path containing a color used on exactly one of its edges. The \emph{conflict-free connection number} of a connected graph $G$, denoted by $cfc(G)$, is defined as the smallest number of colors that are needed in order to make $G$ conflict-free connected. In this paper, we determine all trees $T$ of order $n$ for which $cfc(T)=n-t$, where $t\geq 1$ and $n\geq 2t+2 $. Then we prove that $1\leq cfc(G)\leq n-1$ for a connected graph $G$, and characterize the graphs $G$ with $cfc(G)=1,n-4,n-3,n-2,n-1$, respectively. Finally, we get the Nordhaus-Gaddum-type theorem for the conflict-free connection number of graphs, and prove that if $G$ and $\overline{G}$ are connected, then $4\leq cfc(G)+cfc(\overline{G})\leq n$ and $4\leq cfc(G)\cdot cfc(\overline{G})\leq2(n-2)$, and moreover, $cfc(G)+cfc(\overline{G})=n$ or $cfc(G)\cdot cfc(\overline{G})=2(n-2)$ if and only if one of $G$ and $\overline{G}$ is a tree with maximum degree $n-2$ or a $P_5$, and the lower bounds are sharp.


title: A Cloud-based Service for Real-Time Performance Evaluation of NoSQL Databases

abstract: We have created a cloud-based service that allows the end users to run tests on multiple different databases to find which databases are most suitable for their project. From our research, we could not find another application that enables the user to test several databases to gauge the difference between them. This application allows the user to choose which type of test to perform and which databases to target. The application also displays the results of different tests that were run by other users previously. There is also a map to show the location where all the tests are run to give the user an estimate of the location. Unlike the orthodox static tests and reports conducted to evaluate NoSQL databases, we have created a web application to run and analyze these tests in real time. This web application evaluates the performance of several NoSQL databases. The databases covered are MongoDB, DynamoDB, CouchDB, and Firebase. The web service is accessible from: nosqldb.nextproject.ca.


title: Identification and isotropy characterization of deformed random fields through excursion sets

abstract: A deterministic application $\theta\,:\,\mathbb{R}^2\rightarrow\mathbb{R}^2$ deforms bijectively and regularly the plane and allows to build a deformed random field $X\circ\theta\,:\,\mathbb{R}^2\rightarrow\mathbb{R}$ from a regular, stationary and isotropic random field $X\,:\,\mathbb{R}^2\rightarrow\mathbb{R}$. The deformed field $X\circ\theta$ is in general not isotropic, however we give an explicit characterization of the deformations $\theta$ that preserve the isotropy. Further assuming that $X$ is Gaussian, we introduce a weak form of isotropy of the field $X\circ\theta$, defined by an invariance property of the mean Euler characteristic of some of its excursion sets. Deformed fields satisfying this property are proved to be strictly isotropic. Besides, assuming that the mean Euler characteristic of excursions sets of $X\circ\theta$ over some basic domains is known, we are able to identify $\theta$.


title: Calibrated Entanglement Entropy

abstract: The Ryu-Takayanagi prescription reduces the problem of calculating entanglement entropy in CFTs to the determination of minimal surfaces in a dual anti-de Sitter geometry. For 3D gravity theories and BTZ black holes, we identify the minimal surfaces as special Lagrangian cycles calibrated by the real part of the holomorphic one-form of a spacelike hypersurface. We show that (generalised) calibrations provide a unified way to determine holographic entanglement entropy from minimal surfaces, which is applicable to warped AdS$_3$ geometries. We briefly discuss generalisations to higher dimensions.


title: Explaining Transition Systems through Program Induction

abstract: Explaining and reasoning about processes which underlie observed black-box phenomena enables the discovery of causal mechanisms, derivation of suitable abstract representations and the formulation of more robust predictions. We propose to learn high level functional programs in order to represent abstract models which capture the invariant structure in the observed data. We introduce the $\pi$-machine (program-induction machine) -- an architecture able to induce interpretable LISP-like programs from observed data traces. We propose an optimisation procedure for program learning based on backpropagation, gradient descent and A* search. We apply the proposed method to three problems: system identification of dynamical systems, explaining the behaviour of a DQN agent and learning by demonstration in a human-robot interaction scenario. Our experimental results show that the $\pi$-machine can efficiently induce interpretable programs from individual data traces.


title: Increasing Papers' Discoverability with Precise Semantic Labeling: the sci.AI Platform

abstract: The number of published findings in biomedicine increases continually. At the same time, specifics of the domain's terminology complicates the task of relevant publications retrieval. In the current research, we investigate influence of terms' variability and ambiguity on a paper's likelihood of being retrieved. We obtained statistics that demonstrate significance of the issue and its challenges, followed by presenting the sci.AI platform, which allows precise terms labeling as a resolution.


title: Behavior of digital sequences through exotic numeration systems

abstract: Many digital functions studied in the literature, e.g., the summatory function of the base-$k$ sum-of-digits function, have a behavior showing some periodic fluctuation. Such functions are usually studied using techniques from analytic number theory or linear algebra. In this paper we develop a method based on exotic numeration systems and we apply it on two examples motivated by the study of generalized Pascal triangles and binomial coefficients of words.


title: Soil Property and Class Maps of the Conterminous US at 100 meter Spatial Resolution based on a Compilation of National Soil Point Observations and Machine Learning

abstract: With growing concern for the depletion of soil resources, conventional soil data must be updated to support spatially explicit human-landscape models. Three US soil point datasetswere combined with a stack of over 200 environmental datasets to generate complete coverage gridded predictions at 100 m spatial resolution of soil properties (percent organic C, total N, bulk density, pH, and percent sand and clay) and US soil taxonomic classes (291 great groups and 78 modified particle size classes) for the conterminous US. Models were built using parallelized random forest and gradient boosting algorithms. Soil property predictions were generated at seven standard soil depths (0, 5, 15, 30, 60, 100 and 200 cm). Prediction probability maps for US soil taxonomic classifications were also generated. Model validation results indicate an out-of-bag classification accuracy of 60 percent for great groups, and 66 percent for modified particle size classes; for soil properties cross-validated R-square ranged from 62 percent for total N to 87 percent for pH. Nine independent validation datasets were used to assess prediction accuracies for soil class models and results ranged between 24-58 percent and 24-93 percent for great group and modified particle size class prediction accuracies, respectively. The hybrid "SoilGrids+" modeling system that incorporates remote sensing data, local predictions of soil properties, conventional soil polygon maps, and machine learning opens the possibility for updating conventional soil survey data with machine learning technology to make soil information easier to integrate with spatially explicit models, compared to multi-component map units.


title: Computational simulation of bone remodelling post reverse total shoulder arthroplasty

abstract: Bone is a living material. It adapts, in an optimal sense, to loading by changing its density and trabeculae architecture - a process termed remodelling. Implanted orthopaedic devices can significantly alter the loading on the surrounding bone, which can have a detrimental impact on bone ingrowth that is critical to ensure secure implant fixation. In this contribution, a computational model that accounts for bone remodelling is developed to elucidate the response of bone following a reverse shoulder procedure for rotator cuff deficient patients. The physical process of remodelling is modelled using continuum scale, open system thermodynamics whereby the density of bone evolves isotropically in response to the loading it experiences. The fully-nonlinear continuum theory is solved approximately using the finite element method. The code developed to model the reverse shoulder procedure is validated using a series of benchmark problems.


title: New proofs of Perelman's theorem on shrinking Breathers in Ricci flow

abstract: We give two new proofs of Perelman's theorem that shrinking breathers of Ricci flow on closed manifolds are gradient Ricci solitons, using the fact that the singularity models of type I solutions are shrinking gradient Ricci solitons and the fact that non-collapsed type I ancient solutions have rescaled limits being shrinking gradient Ricci solitons.


title: Quantum K\"ahlerian Lie groups from multiplicative unitaries

abstract: We show that the deformation theory of Fr\'echet algebras for actions of K\"ahlerian Lie groups developed by two of us, leads in a natural way to examples of non-compact locally compact quantum groups. This is achieved by constructing a manageable multiplicative unitary out of the Fr\'echet deformation of $C_0(G)$ for the action $\lambda\otimes \rho$ of $G\times G$ and the undeformed coproduct. We also prove that these quantum groups are isomorphic to those constructed out of the unitary dual $2$-cocycle discovered by Neshveyev and Tuset and associated with Bieliavsky's covariant $\star$-product, via the De Commer's results.


title: A dynamical modeling to study the adaptive immune system and the influence of antibodies in the immune memory

abstract: Immunological systems have been an abundant inspiration to contemporary computer scientists. Problem solving strategies, stemming from known immune system phenomena, have been successfully applied to challenging problems of modern computing (MONROY, SAAB, GODINEZ, 2004). Simulation systems and mathematical modeling are also beginning use to answer more complex immunological questions as immune memory process and duration of vaccines, where the regulation mechanisms are not still known sufficiently (LundegaarD, Lund, Kesmir, Brunak, Nielsen, 2007).In this article we studied in machina a approach to simulate the process of antigenic mutation and its implications for the process of memory. Our results have suggested that the durability of the immune memory is affected by the process of antigenic mutation and by populations of soluble antibodies in the blood. The results also strongly suggest that the decrease of the production of antibodies favors the global maintenance of immune memory.


title: The random field Blume-Capel model revisited

abstract: We have revisited the mean-field treatment for the Blume-Capel model under the presence of a discrete random magnetic field as introduced by Kaufman and Kanner. The magnetic field ($H$) versus temperature ($T$) phase diagrams for given values of the crystal field $D$ were recovered in accordance to Kaufman and Kanner original work. However, our main goal in the present work was to investigate the distinct structures of the crystal field versus temperature phase diagrams as the random magnetic field is varied because similar models have presented reentrant phenomenon due to randomness. Following previous works we have classified the distinct phase diagrams according to five different topologies. The topological structure of the phase diagrams is maintained for both $H-T$ and $D-T$ cases. Althoug the phase diagrams exhibit a richness of multicritical phenomena we did not found any reentrant effect as have been seen in similar models.


title: Mass differences and neutron pairing in Ca, Sn and Pb isotopes

abstract: Various estimates of the even-odd effect of the mass shell of atomic nuclei are considered. Based on the experimental mass values of the Ca, Sn, and Pb isotopes, the dependence of the energy gap on the neutrons number is traced and the relationship of this characteristic to the properties of external neutron subshells is shown. In nuclei with closed proton shells, effects directly related to neutron pairing and effects of nucleon shells are discussed.


title: Crystallites in Color Glass Beads of the 19th Century and Their Influence on Fatal Deterioration of Glass

abstract: Glass corrosion is a crucial problem in keeping and conservation of beadworks in museums. All kinds of glass beads undergo deterioration but blue-green lead-potassium glass beads of the 19th century are subjected to the destruction to the greatest extent. Blue-green lead-potassium glass beads of the 19th century obtained from exhibits kept in Russian museums were studied with the purpose to determine the causes of the observed phenomenon. For the comparison, yellow lead beads of the 19th century were also explored. Both kinds of beads contain Sb but yellow ones are stable. Using scanning electron microscopy, energy dispersive X-ray microspectrometry, electron backscatter diffraction, transmission electron microscopy and X-ray powder analysis, we have registered the presence of crystallites of orthorhombic KSbOSiO$_4$ and cubic Pb$_2$Sb$_{1.5}$Fe$_{0.5}$O$_{6.5}$ in glass matrix of blue-green and yellow beads, respectively. Both compounds form at rather high temperatures obviously during glass melting and/or melt cooling. We suppose that the crystallites generate internal tensile strain in glass during its cooling which causes formation of multiple microcracks in inner domains of blue-green beads. We suggest that the deterioration degree depends on quantity of the precipitates, their sizes and their temperature coefficients of linear expansion. In blue-green beads, the crystallites are distributed in their sizes from $\sim\,$200 nm to several tens of $\mu$m and tend to gather in large colonies. The sizes of crystallites in yellow beads are several hundreds of nm and their clusters contain few crystallites. This explains the difference in corrosion of these kinds of beads containing crystals of Sb compounds.


title: Exact adaptive confidence intervals for linear regression coefficients

abstract: We propose an adaptive confidence interval procedure (CIP) for the coefficients in the normal linear regression model. This procedure has a frequentist coverage rate that is constant as a function of the model parameters, yet provides smaller intervals than the usual interval procedure, on average across regression coefficients. The proposed procedure is obtained by defining a class of CIPs that all have exact $1-\alpha$ frequentist coverage, and then selecting from this class the procedure that minimizes a prior expected interval width. Such a procedure may be described as "frequentist, assisted by Bayes" or FAB. We describe an adaptive approach for estimating the prior distribution from the data so that exact non-asymptotic $1-\alpha$ coverage is maintained. Additionally, in a "$p$ growing with $n$" asymptotic scenario, this adaptive FAB procedure is asymptotically Bayes-optimal among $1-\alpha$ frequentist CIPs.


title: TRILEX and $GW$+EDMFT approach to $d$-wave superconductivity in the Hubbard model

abstract: We generalize the recently introduced TRILEX approach (TRiply Irreducible Local EXpansion) to superconducting phases. The method treats simultaneously Mott and spin-fluctuation physics using an Eliashberg theory supplemented by local vertex corrections determined by a self-consistent quantum impurity model. We show that, in the two-dimensional Hubbard model, at strong coupling, TRILEX yields a $d$-wave superconducting dome as a function of doping. Contrary to the standard cluster dynamical mean field theory (DMFT) approaches, TRILEX can capture $d$-wave pairing using only a single-site effective impurity model. We also systematically explore the dependence of the superconducting temperature on the bare dispersion at weak coupling, which shows a clear link between strong antiferromagnetic (AF) correlations and the onset of superconductivity. We identify a combination of hopping amplitudes particularly favorable to superconductivity at intermediate doping. Finally, we study within $GW$+EDMFT the low-temperature $d$-wave superconducting phase at strong coupling in a region of parameter space with reduced AF fluctuations.


title: A Note on Uniform Integrability of Random Variables in a Probability Space and Sublinear Expectation Space

abstract: In this note we discuss uniform integrability of random variables. In a probability space, we introduce two new notions on uniform integrability of random variables, and prove that they are equivalent to the classic one. In a sublinear expectation space, we give de La Vall\'ee Poussin criterion for the uniform integrability of random variables and do some other discussions.


title: Geometrical frustration yields fiber formation in self-assembly

abstract: Controlling the self-assembly of supramolecular structures is vital for living cells, and a central challenge for engineering at the nano- and microscales. Nevertheless, even particles without optimized shapes can robustly form well-defined morphologies. This is the case in numerous medical conditions where normally soluble proteins aggregate into fibers. Beyond the diversity of molecular mechanisms involved, we propose that fibers generically arise from the aggregation of irregular particles with short-range interactions. Using a minimal model of ill-fitting, sticky particles, we demonstrate robust fiber formation for a variety of particle shapes and aggregation conditions. Geometrical frustration plays a crucial role in this process, and accounts for the range of parameters in which fibers form as well as for their metastable character.


title: Quantum Bianchi identities and characteristic classes via DG categories

abstract: We show how DG categories arise naturally in noncommutative differential geometry and use them to derive noncommutative analogues of the Bianchi identities for the curvature of a connection. We also give a derivation of formulae for characteristic classes in noncommutative geometry following Chern's original derivation, rather than using cyclic cohomology. We show that a related DG category for extendable bimodule connections is a monoidal tensor category and in the metric compatible case give an analogue of a classical antisymmetry of the Riemann tensor. The monoidal structure implies the existence of a cup product on noncommutative sheaf cohomology. Another application is to prove that the curvature of a line module reduces to a 2-form on the base algebra. We also extend our geometric approach to Dirac operators. We illustrate the theory on the q-sphere, the permutation group S_3 and the bicrossproduct model quantum spacetime with algebra [r,t]=\lambda r.


title: A Fluid-Structure Interaction Solver for Compressible Flows with Applications in Blast Loading on Thin Elastic Structures

abstract: We report development and application of a fluid-structure interaction (FSI) solver for compressible flows with large-scale flow-induced deformation of the structure. The FSI solver utilizes partitioned approach to strongly couple a sharp-interface immersed boundary method based flow solver with an open-source finite-element structure dynamics solver. The flow solver is based on a higher-order finite-difference method on Cartesian grid and employs ghost-cell methodology to impose boundary conditions on the immersed boundary. A higher-order accuracy near the immersed boundary is achieved by combining the ghost-cell approach with a weighted least-square error method based on a higher-order approximate polynomial. The second order spatial accuracy of the flow solver is established by performing a grid refinement study. The structure solver is validated with a canonical elastostatics problem. The FSI solver is validated with published measurements and simulations for the large-scale deformation of a thin elastic steel panel subjected to blast loading in a shock tube. The solver correctly predicts oscillating behavior of the tip of the panel with reasonable fidelity and computed shock wave propagation is qualitatively consistent with the published results. In order to demonstrate the fidelity of the solver and to investigate coupled physics of the shock-structure interaction for a thin elastic plate, we employ the solver for simulating 6.4 kg TNT blast loading on the thin elastic plate. The initial conditions of the blast are taken from field tests reported in the literature. Using numerical schlieren, the shock front propagation, Mach reflection and vortex shedding at the tip of the plate are visualized during the shock wave impact on the plate. We discuss coupling between the non-linear dynamics of the plate and blast loading.


title: Identification of Gamma-Rays and Neutrinos from the Cygnus-X Complex Considering Radio Gamma Correlation

abstract: The abundance of accelerators and the ambient conditions make Cygnus X a natural laboratory for studying the life cycle of cosmic-rays (CRs). This naturally makes the Cygnus X complex a highly interesting source in neutrino astronomy, in particular concerning a possible detection with the IceCube Neutrino Observatory, which has a good view of the northern hemisphere. In this paper, we model the multiwavelength spectrum of the Cygnus, for the first time using a broad data set from radio, MeV (COMPTEL), GeV (Fermi), TeV (Argo) and 10s of TeV (Milagro) energies. The modeling is performed assuming a leptohadronic model. We solve the steady-state transport equation for leptons and hadrons injected homogeneously in the region and test the role of diffusive transport and energy loss by radiation and interaction. The result shows that diffusion loss plays a significant role in Cygnus X and always exceeds the advection loss as well as almost all other loss processes. The best-fit parameters we find are a magnetic field of $B=8.9\times10^{-6}$ G, a target density of $N_t=19.4$ cm$^{-3}$, a cosmic ray spectral index of $\alpha=2.37$ and neutral gas distribution over a depth of 116 pc. We find that the fit describes the data up to TeV energies well, while the Milagro data are underestimated. This transport model with a broad multiwavelength fit provides a neutrino flux which approaches the sensitivity of IceCube at very high energies ($>$ 50 TeV). In the future, the flux sensitivity of IceCube will be improved. With this rather pessimistic model, leaving out the influence of possible strong, high-energy point sources, we already expect the flux in the Cygnus X region to suffice for IceCube to measure a significant neutrino flux in the next decade.


title: A Colonel Blotto Game for Interdependence-Aware Cyber-Physical Systems Security in Smart Cities

abstract: Smart cities must integrate a number of interdependent cyber-physical systems that operate in a coordinated manner to improve the well-being of the city's residents. A cyber-physical system (CPS) is a system of computational elements controlling physical entities. Large-scale CPSs are more vulnerable to attacks due to the cyber-physical interdependencies that can lead to cascading failures which can have a significant detrimental effect on a city. In this paper, a novel approach is proposed for analyzing the problem of allocating security resources, such as firewalls and anti-malware, over the various cyber components of an interdependent CPS to protect the system against imminent attacks. The problem is formulated as a Colonel Blotto game in which the attacker seeks to allocate its resources to compromise the CPS, while the defender chooses how to distribute its resources to defend against potential attacks. To evaluate the effects of defense and attack, various CPS factors are considered including human-CPS interactions as well as physical and topological characteristics of a CPS such as flow and capacity of interconnections and minimum path algorithms. Results show that, for the case in which the attacker is not aware of the CPS interdependencies, the defender can have a higher payoff, compared to the case in which the attacker has complete information. The results also show that, in the case of more symmetric nodes, due to interdependencies, the defender achieves its highest payoff at the equilibrium compared to the case with independent, asymmetric nodes.


title: Distributed Precoding Systems in Multi-Gateway Multibeam Satellites: Regularization and Coarse Beamforming

abstract: This paper deals with the problem of beamforming design in a multibeam satellite, which is shared by different groups of terminals -clusters-, each served by an Earth station or gateway. Each gateway precodes the symbols addressed to its respective users; the design follows an MMSE criterion, and a regularization factor judiciously chosen allows to account for the presence of mutually interfering clusters, extending more classical results applicable to one centralized station. More importantly, channel statistics can be used instead of instantaneous channel state information, avoiding the exchange of information among gateways through backhaul links. The on-board satellite beamforming weights are designed to exploit the degrees of freedom of the satellite antennas to minimize the noise impact and the interference to some specific users. On-ground beamforming results are provided as a reference to compare the joint performance of MMSE precoders and on-board beamforming network. A non-adaptive design complements the results and makes them more amenable to practical use by designing a coarse beamforming network.


title: On random stable partitions

abstract: The stable roommates problem does not necessarily have a solution, i.e. a stable matching. We had found that, for the uniformly random instance, the expected number of solutions converges to $e^{1/2}$ as $n$, the number of members, grows, and with Rob Irving we proved that the limiting probability of solvability is $e^{1/2}/2$, at most. Stephan Mertens's extensive numerics compelled him to conjecture that this probability is of order $n^{-1/4}$. Jimmy Tan introduced a notion of a stable cyclic partition, and proved existence of such a partition for every system of members' preferences, discovering that presence of odd cycles in a stable partition is equivalent to absence of a stable matching. In this paper we show that the expected number of stable partitions with odd cycles grows as $n^{1/4}$. However the standard deviation of that number is of order $n^{3/8}\gg n^{1/4}$, too large to conclude that the odd cycles exist with high probability (whp). Still, as a byproduct, we show that whp the fraction of members with more than one stable "predecessor" is of order $n^{-1/4}$. Furthermore, whp the average rank of a predecessor in every stable partition is of order $n^{1/2}$. The likely size of the largest stable matching is $n/2-O(n^{1/4+o(1)})$, and the likely number of pairs of unmatched members blocking the optimal complete matching is $O(n^{3/4+o(1)})$.


title: An impossibility theorem for Parameter Independent hidden variable theories

abstract: Recently, Roger Colbeck and Renato Renner (C&R) have claimed that '[n]o extension of quantum theory can have improved predictive power'. If correct, this is a spectacular impossibility theorem for hidden variable theories, which is more general than the theorems of Bell and Leggett. C&R's claim essentially means that in any hidden variable theory that is compatible with quantum-mechanical predictions, probabilities of measurement outcomes are independent of these hidden variables. On closer inspection, however, the generality and validity of the claim can be contested. First, it is based on an assumption called 'Freedom of Choice'. As the name suggests, this assumption involves the independence of an experimenter's choice of measurement settings. But in the way C&R define this assumption, a no-signalling condition is surreptitiously presupposed, making the assumption less innocent than it sounds. When using this definition, any hidden variable theory violating Parameter Independence, such as Bohmian Mechanics, is immediately shown to be incompatible with quantum-mechanical predictions. Also, the argument of C&R is hard to follow and their mathematical derivation contains several gaps, some of which cannot be closed in the way they suggest. We shall show that these gaps can be filled. The issue with the 'Freedom of Choice' assumption can be circumvented by explicitly assuming Parameter Independence. This makes the result less general, but better founded. We then obtain an impossibility theorem for hidden variable theories satisfying Parameter Independence only. So, while quantum mechanics itself satisfies Parameter Independence, if a variable is added that changes the outcome probabilities, however slightly, Parameter Independence must be violated.


title: A Finite Difference Scheme based on Cubic Trigonometric B-splines for Time Fractional Diffusion-wave Equation

abstract: In this paper, we propose an efficient numerical scheme for the approximate solution of the time fractional diffusion-wave equation with reaction term based on cubic trigonometric basis functions. The time fractional derivative is approximated by the usual finite difference formulation and the derivative in space is discretized using cubic trigonometric B-spline functions. A stability analysis of the scheme is conducted to confirm that the scheme does not amplify errors. Computational experiments are also performed to further establish the accuracy and validity of the proposed scheme. The results obtained are compared with a finite difference schemes based on the Hermite formula and radial basis functions. It is found that our numerical approach performs superior to the existing methods due to its simple implementation, straight forward interpolation and very less computational cost.


title: Counting the number of non-zero coefficients in rows of generalized Pascal triangles

abstract: This paper is about counting the number of distinct (scattered) subwords occurring in a given word. More precisely, we consider the generalization of the Pascal triangle to binomial coefficients of words and the sequence $(S(n))_{n\ge 0}$ counting the number of positive entries on each row. By introducing a convenient tree structure, we provide a recurrence relation for $(S(n))_{n\ge 0}$. This leads to a connection with the $2$-regular Stern-Brocot sequence and the sequence of denominators occurring in the Farey tree. Then we extend our construction to the Zeckendorf numeration system based on the Fibonacci sequence. Again our tree structure permits us to obtain recurrence relations for and the F-regularity of the corresponding sequence.


title: Dissipative systems entropy

abstract: In this paper, we introduce the generalized phase space $\left( \vec{r},\vec{v},\dot{\vec{v}},\ddot{\vec{v}},... \right)$, which expands the known phase space $\left( \vec{r},\vec{v} \right)$. The fact is that the introduced space is the infinity dimensional phase space. The paper shows that dissipative systems in generalized phase space may be considered as conservative systems. It is shown that, in infinity dimensional phase space, the entropy is a constant value. It is shown that the transition to finite dimensional phase space leads to dissipations and change of the entropy. The paper contains the rigorous mathematical result.


title: Observation of the motional Stark shift in low magnetic fields

abstract: We report on the observation of the motional Stark effect of highly excited $^{87}$Rb Rydberg atoms moving in the presence of a weak homogeneous magnetic field in a vapor cell. Employing electromagnetically induced transparency for spectroscopy of an atomic vapor, we observe the velocity-, quantum state- and magnetic field-dependent transition frequencies between the ground and Rydberg excited states. For atoms moving at velocities around $400\,\mathrm{m/s}$, the principal quantum number $n=100$ of the valence electron, and a magnetic field of $B=100\,\mathrm{G}$, we measure a motional Stark shift of $\sim10\,\mathrm{MHz}$. Our experimental results are supported by numerical calculations based on a diagonalization of the effective Hamiltonian governing the valence electron of $^{87}$Rb in the presence of crossed electric and magnetic fields.


title: Limiting Majoron self-interactions from Gravitational waves experiments

abstract: We show how Majoron models may be tested/limited in gravitational waves experiments. In particular, the Majoron self-interaction potential may induce a first order phase transition, producing gravitational waves from bubble collisions. We dubbed such a new scenario {\it violent Majoron model}, because it would be associated to a violent phase transition in the early Universe. Sphaleron constraints can be avoided if the global $U(1)_{B-L}$ is broken at scales lower than the electroweak scale, provided that the B-L spontaneously breaking scale is lower than $10\, {\rm TeV}$ in order to satisfy the cosmological mass density bound. The possibility of a sub-electroweak phase transition is practically unconstrained by cosmological bounds and it may be detected within the sensitivity of next generation of gravitational waves experiments: eLISA, DECIGO and BBO. We also comment on the possible detection in CEPC collider, where Majorons's production can be observed from Higgs' portals in missing transverse energy channels.


title: Quantum Monte Carlo calculation of the bridge reaction barrier to hydrogen dissociation on Pt(111)

abstract: The yardstick of new first-principles approaches to key points on reaction paths at metal surfaces is chemical accuracy compared to reliable experiment. By this we mean that such values as the activation barrier are required to within 1 kcal/mol. Quantum Monte Carlo (QMC) is a promising (albeit lengthy) first-principles method for this and we are now beyond the dawn of QMC benchmarks for these systems, since hydrogen dissociation on Cu(111) has been studied with quite adequate accuracy in two improving QMC studies Hoggan ArXiv 2015, K. Doblhoff-Deir et al JCTC (2017). Pt and Cu require the use of pseudo-potentials in these calculations and we show that those of Pt are less problematic than those for Cu, particularly for QMC work. In this letter, we determine the activation barrier to hydrogen dissociation on Pt(111) using the bridging geometry. This is found to agree to better than chemical accuracy with recent Specific Reaction Parameter (SRP-DFT) work, which is itself within chemical accuracy of measurement (Kroes et al Cplett 2017). The results suggest the dissociation barrier for hydrogen on Pt(111) is close to 5.4 (QMC) and 6.2 (SRP-DFT) kcal/mol with a QMC standard error of 1.08 kcal/mol. This is encouraging for establishing less well-known benchmark values of industrial reaction barriers on Pt(111).


title: A variational principle for problems with a hint of convexity

abstract: A variational principle is introduced to provide a new formulation and resolution for several boundary value problems with a variational structure. This principle allows one to deal with problems well beyond the weakly compact structure. As a result, we study several super-critical semilinear Elliptic problems.


title: Fast Solution of the Linearized Poisson-Boltzmann Equation with nonaffine Parametrized Boundary Conditions Using the Reduced Basis Method

abstract: The Poisson-Boltzmann equation (PBE) is a nonlinear elliptic PDE that arises in biomolecular modeling and is a fundamental tool for structural biology. It is used to calculate electrostatic potentials around an ensemble of fixed charges immersed in an ionic solution. Efficient numerical computation of the PBE yields a high number of degrees of freedom in the resultant algebraic system of equations. Coupled with the fact that in most cases the PBE requires to be solved multiple times for a large number of system configurations, this poses great computational challenges to conventional numerical techniques. To accelerate such computations, we here present the reduced basis method (RBM) which greatly reduces this computational complexity by constructing a reduced order model of typically low dimension. In this study, we employ a simple version of the PBE for proof of concept and discretize the linearized PBE (LPBE) with a centered finite difference scheme. The resultant linear system is solved by the aggregation-based algebraic multigrid method at different samples of ionic strength on a three-dimensional Cartesian grid. The discretized LPBE, which we call the high-fidelity full order model (FOM), yields solution as accurate as other LPBE solvers. We then apply the RBM to FOM. The discrete empirical interpolation method (DEIM) is applied to the Dirichlet boundary conditions which are nonaffine with the parameter, to reduce the complexity of the reduced order model (ROM). From the numerical results, we notice that the RBM reduces the model order from $\mathcal{N} = 2\times 10^{6}$ to $N = 6$ at an accuracy of $10^{-9}$ and reduces computational time by a factor of approximately $7,600$. DEIM, on the other hand, is also used in the offline-online phase of solving the ROM for different values of parameters which provides a speed-up of $20$ for a single iteration of the greedy algorithm.


title: Bounding Cache Miss Costs of Multithreaded Computations Under General Schedulers

abstract: We analyze the caching overhead incurred by a class of multithreaded algorithms when scheduled by an arbitrary scheduler. We obtain bounds that match or improve upon the well-known $O(Q+S \cdot (M/B))$ caching cost for the randomized work stealing (RWS) scheduler, where $S$ is the number of steals, $Q$ is the sequential caching cost, and $M$ and $B$ are the cache size and block (or cache line) size respectively.


title: Dynamics of a spherically symmetric inhomogeneous coupled dark energy model with coupling term proportional to non relatvistic matter

abstract: Quasi--local scalar variables approach is applied to a spherically symmetric inhomogeneous Lema\^\i tre--Tolman--Bondi metric containing a mixture of non-relativistic cold dark matter and coupled dark energy with constant equation of state. The quasi--local coupling term considered is proportional to the quasi--local cold dark matter energy density and a quasi--local Hubble factor-like scalar via a coupling constant $\alpha$. The autonomous numerical system obtained from the evolution equations is classified for different choices of the free parameters: the adiabatic constant of the dark energy $w$ and $\alpha$. The presence of a past attractor in a non-physical region of the energy densities phase-space of the system makes the coupling term non physical when the energy flows from the matter to the dark energy in order to avoid negative values of the dark energy density in the past. On the other hand, if the energy flux goes from dark energy to dark matter, the past attractor lays in a physical region. The system is also numerically solved for some interesting initial profiles leading to different configurations: an ever expanding mixture, a scenario where the dark energy is completely consumed by the non-relativistic matter by means of the coupling term, a scenario where the dark energy disappears in the inner layers while the outer layers expand as a mixture of both sources, and, finally, a structure formation toy model scenario, where the inner shells containing the mixture collapse while the outer shells expand.


title: A natural linear equation in affine geometry: the affine quasi-Einstein equation

abstract: We study the affine quasi-Einstein equation, a second order linear homogeneous equation, which is invariantly defined on any affine manifold. We prove that the space of solutions is finite-dimensional, and its dimension is a strongly projective invariant. Moreover the maximal dimension is shown to be achieved if and only if the manifold is strongly projectively flat.


title: Classifying Galois groups of small iterates via rational points

abstract: We establish several surjectivity theorems regarding the Galois groups of small iterates of $\phi_c(x)=x^2+c$ for $c\in\mathbb{Q}$. To do this, we use explicit techniques from the theory of rational points on curves, including the method of Chabauty-Coleman and the Mordell-Weil sieve. For example, we succeed in finding all rational points on a hyperelliptic curve of genus $7$, with rank $5$ Jacobian, whose points parametrize quadratic polynomials with a "newly small" Galois group at the fifth stage of iteration.


title: The Rafita asteroid family

abstract: The Rafita asteroid family is an S-type group located in the middle main belt, on the right side of the 3J:-1A mean-motion resonance. The proximity of this resonance to the family left side in semi-major axis caused many former family members to be lost. As a consequence, the family shape in the $(a,1/D)$ domain is quite asymmetrical, with a preponderance of objects on the right side of the distribution. The Rafita family is also characterized by a leptokurtic distribution in inclination, which allows the use of methods of family age estimation recently introduced for other leptokurtic families such as Astrid, Hansa, Gallia, and Barcelona. In this work we propose a new method based on the behavior of an asymmetry coefficient function of the distribution in the $(a,1/D)$ plane to date incomplete asteroid families such as Rafita. By monitoring the time behavior of this coefficient for asteroids simulating the initial conditions at the time of the family formation, we were able to estimate that the Rafita family should have an age of $490\pm200$ Myr, in good agreement with results from independent methods such as Monte Carlo simulations of Yarkovsky and Yorp dynamical induced evolution and the time behaviour of the kurtosis of the $\sin{(i)}$ distribution. Asteroids from the Rafita family can reach orbits similar to 8\% of the currently known near Earth objects. $\simeq$1\% of the simulated objects are present in NEO-space during the final 10 Myr of the simulation, and thus would be comparable to objects in the present-day NEO population.


title: One year of monitoring the Vela pulsar using a Phased Array Feed

abstract: We have observed the Vela pulsar for one year using a Phased Array Feed (PAF) receiver on the 12-metre antenna of the Parkes Test-Bed Facility. These observations have allowed us to investigate the stability of the PAF beam-weights over time, to demonstrate that pulsars can be timed over long periods using PAF technology and to detect and study the most recent glitch event that occurred on 12 December 2016. The beam-weights are shown to be stable to 1% on time scales on the order of three weeks. We discuss the implications of this for monitoring pulsars using PAFs on single dish telescopes.


title: The influence of large scale magnetic field in the structure of supercritical accretion flow with outflow

abstract: We present the effects of ordered large scale magnetic field on the structure of supercritical accretion flow in the presence of outflow. In the cylindrical coordinates ($r, \varphi,z$), We write the 1.5 dimensional, the steady state ($\frac{\partial}{\partial t}=0$) and axisymmetric ($\frac{\partial}{\partial \varphi}=0$) inflow-outflow equations by using the self similar solutions. Also a model for radiation pressure supported accretion flow threaded by both toroidal and vertical components of magnetic field has been formulated. For studying the outflows, we adopt a radius dependent mass accretion rate as $\dot{M}=\dot{M}_{out}{(\frac{r}{r_{out}})^{s+\frac{1}{2}}}$ with $s=\frac{1}{2}$. Also by following the previous works, we have considered the interchange of mass, radial and angular momentum and the energy between inflow and outflow. we have found numerically that two components of magnetic field have the opposite effects on the thickness of the disc and similar effects on the radial and angular velocities of the flow. We have found that the existence of the toroidal component of magnetic field will lead to increasing of radial and azimuthal velocities as well as the relative thickness of the disc, which is increased . Moreover, the thickness of the disc will decrease when the vertical component of magnetic field becomes important in magnetized flow. The solutions indicated that the mass inflow rate, the specific energy of outflow affect strongly on the advection parameter.We have shown that by increasing the two components of magnetic field, the temperature of the accretion flow will decrease significantly. On the other hand we have shown the bolometric luminosity of the slim discs for high values of $\dot{m} (\dot{m}>>1)$ is not sensitive to mass accretion rate and is kept constant ($L \approx 10 L_{E}$).


title: New nanostructures of carbon: Quasifullerenes Cn-q (n=20,42,48,60)

abstract: Based on the third allotropic form of carbon (Fullerenes) through theoretical study have been predicted structures described as non-classical fullerenes. We have studied novel allotropic carbon structures with a closed cage configuration that have been predicted for the first time, by using DFT at the B3LYP level. Such carbon Cn-q structures (where, n=20, 42, 48 and 60), combine states of hybridization sp1 and sp2, for the formation of bonds. A comparative analysis of quasi-fullerenes with respect to their isomers of greater stability was also performed. Chemical stability was evaluated with the criteria of aromaticity through the different rings that build the systems. The results show new isomerism of carbon nanostructures with interesting chemical properties such as hardness, chemical potential and HOMO-LUMO gaps. We also studied thermal stability with Lagrangian molecular dynamics method using Atom- Center Density propagation (ADMP) method.


title: Spin mixing and protection of ferromagnetism in a spinor dipolar condensate

abstract: We study spin mixing dynamics in a chromium dipolar Bose-Einstein Condensate, after tilting the atomic spins by an angle $\theta$ with respect to the magnetic field. Spin mixing is triggered by dipolar coupling, but, once dynamics has started, it is mostly driven by contact interactions. For the particular case $\theta=\pi/2$, an external spin-orbit coupling term induced by a magnetic gradient is required to enable the dynamics. Then the initial ferromagnetic character of the gas is locally preserved, an unexpected feature that we attribute to large spin-dependent contact interactions.


title: Discovery of the first quadruple gravitationally lensed quasar candidate with Pan-STARRS

abstract: We report the serendipitous discovery of the first gravitationally lensed quasar candidate from Pan-STARRS. The grizy images reveal four point-like images with magnitudes between 14.9 mag and 18.1 mag. The colors of the point sources are similar, and they are more consistent with quasars than with stars or galaxies. The lensing galaxy is detected in the izy bands, with an inferred photometric redshift of ~0.6, lower than that of the point sources. We successfully model the system with a singular isothermal ellipsoid with shear, using the relative positions of the five objects as constraints. While the brightness ranking of the point sources is consistent with that of the model, we find discrepancies between the model-predicted and observed fluxes, likely due to microlensing by stars and millilensing due to the dark matter substructure. In order to fully confirm the gravitational lens nature of this system and add it to the small but growing number of the powerful probes of cosmology and astrophysics represented by quadruply lensed quasars, we require further spectroscopy and high-resolution imaging.


title: Efficient and principled score estimation with Nystr\"om kernel exponential families

abstract: We propose a fast method with statistical guarantees for learning an exponential family density model where the natural parameter is in a reproducing kernel Hilbert space, and may be infinite-dimensional. The model is learned by fitting the derivative of the log density, the score, thus avoiding the need to compute a normalization constant. Our approach improves the computational efficiency of an earlier solution by using a low-rank, Nystr\"om-like solution. The new solution retains the consistency and convergence rates of the full-rank solution (exactly in Fisher distance, and nearly in other distances), with guarantees on the degree of cost and storage reduction. We evaluate the method in experiments on density estimation and in the construction of an adaptive Hamiltonian Monte Carlo sampler. Compared to an existing score learning approach using a denoising autoencoder, our estimator is empirically more data-efficient when estimating the score, runs faster, and has fewer parameters (which can be tuned in a principled and interpretable way), in addition to providing statistical guarantees.


title: Photonic topological pumping through the edges of a dynamical four-dimensional quantum Hall system

abstract: When a two-dimensional electron gas is exposed to a perpendicular magnetic field and an in-plane electric field, its conductance becomes quantized in the transverse in-plane direction: this is known as the quantum Hall (QH) effect. This effect is a result of the nontrivial topology of the system's electronic band structure, where an integer topological invariant known as the first Chern number leads to the quantization of the Hall conductance. Interestingly, it was shown that the QH effect can be generalized mathematically to four spatial dimensions (4D), but this effect has never been realized for the obvious reason that experimental systems are bound to three spatial dimensions. In this work, we harness the high tunability and control offered by photonic waveguide arrays to experimentally realize a dynamically-generated 4D QH system using a 2D array of coupled optical waveguides. The inter-waveguide separation is constructed such that the propagation of light along the device samples over higher-dimensional momenta in the directions orthogonal to the two physical dimensions, thus realizing a 2D topological pump. As a result, the device's band structure is associated with 4D topological invariants known as second Chern numbers which support a quantized bulk Hall response with a 4D symmetry. In a finite-sized system, the 4D topological bulk response is carried by localized edges modes that cross the sample as a function of of the modulated auxiliary momenta. We directly observe this crossing through photon pumping from edge-to-edge and corner-to-corner of our system. These are equivalent to the pumping of charge across a 4D system from one 3D hypersurface to the opposite one and from one 2D hyperedge to another, and serve as first experimental realization of higher-dimensional topological physics.

categories: quant-ph cond-mat.mes-hall cond-mat.quant-gas physics.optics


title: Efficient Coalgebraic Partition Refinement

abstract: We present a generic partition refinement algorithm that quotients coalgebraic systems by behavioural equivalence, an important task in reactive verification; coalgebraic generality implies in particular that we cover not only classical relational systems but also various forms of weighted systems. Under assumptions on the type functor that allow representing its finite coalgebras in terms of nodes and edges, our algorithm runs in time $\mathcal{O}(m\cdot \log n)$ where $n$ and $m$ are the numbers of nodes and edges, respectively. Instances of our generic algorithm thus match the runtime of the best known algorithms for unlabelled transition systems, Markov chains, and deterministic automata (with fixed alphabets), and improve the best known algorithms for Segala systems.


title: Lifting of Modular Forms

abstract: The existence and construction of vector-valued modular forms (vvmf) for any arbitrary Fuchsian group $\mathrm{G}$, for any representation $\rho:\mathrm{G} \longrightarrow \mathrm{GL}_{d}(\mathbb{C})$ of finite image can be established by lifting scalar-valued modular forms of the finite index subgroup $Ker(\rho)$ of $\mathrm{G}$. In this article vvmf are explicitly constructed for any admissible multiplier (representation) $\rho$, see Section 3 for the definition of admissible multiplier. In other words, the following question has been partially answered: For which representations $\rho$ of a given $\mathrm{G}$, is there a vvmf with at least one nonzero component ?


title: Sharp weighted estimates involving one supremum

abstract: In this note, we study the sharp weighted estimate involving one supremum. In particular, we give a positive answer to an open question raised by Lerner and Moen \cite{LM}. We also extend the result to rough homogeneous singular integral operators.


title: A Short Proof for a Lower Bound on the Zero Forcing Number

abstract: We provide a short proof of a conjecture of Davila and Kenter concerning a lower bound on the zero forcing number $Z(G)$ of a graph $G$. More specifically, we show that $Z(G)\geq (g-2)(\delta-2)+2$ for every graph $G$ of girth $g$ at least $3$ and minimum degree $\delta$ at least $2$.


title: A Bogomolov unobstructedness theorem for log-symplectic manifolds in general position

abstract: We consider compact K\"ahlerian manifolds $X$ of even dimension 4 or more, endowed with a log-symplectic holomorphic Poisson structure $\Pi$ which is sufficiently general, in a precise linear sense, with respect to its (normal-crossing) degeneracy divisor $D(\Pi)$. We prove that $(X, \Pi)$ has unobsrtuced deformations, that the tangent space to its deformation space can be identified in terms of the mixed Hodge structure on $H^2$ of the open symplectic manifold $X\setminus D(\Pi)$, and in fact coincides with this $H^2$ provided the Hodge number $h^{2,0}_X=0$, and finally that the degeneracy locus $D(\Pi)$ deforms locally trivially under deformations of $(X, \Pi)$. It has been pointed out that the general position hypothesis in the original paper is not strong enough and this is corrected in an appended erratum/corrigendum to the revised version.


title: Far-infrared metallicity diagnostics: Application to local ultraluminous infrared galaxies

abstract: The abundance of metals in galaxies is a key parameter which permits to distinguish between different galaxy formation and evolution models. Most of the metallicity determinations are based on optical line ratios. However, the optical spectral range is subject to dust extinction and, for high-z objects (z > 3), some of the lines used in optical metallicity diagnostics are shifted to wavelengths not accessible to ground based observatories. For this reason, we explore metallicity diagnostics using far-infrared (IR) line ratios which can provide a suitable alternative in such situations. To investigate these far-IR line ratios, we modeled the emission of a starburst with the photoionization code CLOUDY. The most sensitive far-IR ratios to measure metallicities are the [OIII]52$\mu$m and 88$\mu$m to [NIII]57$\mu$m ratios. We show that this ratio produces robust metallicities in the presence of an AGN and is insensitive to changes in the age of the ionizing stellar. Another metallicity sensitive ratio is the [OIII]88$\mu$m/[NII]122$\mu$m ratio, although it depends on the ionization parameter. We propose various mid- and far-IR line ratios to break this dependency. Finally, we apply these far-IR diagnostics to a sample of 19 local ultraluminous IR galaxies (ULIRGs) observed with Herschel and Spitzer. We find that the gas-phase metallicity in these local ULIRGs is in the range 0.7 < Z_gas/Z_sun < 1.5, which corresponds to 8.5 < 12 + log (O/H) < 8.9. The inferred metallicities agree well with previous estimates for local ULIRGs and this confirms that they lie below the local mass-metallicity relation.


title: Period proliferation in periodic states in cyclically sheared jammed solids

abstract: Athermal disordered systems can exhibit a remarkable response to an applied oscillatory shear: after a relatively few shearing cycles, the system falls into a configuration that had already been visited in a previous cycle. After this point the system repeats its dynamics periodically despite undergoing many particle rearrangements during each cycle. We study the behavior of orbits as we approach the jamming point in simulations of jammed particles subject to oscillatory shear at fixed pressure and zero temperature. As the pressure is lowered, we find that it becomes more common for the system to find periodic states where it takes multiple cycles before returning to a previously visited state. Thus, there is a proliferation of longer periods as the jamming point is approached.


title: Her2 Challenge Contest: A Detailed Assessment of Automated Her2 Scoring Algorithms in Whole Slide Images of Breast Cancer Tissues

abstract: Evaluating expression of the Human epidermal growth factor receptor 2 (Her2) by visual examination of immunohistochemistry (IHC) on invasive breast cancer (BCa) is a key part of the diagnostic assessment of BCa due to its recognised importance as a predictive and prognostic marker in clinical practice. However, visual scoring of Her2 is subjective and consequently prone to inter-observer variability. Given the prognostic and therapeutic implications of Her2 scoring, a more objective method is required. In this paper, we report on a recent automated Her2 scoring contest, held in conjunction with the annual PathSoc meeting held in Nottingham in June 2016, aimed at systematically comparing and advancing the state-of-the-art Artificial Intelligence (AI) based automated methods for Her2 scoring. The contest dataset comprised of digitised whole slide images (WSI) of sections from 86 cases of invasive breast carcinoma stained with both Haematoxylin & Eosin (H&E) and IHC for Her2. The contesting algorithms automatically predicted scores of the IHC slides for an unseen subset of the dataset and the predicted scores were compared with the 'ground truth' (a consensus score from at least two experts). We also report on a simple Man vs Machine contest for the scoring of Her2 and show that the automated methods could beat the pathology experts on this contest dataset. This paper presents a benchmark for comparing the performance of automated algorithms for scoring of Her2. It also demonstrates the enormous potential of automated algorithms in assisting the pathologist with objective IHC scoring.


title: New explicit expressions for Dirac bilinears

abstract: We derive new explicit expressions for the Dirac bilinears based on a generic representation of the massive Dirac spinors with canonical polarization. These bilinears depend on a direction $n$ in Minkowski space which specifies the form of dynamics. We argue that such a dependence is unavoidable in a relativistic theory with spin, since it originates from Wigner rotation effects. Contrary to most of the expressions found in the literature, our ones are valid for all momenta and canonical polarizations of the spinors. As a by-product, we also obtain a generic explicit expression for the covariant spin vector.


title: Exploring 4D Quantum Hall Physics with a 2D Topological Charge Pump

abstract: The discovery of topological states of matter has profoundly augmented our understanding of phase transitions in physical systems. Instead of local order parameters, topological phases are described by global topological invariants and are therefore robust against perturbations. A prominent example thereof is the two-dimensional integer quantum Hall effect. It is characterized by the first Chern number which manifests in the quantized Hall response induced by an external electric field. Generalizing the quantum Hall effect to four-dimensional systems leads to the appearance of a novel non-linear Hall response that is quantized as well, but described by a 4D topological invariant - the second Chern number. Here, we report on the first observation of a bulk response with intrinsic 4D topology and the measurement of the associated second Chern number. By implementing a 2D topological charge pump with ultracold bosonic atoms in an angled optical superlattice, we realize a dynamical version of the 4D integer quantum Hall effect. Using a small atom cloud as a local probe, we fully characterize the non-linear response of the system by in-situ imaging and site-resolved band mapping. Our findings pave the way to experimentally probe higher-dimensional quantum Hall systems, where new topological phases with exotic excitations are predicted.

categories: cond-mat.quant-gas cond-mat.mes-hall physics.optics quant-ph


title: Coherent many-body spin dynamics in a long-range interacting Ising chain

abstract: Coherent many-body quantum dynamics lies at the heart of quantum simulation and quantum computation. Both require coherent evolution in the exponentially large Hilbert space of an interacting many-body system. To date, trapped ions have defined the state of the art in terms of achievable coherence times in interacting spin chains. Here, we establish an alternative platform by reporting on the observation of coherent, fully interaction-driven quantum revivals of the magnetization in Rydberg-dressed Ising spin chains of atoms trapped in an optical lattice. We identify partial many-body revivals at up to about ten times the characteristic time scale set by the interactions. At the same time, single-site-resolved correlation measurements link the magnetization dynamics with inter-spin correlations appearing at different distances during the evolution. These results mark an enabling step towards the implementation of Rydberg atom based quantum annealers, quantum simulations of higher dimensional complex magnetic Hamiltonians, and itinerant long-range interacting quantum matter.


title: The inverse scattering problem for orthotropic media in Polarization-sensitive Optical Coherence Tomography

abstract: In this paper we provide a mathematical model for imaging an anisotropic, orthotropic medium with Polarization-Sensitive Optical Coherence Tomography (PS-OCT). The imaging problem is formulated as an inverse scattering problem in three dimensions for reconstructing the electrical susceptibility of the medium using Maxwell's equations. Our reconstruction method is based on the second-order Born-approximation of the electric field.


title: Classification of Aerial Photogrammetric 3D Point Clouds

abstract: We present a powerful method to extract per-point semantic class labels from aerialphotogrammetry data. Labeling this kind of data is important for tasks such as environmental modelling, object classification and scene understanding. Unlike previous point cloud classification methods that rely exclusively on geometric features, we show that incorporating color information yields a significant increase in accuracy in detecting semantic classes. We test our classification method on three real-world photogrammetry datasets that were generated with Pix4Dmapper Pro, and with varying point densities. We show that off-the-shelf machine learning techniques coupled with our new features allow us to train highly accurate classifiers that generalize well to unseen data, processing point clouds containing 10 million points in less than 3 minutes on a desktop computer.


title: Extended stirling polynomials of the second kind and extended Bell polynomials

abstract: Recently, several authors have studied the Stirling numbers of the second kind and Bell polynomials. In this paper, we study the extended Stirling polynomials of the second kind and the extended Bell polynomials associated with the Stirling numbers of the second kind. In addition, we note that the extended Bell polynomials can be expressed in terms of the moments of the Poisson random variable with parameter


title: On the failure of lower square function estimates in the non-homogeneous weighted setting

abstract: We show that the classical $A_{\infty}$ condition is not sufficient for a lower square function estimate in the non-homogeneous weighted $L^2$ space. We also show that under the martingale $A_2$ condition, an estimate holds true, but the optimal power of the characteristic jumps from $1 / 2$ to $1$ even when considering the classical $A_2$ characteristic. This is in a sharp contrast to known estimates in the dyadic homogeneous setting as well as the recent positive results in this direction on the discrete timenon-homogeneous martingale transforms. Last, we give a sharp $A_{\infty}$ estimate for the $n$-adic homogeneous case, growing with $n$.


title: Can Dips of Boyajian's Star be Explained by Circumsolar Rings?

abstract: Could the dips of "Boyajian's Star" (KIC 8462852) have been caused by matter in our Solar System? The interval between periods of deep dips is nearly twice the orbital period of the Kepler satellite. I consider heliocentric obscuring rings in the outer Solar System that graze the line of sight to the star once per orbit of Kepler. The hypothesis predicts that future dips may be observed from Earth during windows separated by a year, although their detailed structure depends on the distribution of particles along the ring. Dips observed at terrestrial sites separated by > 1000 km in a direction perpendicular to the Earth's orbital motion may be uncorrelated.


title: Detecting Adversarial Image Examples in Deep Networks with Adaptive Noise Reduction

abstract: Recently, many studies have demonstrated deep neural network (DNN) classifiers can be fooled by the adversarial example, which is crafted via introducing some perturbations into an original sample. Accordingly, some powerful defense techniques were proposed. However, existing defense techniques often require modifying the target model or depend on the prior knowledge of attacks. In this paper, we propose a straightforward method for detecting adversarial image examples, which can be directly deployed into unmodified off-the-shelf DNN models. We consider the perturbation to images as a kind of noise and introduce two classic image processing techniques, scalar quantization and smoothing spatial filter, to reduce its effect. The image entropy is employed as a metric to implement an adaptive noise reduction for different kinds of images. Consequently, the adversarial example can be effectively detected by comparing the classification results of a given sample and its denoised version, without referring to any prior knowledge of attacks. More than 20,000 adversarial examples against some state-of-the-art DNN models are used to evaluate the proposed method, which are crafted with different attack techniques. The experiments show that our detection method can achieve a high overall F1 score of 96.39% and certainly raises the bar for defense-aware attacks.


title: Perfect Edge Domination: Hard and Solvable Cases

abstract: Let $G$ be an undirected graph. An edge of $G$ dominates itself and all edges adjacent to it. A subset $E'$ of edges of $G$ is an edge dominating set of $G$, if every edge of the graph is dominated by some edge of $E'$. We say that $E'$ is a perfect edge dominating set of $G$, if every edge not in $E'$ is dominated by exactly one edge of $E'$. The perfect edge dominating problem is to determine a least cardinality perfect edge dominating set of $G$. For this problem, we describe two NP-completeness proofs, for the classes of claw-free graphs of degree at most 3, and for bounded degree graphs, of maximum degree at most $d \geq 3$ and large girth. In contrast, we prove that the problem admits an $O(n)$ time solution, for cubic claw-free graphs. In addition, we prove a complexity dichotomy theorem for the perfect edge domination problem, based on the results described in the paper. Finally, we describe a linear time algorithm for finding a minimum weight perfect edge dominating set of a $P_5$-free graph. The algorithm is robust, in the sense that, given an arbitrary graph $G$, either it computes a minimum weight perfect edge dominating set of $G$, or it exhibits an induced subgraph of $G$, isomorphic to a $P_5$.


title: An evolutionary strategy for DeltaE - E identification

abstract: In this article we present an automatic method for charge and mass identification of charged nuclear fragments produced in heavy ion collisions at intermediate energies. The algorithm combines a generative model of DeltaE - E relation and a Covariance Matrix Adaptation Evolutionary Strategy (CMA-ES). The CMA-ES is a stochastic and derivative-free method employed to search parameter space of the model by means of a fitness function. The article describes details of the method along with results of an application on simulated labeled data.


title: A finite element implementation of the isotropic exponentiated Hencky-logarithmic model and simulation of the eversion of elastic tubes

abstract: We investigate a finite element formulation of the exponentiated Hencky-logarithmic model whose strain energy function is given by \[ W_\mathrm{eH}(\boldsymbol{F}) = \dfrac{\mu}{k}\, e^{\displaystyle k \left\lVert\mbox{dev}_n \log\boldsymbol{U}\right\rVert^2} + \dfrac{\kappa}{2 \hat{k}}\, e^{\displaystyle \hat{k} [\mbox{tr} (\log\boldsymbol{U})]^2 }\,, \] where $\mu>0$ is the (infinitesimal) shear modulus, $\kappa>0$ is the (infinitesimal) bulk modulus, $k$ and $\hat{k}$ are additional dimensionless material parameters, $\boldsymbol{U}=\sqrt{\boldsymbol{F}^T\boldsymbol{F}}$ and $\boldsymbol{V}=\sqrt{\boldsymbol{F}\boldsymbol{F}^T}$ are the right and left stretch tensor corresponding to the deformation gradient $\boldsymbol{F}$, $\log$ denotes the principal matrix logarithm on the set of positive definite symmetric matrices, $\mbox{dev}_n \boldsymbol{X} = \boldsymbol{X}-\frac{\mbox{tr} \boldsymbol{X}}{n}\boldsymbol{1}$ and $\lVert \boldsymbol{X} \rVert = \sqrt{\mbox{tr}\boldsymbol{X}^T\boldsymbol{X}}$ are the deviatoric part and the Frobenius matrix norm of an $n\times n$-matrix $\boldsymbol{X}$, respectively, and $\mbox{tr}$ denotes the trace operator. To do so, the equivalent different forms of the constitutive equation are recast in terms of the principal logarithmic stretches by use of the spectral decomposition together with the undergoing properties. We show the capability of our approach with a number of relevant examples, including the challenging "eversion of elastic tubes" problem.


title: Monitoring Information Quality within Web Service Composition and Execution

abstract: The composition of web services is a promising approach enabling flexible and loose integration of business applications. Numerous approaches related to web services composition have been developed usually following three main phases: the service discovery is based on the semantic description of advertised services, i.e. the functionality of the service, meanwhile the service selection is based on non- functional quality dimensions of service, and finally the service composition aims to support an underlying process. Most of those approaches explore techniques of static or dynamic design for an optimal service composition. One important aspect so far is mostly neglected, focusing on the output produced of composite web services. In this paper, in contrast to many prominent approaches we introduce a data quality perspective on web services. Based on a data quality management approach, we propose a framework for analyzing data produced by the composite service execution. Utilising process information together with data in service logs, our approach allows identifying problems in service composition and execution. Analyzing the service execution history our approach helps to improve common approaches of service selection and composition.


title: Valley and spin splittings in PbSe nanowires

abstract: We use an empirical tight-binding approach to calculate electron and hole states in [111]-grown PbSe nanowires. We show that the valley-orbit and spin-orbit splittings are very sensitive to the atomic arrangement within the nanowire elementary cell and differ for [111]-nanowires with microscopic $D_{3d}$, $C_{2h}$ and $D_{3}$ symmetries. For the nanowire diameter below 4 nm the valley-orbit splittings become comparable with the confinement energies and the $\boldsymbol{k}\cdot\boldsymbol{p}$ method is inapplicable. Nanowires with the $D_{3}$ point symmetry having no inversion center exhibit giant spin splitting $E = \alpha k_z$, linear in one-dimensional wave vector $k_z$, with the constant $\alpha$ up to 1 eV$\cdot$\AA.


title: A Nitsche Method for Elliptic Problems on Composite Surfaces

abstract: We develop a finite element method for elliptic partial differential equations on so called composite surfaces that are built up out of a finite number of surfaces with boundaries that fit together nicely in the sense that the intersection between any two surfaces in the composite surface is either empty, a point, or a curve segment, called an interface curve. Note that several surfaces can intersect along the same interface curve. On the composite surface we consider a broken finite element space which consists of a continuous finite element space at each subsurface without continuity requirements across the interface curves. We derive a Nitsche type formulation in this general setting and by assuming only that a certain inverse inequality and an approximation property hold we can derive stability and error estimates in the case when the geometry is exactly represented. We discuss several different realizations, including so called cut meshes, of the method. Finally, we present numerical examples.


title: Noise and diffusion of a vibrated self-propelled granular particle

abstract: Granular materials are an important physical realization of active matter. In vibration-fluidized granular matter, both diffusion and self-propulsion derive from the same collisional forcing, unlike many other active systems where there is a clean separation between the origin of single-particle mobility and the coupling to noise. Here we present experimental studies of single-particle motion in a vibrated granular monolayer, along with theoretical analysis that compares grain motion at short and long time scales to the assumptions and predictions, respectively, of the active Brownian particle (ABP) model. The results demonstrate that despite the unique relation between noise and propulsion, granular media do show the generic features predicted by the ABP model and indicate that this is a valid framework to predict collective phenomena. Additionally, our scheme of analysis for validating the inputs and outputs of the model can be applied to other granular and non-granular systems.


title: Better Text Understanding Through Image-To-Text Transfer

abstract: Generic text embeddings are successfully used in a variety of tasks. However, they are often learnt by capturing the co-occurrence structure from pure text corpora, resulting in limitations of their ability to generalize. In this paper, we explore models that incorporate visual information into the text representation. Based on comprehensive ablation studies, we propose a conceptually simple, yet well performing architecture. It outperforms previous multimodal approaches on a set of well established benchmarks. We also improve the state-of-the-art results for image-related text datasets, using orders of magnitude less data.


title: The Brezis-Nirenberg Result for the Fractional Elliptic Problem with Singular Potential

abstract: In this paper, we are concerned with the following type of fractional problems: $$ \begin{cases}\dis (-\Delta)^{s} u-\mu\frac{u}{|x|^{2s}}-\lambda u=|u|^{2^*_{s}-2}u+f(x,u), &\text{in} \Omega,\ \ \, u=0\,&\text{in} \R^N\backslash\Omega \end{cases} \eqno {(*)} $$ where $s\in (0,1)$, $2^*_{s}=2N/(N-2s)$ is the critical Sobolev exponent, $f(x,u)$ is a lower order perturbation of critical Sobolev nonlinearity. We obtain the existence of the solution for (*) through variational methods. In particular we derive a Br\'ezis-Nirenberg type result when $f(x,u)=0$.


title: An X-ray/SDSS sample (II): outflowing gas plasma properties

abstract: Galaxy-scale outflows are nowadays observed in many active galactic nuclei (AGNs); however, their characterisation in terms of (multi-) phase nature, amount of flowing material, effects on the host galaxy, is still unsettled. In particular, ionized gas mass outflow rate and related energetics are still affected by many sources of uncertainties. In this respect, outflowing gas plasma conditions, being largely unknown, play a crucial role. Taking advantage of the spectroscopic analysis results we obtained studying the X-ray/SDSS sample of 563 AGNs at z $<0.8$ presented in our companion paper, we analyse stacked spectra and sub-samples of sources with high signal-to-noise temperature- and density-sensitive emission lines to derive the plasma properties of the outflowing ionized gas component. For these sources, we also study in detail various diagnostic diagrams to infer information about outflowing gas ionization mechanisms. We derive, for the first time, median values for electron temperature and density of outflowing gas from medium-size samples ($\sim 30$ targets) and stacked spectra of AGNs. Evidences of shock excitation are found for outflowing gas. We measure electron temperatures of the order of $\sim 1.7\times10^4$ K and densities of $\sim 1200$ cm$^{-3}$ for faint and moderately luminous AGNs (intrinsic X-ray luminosity $40.5<log(L_X)<44$ in the 2-10 keV band). We caution that the usually assumed electron density ($N_e=100$ cm$^{-3}$) in ejected material might result in relevant overestimates of flow mass rates and energetics and, as a consequence, of the effects of AGN-driven outflows on the host galaxy.


title: A Derandomized Algorithm for RP-ADMM with Symmetric Gauss-Seidel Method

abstract: For multi-block alternating direction method of multipliers(ADMM), where the objective function can be decomposed into multiple block components, we show that with block symmetric Gauss-Seidel iteration, the algorithm will converge quickly. The method will apply a block symmetric Gauss-Seidel iteration in the primal update and a linear correction that can be derived in view of Richard iteration. We also establish the linear convergence rate for linear systems.


title: Teleportation for septuagenarians

abstract: An introduction to the theory of teleportation.


title: Exponential error rates of SDP for block models: Beyond Grothendieck's inequality

abstract: In this paper we consider the cluster estimation problem under the Stochastic Block Model. We show that the semidefinite programming (SDP) formulation for this problem achieves an error rate that decays exponentially in the signal-to-noise ratio. The error bound implies weak recovery in the sparse graph regime with bounded expected degrees, as well as exact recovery in the dense regime. An immediate corollary of our results yields error bounds under the Censored Block Model. Moreover, these error bounds are robust, continuing to hold under heterogeneous edge probabilities and a form of the so-called monotone attack. Significantly, this error rate is achieved by the SDP solution itself without any further pre- or post-processing, and improves upon existing polynomially-decaying error bounds proved using the Grothendieck\textquoteright s inequality. Our analysis has two key ingredients: (i) showing that the graph has a well-behaved spectrum, even in the sparse regime, after discounting an exponentially small number of edges, and (ii) an order-statistics argument that governs the final error rate. Both arguments highlight the implicit regularization effect of the SDP formulation.

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


title: A Decidable Confluence Test for Cognitive Models in ACT-R

abstract: Computational cognitive modeling investigates human cognition by building detailed computational models for cognitive processes. Adaptive Control of Thought - Rational (ACT-R) is a rule-based cognitive architecture that offers a widely employed framework to build such models. There is a sound and complete embedding of ACT-R in Constraint Handling Rules (CHR). Therefore analysis techniques from CHR can be used to reason about computational properties of ACT-R models. For example, confluence is the property that a program yields the same result for the same input regardless of the rules that are applied. In ACT-R models, there are often cognitive processes that should always yield the same result while others e.g. implement strategies to solve a problem that could yield different results. In this paper, a decidable confluence criterion for ACT-R is presented. It allows to identify ACT-R rules that are not confluent. Thereby, the modeler can check if his model has the desired behavior. The sound and complete translation of ACT-R to CHR from prior work is used to come up with a suitable invariant-based confluence criterion from the CHR literature. Proper invariants for translated ACT-R models are identified and proven to be decidable. The presented method coincides with confluence of the original ACT-R models.


title: Unifying and Generalizing Methods for Removing Unwanted Variation Based on Negative Controls

abstract: Unwanted variation, including hidden confounding, is a well-known problem in many fields, particularly large-scale gene expression studies. Recent proposals to use control genes --- genes assumed to be unassociated with the covariates of interest --- have led to new methods to deal with this problem. Going by the moniker Removing Unwanted Variation (RUV), there are many versions --- RUV1, RUV2, RUV4, RUVinv, RUVrinv, RUVfun. In this paper, we introduce a general framework, RUV*, that both unites and generalizes these approaches. This unifying framework helps clarify connections between existing methods. In particular we provide conditions under which RUV2 and RUV4 are equivalent. The RUV* framework also preserves an advantage of RUV approaches --- their modularity --- which facilitates the development of novel methods based on existing matrix imputation algorithms. We illustrate this by implementing RUVB, a version of RUV* based on Bayesian factor analysis. In realistic simulations based on real data we found that RUVB is competitive with existing methods in terms of both power and calibration, although we also highlight the challenges of providing consistently reliable calibration among data sets.


title: Information Theoretic Principles of Universal Discrete Denoising

abstract: Today, the internet makes tremendous amounts of data widely available. Often, the same information is behind multiple different available data sets. This lends growing importance to latent variable models that try to learn the hidden information from the available imperfect versions. For example, social media platforms can contain an abundance of pictures of the same person or object, yet all of which are taken from different perspectives. In a simplified scenario, one may consider pictures taken from the same perspective, which are distorted by noise. This latter application allows for a rigorous mathematical treatment, which is the content of this contribution. We apply a recently developed method of dependent component analysis to image denoising when multiple distorted copies of one and the same image are available, each being corrupted by a different and unknown noise process. In a simplified scenario, we assume that the distorted image is corrupted by noise that acts independently on each pixel. We answer completely the question of how to perform optimal denoising, when at least three distorted copies are available: First we define optimality of an algorithm in the presented scenario, and then we describe an aymptotically optimal universal discrete denoising algorithm (UDDA). In the case of binary data and binary symmetric noise, we develop a simplified variant of the algorithm, dubbed BUDDA, which we prove to attain universal denoising uniformly.