title: Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements

abstract: We propose a method to mitigate heavy-tailed distributions in fermion Quantum Monte Carlo simulations originating from zeros of the fermion determinant. In this case the second moment of the observables might be not well defined, and we show that by merely changing the synchronization between local updates and computation of observables, one can reduce the prefactor of the heavy-tailed distribution, thus substantially suppressing statistical fluctuations of observables. We also show that the average, or the first moment, is well defined and hence is independent on our measuring scheme. The method is especially suitable for local observables similar to e.g. double occupancy, where the resulting speedup can reach two orders of magnitude. For observables, containing spatial correlators, the speedup is more moderate, but still ranges between five and ten. Our results are independent on the nature of the auxiliary field, discrete or continuous, and pave the way to improve measurement strategies for Hybrid Monte Carlo simulations.


title: Equivariant Quantum Graph Circuits

abstract: We investigate quantum circuits for graph representation learning, and propose equivariant quantum graph circuits (EQGCs), as a class of parameterized quantum circuits with strong relational inductive bias for learning over graph-structured data. Conceptually, EQGCs serve as a unifying framework for quantum graph representation learning, allowing us to define several interesting subclasses which subsume existing proposals. In terms of the representation power, we prove that the studied subclasses of EQGCs are universal approximators for functions over the bounded graph domain. This theoretical perspective on quantum graph machine learning methods opens many directions for further work, and could lead to models with capabilities beyond those of classical approaches. We empirically verify the expressive power of EQGCs through a dedicated experiment on synthetic data, and additionally observe that the performance of EQGCs scales well with the depth of the model and does not suffer from barren plateu issues.


title: The existence of partitioned balanced tournament designs

abstract: E. R. Lamken prove in [5] that there exists a partitioned balanced tournament design of side $n$, PBTD($n$), for $n$ a positive integer, $n \ge 5$, except possibly for $n \in \{9,11,15\}$. In this article, we show the existence of PBTD($n$) for $n \in \{9,11,15\}$. As a consequence, the existence of PBTD($n$) has been completely determined.


title: System-Level Analysis of Full-Duplex Self-Backhauled Millimeter Wave Networks

abstract: Integrated access and backhaul (IAB) facilitates cost-effective deployment of millimeter wave(mmWave) cellular networks through multihop self-backhauling. Full-duplex (FD) technology, particularly for mmWave systems, is a potential means to overcome latency and throughput challenges faced by IAB networks. We derive practical and tractable throughput and latency constraints using queueing theory and formulate a network utility maximization problem to evaluate both FD-IAB and half-duplex(HD)-IAB networks. We use this to characterize the network-level improvements seen when upgrading from conventional HD IAB nodes to FD ones by deriving closed-form expressions for (i) latency gain of FD-IAB over HD-IAB and (ii) the maximum number of hops that a HD- and FD-IAB network can support while satisfying latency and throughput targets. Extensive simulations illustrate that FD-IAB can facilitate reduced latency, higher throughput, deeper networks, and fairer service. Compared to HD-IAB,FD-IAB can improve throughput by 8x and reduce latency by 4x for a fourth-hop user. In fact, upgrading IAB nodes with FD capability can allow the network to support latency and throughput targets that its HD counterpart fundamentally cannot meet. The gains are more profound for users further from the donor and can be achieved even when residual self-interference is significantly above the noise floor.


title: Slope of Dry Granular Materials Surface is Generally Curved

abstract: As if it has become a consensus that the slope of the granular heap is straight, and from that, an angle of repose is defined. However, closer inspection shows that the slopes are not truly straight, instead, more often they have convex shapes, although with small convexity. Similar curvature is observed on the sand dunes' surfaces. We derive equations together with experiments on 10 types of granular materials to show that, in general, the slope of the granular heap is convex while the straight slope is a critical state that rarely occurs. The angle of repose is allowed to vary within a narrow range around a critical value. If the angle of repose is slightly greater (smaller) than the critical value, the slope is convex (concave). The model also answers why the surfaces of granular heaps, including sand dunes, are more often convex rather than concave or straight, and advocate the simultaneous estimation of the coefficient of friction on the grain and the type of packing in the heap. The results might open new understandings related to granular material properties as well as sand transport and geological records in several deserts.


title: Viterbi Decoding of CRES Signals in Project 8

abstract: Cyclotron Radiation Emission Spectroscopy (CRES) is a modern approach for determining charged particle energies via high-precision frequency measurements of the emitted cyclotron radiation. For CRES experiments with gas within the fiducial volume, signal and noise dynamics can be modelled by a hidden Markov model. We introduce a novel application of the Viterbi algorithm in order to derive informational limits on the optimal detection of cyclotron radiation signals in this class of gas-filled CRES experiments, thereby providing concrete limits from which future reconstruction algorithms, as well as detector designs, can be constrained. The validity of the resultant decision rules is confirmed using both Monte Carlo and Project 8 data.


title: Mass changes the diffusion coefficient of particles with ligand-receptor contacts in the overdamped limit

abstract: Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands, that bind reversibly to surface receptors -- heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.


title: The Many Faces of Anger: A Multicultural Video Dataset of Negative Emotions in the Wild (MFA-Wild)

abstract: The portrayal of negative emotions such as anger can vary widely between cultures and contexts, depending on the acceptability of expressing full-blown emotions rather than suppression to maintain harmony. The majority of emotional datasets collect data under the broad label ``anger", but social signals can range from annoyed, contemptuous, angry, furious, hateful, and more. In this work, we curated the first in-the-wild multicultural video dataset of emotions, and deeply explored anger-related emotional expressions by asking culture-fluent annotators to label the videos with 6 labels and 13 emojis in a multi-label framework. We provide a baseline multi-label classifier on our dataset, and show how emojis can be effectively used as a language-agnostic tool for annotation.


title: On Markov chain approximations for computing boundary crossing probabilities of diffusion processes

abstract: We propose a discrete time discrete space Markov chain approximation with a Brownian bridge correction for computing curvilinear boundary crossing probabilities of a general diffusion process on a finite time interval. For broad classes of curvilinear boundaries and diffusion processes, we prove the convergence of the constructed approximations in the form of products of the respective substochastic matrices to the boundary crossing probabilities for the process as the time grid used to construct the Markov chains is getting finer. Numerical results indicate that the convergence rate for the proposed approximation with the Brownian bridge correction is $O(n^{-2})$ in the case of $C^2$-boundaries and a uniform time grid with $n$ steps.


title: Performance Enhancement of Tree-based Friends-of-friend Galaxy-finder for High-resolution Simulations of Galaxy Formation

abstract: Cosmological simulations are useful tools for studying the evolution of galaxies, and it is critical to accurately identify galaxies and their halos from raw simulation data. The friends-of-friend (FoF) algorithm has been widely adopted for this purpose because of its simplicity and expandability to higher dimensions. However, it is cost-inefficient when applied to high-resolution simulations because standard FoF implementation leads to too many distance calculations in dense regions. We confirm this through our exercise of applying the 6-dimensional (6D) FoF galaxy finder code, VELOCIraptor (Elahi et al. 2019), on the NewHorizon simulation (Dubois et al. 2021). The high particle resolution of NewHorizon ($M_{\rm star} \sim 10^4 M_{\odot}$) allows a large central number density ($10^{6}\,{\rm kpc}^{-3}$) for typical galaxies, resulting in a few days to weeks of galaxy searches for just one snapshot. Even worse, we observed a significant decrease in the FoF performance in the high-dimensional 6D searches: "the curse of dimensionality" problem. To overcome these issues, we have developed several implementations that can be readily applied to any tree-based FoF code. They include limiting visits to tree nodes, re-ordering the list of particles for searching neighbor particles, and altering the tree structure. Compared to the run with the original code, the new run with these implementations results in the identical galaxy detection with the ideal performance, $O(N \log{N})$, $N$ being the number of particles in a galaxy -- with a speed gain of a factor of 2700 in 3D or 12 in 6D FoF search.


title: Asymptotic nonlocality in non-Abelian gauge theories

abstract: Asymptotically nonlocal field theories represent a sequence of higher-derivative theories whose limit point is a ghost-free, infinite-derivative theory. Here, we extend previous work on pure scalar and Abelian gauge theories to asymptotically nonlocal non-Abelian theories. In particular, we confirm that there is a limit in which the Lee-Wick spectrum can be decoupled, but where the hierarchy problem is resolved via an emergent nonlocal scale that regulates loop diagrams and that is hierarchically smaller than the lightest Lee-Wick resonance.


title: Further metrics for flash-radiography source spots

abstract: Flash radiography of large hydrodynamic experiments driven by high explosives is a diagnostic method in use at many laboratories. A major contributor to the practical resolution obtained by this point-projection technique is the size of the source spot. A variety of techniques for measuring spot size are in use at the different laboratories, as are different definitions of spot size. Additional metrics in common use at Los Alamos National Laboratory are considered here, and compared with those previously discussed in J. Opt. Soc.


title: $U(1)$ local strings in generalized hybrid metric-Palatini gravity

abstract: In this work we will explore $U(1)$ local cosmic string solutions in the context of the generalized hybrid metric-Palatini theory of gravity in its scalar-tensor representation. Using a general static cylindrically symmetric metric to find the dynamical equations for this particular case, we will simplify the equations by imposing boost invariance along $t$ and $z$ directions. The strings properties are determined by both scalar fields and by the effective potential, function of the scalar fields. While for some forms of the potential, the dynamical equations can be solved exactly, for more general forms of the potential the solutions are found numerically. Several stable string configurations were found, whose basic parameters depend essentially on the effective field potential, and on the boundary conditions.


title: From the simplex to the sphere: Faster constrained optimization using the Hadamard parametrization

abstract: The standard simplex in R^n, also known as the probability simplex, is the set of nonnegative vectors whose entries sum up to 1. They frequently appear as constraints in optimization problems that arise in machine learning, statistics, data science, operations research, and beyond. We convert the standard simplex to the unit sphere and thus transform the corresponding constrained optimization problem into an optimization problem on a simple, smooth manifold. We show that KKT points and strict-saddle points of the minimization problem on the standard simplex all correspond to those of the transformed problem, and vice versa. So, solving one problem is equivalent to solving the other problem. Then, we propose several simple, efficient, and projection-free algorithms using the manifold structure. The equivalence and the proposed algorithm can be extended to optimization problems with unit simplex, weighted probability simplex, or `1-norm sphere constraints. Numerical experiments between the new algorithms and existing ones show the advantages of the new approach


title: Handling missing data when estimating causal effects with Targeted Maximum Likelihood Estimation

abstract: Targeted Maximum Likelihood Estimation (TMLE) is increasingly used for doubly robust causal inference, but how missing data should be handled when using TMLE with data-adaptive approaches is unclear. Based on the Victorian Adolescent Health Cohort Study, we conducted a simulation study to evaluate eight missing data methods in this context: complete-case analysis, extended TMLE incorporating outcome-missingness model, missing covariate missing indicator method, five multiple imputation (MI) approaches using parametric or machine-learning models. Six scenarios were considered, varying in exposure/outcome generation models (presence of confounder-confounder interactions) and missingness mechanisms (whether outcome influenced missingness in other variables and presence of interaction/non-linear terms in missingness models). Complete-case analysis and extended TMLE had small biases when outcome did not influence missingness in other variables. Parametric MI without interactions had large bias when exposure/outcome generation models included interactions. Parametric MI including interactions performed best in bias and variance reduction across all settings, except when missingness models included a non-linear term. When choosing a method to handle missing data in the context of TMLE, researchers must consider the missingness mechanism and, for MI, compatibility with the analysis method. In many settings, a parametric MI approach that incorporates interactions and non-linearities is expected to perform well.


title: Predicting the orientation of adsorbed proteins steered with electric fields using a simple electrostatic model

abstract: Under the most common experimental conditions, the adsorption of proteins to solid surfaces is an spontaneous process that leads to a rather compact layer of randomly oriented molecules. Due to the importance of this process for the development of catalytic surfaces, a number of existing computational and experimental approaches try to predict and control the orientation of such molecules. However, and despite their own advantages, these tend to be either too expensive computationally, or oversimplified, undermining their ability to predict the most appropriate experimental conditions to maximize the catalytic activity of adsorbed proteins. To address this current need, we present an efficient computational approach to model the behavior of proteins near surfaces in the presence of an external electric field, based on continuum electrostatics. Our model can not only estimate the overall affinity of the protein with the surface, but also their most likely orientation as a function of the potential applied. In this way, a rational selection of the potential can be performed to maximize the accessibility of the protein's active site to the solvent. The model relies on the Poisson-Boltzmann equation and was implemented in an extension of the code PyGBe that includes an external electric field, and renders the electrostatic component of the solvation free energy. To demonstrate the feasibility of this technique, we investigate the adsorption of trypsin onto a carbon electrode under potentiostatic conditions both numerically and experimentally. We found that even though the adsorption process is largely dominated by hydrophobic effects, the orientation of trypsin can be controlled through an external potential, influencing the position of the active sites, and resulting in an important change in the catalytic activity of the surface.


title: The Lagrange-D'Alembert Principle in Banach Space

abstract: The Lagrange-D'Alembert Principle is one of the fundamental tools of classical mechanics. We generalize this principle to mechanics-like ODE in Banach spaces. As an application we discuss geodesics in infinite dimensional manifolds and a random ODE with nonholonomic constraint.


title: Skeletal Graph Self-Attention: Embedding a Skeleton Inductive Bias into Sign Language Production

abstract: Recent approaches to Sign Language Production (SLP) have adopted spoken language Neural Machine Translation (NMT) architectures, applied without sign-specific modifications. In addition, these works represent sign language as a sequence of skeleton pose vectors, projected to an abstract representation with no inherent skeletal structure. In this paper, we represent sign language sequences as a skeletal graph structure, with joints as nodes and both spatial and temporal connections as edges. To operate on this graphical structure, we propose Skeletal Graph Self-Attention (SGSA), a novel graphical attention layer that embeds a skeleton inductive bias into the SLP model. Retaining the skeletal feature representation throughout, we directly apply a spatio-temporal adjacency matrix into the self-attention formulation. This provides structure and context to each skeletal joint that is not possible when using a non-graphical abstract representation, enabling fluid and expressive sign language production. We evaluate our Skeletal Graph Self-Attention architecture on the challenging RWTH-PHOENIX-Weather-2014T(PHOENIX14T) dataset, achieving state-of-the-art back translation performance with an 8% and 7% improvement over competing methods for the dev and test sets.


title: Simulation-Based Inference of Strong Gravitational Lensing Parameters

abstract: In the coming years, a new generation of sky surveys, in particular, Euclid Space Telescope (2022), and the Rubin Observatory's Legacy Survey of Space and Time (LSST, 2023) will discover more than 200,000 new strong gravitational lenses, which represents an increase of more than two orders of magnitude compared to currently known sample sizes. Accurate and fast analysis of such large volumes of data under a statistical framework is therefore crucial for all sciences enabled by strong lensing. Here, we report on the application of simulation-based inference methods, in particular, density estimation techniques, to the predictions of the set of parameters of strong lensing systems from neural networks. This allows us to explicitly impose desired priors on lensing parameters, while guaranteeing convergence to the optimal posterior in the limit of perfect performance.


title: Limiting Betti distributions of Hilbert schemes on $n$ points

abstract: Hausel and Rodriguez-Villegas recently observed that work of G\"ottsche, combined with a classical result of Erd\H{o}s and Lehner on integer partitions, implies that the limiting Betti distribution for the Hilbert schemes $(\mathbb{C}^2)^{[n]}$ on $n$ points, as $n\rightarrow +\infty,$ is a \textit{Gumbel distribution}. In view of this example, they ask for further such Betti distributions. We answer this question for the quasihomogeneous Hilbert schemes $((\mathbb{C}^2)^{[n]})^{T_{\alpha,\beta}}$ that are cut out by torus actions. We prove that their limiting distributions are also of Gumbel type. To obtain this result, we combine work of Buryak, Feigin, and Nakajima on these Hilbert schemes with our generalization of the result of Erd\H{o}s and Lehner, which gives the distribution of the number of parts in partitions that are multiples of a fixed integer $A\geq 2.$ Furthermore, if $p_k(A;n)$ denotes the number of partitions of $n$ with exactly $k$ parts that are multiples of $A$, then we obtain the asymptotic $$ p_k(A,n)\sim \frac{24^{\frac k2-\frac14}(n-Ak)^{\frac k2-\frac34}}{\sqrt2\left(1-\frac1A\right)^{\frac k2-\frac14}k!A^{k+\frac12}(2\pi)^k}e^{2\pi\sqrt{\frac1{6}\left(1-\frac1A\right)(n-Ak)}}, $$ a result which is of independent interest.


title: Long-Range Thermal 3D Perception in Low Contrast Environments

abstract: This report discusses the results of SBIR Phase I effort to prove the feasibility of dramatic improvement of the microbolometer-based Long Wave Infrared (LWIR) detectors sensitivity, especially for the 3D measurements. The resulting low SWaP-C thermal depth-sensing system will enable the situational awareness of Autonomous Air Vehicles for Advanced Air Mobility (AAM). It will provide robust 3D information of the surrounding environment, including low-contrast static and moving objects, at far distances in degraded visual conditions and GPS-denied areas. Our multi-sensor 3D perception enabled by COTS uncooled thermal sensors mitigates major weakness of LWIR sensors - low contrast by increasing the system sensitivity over an order of magnitude. There were no available thermal image sets suitable for evaluating this technology, making datasets acquisition our first goal. We discuss the design and construction of the prototype system with sixteen 640pix x 512pix LWIR detectors, camera calibration to subpixel resolution, capture, and process synchronized image. The results show the 3.84x contrast increase for intrascene-only data and an additional 5.5x - with the interscene accumulation, reaching system noise-equivalent temperature difference (NETD) of 1.9 mK with the 40 mK sensors.


title: Expected value of letters of permutations with a given number of $k$-cycles

abstract: In this paper, we study permutations $\pi \in S_n$ with exactly $m$ transpositions. In particular, we are interested in the expected value of $\pi(1)$ when such permutations are chosen uniformly at random. When $n$ is even, this expected value is approximated closely by $(n+1)/2$, with an error term that is related to the number isometries of the $(n/2-m)$-dimensional hypercube that move every face. Furthermore, when $k \mid n$, this construction generalizes to allow us to compute the expected value of $\pi(1)$ for permutations with exactly $m$ $k$-cycles. In this case, the expected value has an error term which is related instead to the number derangements of the generalized symmetric group $S(k,n/k-m)$. When $k$ does not divide $n$, the expected value of $\pi(1)$ is precisely $(n+1)/2$. Indirectly, this suggests the existence of a reversible algorithm to insert a letter into a permutation which preserves the number of $k$-cycles, which we construct.


title: RamBoAttack: A Robust Query Efficient Deep Neural Network Decision Exploit

abstract: Machine learning models are critically susceptible to evasion attacks from adversarial examples. Generally, adversarial examples, modified inputs deceptively similar to the original input, are constructed under whitebox settings by adversaries with full access to the model. However, recent attacks have shown a remarkable reduction in query numbers to craft adversarial examples using blackbox attacks. Particularly, alarming is the ability to exploit the classification decision from the access interface of a trained model provided by a growing number of Machine Learning as a Service providers including Google, Microsoft, IBM and used by a plethora of applications incorporating these models. The ability of an adversary to exploit only the predicted label from a model to craft adversarial examples is distinguished as a decision-based attack. In our study, we first deep dive into recent state-of-the-art decision-based attacks in ICLR and SP to highlight the costly nature of discovering low distortion adversarial employing gradient estimation methods. We develop a robust query efficient attack capable of avoiding entrapment in a local minimum and misdirection from noisy gradients seen in gradient estimation methods. The attack method we propose, RamBoAttack, exploits the notion of Randomized Block Coordinate Descent to explore the hidden classifier manifold, targeting perturbations to manipulate only localized input features to address the issues of gradient estimation methods. Importantly, the RamBoAttack is more robust to the different sample inputs available to an adversary and the targeted class. Overall, for a given target class, RamBoAttack is demonstrated to be more robust at achieving a lower distortion within a given query budget. We curate our extensive results using the large-scale high-resolution ImageNet dataset and open-source our attack, test samples and artifacts on GitHub.


title: Relaxing inflation models with non-minimal coupling: A general study

abstract: We investigate the predictions of inflation models with a non-minimal coupling to gravity for inflationary observables such as the spectral index and tensor-to-scalar ratio in a general setting. We argue that, depending on the relation between the Jordan frame potential and a function characterizing the non-minimal coupling, one can classify the model into three categories each of which gives distinctive predictions for the inflationary observables. We derive general predictions for each class and also investigate some explicit models to discuss how the general features arise. Our results would be useful to design an inflation model consistent with observational constraints with a non-minimal coupling to gravity.


title: Singular jets in compound drop impact

abstract: Compound drop impacting on a solid surface is of considerable importance in industrial applications, such as combustion, food industry, and drug encapsulation. An intriguing phenomenon associated with this process is the occurrence of singular jets that are up to dozens of times faster than the impact velocity. These jets break into micro-droplets, which can produce aerosols and affect the quality of printing technologies. Here, we investigate experimentally and numerically the jetting process after a coaxial water-in-oil compound drop impacts on a glass substrate with different releasing heights and volumetric ratios. After impact, the water core spreads and retracts, giving rise to a vertical jet initially made of oil. For certain values of the impacting velocity, high speed and very thin jets are observed, the so-called singular jets. Depending on the volumetric ratio, one or two velocity peaks can be observed when varying the impact velocity, triggered by the contraction dynamics of a deep and cylindrical cavity. The self-similar time-evolution of the collapse for the first singularity regime follows a 1/2 power law in time, which can be derived from bubble pinch-off. In contrast, the collapse at the second peak follows a 2/3 power law, which can be accounted for by a balance between inertial and capillary forces.


title: Well-posedness for the free boundary hard phase model in general relativity

abstract: The hard phase model describes a relativistic barotropic and irrotational fluid with sound speed equal to the speed of light. In the framework of general relativity, the fluid, as a matter field, affects the geometry of the background spacetime. Therefore the motion of the fluid must be coupled to the Einstein equations which describe the structure of the underlying spacetime. In this work we prove a priori estimates and well-posedness in Sobolev spaces for this model with free boundary. Estimates for the curvature are derived using the Bianchi equations in a frame that is parallel transported by the fluid velocity. The fluid velocity is also decomposed with respect to this parallel frame, and its components are estimated using a coupled interior-boundary system of wave equations.


title: SmartCon: Deep Probabilistic Learning Based Intelligent Link-Configuration in Narrowband-IoT Towards 5G and B5G

abstract: To enhance the coverage and transmission reliability, repetitions adopted by Narrowband Internet of Things (NB-IoT) allow repeating transmissions several times. However, this results in a waste of radio resources when the signal strength is high. In addition, in low signal quality, the selection of a higher modulation and coding scheme (MCS) level leads to a huge packet loss in the network. Moreover, the number of physical resource blocks (PRBs) per-user needs to be chosen dynamically, such that the utilization of radio resources can be improved on per-user basis. Therefore, in NB-IoT systems, dynamic adaptation of repetitions, MCS, and radio resources, known as auto link-configuration, is crucial. Accordingly, in this paper, we propose SmartCon which is a Generative Adversarial Network (GAN)-based deep learning approach for auto link-configuration during uplink or downlink scheduling, such that the packet loss rate is significantly reduced in NB-IoT networks. For the training purpose of the GAN, we use a Multi-Armed Bandit (MAB)-based reinforcement learning mechanism that intelligently tunes its output depending on the present network condition. The performance of SmartCon is thoroughly evaluated through simulations where it is shown to significantly improve the performance of NB-IoT systems compared to baseline schemes.


title: Polarization rotator Bragg grating assisted wavelength selective polarization alignment

abstract: A device performing both the polarization-alignment and the wavelength selection is demonstrated by fabricated device. Waveguides incorporating Bragg gratings for polarization rotation and mode conversion having different widths are placed near to each other. For the polarization rotation, experimental results show that the field strength ratio in the two waveguides defines the diffraction coupling coefficient and excitation strength of unwanted modes. The excitation of unwanted mode is reduced, when a large radius short curved waveguide is used for input/output waveguides as shown in the experiment. Polarization aligned wavelength peaks were observed at the backward drop port in the experiment.


title: Solving linear Bayesian inverse problems using a fractional total variation-Gaussian (FTG) prior and transport map

abstract: The Bayesian inference is widely used in many scientific and engineering problems, especially in the linear inverse problems in infinite-dimensional setting where the unknowns are functions. In such problems, choosing an appropriate prior distribution is an important task. In particular, when the function to infer has much detail information, such as many sharp jumps, corners, and the discontinuous and nonsmooth oscillation, the so-called total variation-Gaussian (TG) prior is proposed in function space to address it. However, the TG prior is easy to lead the blocky (staircase) effect in numerical results. In this work, we present a fractional order-TG (FTG) hybrid prior to deal with such problems, where the fractional order total variation (FTV) term is used to capture the detail information of the unknowns and simultaneously uses the Gaussian measure to ensure that it results in a well-defined posterior measure. For the numerical implementations of linear inverse problems in function spaces, we also propose an efficient independence sampler based on a transport map, which uses a proposal distribution derived from a diagonal map, and the acceptance probability associated to the proposal is independent of discretization dimensionality. And in order to take full advantage of the transport map, the hierarchical Bayesian framework is applied to flexibly determine the regularization parameter. Finally we provide some numerical examples to demonstrate the performance of the FTG prior and the efficiency and robustness of the proposed independence sampler method.


title: From Noether's Theorem to Bremsstrahlung: a pedagogical introduction to large gauge transformations and classical soft theorems

abstract: Electromagnetism contains an infinite dimensional symmetry group of large gauge transformations. This gives rise to an infinite number of conserved quantities called "soft charges" via Noether's theorem. When charged particles scatter, the conservation of soft charge constrains the overall amount of radiation emitted per angle. Here we describe the physical consequences of soft charge conservation and give fresh accounts of the roles of spacelike and timelike infinity in these conservation laws. We conclude by exploring the possibility of creating a dual boundary theory of electromagnetism.


title: Image-to-Image Translation-based Data Augmentation for Robust EV Charging Inlet Detection

abstract: This work addresses the task of electric vehicle (EV) charging inlet detection for autonomous EV charging robots. Recently, automated EV charging systems have received huge attention to improve users' experience and to efficiently utilize charging infrastructures and parking lots. However, most related works have focused on system design, robot control, planning, and manipulation. Towards robust EV charging inlet detection, we propose a new dataset (EVCI dataset) and a novel data augmentation method that is based on image-to-image translation where typical image-to-image translation methods synthesize a new image in a different domain given an image. To the best of our knowledge, the EVCI dataset is the first EV charging inlet dataset. For the data augmentation method, we focus on being able to control synthesized images' captured environments (e.g., time, lighting) in an intuitive way. To achieve this, we first propose the environment guide vector that humans can intuitively interpret. We then propose a novel image-to-image translation network that translates a given image towards the environment described by the vector. Accordingly, it aims to synthesize a new image that has the same content as the given image while looking like captured in the provided environment by the environment guide vector. Lastly, we train a detection method using the augmented dataset. Through experiments on the EVCI dataset, we demonstrate that the proposed method outperforms the state-of-the-art methods. We also show that the proposed method is able to control synthesized images using an image and environment guide vectors.


title: LCTR: On Awakening the Local Continuity of Transformer for Weakly Supervised Object Localization

abstract: Weakly supervised object localization (WSOL) aims to learn object localizer solely by using image-level labels. The convolution neural network (CNN) based techniques often result in highlighting the most discriminative part of objects while ignoring the entire object extent. Recently, the transformer architecture has been deployed to WSOL to capture the long-range feature dependencies with self-attention mechanism and multilayer perceptron structure. Nevertheless, transformers lack the locality inductive bias inherent to CNNs and therefore may deteriorate local feature details in WSOL. In this paper, we propose a novel framework built upon the transformer, termed LCTR (Local Continuity TRansformer), which targets at enhancing the local perception capability of global features among long-range feature dependencies. To this end, we propose a relational patch-attention module (RPAM), which considers cross-patch information on a global basis. We further design a cue digging module (CDM), which utilizes local features to guide the learning trend of the model for highlighting the weak local responses. Finally, comprehensive experiments are carried out on two widely used datasets, ie, CUB-200-2011 and ILSVRC, to verify the effectiveness of our method.


title: Fluctuations of Power versus Energy for Random Fields Near a Perfectly Conducting Boundary

abstract: The standard deviations of the energy and Poynting power densities for an isotropic random field near a perfectly conducting planar boundary are characterized, based on quartic plane-wave expansions. For normal and transverse components, different rates of decay exist as a function of electrical distance from the boundary. At large distances, the envelopes for the power are more strongly damped than for the energy, both showing inverse power law decay. The decay for the standard deviation is generally one order faster than for the corresponding mean. For the normally directed power flux, its standard deviation near the boundary increases linearly with distance. The relative uncertainty of the scalar power is much smaller than for the Poynting power. Poynting's theorem for standard deviations is obtained and demonstrates larger standard deviations of the energy imbalance and power flux than their mean values.


title: Gillespie algorithms for stochastic multiagent dynamics in populations and network

abstract: Many multiagent dynamics, including various collective dynamics occurring on networks, can be modeled as a stochastic process in which the agents in the system change their state over time in interaction with each other. The Gillespie algorithms are popular algorithms that exactly simulate such stochastic multiagent dynamics when each state change is driven by a discrete event, the dynamics is defined in continuous time, and the stochastic law of event occurrence is governed by independent Poisson processes. In the first main part of this volume, we provide a tutorial on the Gillespie algorithms focusing on simulation of social multiagent dynamics occurring in populations and networks. We do not assume advanced knowledge of mathematics (or computer science or physics). We clarify why one should use the continuous-time models and the Gillespie algorithms in many cases, instead of easier-to-understand discrete-time models. In the remainder of this volume, we review recent extensions of the Gillespie algorithms aiming to add more reality to the model (i.e., non-Poissonian cases) or to speed up the simulations.


title: Motion by crystalline-like mean curvature: a survey

abstract: We consider a class of anisotropic curvature flows called a crystalline curvature flow. We present a survey on this class of flows with special emphasis on the well-posedness of its initial value problem.


title: 3D Scene Understanding at Urban Intersection using Stereo Vision and Digital Map

abstract: The driving behavior at urban intersections is very complex. It is thus crucial for autonomous vehicles to comprehensively understand challenging urban traffic scenes in order to navigate intersections and prevent accidents. In this paper, we introduce a stereo vision and 3D digital map based approach to spatially and temporally analyze the traffic situation at urban intersections. Stereo vision is used to detect, classify and track obstacles, while a 3D digital map is used to improve ego-localization and provide context in terms of road-layout information. A probabilistic approach that temporally integrates these geometric, semantic, dynamic and contextual cues is presented. We qualitatively and quantitatively evaluate our proposed technique on real traffic data collected at an urban canyon in Tokyo to demonstrate the efficacy of the system in providing comprehensive awareness of the traffic surroundings.


title: On the Dilution of Precision for Time Difference of Arrival with Station Deployment

abstract: The paper aims to reveal the relationship between the performance of moving object tracking algorithms and the tracking anchors (station) deployment. The Dilution of Precision (DoP) for Time difference of arrival (TDoA) technique with respect to anchor deployment is studied. Linear estimator and non-linear estimator are used for TDoA algorithms. The research findings are: for the linear estimator, the DoP attain a lower value when other anchors are scattered around a central anchor; for the nonlinear estimator, the DoP is optimal when the anchors are scattered around the target tag. Experiments on both of the algorithms are conducted, targeting the location precision related to the anchors' deployment, with practical situations for tracking moving objects integrated with a Kalman Filter (KF) in an Ultra-Wide Band (UWB) based real-time localization system. The work provides a guideline for deploying anchors in UWB-based tracking systems.


title: Quandle Coloring Quivers of (p, 2)-Torus Links

abstract: A quandle coloring quiver is a quiver structure, introduced by Karina Cho and Sam Nelson, which is defined on the set of quandle colorings of an oriented knot or link by a finite quandle. We study quandle coloring quivers of (p, 2)-torus knots and links with respect to dihedral quandles.


title: IFR-Explore: Learning Inter-object Functional Relationships in 3D Indoor Scenes

abstract: Building embodied intelligent agents that can interact with 3D indoor environments has received increasing research attention in recent years. While most works focus on single-object or agent-object visual functionality and affordances, our work proposes to study a new kind of visual relationship that is also important to perceive and model -- inter-object functional relationships (e.g., a switch on the wall turns on or off the light, a remote control operates the TV). Humans often spend little or no effort to infer these relationships, even when entering a new room, by using our strong prior knowledge (e.g., we know that buttons control electrical devices) or using only a few exploratory interactions in cases of uncertainty (e.g., multiple switches and lights in the same room). In this paper, we take the first step in building AI system learning inter-object functional relationships in 3D indoor environments with key technical contributions of modeling prior knowledge by training over large-scale scenes and designing interactive policies for effectively exploring the training scenes and quickly adapting to novel test scenes. We create a new benchmark based on the AI2Thor and PartNet datasets and perform extensive experiments that prove the effectiveness of our proposed method. Results show that our model successfully learns priors and fast-interactive-adaptation strategies for exploring inter-object functional relationships in complex 3D scenes. Several ablation studies further validate the usefulness of each proposed module.


title: Zero-Shot Uncertainty-Aware Deployment of Simulation Trained Policies on Real-World Robots

abstract: While deep reinforcement learning (RL) agents have demonstrated incredible potential in attaining dexterous behaviours for robotics, they tend to make errors when deployed in the real world due to mismatches between the training and execution environments. In contrast, the classical robotics community have developed a range of controllers that can safely operate across most states in the real world given their explicit derivation. These controllers however lack the dexterity required for complex tasks given limitations in analytical modelling and approximations. In this paper, we propose Bayesian Controller Fusion (BCF), a novel uncertainty-aware deployment strategy that combines the strengths of deep RL policies and traditional handcrafted controllers. In this framework, we can perform zero-shot sim-to-real transfer, where our uncertainty based formulation allows the robot to reliably act within out-of-distribution states by leveraging the handcrafted controller while gaining the dexterity of the learned system otherwise. We show promising results on two real-world continuous control tasks, where BCF outperforms both the standalone policy and controller, surpassing what either can achieve independently. A supplementary video demonstrating our system is provided at https://bit.ly/bcf_deploy.


title: Representing 3D Shapes with Probabilistic Directed Distance Fields

abstract: Differentiable rendering is an essential operation in modern vision, allowing inverse graphics approaches to 3D understanding to be utilized in modern machine learning frameworks. Explicit shape representations (voxels, point clouds, or meshes), while relatively easily rendered, often suffer from limited geometric fidelity or topological constraints. On the other hand, implicit representations (occupancy, distance, or radiance fields) preserve greater fidelity, but suffer from complex or inefficient rendering processes, limiting scalability. In this work, we endeavour to address both shortcomings with a novel shape representation that allows fast differentiable rendering within an implicit architecture. Building on implicit distance representations, we define Directed Distance Fields (DDFs), which map an oriented point (position and direction) to surface visibility and depth. Such a field can render a depth map with a single forward pass per pixel, enable differential surface geometry extraction (e.g., surface normals and curvatures) via network derivatives, be easily composed, and permit extraction of classical unsigned distance fields. Using probabilistic DDFs (PDDFs), we show how to model inherent discontinuities in the underlying field. Finally, we apply our method to fitting single shapes, unpaired 3D-aware generative image modelling, and single-image 3D reconstruction tasks, showcasing strong performance with simple architectural components via the versatility of our representation.


title: Self-Ensemling for 3D Point Cloud Domain Adaption

abstract: Recently 3D point cloud learning has been a hot topic in computer vision and autonomous driving. Due to the fact that it is difficult to manually annotate a qualitative large-scale 3D point cloud dataset, unsupervised domain adaptation (UDA) is popular in 3D point cloud learning which aims to transfer the learned knowledge from the labeled source domain to the unlabeled target domain. However, the generalization and reconstruction errors caused by domain shift with simply-learned model are inevitable which substantially hinder the model's capability from learning good representations. To address these issues, we propose an end-to-end self-ensembling network (SEN) for 3D point cloud domain adaption tasks. Generally, our SEN resorts to the advantages of Mean Teacher and semi-supervised learning, and introduces a soft classification loss and a consistency loss, aiming to achieve consistent generalization and accurate reconstruction. In SEN, a student network is kept in a collaborative manner with supervised learning and self-supervised learning, and a teacher network conducts temporal consistency to learn useful representations and ensure the quality of point clouds reconstruction. Extensive experiments on several 3D point cloud UDA benchmarks show that our SEN outperforms the state-of-the-art methods on both classification and segmentation tasks. Moreover, further analysis demonstrates that our SEN also achieves better reconstruction results.


title: Realized GARCH, CBOE VIX, and the Volatility Risk Premium

abstract: We show that the Realized GARCH model yields close-form expression for both the Volatility Index (VIX) and the volatility risk premium (VRP). The Realized GARCH model is driven by two shocks, a return shock and a volatility shock, and these are natural state variables in the stochastic discount factor (SDF). The volatility shock endows the exponentially affine SDF with a compensation for volatility risk. This leads to dissimilar dynamic properties under the physical and risk-neutral measures that can explain time-variation in the VRP. In an empirical application with the S&P 500 returns, the VIX, and the VRP, we find that the Realized GARCH model significantly outperforms conventional GARCH models.


title: Surrogate-based cross-correlation for particle image velocimetry

abstract: This paper presents a novel surrogate-based cross-correlation (SBCC) framework to improve the correlation performance for practical particle image velocimetry~(PIV). The basic idea is that an optimized surrogate filter/image, replacing one raw image, will produce a more accurate and robust correlation signal. Specifically, the surrogate image is encouraged to generate perfect Gaussian-shaped correlation map to tracking particles (PIV image pair) while producing zero responses to image noise (context images). And the problem is formularized with an objective function composed of surrogate loss and consistency loss. As a result, the closed-form solution provides an efficient multivariate operator that could consider other negative context images. Compared with the state-of-the-art baseline methods (background subtraction, robust phase correlation, etc.), our SBCC method exhibits significant performance improvement (accuracy and robustness) on the synthetic dataset and several challenging experimental PIV cases. Besides, our implementation with experimental details (\url{https://github.com/yongleex/SBCC}) is also available for interested researchers.


title: Inferring Invariants with Quantifier Alternations: Taming the Search Space Explosion

abstract: We present a PDR/IC3 algorithm for finding inductive invariants with quantifier alternations. We tackle scalability issues that arise due to the large search space of quantified invariants by combining a breadth-first search strategy and a new syntactic form for quantifier-free bodies. The breadth-first strategy prevents inductive generalization from getting stuck in regions of the search space that are expensive to search and focuses instead on lemmas that are easy to discover. The new syntactic form is well-suited to lemmas with quantifier alternations by allowing both limited conjunction and disjunction in the quantifier-free body, while carefully controlling the size of the search space. Combining the breadth-first strategy with the new syntactic form results in useful inductive bias by prioritizing lemmas according to: (i) well-defined syntactic metrics for simple quantifier structures and quantifier-free bodies, and (ii) the empirically useful heuristic of preferring lemmas that are fast to discover. On a benchmark suite of primarily distributed protocols and complex Paxos variants, we demonstrate that our algorithm can solve more of the most complicated examples than state-of-the-art techniques.


title: Thick embeddings of graphs into symmetric spaces via coarse geometry

abstract: We prove estimates for the optimal volume of thick embeddings of finite graphs into symmetric spaces, generalising results of Kolmogorov-Barzdin and Gromov-Guth for embeddings into Euclidean spaces. We distinguish two very different behaviours depending on the rank of the non-compact factor. For rank at least 2, we construct thick embeddings of $N$-vertex graphs with volume $CN\ln(1+N)$ and prove that this is optimal. For rank at most $1$ we prove lower bounds of the form $cN^a$ for some (explicit) $a>1$ which depends on the dimension of the Euclidean factor and the conformal dimension of the boundary of the non-compact factor. The main tool is a coarse geometric analogue of a thick embedding called a coarse wiring, with the key property that the minimal volume of a thick embedding is comparable to the ``minimal volume'' of a coarse wiring for symmetric spaces of dimension at least $3$. In the appendix it is proved that for each $k\geq 3$ every bounded degree graph admits a coarse wiring into $\mathbb{R}^k$ with volume at most $CN^{1+\frac{1}{k-1}}$. As a corollary, the same upper bound holds for real hyperbolic space of dimension $k+1$ and in both cases this result is optimal.


title: Eigenfunction and eigenmode-spacing statistics in chaotic photonic crystal graphs

abstract: The statistical properties of wave chaotic systems of varying dimensionalities and realizations have been studied extensively. These systems are commonly characterized by the statistics of the eigenmode-spacings and the statistics of the eigenfunctions. Here, we propose photonic crystal (PC) defect waveguide graphs as a new physical setting for chaotic graph studies. Photonic crystal waveguides possess a dispersion relation for the propagating modes which is engineerable. Graphs constructed by joining these waveguides possess junctions and bends with distinct scattering properties. We present numerically determined statistical properties of an ensemble of such PC-graphs including both eigenfunction amplitude and eigenmode-spacing studies. Our proposed system is compatible with silicon nanophotonic technology and opens chaotic graph studies to a new community of researchers.


title: Are We There Yet? Timing and Floating-Point Attacks on Differential Privacy Systems

abstract: Differential privacy is a de facto privacy framework that has seen adoption in practice via a number of mature software platforms. Implementation of differentially private (DP) mechanisms has to be done carefully to ensure end-to-end security guarantees. In this paper we study two implementation flaws in the noise generation commonly used in DP systems. First we examine the Gaussian mechanism's susceptibility to a floating-point representation attack. The premise of this first vulnerability is similar to the one carried out by Mironov in 2011 against the Laplace mechanism. Our experiments show attack's success against DP algorithms, including deep learning models trained using differentially-private stochastic gradient descent. In the second part of the paper we study discrete counterparts of the Laplace and Gaussian mechanisms that were previously proposed to alleviate the shortcomings of floating-point representation of real numbers. We show that such implementations unfortunately suffer from another side channel: a novel timing attack. An observer that can measure the time to draw (discrete) Laplace or Gaussian noise can predict the noise magnitude, which can then be used to recover sensitive attributes. This attack invalidates differential privacy guarantees of systems implementing such mechanisms. We demonstrate that several commonly used, state-of-the-art implementations of differential privacy are susceptible to these attacks. We report success rates up to 92.56% for floating point attacks on DP-SGD, and up to 99.65% for end-to-end timing attacks on private sum protected with discrete Laplace. Finally, we evaluate and suggest partial mitigations.


title: Option Pricing with State-dependent Pricing Kernel

abstract: We introduce a new volatility model for option pricing that combines Markov switching with the Realized GARCH framework. This leads to a novel pricing kernel with a state-dependent variance risk premium and a pricing formula for European options, which is derived with an analytical approximation method. We apply the Markov switching Realized GARCH model to S&P 500 index options from 1990 to 2019 and find that investors' aversion to volatility-specific risk is time-varying. The proposed framework outperforms competing models and reduces (in-sample and out-of-sample) option pricing errors by 15% or more.


title: A Large-Scale Benchmark for the Incompressible Navier-Stokes Equations

abstract: We introduce a collection of benchmark problems in 2D and 3D (geometry description and boundary conditions), including simple cases with known analytic solution, classical experimental setups, and complex geometries with fabricated solutions for evaluation of numerical schemes for incompressible Navier-Stokes equations in laminar flow regime. We compare the performance of a representative selection of most broadly used algorithms for Navier-Stokes equations on this set of problems. Where applicable, we compare the most common spatial discretization choices (unstructured triangle/tetrahedral meshes and structured or semi-structured quadrilateral/hexahedral meshes). The study shows that while the type of spatial discretization used has a minor impact on the accuracy of the solutions, the choice of time integration method, spatial discretization order, and the choice of solving the coupled equations or reducing them to simpler subproblems have very different properties. Methods that are directly solving the original equations tend to be more accurate than splitting approaches for the same number of degrees of freedom, but numerical or computational difficulty arise when they are scaled to larger problem sizes. Low-order splitting methods are less accurate, but scale more easily to large problems, while higher-order splitting methods are accurate but require dense time discretizations to be stable. We release the description of the experiments and an implementation of our benchmark, which we believe will enable statistically significant comparisons with the state of the art as new approaches for solving the incompressible Navier-Stokes equations are introduced.


title: Robustness Certificates for Implicit Neural Networks: A Mixed Monotone Contractive Approach

abstract: Implicit neural networks are a general class of learning models that replace the layers in traditional feedforward models with implicit algebraic equations. Compared to traditional learning models, implicit networks offer competitive performance and reduced memory consumption. However, they can remain brittle with respect to input adversarial perturbations. This paper proposes a theoretical and computational framework for robustness verification of implicit neural networks; our framework blends together mixed monotone systems theory and contraction theory. First, given an implicit neural network, we introduce a related embedded network and show that, given an $\ell_\infty$-norm box constraint on the input, the embedded network provides an $\ell_\infty$-norm box overapproximation for the output of the given network. Second, using $\ell_{\infty}$-matrix measures, we propose sufficient conditions for well-posedness of both the original and embedded system and design an iterative algorithm to compute the $\ell_{\infty}$-norm box robustness margins for reachability and classification problems. Third, of independent value, we propose a novel relative classifier variable that leads to tighter bounds on the certified adversarial robustness in classification problems. Finally, we perform numerical simulations on a Non-Euclidean Monotone Operator Network (NEMON) trained on the MNIST dataset. In these simulations, we compare the accuracy and run time of our mixed monotone contractive approach with the existing robustness verification approaches in the literature for estimating the certified adversarial robustness.


title: Adaptive projected SOR algorithms for nonnegative quadratic programming

abstract: The optimal value of the projected successive overrelaxation (PSOR) method for nonnegative quadratic programming problems is problem-dependent. We present a novel adaptive PSOR algorithm that adaptively controls the relaxation parameter using the Wolfe conditions. The method and its variants can be applied to various problems without requiring a specific assumption regarding the matrix defining the objective function, and the cost for updating the parameter is negligible in the whole iteration. Numerical experiments show that the proposed methods often perform comparably to (or sometimes superior to) the PSOR method with a nearly optimal relaxation parameter.


title: First application of mass measurement with the Rare-RI Ring reveals the solar r-process abundance trend at A=122 and A=123

abstract: The Rare-RI Ring (R3) is a recently commissioned cyclotron-like storage ring mass spectrometer dedicated to mass measurements of exotic nuclei far from stability at Radioactive Isotope Beam Factory (RIBF) in RIKEN. The first application of mass measurement using the R3 mass spectrometer at RIBF is reported. Rare isotopes produced at RIBF, $^{127}$Sn, $^{126}$In, $^{125}$Cd, $^{124}$Ag, $^{123}$Pd, were injected in R3. Masses of $^{126}$In, $^{125}$Cd, and $^{123}$Pd were measured whereby the mass uncertainty of $^{123}$Pd was improved. This is the first reported measurement with a new storage ring mass spectrometery technique realized at a heavy-ion cyclotron and employing individual injection of the pre-identified rare nuclei. The latter is essential for the future mass measurements of the rarest isotopes produced at RIBF. The impact of the new $^{123}$Pd result on the solar $r$-process abundances in a neutron star merger event is investigated by performing reaction network calculations of 20 trajectories with varying electron fraction $Y_e$. It is found that the neutron capture cross section on $^{123}$Pd increases by a factor of 2.2 and $\beta$-delayed neutron emission probability, $P_\mathrm{1n}$, of $^{123}$Rh increases by 14\%. The neutron capture cross section on $^{122}$Pd decreases by a factor of 2.6 leading to pileup of material at $A=122$, thus reproducing the trend of the solar $r$-process abundances. The trend of the two-neutron separation energies (S$_\mathrm{2n}$) was investigated for the Pd isotopic chain. The new mass measurement with improved uncertainty excludes large changes of the S$_\mathrm{2n}$ value at $N=77$. Such large increase of the S$_\mathrm{2n}$ values before $N=82$ was proposed as an alternative to the quenching of the $N=82$ shell gap to reproduce $r$-process abundances in the mass region of $A=112-124$.


title: Building Autocorrelation-Aware Representations for Fine-Scale Spatiotemporal Prediction

abstract: Many scientific prediction problems have spatiotemporal data- and modeling-related challenges in handling complex variations in space and time using only sparse and unevenly distributed observations. This paper presents a novel deep learning architecture, Deep learning predictions for LocATion-dependent Time-sEries data (DeepLATTE), that explicitly incorporates theories of spatial statistics into neural networks to address these challenges. In addition to a feature selection module and a spatiotemporal learning module, DeepLATTE contains an autocorrelation-guided semi-supervised learning strategy to enforce both local autocorrelation patterns and global autocorrelation trends of the predictions in the learned spatiotemporal embedding space to be consistent with the observed data, overcoming the limitation of sparse and unevenly distributed observations. During the training process, both supervised and semi-supervised losses guide the updates of the entire network to: 1) prevent overfitting, 2) refine feature selection, 3) learn useful spatiotemporal representations, and 4) improve overall prediction. We conduct a demonstration of DeepLATTE using publicly available data for an important public health topic, air quality prediction, in a well-studied, complex physical environment - Los Angeles. The experiment demonstrates that the proposed approach provides accurate fine-spatial-scale air quality predictions and reveals the critical environmental factors affecting the results.


title: The Design and Performance of Charged Particle Detector onboard the GECAM Mission

abstract: The Gravitational Wave highly energetic Electromagnetic Counterpart All-sky Monitor (GECAM) is dedicated to detecting gravitational wave gamma-ray bursts. It is capable of all-sky monitoring over and discovering gamma-ray bursts and new radiation phenomena. GECAM consists of two microsatellites, each equipped with 8 charged particle detectors (CPDs) and 25 gamma-ray detectors (GRDs). The CPD is used to measure charged particles in the space environment, monitor energy and flow intensity changes, and identify between gamma-ray bursts and space charged particle events in conjunction with GRD. CPD uses plastic scintillator as the sensitive material for detection, silicon photomultiplier (SiPM) array as the optically readable device, and the inlaid Am-241 radioactive source as the onboard calibration means. In this paper, we will present the working principle, physical design, functional implementation and preliminary performance test results of the CPD.


title: Mathematical models of COVID-19 spread

abstract: The paper presents classification and analysis of the mathematical models of COVID-19 spread in different groups of populations such as the family, school, office (3-100 people), neighborhood (100-5000 people), city, region (0.5-15 million people), country, continent and the world. The classification covers the main types of models including time-series, differential, imitation ones, and their combinations. The time-series models are built from analysis of the time series derived using filtration, regression and network methods (Section 2). The differential models include those derived from systems of ordinary and stochastic differential equations as well as partial-derivative equations (Section 3). The imitation models include cellular automata and agent-based models (Section 4). The fourth group in the classification is combinations of nonlinear Markov chains and optimal control theory, derived within the framework of the mean-field game theory. Due to the novelty of the disease and the difficulties it causes, the parameters of most models are, as a rule, unknown, which necessitates one to solve the inverse problem, so the paper also analyses the main algorithms to solve the inverse problem such as stochastic optimization; nature-like algorithms (genetic; differential evolution; particle swarm, etc.); the understanding method; big-data analysis, and machine learning.


title: The DP Color Function of Clique-Gluings of Graphs

abstract: DP-coloring (also called correspondence coloring) is a generalization of list coloring that has been widely studied in recent years after its introduction by Dvo\v{r}\'{a}k and Postle in 2015. As the analogue of the chromatic polynomial of a graph $G$, $P(G,m)$, the DP color function of $G$, denoted by $P_{DP}(G,m)$, counts the minimum number of DP-colorings over all possible $m$-fold covers. Formulas for chromatic polynomials of clique-gluings of graphs, a fundamental graph operation, are well-known, but the effect of such gluings on the DP color function is not well understood. In this paper we study the DP color function of $K_p$-gluings of graphs. Recently, Becker et. al. asked whether $P_{DP}(G,m) \leq (\prod_{i=1}^n P_{DP}(G_i,m))/\left( \prod_{i=0}^{p-1} (m-i) \right)^{n-1}$ whenever $m \geq p$, where the expression on the right is the DP-coloring analogue of the corresponding chromatic polynomial formula for a $K_p$-gluing, $G$, of $G_1, \ldots, G_n$. Becker et. al. showed this inequality holds when $p=1$. In this paper we show this inequality holds for edge-gluings ($p=2$). On the other hand, we show it does not hold for triangle-gluings ($p=3$), which also answers a question of Dong and Yang (2021). Finally, we show a relaxed version, based on a class of $m$-fold covers that we conjecture would yield the fewest DP-colorings for a given graph, of the inequality holds when $p \geq 3$.


title: Connections between nonlocal operators: from vector calculus identities to a fractional Helmholtz decomposition

abstract: Nonlocal vector calculus, which is based on the nonlocal forms of gradient, divergence, and Laplace operators in multiple dimensions, has shown promising applications in fields such as hydrology, mechanics, and image processing. In this work, we study the analytical underpinnings of these operators. We rigorously treat compositions of nonlocal operators, prove nonlocal vector calculus identities, and connect weighted and unweighted variational frameworks. We combine these results to obtain a weighted fractional Helmholtz decomposition which is valid for sufficiently smooth vector fields. Our approach identifies the function spaces in which the stated identities and decompositions hold, providing a rigorous foundation to the nonlocal vector calculus identities that can serve as tools for nonlocal modeling in higher dimensions.


title: Potential Black Hole Seeding of the Spiral Galaxy NGC 4424 via an Infalling Star Cluster

abstract: Galaxies can grow through their mutual gravitational attraction and subsequent union. While orbiting a regular high-surface-brightness galaxy, the body of a low-mass galaxy can be stripped away. However, the stellar heart of the infalling galaxy, if represented by a tightly-bound nuclear star cluster, is more resilient. From archival Hubble Space Telescope images, we have discovered a red, tidally-stretched star cluster positioned ~5 arcseconds (~400 pc in projection) from, and pointing toward the center of, the post-merger spiral galaxy NGC 4424. The star cluster, which we refer to as `Nikhuli', has a near-infrared luminosity of (6.88+/-1.85)x10^6 L_{solar,F160W} and likely represents the nucleus of a captured/wedded galaxy. Moreover, from our Chandra X-ray Observatory image, Nikhuli is seen to contain a high-energy X-ray point source, with L_{0.5-8 keV} = 6.31^{+7.50}_{-3.77}x10^{38} erg/s (90% confidence). We argue that this is more likely to be an active massive black hole than an X-ray binary. Lacking an outward-pointing comet-like appearance, the stellar structure of Nikhuli favors infall rather than the ejection from a gravitational-wave recoil event. A minor merger with a low-mass early-type galaxy may have sown a massive black hole, aided an X-shaped pseudobulge, and be sewing a small bulge. The stellar mass and the velocity dispersion of NGC 4424 predict a central black hole of (0.6-1.0)x10^5 M_solar, similar to the expected intermediate-mass black hole in Nikhuli, and suggestive of a black hole supply mechanism for bulgeless late-type galaxies. We may potentially be witnessing black hole seeding by capture and sinking, with a nuclear star cluster the delivery vehicle.


title: Multivariate double truncated expectation and covariance risk measures for elliptical distributions

abstract: The main objective of this work is to calculate the multivariate double truncated expectation (MDTE) and covariance (MDTCov) for elliptical distributions. We also consider double truncated expectation (DTE) and variance (DTV) for univariate elliptical distributions. The exact expressions of MDTE and MDTCov are derived for some special cases of the family, such as normal, student-$t$, logistic, Laplace and Pearson type VII distributions. As numerical illustration, the DTE, DTV, MDTE and MDTCov for normal distribution are computed in details. Finally, we discuss MDTE and MDTCov of three industry segments' (Banks, Insurance, Financial and Credit Service) stock return in London stock exchange.


title: Open-Access Data and Toolbox for Tracking COVID-19 Impact on Power Systems

abstract: Intervention policies against COVID-19 have caused large-scale disruptions globally, and led to a series of pattern changes in the power system operation. Analyzing these pandemic-induced patterns is imperative to identify the potential risks and impacts of this extreme event. With this purpose, we developed an open-access data hub (COVID-EMDA+), an open-source toolbox (CoVEMDA), and a few evaluation methods to explore what the U.S. power systems are experiencing during COVID-19. These resources could be broadly used for research, public policy, and educational purposes. Technically, our data hub harmonizes a variety of raw data such as generation mix, demand profiles, electricity price, weather observations, mobility, confirmed cases and deaths. Typical methods are reformulated and standardized in our toolbox, including baseline estimation, regression analysis, and scientific visualization. Here the fluctuation index and probabilistic baseline are proposed for the first time to consider data fluctuation and estimation uncertainty. Based on these, we conduct three empirical studies on the U.S. power systems, and share new solutions and unexpected findings to address the issues of public concerns. This conveys a more complete picture of the pandemic's impacts, and also opens up several attractive topics for future work. Python, Matlab source codes, and user manuals are all publicly shared on a Github repository.


title: PMFL: Partial Meta-Federated Learning for heterogeneous tasks and its applications on real-world medical records

abstract: Federated machine learning is a versatile and flexible tool to utilize distributed data from different sources, especially when communication technology develops rapidly and an unprecedented amount of data could be collected on mobile devices nowadays. Federated learning method exploits not only the data but the computational power of all devices in the network to achieve more efficient model training. Nevertheless, while most traditional federated learning methods work well for homogeneous data and tasks, adapting the method to a different heterogeneous data and task distribution is challenging. This limitation has constrained the applications of federated learning in real-world contexts, especially in healthcare settings. Inspired by the fundamental idea of meta-learning, in this study we propose a new algorithm, which is an integration of federated learning and meta-learning, to tackle this issue. In addition, owing to the advantage of transfer learning for model generalization, we further improve our algorithm by introducing partial parameter sharing. We name this method partial meta-federated learning (PMFL). Finally, we apply the algorithms to two medical datasets. We show that our algorithm could obtain the fastest training speed and achieve the best performance when dealing with heterogeneous medical datasets.


title: Dynamic hardware system for cascade SVM classification of melanoma

abstract: Melanoma is the most dangerous form of skin cancer, which is responsible for the majority of skin cancer-related deaths. Early diagnosis of melanoma can significantly reduce mortality rates and treatment costs. Therefore, skin cancer specialists are using image-based diagnostic tools for detecting melanoma earlier. We aim to develop a handheld device featured with low cost and high performance to enhance early detection of melanoma at the primary healthcare. But, developing this device is very challenging due to the complicated computations required by the embedded diagnosis system. Thus, we aim to exploit the recent hardware technology in reconfigurable computing to achieve a high-performance embedded system at low cost. Support vector machine (SVM) is a common classifier that shows high accuracy for classifying melanoma within the diagnosis system and is considered as the most compute-intensive task in the system. In this paper, we propose a dynamic hardware system for implementing a cascade SVM classifier on FPGA for early melanoma detection. A multi-core architecture is proposed to implement a two-stage cascade classifier using two classifiers with accuracies of 98% and 73%. The hardware implementation results were optimized by using the dynamic partial reconfiguration technology, where very low resource utilization of 1% slices and power consumption of 1.5 W were achieved. Consequently, the implemented dynamic hardware system meets vital embedded system constraints of high performance and low cost, resource utilization, and power consumption, while achieving efficient classification with high accuracy.


title: Hot neutron stars and their equation of state

abstract: A set of microscopic, covariant density-functional, and non-relativistic Skyrme-type equations of state is employed to study the structure of purely nucleonic neutron stars at finite temperature. After examining the agreement with presently available astrophysical observational constraints, we find that the magnitude of thermal effects depends on the nucleon effective mass as well as on the stiffness of the cold equation of state. We evidence a fairly small but model-dependent effect of finite temperature on stellar stability that is correlated with the relative thermal pressure inside the star.


title: Attention-based Transformation from Latent Features to Point Clouds

abstract: In point cloud generation and completion, previous methods for transforming latent features to point clouds are generally based on fully connected layers (FC-based) or folding operations (Folding-based). However, point clouds generated by FC-based methods are usually troubled by outliers and rough surfaces. For folding-based methods, their data flow is large, convergence speed is slow, and they are also hard to handle the generation of non-smooth surfaces. In this work, we propose AXform, an attention-based method to transform latent features to point clouds. AXform first generates points in an interim space, using a fully connected layer. These interim points are then aggregated to generate the target point cloud. AXform takes both parameter sharing and data flow into account, which makes it has fewer outliers, fewer network parameters, and a faster convergence speed. The points generated by AXform do not have the strong 2-manifold constraint, which improves the generation of non-smooth surfaces. When AXform is expanded to multiple branches for local generations, the centripetal constraint makes it has properties of self-clustering and space consistency, which further enables unsupervised semantic segmentation. We also adopt this scheme and design AXformNet for point cloud completion. Considerable experiments on different datasets show that our methods achieve state-of-the-art results.


title: Efficient training of artificial neural network surrogates for a collisional-radiative model through adaptive parameter space sampling

abstract: Reliable plasma transport modeling for magnetic confinement fusion depends on accurately resolving the ion charge state distribution and radiative power losses of the plasma. These quantities can be obtained from solutions of a collisional-radiative (CR) model at each time step within a plasma transport simulation. However, even compact, approximate CR models can be computationally onerous to evaluate, and in-situ evaluations of these models within a coupled plasma transport code can lead to a rigid bottleneck. A way to bypass this bottleneck is to deploy artificial neural network surrogates for rapid evaluations of the necessary plasma quantities. However, one issue with training an accurate artificial neural network surrogate is the reliance on a sufficiently large and representative data set for both training and validation, which can be time-consuming to generate. In this study we further explore a data-driven active learning and training scheme to allow autonomous adaptive sampling of the problem parameter space that ensures a sufficiently large and meaningful set of training data assembled for the surrogate training. Our numerical experiments show that in order to produce a comparably accurate CR surrogate, the proposed approach requires a total number of data samples that is an order-of-magnitude smaller than a conventional approach.


title: Born Machines for Periodic and Open XY Quantum Spin Chains

abstract: Quantum phase transitions are ubiquitous in quantum many body systems. The quantum fluctuations that occur at very low temperatures are known to be responsible for driving the system across different phases as a function of an external control parameter. The XY Hamiltonian with a transverse field is a basic model that manifests two distinct quantum phase transitions, including spontaneous $Z_2$ symmetry breaking from an ordered to a disordered state. While programmable quantum devices have shown great success in investigating the various exotic quantum phases of matter, in parallel, the quest for harnessing machine learning tools in learning quantum phases of matter is ongoing. In this paper, we present a numerical study of the power of a quantum-inspired generative model known as the Born machine in learning quantum phases of matter. Data obtained from the system under open and periodic boundary conditions is considered. Our results indicate that a Born machine based on matrix product states can successfully capture the quantum state across various phases of the XY Hamiltonian and close to a critical point, despite the existence of long-range correlations. We further impose boundary conditions on the Born machine and show that matching the boundary condition of the Born machine and that of the training data improves performance when limited data is available and a small bond dimension is employed.


title: Novel Multiferroic Nanocomposite with High Pressure-Modulated Magnetoelectric Coupling

abstract: In this work, we have designed obtained a novel multiferroic nanocomposite using the high-pressure torsion (HPT) method. The crystal structure, phase composition, morphology, ferromagnetic (FM) and ferroelectric (FE) properties of the initial powders and "multiferroic-ferromagnetic" nanocomposites have been studied comprehensively. The initial powders and their composites show the perovskite and spinel crystalline phases for the FE and FM fractions, respectively. After HPT, the particle sizes of the initial powders are decreased significantly. It is shown that the novel nanocomposite consists of exchange-interacting FE and FM phases and demonstrates improved magnetic and electrical properties in low fields at room temperature. A colossal increase in residual polarization with an increase in external high-pressure is found in new composite. The obtained results make it possible to consider the novel nanocomposite as a new functional material for its use both in electronic devices for monitoring ultra-high-pressure and in integrated circuits of high-speed computing nanosystems with low switching energy. The HPT method is a promising method for obtaining new heterophase multiferroic nanosystems.


title: Multimodal Interactions Using Pretrained Unimodal Models for SIMMC 2.0

abstract: This paper presents our work on the Situated Interactive MultiModal Conversations 2.0 challenge held at Dialog State Tracking Challenge 10. SIMMC 2.0 includes 4 subtasks, and we introduce our multimodal approaches for the subtask \#1, \#2 and the generation of subtask \#4. SIMMC 2.0 dataset is a multimodal dataset containing image and text information, which is more challenging than the problem of only text-based conversations because it must be solved by understanding the relationship between image and text. Therefore, since there is a limit to solving only text models such as BERT or GPT2, we propose a multimodal model combining image and text. We first pretrain the multimodal model to understand the relationship between image and text, then finetune our model for each task. We achieve the 3rd best performance in subtask \#1, \#2 and a runner-up in the generation of subtask \#4. The source code is available at https://github.com/rungjoo/simmc2.0.


title: FaceFormer: Speech-Driven 3D Facial Animation with Transformers

abstract: Speech-driven 3D facial animation is challenging due to the complex geometry of human faces and the limited availability of 3D audio-visual data. Prior works typically focus on learning phoneme-level features of short audio windows with limited context, occasionally resulting in inaccurate lip movements. To tackle this limitation, we propose a Transformer-based autoregressive model, FaceFormer, which encodes the long-term audio context and autoregressively predicts a sequence of animated 3D face meshes. To cope with the data scarcity issue, we integrate the self-supervised pre-trained speech representations. Also, we devise two biased attention mechanisms well suited to this specific task, including the biased cross-modal multi-head (MH) attention and the biased causal MH self-attention with a periodic positional encoding strategy. The former effectively aligns the audio-motion modalities, whereas the latter offers abilities to generalize to longer audio sequences. Extensive experiments and a perceptual user study show that our approach outperforms the existing state-of-the-arts. The code will be made available.


title: Tensor Charges and their Impact on Physics Beyond the Standard Model

abstract: The nucleon tensor charge, $g_T$, is an important quantity in the search for beyond the Standard Model tensor interactions in neutron and nuclear $\beta$-decays as well as the contribution of the quark electric dipole moment (EDM) to the neutron EDM. We present results from the QCDSF/UKQCD/CSSM collaboration for the tensor charge, $g_T$, using lattice QCD methods and the Feynman-Hellmann theorem. We use a flavour symmetry breaking method to systematically approach the physical quark mass using ensembles that span three lattice spacings.


title: Segmenting Time Series via Self-Normalization

abstract: We propose a novel and unified framework for change-point estimation in multivariate time series. The proposed method is fully nonparametric, enjoys effortless tuning and is robust to temporal dependence. One salient and distinct feature of the proposed method is its versatility, where it allows change-point detection for a broad class of parameters (such as mean, variance, correlation and quantile) in a unified fashion. At the core of our method, we couple the self-normalization (SN) based tests with a novel nested local-window segmentation algorithm, which seems new in the growing literature of change-point analysis. Due to the presence of an inconsistent long-run variance estimator in the SN test, non-standard theoretical arguments are further developed to derive the consistency and convergence rate of the proposed SN-based change-point detection method. Extensive numerical experiments and relevant real data analysis are conducted to illustrate the effectiveness and broad applicability of our proposed method in comparison with state-of-the-art approaches in the literature.


title: Dispersive qubit readout with machine learning

abstract: Open quantum systems can undergo dissipative phase transitions, and their critical behavior can be exploited in sensing applications. For example, it can be used to enhance the fidelity of superconducting qubit readout measurements, a central problem toward the creation of reliable quantum hardware. A recently introduced measurement protocol, named ``critical parametric quantum sensing'', uses the parametric (two-photon driven) Kerr resonator's driven-dissipative phase transition to reach single-qubit detection fidelity of 99.9\% [arXiv:2107.04503]. In this work, we improve upon the previous protocol by using machine learning-based classification algorithms to \textit{efficiently and rapidly} extract information from this critical dynamics, which has so far been neglected to focus only on stationary properties. These classification algorithms are applied to the time series data of weak quantum measurements (homodyne detection) of a circuit-QED implementation of the Kerr resonator coupled to a superconducting qubit. This demonstrates how machine learning methods enable a faster and more reliable measurement protocol in critical open quantum systems.


title: Neural quantum states for supersymmetric quantum gauge theories

abstract: Supersymmetric quantum gauge theories are important mathematical tools in high energy physics. As an example, supersymmetric matrix models can be used as a holographic description of quantum black holes. The wave function of such supersymmetric gauge theories is not known and it is challenging to obtain with traditional techniques. We employ a neural quantum state ansatz for the wave function of a supersymmetric matrix model and use a variational quantum Monte Carlo approach to discover the ground state of the system. We discuss the difficulty of including bosonic particles and fermionic particles, as well as gauge degrees of freedom.


title: Controlling Valley-Polarisation in Graphene via Tailored Light Pulses

abstract: Analogous to charge and spin, electrons in solids endows an additional degree of freedom: the valley pseudospin. Two-dimensional hexagonal materials such as graphene exhibit two valleys, labelled as $\mathbf{K}$ and $\mathbf{K}^{\prime}$. These two valleys have the potential to realise logical operations in two-dimensional materials. Obtaining the desired control over valley polarisation between the two valleys is a prerequisite for the logical operations. Recently, it was shown that two counter-rotating circularly polarised laser pulses can induce a significant valley-polarisation in graphene. The main focus of the present work is to optimise the valley polarisation in monolayer graphene by controlling different laser parameters, such as wavelength, intensity ratio, frequency ratio and sub-cycle phase in two counter-rotating circularly polarised laser setup. Moreover, an alternate approach, based on single or few-cycle linearly polarised laser pulse, is also explored to induce significant valley polarisation in graphene. Our work could help experimentalists to choose a suitable method with optimised parameter space to obtain the desired control over valley polarisation in monolayer graphene.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el physics.optics


title: Uncertainty, Edge, and Reverse-Attention Guided Generative Adversarial Network for Automatic Building Detection in Remotely Sensed Images

abstract: Despite recent advances in deep-learning based semantic segmentation, automatic building detection from remotely sensed imagery is still a challenging problem owing to large variability in the appearance of buildings across the globe. The errors occur mostly around the boundaries of the building footprints, in shadow areas, and when detecting buildings whose exterior surfaces have reflectivity properties that are very similar to those of the surrounding regions. To overcome these problems, we propose a generative adversarial network based segmentation framework with uncertainty attention unit and refinement module embedded in the generator. The refinement module, composed of edge and reverse attention units, is designed to refine the predicted building map. The edge attention enhances the boundary features to estimate building boundaries with greater precision, and the reverse attention allows the network to explore the features missing in the previously estimated regions. The uncertainty attention unit assists the network in resolving uncertainties in classification. As a measure of the power of our approach, as of December 4, 2021, it ranks at the second place on DeepGlobe's public leaderboard despite the fact that main focus of our approach -- refinement of the building edges -- does not align exactly with the metrics used for leaderboard rankings. Our overall F1-score on DeepGlobe's challenging dataset is 0.745. We also report improvements on the previous-best results for the challenging INRIA Validation Dataset for which our network achieves an overall IoU of 81.28% and an overall accuracy of 97.03%. Along the same lines, for the official INRIA Test Dataset, our network scores 77.86% and 96.41% in overall IoU and accuracy.


title: Crystal Electic Field and Possible Coupling with Phonon in Kondo lattice CeCuGa3

abstract: We investigate the magnetic and crystal electric field (CEF) states of the Kondo lattice system CeCuGa3 by muon spin relaxation (muSR), neutron diffraction, and inelastic neutron scattering (INS) measurements. A noncentrosymmetric BaNiSn3-type tetragonal crystal structure (space group I4mm) is inferred from x-ray as well as from neutron powder diffraction. The low-temperature magnetic susceptibility and heat capacity data show an anomaly near 2.3 - 2.5~K associated with long range magnetic ordering, which is further confirmed by muSR and neutron diffraction data. The neutron powder diffraction collected at 1.7 K shows the presence of magnetic Bragg peaks indexed by an incommensurate magnetic propagation vector k = (0.148, 0.148, 0) and the magnetic structure is best described by a longitudinal spin density wave with ordered moments lying in ab-plane. An analysis of the INS data based on a CEF model reveals the presence of two magnetic excitations near 4.5 meV and 6.9 meV. The magnetic heat capacity data suggest an overall CEF splitting of 20.7 meV, however the excitation between 20 and 30 meV is very broad and weak in our INS data, but could provide an evidence of CEF level in this energy range in agreement with the magnetic entropy. Our analysis of INS data based on the CEF-phonon model indicates that the two excitations at 4.5 meV and 6.9 meV have their origin in CEF-phonon coupling (i.e. splitting of one CEF peak into two peaks, called vibron), with an overall splitting of 28.16 meV, similar to the case of CeCuAl3 and CeAuAl3.


title: An Integrative Analysis of the HD 219134 Planetary System and the Inner Solar System: Extending DYNAMITE with Enhanced Orbital Dynamical Stability Criteria

abstract: Planetary architectures remain unexplored for the vast majority of exoplanetary systems, even among the closest ones, with potentially hundreds of planets still ``hidden" from our knowledge. DYNAMITE is a powerful software package that can predict the presence and properties of these yet undiscovered planets. We have significantly expanded the integrative capabilities of DYNAMITE, which now allows for (i) planets of unknown inclinations alongside planets of known inclinations, (ii) population statistics and model distributions for the eccentricity of planetary orbits, and (iii) three different dynamical stability criteria. We demonstrate the new capabilities with a study of the HD 219134 exoplanet system consisting of four confirmed planets and two likely candidates, where five of the likely planets are Neptune-size or below with orbital periods less than 100 days. By integrating the known data for the HD 219134 planetary system with contextual and statistical exoplanet population information, we tested different system architecture hypotheses to determine their likely dynamical stability. Our results provide support for the planet candidates, and we predict at least two additional planets in this system. We also deploy DYNAMITE on analogs of the inner Solar System by excluding Venus or Earth from the input parameters to test DYNAMITE's predictive power. Our analysis finds the system remains stable while also recovering the excluded planets, demonstrating the increasing capability of DYNAMITE to accurately and precisely model the parameters of additional planets in multi-planet systems.


title: There exist infinitely many kinds of partial separability/entanglement

abstract: In tri-partite systems, there are three basic biseparability, $A$-$BC$, $B$-$CA$ and $C$-$AB$ biseparability according to bipartitions of local systems. We begin with three convex sets consisting of these basic biseparable states in the three qubit system, and consider arbitrary iterations of intersections and/or convex hulls of them to get convex cones. One natural way to classify tri-partite states is to consider those convex sets to which they belong or do not belong. This is especially useful to classify partial entanglement of mixed states. We show that the lattice generated by those three basic convex sets with respect to convex hull and intersection has infinitely many mutually distinct members, to see that there are infinitely many kinds of three qubit partial entanglement. To do this, we consider an increasing chain of convex sets in the lattice and exhibit three qubit Greenberger-Horne-Zeilinger diagonal states distinguishing those convex sets in the chain.


title: Visualization of reaction chemistry in W-KClO4-BaCrO4 delay mixtures via a Sestak-Berggren model based isoconversional method

abstract: The combustion delay mixture of tungsten (W), potassium perchlorate (KClO4), and barium chromate (BaCrO4), also known as the WKB mixture, has long been considered to be an integral part of military-grade ammunition. Despite its long history, however, their progressive reaction dynamics remains a question mark, especially due to the complex nature of their combustion reaction. As opposed to a one-step oxidation commonly observed in conventional combustions, the WKB mixture is associated with a multibody reaction between its solid-state components. To this end, the emergence of three combustion peaks, which we corresponded with disparate chemical reactions, was observed using thermogravimetric analysis on two separate WKB mixtures with differing mixture ratios. We applied the stepwise isoconversional method on each of the peaks to match the combustion chemistry it represents to the Sestak-Berggren model and computed the conceptual activation energy. Further plotting the logarithmic pre-exponential factor as a function of the reaction progress, we demonstrate a method of using the plot as an intuitive tool to understand the dynamics of individual reactions that compose multi-step chemical reactions. Our study provides a systematic approach in visualizing the reaction chemistry, thereby strengthening the analytical arsenal against reaction dynamics of combustion compounds in general.


title: Tradeoffs Between Contrastive and Supervised Learning: An Empirical Study

abstract: Contrastive learning has made considerable progress in computer vision, outperforming supervised pretraining on a range of downstream datasets. However, is contrastive learning the better choice in all situations? We demonstrate two cases where it is not. First, under sufficiently small pretraining budgets, supervised pretraining on ImageNet consistently outperforms a comparable contrastive model on eight diverse image classification datasets. This suggests that the common practice of comparing pretraining approaches at hundreds or thousands of epochs may not produce actionable insights for those with more limited compute budgets. Second, even with larger pretraining budgets we identify tasks where supervised learning prevails, perhaps because the object-centric bias of supervised pretraining makes the model more resilient to common corruptions and spurious foreground-background correlations. These results underscore the need to characterize tradeoffs of different pretraining objectives across a wider range of contexts and training regimes.