title: Efficient Algorithms for Smooth Minimax Optimization

abstract: This paper studies first order methods for solving smooth minimax optimization problems $\min_x \max_y g(x,y)$ where $g(\cdot,\cdot)$ is smooth and $g(x,\cdot)$ is concave for each $x$. In terms of $g(\cdot,y)$, we consider two settings -- strongly convex and nonconvex -- and improve upon the best known rates in both. For strongly-convex $g(\cdot, y),\ \forall y$, we propose a new algorithm combining Mirror-Prox and Nesterov's AGD, and show that it can find global optimum in $\tilde{O}(1/k^2)$ iterations, improving over current state-of-the-art rate of $O(1/k)$. We use this result along with an inexact proximal point method to provide $\tilde{O}(1/k^{1/3})$ rate for finding stationary points in the nonconvex setting where $g(\cdot, y)$ can be nonconvex. This improves over current best-known rate of $O(1/k^{1/5})$. Finally, we instantiate our result for finite nonconvex minimax problems, i.e., $\min_x \max_{1\leq i\leq m} f_i(x)$, with nonconvex $f_i(\cdot)$, to obtain convergence rate of $O(m(\log m)^{3/2}/k^{1/3})$ total gradient evaluations for finding a stationary point.


title: On the structure of modules over walled Brauer algebra via normal form and random walks

abstract: We analyze cyclic cell modules over walled Brauer algebra in terms of a certain normal form. The latter allows us to decompose the algebra into the generating set and annihilator ideal of a certain cyclic vector. In addition, we show that the numbers of reduced basis monomials of given length coincide with those for the symmetric group. For the semisimple case we utilize the theory of differential posets to calculate the dimensions of modules in terms of the paths in Bratelli diagram. It turns out that the number of primitive idempotents is the same as for the symmetric group.


title: Future High Energy Frontier Colliders

abstract: Colliders have been at the forefront of scientic discoveries in high-energy particle physics since the inception of the colliding beams method in the middle of the 20th century. The field of accelerators is very dynamic and many innovative concepts are currently being considered such future facilities as Higgs factories and energy frontier colliders beyond the LHC. Here we briefly overview leading proposals and studies towards the next generation colliders and discuss their major challenges as well as directions of corresponding accelerator R and D programs needed to address their cost and performance risks.


title: Constraints on anisotropic RG flows from holographic entanglement entropy

abstract: In the context of the gauge/gravity duality, using the proposed candidate $c$-function, which is derived from the entanglement entropy of a strip-shaped region, we investigate the RG flow for $d+1$-dimensional quantum field theories with broken Lorentz and rotational symmetries in the IR, but preserved conformal invariance in the UV boundary. We examine conditions of monotonicity of the $c$-function for holographic anisotropic theories dual to the Einstein gravity via the constraints imposed by the null energy conditions. We consider near UV and IR behaviors and identify the sufficient conditions that guarantee the $c$ function decreases monotonically along the RG flows.


title: Learning algebraic decompositions using Prony structures

abstract: We propose an algebraic framework generalizing several variants of Prony's method and explaining their relations. This includes Hankel and Toeplitz variants of Prony's method for the decomposition of multivariate exponential sums, polynomials (with respect to the monomial and Chebyshev bases), Gau{\ss}ian sums, spherical harmonic sums, taking also into account whether they have their support on an algebraic set.


title: Influence of scattering on optical response of superconductivity

abstract: By using the gauge-invariant kinetic equation, we analytically investigate the influence of the scattering on the optical properties of superconductors in the normal-skin-effect region. Both linear and second-order responses are studied under a multi-cycle terahertz pulse. In the linear regime, we reveal that the optical absorption $\sigma_{1s}(\omega)$, which origins from the scattering, exhibits a crossover point at $\omega=2|\Delta|$. Particularly, it is further shown that when $\omega<2|\Delta|$, $\sigma_{1s}(\omega)$ from the scattering always exhibits a finite value even at low temperature, in contrast to the vanishing $\sigma_{1s}(\omega)$ in the anomalous-skin-effect region as the Mattis-Bardeen theory revealed. In the second-order regime, responses of the Higgs mode during and after the optical pulse are studied. During the pulse, we show that the scattering causes a phase shift in the second-order response of the Higgs mode. Particularly, this phase shift exhibits a significant $\pi$-jump at $\omega=|\Delta|$, which provides a very clear feature for the experimental detection. After the pulse, by studying the damping of the Higgs-mode excitation, we reveal a relaxation mechanism from the elastic scattering, which shows a monotonic enhancement with the increase of the impurity density.


title: Learning to Rank Broad and Narrow Queries in E-Commerce

abstract: Search is a prominent channel for discovering products on an e-commerce platform. Ranking products retrieved from search becomes crucial to address customer's need and optimize for business metrics. While learning to Rank (LETOR) models have been extensively studied and have demonstrated efficacy in the context of web search; it is a relatively new research area to be explored in the e-commerce. In this paper, we present a framework for building LETOR model for an e-commerce platform. We analyze user queries and propose a mechanism to segment queries between broad and narrow based on user's intent. We discuss different types of features - query, product and query-product and discuss challenges in using them. We show that sparsity in product features can be tackled through a denoising auto-encoder while skip-gram based word embeddings help solve the query-product sparsity issues. We also present various target metrics that can be employed for evaluating search results and compare their robustness. Further, we build and compare performances of both pointwise and pairwise LETOR models on fashion category data set. We also build and compare distinct models for broad and narrow queries, analyze feature importance across these and show that these specialized models perform better than a combined model in the fashion world.


title: A filtering problem with uncertainty in observation

abstract: This paper is concerned with a generalized Kalman-Bucy filtering model and corresponding robust problem under model uncertainty. We find that this robust problem is equivalent to considering an estimate problem under some sublinear operator. Therefore, we turn to obtaining the minimum mean square estimator under a sublinear operator. By Girsanov theorem and minimax theorem, we obtain the optimal estimator $\hat{x}_{t}$ of the signal process $x_{t}$ for given time $t\in\lbrack0,T]$.


title: Adaptive particle-based approximations of the Gibbs posterior for inverse problems

abstract: In this work, we adopt a general framework based on the Gibbs posterior to update belief distributions for inverse problems governed by partial differential equations (PDEs). The Gibbs posterior formulation is a generalization of standard Bayesian inference that only relies on a loss function connecting the unknown parameters to the data. It is particularly useful when the true data generating mechanism (or noise distribution) is unknown or difficult to specify. The Gibbs posterior coincides with Bayesian updating when a true likelihood function is known and the loss function corresponds to the negative log-likelihood, yet provides subjective inference in more general settings. We employ a sequential Monte Carlo (SMC) approach to approximate the Gibbs posterior using particles. To manage the computational cost of propagating increasing numbers of particles through the loss function, we employ a recently developed local reduced basis method to build an efficient surrogate loss function that is used in the Gibbs update formula in place of the true loss. We derive error bounds for our approximation and propose an adaptive approach to construct the surrogate model in an efficient manner. We demonstrate the efficiency of our approach through several numerical examples.


title: Forecasting high-dimensional dynamics exploiting suboptimal embeddings

abstract: Delay embedding---a method for reconstructing dynamical systems by delay coordinates---is widely used to forecast nonlinear time series as a model-free approach. When multivariate time series are observed, several existing frameworks can be applied to yield a single forecast combining multiple forecasts derived from various embeddings. However, the performance of these frameworks is not always satisfactory because they randomly select embeddings or use brute force and do not consider the diversity of the embeddings to combine. Herein, we develop a forecasting framework that overcomes these existing problems. The framework exploits various "suboptimal embeddings" obtained by minimizing the in-sample error via combinatorial optimization. The framework achieves the best results among existing frameworks for sample toy datasets and a real-world flood dataset. We show that the framework is applicable to a wide range of data lengths and dimensions. Therefore, the framework can be applied to various fields such as neuroscience, ecology, finance, fluid dynamics, weather, and disaster prevention.


title: Security analysis of practical continuous-variable quantum key distribution systems under laser seeding attack

abstract: Here, we investigate the security of the practical one-way CVQKD and CV-MDI-QKD systems under laser seeding attack. In particular, Eve can inject a suitable light into the laser diodes of the light source modules in the two kinds of practical CVQKD systems, which results in the increased intensity of the generated optical signal. The parameter estimation under the attack shows that the secret key rates of these two schemes may be overestimated, which opens a security loophole for Eve to successfully perform an intercept-resend attack on these systems. To close this loophole, we propose a real-time monitoring scheme to precisely evaluate the secret key rates of these schemes. The analysis results indicate the implementation of the proposed monitoring scheme can effectively resist this potential attack.


title: Global existence in critical spaces for non Newtonian compressible viscoelastic flows

abstract: We are interested in the multi-dimentional compressible viscoelastic flows of Oldroyd type, which is one of non-Newtonian fluids exhibiting the elastic behavior. In order to capture the damping effect of the additional deformation tensor, to the best of our knowledge, the "div-curl" structural condition plays a key role in previous efforts. Our aim of this paper is to remove the structural condition and prove a global existence of strong solutions to compressible viscoelastic flows in critical spaces. The new ingredient lies in the introduction of effective flux $(\theta,\mathcal{G})$, which enables us to capture the dissipation arising from \textit{combination} of density and deformation tensor. In absence of compatible conditions, the partial dissipation is found in non-Newtonian compressible fluids, which is weaker than that of usual Navier-Stokes equations.


title: An order-disorder phase transition in black-hole star clusters

abstract: The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.


title: Three-pion contribution to hadronic vacuum polarization

abstract: We address the contribution of the $3\pi$ channel to hadronic vacuum polarization (HVP) using a dispersive representation of the $e^+e^-\to 3\pi$ amplitude. This channel gives the second-largest individual contribution to the total HVP integral in the anomalous magnetic moment of the muon $(g-2)_\mu$, both to its absolute value and uncertainty. It is largely dominated by the narrow resonances $\omega$ and $\phi$, but not to the extent that the off-peak regions were negligible, so that at the level of accuracy relevant for $(g-2)_\mu$ an analysis of the available data as model independent as possible becomes critical. Here, we provide such an analysis based on a global fit function using analyticity and unitarity of the underlying $\gamma^*\to3\pi$ amplitude and its normalization from a chiral low-energy theorem, which, in particular, allows us to check the internal consistency of the various $e^+e^-\to 3\pi$ data sets. Overall, we obtain $a_\mu^{3\pi}|_{\leq 1.8\,\text{GeV}}=46.2(6)(6)\times 10^{-10}$ as our best estimate for the total $3\pi$ contribution consistent with all (low-energy) constraints from QCD. In combination with a recent dispersive analysis imposing the same constraints on the $2\pi$ channel below $1\,\text{GeV}$, this covers nearly $80\%$ of the total HVP contribution, leading to $a_\mu^\text{HVP}=692.3(3.3)\times 10^{-10}$ when the remainder is taken from the literature, and thus reaffirming the $(g-2)_\mu$ anomaly at the level of at least $3.4\sigma$. As side products, we find for the vacuum-polarization-subtracted masses $M_\omega=782.63(3)(1)\,\text{MeV}$ and $M_\phi=1019.20(2)(1)\,\text{MeV}$, confirming the tension to the $\omega$ mass as extracted from the $2\pi$ channel.


title: Neutron star mergers and rare core-collapse supernovae as sources of r-process enrichment in simulated galaxies

abstract: We use cosmological, magnetohydrodynamical simulations of Milky Way-mass galaxies from the Auriga project to study their enrichment with rapid neutron capture (r-process) elements. We implement a variety of enrichment models from both binary neutron star mergers and rare core-collapse supernovae. We focus on the abundances of (extremely) metal-poor stars, most of which were formed during the first ~Gyr of the Universe in external galaxies and later accreted onto the main galaxy. We find that the majority of metal-poor stars are r-process enriched in all our enrichment models. Neutron star merger models result in a median r-process abundance ratio which increases with metallicity, whereas the median trend in rare core-collapse supernova models is approximately flat. The scatter in r-process abundance increases for models with longer delay times or lower rates of r-process producing events. Our results are nearly perfectly converged, in part due to the mixing of gas between mesh cells in the simulations. Additionally, different Milky Way-mass galaxies show only small variation in their respective r-process abundance ratios. Current (sparse and potentially biased) observations of metal-poor stars in the Milky Way seem to prefer rare core-collapse supernovae over neutron star mergers as the dominant source of r-process elements at low metallicity, but we discuss possible caveats to our models. Dwarf galaxies which experience a single r-process event early in their history show highly enhanced r-process abundances at low metallicity, which is seen both in observations and in our simulations. We also find that the elements produced in a single event are mixed with ~10^8 Msun of gas relatively quickly, distributing the r-process elements over a large region.


title: Topological states on fractal lattices

abstract: We investigate the fate of topological states on fractal lattices. Focusing on a spinless chiral p-wave paired superconductor, we find that this model supports two qualitatively distinct phases when defined on a Sierpinski gasket. While the trivial phase is characterized by a self-similar spectrum with infinitely many gaps and extended eigenstates, the novel "topological" phase has a gapless spectrum and hosts chiral states propagating along edges of the graph. Besides employing theoretical probes such as the real-space Chern number, inverse participation ratio, and energy-level statistics in the presence of disorder, we develop a simple physical picture capturing the essential features of the model on the gasket. Extending this picture to other fractal lattices and topological states, we show that the p+ip state admits a gapped topological phase on the Sierpinski carpet and that a higher-order topological insulator placed on this lattice hosts gapless modes localized on corners.


title: Observational Properties of Ultra-Diffuse Galaxies in the Field: Field-UDGs are Predominantly Blue and Starforming

abstract: While we have learned much about Ultra-Diffuse Galaxies (UDGs) in groups and clusters, relatively little is known about them in less-dense environments. More isolated UDGs are important for our understanding of UDG formation scenarios because they form via secular mechanisms, allowing us to determine the relative importance of environmentally-driven formation in groups and clusters. We have used the public Kilo-Degree Survey (KiDS) together with the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the abundance and properties of UDGs in the field, targeting sources with low surface brightness (24.0$\leq$\bar{\mu}_{e,r}}$\leq$\26.5) and large apparent sizes (3.0\arcsec$\leq$\bar{r}_{e,r}}$\leq$8.0\arcsec). Accounting for several sources of interlopers in our selection based on canonical scaling relations, and using an empirical UDG model based on measurements from the literature, we show that a scenario in which cluster-like red sequence UDGs occupy a significant number of field galaxies is unlikely, with most field UDGs being significantly bluer and showing signs of localised star formation. An immediate conclusion is that UDGs are much more efficiently quenched in high-density environments. We estimate an upper-limit on the total field abundance of UDGs of 8$\pm$3$\times10^{-3}$cMpc$^{-3}$ within our selection range. We also compare the total field abundance of UDGs to a measurement of the abundance of HI-rich UDGs from the literature, suggesting that they occupy at least one-fifth of the overall UDG population. The mass formation efficiency of UDGs implied by this upper-limit is similar to what is measured in groups and clusters.


title: The Fundamental Relation between Halo Mass and Galaxy Group Properties

abstract: We explore the interrelationships between the galaxy group halo mass and various observable group properties. We propose a simple scenario that describes the evolution of the central galaxies and their host dark matter halos. Star formation quenching is one key process in this scenario, which leads to the different assembly histories of blue groups (group with a blue central) and red groups (group with a red central). For blue groups, both the central galaxy and the halo continue to grow their mass. For red groups, the central galaxy has been quenched and its stellar mass remains about constant, while its halo continues to grow by merging smaller halos. From this simple scenario, we speculate about the driving properties that should strongly correlate with the group halo mass. We then apply the machine learning algorithm the Random Forest (RF) regressor to blue groups and red groups separately in the semianalytical model L-GALAXIES to explore these nonlinear multicorrelations and to verify the scenario as proposed above. Remarkably, the results given by the RF regressor are fully consistent with the prediction from our simple scenario and hence provide strong support for it. As a consequence, the group halo mass can be more accurately determined from observable galaxy properties by the RF regressor with a 50% reduction in error. A halo mass more accurately determined in this way also enables more accurate investigations on the galaxy-halo connection and other important related issues, including galactic conformity and the effect of halo assembly bias on galaxy assembly.


title: Distinguishing black holes from horizonless objects through the excitation of resonances during inspiral

abstract: How well is the vacuum Kerr geometry a good description of the dark, compact objects in our universe? Precision measurements of accreting matter in the deep infrared and gravitational-wave measurements of coalescing objects are finally providing answers to this question. Here, we study the possibility of resonant excitation of the modes of the central object -- taken to be very compact but horizonless -- during an extreme-mass-ratio inspiral. We show that for very compact objects resonances are indeed excited. However, the impact of such excitation on the phase of the gravitational-wave signal is negligible, since resonances are crossed very quickly during inspiral.


title: Novel color superconducting phases of $\cal{N}$ = 4 super Yang-Mills at strong coupling

abstract: We revisit the large-$N_c$ phase diagram of $\cal{N}$ = 4 super Yang-Mills theory at finite $R$-charge density and strong coupling, by means of the AdS/CFT correspondence. We conjecture new phases that result from a black hole shedding some of its charge through the nucleation of probe color D3-branes that remain at a finite distance from the black hole when the dual field theory lives on a sphere. In the corresponding ground states the color group is partially Higgsed, so these phases can be identified as having a type of color superconductivity. The new phases would appear at intermediate values of the $R$-charge chemical potential and we expect them to be metastable but long-lived in the large-$N_c$ limit.


title: The Dependence of AGN Activity on Environment in SDSS

abstract: Environment is one of the key external drivers of the galaxies, while active galactic nucleus (AGN) is one of the key internal drivers. Both of them play fundamental roles in regulating the formation and evolution of galaxies. We explore the interrelationship between environment and AGN in SDSS. At a given stellar mass, the specific star formation rate distribution of the AGN host galaxies remains unchanged with over-density, with the peak of the distribution around the Green Valley. We show that, at a given stellar mass, the AGN fraction that has been commonly used in previous studies (defined as the number of AGNs relative to all galaxies including passive and star forming ones) does decrease with increasing over-density for satellites. This is largely due to the fact that the fraction of passive galaxies strongly depends on environment. In order to investigate the intrinsic correlation between AGN and environment, especially under the assumption that AGN feedback is responsible for star formation quenching, the AGN fraction should be defined as the number of AGNs relative to the star-forming galaxies only. With the new definition, we find little dependence of AGN fraction on over-density, central/satellite, and group halo mass. There is only marginal evidence that AGN may prefer denser regions, which is possibly due to more frequent interaction of galaxies or higher merger rate in groups. Our results support the scenario that internal secular evolution is the predominant mechanism of triggering AGN activity, while external environment related processes only play a minor role.


title: The radial acceleration relation and dark baryons in MOND

abstract: Recent observations of rotationally supported galaxies show a tight correlation between the observed radial acceleration at every radius and the Newtonian acceleration generated by the baryonic mass distribution, the so-called radial acceleration relation (RAR). The rotation curves (RCs) of the SPARC sample of disk galaxies with different morphologies, masses, sizes and gas fractions are investigated in the context of modified Newtonian dynamics (MOND). We include the effect of cold dark baryons by scaling the measured mass in the atomic form by a factor of $c$ in the mass budget of galaxies. In addition to the standard interpolating function, we also fit the RCs and the RAR with the empirical RAR-inspired interpolating function. Slightly better fits for about $47\%$ of galaxies in our sample are achieved in the presence of dark baryons ($c>1$) with the mean value of $c = 2.4\pm 1.3$. Although the MOND fits are not significantly improved by including dark baryons, it results in a decrease in the characteristic acceleration $g_\dag$ by $40\%$. We find no correlation between the MOND critical acceleration $a_0$ and the central surface brightness of the stellar disk, $\mu_{3.6}$. This supports $a_0$ being a universal constant for all galaxies.


title: Bayesian hierarchical inference of asteroseismic inclination angles

abstract: The stellar inclination angle-the angle between the rotation axis of a star and our line of sight-provides valuable information in many different areas, from the characterisation of the geometry of exoplanetary and eclipsing binary systems, to the formation and evolution of those systems. We propose a method based on asteroseismology and a Bayesian hierarchical scheme for extracting the inclination angle of a single star. This hierarchical method therefore provides a means to both accurately and robustly extract inclination angles from red giant stars. We successfully apply this technique to an artificial dataset with an underlying isotropic inclination angle distribution to verify the method. We also apply this technique to 123 red giant stars observed with $\textit{Kepler}$. We also show the need for a selection function to account for possible population-level biases, that are not present in individual star-by-star cases, in order to extend the hierarchical method towards inferring underlying population inclination angle distributions.


title: Temperature and density dependent cooling function for H$_2$ with updated H$_2$/H collisional rates

abstract: The energy transfer among the components in a gas determines its fate. Especially at low temperatures, inelastic collisions drive the cooling and the heating mechanisms. In the early Universe as well as in zero- or low- metallicity environments the major contribution comes from the collisions among atomic and molecular hydrogen, also in its deuterated version. The present work shows some updated calculations of the H$_2$ cooling function based on novel collisional data which explicitely take into account the reactive pathway at low temperatures. Deviations from previous calculations are discussed and a multivariate data analysis is performed to provide a fit depending on both the gas temperature and the density of the gas.


title: The Southern Photometric Local Universe Survey (S-PLUS): improved SEDs, morphologies and redshifts with 12 optical filters

abstract: The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

categories: astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.IM astro-ph.SR


title: Locality & Entanglement in Table-Top Testing of the Quantum Nature of Linearized Gravity

abstract: This paper points out the importance of the assumption of locality of physical interactions, and the concomitant necessity of propagation of an entity (in this case, off-shell quanta - virtual gravitons) between two non-relativistic test masses in unveiling the quantum nature of linearized gravity through a laboratory experiment. At the outset, we will argue that observing the quantum nature of a system is not limited to evidencing $O\left(\hbar\right)$ corrections to a classical theory: it instead hinges upon verifying tasks that a classical system cannot accomplish. We explain the background concepts needed from quantum field theory and quantum information theory to fully appreciate the previously proposed table-top experiments: namely forces arising through the exchange of virtual (off-shell) quanta, as well as Local Operations and Classical Communication (LOCC) and entanglement witnesses. We clarify the key assumption inherent in our evidencing experiment, namely the locality of physical interactions, which is a generic feature of interacting systems of quantum fields around us, and naturally incorporates micro-causality in the description of our experiment. We also present the types of states the matter field must inhabit, putting the experiment on firm relativistic quantum field theoretic grounds. At the end we use a non-local theory of gravity to illustrate how our mechanism may still be used to detect the qualitatively quantum nature of a force when the scale of non-locality is finite. We find that the scale of non-locality, including the entanglement entropy production in local/ non-local gravity, may be revealed from the results of our experiment.


title: Structure of Quantum Entanglement at a Finite Temperature Critical Point

abstract: We propose a scheme to characterize long-range quantum entanglement close to a finite temperature critical point using tripartite entanglement negativity. As an application, we study a model with mean-field Ising critical exponents and find that the tripartite negativity does not exhibit any singularity in the thermodynamic limit across the transition. This indicates that the long-distance critical fluctuations are completely classical, allowing one to define a `quantum correlation length' that remains finite at the transition despite a divergent physical correlation length. Motivated by our model, we also study mixed state entanglement in tight-binding models of bosons with $U(1)$ and time-reversal symmetry. By employing Glauber-Sudarshan `P-representation', we find a surprising result that such states have zero entanglement.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: Using machine learning to understand physics graduate school admissions

abstract: Among all of the first-year graduate students enrolled in doctoral-granting physics departments, the percentage of female and racial minority students has remained unchanged for the past 20 years. The current graduate program admissions process can create challenges for achieving diversity goals in physics. In this paper, we will investigate how the various aspects of a prospective student's application to a physics doctoral program affect the likelihood the applicant will be admitted. Admissions data was collected from a large, Midwestern public research university that has a decentralized admissions process and included applicants' undergraduate GPAs and institutions, research interests, and GRE scores. Because the collected data varied in scale, we used supervised machine learning algorithms to create models that predict who was admitted into the PhD program. We find that using only the applicant's undergraduate GPA and physics GRE score, we are able to predict with 75% accuracy who will be admitted to the program.


title: Characterization of Sobolev spaces on the sphere

abstract: We prove a characterization of the Sobolev spaces $H^\alpha$ on the unit sphere $\mathbb{S}^{d-1}$, where the smoothness index $\alpha$ is any positive real number and $d\geq 2$. This characterization does not use differentiation and it is given in terms of $([\alpha/2]+1)$-multidimensional square functions $S_\alpha$. For $[\alpha/2]=0,$ a function $f\in L^2(\mathbb{S}^{d-1})$ belongs to $H^\alpha(\mathbb{S}^{d-1})$ if and only if $S_\alpha (f)\in L^2(\mathbb{S}^{d-1})$. If $n=[\alpha/2]>0$, the membership of $f$ is equivalent to the existence of $g_1,\cdots,g_n$ in $L^2(\mathbb{S}^{d-1})$ such that $S_\alpha(f,g_1,\ldots,g_n)\in L^2(\mathbb{S}^{d-1})$ and in this case, $g_j=T_j((-\Delta_S)^j f)$, where $T_j$ is a zonal Fourier multiplier in the sphere and $\Delta_S$ is the Laplace-Beltrami operator. The square functions $S_\alpha $ are based on averaging operators over euclidean balls (caps) in the sphere that may be viewed as zonal multipliers. The results in the paper are in the spirit of the characterization of fractional Sobolev spaces given in $\mathbb{R}^d$ proved in \cite{AMV}. The development of the theory is fully based on zonal Fourier multipliers and special functions.


title: A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity

abstract: A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other discontinuities. This paper investigates discontinuous Galerkin methods for the solution of elliptic problems in numerical relativity. We present a novel hp-adaptive numerical scheme for curvilinear and non-conforming meshes. It uses a multigrid preconditioner with a Chebyshev or Schwarz smoother to create a very scalable discontinuous Galerkin code on generic domains. The code employs compactification to move the outer boundary near spatial infinity. We explore the properties of the code on some test problems, including one mimicking Neutron stars with phase transitions. We also apply it to construct initial data for two or three black holes.


title: Target Flow-Pressure Operating Range for Designing a Failing Fontan Cavopulmonary Support Device

abstract: Fontan operation as the current standard of care for the palliation of single ventricle defects results in significant late complications. Using a mechanical circulatory device for the right circulation to serve the function of the missing subpulmonary ventricle could potentially stabilize the failing Fontan circulation. This study aims to elucidate the hydraulic operating regions that should be targeted for designing cavopulmonary blood pumps. By integrating numerical analysis and available clinical information, the interaction of the cavopulmonary support via the IVC and full assist configurations with a wide range of simulated adult failing scenarios was investigated; with IVC and full assist corresponding to the inferior venous return or the entire venous return, respectively, being routed through the device. We identified the desired hydraulic operating regions for a cavopulmonary assist device by clustering all head pressures and corresponding pump flows that result in hemodynamic improvement for each simulated failing Fontan physiology. Results show that IVC support can produce beneficial hemodynamics in only a small fraction of failing Fontan scenarios. Cavopulmonary assist device could increase cardiac index by 35% and decrease the inferior vena cava pressure by 45% depending on the patient's pre-support hemodynamic state and surgical configuration of the cavopulmonary assist device (IVC or full support). The desired flow-pressure operating regions we identified can serve as the performance criteria for designing cavopulmonary assist devices as well as evaluating off-label use of commercially available left-side blood pumps for failing Fontan cavopulmonary support.


title: The interplay between SM precision, BSM physics and measurements of $\alpha_{\rm{had}}$ in $\mu$-$e$ scattering

abstract: Muon electron scattering experiments such as the proposed MUonE experiment, offer an opportunity for an improved measurement of the Leading Order hadronic running of $\alpha$, denoted $\Delta \alpha_{\rm{had}}$. Such a measurement could be utilized to reduce the theoretical uncertainty on the prediction of the anomalous magnetic moment of the muon, $g-2$. Currently there is a discrepancy between theory and data for this observable which could potentially be explained by Beyond the Standard Model (BSM) physics. Here we investigate the possible impact of missing Standard Model (SM) higher order corrections and BSM physics on the proposed measurement of $\Delta \alpha_{\rm{had}}$. In principle either could be indirectly fitted into $\Delta \alpha_{\rm{had}}$, causing inconsistencies if used in a $g-2$ application. The literature suggests a target of 10 ppm on the cross section for the theoretical accuracy. We assess the validity of this target in detail using a variety of methods, finding that a 1 ppm target is a more conservative estimate to ensure missing higher orders do not dominate the theoretical uncertainty. For the potential BSM contributions we study various models which contribute first at tree- and loop-level. Of particular interest is the impact from dark photon models, which can potentially affect the measurement of $\Delta \alpha_{\rm{had}}$ at the desired accuracy. At loop-level there exists in general a kinematic suppression adequate to reduce the BSM contributions to a level which can be neglected for the extraction of $\Delta \alpha_{\rm{had}}$.


title: Cell-scale degradation of peritumoural extracellular matrix fibre network and its role within tissue-scale cancer invasion

abstract: Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. Occurring over many different temporal and spatial scales, the first stage of invasion is the secretion of matrix degrading enzymes (MDEs) by the cancer cells that consequently degrade the surrounding extracellular matrix (ECM). This process is vital for creating space in which the cancer cells can progress and it is driven by the activities of specific matrix metalloproteinases (MMPs). In this paper, we consider the key role of two MMPs by developing further the novel two-part multiscale model introduced in [33] to better relate at micro-scale the two micro-scale activities that were considered there, namely, the micro-dynamics concerning the continuous rearrangement of the naturally oriented ECM fibres within the bulk of the tumour and MDEs proteolytic micro-dynamics that take place in an appropriate cell-scale neighbourhood of the tumour boundary. Focussing primarily on the activities of the membrane-tethered MT1-MMP and the soluble MMP-2 with the fibrous ECM phase, in this work we investigate the MT1-MMP/MMP-2 cascade and its overall effect on tumour progression. To that end, we will propose a new multiscale modelling framework by considering the degradation of the ECM fibres not only to take place at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour.


title: Election predictions are arbitrage-free: response to Taleb

abstract: Taleb (2018) claimed a novel approach to evaluating the quality of probabilistic election forecasts via no-arbitrage pricing techniques and argued that popular forecasts of the 2016 U.S. Presidential election had violated arbitrage boundaries. We show that under mild assumptions all such political forecasts are arbitrage-free and that the heuristic that Taleb's argument was based on is false.


title: SVM Enhanced Frenet Frame Planner For Safe Navigation Amidst Moving Agents

abstract: This paper proposes an SVM Enhanced Trajectory Planner for dynamic scenes, typically those encountered in on road settings. Frenet frame based trajectory generation is popular in the context of autonomous driving both in research and industry. We incorporate a safety based maximal margin criteria using a SVM layer that generates control points that are maximally separated from all dynamic obstacles in the scene. A kinematically consistent trajectory generator then computes a path through these waypoints. We showcase through simulations as well as real world experiments on a self driving car that the SVM enhanced planner provides for a larger offset with dynamic obstacles than the regular Frenet frame based trajectory generation. Thereby, the authors argue that such a formulation is inherently suited for navigation amongst pedestrians. We assume the availability of an intent or trajectory prediction module that predicts the future trajectories of all dynamic actors in the scene.


title: The Information Processing Factory: Organization, Terminology, and Definitions

abstract: The Information Processing Factory (IPF) project has recently introduced the abstraction of complex architectures as self-aware information processing factories. These factories consist of a set of highly configurable resources, e.g., processing elements and interconnects, whose use is monitored, planned, and configured during runtime. Managing a factory involves multiple facets, such as efficiency, availability, reliability, integrity, and timing. IPF conquers the complexity of managing facets in digital systems by hierarchically decomposing the challenges and addressing them with different co-existing entities in the factory. This paper introduces the organization, terminology, and definitions of IPF.


title: Peculiar Index Relations, 2D TQFT, and Universality of SUSY Enhancement

abstract: We study certain exactly marginal gaugings involving arbitrary numbers of Argyres-Douglas (AD) theories and show that the resulting Schur indices are related to those of certain Lagrangian theories of class $\mathcal{S}$ via simple transformations. By writing these quantities in the language of 2D topological quantum field theory (TQFT), we easily read off the $S$-duality action on the flavor symmetries of the AD quivers and also find expressions for the Schur indices of various classes of exotic AD theories appearing in different decoupling limits. The TQFT expressions for these latter theories are related by simple transformations to the corresponding quantities for certain well-known isolated theories with regular punctures (e.g., the Minahan-Nemeschansky $E_6$ theory and various generalizations). We then reinterpret the TQFT expressions for the indices of our AD theories in terms of the topology of the corresponding 3D mirror quivers, and we show that our isolated AD theories generically admit renormalization group (RG) flows to interacting superconformal field theories (SCFTs) with thirty-two (Poincar\'e plus special) supercharges. Motivated by these examples, we argue that, in a sense we make precise, the existence of RG flows to interacting SCFTs with thirty-two supercharges is generic in a far larger class of 4D $\mathcal{N}=2$ SCFTs arising from compactifications of the 6D $(2,0)$ theory on surfaces with irregular singularities.


title: Dark energy explained by an inadequate fitting of the FLRW metric

abstract: Approximating a real manifold by an idealized one requires to calibrate the parameters characterizing the idealized manifold in function of the real one. This calibration is a purely conventional process and can generally be done in several ways, leading to different fittings. In practice, however, all possible fittings cannot be considered as representative of the real manifold. Approximating the real metric of the universe by the FLRW metric would be adequate only if both corresponding structures, defined by the space-time interval, are equivalent on large scales. This requirement puts some constraints on what would be a representative FLRW metric. We show that the way how measurements on SNIa are interpreted to determine the evolution of the scale factor implicitly define the calibration process, and that this one is compatible with the aforementioned constraints. On a theoretical point of view, this indicates that the as fitted FLRW metric would indeed be representative of the real one. On a practical point of view, however, we show that a bias in the measurements could invalidate this conclusion. The bias comes from the fact that SNIa are not randomly distributed over space, but are probably mostly located in regions were matter is largely present, i.e., in overdense regions. We explain how this bias could account for the apparent accelerated expansion of the universe, without needing to introduce the dark energy assumption. We show in particular that this bias leads to an inadequate fitting of the FLRW metric, resulting in the appearance of a new term in the evolution equation of the related scale factor, being equivalent to the cosmological constant.


title: Covering graphs with convex sets and partitioning graphs into convex sets

abstract: We present some complexity results concerning the problems of covering a graph with $p$ convex sets and of partitioning a graph into $p$ convex sets. The following convexities are considered: digital convexity, monophonic convexity, $P_3$-convexity, and $P_3^*$-convexity.


title: Minimum Power to Maintain a Nonequilibrium Distribution of a Markov Chain

abstract: Biological systems use energy to maintain non-equilibrium distributions for long times, e.g. of chemical concentrations or protein conformations. What are the fundamental limits of the power used to "hold" a stochastic system in a desired distribution over states? We study the setting of an uncontrolled Markov chain $Q$ altered into a controlled chain $P$ having a desired stationary distribution. Thermodynamics considerations lead to an appropriately defined Kullback-Leibler (KL) divergence rate $D(P||Q)$ as the cost of control, a setting introduced by Todorov, corresponding to a Markov decision process with mean log loss action cost. The optimal controlled chain $P^*$ minimizes the KL divergence rate $D(\cdot||Q)$ subject to a stationary distribution constraint, and the minimal KL divergence rate lower bounds the power used. While this optimization problem is familiar from the large deviations literature, we offer a novel interpretation as a minimum "holding cost" and compute the minimizer $P^*$ more explicitly than previously available. We state a version of our results for both discrete- and continuous-time Markov chains, and find nice expressions for the important case of a reversible uncontrolled chain $Q$, for a two-state chain, and for birth-and-death processes.

categories: cond-mat.stat-mech cs.IT math.IT physics.bio-ph


title: Could face-centered cubic titanium in cold-rolled commercially-pure titanium only be a Ti-hydride?

abstract: A face-centered cubic (FCC) phase in electro-polished specimens for transmission electron microscopy of commercially pure titanium has sometimes been reported. Here, a combination of atom-probe tomography, scanning transmission electron microscopy and low-loss electron energy loss spectroscopy is employed to study both the crystal structural and chemical composition of this FCC phase. Our results prove that the FCC phase is actually a TiHx (x>1) hydride, and not a new allotrope of Ti, in agreement with previous reports. The formation of the hydride is discussed.


title: Constructive description of analytic Besov spaces in strictly pseudoconvex domains

abstract: We use the method of pseudoanalytic continuation to obtain a characterization of spaces of holomorphic functions with boundary values in Besov spaces in terms of polynomial approximations.


title: Strongly bent double-stranded DNA: reconciling theory and experiment

abstract: The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the polymer makes the harmonic deformation energetically unfavorable: in this strong bending regime the energy of the polymer varies linearly with the average bending angle as the system follows the convex hull of the deformation energy function. For double-stranded DNA, the effective bending deformation energy becomes non-convex for bends greater than about 2 degrees per base-pair, equivalent to the curvature of a closed circular loop of about 160 base pairs. A simple equation is derived for the polymer loop energy that covers both the weak and strong bending regimes. The theory shows quantitative agreement with recent DNA cyclization experiments on short DNA fragments, while maintaining the expected agreement with experiment in the weak bending regime. Counter-intuitively, cyclization probability (j-factor) of very short DNA loops is predicted to increase with decreasing loop length; the j-factor reaches its minimum for loops of 45 base pairs. Atomistic simulations reveal that the attractive component of the short-range Lennard-Jones interaction between the backbone atoms can explain the underlying non-convexity of the DNA effective bending energy, leading to the linear bending regime. Applicability of the theory to protein-DNA complexes, including the nucleosome, is discussed.


title: Protecting Privacy of Users in Brain-Computer Interface Applications

abstract: Machine learning (ML) is revolutionizing research and industry. Many ML applications rely on the use of large amounts of personal data for training and inference. Among the most intimate exploited data sources is electroencephalogram (EEG) data, a kind of data that is so rich with information that application developers can easily gain knowledge beyond the professed scope from unprotected EEG signals, including passwords, ATM PINs, and other intimate data. The challenge we address is how to engage in meaningful ML with EEG data while protecting the privacy of users. Hence, we propose cryptographic protocols based on Secure Multiparty Computation (SMC) to perform linear regression over EEG signals from many users in a fully privacy-preserving (PP) fashion, i.e.~such that each individual's EEG signals are not revealed to anyone else. To illustrate the potential of our secure framework, we show how it allows estimating the drowsiness of drivers from their EEG signals as would be possible in the unencrypted case, and at a very reasonable computational cost. Our solution is the first application of commodity-based SMC to EEG data, as well as the largest documented experiment of secret sharing based SMC in general, namely with 15 players involved in all the computations.


title: Deuteron-alpha scattering: separable vs nonseparable Faddeev approach

abstract: {\bf Background} Deuteron induced reactions are widely used to probe nuclear structure and astrophysical information. Those (d,p) reactions may be viewed as three-body reactions and described with Faddeev techniques. {\bf Purpose} Faddeev-AGS equations in momentum space have a long tradition of utilizing separable interactions in order to arrive at sets of coupled integral equations in one variable. However, it needs to be demonstrated that observables calculated based on separable interactions agree exactly with those based on nonseparable forces. {\bf Methods} Momentum space AGS equations are solved with separable and nonseparable forces as coupled integral equations. {\bf Results} Deuteron-alpha scattering is calculated via momentum space AGS equations using the CD-Bonn neutron-proton force and a Woods-Saxon type neutron(proton)-$^4$He force, for which the Pauli-forbidden S-wave bound state is projected out. Elastic as well as breakup observables are calculated and compared to results in which the interactions in the two-body sub-systems are represented by separable interactions derived in the Ernst-Shakin-Thaler (EST) framework. {\bf Conclusions} We find that the calculations based on the separable representation of the interactions and the original interactions give results that are in excellent agreement. Specifically, integrated cross sections and angular distributions for elastic scattering agree within $\approx$ 1\%, which is well below typical experimental errors. In addition, the five-fold differential cross sections corresponding to breakup of the deuteron agree extremely well.


title: Reverse engineering neural networks from many partial recordings

abstract: Much of neuroscience aims at reverse engineering the brain, but we only record a small number of neurons at a time. We do not currently know if reverse engineering the brain requires us to simultaneously record most neurons or if multiple recordings from smaller subsets suffice. This is made even more important by the development of novel techniques that allow recording from selected subsets of neurons, e.g. using optical techniques. To get at this question, we analyze a neural network, trained on the MNIST dataset, using only partial recordings and characterize the dependency of the quality of our reverse engineering on the number of simultaneously recorded "neurons". We find that reverse engineering of the nonlinear neural network is meaningfully possible if a sufficiently large number of neurons is simultaneously recorded but that this number can be considerably smaller than the number of neurons. Moreover, recording many times from small random subsets of neurons yields surprisingly good performance. Application in neuroscience suggests to approximate the I/O function of an actual neural system, we need to record from a much larger number of neurons. The kind of scaling analysis we perform here can, and arguably should be used to calibrate approaches that can dramatically scale up the size of recorded data sets in neuroscience.


title: Hyper-Molecules: on the Representation and Recovery of Dynamical Structures, with Application to Flexible Macro-Molecular Structures in Cryo-EM

abstract: Cryo-electron microscopy (cryo-EM), the subject of the 2017 Nobel Prize in Chemistry, is a technology for determining the 3-D structure of macromolecules from many noisy 2-D projections of instances of these macromolecules, whose orientations and positions are unknown. The molecular structures are not rigid objects, but flexible objects involved in dynamical processes. The different conformations are exhibited by different instances of the macromolecule observed in a cryo-EM experiment, each of which is recorded as a particle image. The range of conformations and the conformation of each particle are not known a priori; one of the great promises of cryo-EM is to map this conformation space. Remarkable progress has been made in determining rigid structures from homogeneous samples of molecules in spite of the unknown orientation of each particle image and significant progress has been made in recovering a few distinct states from mixtures of rather distinct conformations, but more complex heterogeneous samples remain a major challenge. We introduce the ``hyper-molecule'' framework for modeling structures across different states of heterogeneous molecules, including continuums of states. The key idea behind this framework is representing heterogeneous macromolecules as high-dimensional objects, with the additional dimensions representing the conformation space. This idea is then refined to model properties such as localized heterogeneity. In addition, we introduce an algorithmic framework for recovering such maps of heterogeneous objects from experimental data using a Bayesian formulation of the problem and Markov chain Monte Carlo (MCMC) algorithms to address the computational challenges in recovering these high dimensional hyper-molecules. We demonstrate these ideas in a prototype applied to synthetic data.


title: Measurement of the top quark Yukawa coupling from $\mathrm{t\bar{t}}$ kinematic distributions in the lepton+jets final state in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: Results are presented for an extraction of the top quark Yukawa coupling from top quark-antiquark ($\mathrm{t\bar{t}}$) kinematic distributions in the lepton plus jets final state in proton-proton collisions, based on data collected by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.8 fb$^{-1}$. Corrections from weak boson exchange, including Higgs bosons, between the top quarks can produce large distortions of differential distributions near the energy threshold of $\mathrm{t\bar{t}}$ production. Therefore, precise measurements of these distributions are sensitive to the Yukawa coupling. Top quark events are reconstructed with at least three jets in the final state, and a novel technique is introduced to reconstruct the $\mathrm{t\bar{t}}$ system for events with one missing jet. This technique enhances the experimental sensitivity in the low invariant mass region, $M_\mathrm{t\bar{t}}$. The data yields in $M_\mathrm{t\bar{t}}$, the rapidity difference $|y_\mathrm{t}-y_{\mathrm{\bar{t}}}|$, and the number of reconstructed jets are compared with distributions representing different Yukawa couplings. These comparisons are used to measure the ratio of the top quark Yukawa coupling to its standard model predicted value to be 1.07 $^{+0.34}_{-0.43}$ with an upper limit of 1.67 at the 95% confidence level.


title: Combating the Filter Bubble: Designing for Serendipity in a University Course Recommendation System

abstract: Collaborative filtering based algorithms, including Recurrent Neural Networks (RNN), tend towards predicting a perpetuation of past observed behavior. In a recommendation context, this can lead to an overly narrow set of suggestions lacking in serendipity and inadvertently placing the user in what is known as a "filter bubble." In this paper, we grapple with the issue of the filter bubble in the context of a course recommendation system in production at a public university. Most universities in the United States encourage students to explore developing interests while simultaneously advising them to adhere to course taking norms which progress them towards graduation. These competing objectives, and the stakes involved for students, make this context a particularly meaningful one for investigating real-world recommendation strategies. We introduce a novel modification to the skip-gram model applied to nine years of historic course enrollment sequences to learn course vector representations used to diversify recommendations based on similarity to a student's specified favorite course. This model, which we call multifactor2vec, is intended to improve the semantics of the primary token embedding by also learning embeddings of potentially conflated factors of the token (e.g., instructor). Our offline testing found this model improved accuracy and recall on our course similarity and analogy validation sets over a standard skip-gram. Incorporating course catalog description text resulted in further improvements. We compare the performance of these models to the system's existing RNN-based recommendations with a user study of undergraduates (N = 70) rating six characteristics of their course recommendations. Results of the user study show a dramatic lack of novelty in RNN recommendations and depict the characteristic trade-offs that make serendipity difficult to achieve.


title: Geolocation of Multiple Noncooperative Emitters Using Received Signal Strength: Sparsity, Resolution, and Detectability

abstract: In this paper we investigate the problem of locating multiple non-cooperative radio frequency (RF) emitters using only received signal strength (RSS) data. We assume that the number of emitters is unknown and that individual emitters cannot be distinguished in the RSS data. Moreover, we assume that the environment in which the data has been collected has not been mapped or "fingerprinted" by the prior collection of RSS data. Our primary interest is the limiting resolution that can be obtained by this type of data, and the lowest power emitters that can be detected, as a function of noise level, sensor geometry, and other variables. We formulate the recovery problem as one of sparse approximation or compressed sensing, and investigate an appropriate recovery algorithm for this setting, and use it to illustrate our conclusions. We also include a reconstruction based on sampled data we collected, to illustrate the reasonableness of our parameter choices and conclusions.


title: Incompressible image registration using divergence-conforming B-splines

abstract: Anatomically plausible image registration often requires volumetric preservation. Previous approaches to incompressible image registration have exploited relaxed constraints, ad hoc optimisation methods or practically intractable computational schemes. Divergence-free velocity fields have been used to achieve incompressibility in the continuous domain, although, after discretisation, no guarantees have been provided. In this paper, we introduce stationary velocity fields (SVFs) parameterised by divergence-conforming B-splines in the context of image registration. We demonstrate that sparse linear constraints on the parameters of such divergence-conforming B-Splines SVFs lead to being exactly divergence-free at any point of the continuous spatial domain. In contrast to previous approaches, our framework can easily take advantage of modern solvers for constrained optimisation, symmetric registration approaches, arbitrary image similarity and additional regularisation terms. We study the numerical incompressibility error for the transformation in the case of an Euler integration, which gives theoretical insights on the improved accuracy error over previous methods. We evaluate the proposed framework using synthetically deformed multimodal brain images, and the STACOM11 myocardial tracking challenge. Accuracy measurements demonstrate that our method compares favourably with state-of-the-art methods whilst achieving volume preservation.


title: Quaternionic Analysis, Representation Theory and Physics II

abstract: We develop further quaternionic analysis introducing left and right doubly regular functions. We derive Cauchy-Fueter type formulas for these doubly regular functions that can be regarded as another counterpart of Cauchy's integral formula for the second order pole, in addition to the one studied in the first paper with the same title. We also realize the doubly regular functions as a subspace of the quaternionic-valued functions satisfying a Euclidean version of Maxwell's equations for the electromagnetic field. Then we return to the study of the original quaternionic analogue of Cauchy's second order pole formula and its relation to the polarization of vacuum. We find the decomposition of the space of quaternionic-valued functions into irreducible components that include the spaces of doubly left and right regular functions. Using this decomposition, we show that a regularization of the vacuum polarization diagram is achieved by subtracting the component corresponding to the one-dimensional subrepresentation of the conformal group. After the regularization, the vacuum polarization diagram is identified with a certain second order differential operator which yields a quaternionic version of Maxwell equations. Next, we introduce two types of quaternionic algebras consisting of spaces of scalar-valued and quaternionic-valued functions. We emphasize that these algebra structures are invariant under the action of the conformal Lie algebra. This uses techniques from our study of the vacuum polarization diagram. These algebras are not associative, but we can define an infinite family of n-multiplications, and we conjecture that they have structures of weak cyclic A-infinity algebras. We also conjecture the relation between the multiplication operations of the scalar and non-scalar quaternionic algebras with the n-photon Feynman diagrams in the scalar and ordinary conformal QED.


title: Full color generation with Fano-type resonant HfO$_2$ nanopillars designed by a deep-learning approach

abstract: In contrast to lossy plasmonic metasurfaces (MSs), wideband dielectric MSs comprising of subwavelength nanostructures supporting Mie resonances are of great interest in the visible wavelength range. Here, for the first time to our knowledge, we experimentally demonstrate a reflective MS consisting of a square-lattice array of Hafnia (HfO$_2$) nanopillars to generate a wide color gamut. To design and optimize these MSs, we use a deep-learning algorithm based on a dimensionality reduction technique. Good agreement is observed between simulation and experimental results in yielding vivid and high-quality colors. We envision that these structures not only empower the high-resolution digital displays and sensitive colorimetric biosensors but also can be applied to on-demand applications of beaming in a wide wavelength range down to deep ultraviolet.


title: Quantum Thermodynamics: An introduction to the thermodynamics of quantum information

abstract: This book provides an introduction to the emerging field of quantum thermodynamics, with particular focus on its relation to quantum information and its implications for quantum computers and next generation quantum technologies. The text, aimed at graduate level physics students with a working knowledge of quantum mechanics and statistical physics, provides a brief overview of the development of classical thermodynamics and its quantum formulation in Chapter 1. Chapter 2 then explores typical thermodynamic settings, such as cycles and work extraction protocols, when the working material is genuinely quantum. Finally, Chapter 3 explores the thermodynamics of quantum information processing and introduces the reader to some more state-of-the-art topics in this exciting and rapidly developing research field.

categories: quant-ph cond-mat.mes-hall cond-mat.quant-gas cond-mat.stat-mech


title: Bounded Displacement Non-Equivalence In Substitution Tilings

abstract: In the study of aperiodic order and mathematical models of quasicrystals, questions regarding equivalence relations on Delone sets naturally arise. This work is dedicated to the bounded displacement (BD) equivalence relation, and especially to results concerning instances of non-equivalence. We present a general condition for two Delone sets to be BD non-equivalent, and apply our result to Delone sets associated with tilings of Euclidean space. First we consider substitution tilings, and exhibit a substitution matrix associated with two distinct substitution rules. The first rule generates only periodic tilings, while the second generates tilings for which any associated Delone set is non-equivalent to any lattice in space. As an extension of this result, we introduce arbitrarily many distinct substitution rules associated with a single matrix, with the property that Delone sets generated by distinct rules are non-equivalent. We then turn to the study of mixed substitution tilings, and present a mixed substitution system that generates representatives of continuously many distinct BD equivalence classes.


title: Rod in a train: a mechanical problem of H.Whitney, or Much Ado About Nothing

abstract: In 1941 a mechanical problem about a rod in a moving train (there is a initial position such that rod does not touch the floor while train is moving) was published by R.Courant and H.Robbins in their popular book "What is mathematics?" and attributed to H.Whitney. Many mathematicians, including G.E.Littlewood, A.Broman, T.Poston, I.Stewart, V.Arnold, commented on this problem and its solution based on a continuity argument, and created a lot of confusion. In this paper we follow these developments and discuss at what extent the objections were justified. (In Russian)


title: Multi-object Tracking in Unknown Detection Probability with the PMBM Filter

abstract: This paper focuses on the joint multi-object tracking (MOT) and the estimate of detection probability with the \emph{Poisson multi-Bernoulli mixture} (PMBM) filter. In a majority of multi-object scenarios, the knowledge of detection probability is usually uncertain, which is often estimated offline from the training data. In such cases, online filtering is not allowed or believable, otherwise, significate parameter mismatch will result in biased estimates (state and cardinality of objects). Consequently, the ability of adaptively estimating the detection probability is essential in practice. In this paper, we detail how the detection probability can be estimated accompanied with the state estimates. Besides, closed-form solutions to the proposed method are derived by approximating the intensity of Poisson random finite set (RFS) to a Beta-Gaussian mixture form and density of Bernoulli RFS to a single Beta-Gaussian form. Simulation results show the effectiveness and superiority of the proposed method.


title: Approximate Similarity Search Under Edit Distance Using Locality-Sensitive Hashing

abstract: Edit distance similarity search, also called approximate pattern matching, is a fundamental problem with widespread database applications. The goal of the problem is to preprocess $n$ strings of length $d$, to quickly answer queries $q$ of the form: if there is a database string within edit distance $r$ of $q$, return a database string within edit distance $cr$ of $q$. Previous approaches to this problem either rely on very large (superconstant) approximation ratios $c$, or very small search radii $r$. Outside of a narrow parameter range, these solutions are not competitive with trivially searching through all $n$ strings. In this work give a simple and easy-to-implement hash function that can quickly answer queries for a wide range of parameters. Specifically, our strategy can answer queries in time $\tilde{O}(d3^rn^{1/c})$. The best known practical results require $c \gg r$ to achieve any correctness guarantee; meanwhile, the best known theoretical results are very involved and difficult to implement, and require query time at least $24^r$. Our results significantly broaden the range of parameters for which we can achieve nontrivial bounds, while retaining the practicality of a locality-sensitive hash function. We also show how to apply our ideas to the closely-related Approximate Nearest Neighbor problem for edit distance, obtaining similar time bounds.


title: The Derrida--Retaux conjecture on recursive models

abstract: We are interested in the nearly supercritical regime in a family of max-type recursive models studied by Collet, Eckman, Glaser and Martin and by Derrida and Retaux, and prove that under a suitable integrability assumption on the initial distribution, the free energy vanishes at the transition with an essential singularity with exponent $\tfrac12$. This gives a weaker answer to a conjecture of Derrida and Retaux. Other behaviours are obtained when the integrability condition is not satisfied.


title: Continuous Integration Theater

abstract: Background: Continuous Integration (CI) systems are now the bedrock of several software development practices. Several tools such as TravisCI, CircleCI, and Hudson, that implement CI practices, are commonly adopted by software engineers. However, the way that software engineers use these tools could lead to what we call "Continuous Integration Theater", a situation in which software engineers do not employ these tools effectively, leading to unhealthy CI practices. Aims: The goal of this paper is to make sense of how commonplace are these unhealthy continuous integration practices being employed in practice. Method: By inspecting 1,270 open-source projects that use TravisCI, the most used CI service, we quantitatively studied how common is to use CI (1) with infrequent commits, (2) in a software project with poor test coverage, (3) with builds that stay broken for long periods, and (4) with builds that take too long to run. Results: We observed that 748 ($sim$60%) projects face infrequent commits, which essentially makes the merging process harder. Moreover, we were able to find code coverage information for 51 projects. The average code coverage was 78%, although Ruby projects have a higher code coverage than Java projects (86% and 63%, respectively). However, some projects with very small coverage ($sim$4%) were found. Still, we observed that 85% of the studied projects have at least one broken build that take more than four days to be fixed. Interestingly, very small projects (up to 1,000 lines of code) are the ones that take the longest to fix broken builds. Finally, we noted that, for the majority of the studied projects, the build is executed under the 10 minutes rule of thumb. Conclusions: Our results are important to an increasing community of software engineers that employ CI practices on daily basis but may not be aware of bad practices that are eventually employed.


title: Minimizing the numbers of cliques and cycles of fixed size in an $F$-saturated graph

abstract: This paper considers two important questions in the well-studied theory of graphs that are $F$-saturated. A graph $G$ is called $F$-saturated if $G$ does not contain a subgraph isomorphic to $F$, but the addition of any edge creates a copy of $F$. We first resolve a fundamental question of minimizing the number of cliques of size $r$ in a $K_s$-saturated graph for all sufficiently large numbers of vertices, confirming a conjecture of Kritschgau, Methuku, Tait, and Timmons. We also go further and prove a corresponding stability result. Next we minimize the number of cycles of length $r$ in a $K_s$-saturated graph for all sufficiently large numbers of vertices, and classify the extremal graphs for most values of $r$, answering another question of Kritschgau, Methuku, Tait, and Timmons for most $r$. We then move on to a central and longstanding conjecture in graph saturation made by Tuza, which states that for every graph $F$, the limit $\lim_{n \rightarrow \infty} \frac{\sat(n, F)}{n}$ exists, where $\sat(n, F)$ denotes the minimum number of edges in an $n$-vertex $F$-saturated graph. Pikhurko made progress in the negative direction by considering families of graphs instead of a single graph, and proved that there exists a graph family $\mathcal{F}$ of size $4$ for which $\lim_{n \rightarrow \infty} \frac{\sat(n, \mathcal{F})}{n}$ does not exist (for a family of graphs $\mathcal{F}$, a graph $G$ is called $\mathcal{F}$-saturated if $G$ does not contain a copy of any graph in $\mathcal{F}$, but the addition of any edge creates a copy of a graph in $\mathcal{F}$, and $\sat(n, \mathcal{F})$ is defined similarly). We make the first improvement in 15 years by showing that there exist infinitely many graph families of size $3$ where this limit does not exist. Our construction also extends to the generalized saturation problem when we minimize the number of fixed-size cliques.


title: Zero Spacings of Paraorthogonal Polynomials on the Unit Circle

abstract: We prove some new results about the spacing between neighboring zeros of paraorthogonal polynomials on the unit circle. Our methods also provide new proofs of some existing results. The main tool we will use is a formula for the phase of the appropriate Blaschke product at points on the unit circle.


title: Limits of Sparse Configuration Models and Beyond: Graphexes and Multi-Graphexes

abstract: We investigate structural properties of large, sparse random graphs through the lens of "sampling convergence" (Borgs et. al. (2017)). Sampling convergence generalizes left convergence to sparse graphs, and describes the limit in terms of a "graphex". We introduce a notion of sampling convergence for sequences of multigraphs, and establish the graphex limit for the configuration model, a preferential attachment model, the generalized random graph, and a bipartite variant of the configuration model. The results for the configuration model, preferential attachment model and bipartite configuration model provide necessary and sufficient conditions for these random graph models to converge. The limit for the configuration model and the preferential attachment model is an augmented version of an exchangeable random graph model introduced by Caron and Fox (2017).


title: A Different Angle on Quantum Uncertainty (Measure Angle)

abstract: The uncertainty associated with probing the quantum state is expressed as the effective abundance (measure) of possibilities for its collapse. New kinds of uncertainty limits entailed by quantum description of the physical system arise in this manner.


title: MIDI-Sandwich: Multi-model Multi-task Hierarchical Conditional VAE-GAN networks for Symbolic Single-track Music Generation

abstract: Most existing neural network models for music generation explore how to generate music bars, then directly splice the music bars into a song. However, these methods do not explore the relationship between the bars, and the connected song as a whole has no musical form structure and sense of musical direction. To address this issue, we propose a Multi-model Multi-task Hierarchical Conditional VAE-GAN (Variational Autoencoder-Generative adversarial networks) networks, named MIDI-Sandwich, which combines musical knowledge, such as musical form, tonic, and melodic motion. The MIDI-Sandwich has two submodels: Hierarchical Conditional Variational Autoencoder (HCVAE) and Hierarchical Conditional Generative Adversarial Network (HCGAN). The HCVAE uses hierarchical structure. The underlying layer of HCVAE uses Local Conditional Variational Autoencoder (L-CVAE) to generate a music bar which is pre-specified by the First and Last Notes (FLN). The upper layer of HCVAE uses Global Variational Autoencoder(G-VAE) to analyze the latent vector sequence generated by the L-CVAE encoder, to explore the musical relationship between the bars, and to produce the song pieced together by multiple music bars generated by the L-CVAE decoder, which makes the song both have musical structure and sense of direction. At the same time, the HCVAE shares a part of itself with the HCGAN to further improve the performance of the generated music. The MIDI-Sandwich is validated on the Nottingham dataset and is able to generate a single-track melody sequence (17x8 beats), which is superior to the length of most of the generated models (8 to 32 beats). Meanwhile, by referring to the experimental methods of many classical kinds of literature, the quality evaluation of the generated music is performed. The above experiments prove the validity of the model.


title: Topological Entanglement Entropy in \(d\)-dimensions for Abelian Higher Gauge Theories

abstract: We compute the topological entanglement entropy for a large set of lattice models in $d$-dimensions. It is well known that many such quantum systems can be constructed out of lattice gauge models. For dimensionality higher than two, there are generalizations going beyond gauge theories, which are called higher gauge theories and rely on higher-order generalizations of groups. Our main concern is a large class of $d$-dimensional quantum systems derived from Abelian higher gauge theories. In this paper, we derive a general formula for the bipartition entanglement entropy for this class of models, and from it we extract both the area law and the sub-leading terms, which explicitly depend on the topology of the entangling surface. We show that the entanglement entropy $S_A$ in a sub-region $A$ is proportional to $\log(GSD_{\tilde{A}})$, where \(GSD_{\tilde{A}}\) is the ground state degeneracy of a particular restriction of the full model to \(A\). The quantity $GSD_{\tilde{A}}$ can be further divided into a contribution that scales with the size of the boundary $\partial A$ and a term which depends on the topology of $\partial A$. There is also a topological contribution coming from $A$ itself, that may be non-zero when $A$ has a non-trivial homology. We present some examples and discuss how the topology of $A$ affects the topological entropy. Our formalism allows us to do most of the calculation for arbitrary dimension $d$. The result is in agreement with entanglement calculations for known topological models.


title: Intermittency of dynamical phases in a quantum spin glass

abstract: Answering the question of existence of efficient quantum algorithms for NP-hard problems require deep theoretical understanding of the properties of the low-energy eigenstates and long-time coherent dynamics in quantum spin glasses. We discovered and described analytically the property of asymptotic orthogonality resulting in a new type of structure in quantum spin glass. Its eigen-spectrum is split into the alternating sequence of bands formed by quantum states of two distinct types ($x$ and $z$). Those of $z$-type are non-ergodic extended eigenstates (NEE) in the basis of $\{\sigma_z\}$ operators that inherit the structure of the classical spin glass with exponentially long decay times of Edwards Anderson order parameter at any finite value of transverse field $B_{\perp}$. Those of $x$-type form narrow bands of NEEs that conserve the integer-valued $x$-magnetization. Quantum evolution within a given band of each type is described by a Hamiltonian that belongs to either the ensemble of Preferred Basis Levi matrices ($z$-type) or Gaussian Orthogonal ensemble ($x$-type). We characterize the non-equilibrium dynamics using fractal dimension $D$ that depends on energy density (temperature) and plays a role of thermodynamic potential: $D=0$ in MBL phase, $0<D<1$ in NEE phase, $D\rightarrow 1$ in ergodic phase in infinite temperature limit. MBL states coexist with NEEs in the same range of energies even at very large $B_{\perp}$. Bands of NEE states can be used for new quantum search-like algorithms of population transfer in the low-energy part of spin-configuration space. Remarkably, the intermitted structure of the eigenspectrum emerges in quantum version of a statistically featureless Random Energy Model and is expected to exist in a class of paractically important NP-hard problems that unlike REM can be implemented on a computer with polynomial resources.


title: Going beyond copper: wafer-scale synthesis of graphene on sapphire

abstract: The adoption of graphene in electronics, optoelectronics and photonics is hindered by the difficulty in obtaining high quality material on technologically-relevant substrates, over wafer-scale sizes and with metal contamination levels compatible with industrial requirements. To date, the direct growth of graphene on insulating substrates has proved to be challenging, usually requiring metal-catalysts or yielding defective graphene. In this work, we demonstrate a metal-free approach implemented in commercially available reactors to obtain high-quality monolayer graphene on c-plane sapphire substrates via chemical vapour deposition (CVD). We identify via low energy electron diffraction (LEED), low energy electron microscopy (LEEM) and scanning tunneling microscopy (STM) measurements the Al-rich reconstruction root31R9 of sapphire to be crucial for obtaining epitaxial graphene. Raman spectroscopy and electrical transport measurements reveal high-quality graphene with mobilities consistently above 2000 cm2/Vs. We scale up the process to 4-inch and 6-inch wafer sizes and demonstrate that metal contamination levels are within the limits for back-end-of-line (BEOL) integration. The growth process introduced here establishes a method for the synthesis of wafer-scale graphene films on a technologically viable basis.


title: CS563-QA: A Collection for Evaluating Question Answering Systems

abstract: Question Answering (QA) is a challenging topic since it requires tackling the various difficulties of natural language understanding. Since evaluation is important not only for identifying the strong and weak points of the various techniques for QA, but also for facilitating the inception of new methods and techniques, in this paper we present a collection for evaluating QA methods over free text that we have created. Although it is a small collection, it contains cases of increasing difficulty, therefore it has an educational value and it can be used for rapid evaluation of QA systems.


title: RF current condensation in magnetic islands and associated hysteresis phenomena

abstract: The nonlinear RF current condensation effect suggests that magnetic islands might be well controlled with broader deposition profiles than previously thought possible. To assess this possibility, a simplified energy deposition model in a symmetrised 1D slab geometry is constructed. By limiting the RF wave power that can be absorbed through damping, this model describes also the predicted hysteresis phenomena. Compared to the linear model, the nonlinear effects lead to larger temperature variations, narrower deposition widths, and more robust island stabilisation. Although, in certain regimes, the island centre can be disadvantageously shaded because of the nonlinear effects, in general, the RF condensation effect can take place, with current preferentially generated, advantageously, close to the island centre.


title: A correction to Kallenberg's theorem for jointly exchangeable random measures

abstract: Kallenberg (2005) provided a necessary and sufficient condition for the local finiteness of a jointly exchangeable random measure on $\R_+^2$. Here we note an additional condition that was missing in Kallenberg's theorem, but was implicitly used in the proof. We also provide a counter-example when the additional condition does not hold.


title: Validity of the Einstein Hole Argument

abstract: Arguing from his "hole" thought experiment, Einstein became convinced that, in cases in which the energy-momentum-tensor source vanishes in a spacetime hole, a solution to his general relativistic field equation cannot be uniquely determined by that source. After reviewing the definition of active diffeomorphisms, this paper uses them to outline a mathematical proof of Einstein's result. The relativistic field equation is shown to have multiple solutions, just as Einstein thought. But these multiple solutions can be distinguished by the different physical meaning that each metric solution attaches to the local coordinates used to write it. Thus the hole argument, while formally correct, does not prohibit the subsequent rejection of spurious solutions and the selection of a physically unique metric. This conclusion is illustrated using the Schwarzschild metric. It is suggested that the Einstein hole argument therefore cannot be used to argue against substantivalism.


title: E-Sports Talent Scouting Based on Multimodal Twitch Stream Data

abstract: We propose and investigate feasibility of a novel task that consists in finding e-sports talent using multimodal Twitch chat and video stream data. In that, we focus on predicting the ranks of Counter-Strike: Global Offensive (CS:GO) gamers who broadcast their games on Twitch. During January 2019-April 2019, we have built two Twitch stream collections: One for 425 publicly ranked CS:GO gamers and one for 9,928 unranked CS:GO gamers. We extract neural features from video, audio and text chat data and estimate modality-specific probabilities for a gamer to be top-ranked during the data collection time-frame. A hierarchical Bayesian model is then used to pool the evidence across modalities and generate estimates of intrinsic skill for each gamer. Our modeling is validated through correlating the intrinsic skill predictions with May 2019 ranks of the publicly profiled gamers.


title: Harnessing Spectral Singularities in Non-Hermitian Cylindrical Structures

abstract: Non-Hermitian systems characterized by suitable spatial distributions of gain and loss can exhibit "spectral singularities" in the form of zero-width resonances associated to real-frequency poles in the scattering operator. Here, we study this intriguing phenomenon in connection with cylindrical geometries, and explore possible applications to controlling and tailoring in unconventional ways the scattering response of sub-wavelength and wavelength-sized objects. Among the possible implications and applications, we illustrate the additional degrees of freedom available in the scattering-absorption-extinction tradeoff, and address the engineering of zero-forward-scattering, transverse scattering, and gain-controlled reconfigurability of the scattering pattern, also paying attention to stability issues. Our results may open up new vistas in active and reconfigurable nanophotonics platforms.


title: A Short Proof of the Toughness of Delaunay Triangulations

abstract: We present a self-contained short proof of the seminal result of Dillencourt (SoCG 1987 and DCG 1990) that Delaunay triangulations, of planar point sets in general position, are 1-tough. An important implication of this result is that Delaunay triangulations have perfect matchings. Another implication of our result is a proof of the conjecture of Aichholzer et al. (2010) that at least $n$ points are required to block any $n$-vertex Delaunay triangulation


title: Impacts of radiative corrections on measurements of lepton flavour universality in $B \to D \ell \nu_{\ell}$ decays

abstract: Radiative corrections to $B \to D \ell \nu_{\ell}$ decays can have an impact on predictions and measurements of the lepton universality ratios $\mathcal{R}(D^+)$ and $\mathcal{R}(D^0)$. These proceedings summarise a study on the comparison between recent calculations of soft-photon corrections on these ratios and the corrections simulated by the PHOTOS package. Also the impact of Coulomb interactions, not simulated in PHOTOS, is discussed. Using pseudo-experiments, the effect of high-energy photon emission is studied in an LHCb-like environment, showing a bias of up to 7\% on measurements of $\mathcal{R}(D)$.


title: Learning from satisfying assignments under continuous distributions

abstract: What kinds of functions are learnable from their satisfying assignments? Motivated by this simple question, we extend the framework of De, Diakonikolas, and Servedio [DDS15], which studied the learnability of probability distributions over $\{0,1\}^n$ defined by the set of satisfying assignments to "low-complexity" Boolean functions, to Boolean-valued functions defined over continuous domains. In our learning scenario there is a known "background distribution" $\mathcal{D}$ over $\mathbb{R}^n$ (such as a known normal distribution or a known log-concave distribution) and the learner is given i.i.d. samples drawn from a target distribution $\mathcal{D}_f$, where $\mathcal{D}_f$ is $\mathcal{D}$ restricted to the satisfying assignments of an unknown low-complexity Boolean-valued function $f$. The problem is to learn an approximation $\mathcal{D}'$ of the target distribution $\mathcal{D}_f$ which has small error as measured in total variation distance. We give a range of efficient algorithms and hardness results for this problem, focusing on the case when $f$ is a low-degree polynomial threshold function (PTF). When the background distribution $\mathcal{D}$ is log-concave, we show that this learning problem is efficiently solvable for degree-1 PTFs (i.e.,~linear threshold functions) but not for degree-2 PTFs. In contrast, when $\mathcal{D}$ is a normal distribution, we show that this learning problem is efficiently solvable for degree-2 PTFs but not for degree-4 PTFs. Our hardness results rely on standard assumptions about secure signature schemes.


title: Introducing Astrocytes on a Neuromorphic Processor: Synchronization, Local Plasticity and Edge of Chaos

abstract: While there is still a lot to learn about astrocytes and their neuromodulatory role in the spatial and temporal integration of neuronal activity, their introduction to neuromorphic hardware is timely, facilitating their computational exploration in basic science questions as well as their exploitation in real-world applications. Here, we present an astrocytic module that enables the development of a spiking Neuronal-Astrocytic Network (SNAN) into Intel's Loihi neuromorphic chip. The basis of the Loihi module is an end-to-end biophysically plausible compartmental model of an astrocyte that simulates the intracellular activity in response to the synaptic activity in space and time. To demonstrate the functional role of astrocytes in SNAN, we describe how an astrocyte may sense and induce activity-dependent neuronal synchronization, switch on and off spike-time-dependent plasticity (STDP) to introduce single-shot learning, and monitor the transition between ordered and chaotic activity at the synaptic space. Our module may serve as an extension for neuromorphic hardware, by either replicating or exploring the distinct computational roles that astrocytes have in forming biological intelligence.


title: Renormalization in Nonminimal Lorentz-Violating Field Theory

abstract: We provide the first step towards renormalization in a nonminimal Lorentz-violating model consisting of normal scalars and modified fermions with mass dimension five operators. We compute the radiative corrections corresponding to the scalar self-energy, and we show that some divergencies are improved and in the scalar sector they are finite. The pole mass of the scalar two-point function is found and shown to lead to modifications of asymptotic states.


title: On internal structure of smaller domains in domain coarsening dynamics of spontaneous Z_2-symmetry breaking in two dimensions

abstract: The internal structure of domains smaller than the characteristic size in domain coarsening dynamics of $Z_2$ symmetry breaking is evaluated theoretically for different phase ordering systems in two dimensions. In the previous works on (non-) conserved Ising systems and binary superfluids, the statistical properties of smaller domains are analyzed by assuming that the contribution from its internal structure is negligible. It is shown that this assumption is justified analytically with respect to the statistical quantities, such as domain area, domain-wall length and superfluid circulation, according to the empirical dynamic scaling law for the smaller domains.


title: Evolving the Hearthstone Meta

abstract: Balancing an ever growing strategic game of high complexity, such as Hearthstone is a complex task. The target of making strategies diverse and customizable results in a delicate intricate system. Tuning over 2000 cards to generate the desired outcome without disrupting the existing environment becomes a laborious challenge. In this paper, we discuss the impacts that changes to existing cards can have on strategy in Hearthstone. By analyzing the win rate on match-ups across different decks, being played by different strategies, we propose to compare their performance before and after changes are made to improve or worsen different cards. Then, using an evolutionary algorithm, we search for a combination of changes to the card attributes that cause the decks to approach equal, 50% win rates. We then expand our evolutionary algorithm to a multi-objective solution to search for this result, while making the minimum amount of changes, and as a consequence disruption, to the existing cards. Lastly, we propose and evaluate metrics to serve as heuristics with which to decide which cards to target with balance changes.


title: Efficient Circuit Simulation in MapReduce

abstract: The MapReduce framework has firmly established itself as one of the most widely used parallel computing platforms for processing big data on tera- and peta-byte scale. Approaching it from a theoretical standpoint has proved to be notoriously difficult, however. In continuation of Goodrich et al.'s early efforts, explicitly espousing the goal of putting the MapReduce framework on footing equal to that of long-established models such as the PRAM, we investigate the obvious complexity question of how the computational power of MapReduce algorithms compares to that of combinational Boolean circuits commonly used for parallel computations. Relying on the standard MapReduce model introduced by Karloff et al. a decade ago, we develop an intricate simulation technique to show that any problem in NC (i.e., a problem solved by a logspace-uniform family of Boolean circuits of polynomial size and a depth polylogarithmic in the input size) can be solved by a MapReduce computation in O(T(n)/ log n) rounds, where n is the input size and T(n) is the depth of the witnessing circuit family. Thus, we are able to closely relate the standard, uniform NC hierarchy modeling parallel computations to the deterministic MapReduce hierarchy DMRC by proving that NC^(i+1) is contained in DMRC^i for all natural i, including 0. Besides the theoretical significance, this result that has important applied aspects as well. In particular, we show for all problems in NC^1---many practically relevant ones such as integer multiplication and division, the parity function, and recognizing balanced strings of parentheses being among these---how to solve them in a constant number of deterministic MapReduce rounds.