title: Traversable Wormhole in Einstein 3-Form Theory With Self-Interacting Potential

abstract: We numerically construct a symmetric wormhole solution in pure Einstein gravity supported by a massive $3$-form field with a potential that contains a quartic self-interaction term. The wormhole spacetimes have only a single throat and they are everywhere regular and asymptotically flat. Furthermore, their mass and throat circumference increase almost linearly as the coefficient of the quartic self-interaction term $\Lambda$ increases. The amount of violation of the null energy condition (NEC) is proportional to the magnitude of $3$-form, thus the NEC is less violated as $\Lambda$ increases, since the magnitude of $3$-form decreases with $\Lambda$. In addition, we investigate the geodesics of particles moving around the wormhole. The unstable photon orbit is located at the throat. We also find that the wormhole can cast a shadow whose apparent size is smaller than that cast by the Schwarzschild black hole, but reduces to it when $\Lambda$ acquires a large value. The behavior of the innermost stable circular orbit around this wormhole is also discussed. The results of this paper hint toward the possibility that the 3-form wormholes could be potential black hole mimickers, as long as $\Lambda$ is sufficiently large, precisely when NEC is weakly violated.


title: Reimagining the water snowline

abstract: Water is a molecule that is tightly related to many facets of star and planet formation. Water's abundance and distribution, especially the location of it's snowline has thus been the subject of much study. While water is seen to be abundant in the inner region of proto-planetary disks in infrared spectroscopy, detections of water in the disk in the sub-millimeter are rare, with only one detection towards AS 205. Here we put the multitude of non-detections and the single detection into context of recent physico-chemical models. We find that the 321.2257 GHz (10(2,9) - 9(3,6)) line detection towards AS 205 is inconsistent with a normal inner disk temperature structure and that the observed line must be masing. Furthermore, the emitting area derived from the line width, together with published analyses on water in disks around T-Tauri stars implies that the water snowline in the disk surface is at the same location as the snowline in the mid-plane. We propose that this is caused by vertical mixing continuously sequestering water from the warm surface layers into the cold disk midplane.


title: The regularity of almost all edge ideals

abstract: A fruitful contemporary paradigm in graph theory is that almost all graphs that do not contain a certain subgraph have common structural characteristics. The "almost" is crucial, without it there is no structure. In this paper we transfer this paradigm to commutative algebra and make use of deep graph theoretic results. A key tool are the critical graphs introduced by Balogh and Butterfield. We consider edge ideals $I_G$ of graphs and their Betti numbers. The numbers of the form $\beta_{i,2i+2}$ constitute the "main diagonal" of the Betti table. It is well known that any Betti number $\beta_{i,j}(I_G)$ below (or equivalently, to the left of) this diagonal is always zero. We identify a certain "parabola" inside the Betti table and call parabolic Betti numbers the entries of the Betti table bounded on the left by the main diagonal and on the right by this parabola. Let $\beta_{i,j}$ be a parabolic Betti number on the $r$-th row of the Betti table, for $r\ge3$. Our main results state that almost all graphs $G$ with $\beta_{i,j}(I_G)=0$ can be partitioned into $r-2$ cliques and one independent set, and in particular for almost all graphs $G$ with $\beta_{i,j}(I_G)=0$ the regularity of $I_G$ is $r-1$.


title: Flavors of Magnetic Noise in Quantum Materials

abstract: The complexity of electronic band structures in quantum materials offers new charge-neutral degrees of freedom stable for transport, a promising example being the valley (axial) degree of freedom in Weyl semimetals (WSMs). A noninvasive probe of their transport properties is possible by exploiting the frequency dependence of the magnetic noise generated in the vicinity of the material. In this work, we investigate the magnetic noise generically associated with diffusive transport using a systematic Langevin approach. Taking a minimal model of magnetic WSMs for demonstration, we show that thermal fluctuations of the charge current, the valley current, and the magnetic order can give rise to magnetic noise with distinctively different spectral characters, which provide a theoretical guidance to separate their contributions. Our approach is extendable to the study of magnetic noise and its spectral features arising from other transport degrees of freedom in quantum materials.


title: When does the zero fiber of the moment map have rational singularities?

abstract: Let $G$ be a complex reductive group and $V$ a $G$-module. There is a natural moment mapping $\mu\colon V\oplus V^*\to\mathfrak{g}^*$ and we denote $\mu^{-1}(0)$ (the shell) by $N_V$. We use invariant theory and results of Musta\c{t}\u{a} [Mus01] to find criteria for $N_V$ to have rational singularities and for the categorical quotient $N_V /\!\!/ G$ to have symplectic singularities, the latter results improving upon [HSS20]. It turns out that for ``most'' $G$-modules $V$, the shell $N_V$ has rational singularities. For the case of direct sums of classical representations of the classical groups, $N_V$ has rational singularities and $N_V /\!\!/ G$ has symplectic singularities if $N_V$ is a reduced and irreducible complete intersection. Another important special case is $V=p\,\mathfrak{g}$ (the direct sum of $p$ copies of the Lie algebra of $G$) where $p\geq 2$. We show that $N_V$ has rational singularities and that $N_V /\!\!/ G$ has symplectic singularities, improving upon results of [Bud19], [AA16], [Kap19] and [GH20]. Let $\pi=\pi_1(\Sigma)$ where $\Sigma$ is a closed Riemann surface of genus $p\geq 2$. Let $G$ be semisimple and let $\operatorname{Hom}(\pi,G)$ and $\mathscr X\!(\pi,G)$ be the corresponding representation variety and character variety. We show that $\operatorname{Hom}(\pi,G)$ is a complete intersection with rational singularities and that $\mathscr X\!(\pi,G)$ has symplectic singularities. If $p>2$ or $G$ contains no simple factor of rank $1$, then the singularities of $\operatorname{Hom}(\pi,G)$ and $\mathscr X\!(\pi,G)$ are in codimension at least four and $\operatorname{Hom}(\pi,G)$ is locally factorial. If, in addition, $G$ is simply connected, then $\mathscr X\!(\pi,G)$ is locally factorial.


title: Synthesizing Pareto-Optimal Interpretations for Black-Box Models

abstract: We present a new multi-objective optimization approach for synthesizing interpretations that "explain" the behavior of black-box machine learning models. Constructing human-understandable interpretations for black-box models often requires balancing conflicting objectives. A simple interpretation may be easier to understand for humans while being less precise in its predictions vis-a-vis a complex interpretation. Existing methods for synthesizing interpretations use a single objective function and are often optimized for a single class of interpretations. In contrast, we provide a more general and multi-objective synthesis framework that allows users to choose (1) the class of syntactic templates from which an interpretation should be synthesized, and (2) quantitative measures on both the correctness and explainability of an interpretation. For a given black-box, our approach yields a set of Pareto-optimal interpretations with respect to the correctness and explainability measures. We show that the underlying multi-objective optimization problem can be solved via a reduction to quantitative constraint solving, such as weighted maximum satisfiability. To demonstrate the benefits of our approach, we have applied it to synthesize interpretations for black-box neural-network classifiers. Our experiments show that there often exists a rich and varied set of choices for interpretations that are missed by existing approaches.


title: Determinant of the finite volume Laplacian

abstract: The finite volume Laplacian can be defined in all dimensions and is a natural way to approximate the operator on a simplicial mesh. In the most general setting, its definition with orthogonal duals may require that not all volumes are positive; an example is the case corresponding to two-dimensional finite elements on a non-Delaunay triangulation. Nonetheless, in many cases two- and three-dimensional Laplacians can be shown to be negative semidefinite with a kernel consisting of constants. This work generalizes work in two dimensions that gives a geometric description of the Laplacian determinant; in particular, it relates the Laplacian determinant on a simplex in any dimension to certain volume quantities derived from the simplex geometry.


title: About the solution of the numerical instability for topological solitons with long range interaction

abstract: The computations of solutions of the field equations in the Model of Topological Particles, formulated with a scalar SU(2)-field, have shown instabilities leading to discrepancies between the numerical and analytical solutions. We identify the origin of these deviations in misalignments of the rotational axes corresponding to the SU(2) elements. The system of a single soliton we use as an example to show that a constraint suppressing the wave-like disturbances is able to lead to excellent agreement between the result of the numerical minimisation procedure and the analytical solution.


title: The physics potential of a reactor neutrino experiment with Skipper-CCDs: Searching for new physics with light mediators

abstract: We explore the sensitivity to new physics of the recently proposed vIOLETA experiment: a 10 kg Skipper Charged Coupled Device detector deployed 12 meters away from a commercial nuclear reactor core. We investigate two broad classes of models which benefit from the very low energy recoil threshold of these detectors, namely neutrino magnetic moments and light mediators coupled to neutrinos and quarks or electrons. We find that this experimental setup is very sensitive to light, weakly coupled new physics, and in particular that it could probe potential explanations of the event excess observed in XENON1T. We also provide a detailed study on the dependence of the sensitivity on the experimental setup assumptions and on the neutrino flux systematic uncertainties.


title: On the dimension of exceptional parameters for nonlinear projections, and the discretized Elekes-R\'onyai theorem

abstract: We consider four related problems. (1) Obtaining dimension estimates for the set of exceptional vantage points for the pinned Falconer distance problem. (2) Nonlinear projection theorems, in the spirit of Kaufman, Bourgain, and Shmerkin. (3) The parallelizability of planar $d$-webs. (4) The Elekes-R\'onyai theorem on expanding polynomials. Given a Borel set $A$ in the plane, we study the set of exceptional vantage points, for which the pinned distance $\Delta_p(A)$ has small dimension, that is, close to $(\dim A)/2$. We show that if this set has positive dimension, then it must have very special structure. This result follows from a more general single-scale nonlinear projection theorem, which says that if $\phi_1,\phi_2,\phi_3$ are three smooth functions whose associated 3-web has non-vanishing Blaschke curvature, and if $A$ is a $(\delta,\alpha)_2$-set in the sense of Katz and Tao, then at least one of the images $\phi_i(A)$ must have measure much larger than $|A|^{1/2}$, where $|A|$ stands for the measure of $A$. We prove analogous results for $d$ smooth functions $\phi_1,\ldots,\phi_d$, whose associated $d$-web is not parallelizable. We use similar tools to characterize when bivariate real analytic functions are "dimension expanding" when applied to a Cartesian product: if $P$ is a bivariate real analytic function, then $P$ is either locally of the form $h(a(x) + b(y))$, or $P(A,B)$ has dimension at least $\alpha+c$ whenever $A$ and $B$ are Borel sets with Hausdorff dimension $\alpha$. Again, this follows from a single-scale estimate, which is an analogue of the Elekes-R\'onyai theorem in the setting of the Katz-Tao discretized ring conjecture.


title: The X-ray binary population in the Galactic Center revealed through multi-decade observations

abstract: We present an investigation of the quiescent and transient X-ray binaries (XRBs) of the Galactic Center (GC). We extended our Chandra analysis of the non-thermal X-ray sources, located in the central parsec, from Hailey et al. (2018), using an additional 4.6 Msec of ACIS-S data obtained in 2012-2018. The individual Chandra spectra of the 12 sources fit to an absorbed power-law model with a mean photon index $\Gamma$~2 and show no Fe emission lines. Long-term variability was detected from nine of them, confirming that a majority are quiescent XRBs. Frequent X-ray monitoring of the GC revealed that the 12 non-thermal X-ray sources, as well as four X-ray transients have shown at most a single outburst over the last two decades. They are distinct from the six known neutron star LMXBs in the GC, which have all undergone multiple outbursts with <~ 5 year recurrence time on average. Based on the outburst history data of the broader population of X-ray transients, we conclude that the 16 sources represent a population of ~240-630 tightly-bound BH-LMXBs with ~4-12 hour orbital periods, consistent with the stellar/binary dynamics modelling in the vicinity of Sgr A*. The distribution of the 16 BH-LMXB candidates is disk-like (at 87% CL) and aligned with the nuclear star cluster. Our results have implications for XRB formation and the rate of gravitational wave events in other galactic nuclei.


title: Federated Asymptotics: a model to compare federated learning algorithms

abstract: We propose an asymptotic framework to analyze the performance of (personalized) federated learning algorithms. In this new framework, we formulate federated learning as a multi-criterion objective, where the goal is to minimize each client's loss using information from all of the clients. We analyze a linear regression model where, for a given client, we may theoretically compare the performance of various algorithms in the high-dimensional asymptotic limit. This asymptotic multi-criterion approach naturally models the high-dimensional, many-device nature of federated learning. These tools make fairly precise predictions about the benefits of personalization and information sharing in federated scenarios -- at least in our (stylized) model -- including that Federated Averaging with simple client fine-tuning achieves the same asymptotic risk as the more intricate meta-learning and proximal-regularized approaches and outperforming Federated Averaging without personalization. We evaluate these predictions on federated versions of the EMNIST, CIFAR-100, Shakespeare, and Stack Overflow datasets, where the experiments corroborate the theoretical predictions, suggesting such frameworks may provide a useful guide to practical algorithmic development.


title: A Matrix Model with BMS$_3$ Constraints

abstract: Our goal is to find a matrix model with $BMS_3$ constraints built in. These constraints are imposed through Loop equations. We solve them using a free field realisation of the algebra and write down the partition function in eigenvalue form. We comment on the nature of this partition function and it's relations with integrable systems.


title: Iterative learning control with discrete-time nonlinear nonminimum phase models via stable inversion

abstract: Output reference tracking can be improved by iteratively learning from past data to inform the design of feedforward control inputs for subsequent tracking attempts. This process is called iterative learning control (ILC). This article develops a method to apply ILC to systems with nonlinear discrete-time dynamical models with unstable inverses (i.e. discrete-time nonlinear non-minimum phase models). This class of systems includes piezoactuators, electric power converters, and manipulators with flexible links, which may be found in nanopositioning stages, rolling mills, and robotic arms, respectively. As these devices may be required to execute fine transient reference tracking tasks repetitively in contexts such as manufacturing, they may benefit from ILC. Specifically, this article facilitates ILC of such systems by presenting a new ILC synthesis framework that allows combination of the principles of Newton's root finding algorithm with stable inversion, a technique for generating stable trajectories from unstable models. The new framework, called Invert-Linearize ILC (ILILC), is validated in simulation on a cart-and-pendulum system with model error, process noise, and measurement noise. Where preexisting Newton-based ILC diverges, ILILC with stable inversion converges, and does so in less than one third the number of trials necessary for the convergence of a gradient-descent-based ILC technique used as a benchmark.


title: Heterotic String Model Building with Monad Bundles and Reinforcement Learning

abstract: We use reinforcement learning as a means of constructing string compactifications with prescribed properties. Specifically, we study heterotic SO(10) GUT models on Calabi-Yau three-folds with monad bundles, in search of phenomenologically promising examples. Due to the vast number of bundles and the sparseness of viable choices, methods based on systematic scanning are not suitable for this class of models. By focusing on two specific manifolds with Picard numbers two and three, we show that reinforcement learning can be used successfully to explore monad bundles. Training can be accomplished with minimal computing resources and leads to highly efficient policy networks. They produce phenomenologically promising states for nearly 100% of episodes and within a small number of steps. In this way, hundreds of new candidate standard models are found.


title: An Archival Search for Very-High-Energy Counterparts to Sub-Threshold Neutron-Star Merger Candidates

abstract: The recent discovery of electromagnetic signals in coincidence with gravitational waves from neutron-star mergers has solidified the importance of multimessenger campaigns for studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as the LIGO/Virgo gravitational wave detectors and the IceCube neutrino observatory, record many candidate signals that fall short of the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by these candidates may, when combined with time-coincident gamma-ray observations, lead to significant detections. In this contribution, I describe our use of sub-threshold binary neutron star merger candidates identified in Advanced LIGO's first observing run (O1) to search for transient events in very-high-energy gamma rays using archival observations from the VERITAS imaging atmospheric Cherenkov telescope array. I describe the promise of this technique for future joint sub-threshold searches.


title: Peak Sidelobe Level and Peak Crosscorrelation of Golay-Rudin-Shapiro Sequences

abstract: Sequences with low aperiodic autocorrelation and crosscorrelation are used in communications and remote sensing. Golay and Shapiro independently devised a recursive construction that produces families of complementary pairs of binary sequences. In the simplest case, the construction produces the Rudin-Shapiro sequences, and in general it produces what we call Golay-Rudin-Shapiro sequences. Calculations by Littlewood show that the Rudin-Shapiro sequences have low mean square autocorrelation. A sequence's peak sidelobe level is its largest magnitude of autocorrelation over all nonzero shifts. H{\o}holdt, Jensen, and Justesen showed that there is some undetermined positive constant $A$ such that the peak sidelobe level of a Rudin-Shapiro sequence of length $2^n$ is bounded above by $A(1.842626\ldots)^n$, where $1.842626\ldots$ is the positive real root of $X^4-3 X-6$. We show that the peak sidelobe level is bounded above by $5(1.658967\ldots)^{n-4}$, where $1.658967\ldots$ is the real root of $X^3+X^2-2 X-4$. Any exponential bound with lower base will fail to be true for almost all $n$, and any bound with the same base but a lower constant prefactor will fail to be true for at least one $n$. We provide a similar bound on the peak crosscorrelation (largest magnitude of crosscorrelation over all shifts) between the sequences in each Rudin-Shapiro pair. The methods that we use generalize to all families of complementary pairs produced by the Golay-Rudin-Shapiro recursion, for which we obtain bounds on the peak sidelobe level and peak crosscorrelation with the same exponential growth rate as we obtain for the original Rudin-Shapiro sequences.

categories: cs.IT cs.DM eess.SP math.CO math.IT math.NT


title: Useful relations for the analysis of stellar scintillation at the entrance pupil of a telescope

abstract: The development of new techniques for characterizing atmospheric optical turbulence (OT) has become an active topic of research again in recent years. In order to facilitate these studies, we reconsidered known theoretical results and obtained some new practically useful conclusions. We introduce a dimensionless Fresnel filter, which allows us to approximate a polychromatic weighting function (WF) by a monochromatic one with a typical precision of several percent. A so-called dimensionless WF can be easily scaled for a receiving aperture of any size. For the case of a circular aperture and monochromatic radiation, an analytical expression for the WF was found. The WFs for a square aperture and for a circular aperture match with relative difference less than 0.01 if the circular aperture diameter is 1.15 times larger than the square aperture side. A linear digital filter can be applied to the scintillation signal from an image detector. As an example of digital filtering, we considered the power law filter $\propto f^{5/3}$ with the WF being constant in a wide range of altitudes. We discuss the main limitations of this approach for measuring OT integral: finite pixel size, aliasing, and finite image detector size.


title: Sustaining Efficiency at Elevated Power Densities with InGaAs Air Bridge Cells

abstract: Here we investigate the use of single-junction InGaAs airbridge cells (ABCs) at elevated power densities. Such conditions are relevant to many thermophotovoltaic (TPV) applications, ranging from space to on-demand renewable electricity, and require effective management of heat and charge carriers. Experimental characterization of an InGaAs ABC with varying emitter and cell temperature is used to develop a predictive device model where carrier lifetimes and series resistances are the sole fitting parameters. The utility of this model is demonstrated through its use in identifying near-term opportunities for improving performance at elevated power densities, and for designing a thermal management strategy that maximizes overall power output. This model shows that an InGaAs ABC with material quality that leads to the longest reported carrier lifetimes can attain efficiencies exceeding 40% at 0.5 W/cm2, even when considering the power necessary to cool the cells.


title: The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites

abstract: A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $\delta \rho_{max}$=30\%(50\%) with respect to underdense regions, and $\delta \rho_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers.


title: QoT Assessment of the Optical Spectrum as a Service in Disaggregated Network Scenarios

abstract: The potential to operate third-party terminals over multi-domain transparent optical networks attracts operators and customers to implement Optical Spectrum as a Service (OSaaS). As infrastructure information cannot always be shared with OSaaS end customers, alternatives to off-line Quality of Transmission (QoT) estimation tools are required to assess the performance of the spectrum slot in order to estimate achievable throughput. In this paper, commercially available sliceable coherent transceivers are used to assess the Generalized Signal To Noise Ratio (GSNR) based QoT of the OSaaS in a live production network for both, narrow-band and wide-band OSaaS configurations. Extended channel probing based on symbol rate variability is combined with spectral sweeping and operation regime detection to characterize OSaaS implementations on 17 links with different underlying infrastructure configurations in order to maximize capacity and increase service margins in a low-margin operation regime. We achieve 0.05 dB estimation accuracy in GSNR for a wide-band spectrum services and 0.32 dB accuracy for narrow-band spectrum services. Based on the GSNR profile, spectral misalignment, spectral ripple, and operation regime are detected and service margin improvements are demonstrated. Finally, we discuss the network optimization perspective based on acquired data from channel probing and propose use-cases for continuous channel probing in transparent optical networks.


title: Clustering augmented Self-Supervised Learning: Anapplication to Land Cover Mapping

abstract: Collecting large annotated datasets in Remote Sensing is often expensive and thus can become a major obstacle for training advanced machine learning models. Common techniques of addressing this issue, based on the underlying idea of pre-training the Deep Neural Networks (DNN) on freely available large datasets, cannot be used for Remote Sensing due to the unavailability of such large-scale labeled datasets and the heterogeneity of data sources caused by the varying spatial and spectral resolution of different sensors. Self-supervised learning is an alternative approach that learns feature representation from unlabeled images without using any human annotations. In this paper, we introduce a new method for land cover mapping by using a clustering based pretext task for self-supervised learning. We demonstrate the effectiveness of the method on two societally relevant applications from the aspect of segmentation performance, discriminative feature representation learning and the underlying cluster structure. We also show the effectiveness of the active sampling using the clusters obtained from our method in improving the mapping accuracy given a limited budget of annotating.


title: First-Order Theory of Probabilistic Independence and Single-Letter Characterizations of Capacity Regions

abstract: We consider the first-order theory of random variables with the probabilistic independence relation, which concerns statements consisting of random variables, the probabilistic independence symbol, logical operators, and existential and universal quantifiers. Although probabilistic independence is the only non-logical relation included, this theory is surprisingly expressive, and is able to interpret the true first-order arithmetic over natural numbers (and hence is undecidable). We also construct a single-letter characterization of the capacity region for a general class of multiuser coding settings (including broadcast channel, interference channel and relay channel), using a first-order formula. We then introduce the linear entropy hierarchy to classify single-letter characterizations according to their complexity.


title: Role of New Kernel Function in Complexity Analysis of an Interior Point Algorithm for Semi definite Linear Complementarity Problem

abstract: In this paper, we introduce a new kernel function which differs from previous functions, and play an important role for generating a new design of primal-dual interior point algorithms for semidefinite linear complementarity problem. Its properties, allow us a great simplicity for the analysis of interior-point method, therefore the complexity of large-update primal-dual interior point is the best so far. Numerical tests have shown that the use of this function gave a big improvement in the results concerning the time and the number of iterations. so is well promising and perform well enough in practice in comparison with some other existing results in the literature.


title: Liquid Transfer for Viscoelastic Solutions

abstract: Viscoelastic liquid transfer from one surface to another is a process that finds applications in many technologies, primarily in printing. Here, cylindrical shaped capillary bridges pinned between two parallel disks are considered. Specifically, the effects of polymer mass fraction, solution viscosity, disk diameter, initial aspect ratio, final aspect ratio, stretching velocity and filling fraction (alike contact angle) are experimentally investigated in uniaxial extensional flow. Both Newtonian and viscoelastic polymer solutions are prepared using polyethylene glycol (PEG) and polyethylene oxide (PEO), with a wide variety of mass fractions. The results show that the increase in polymer mass fraction and solvent viscosity reduces the liquid transfer to the top surface. Moreover, the increase in the initial and final stretching height of the capillary bridge also decreases the liquid transfer, for both Newtonian and viscoelastic solutions. Finally, the shape of the capillary bridge is varied by changing the liquid volume. Now, Newtonian and viscoelastic solutions exhibit opposite behaviors for the liquid transfer. These findings are discussed in terms of interfacial shape instability and gravitational drainage.


title: Delocalized states in three-terminal superconductor-semiconductor nanowire devices

abstract: We fabricate three-terminal hybrid devices with a nanowire segment proximitized by a superconductor, and with two tunnel probe contacts on either side of that segment. We perform simultaneous tunneling measurements on both sides. We identify some states as delocalized above-gap states observed on both ends, and some states as localized near one of the tunnel barriers. Delocalized states can be traced from zero to finite magnetic fields beyond 0.5 T. In the parameter regime of delocalized states, we search for correlated subgap resonances required by the Majorana zero mode hypothesis. While both sides exhibit ubiquitous low-energy features at high fields, no correlation is inferred. Simulations using a one-dimensional effective model suggest that delocalized states may belong to lower one-dimensional subbands, while the localized states originate from higher subbands. To avoid localization in higher subbands, disorder may need to be further reduced to realize Majorana zero modes.


title: Tidally induced multipole moments of a nonrotating black hole vanish to all post-Newtonian orders

abstract: The tidal Love numbers of a black hole vanish, and this is often taken to imply that the hole's tidally induced multipole moments vanish also. An obstacle to establishing a link between these statements is that the multipole moments of individual bodies are not defined in general relativity, when the bodies are subjected to a mutual gravitational interaction. In a previous publication [Phys. Rev. D 103, 064023 (2021)] I promoted the view that individual multipole moments can be defined when the mutual interaction is sufficiently weak to be described by a post-Newtonian expansion. In this view, a compact body is perceived far away as a skeletonized post-Newtonian object with a multipole structure, and the multipole moments can then be related to the body's Love numbers. I expand on this view, and demonstrate that all static, tidally induced, mass multipole moments of a nonrotating black hole vanish to all post-Newtonian orders. The proof rests on a perturbative solution to the Einstein-Maxwell equations that describes an electrically charged particle placed in the presence of a charged black hole. The gravitational attraction between particle and black hole is balanced by electrostatic repulsion, and the system is in an equilibrium state. The particle provides a tidal environment to the black hole, and the multipole moments vanish for this environment. I argue that the vanishing is robust, and applies to all slowly-varying tidal environments. The black hole's charge can be as small as desired (though not identically zero); by continuity, the multipole moments of an electrically neutral black hole will continue to vanish.


title: TUM-VIE: The TUM Stereo Visual-Inertial Event Dataset

abstract: Event cameras are bio-inspired vision sensors which measure per pixel brightness changes. They offer numerous benefits over traditional, frame-based cameras, including low latency, high dynamic range, high temporal resolution and low power consumption. Thus, these sensors are suited for robotics and virtual reality applications. To foster the development of 3D perception and navigation algorithms with event cameras, we present the TUM-VIE dataset. It consists of a large variety of handheld and head-mounted sequences in indoor and outdoor environments, including rapid motion during sports and high dynamic range scenarios. The dataset contains stereo event data, stereo grayscale frames at 20Hz as well as IMU data at 200Hz. Timestamps between all sensors are synchronized in hardware. The event cameras contain a large sensor of 1280x720 pixels, which is significantly larger than the sensors used in existing stereo event datasets (at least by a factor of ten). We provide ground truth poses from a motion capture system at 120Hz during the beginning and end of each sequence, which can be used for trajectory evaluation. TUM-VIE includes challenging sequences where state-of-the art visual SLAM algorithms either fail or result in large drift. Hence, our dataset can help to push the boundary of future research on event-based visual-inertial perception algorithms.


title: Weakly Supervised Classification Using Group-Level Labels

abstract: In many applications, finding adequate labeled data to train predictive models is a major challenge. In this work, we propose methods to use group-level binary labels as weak supervision to train instance-level binary classification models. Aggregate labels are common in several domains where annotating on a group-level might be cheaper or might be the only way to provide annotated data without infringing on privacy. We model group-level labels as Class Conditional Noisy (CCN) labels for individual instances and use the noisy labels to regularize predictions of the model trained on the strongly-labeled instances. Our experiments on real-world application of land cover mapping shows the utility of the proposed method in leveraging group-level labels, both in the presence and absence of class imbalance.


title: Volatility of Sodium in Carbonaceous Chondrites at Temperatures Consistent with Low-Perihelia Asteroids

abstract: Solar system bodies with surface and sub-surface volatiles will show observational evidence of activity when they reach a temperature where those volatiles change from solid to gas and are released. This is most frequently seen in comets, where activity is driven by the sublimation of water, carbon dioxide, or carbon monoxide ices. However, some bodies (notably the asteroid (3200) Phaethon) show initiation of activity at very small heliocentric distances, long after they have reached the sublimation temperatures of these ices. We investigate whether the sodium present in the mineral matrix could act as the volatile element responsible for this activity. We conduct theoretical modeling which indicates that sodium has the potential to sublimate in the conditions that Phaethon experiences, depending on the mineral phase it is held in. To test this, we then exposed samples of the carbonaceous chondrite Allende to varying heating events similar to what would be experienced by low perihelion asteroids. We measured the change in sodium present in each sample, and find that the highest temperature samples show a significant loss of sodium from specific mineral phases over a single heating event, comparable to a day on the surface of Phaethon. Under specific thermal histories possible for Phaethon, this outgassing could be sufficient to explain this object's observed activity. This effect would also be expected to be observed for other low-perihelia asteroids as well, and may act as a critical step in the process of disrupting small low-albedo asteroids.


title: Andreev spectroscopy of nonhelical spin textures in topological insulators

abstract: We study how nonhelical spin textures affect the proximity-induced superconductivity of topological insulator (TI)-superconductor (SC) interface states. In particular we calculate the conductance of lateral heterojunctions which comprise a TI surface that is only partially covered by a superconducting material. Interface potentials at the TI-SC interface may lead to a Fermi velocity and spin texture mismatch between the two regions of the lateral heterojunction. By enforcing the hermiticity of the total Hamiltonian, we derive the boundary conditions and calculate the conductance of the structure in both the normal and superconducting state. The total Andreev conductance is calculated for both $s$-wave and spin-triplet parent SCs, and for several examples of nonhelical spin textures which lead to different Fermi surface mismatches between the two planar regions of the heterojunction. We find that for spin-triplet SCs, nonzero conductance signatures only appear for certain combinations of nonhelical spin textures and parent superconducting material.


title: Localized Phonon Densities of States at Grain Boundaries in Silicon

abstract: Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope it is of interest to explore whether defects such as dislocations or grain boundaries will result in measurable changes of the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. Since these boundaries are modeled by supercells with up to 160 atoms, the density of states was calculated by taking the Fourier transform of the velocity-velocity autocorrelation function from molecular dynamics simulations based on new supercells doubled in each direction. In select cases the results were checked on the original supercells with fewer atoms by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary the height of the optic phonon peak in the density of states at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths. It also means that changes in the density of states intrinsic to grain boundaries are unlikely to affect thermal conductivity at ambient temperatures, which are most likely dominated by the scattering of acoustic phonons.


title: A note on inverse results of random walks in Abelian groups

abstract: In this short note we give various near optimal characterizations of random walks over finite Abelian groups with large maximum discrepancy from the uniform measure. We also provide several interesting connections to existing results in the literature.


title: Augmenting control arms with Real-World Data for cancer trials: Hybrid control arm methods and considerations

abstract: Randomized controlled trials (RCTs) are the gold standard for assessing drug safety and efficacy. However, RCTs have some drawbacks which have led to the use of single-arm studies to make certain internal drug development and regulatory decisions, particularly in oncology. Hybrid controlled trials with real-world data (RWD), in which the control arm is composed of both trial and real-world patients, have the potential to help address some of the shortcomings of both RCTs and single-arm studies in particular situations, such as when a disease has low prevalence or when the standard of care to be used in the control arm is ineffective or highly toxic and an experimental therapy shows early promise. This paper discusses why it may be beneficial to consider hybrid controlled trials with RWD, what such a design entails, when it may be appropriate, and how to conduct the analyses. We propose a novel two-step borrowing method for the construction of hybrid control arms. We use simulations to demonstrate the operating characteristics of dynamic and static borrowing methods, and highlight the trade-offs and analytic decisions that study teams will need to address when designing a hybrid study.


title: Fermions on a torus knot

abstract: In this work, we investigate the effects of a nontrivial topology (and geometry) of a system considering \textit{interacting} and \textit{noninteracting} particle modes, which are restricted to follow a closed path over the torus surface. In order to present a prominent thermodynamical investigation of this system configuration, we carry out a detailed analysis using statistical mechanics within the grand canonical ensemble approach to deal with \textit{noninteracting} fermions. In an analytical manner, we study the following thermodynamic functions in such context: the Helmholtz free energy, the mean energy, the magnetization and the susceptibility. Further, we take into account the behavior of Fermi energy of the thermodynamic system. Finally, we briefly outline how to proceed in case of \textit{interacting} fermions.


title: Generative Relation Linking for Question Answering over Knowledge Bases

abstract: Relation linking is essential to enable question answering over knowledge bases. Although there are various efforts to improve relation linking performance, the current state-of-the-art methods do not achieve optimal results, therefore, negatively impacting the overall end-to-end question answering performance. In this work, we propose a novel approach for relation linking framing it as a generative problem facilitating the use of pre-trained sequence-to-sequence models. We extend such sequence-to-sequence models with the idea of infusing structured data from the target knowledge base, primarily to enable these models to handle the nuances of the knowledge base. Moreover, we train the model with the aim to generate a structured output consisting of a list of argument-relation pairs, enabling a knowledge validation step. We compared our method against the existing relation linking systems on four different datasets derived from DBpedia and Wikidata. Our method reports large improvements over the state-of-the-art while using a much simpler model that can be easily adapted to different knowledge bases.


title: Sterile neutrinos in $\Lambda_b^0\to (\Lambda_c^+,p^+)\ell^-_1\ell^-_2\ell^+_3\nu$ decays

abstract: We study lepton number violating and lepton number conserving semi-leptonic decays of heavy baryons $\Lambda_b^0$ to three charged leptons and a neutrino. The decays occur through two intermediate quasi-degenerate GeV-scale sterile neutrinos of either Majorana or Dirac type that can be on-shell. Interference between the intermediate heavy neutrinos leads to CP violation in the final states. Effect of neutrino oscillations between the heavy states are also considered in observables of interests, \emph{i.e.,} branching ratio, and the CP-asymmetry. Given the present constraints on the heavy-to-light mixing elements $|V_{eN}|$ and $|V_{\mu N}|$, CP-averaged branching ratio of $\Lambda_b^0\to (\Lambda_c,p)\mu\mu e\nu$ with intermediate Majorana neutrinos is almost two orders of magnitude larger than for the same with Dirac neutrinos, whereas, CP-averaged branching ratio of $\Lambda_b^0\to (\Lambda_c,p)ee\mu\nu$ is of the same order of magnitude for both Majorana and Dirac neutrino induced decays. CP-violation is found to be appreciable when the neutrino mass difference is comparable with the average decay widths.


title: Classification of Common Waveforms Including a Watchdog for Unknown Signals

abstract: In this paper, we examine the use of a deep multi-layer perceptron model architecture to classify received signal samples as coming from one of four common waveforms, Single Carrier (SC), Single-Carrier Frequency Division Multiple Access (SC-FDMA), Orthogonal Frequency Division Multiplexing (OFDM), and Linear Frequency Modulation (LFM), used in communication and radar networks. Synchronization of the signals is not needed as we assume there is an unknown and uncompensated time and frequency offset. An autoencoder with a deep CNN architecture is also examined to create a new fifth classification category of an unknown waveform type. This is accomplished by calculating a minimum and maximum threshold values from the root mean square error (RMSE) of the radar and communication waveforms. The classifier and autoencoder work together to monitor a spectrum area to identify the common waveforms inside the area of operation along with detecting unknown waveforms. Results from testing showed the classifier had 100\% classification rate above 0 dB with accuracy of 83.2\% and 94.7\% at -10 dB and -5 dB, respectively, with signal impairments present. Results for the anomaly detector showed 85.3\% accuracy at 0 dB with 100\% at SNR greater than 0 dB with signal impairments present when using a high-value Fast Fourier Transform (FFT) size. Accurate detection rates decline as additional noise is introduced to the signals, with 78.1\% at -5 dB and 56.5\% at -10 dB. However, these low rates seen can be potentially mitigated by using even higher FFT sizes also shown in our results.


title: Detecting changes in covariance via random matrix theory

abstract: A novel method is proposed for detecting changes in the covariance structure of moderate dimensional time series. This non-linear test statistic has a number of useful properties. Most importantly, it is independent of the underlying structure of the covariance matrix. We discuss how results from Random Matrix Theory, can be used to study the behaviour of our test statistic in a moderate dimensional setting (i.e. the number of variables is comparable to the length of the data). In particular, we demonstrate that the test statistic converges point wise to a normal distribution under the null hypothesis. We evaluate the performance of the proposed approach on a range of simulated datasets and find that it outperforms a range of alternative recently proposed methods. Finally, we use our approach to study changes in the amount of water on the surface of a plot of soil which feeds into model development for degradation of surface piping.


title: Astrometric Constraints on a Massive Black Hole Binary in NGC\,4472

abstract: In an EHT study of a Jy-level target, Safarzadeh et al. (2019) show how astrometric monitoring could constrain massive black hole binaries with the wide separations that make them long-lived against gravitational wave losses, and with the small mass ratios expected from merged satellite galaxies. With this ngVLA study, we show how such frontier topics could be explored for the more numerous mJy-level targets, such as NGC\,4472. We also discuss how ngVLA astrometric monitoring could test the upper limits from pulsar timing arrays on gravitational waves from NGC\,4472.


title: Density control of interacting agent systems

abstract: We consider the problem of controlling the group behavior of a large number of dynamic systems that are constantly interacting with each other. These systems are assumed to have identical dynamics (e.g., birds flock, robot swarm) and their group behavior can be modeled by a distribution. Thus, this problem can be viewed as an optimal control problem over the space of distributions. We propose a novel algorithm to compute a feedback control strategy so that, when adopted by the agents, the distribution of them would be transformed from an initial one to a target one over a finite time window. Our method is built on optimal transport theory but differs significantly from existing work in this area in that our method models the interactions among agents explicitly. From an algorithmic point of view, our algorithm is based on a generalized version of the proximal gradient descent algorithm and has a convergence guarantee with a sublinear rate. We further extend our framework to account for the scenarios where the agents are from multiple species. In the linear quadratic setting, the solution is characterized by coupled Riccati equations which can be solved in closed-form. Finally, several numerical examples are presented to illustrate our framework.


title: Mass-Loading and Non-Spherical Divergence in Hot Galactic Winds: Implications for X-ray Observations

abstract: Cool clouds are expected to be destroyed and incorporated into hot supernova-driven galactic winds. The mass-loading of a wind by the cool medium modifies the bulk velocity, temperature, density, entropy, and abundance profiles of the hot phase relative to an un-mass-loaded outflow. We provide general equations and limits for this physics that can be used to infer the rate of cool gas entrainment from X-ray observations, accounting for non-spherical expansion. In general, mass-loading flattens the density and temperature profiles, decreases the velocity and increases the entropy if the Mach number is above a critical value. We first apply this model to a recent high-resolution galactic outflow simulation where the mass-loading can be directly inferred. We show that the temperature, entropy, and composition profiles are well-matched, providing evidence that this physics sets the bulk hot gas profiles. We then model the diffuse X-ray emission from the local starburst M82. The non-spherical (more cylindrical) outflow geometry is directly taken from the observed X-ray surface brightness profile. These models imply a total mass-loading rate that is about equal to that injected in the starburst, $\simeq 10$ M$_\odot$ yr$^{-1}$, and they predict an asymptotic hot wind velocity of $\sim 1000\,{\rm km \ s^{-1}}$ that is $\sim1.5-2$ times smaller than previous predictions. We also show how the observed entropy profile can be used to constrain the outflow velocity, making predictions for future missions like XRISM. We argue that the observed X-ray limb-brightening may be explained by mass-loading at the outflow's edges.


title: IsoScore: Measuring the Uniformity of Embedding Space Utilization

abstract: The recent success of distributed word representations has led to an increased interest in analyzing the properties of their spatial distribution. Several studies have suggested that contextualized word embedding models do not isotropically project tokens into vector space. However, current methods designed to measure isotropy, such as average random cosine similarity and the partition score, have not been thoroughly analyzed and are not appropriate for measuring isotropy. We propose IsoScore: a novel tool that quantifies the degree to which a point cloud uniformly utilizes the ambient vector space. Using rigorously designed tests, we demonstrate that IsoScore is the only tool available in the literature that accurately measures how uniformly distributed variance is across dimensions in vector space. Additionally, we use IsoScore to challenge a number of recent conclusions in the NLP literature that have been derived using brittle metrics of isotropy. We caution future studies from using existing tools to measure isotropy in contextualized embedding space as resulting conclusions will be misleading or altogether inaccurate.


title: An NLP approach to quantify dynamic salience of predefined topics in a text corpus

abstract: The proliferation of news media available online simultaneously presents a valuable resource and significant challenge to analysts aiming to profile and understand social and cultural trends in a geographic location of interest. While an abundance of news reports documenting significant events, trends, and responses provides a more democratized picture of the social characteristics of a location, making sense of an entire corpus to extract significant trends is a steep challenge for any one analyst or team. Here, we present an approach using natural language processing techniques that seeks to quantify how a set of pre-defined topics of interest change over time across a large corpus of text. We found that, given a predefined topic, we can identify and rank sets of terms, or n-grams, that map to those topics and have usage patterns that deviate from a normal baseline. Emergence, disappearance, or significant variations in n-gram usage present a ground-up picture of a topic's dynamic salience within a corpus of interest.


title: An Efficient Parallel Algorithm for finding Bridges in a Dense Graph

abstract: This paper presents a simple and efficient approach for finding the bridges and failure points in a densely connected network mapped as a graph. The algorithm presented here is a parallel algorithm which works in a distributed environment. The main idea of our algorithm is to generate a sparse certificate for a graph and finds bridges using a simple DFS (Depth First Search). We first decompose the graph into independent and minimal subgraphs using a minimum spanning forest algorithm. To identify the bridges in the graph network, we convert these subgraphs into a single compressed graph and use a DFS approach to find bridges. The approach presented here is optimized for the use cases of dense graphs and gives the time complexity of O(E/M + Vlog(M)), for a given graph G(V,E) running on M machines.


title: Issues with Positivity-Preserving Patankar-type Schemes

abstract: Patankar-type schemes are linearly implicit time integration methods designed to be unconditionally positivity-preserving. However, there are only little results on their stability or robustness. We suggest two approaches to analyze the performance and robustness of these methods. In particular, we demonstrate problematic behaviors of these methods that, even on very simple linear problems, can lead to undesired oscillations and order reduction for vanishing initial condition. Finally, we demonstrate in numerical simulations that our theoretical results for linear problems apply analogously to nonlinear stiff problems.


title: Elephants or Goldfish? An Empirical Analysis of Carrier Reciprocity in Dynamic Freight Markets

abstract: Dynamic macroeconomic conditions and non-binding truckload freight contracts enable both shippers and carriers to behave opportunistically. We present an empirical analysis of carrier reciprocity in the US truckload transportation sector to demonstrate whether consistent performance and fair pricing by shippers when markets are in their favor result in maintained primary carrier tender acceptance when markets turn. The results suggest carriers have short memories; they do not remember shippers' previous period pricing, tendering behavior, or performance when making freight acceptance decisions. However, carriers appear to be myopic and respond to shippers' current market period behaviors, ostensibly without regard to shippers' previous behaviors.


title: Lights Out On A Random Graph

abstract: We consider the generalized game Lights Out played on a graph and investigate the following question: for a given positive integer $n$, what is the probability that a graph chosen uniformly at random from the set of graphs with $n$ vertices yields a universally solvable game of Lights Out? When $n \leq 11$, we compute this probability exactly by determining if the game is universally solvable for each graph with $n$ vertices. We approximate this probability for each positive integer $n$ with $n \leq 100$ by applying a Monte Carlo method using 1,000,000 trials. We also perform the analogous computations for connected graphs.


title: Energy Balance at Interplanetary Shocks: In-situ Measurement of the Fraction in Supra-thermal and Energetic Particles with ACE and Wind

abstract: The acceleration of charged particles by interplanetary shocks can drain a non-negligible fraction of the upstream ram pressure. For a sample of shocks observed in-situ at 1 AU by the ACE and Wind spacecraft, time-series of the non-Maxwellian components (supra-thermal and higher-energy) of the ion and electron energy spectra were acquired for each event. These were averaged for one hour before and after the time of the shock passage to determine their partial pressure. Using the MHD Rankine-Hugoniot jump conditions, we find that the fraction of the total upstream energy flux transferred to non-Maxwellian downstream particles is typically about 2-16%. Notably, our sample shows that neither the fast magnetosonic Mach number nor the angle between the shock normal and average upstream magnetic field are correlated with non-Maxwellian particle pressure.


title: Uncover band topology via quantized drift in two-dimensional Bloch oscillations

abstract: We propose to measure band topology via quantized drift of Bloch oscillations in a two-dimensional Harper-Hofstadter lattice subjected to tilted fields in both directions. When the difference between the two tilted fields is large, Bloch oscillations uniformly sample all momenta, and hence the displacement in each direction tends to be quantized at multiples of the overall period, regardless of any momentum of initial state. The quantized displacement is related to a reduced Chern number defined as a line integral of Berry curvature in each direction, providing an almost perfect measurement of Chern number. Our scheme can apply to detect Chern number and topological phase transitions not only for the energy-separable band, but also for energy-inseparable bands which cannot be achieved by conventional Thouless pumping or integer quantum Hall effect.


title: Principal bundle groupoids, their gauge group and their nerve

abstract: We consider groupoids in the category of principal bundles, which we call principal bundles (PB) groupoids. Inspired by work by Th. Nikolaus and K. Waldorf, we generalise bundle gerbes over manifolds to bundle gerbes over groupoids and discuss a functorial correspondence between PB groupoids and bundle gerbes over groupoids. From a PB groupoid over a fibre product groupoid, we build a bundle gerbe over another fibre product groupoid. Conversely, from a bundle gerbe over a Lie groupoid, we build a PB groupoid. It has a trivial base and from any PB groupoid with trivial base, we build a bundle gerbe over a Lie groupoid. In that case, the resulting bundle gerbe is isomorphic as a groupoid to a partial quotient of the PB groupoid. We describe the nerves of PB groupoids and their partial quotients, which are simplicial objects in the category of principal bundles. Applying this construction enables us to define the inner transformation group of the nerve of a partial quotient groupoid and to describe the transformations of the corresponding bundle gerbe.


title: Scene Designer: a Unified Model for Scene Search and Synthesis from Sketch

abstract: Scene Designer is a novel method for searching and generating images using free-hand sketches of scene compositions; i.e. drawings that describe both the appearance and relative positions of objects. Our core contribution is a single unified model to learn both a cross-modal search embedding for matching sketched compositions to images, and an object embedding for layout synthesis. We show that a graph neural network (GNN) followed by Transformer under our novel contrastive learning setting is required to allow learning correlations between object type, appearance and arrangement, driving a mask generation module that synthesises coherent scene layouts, whilst also delivering state of the art sketch based visual search of scenes.


title: Private Delivery Networks -- Extended Abstract

abstract: The past decade has seen tremendous shifts in how people live, work, and buy goods, with an increased reliance on e-commerce and deliveries. Purchase histories generated through e-commerce can be highly personal, revealing identifying information about individuals and households. Constructing profiles from these data allows for the targeting of individuals and communities through practices such as targeted marketing and information campaigns. Furthermore, when purchase profiles are connected with delivery addresses, these data can measure the demographics of a local community and allow for individualized targeting to reach beyond the digital realm to the physical one. Events that accelerated shifts towards e-commerce, such as an infectious disease epidemic, have also widened equity gaps. This work is about alternative e-commerce delivery network models that address both rising privacy and wealth inequality concerns. This includes strategies that mask and add noise to purchase histories, and allow people to "buy privacy" through charitable contributions.


title: Card guessing and the birthday problem for sampling without replacement

abstract: Consider a uniformly random deck consisting of cards labelled by numbers from $1$ through $n$, possibly with repeats. A guesser guesses the top card, after which it is revealed and removed and the game continues. What is the expected number of correct guesses under the best and worst strategies? We establish sharp asymptotics for both strategies. For the worst case, this answers a recent question of Diaconis, Graham, He and Spiro, who found the correct order. As part of the proof, we study the birthday problem for sampling without replacement using Stein's method.


title: Stochastic Optimization under Distributional Drift

abstract: We consider the problem of minimizing a convex function that is evolving according to unknown and possibly stochastic dynamics, which may depend jointly on time and on the decision variable itself. Such problems abound in the machine learning and signal processing literature, under the names of concept drift, stochastic tracking, and performative prediction. We provide novel non-asymptotic convergence guarantees for stochastic algorithms with iterate averaging, focusing on bounds valid both in expectation and with high probability. The efficiency estimates we obtain clearly decouple the contributions of optimization error, gradient noise, and time drift. Notably, we identify a low drift-to-noise regime in which the tracking efficiency of the proximal stochastic gradient method benefits significantly from a step decay schedule. Numerical experiments illustrate our results.


title: Task-Oriented Multi-User Semantic Communications for VQA Task

abstract: Semantic communications focus on the transmission of semantic features. In this letter, we consider a task-oriented multi-user semantic communication system for multimodal data transmission. Particularly, partial users transmit images while the others transmit texts to inquiry the information about the images. To exploit the correlation among the multimodal data from multiple users, we propose a deep neural network enabled semantic communication system, named MU-DeepSC, to execute the visual question answering (VQA) task as an example. Specifically, the transceiver for MU-DeepSC is designed and optimized jointly to capture the features from the correlated multimodal data for task-oriented transmission. Simulation results demonstrate that the proposed MU-DeepSC is more robust to channel variations than the traditional communication systems, especially in the low signal-to-noise (SNR) regime.


title: Degeneracies Between Self-interacting Dark Matter and Supernova Feedback as cusp-core transformation mechanisms

abstract: We present a suite of 16 high-resolution hydrodynamic simulations of an isolated dwarf galaxy (gaseous and stellar disk plus a stellar bulge) within an initially cuspy dark matter (DM) halo, including self-interactions between the DM particles (SIDM); as well as stochastic star formation and subsequent supernova feedback (SNF), implemented using the stellar feedback model SMUGGLE. The simulations start from identical initial conditions and we regulate the strength of SIDM and SNF by systematically varying the SIDM momentum transfer cross section and the gas density threshold for star formation. The DM halo forms a constant density core of similar size and shape for several combinations of those two parameters. Haloes with cores that are formed due to SIDM (adiabatic cusp-core transformation) have velocity dispersion profiles which are closer to isothermal than those of haloes with cores that are formed due to SNF in simulations with bursty star formation (impulsive cusp-core transformation). Impulsive SNF can generate positive stellar age gradients and increase random motion in the gas at the centre of the galaxy. Simulated galaxies in haloes with cores that were formed adiabatically are spatially more extended, with stellar metallicity gradients that are shallower (at late times) than those of galaxies in other simulations. Such observable properties of the gas and the stars, which indicate either an adiabatic or an impulsive evolution of the gravitational potential, may be used to determine whether observed cores in DM haloes are formed through self-interactions between the DM particles or in response to impulsive SNF.


title: Approximating the Permanent with Deep Rejection Sampling

abstract: We present a randomized approximation scheme for the permanent of a matrix with nonnegative entries. Our scheme extends a recursive rejection sampling method of Huber and Law (SODA 2008) by replacing the upper bound for the permanent with a linear combination of the subproblem bounds at a moderately large depth of the recursion tree. This method, we call deep rejection sampling, is empirically shown to outperform the basic, depth-zero variant, as well as a related method by Kuck et al. (NeurIPS 2019). We analyze the expected running time of the scheme on random $(0, 1)$-matrices where each entry is independently $1$ with probability $p$. Our bound is superior to a previous one for $p$ less than $1/5$, matching another bound that was known to hold when every row and column has density exactly $p$.


title: Structure of $^{43}$P and $^{42}$Si in a two-level shape-coexistence model

abstract: Exclusive cross sections for the $^{43}$P$(-1p)^{42}$Si reaction to the lowest $0^+$ and $2^+$ states, measured at NSCL with GRETINA and the S800, are interpreted in terms of a two-level mixing (collective) model of oblate and prolate co-existing shapes. Using the formalism developed for deformed nuclei we calculate the spectroscopic amplitudes and exclusive cross-sections in the strong coupling limit, where for $^{43}$P the schematic wavefunction includes the coupling of the Nilsson [211]$\frac{1}{2}$ proton orbit. Good agreement with the experimental data is obtained when the amplitude of the oblate configuration is $\gtrsim$ 80\%, suggesting that both nuclei are predominantly oblate, in line with theoretical expectations.


title: Dynamical stability in the vicinity of Saturnian small moons. The cases of Aegaeon, Methone, Anthe and Pallene

abstract: In this work we analyze the orbital evolution and the dynamical stability in the vicinity of the small Saturnian moons Aegaeon, Methone, Anthe and Pallene. We numerically resolve the exact equations of motions to investigate the orbital motion of thousands of test particles within and near to the domain of the 7/6, 14/15, 10/11 mean motion resonances of Aegaeon, Methone and Anthe with Mimas, respectively. We show that, for massless small moons, the orbits of particles initially restricted to the resonant domains remain stable for at least $10^4$ yr. We also conduct numerical simulations considering Aegaeon, Methone, Anthe and Pallene as massive bodies. The results show that most particles undergo significant perturbations in their orbital motions, ultimately destabilizing in timescales of a few hundreds of years or even less through collisions with the four small moons. In addition, we also simulate the orbital evolution of test particles initially distributed in form of arc around Aegaeon, Methone and Anthe. We show that the initial arcs are dynamically eroded on timescales of hundreds of years, allowing us to constraint the timescales for which gravitational forces operate to remove particles from the observed arcs.


title: A Game-Theoretic Approach to Self-Stabilization with Selfish Agents

abstract: Self-stabilization is an excellent approach for adding fault tolerance to a distributed multi-agent system. However, two properties of self-stabilization theory, convergence and closure, may not be satisfied if agents are selfish. To guarantee convergence, we formulate the problem as a stochastic Bayesian game and introduce probabilistic self-stabilization to adjust the probabilities of rules with behavior strategies. This satisfies agents' self-interests such that no agent deviates the rules. To guarantee closure in the presence of selfish agents, we propose fault-containment as a method to constrain legitimate configurations of the self-stabilizing system to be Nash equilibria. We also assume selfish agents as capable of performing unauthorized actions at any time, which threatens both properties, and present a stepwise solution to handle it. As a case study, we consider the problem of distributed clustering and propose five self-stabilizing algorithms for forming clusters. Simulation results show that our algorithms react correctly to rule deviations and outperform comparable schemes in terms of fairness and stabilization time.


title: Neutron Scattering Signature of Phonon Renormalization in Nickel (II) Oxide

abstract: The physics of mutual interaction of phonon quasiparticles with electronic spin degrees of freedom, leading to unusual transport phenomena of spin and heat, has been a subject of continuing interests for decades. Despite its pivotal role in transport processes, the effect of spin-phonon coupling on the phonon system, especially acoustic phonon properties, has so far been elusive. By means of inelastic neutron scattering and first-principles calculations, anomalous scattering spectral intensity from acoustic phonons was identified in the exemplary collinear antiferromagnetic nickel (II) oxide, unveiling strong spin-lattice correlations that renormalize the polarization of acoustic phonon. In particular, a clear magnetic scattering signature of the measured neutron scattering intensity from acoustic phonons is demonstrated by its momentum transfer and temperature dependences. The anomalous scattering intensity is successfully modeled with a modified magneto-vibrational scattering cross section, suggesting the presence of spin precession driven by phonon. The renormalization of phonon eigenvector is indicated by the observed "geometry-forbidden" neutron scattering intensity from transverse acoustic phonon. Importantly, the eigenvector renormalization cannot be explained by magnetostriction but instead, it could result from the coupling between phonon and local magnetization of ions.


title: Fano resonance and incoherent interlayer excitons in molecular van der Waals heterostructures

abstract: Complex van der Waals heterostructures from layered molecular stacks are promising optoelectronic materials offering means to efficient, modular charge separation and collection layers. The effect of stacking in the electrodynamics of such hybrid organic-inorganic two-dimensional materials remains largely unexplored, whereby molecular scale engineering could lead to advanced optical phenomena. For instance, tunable Fano engineering could make possible on-demand transparent conducting layers or photoactive elements, and passive cooling. We employ an adapted Gersten-Nitzan model and real time time-dependent density functional tight-binding to study the optoelectronics of self-assembled monolayers on graphene nanoribbons. We find Fano resonances that cause electromagnetic induced opacity and transparency, and reveal an additional incoherent process leading to interlayer exciton formation with a characteristic charge transfer rate. These results showcase hybrid van der Waals heterostructures as paradigmatic 2D optoelectronic stacks, featuring tunable Fano optics and unconventional charge transfer channels.


title: Bell-state generation for spin qubits via dissipative coupling

abstract: We theoretically investigate the dynamics of two spin qubits interacting with a magnetic medium. A systematic formal framework for this qubit-magnet hybrid system is developed in terms of the steady-state properties of the magnetic medium. Focusing on the induced dissipative coupling between the spin qubits, we show how a sizable long-lived entanglement can be established via the magnetic environment, in the absence of any coherent coupling. Moreover, we demonstrate that maximally-entangled two-qubit states (Bell states) can be achieved in this scheme when complemented by proper postselection. In this situation, the time evolution of the entanglement is governed by a non-Hermitian Hamiltonian, where dynamical phases are separated by an exceptional point. The resultant Bell state is robust against weak random perturbations and does not require the preparation of a particular initial state. Our study may find applications in quantum information science, quantum spintronics, and for sensing of nonlocal quantum correlations.


title: The Visibility Center of a Simple Polygon

abstract: We introduce the \emph{visibility center} of a set of points inside a polygon -- a point $c_V$ such that the maximum geodesic distance from $c_V$ to see any point in the set is minimized. For a simple polygon of $n$ vertices and a set of $m$ points inside it, we give an $O((n+m) \log {(n+m)})$ time algorithm to find the visibility center. We find the visibility center of \emph{all} points in a simple polygon in $O(n \log n)$ time. Our algorithm reduces the visibility center problem to the problem of finding the geodesic center of a set of half-polygons inside a polygon, which is of independent interest. We give an $O((n+k) \log (n+k))$ time algorithm for this problem, where $k$ is the number of half-polygons.


title: The strong CP problem and higher dimensional gauge theories

abstract: We discuss a natural scenario to solve the strong CP problem in the framework of the higher dimensional gauge theory. An axion-like field $A_y$ has been built-in as the extra-space component of the higher dimensional gauge field. The coupling of $A_y$ with gluons is attributed to the radiatively induced "Chern-Simons" (CS) term. We adopt a toy model with some unknown gauge symmetry U(1)$_X$. The CS term is obtained in two ways: first by a concrete 1-loop calculation and next by use of the Fujikawa's method to deal with the chiral anomaly in 4D space-time. The obtained results are identical, which implies that the radiative correction to the CS term is "1-loop exact" and is also free from UV-divergence even though the theory itself is non-renormalizable. As a novel feature of this scenario, such obtained CS term is no longer linear in the field $A_y$ as in the usually discussed CS term in 5D space-time but a periodic function of $A_y$, since $A_y$ has a physical meaning as the Wilson-loop phase. We argue how such novel feature of this scenario causes the modification of the ordinary solutions of the strong CP problem based on the axion fields.


title: CaraNet: Context Axial Reverse Attention Network for Segmentation of Small Medical Objects

abstract: Segmenting medical images accurately and reliably is important for disease diagnosis and treatment. It is a challenging task because of the wide variety of objects' sizes, shapes, and scanning modalities. Recently, many convolutional neural networks (CNN) have been designed for segmentation tasks and achieved great success. Few studies, however, have fully considered the sizes of objects, and thus most demonstrate poor performance for small objects segmentation. This can have a significant impact on the early detection of diseases. This paper proposes a Context Axial Reserve Attention Network (CaraNet) to improve the segmentation performance on small objects compared with several recent state-of-the-art models. We test our CaraNet on brain tumor (BraTS 2018) and polyp (Kvasir-SEG, CVC-ColonDB, CVC-ClinicDB, CVC-300, and ETIS-LaribPolypDB) segmentation datasets. Our CaraNet achieves the top-rank mean Dice segmentation accuracy, and results show a distinct advantage of CaraNet in the segmentation of small medical objects.


title: Coherent Ising Machines with Optical Error Correction Circuits

abstract: We propose a network of open-dissipative quantum oscillators with optical error correction circuits. In the proposed network, the squeezed/anti-squeezed vacuum states of the constituent optical parametric oscillators below the threshold establish quantum correlations through optical mutual coupling, while collective symmetry breaking is induced above the threshold as a decision-making process. This initial search process is followed by a chaotic solution search step facilitated by the optical error correction feedback. As an optical hardware technology, the proposed coherent Ising machine (CIM) has several unique features, such as programmable all-to-all Ising coupling in the optical domain, directional coupling $(J_{ij} \neq J_{ji})$ induced chaotic behavior, and low power operation at room temperature. We study the performance of the proposed CIMs and investigate how the performance scales with different problem sizes. The quantum theory of the proposed CIMs can be used as a heuristic algorithm and efficiently implemented on existing digital platforms. This particular algorithm is derived from the truncated Wigner stochastic differential equation. We find that the various CIMs discussed are effective at solving many problem types, however the optimal algorithm is different depending on the instance. We also find that the proposed optical implementations have the potential for low energy consumption when implemented optically on a thin film LiNbO3 platform.


title: Mixed Modes and Asteroseismic Surface Effects: II. Subgiant Systematics

abstract: Models of solar-like oscillators yield acoustic modes at different frequencies than would be seen in actual stars possessing identical interior structure, due to modelling error near the surface. This asteroseismic "surface term" must be corrected when mode frequencies are used to infer stellar structure. Subgiants exhibit oscillations of mixed acoustic ($p$-mode) and gravity ($g$-mode) character, which defy description by the traditional $p$-mode asymptotic relation. Since nonparametric diagnostics of the surface term rely on this description, they cannot be applied to subgiants directly. In Paper I, we generalised such nonparametric methods to mixed modes, and showed that traditional surface-term corrections only account for mixed-mode coupling to, at best, first order in a perturbative expansion. Here, we apply those results, modelling subgiants using asteroseismic data. We demonstrate that, for grid-based inference of subgiant properties using individual mode frequencies, neglecting higher-order effects of mode coupling in the surface term results in significant systematic differences in the inferred stellar masses, and measurable systematics in other fundamental properties. While these systematics are smaller than those resulting from other choices of model construction, they persist for both parametric and nonparametric formulations of the surface term. This suggests that mode coupling should be fully accounted for when correcting for the surface term in seismic modelling with mixed modes, irrespective of the choice of correction used. The inferred properties of subgiants, in particular masses and ages, also depend on the choice of surface-term correction, in a different manner from both main-sequence and red giant stars.


title: New method to achieve the proper polarization state for a vector vortex coronagraph

abstract: The vector vortex coronagraph (VVC) performance in the laboratory and in ground-based observatories has earned it a spot on the NASA mission concepts HabEx and LUVOIR. The VVC induces a phase ramp through the manipulation of the polarization state. Left- and right-circular polarizations get imprinted a phase ramp of opposite signs, which prevents model-based focal plane wavefront sensing and control strategies in natural light. We thus have to work with a polarization state than ensures circularly polarized light at the VVC mask. However, achieving this polarization state can be non trivial if there are optics that add phase retardance of any kind between the circular polarizer and the focal plane mask. Here we present the method currently used at the Caltech high contrast spectroscopy testbed (HCST) to achieve the proper circular polarization state for a VVC, which only uses the deformable mirror and appropriate rotation of the circular polarizer and analyzer optics. At HCST we achieve raw contrast levels of \tentoe~for broadband light with a VVC.


title: Density Sharpening: Principles and Applications to Discrete Data Analysis

abstract: This article introduces a general statistical modeling principle called "Density Sharpening" and applies it to the analysis of discrete count data. The underlying foundation is based on a new theory of nonparametric approximation and smoothing methods for discrete distributions which play a useful role in explaining and uniting a large class of applied statistical methods. The proposed modeling framework is illustrated using several real applications, from seismology to healthcare to physics.


title: The existence of $\mathbb{F}_q$-primitive points on curves using freeness

abstract: Let $\mathcal C_Q$ be the cyclic group of order $Q$, $n$ a divisor of $Q$ and $r$ a divisor of $Q/n$. We introduce the set of $(r,n)$-free elements of $\mathcal C_Q$ and derive a lower bound for the the number of elements $\theta \in \mathbb F_q$ for which $f(\theta)$ is $(r,n)$-free and $F(\theta)$ is $(R,N)$-free, where $ f, F \in \mathbb F_q[x]$. As an application, we consider the existence of $\mathbb F_q$-primitive points on curves like $y^n=f(x)$ and find, in particular, all the odd prime powers $q$ for which the elliptic curves $y^2=x^3 \pm x$ contain an $\mathbb F_q$-primitive point.


title: Reusable Templates and Guides For Documenting Datasets and Models for Natural Language Processing and Generation: A Case Study of the HuggingFace and GEM Data and Model Cards

abstract: Developing documentation guidelines and easy-to-use templates for datasets and models is a challenging task, especially given the variety of backgrounds, skills, and incentives of the people involved in the building of natural language processing (NLP) tools. Nevertheless, the adoption of standard documentation practices across the field of NLP promotes more accessible and detailed descriptions of NLP datasets and models, while supporting researchers and developers in reflecting on their work. To help with the standardization of documentation, we present two case studies of efforts that aim to develop reusable documentation templates -- the HuggingFace data card, a general purpose card for datasets in NLP, and the GEM benchmark data and model cards with a focus on natural language generation. We describe our process for developing these templates, including the identification of relevant stakeholder groups, the definition of a set of guiding principles, the use of existing templates as our foundation, and iterative revisions based on feedback.


title: BN-NAS: Neural Architecture Search with Batch Normalization

abstract: We present BN-NAS, neural architecture search with Batch Normalization (BN-NAS), to accelerate neural architecture search (NAS). BN-NAS can significantly reduce the time required by model training and evaluation in NAS. Specifically, for fast evaluation, we propose a BN-based indicator for predicting subnet performance at a very early training stage. The BN-based indicator further facilitates us to improve the training efficiency by only training the BN parameters during the supernet training. This is based on our observation that training the whole supernet is not necessary while training only BN parameters accelerates network convergence for network architecture search. Extensive experiments show that our method can significantly shorten the time of training supernet by more than 10 times and shorten the time of evaluating subnets by more than 600,000 times without losing accuracy.


title: Convolutive Prediction for Monaural Speech Dereverberation and Noisy-Reverberant Speaker Separation

abstract: A promising approach for speech dereverberation is based on supervised learning, where a deep neural network (DNN) is trained to predict the direct sound from noisy-reverberant speech. This data-driven approach is based on leveraging prior knowledge of clean speech patterns and seldom explicitly exploits the linear-filter structure in reverberation, i.e., that reverberation results from a linear convolution between a room impulse response (RIR) and a dry source signal. In this work, we propose to exploit this linear-filter structure within a deep learning based monaural speech dereverberation framework. The key idea is to first estimate the direct-path signal of the target speaker using a DNN and then identify signals that are decayed and delayed copies of the estimated direct-path signal, as these can be reliably considered as reverberation. They can be either directly removed for dereverberation, or used as extra features for another DNN to perform better dereverberation. To identify the copies, we estimate the underlying filter (or RIR) by efficiently solving a linear regression problem per frequency in the time-frequency domain. We then modify the proposed algorithm for speaker separation in reverberant and noisy-reverberant conditions. State-of-the-art speech dereverberation and speaker separation results are obtained on the REVERB, SMS-WSJ, and WHAMR! datasets.


title: Precision and accuracy of acoustic gunshot location in an urban environment

abstract: The muzzle blast caused by the discharge of a firearm generates a loud, impulsive sound that propagates away from the shooter in all directions. The location of the source can be computed from time-of-arrival measurements of the muzzle blast on multiple acoustic sensors at known locations, a technique known as multilateration. The multilateration problem is considerably simplified by assuming straight-line propagation in a homogeneous medium, a model for which there are multiple published solutions. Live-fire tests of the ShotSpotter gunshot location system in Pittsburgh, PA were analyzed off-line under several algorithms and geometric constraints to evaluate the accuracy of acoustic multilateration in a forensic context. Best results were obtained using the algorithm due to Mathias, Leonari and Galati under a two-dimensional geometric constraint. Multilateration on random subsets of the participating sensor array show that 96% of shots can be located to an accuracy of 15 m or better when six or more sensors participate in the solution.


title: PnP-3D: A Plug-and-Play for 3D Point Clouds

abstract: With the help of the deep learning paradigm, many point cloud networks have been invented for visual analysis. However, there is great potential for development of these networks since the given information of point cloud data has not been fully exploited. To improve the effectiveness of existing networks in analyzing point cloud data, we propose a plug-and-play module, PnP-3D, aiming to refine the fundamental point cloud feature representations by involving more local context and global bilinear response from explicit 3D space and implicit feature space. To thoroughly evaluate our approach, we conduct experiments on three standard point cloud analysis tasks, including classification, semantic segmentation, and object detection, where we select three state-of-the-art networks from each task for evaluation. Serving as a plug-and-play module, PnP-3D can significantly boost the performances of established networks. In addition to achieving state-of-the-art results on four widely used point cloud benchmarks, we present comprehensive ablation studies and visualizations to demonstrate our approach's advantages. The code will be available at https://github.com/ShiQiu0419/pnp-3d.


title: DRB-GAN: A Dynamic ResBlock Generative Adversarial Network for Artistic Style Transfer

abstract: The paper proposes a Dynamic ResBlock Generative Adversarial Network (DRB-GAN) for artistic style transfer. The style code is modeled as the shared parameters for Dynamic ResBlocks connecting both the style encoding network and the style transfer network. In the style encoding network, a style class-aware attention mechanism is used to attend the style feature representation for generating the style codes. In the style transfer network, multiple Dynamic ResBlocks are designed to integrate the style code and the extracted CNN semantic feature and then feed into the spatial window Layer-Instance Normalization (SW-LIN) decoder, which enables high-quality synthetic images with artistic style transfer. Moreover, the style collection conditional discriminator is designed to equip our DRB-GAN model with abilities for both arbitrary style transfer and collection style transfer during the training stage. No matter for arbitrary style transfer or collection style transfer, extensive experiments strongly demonstrate that our proposed DRB-GAN outperforms state-of-the-art methods and exhibits its superior performance in terms of visual quality and efficiency. Our source code is available at \color{magenta}{\url{https://github.com/xuwenju123/DRB-GAN}}.


title: InfoGram and Admissible Machine Learning

abstract: We have entered a new era of machine learning (ML), where the most accurate algorithm with superior predictive power may not even be deployable, unless it is admissible under the regulatory constraints. This has led to great interest in developing fair, transparent and trustworthy ML methods. The purpose of this article is to introduce a new information-theoretic learning framework (admissible machine learning) and algorithmic risk-management tools (InfoGram, L-features, ALFA-testing) that can guide an analyst to redesign off-the-shelf ML methods to be regulatory compliant, while maintaining good prediction accuracy. We have illustrated our approach using several real-data examples from financial sectors, biomedical research, marketing campaigns, and the criminal justice system.


title: The infuence of compressive lattice deformations on the zone-center energy band properties of zincblende GaN and InN. Hybrid density functional results

abstract: We have investigated the effects of hydrostatic pressure and compressive biaxial strain on the $\Gamma$-point energy states of GaN and InN with zincblende crystal structure via first-principles DFT+HSE06 computation. To correctly reproduce accepted experimental values of the energy band gap of both compounds, the procedure includes modified exchange-correlation fractions in the HSE hybrid functional, changing from the standard value $\alpha=0.25$ to $\alpha=0.43$ (GaN) and $\alpha=0.40$ (InN). Within this environment, the work reports on the variation of conduction and valence band edges as functions of the lattice deformation. In addition, we present fitted expressions describing the change in the band gap and the spin-orbit splitting energy due to unit cell size modification. All these parameters are main input quantities in the description of electronic and optical properties of InN/GaN-based heterostructures.


title: A Hamiltonian model for the macroscopic Maxwell equations using exterior calculus

abstract: A Hamiltonian field theory for the macroscopic Maxwell equations with fully general polarization and magnetization is stated in the language of differential forms. The precise procedure for translating the vector calculus formulation into differential forms is discussed in detail. We choose to distinguish between straight and twisted differential forms so that all integrals be taken over densities (i.e. twisted top forms). This ensures that the duality pairings, which are stated as integrals over densities, are orientation independent. The relationship between functional differentiation with respect to vector fields and with respect to differential forms is established using the chain rule. The theory is developed such that the Poisson bracket is metric and orientation independent with all metric dependence contained in the Hamiltonian. As is typically seen in the exterior calculus formulation of Maxwell's equations, the Hodge star operator plays a key role in modeling the constitutive relations. As a demonstration of the kind of constitutive models this theory accommodates, the paper concludes with several examples.


title: Learning to Cluster via Same-Cluster Queries

abstract: We study the problem of learning to cluster data points using an oracle which can answer same-cluster queries. Different from previous approaches, we do not assume that the total number of clusters is known at the beginning and do not require that the true clusters are consistent with a predefined objective function such as the K-means. These relaxations are critical from the practical perspective and, meanwhile, make the problem more challenging. We propose two algorithms with provable theoretical guarantees and verify their effectiveness via an extensive set of experiments on both synthetic and real-world data.


title: Optical representation of thermal nuclear fluctuation effect on band-gap renormalization

abstract: The bandgap of insulating materials is renormalized in various ways by the electron-phonon interaction owing to the dynamical and quantum fluctuations of nuclei. These fluctuation effects are considered in the perturbative Allen-Heine-Cardona theory using the formulae for the Fan-Migdal and Debye-Waller terms. However, the material dependence is not clear in the formulae. Thus, in this study, we focus on the analytical form of the Debye-Waller term and find that the term can be reformulated using the momentum matrix. In addition, the optical selection rule is found to play a role. For diamond-type materials, the Debye-Waller term can be approximately decomposed into a product of the optical transition energy, the mean square displacement of nuclei, and the dipole transition probability. The decomposition can also be applied with an additional approximation to zinc-blende-type materials, as revealed by our first-principles calculation. The magnitudes of the Debye-Waller term of several materials can thus be estimated using basic physical quantities prior to performing the calculation of the electron-phonon interaction.


title: A mimetic discretization of the macroscopic Maxwell equations in Hamiltonian form

abstract: A mimetic spectral element discretization, utilizing a novel Galerkin projection Hodge star operator, of the macroscopic Maxwell equations in Hamiltonian form is presented. The idea of splitting purely topological and metric dependent quantities is natural in the Hamiltonian modeling framework as the Poisson bracket is metric free with the Hamiltonian containing all metric information. This idea may be incorporated into the mimetic spectral element method by directly discretizing the Poincar\'e duality structure. This "split exterior calculus mimetic spectral element method" yields spatially discretized Maxwell's equations which are Hamiltonian and exactly and strongly conserve Gauss's laws. Moreover, the new discrete Hodge star operator is itself of interest as a partition of the purely topological and metric dependent portions of the Hodge star operator. As a simple test case, the numerical results of applying this method to a one-dimensional version of Maxwell's equations are given.

categories: physics.comp-ph cs.NA math.NA physics.class-ph


title: BOBCAT: Bilevel Optimization-Based Computerized Adaptive Testing

abstract: Computerized adaptive testing (CAT) refers to a form of tests that are personalized to every student/test taker. CAT methods adaptively select the next most informative question/item for each student given their responses to previous questions, effectively reducing test length. Existing CAT methods use item response theory (IRT) models to relate student ability to their responses to questions and static question selection algorithms designed to reduce the ability estimation error as quickly as possible; therefore, these algorithms cannot improve by learning from large-scale student response data. In this paper, we propose BOBCAT, a Bilevel Optimization-Based framework for CAT to directly learn a data-driven question selection algorithm from training data. BOBCAT is agnostic to the underlying student response model and is computationally efficient during the adaptive testing process. Through extensive experiments on five real-world student response datasets, we show that BOBCAT outperforms existing CAT methods (sometimes significantly) at reducing test length.