title: Molecular dynamics studies of interactions between Arg9(nona-arginine) and a DOPC/DOPG(4:1) membrane

abstract: It has been known that the uptake mechanisms of cell-penetrating peptides(CPPs) depend on the experimental conditions such as concentration of peptides, lipid composition, temperature, etc. In this study we investigate the temperature dependence of the penetration of Arg9s into a DOPC/DOPG(4:1) membrane using molecular dynamics(MD) simulations at two different temperatures, T = 310 K and T = 288 K. Although it is difficult to identify the temperature dependence because of having only one single simulation at each temperature and no evidence of translocation of Arg9s across the membrane at both temperatures, our simulations suggest that followings are strongly correlated with the penetration of Arg9s: a number of water molecules coordinated by Arg9s, electrostatic energy between Arg9s and the lipids molecules. We also present how Arg9s change a bending rigidity of the membrane and how a collective behavior between Arg9s enhances the penetration and the membrane bending. Our analyses can be applicable to any cell-penetrating peptides(CPPs) to investigate their interactions with various membranes using MD simulations.


title: Open quantum systems decay across time

abstract: The description of an open quantum system's decay almost always requires several approximations as to remain tractable. Here, we first revisit the meaning, domain and seeming contradictions of a few of the most widely used of such approximations: semi-group Markovianity, linear response theory, Wigner--Weisskopf and rotating-wave approximation. Secondly, we derive an effective time-dependent decay theory and corresponding generalized quantum regression relations for an open quantum system linearly coupled to an environment. This theory covers all timescales, and subsumes the Markovian and linear-response results as limiting cases. Finally, we apply our theory to the phenomenon of quantum friction.


title: From Probability to Consilience: How Explanatory Values Implement Bayesian Reasoning

abstract: Recent work in cognitive science has uncovered a diversity of explanatory values, or dimensions along which we judge explanations as better or worse. We propose a Bayesian account of how these values fit together to guide explanation. The resulting taxonomy provides a set of predictors for which explanations people prefer and shows how core values from psychology, statistics, and the philosophy of science emerge from a common mathematical framework. In addition to operationalizing the explanatory virtues associated with, for example, scientific argument-making, this framework also enables us to reinterpret the explanatory vices that drive conspiracy theories, delusions, and extremist ideologies.


title: Pseudo-ternary LiBH4-LiCl-P2S5 system as structurally disordered bulk electrolyte for all-solid-state lithium batteries

abstract: The properties of the mixed system LiBH4 LiCl P2S5 are studied with respect to all-solid-state batteries. The studied material undergoes an amorphization upon heating above 601C, accompanied with increased Li+ conductivity beneficial for battery electrolyte applications. The measured ionic conductivity is 10-3 Scm-1 at room temperature with an activation energy of 0.40(2) eV after amorphization. Structural analysis and characterization of the material suggest that BH4 groups and PS4 may belong to the same molecular structure, where Cl ions interplay to accommodate the structural unit. Thanks to its conductivity, ductility and electrochemical stability (up to 5 V, Au vs. Li+/Li), this new electrolyte is successfully tested in battery cells operated with a cathode material (layered TiS2, theo. capacity 239 mAh g-1) and Li anode resulting in 93% capacity retention (10 cycles) and notable cycling stability under the current density 12 mA g-1 (0.05C-rate) at 501C. Further advanced characterisation by means of operando synchrotron X-ray diffraction in transmission mode contributes explicitly to a better understanding of the (de)lithiation processes of solid-state battery electrodes operated at moderate temperatures.


title: Optimizing Neural Networks via Koopman Operator Theory

abstract: Koopman operator theory, a powerful framework for discovering the underlying dynamics of nonlinear dynamical systems, was recently shown to be intimately connected with neural network training. In this work, we take the first steps in making use of this connection. As Koopman operator theory is a linear theory, a successful implementation of it in evolving network weights and biases offers the promise of accelerated training, especially in the context of deep networks, where optimization is inherently a non-convex problem. We show that Koopman operator theoretic methods allow for accurate predictions of weights and biases of feedforward, fully connected deep networks over a non-trivial range of training time. During this window, we find that our approach is >10x faster than various gradient descent based methods (e.g. Adam, Adadelta, Adagrad), in line with our complexity analysis. We end by highlighting open questions in this exciting intersection between dynamical systems and neural network theory. We highlight additional methods by which our results could be expanded to broader classes of networks and larger training intervals, which shall be the focus of future work.


title: First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

abstract: Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as high-$T_c$ superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of $2.88\unicode{xb0}$ and simulate how its electronic states evolve when doping level and out-of-plane electric field are gate-tuned. Our calculations show that, at doping levels between two electrons and four holes per moir\'{e} unit cell, a ferromagnetic solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this ferromagnetic state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moir\'{e} unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. Moreover, upon increasing the out-of-plane electric field at zero doping, a transition from the ferromagnetic state to the nonmagnetic one is seen. We also analyze the interlayer bonding character due to the flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the ferromagnetic state at a certain doping level.


title: A Self Consistent Field Formulation of Excited State Mean Field Theory

abstract: We show that, as in Hartree Fock theory, the orbitals for excited state mean field theory can be optimized via a self-consistent one-electron equation in which electron-electron repulsion is accounted for through mean field operators. In addition to showing that this excited state ansatz is sufficiently close to a mean field product state to admit a one-electron formulation, this approach brings the orbital optimization speed to within roughly a factor of two of ground state mean field theory. The approach parallels Hartree Fock theory in multiple ways, including the presence of a commutator condition, a one-electron mean-field working equation, and acceleration via direct inversion in the iterative subspace. When combined with a configuration interaction singles Davidson solver for the excitation coefficients, the self consistent field formulation dramatically reduces the cost of the theory compared to previous approaches based on quasi-Newton descent.


title: The Trouble with Water: Condensation, Circulation and Climate

abstract: This article discussesl a few of the problems that arise in geophysical fluid dynamics and climate that are associated with the presence of moisture in the air, its condensation and release of latent heat. Our main focus is Earth's atmosphere but we also discuss how these problems might manifest themselves on other planetary bodies, with particular attention to Titan where methane takes on the role of water. GFD has traditionally been concerned with understanding the very basic problems that lie at the foundation of dynamical meteorology and ocean\-ography. Conventionally, and a little ironically, the subject mainly considers `dry' fluids, meaning it does not concern itself overly much with phase changes. The subject is often regarded as dry in another way, because it does not consider problems perceived as relevant to the real world, such as clouds or rainfall, which have typically been the province of complicated numerical models. Those models often rely on parameterizations of unresolved processes, parameterizations that may work very well but that often have a semi-empirical basis. The consequent dichotomy between the foundations and the applications prevents progress being made that has both a secure basis in scientific understanding and a relevance to the Earth's climate, especially where moisture is concerned. The dichotomy also inhibits progress in understanding the climate of other planets, where observations are insufficient to tune the parameterizations that weather and climate models for Earth rely upon, and a more fundamental approach is called for. Here we discuss four diverse examples of the problems with moisture: the determination of relative humidity and cloudiness; the transport of water vapor and its possible change under global warming; the moist shallow water equations and the Madden-Julian Oscillation; and the hydrology cycle on other planetary bodies.


title: Large optical tunability from charge density waves in 1T-TaS$_2$ under incoherent illumination

abstract: Strongly correlated materials possess a complex energy landscape and host many interesting physical phenomena, including charge density waves (CDWs). CDWs have been observed and extensively studied in many materials since their first discovery in 1972. Yet, they present ample opportunities for discovery. Here, we report a large tunability in the optical response of a quasi-2D CDW material, 1T-TaS$_2$, upon incoherent light illumination at room temperature. We show that the observed tunability is a consequence of light-induced rearrangement of CDW stacking across the layers of 1T-TaS$_2$. Our model, based on this hypothesis, agrees reasonably well with experiments suggesting that the interdomain CDW interaction is a vital knob to control the phase of strongly correlated materials.


title: Mapping the co-evolution of artificial intelligence, robotics, and the internet of things over 20 years (1998-2017)

abstract: Understanding the emergence, co-evolution, and convergence of science and technology (S&T) areas offers competitive intelligence for researchers, managers, policy makers, and others. The resulting data-driven decision support helps set proper research and development (R&D) priorities; develop future S&T investment strategies; monitor key authors, organizations, or countries; perform effective research program assessment; and implement cutting-edge education/training efforts. This paper presents new funding, publication, and scholarly network metrics and visualizations that were validated via expert surveys. The metrics and visualizations exemplify the emergence and convergence of three areas of strategic interest: artificial intelligence (AI), robotics, and internet of things (IoT) over the last 20 years (1998-2017). For 32,716 publications and 4,497 NSF awards, we identify their conceptual space (using the UCSD map of science), geospatial network, and co-evolution landscape. The findings demonstrate how the transition of knowledge (through cross-discipline publications and citations) and the emergence of new concepts (through term bursting) create a tangible potential for interdisciplinary research and new disciplines.


title: Online adaptation in robots as biological development provides phenotypic plasticity

abstract: The ability of responding to environmental stimuli with appropriate actions is a property shared by all living organisms, and it is also sought in the design of robotic systems. Phenotypic plasticity provides a way for achieving this property as it characterises those organisms that, from one genotype, can express different phenotypes in response to different environments, without involving genetic modifications. In this work we study phenotypic plasticity in robots that are equipped with online sensor adaptation. We show that Boolean network controlled robots can attain navigation with collision avoidance by adapting the coupling between proximity sensors and their controlling network without changing its structure. In other terms, these robots, while being characterised by one genotype (i.e. the network) can express a phenotype among many that is suited for the specific environment. We also show that the dynamical regime that makes it possible to attain the best overall performance is the critical one, bringing further evidence to the hypothesis that natural and artificial systems capable of optimally balancing robustness and adaptivity are critical.


title: How to Spread a Rumor: Call Your Neighbors or Take a Walk?

abstract: We study the problem of randomized information dissemination in networks. We compare the now standard PUSH-PULL protocol, with agent-based alternatives where information is disseminated by a collection of agents performing independent random walks. In the VISIT-EXCHANGE protocol, both nodes and agents store information, and each time an agent visits a node, the two exchange all the information they have. In the MEET-EXCHANGE protocol, only the agents store information, and exchange their information with each agent they meet. We consider the broadcast time of a single piece of information in an $n$-node graph for the above three protocols, assuming a linear number of agents that start from the stationary distribution. We observe that there are graphs on which the agent-based protocols are significantly faster than PUSH-PULL, and graphs where the converse is true. We attribute the good performance of agent-based algorithms to their inherently fair bandwidth utilization, and conclude that, in certain settings, agent-based information dissemination, separately or in combination with PUSH-PULL, can significantly improve the broadcast time. The graphs considered above are highly non-regular. Our main technical result is that on any regular graph of at least logarithmic degree, PUSH-PULL and VISIT-EXCHANGE have the same asymptotic broadcast time. The proof uses a novel coupling argument which relates the random choices of vertices in PUSH-PULL with the random walks in VISIT-EXCHANGE. Further, we show that the broadcast time of MEET-EXCHANGE is asymptotically at least as large as the other two's on all regular graphs, and strictly larger on some regular graphs. As far as we know, this is the first systematic and thorough comparison of the running times of these very natural information dissemination protocols.


title: A machine learning approach to Bayesian parameter estimation

abstract: Bayesian estimation is a powerful theoretical paradigm for the operation of quantum sensors. However, the Bayesian method for statistical inference generally suffers from demanding calibration requirements that have so far restricted its use to proof-of-principle experiments. In this theoretical study, we formulate parameter estimation as a classification task and use artificial neural networks to efficiently perform Bayesian estimation. We show that the network's posterior distribution is centered at the true (unknown) value of the parameter within an uncertainty given by the inverse Fisher information, representing the ultimate sensitivity limit for the given apparatus. When only a limited number of calibration measurements are available, our machine-learning based procedure outperforms standard calibration methods. Thus, our work paves the way for Bayesian quantum sensors which can benefit from efficient optimization methods, such as in adaptive schemes, and take advantage of complex non-classical states. These capabilities can significantly enhance the sensitivity of future devices.


title: Exploring Context-Aware Conversational Agents in Software Development

abstract: Software development is a complex endeavor that depends on a wide variety of contextual factors involving a large amount of distributed information. This knowledge could include: technology-related tasks, software operating environments and stakeholder requirements. A major roadblock to using this knowledge in software development is that most of this information is implicit and captured in the developers' minds (tacit) or spread through volumes of documentation. Developers, as they work often have to maintain mental models of these tasks as they produce the software. As a result, context can be easily lost or forgotten and developers often use trial-and-error approaches while finishing the project. This study aims at analyzing whether supporting software developers with a chatbot during task execution can improve the overall development experience. The chatbot can assist the developers in executing different tasks based on implicit contextual information. We propose an implementation to explore the viability of using textual chatbots to assist developers automatically and proactively with software development project activities that recur.


title: How Gamification Affects Software Developers: Cautionary Evidence from a Natural Experiment on GitHub

abstract: We examine how the behavior of software developers changes in response to removing gamification elements from GitHub, an online platform for collaborative programming and software development. We find that the unannounced removal of daily activity streak counters from the user interface (from user profile pages) was followed by significant changes in behavior. Long-running streaks of activity were abandoned and became less common. Weekend activity decreased and days in which developers made a single contribution became less common. Synchronization of streaking behavior in the platform's social network also decreased, suggesting that gamification is a powerful channel for social influence. Focusing on a set of software developers that were publicly pursuing a goal to make contributions for 100 days in a row, we find that some of these developers abandon this quest following the removal of the public streak counter. Our findings provide evidence for the significant impact of gamification on the behavior of developers on large collaborative programming and software development platforms. They urge caution: gamification can steer the behavior of software developers in unexpected and unwanted directions.


title: Semilattice ordered algebras with constants

abstract: We continue our studies on semilattice ordered algebras. This time we accept constants in the type of algebras. We investigate identities satisfied by such algebras and describe the free objects in varieties of semilattice ordered algebras with constants.


title: The reciprocal sum of divisors of Mersenne numbers

abstract: We investigate various questions concerning the reciprocal sum of divisors, or prime divisors, of the Mersenne numbers $2^n-1$. Conditional on the Elliott-Halberstam Conjecture and the Generalized Riemann Hypothesis, we determine $\max_{n\le x} \sum_{p \mid 2^n-1} 1/p$ to within $o(1)$ and $\max_{n\le x} \sum_{d\mid 2^n-1}1/d$ to within a factor of $1+o(1)$, as $x\to\infty$. This refines, conditionally, earlier estimates of Erd\H{o}s and Erd\H{o}s-Kiss-Pomerance. Conditionally (only) on GRH, we also determine $\sum 1/d$ to within a factor of $1+o(1)$ where $d$ runs over all numbers dividing $2^n-1$ for some $n\le x$. This conditionally confirms a conjecture of Pomerance and answers a question of Murty-Rosen-Silverman. Finally, we show that both $\sum_{p\mid 2^n-1} 1/p$ and $\sum_{d\mid 2^n-1}1/d$ admit continuous distribution functions in the sense of probabilistic number theory.


title: On tensor rank and commuting matrices

abstract: Obtaining superlinear lower bounds on tensor rank is a major open problem in complexity theory. In this paper we propose a generalization of the approach used by Strassen in the proof of his 3n/2 border rank lower bound. Our approach revolves around a problem on commuting matrices: Given matrices Z_1,...,Z_p of size n and an integer r>n, are there commuting matrices Z'_1,...,Z'_p of size r such that every Z_k is embedded as a submatrix in the top-left corner of Z'_k? As one of our main results, we show that this question always has a positive answer for r larger than rank(T)+n, where T denotes the tensor with slices Z_1,..,Z_p. Taking the contrapositive, if one can show for some specific matrices Z_1,...,Z_p and a specific integer r that this question has a negative answer, this yields the lower bound rank(T) > r-n. There is a little bit of slack in the above rank(T)+n bound, but we also provide a number of exact characterizations of tensor rank and symmetric rank, for ordinary and symmetric tensors, over the fields of real and complex numbers. Each of these characterizations points to a corresponding variation on the above approach. In order to explain how Strassen's theorem fits within this framework we also provide a self-contained proof of his lower bound.


title: Eikonal phase matrix, deflection angle and time delay in effective field theories of gravity

abstract: The eikonal approximation is an ideal tool to extract classical observables in gauge theory and gravity directly from scattering amplitudes. Here we consider effective theories of gravity where in addition to the Einstein-Hilbert term we include non-minimal couplings of the type $R^3$, $R^4$ and $FFR$. In particular, we study the scattering of gravitons and photons of frequency $\omega$ off heavy scalars of mass $m$ in the limit $m\gg \omega \gg |\vec{q}\,|$, where $\vec{q}$ is the momentum transfer. The presence of non-minimal couplings induces helicity-flip processes which survive the eikonal limit, thereby promoting the eikonal phase to an eikonal phase matrix. We obtain the latter from the relevant two-to-two helicity amplitudes that we compute up to one-loop order, and confirm that the leading-order terms in $\omega$ exponentiate \`{a} la Amati, Ciafaloni and Veneziano. From the eigenvalues of the eikonal phase matrix we then extract two physical observables, to 2PM order: the classical deflection angle and Shapiro time delay/advance. Whenever the classical expectation of helicity conservation of the massless scattered particle is violated, i.e. the eigenvalues of the eikonal matrix are non-degenerate, causality violation due to time advance is a generic possibility for small impact parameter. We show that for graviton scattering in the $R^4$ and $FFR$ theories, time advance is circumvented if the couplings of these interactions satisfy certain positivity conditions, while it is unavoidable for graviton scattering in the $R^3$ theory and photon scattering in the $FFR$ theory. The scattering processes we consider mimic the deflection of photons and gravitons off spinless heavy objects such as black~holes.


title: First-principles theory of infrared vibrational spectroscopy in metals and semimetals: application to graphite

abstract: We develop an ab initio method to simulate the infrared vibrational response of metallic systems in the framework of time-dependent density functional perturbation theory. By introducing a generalized frequency-dependent Born effective charge tensor, we show that phonon peaks in the reflectivity of metals can be always described by a Fano function, whose shape is determined by the complex nature of the frequency-dependent effective charges and electronic dielectric tensor. The IR vibrational properties of graphite, chosen as a representative test case to benchmark our method, are found to be accurately reproduced. Our approach offers a first-principle scheme for the prediction and understanding of IR reflectance spectra of metals, that may represent one of the few available tools of investigation of these materials when subjected to extremely high-pressure conditions.


title: RODE-Net: Learning Ordinary Differential Equations with Randomness from Data

abstract: Random ordinary differential equations (RODEs), i.e. ODEs with random parameters, are often used to model complex dynamics. Most existing methods to identify unknown governing RODEs from observed data often rely on strong prior knowledge. Extracting the governing equations from data with less prior knowledge remains a great challenge. In this paper, we propose a deep neural network, called RODE-Net, to tackle such challenge by fitting a symbolic expression of the differential equation and the distribution of parameters simultaneously. To train the RODE-Net, we first estimate the parameters of the unknown RODE using the symbolic networks \cite{long2019pde} by solving a set of deterministic inverse problems based on the measured data, and use a generative adversarial network (GAN) to estimate the true distribution of the RODE's parameters. Then, we use the trained GAN as a regularization to further improve the estimation of the ODE's parameters. The two steps are operated alternatively. Numerical results show that the proposed RODE-Net can well estimate the distribution of model parameters using simulated data and can make reliable predictions. It is worth noting that, GAN serves as a data driven regularization in RODE-Net and is more effective than the $\ell_1$ based regularization that is often used in system identifications.

categories: math.NA cs.LG cs.NA physics.comp-ph stat.ML


title: The Drift Chamber of the MEG II experiment

abstract: The MEG experiment at the Paul Scherrer Institut searches for the charged-Lepton-Flavor-Violating mu+ -> e+ gamma decay. MEG has already set the world best upper limit on the branching ratio: BR<4.2x10^-13 @ 90% C.l. An upgrade (MEG II) of the whole detector has been approved to obtain a substantial increase of sensitivity. Currently MEG II is completing the upgrade of the various detectors, an engineering run and a pre-commissioning run were carried out during 2018 and 2019. The new positron tracker is a unique volume, ultra-light He based cylindrical drift chamber (CDCH), with high granularity: 9 layers of 192 square drift cells, ~6-9 mm wide, consist of ~12000 wires in a full stereo configuration. To ensure the electrostatic stability of the drift cells a new wiring strategy should be developed due to the high wire density (12 wires/cm^2 ), the stringent precision requirements on the wire position and uniformity of the wire mechanical tension (better than 0.5 g) The basic idea is to create multiwire frames, by soldering a set of (16 or 32) wires on 40 um thick custom wire-PCBs. Multiwire frames and PEEK spacers are overlapped alternately along the radius, to set the proper cell width, in each of the twelve sectors defined by the spokes of the rudder wheel shaped end-plates. Despite to the conceptual simplicity of the assembling strategies, the building of the multiwire frames, with the set requirements, imposes a use of an automatic wiring system. The MEG II CDCH is the first cylindrical drift chamber ever designed and built in a modular way and it will allow to track positrons, with a momentum greater than 45 MeV/c, with high efficiency by using a very small amount of material, 1.5x10^-3 X0 . We describe the CDCH design and construction, the wiring phase at INFN-Lecce, the choice of the wires, their mechanical properties, the assembly and sealing at INFN-Pisa and the commissioning.


title: The Neural Tangent Link Between CNN Denoisers and Non-Local Filters

abstract: Convolutional Neural Networks (CNNs) are now a well-established tool for solving computational imaging problems. Modern CNN-based algorithms obtain state-of-the-art performance in diverse image restoration problems. Furthermore, it has been recently shown that, despite being highly overparameterized, networks trained with a single corrupted image can still perform as well as fully trained networks. We introduce a formal link between such networks through their neural tangent kernel (NTK), and well-known non-local filtering techniques, such as non-local means or BM3D. The filtering function associated with a given network architecture can be obtained in closed form without need to train the network, being fully characterized by the random initialization of the network weights. While the NTK theory accurately predicts the filter associated with networks trained using standard gradient descent, our analysis shows that it falls short to explain the behaviour of networks trained using the popular Adam optimizer. The latter achieves a larger change of weights in hidden layers, adapting the non-local filtering function during training. We evaluate our findings via extensive image denoising experiments.


title: Self-supervised Training of Graph Convolutional Networks

abstract: Graph Convolutional Networks (GCNs) have been successfully applied to analyze non-grid data, where the classical convolutional neural networks (CNNs) cannot be directly used. One similarity shared by GCNs and CNNs is the requirement of massive amount of labeled data for network training. In addition, GCNs need the adjacency matrix as input to define the relationship between those non-grid data, which leads to all of data including training, validation and test data typically forms only one graph structures data for training. Furthermore, the adjacency matrix is usually pre-defined and stationary, which makes the data augmentation strategies cannot be employed on the constructed graph structures data to augment the amount of training data. To further improve the learning capacity and model performance under the limited training data, in this paper, we propose two types of self-supervised learning strategies to exploit available information from the input graph structure data itself. Our proposed self-supervised learning strategies are examined on two representative GCN models with three public citation network datasets - Citeseer, Cora and Pubmed. The experimental results demonstrate the generalization ability as well as the portability of our proposed strategies, which can significantly improve the performance of GCNs with the power of self-supervised learning in improving feature learning.


title: Charged particle guiding and beam splitting with auto-ponderomotive potentials on a chip

abstract: We report guiding and manipulation of charged particle beams by means of electrostatic optics based on a principle similar to the electrodynamic Paul trap. We use hundreds of electrodes fabricated on planar substrates and supplied with static voltages to create a ponderomotive potential for charged particles in motion. Shape and strength of the potential can be locally tailored by the electrodes' layout and the applied voltages, enabling the control of charged particle beams within precisely engineered effective potentials. We demonstrate guiding of electrons and ions for a large range of energies (from 20 to 5000 eV) and masses (5E-4 to 131 atomic mass units) as well as electron beam splitting as a proof-of-concept for more complex beam manipulation. Simultaneous confinement of charged particles with different masses is possible, as well as guiding of electrons with energies in the keV regime, and the creation of highly customizable potential landscapes, which is all hard to impossible with conventional electrodynamic Paul traps.


title: The Panchromatic Afterglow of GW170817: The full uniform dataset, modeling, comparison with previous results and implications

abstract: We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.


title: The non-linear dual gravity equation of motion in eleven dimensions

abstract: We derive the non-linear dual graviton equation of motion in eleven dimensions in the context of E theory.


title: The First Light Curve Solutions and Period Changes of EW Psc and HN Psc

abstract: The first light curve solutions of the stars EW and HN in Pisces constellation are presented. Photometry, and its' periodic changes are calculated and discussed. The analysis of O-C diagram done by MCMC approach in OCFit code and the new ephemeris provided for two binary systems. The light curve solutions obtained. The results show that EW Psc is a near contact eclipsing binary system with a photometric mass ratio q = 0.587, and the fillout factor -0.034 and -0.018 for primary and secondary components, respectively. The solution results also show that the system HN Psc is a weak-contact W UMa eclipsing binary with a photometric mass ratio q = 0.853, and with fillout factor 5.6%. The light curves solutions required the cold spot accounting for the O'Connell effect.


title: The strange critical endpoint and isentropic trajectories in an Extended PNJL Model with Eight Quark Interactions

abstract: In this work, we explore the possible existence of several critical endpoints in the phase diagram of strongly interacting matter using an extended PNJL model with 't Hooft determinant and eight quark interactions in the up, down and strange sectors. Besides, we also study the isentropic trajectories crossing both (light and strange) chiral phase transitions and around the critical endpoint in both the crossover and first-order transition regions.


title: Graphene Nanoribbons Under Axial Compressive and Point Tensile Stresses

abstract: The geometric, electronic and magnetic properties of strained graphene nanoribbons were investigated using spin polarized calculations within the framework of density functional theory. Cases of compressive stress along the longer axis of a nanoribbon and tensile stress at the midpoint and perpendicular to the plane of the nanoribbon were considered. Significant structural changes were observed including the formation of nanoripples. The calculated electronic and magnetic properties strongly depend on the size and shape of nanoribbons. The tunable magnetic properties of strained nanoribbons can be employed for designing magnetic nano-switches.


title: Exchange splitting and exchange-induced non-reciprocal photonic behavior of graphene in CrI3-graphene vdW heterostructures

abstract: The exchange splitting and resulting electromagnetic response of graphene in a monolayer chromium triiodide(CrI3) graphene van der Waals(vdW) heterostructure are considered using a combination of density functional theory and electromagnetic calculations. Although the effective exchange fields are in the hundreds of Tesla, for the equilibrium separation non reciprocal effects are found to be weak compared to those for a comparable external magnetic bias. For non equilibrium separations, nonreciprocal effects can be substantial.


title: Generalization Study of Quantum Neural Network

abstract: Generalization is an important feature of neural network, and there have been many studies on it. Recently, with the development of quantum compu-ting, it brings new opportunities. In this paper, we studied a class of quantum neural network constructed by quantum gate. In this model, we mapped the feature data to a quantum state in Hilbert space firstly, and then implement unitary evolution on it, in the end, we can get the classification result by im-plement measurement on the quantum state. Since all the operations in quan-tum neural networks are unitary, the parameters constitute a hypersphere of Hilbert space. Compared with traditional neural network, the parameter space is flatter. Therefore, it is not easy to fall into local optimum, which means the quantum neural networks have better generalization. In order to validate our proposal, we evaluated our model on three public datasets, the results demonstrated that our model has better generalization than the classical neu-ral network with the same structure.


title: Approximation on compact sets of functions and all derivatives

abstract: In Mergelyan type approximation we uniformly approximate functions on compact sets K by polynomials or rational functions or holomorphic functions on varying open sets containing K. In the present paper we consider analogous approximation, where uniform convergence on K is replaced by uniform approximation on K of all order derivatives.


title: Analytical solution of diffusion probability for a flat potential with a gaussian sink

abstract: We give a very simple method for finding the exact analytical solution for the problem of a particle undergoing diffusive motion on a flat potential in the presence of a gaussian sink function. The diffusion process is modelled by using one dimensional Smoluchowski equation. Our method provides solution in Laplace domain, which is used to derive an analytical expression for time average rate constant. Our solution can be used to analyze several related problems involving diffusion-reaction systems.


title: Analytical solutions of the simple shear problem for certain types of micromorphic continuum models -- including full derivations

abstract: To draw conclusions as regards the stability and modelling limits of the investigated continuum, we consider a family of infinitesimal isotropic generalized continuum models (Mindlin-Eringen micromorphic, relaxed micromorphic continuum, Cosserat, micropolar, microstretch, microstrain, microvoid, indeterminate couple stress, second gradient elasticity, etc.) and solve analytically the simple shear problem of an infinite stripe. A qualitative measure characterizing the different generalized continuum moduli is given by the shear stiffness $\mu^{*}$. This stiffness is in general length-scale dependent. Interesting limit cases are highlighted, which allow to interpret some of the appearing material parameter of the investigated continua.


title: Data-driven learning of non-autonomous systems

abstract: We present a numerical framework for recovering unknown non-autonomous dynamical systems with time-dependent inputs. To circumvent the difficulty presented by the non-autonomous nature of the system, our method transforms the solution state into piecewise integration of the system over a discrete set of time instances. The time-dependent inputs are then locally parameterized by using a proper model, for example, polynomial regression, in the pieces determined by the time instances. This transforms the original system into a piecewise parametric system that is locally time invariant. We then design a deep neural network structure to learn the local models. Once the network model is constructed, it can be iteratively used over time to conduct global system prediction. We provide theoretical analysis of our algorithm and present a number of numerical examples to demonstrate the effectiveness of the method.

categories: eess.SP cs.LG cs.NA math.DS math.NA physics.comp-ph


title: M-regular Decompositions for pushforwards of pluricanonical bundles of pairs to abelian varieties

abstract: We extend the so called Chen-Jiang decomposition for pushforward of pluricanocanical bundles to abelian varieties to the setting of klt pairs. We also provide a geometric application of this decomposition.


title: Analytic Baby Skyrmions at Finite Density

abstract: We introduce a consistent ansatz for the baby Skyrme model in (2+1)-dimensions which is able to reduce the complete set of field equations to just one equation for the profile function in situations in which the baby baryon charge can be arbitrary. Many analytic solutions both with and without the inclusion of the effects of the minimal coupling with the Maxwell field are constructed. Linear stability and other physical properties are discussed. These analytic gauged baby Skyrmions generate a persistent $U(1)$ current which cannot be turned off continuously as it is tied to the topological charge of the baby Skyrmions themselves. In the simplest non-trivial case of a gauged baby Skyrmion, a very important role is played by the Mathieu equation with an effective coupling constant which can be computed explicitly. These configurations are a very suitable arena to test resurgence in a non-integrable context.


title: On secant defectiveness and identifiability of Segre-Veronese varieties

abstract: We give an almost asymptotically sharp bound for the non secant defectiveness and identifiability of Segre-Veronese varieties. We also provide new examples of defective Segre-Veronese varieties.


title: How initial distribution affects symmetry breaking induced by panic in ants: experiment and flee-pheromone model

abstract: Collective escaping is a ubiquitous phenomenon in animal groups. Symmetry breaking caused by panic escape exhibits a shared feature across species that one exit is used more than the other when agents escaping from a closed space with two symmetrically located exists. Intuitively, one exit will be used more by more individuals close to it, namely there is an asymmetric distribution initially. We used ant groups to investigate how initial distribution of colonies would influence symmetry breaking in collective escaping. Surprisingly, there was no positive correlation between symmetry breaking and the asymmetrically initial distribution, which was quite counter-intuitive. In the experiments, a flee stage was observed and accordingly a flee-pheromone model was introduced to depict this special behavior in the early stage of escaping. Simulation results fitted well with the experiment. Furthermore, the flee stage duration was calibrated quantitatively and the model reproduced the observation demonstrated by our previous work. This paper explicitly distinguished two stages in ant panic escaping for the first time, thus enhancing the understanding in escaping behavior of ant colonies.


title: One Step to Efficient Synthetic Data

abstract: A common approach to synthetic data is to sample from a fitted model. We show that under general assumptions, this approach results in a sample with inefficient estimators and whose joint distribution is inconsistent with the true distribution. Motivated by this, we propose a general method of producing synthetic data, which is widely applicable for parametric models, has asymptotically efficient summary statistics, and is both easily implemented and highly computationally efficient. Our approach allows for the construction of both partially synthetic datasets, which preserve certain summary statistics, as well as fully synthetic data which satisfy the strong guarantee of differential privacy (DP), both with the same asymptotic guarantees. We also provide theoretical and empirical evidence that the distribution from our procedure converges to the true distribution. Besides our focus on synthetic data, our procedure can also be used to perform approximate hypothesis tests in the presence of intractable likelihood functions.


title: SQUIRREL: Testing Database Management Systems with Language Validity and Coverage Feedback

abstract: Fuzzing is an increasingly popular technique for verifying software functionalities and finding security vulnerabilities. However, current mutation-based fuzzers cannot effectively test database management systems (DBMSs), which strictly check inputs for valid syntax and semantics. Generation-based testing can guarantee the syntax correctness of the inputs, but it does not utilize any feedback, like code coverage, to guide the path exploration. In this paper, we develop Squirrel, a novel fuzzing framework that considers both language validity and coverage feedback to test DBMSs. We design an intermediate representation (IR) to maintain SQL queries in a structural and informative manner. To generate syntactically correct queries, we perform type-based mutations on IR, including statement insertion, deletion and replacement. To mitigate semantic errors, we analyze each IR to identify the logical dependencies between arguments, and generate queries that satisfy these dependencies. We evaluated Squirrel on four popular DBMSs: SQLite, MySQL, PostgreSQL and MariaDB. Squirrel found 51 bugs in SQLite, 7 in MySQL and 5 in MariaDB. 52 of the bugs are fixed with 12 CVEs assigned. In our experiment, Squirrel achieves 2.4x-243.9x higher semantic correctness than state-of-the-art fuzzers, and explores 2.0x-10.9x more new edges than mutation-based tools. These results show that Squirrel is effective in finding memory errors of database management systems.


title: ExKMC: Expanding Explainable $k$-Means Clustering

abstract: Despite the popularity of explainable AI, there is limited work on effective methods for unsupervised learning. We study algorithms for $k$-means clustering, focusing on a trade-off between explainability and accuracy. Following prior work, we use a small decision tree to partition a dataset into $k$ clusters. This enables us to explain each cluster assignment by a short sequence of single-feature thresholds. While larger trees produce more accurate clusterings, they also require more complex explanations. To allow flexibility, we develop a new explainable $k$-means clustering algorithm, ExKMC, that takes an additional parameter $k' \geq k$ and outputs a decision tree with $k'$ leaves. We use a new surrogate cost to efficiently expand the tree and to label the leaves with one of $k$ clusters. We prove that as $k'$ increases, the surrogate cost is non-increasing, and hence, we trade explainability for accuracy. Empirically, we validate that ExKMC produces a low cost clustering, outperforming both standard decision tree methods and other algorithms for explainable clustering. Implementation of ExKMC available at https://github.com/navefr/ExKMC.


title: Nonequilibrium phase transition in transport through a driven quantum point contact

abstract: We study transport of noninteracting fermions through a periodically driven quantum point contact (QPC) connecting two tight-binding chains. Initially, each chain is prepared in its own equilibrium state, generally with a bias in chemical potentials and temperatures. We examine the heating rate (or, alternatively, energy increase per cycle) in the nonequilibrium time-periodic steady state established after initial transient dynamics. We find that the heating rate vanishes identically when the driving frequency exceeds the bandwidth of the chain. We first establish this fact for a particular type of QPC where the heating rate can be calculated analytically. Then we verify numerically that this nonequilibrium phase transition is present for a generic QPC. Finally, we derive this effect perturbatively in leading order for cases when the QPC Hamiltonian can be considered as a small perturbation. Strikingly, we discover that for certain QPCs the current averaged over the driving cycle also vanishes above the critical frequency, despite a persistent bias. This shows that a driven QPC can act as a frequency-controlled quantum switch.


title: Irrational-slope versions of Thompson's groups $T$ and $V$

abstract: In this paper we consider the $T$- and $V$- versions, $T_{\tau}$ and $V_{\tau}$ , of the irrational slope Thompson group $F_{\tau}$ considered in [3]. We give infinite presentations for these groups and show how they can be represented by tree-pair diagrams similar to those for $T$ and $V$. We also show that $T_{\tau}$ and $V_{\tau}$ have index-2 normal subgroups, unlike their original Thompson counterparts $T$ and $V$. These index-2 subgroups are shown to be simple.


title: Learning Memory-Based Control for Human-Scale Bipedal Locomotion

abstract: Controlling a non-statically stable biped is a difficult problem largely due to the complex hybrid dynamics involved. Recent work has demonstrated the effectiveness of reinforcement learning (RL) for simulation-based training of neural network controllers that successfully transfer to real bipeds. The existing work, however, has primarily used simple memoryless network architectures, even though more sophisticated architectures, such as those including memory, often yield superior performance in other RL domains. In this work, we consider recurrent neural networks (RNNs) for sim-to-real biped locomotion, allowing for policies that learn to use internal memory to model important physical properties. We show that while RNNs are able to significantly outperform memoryless policies in simulation, they do not exhibit superior behavior on the real biped due to overfitting to the simulation physics unless trained using dynamics randomization to prevent overfitting; this leads to consistently better sim-to-real transfer. We also show that RNNs could use their learned memory states to perform online system identification by encoding parameters of the dynamics into memory.


title: Comparison of motivic Chern classes and stable envelopes for cotangent bundles

abstract: We consider a complex smooth projective variety equipped with an action of an algebraic torus with a finite number of fixed points. We compare the motivic Chern classes of Bia{\l}ynicki-Birula cells with the $K$-theoretic stable envelopes of cotangent bundle. We prove that under certain geometric assumptions satisfied e.g. by homogenous spaces these two notions coincide up to normalization.


title: Towards Kac - van de Leur conjecture: locality of superconformal algebras

abstract: We prove locality of superconformal algebras: every pluperfect superconformal algebra is spanned by coefficients of a finite family of mutually local distributions. We also introduce quasi-Poisson algebras and show that they can be used to construct all known simple superconformal algebras.


title: Universal Thermodynamics in the Kitaev Fractional Liquid

abstract: In the Kitaev honeycomb model, the quantum spin fractionalizes into itinerant Majorana and gauge flux spontaneously upon cooling, leading to rich experimental ramifications at finite temperature and an upsurge of research interest. In this work, we employ the exponential tensor renormalization group approach to explore the Kitaev model under various perturbations, including the external fields, Heisenberg, and the off-diagonal couplings that are common in the Kitaev materials. Through large-scale manybody calculations, we find a Kitaev fractional liquid at intermediate temperature that is robust against perturbations. The fractional liquid exhibits universal thermodynamic behaviors, including the fractional thermal entropy, metallic specific heat, and an intermediate-temperature Curie law of magnetic susceptibility. The emergent universal susceptibility behavior, with a modified Curie constant, can be ascribed to the strongly fluctuating $\mathbb{Z}_2$ fluxes as well as the extremely short-ranged and bond-directional spin correlations. With this insight, we revisit the susceptibility measurements of Na$_2$IrO$_3$ and $\alpha$-RuCl$_3$, and find evident signatures of finite-temperature fractionalization and ferromagnetic Kitaev couplings. Moreover, the peculiar spin correlation in the fractional liquid corresponds to a stripy structure factor which rotates in the extended Brillouin zone as the spin component changes. Therefore, our findings encourage future experimental exploration of fractional liquid in the Kitaev materials by thermodynamic measurements and spin-resolved structure factor probes.


title: Emergence of quantum critical charge and spin-state fluctuations near the pressure-induced Mott transition in MnO, FeO, CoO, and NiO

abstract: We perform a comprehensive theoretical study of the pressure-induced evolution of the electronic structure, magnetic state, and phase stability of the late transition metal monoxides MnO, FeO, CoO, and NiO using a fully charge self-consistent DFT+dynamical mean-field theory method. Our results reveal that the pressure-induced Mott insulator-to-metal phase transition in MnO-NiO is accompanied by a simultaneous collapse of local magnetic moments and lattice volume, implying a complex interplay between chemical bonding and electronic correlations. We compute the pressure-induced evolution of relative weights of the different valence states and spin-state configurations. Employing the concept of fluctuating valence in a correlated solid, we demonstrate that in MnO, FeO, and CoO a Mott insulator-metal transition and collapse of the local moments is accompanied by a sharp crossover of the spin-state and valence configurations. Our microscopic explanation of the magnetic collapse differs from the accepted picture and points out a remarkable dynamical coexistence (frustration) of the high-, intermediate-, and low-spin states. In particular, in MnO, the magnetic collapse is found to be driven by the appearance of the intermediate-spin state (IS), competing with the low-spin (LS) state; in FeO, we observe a conventional high-spin to low-spin (HS-LS) crossover. Most interestingly, in CoO, we obtain a remarkable (dynamical) coexistence of the HS and LS states, i.e., a HS-LS frustration, up to high pressure. Our results demonstrate the importance of quantum fluctuations of the valence and spin states for the understanding of quantum criticality of the Mott transitions.


title: Exotic t-structure for partial resolutions of the nilpotent cone

abstract: We define and study an exotic t-structure on the bounded derived category of equivariant coherent sheaves on partial resolutions of the nilpotent cone.


title: Dynamic Longest Common Substring in Polylogarithmic Time

abstract: The longest common substring problem consists in finding a longest string that appears as a (contiguous) substring of two input strings. We consider the dynamic variant of this problem, in which we are to maintain two dynamic strings $S$ and $T$, each of length at most $n$, that undergo substitutions of letters, in order to be able to return a longest common substring after each substitution. Recently, Amir et al. [ESA 2019] presented a solution for this problem that needs only $\tilde{\mathcal{O}}(n^{2/3})$ time per update. This brought the challenge of determining whether there exists a faster solution with polylogarithmic update time, or (as is the case for other dynamic problems), we should expect a polynomial (conditional) lower bound. We answer this question by designing a significantly faster algorithm that processes each substitution in amortized $\log^{\mathcal{O}(1)} n$ time with high probability. Our solution relies on exploiting the local consistency of the parsing of a collection of dynamic strings due to Gawrychowski et al. [SODA 2018], and on maintaining two dynamic trees with labeled bicolored leaves, so that after each update we can report a pair of nodes, one from each tree, of maximum combined weight, which have at least one common leaf-descendant of each color. We complement this with a lower bound of $\Omega(\log n/ \log\log n)$ for the update time of any polynomial-size data structure that maintains the LCS of two dynamic strings, and the same lower bound for the update time of any data structure of size $\tilde{\mathcal{O}}(n)$ that maintains the LCS of a static and a dynamic string. Both lower bounds hold even allowing amortization and randomization.


title: On the Promise of the Stochastic Generalized Gauss-Newton Method for Training DNNs

abstract: Following early work on Hessian-free methods for deep learning, we study a stochastic generalized Gauss-Newton method (SGN) for training DNNs. SGN is a second-order optimization method, with efficient iterations, that we demonstrate to often require substantially fewer iterations than standard SGD to converge. As the name suggests, SGN uses a Gauss-Newton approximation for the Hessian matrix, and, in order to compute an approximate search direction, relies on the conjugate gradient method combined with forward and reverse automatic differentiation. Despite the success of SGD and its first-order variants, and despite Hessian-free methods based on the Gauss-Newton Hessian approximation having been already theoretically proposed as practical methods for training DNNs, we believe that SGN has a lot of undiscovered and yet not fully displayed potential in big mini-batch scenarios. For this setting, we demonstrate that SGN does not only substantially improve over SGD in terms of the number of iterations, but also in terms of runtime. This is made possible by an efficient, easy-to-use and flexible implementation of SGN we propose in the Theano deep learning platform, which, unlike Tensorflow and Pytorch, supports forward automatic differentiation. This enables researchers to further study and improve this promising optimization technique and hopefully reconsider stochastic second-order methods as competitive optimization techniques for training DNNs; we also hope that the promise of SGN may lead to forward automatic differentiation being added to Tensorflow or Pytorch. Our results also show that in big mini-batch scenarios SGN is more robust than SGD with respect to its hyperparameters (we never had to tune its step-size for our benchmarks!), which eases the expensive process of hyperparameter tuning that is instead crucial for the performance of first-order methods.


title: On a $p(\cdot)$-biharmonic problem of Kirchhoff type involving critical growth

abstract: We establish a concentration-compactness principle for the Sobolev space $W^{2,p(\cdot)}(\Omega)\cap W_0^{1,p(\cdot)}(\Omega)$ that is a tool for overcoming the lack of compactness of the critical Sobolev imbedding. Using this result we obtain several existence and multiplicity results for a class of Kirchhoff type problems involving $p(\cdot)$-biharmonic operator and critical growth.


title: Explaining Multi-stage Tasks by Learning Temporal Logic Formulas from Suboptimal Demonstrations

abstract: We present a method for learning multi-stage tasks from demonstrations by learning the logical structure and atomic propositions of a consistent linear temporal logic (LTL) formula. The learner is given successful but potentially suboptimal demonstrations, where the demonstrator is optimizing a cost function while satisfying the LTL formula, and the cost function is uncertain to the learner. Our algorithm uses the Karush-Kuhn-Tucker (KKT) optimality conditions of the demonstrations together with a counterexample-guided falsification strategy to learn the atomic proposition parameters and logical structure of the LTL formula, respectively. We provide theoretical guarantees on the conservativeness of the recovered atomic proposition sets, as well as completeness in the search for finding an LTL formula consistent with the demonstrations. We evaluate our method on high-dimensional nonlinear systems by learning LTL formulas explaining multi-stage tasks on 7-DOF arm and quadrotor systems and show that it outperforms competing methods for learning LTL formulas from positive examples.


title: Hausdorff measure and Assouad dimension of generic self-conformal IFS on the line

abstract: This paper considers self-conformal iterated function systems (IFSs) on the real line whose first level cylinders overlap. In the space of self-conformal IFSs, we show that generically (in topological sense) if the attractor of such a system has Hausdorff dimension less than $1$ then it has zero appropriate dimensional Hausdorff measure and its Assouad dimension is equal to $1$. Our main contribution is in showing that if the cylinders intersect then the IFS generically does not satisfy the weak separation property and hence, we may apply a recent result of Angelevska, K\"aenm\"aki and Troscheit [BLMS, 2020]. This phenomenon holds for transversal families (in particular for the translation family) typically, in the self-similar case, in both topological and in measure theoretical sense, and in the more general self-conformal case in the topological sense.


title: DGSAC: Density Guided Sampling and Consensus

abstract: Robust multiple model fitting plays a crucial role in many computer vision applications. Unlike single model fitting problems, the multi-model fitting has additional challenges. The unknown number of models and the inlier noise scale are the two most important of them, which are in general provided by the user using ground-truth or some other auxiliary information. Mode seeking/ clustering-based approaches crucially depend on the quality of model hypotheses generated. While preference analysis based guided sampling approaches have shown remarkable performance, they operate in a time budget framework, and the user provides the time as a reasonable guess. In this paper, we deviate from the mode seeking and time budget framework. We propose a concept called Kernel Residual Density (KRD) and apply it to various components of a multiple-model fitting pipeline. The Kernel Residual Density act as a key differentiator between inliers and outliers. We use KRD to guide and automatically stop the sampling process. The sampling process stops after generating a set of hypotheses that can explain all the data points. An explanation score is maintained for each data point, which is updated on-the-fly. We propose two model selection algorithms, an optimal quadratic program based, and a greedy. Unlike mode seeking approaches, our model selection algorithms seek to find one representative hypothesis for each genuine structure present in the data. We evaluate our method (dubbed as DGSAC) on a wide variety of tasks like planar segmentation, motion segmentation, vanishing point estimation, plane fitting to 3D point cloud, line, and circle fitting, which shows the effectiveness of our method and its unified nature.


title: On countable elementary free groups

abstract: We prove that if a countable group is elementarily equivalent to a non-abelian free group and all of its abelian subgroups are cyclic, then the group is a union of a chain of regular NTQ groups (i.e., hyperbolic towers).


title: Anatomical Mesh-Based Virtual Fixtures for Surgical Robots

abstract: This paper presents a dynamic constraint formulation to provide protective virtual fixtures of 3D anatomical structures from polygon mesh representations. The proposed approach can anisotropically limit the tool motion of surgical robots without any assumption of the local anatomical shape close to the tool. Using a bounded search strategy and Principle Directed tree, the proposed system can run efficiently at 180 Hz for a mesh object containing 989,376 triangles and 493,460 vertices. The proposed algorithm has been validated in both simulation and skull cutting experiments. The skull cutting experiment setup uses a novel piezoelectric bone cutting tool designed for the da Vinci research kit. The result shows that the virtual fixture assisted teleoperation has statistically significant improvements in the cutting path accuracy and penetration depth control. The code has been made publicly available at https://github.com/mli0603/PolygonMeshVirtualFixture.


title: Assessing Holistic Impacts of Major Events on the Bitcoin Blockchain Network

abstract: As the pioneer of blockchain technology, Bitcoin is the most popular cryptocurrency to date. Given its dramatic price spikes (and crashes) along with the never-ending news from SEC regulations to security breaches, there seems to be a lack of understanding about the dynamics of cryptocurrencies. These dynamics are believed to be affected by various political, security, financial, and regulatory events. In this paper, we present an efficient framework for holistic analysis of cryptocurrency fluctuations by introducing the Impact-Score metric to distinguish event-induced changes from normal variations. We have applied our framework to 16 major worldwide events and the Bitcoin blockchain network (defined as Bitcoin transaction and users, blockchain data, and memory pool data) from 2016-2018. The results show that a majority of the events are correlated with substantial network changes. We observed roughly generalizable correlations between event types (e.g. financial events) and sub-structures of the Bitcoin blockchain network. Subgroups of these events have strongly consistent temporal impacts on specific facets (e.g. activity or fees) of the Bitcoin ecosystem. Furthermore, we demonstrate the robustness of our process by correlating a majority of spikes in network/subnetwork change with major events.


title: Curry-Howard-Lambek Correspondence for Intuitionistic Belief

abstract: This paper introduces a natural deduction calculus for intuitionistic logic of belief $\mathsf{IEL}^{-}$ which is easily turned into a modal $\lambda$-calculus giving a computational semantics for deductions in $\mathsf{IEL}^{-}$. By using that interpretation, it is also proved that $\mathsf{IEL}^{-}$ has good proof-theoretic properties. The correspondence between deductions and typed terms is then extended to a categorical semantics for identity of proofs in $\mathsf{IEL}^{-}$ showing the general structure of such a modality for belief in an intuitionistic framework.


title: Intergalactic filaments spin

abstract: Matter in the Universe is arranged in a cosmic web, with a filament of matter typically connecting each neighbouring galaxy pair, separated by tens of millions of light-years. A quadrupolar pattern of the spin field around filaments is known to influence the spins of galaxies and haloes near them, but it remains unknown whether filaments themselves spin. Here, we measure dark-matter velocities around filaments in cosmological simulations, finding that matter generally rotates around them, much faster than around a randomly located axis. It also exhibits some coherence along the filament. The net rotational component is comparable to, and often dominant over, the known quadrupolar flow. The evidence of net rotations revises previous emphasis on a quadrupolar spin field around filaments. The full picture of rotation in the cosmic web is more complicated and multiscale than a network of spinning filamentary rods, but we argue that filament rotation is substantial enough to be an essential part of the picture. It is likely that the longest coherently rotating objects in the Universe are filaments. Also, we speculate that this rotation could provide a mechanism to generate or amplify intergalactic magnetic fields in filaments.


title: Emergent Multi-Agent Communication in the Deep Learning Era

abstract: The ability to cooperate through language is a defining feature of humans. As the perceptual, motory and planning capabilities of deep artificial networks increase, researchers are studying whether they also can develop a shared language to interact. From a scientific perspective, understanding the conditions under which language evolves in communities of deep agents and its emergent features can shed light on human language evolution. From an applied perspective, endowing deep networks with the ability to solve problems interactively by communicating with each other and with us should make them more flexible and useful in everyday life. This article surveys representative recent language emergence studies from both of these two angles.


title: Learning to Scan: A Deep Reinforcement Learning Approach for Personalized Scanning in CT Imaging

abstract: Computed Tomography (CT) takes X-ray measurements on the subjects to reconstruct tomographic images. As X-ray is radioactive, it is desirable to control the total amount of dose of X-ray for safety concerns. Therefore, we can only select a limited number of measurement angles and assign each of them limited amount of dose. Traditional methods such as compressed sensing usually randomly select the angles and equally distribute the allowed dose on them. In most CT reconstruction models, the emphasize is on designing effective image representations, while much less emphasize is on improving the scanning strategy. The simple scanning strategy of random angle selection and equal dose distribution performs well in general, but they may not be ideal for each individual subject. It is more desirable to design a personalized scanning strategy for each subject to obtain better reconstruction result. In this paper, we propose to use Reinforcement Learning (RL) to learn a personalized scanning policy to select the angles and the dose at each chosen angle for each individual subject. We first formulate the CT scanning process as an MDP, and then use modern deep RL methods to solve it. The learned personalized scanning strategy not only leads to better reconstruction results, but also shows strong generalization to be combined with different reconstruction algorithms.


title: The cardinality of the sublattice of closed ideals of operators between certain classical sequence spaces

abstract: Theorem A and Theorem B of [1] state that for $1<p<\infty$ the lattice of closed ideals of $\mathcal{L}(\ell_p,c_0)$, $\mathcal{L}(\ell_p,\ell_\infty)$ and of $\mathcal{L}(\ell_1,\ell_p)$ are at least of cardinality $2^{\omega}$. Here we show that the cardinality of the lattice of closed ideals of $\mathcal{L}(\ell_p,c_0)$, $\mathcal{L}(\ell_p,\ell_\infty)$ and of $\mathcal{L}(\ell_1,\ell_p)$, is at least $2^{2^\omega}$, and thus equal to it.


title: The Grassmannian VOA

abstract: We study the 3-parametric family of vertex operator algebras based on the unitary Grassmannian coset CFT $\mathfrak{u}(M+N)_k/(\mathfrak{u}(M)_k \times \mathfrak{u}(N)_k)$. This VOA serves as a basic building block for a large class of cosets and generalizes the $\mathcal{W}_\infty$ algebra. We analyze representations and their characters in detail and find surprisingly simple character formulas for the representations in the generic parameter regime that admit an elegant combinatorial formulation. We also discuss truncations of the algebra and give a conjectural formula for the complete set of truncation curves. We develop a theory of gluing for these algebras in order to build more complicated coset and non-coset algebras. We demonstrate the power of this technology with some examples and show in particular that the $\mathcal{N}=2$ supersymmetric Grassmannian can be obtained by gluing three bosonic Grassmannian algebras in a loop. We finally speculate about the tantalizing possibility that this algebra is a specialization of an even larger 4-parametric family of algebras exhibiting pentality symmetry. Specialization of this conjectural family should include both the unitary Grassmannian family as well as the Lagrangian Grassmannian family of VOAs which interpolates between the unitary and the orthosymplectic cosets.


title: A Negative Correlation Strategy for Bracketing in Difference-in-Differences

abstract: The method of difference-in-differences (DID) is widely used to study the causal effect of policy interventions in observational studies. DID employs a before and after comparison of the treated and control units to remove bias due to time-invariant unmeasured confounders under the parallel trends assumption. Estimates from DID, however, will be biased if the outcomes for the treated and control units evolve differently in the absence of treatment, namely if the parallel trends assumption is violated. We propose a general identification strategy that leverages two groups of control units whose outcomes relative to the treated units exhibit a negative correlation, and achieves partial identification of the average treatment effect for the treated. The identified set is of a union bounds form that involves the minimum and maximum operators, which makes the canonical bootstrap generally inconsistent and naive methods overly conservative. By utilizing the directional inconsistency of the bootstrap distribution, we develop a novel bootstrap method to construct uniformly valid confidence intervals for the identified set and parameter of interest when the identified set is of a union bounds form, and we establish the method's theoretical properties. We develop a simple falsification test and sensitivity analysis. We apply the proposed strategy for bracketing to study whether minimum wage laws affect employment levels.


title: Quantum process tomography with unsupervised learning and tensor networks

abstract: The impressive pace of advance of quantum technology calls for robust and scalable techniques for the characterization and validation of quantum hardware. Quantum process tomography, the reconstruction of an unknown quantum channel from measurement data, remains the quintessential primitive to completely characterize quantum devices. However, due to the exponential scaling of the required data and classical post-processing, its range of applicability is typically restricted to one- and two-qubit gates. Here, we present a new technique for performing quantum process tomography that addresses these issues by combining a tensor network representation of the channel with a data-driven optimization inspired by unsupervised machine learning. We demonstrate our technique through synthetically generated data for ideal one- and two-dimensional random quantum circuits of up to 10 qubits, and a noisy 5-qubit circuit, reaching process fidelities above 0.99 using only a limited set of single-qubit measurement samples and input states. Our results go far beyond state-of-the-art, providing a practical and timely tool for benchmarking quantum circuits in current and near-term quantum computers.


title: Equivariant Flows: Exact Likelihood Generative Learning for Symmetric Densities

abstract: Normalizing flows are exact-likelihood generative neural networks which approximately transform samples from a simple prior distribution to samples of the probability distribution of interest. Recent work showed that such generative models can be utilized in statistical mechanics to sample equilibrium states of many-body systems in physics and chemistry. To scale and generalize these results, it is essential that the natural symmetries in the probability density -- in physics defined by the invariances of the target potential -- are built into the flow. We provide a theoretical sufficient criterion showing that the distribution generated by \textit{equivariant} normalizing flows is invariant with respect to these symmetries by design. Furthermore, we propose building blocks for flows which preserve symmetries which are usually found in physical/chemical many-body particle systems. Using benchmark systems motivated from molecular physics, we demonstrate that those symmetry preserving flows can provide better generalization capabilities and sampling efficiency.

categories: stat.ML cs.LG physics.chem-ph physics.comp-ph


title: Models of zero-range interaction for the bosonic trimer at unitarity

abstract: We present the mathematical construction of the physically relevant quantum Hamiltonians for a three-body systems consisting of identical bosons mutually coupled by a two-body interaction of zero range. For a large part of the presentation, infinite scattering length will be considered (the unitarity regime). The subject has several precursors in the mathematical literature. We proceed through an operator-theoretic construction of the self-adjoint extensions of the minimal operator obtained by restricting the free Hamiltonian to wave-functions that vanish in the vicinity of the coincidence hyperplanes: all extensions thus model an interaction precisely supported at the spatial configurations where particles come on top of each other. Among them, we select the physically relevant ones, by implementing in the operator construction the presence of the specific short-scale structure suggested by formal physical arguments that are ubiquitous in the physical literature on zero-range methods. This is done by applying at different stages the self-adjoint extension schemes a la Kre{\u\i}n-Vi\v{s}ik-Birman and a la von Neumann. We produce a class of canonical models for which we also analyse the structure of the negative bound states. Bosonicity and zero range combined together make such canonical models display the typical Thomas and Efimov spectra, i.e., sequence of energy eigenvalues accumulating to both minus infinity and zero. We also discuss a type of regularisation that prevents such spectral instability while retaining an effective short-scale pattern. Beside the operator qualification, we also present the associated energy quadratic forms. We structured our analysis so as to clarify certain steps of the operator-theoretic construction that are notoriously subtle for the correct identification of a domain of self-adjointness.


title: Lyapunov spectra of chaotic recurrent neural networks

abstract: Brains process information through the collective dynamics of large neural networks. Collective chaos was suggested to underlie the complex ongoing dynamics observed in cerebral cortical circuits and determine the impact and processing of incoming information streams. In dissipative systems, chaotic dynamics takes place on a subset of phase space of reduced dimensionality and is organized by a complex tangle of stable, neutral and unstable manifolds. Key topological invariants of this phase space structure such as attractor dimension, and Kolmogorov-Sinai entropy so far remained elusive. Here we calculate the complete Lyapunov spectrum of recurrent neural networks. We show that chaos in these networks is extensive with a size-invariant Lyapunov spectrum and characterized by attractor dimensions much smaller than the number of phase space dimensions. We find that near the onset of chaos, for very intense chaos, and discrete-time dynamics, random matrix theory provides analytical approximations to the full Lyapunov spectrum. We show that a generalized time-reversal symmetry of the network dynamics induces a point-symmetry of the Lyapunov spectrum reminiscent of the symplectic structure of chaotic Hamiltonian systems. Fluctuating input reduces both the entropy rate and the attractor dimension. For trained recurrent networks, we find that Lyapunov spectrum analysis provides a quantification of error propagation and stability achieved. Our methods apply to systems of arbitrary connectivity, and we describe a comprehensive set of controls for the accuracy and convergence of Lyapunov exponents. Our results open a novel avenue for characterizing the complex dynamics of recurrent neural networks and the geometry of the corresponding chaotic attractors. They also highlight the potential of Lyapunov spectrum analysis as a diagnostic for machine learning applications of recurrent networks.


title: Thermodynamic perspective on the field-induced behavior of $\alpha$-RuCl$_3$

abstract: Measurements of the magnetic Gr\"uneisen parameter ($\Gamma_B$) and specific heat on the Kitaev material candidate $\alpha$-RuCl$_3$ are used to access in-plane field- and temperature-dependence of the entropy up to 12 T and down to 1 K. No signatures corresponding to phase transitions are detected beyond the boundary of the magnetically ordered region, but only a shoulder-like anomaly in $\Gamma_B$, involving an entropy increment as small as $10^{-5} R\log 2$. These observations put into question the presence of a thermodynamic phase transition between the purported quantum spin liquid and the field-polarized state of $\alpha$-RuCl$_3$. We show theoretically that at low temperatures $\Gamma_B$ is sensitive to crossings in the lowest excitations within gapped phases, and identify the measured shoulder-like anomaly as being of such origin. Exact diagonalization calculations demonstrate that the shoulder-like anomaly can be reproduced in extended Kitaev models that gain proximity to an additional phase at finite field without entering it. We discuss manifestations of this proximity in other measurements.


title: Infinite co-minimal pairs involving lacunary sequences and generalisations to higher dimensions

abstract: The study of minimal complements in a group or a semigroup was initiated by Nathanson. The notion of minimal complements and being a minimal complement leads to the notion of co-minimal pairs which was considered in a prior work of the authors. In this article, we study which type of subsets in the integers and free abelian groups of higher rank can be a part of a co-minimal pair. We show that a majority of lacunary sequences have this property. From the conditions established, one can show that any infinite subset of any finitely generated abelian group has uncountably many subsets which is a part of a co-minimal pair. Further, the uncountable collection of sets can be chosen so that they satisfy certain algebraic properties.


title: Virasoro hair and entropy for axisymmetric Killing horizons

abstract: We show that the gravitational phase space for the near-horizon region of a bifurcate, axisymmetric Killing horizon in any dimension admits a 2D conformal symmetry algebra with central charges proportional to the area. This extends the construction of [Haco et. al., JHEP 12, 098 (2018)] to generic Killing horizons appearing in solutions of Einstein's equations, and motivates a holographic description in terms of a 2D conformal field theory. The Cardy entropy in such a field theory agrees with the Bekenstein-Hawking entropy of the horizon, suggesting a microscopic interpretation.


title: Targeting SARS-CoV-2 with AI- and HPC-enabled Lead Generation: A First Data Release

abstract: Researchers across the globe are seeking to rapidly repurpose existing drugs or discover new drugs to counter the the novel coronavirus disease (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). One promising approach is to train machine learning (ML) and artificial intelligence (AI) tools to screen large numbers of small molecules. As a contribution to that effort, we are aggregating numerous small molecules from a variety of sources, using high-performance computing (HPC) to computer diverse properties of those molecules, using the computed properties to train ML/AI models, and then using the resulting models for screening. In this first data release, we make available 23 datasets collected from community sources representing over 4.2 B molecules enriched with pre-computed: 1) molecular fingerprints to aid similarity searches, 2) 2D images of molecules to enable exploration and application of image-based deep learning methods, and 3) 2D and 3D molecular descriptors to speed development of machine learning models. This data release encompasses structural information on the 4.2 B molecules and 60 TB of pre-computed data. Future releases will expand the data to include more detailed molecular simulations, computed models, and other products.


title: A unifying perspective on linear continuum equations prevalent in science. Part III: Canonical forms for dynamic equations with moduli that may, or may not, vary with time

abstract: Enlarging on Parts I and II we write more equations in the desired format of the extended abstract theory of composites. We focus on a multitude of full dynamic equations, including equations where the medium is moving or otherwise changing in time. The motivation is that results and methods in the theory of composites then extend to these equations.


title: A unifying perspective on linear continuum equations prevalent in science. Part II: Canonical forms for time-harmonic equations

abstract: Following some past advances, we reformulate a large class of linear continuum science equations in the format of the extended abstract theory of composites so that we can apply this theory to better understand and efficiently solve those equations. Here in part II we elucidate the form for many time-harmonic equations that do not involve higher order gradients.

categories: math-ph math.AP math.MP physics.class-ph


title: Visual Summarization of Lecture Video Segments for Enhanced Navigation

abstract: Lecture videos are an increasingly important learning resource for higher education. However, the challenge of quickly finding the content of interest in a lecture video is an important limitation of this format. This paper introduces visual summarization of lecture video segments to enhance navigation. A lecture video is divided into segments based on the frame-to-frame similarity of content. The user navigates the lecture video content by viewing a single frame visual and textual summary of each segment. The paper presents a novel methodology to generate the visual summary of a lecture video segment by computing similarities between images extracted from the segment and employing a graph-based algorithm to identify the subset of most representative images. The results from this research are integrated into a real-world lecture video management portal called Videopoints. To collect ground truth for evaluation, a survey was conducted where multiple users manually provided visual summaries for 40 lecture video segments. The users also stated whether any images were not selected for the summary because they were similar to other selected images. The graph based algorithm for identifying summary images achieves 78% precision and 72% F1-measure with frequently selected images as the ground truth, and 94% precision and 72% F1-measure with the union of all user selected images as the ground truth. For 98% of algorithm selected visual summary images, at least one user also selected that image for their summary or considered it similar to another image they selected. Over 65% of automatically generated summaries were rated as good or very good by the users on a 4-point scale from poor to very good. Overall, the results establish that the methodology introduced in this paper produces good quality visual summaries that are practically useful for lecture video navigation.


title: K3 metrics

abstract: We provide an explicit construction of Ricci-flat K3 metrics. It employs the technology of D-geometry, which in the case of interest is equivalent to a hyper-K\"ahler quotient. We relate it to the construction of arXiv:1810.10540, and in particular show that it contains the solution to the BPS state counting problem (that of computing the BPS index of a heterotic little string theory compactified on $T^2$) discussed therein, which is the data needed for this second construction of K3 metrics.


title: Revisiting Soltan's argument based on a semi-analytical model for galaxy and black hole evolution

abstract: We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions.


title: Transverse-Energy-Energy Correlations in Deep Inelastic Scattering

abstract: Event shape observables have been widely used for precision QCD studies at various lepton and hadron colliders. We present the most accurate calculation of the transverse-energy-energy correlation event shape variable in deep-inelastic scattering. In the framework of soft-collinear effective theory the cross section is factorized as the convolution of the hard function, beam function, jet function and soft function in the back-to-back limit. A close connection to TMD factorization is established, as the beam function when combined with part of the soft function is identical to the conventional TMD parton distribution function, and the jet function is the second moment of the TMD fragmentation function matching coefficient. We validate our framework by comparing the obtained LO and NLO leading singular distributions to the full QCD calculations in the back-to-back limit. We report the resummed transverse-energy-energy correlation distributions up to N$^3$LL accuracy matched with the NLO cross section for the production of a lepton and two jets. Our work provides a new way to precisely study TMD physics at the future Electron-Ion Collider.


title: Massive Islands

abstract: We comment on the role of the graviton mass in recent calculations of the Page curve using holographic ideas. All reliable calculations of the Page curve in more than 2+1 spacetime dimensions have been performed in systems with massive gravitons. A crucial ingredient in these calculations is the formation of islands, regions that contribute to the entropy of degrees of freedom located elsewhere. While most often simply ignored, it is indeed true that mass of the graviton does not appear to significantly affect the calculations that appeared in the literature. We use the freedom to change the graviton mass to give an extremely simple model of analytically tractable island formation in general dimensions. We do however note that if one attempts to take the limit of zero graviton mass, any contribution from the islands disappears. This raises the question to what extent entanglement islands can play a role in standard massless gravity.


title: Ongoing flyby in the young multiple system UX Tauri

abstract: We present observations of the young multiple system UX Tauri to look for circumstellar disks and for signs of dynamical interactions. We obtained SPHERE/IRDIS deep differential polarization images in the J and H bands. We also used ALMA archival CO data. Large extended spirals are well detected in scattered light coming out of the disk of UX Tau A. The southern spiral forms a bridge between UX Tau A and C. These spirals, including the bridge connecting the two stars, all have a CO (3-2) counterpart seen by ALMA. The disk of UX Tau C is detected in scattered light. It is much smaller than the disk of UX Tau A and has a major axis along a different position angle, suggesting a misalignment. We performed PHANTOM SPH hydrodynamical models to interpret the data. The scattered light spirals, CO emission spirals and velocity patterns of the rotating disks, and the compactness of the disk of UX Tau C all point to a scenario in which UX Tau A has been perturbed very recently (about 1000 years) by the close passage of UX Tau C.


title: Nonergodic Quantum Dynamics from Deformations of Classical Cellular Automata

abstract: Classical reversible cellular automata (CAs), which describe the discrete-time dynamics of classical degrees of freedom in a finite state-space, can exhibit exact, nonthermal quantum eigenstates despite being classically chaotic. We show that every classical CA defines a family of generically non-integrable, periodically-driven (Floquet) quantum dynamics with exact, nonthermal eigenstates. These Floquet dynamics are nonergodic in the sense that certain product states on a periodic classical orbit fail to thermalize, while generic initial states thermalize as expected in a quantum chaotic system. We demonstrate that some signatures of these effects can be probed in quantum simulators based on Rydberg atoms in the blockade regime. These results establish classical CAs as parent models for a class of quantum chaotic systems with rare nonthermal eigenstates.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph quant-ph


title: Detecting isolated stellar-mass black holes in the absence of microlensing parallax effect

abstract: Gravitational microlensing can detect isolated stellar-mass black holes (BHs), which are believed to be the dominant form of Galactic BHs according to population synthesis models. Previous searches for BH events in microlensing data focused on long-timescale events with significant microlensing parallax detections. Here we show that, although BH events preferentially have long timescales, the microlensing parallax amplitudes are so small that in most cases the parallax signals cannot be detected statistically significantly. We then identify OGLE-2006-BLG-044 to be a candidate BH event because of its long timescale and small microlensing parallax. Our findings have implications to future BH searches in microlensing data.


title: UV Freeze-in in Starobinsky Inflation

abstract: In the Starobinsky model of inflation, the observed dark matter abundance can be produced from the direct decay of the inflaton field only in a very narrow spectrum of close-to-conformal scalar fields and spinors of mass $\sim 10^7$ GeV. This spectrum can be, however, significantly broadened in the presence of effective non-renormalizable interactions between the dark and the visible sectors. In particular, we show that UV freeze-in can efficiently generate the right dark matter abundance for a large range of masses spanning from the keV to the PeV scale and arbitrary spin, without significantly altering the heating dynamics. We also consider the contribution of effective interactions to the inflaton decay into dark matter.


title: Luminosity Functions and Host-to-Host Scatter of Dwarf Satellite Systems in the Local Volume

abstract: Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation (SBF) distances to recently cataloged candidate dwarf satellites around 10 massive hosts within $D<12$ Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to $M_g{\sim}-9$ to $-10$, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting a SHMR that matches hydrodynamic simulations, the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function amongst LV hosts. Contrary to recent results, we find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, necessitating a higher normalization of the SHMR around halo masses of $10^{11}$ \msun\ than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy-subhalo connection in the low-mass regime.


title: Radial Distributions of Dwarf Satellite Systems in the Local Volume

abstract: The radial spatial distribution of low-mass satellites around a Milky Way (MW)-like host is an important benchmark for simulations of small-scale structure. The distribution is sensitive to the disruption of subhalos by the central disk and can indicate whether the disruption observed in simulations of MW analogs is artificial (i.e., numeric) or physical in origin. We consider a sample of 12 well-surveyed satellite systems of MW-like hosts in the Local Volume that are complete to $M_V<-9$ and within 150 projected kpc. We investigate the radial distribution of satellites and compare with $\Lambda$CDM cosmological simulations, including big-box cosmological simulations and high resolution zoom in simulations of MW sized halos. We find that the observed satellites are significantly more centrally concentrated than the simulated systems. Several of the observed hosts, including the MW, are $\sim2\sigma$ outliers relative to the simulated hosts in being too concentrated, while none of the observed hosts are less centrally concentrated than the simulations. This result is robust to different ways of measuring the radial concentration. We find that this discrepancy is more significant for bright, $M_V<-12$ satellites, suggestive that this is not the result of observational incompleteness. We argue that the discrepancy is possibly due to artificial disruption in the simulations, but, if so, this has important ramifications for what stellar to halo mass relation is allowed in the low-mass regime by the observed abundance of satellites.