title: A combinatorial model for lane merging

abstract: A two lane road approaches a stoplight. The left lane merges into the right just past the intersection. Vehicles approach the intersection one at a time, with some drivers always choosing the right lane, while others always choose the shorter lane, giving preference to the right lane to break ties. An arrival sequence of vehicles can be represented as a binary string, where the zeros represent drivers always choosing the right lane, and the ones represent drivers choosing the shorter lane. From each arrival sequence we construct a merging path, which is a lattice path determined by the lane chosen by each car. We give closed formulas for the number of merging paths reaching the point $(n,m)$ with exactly $k$ zeros in the arrival sequence, and the expected length of the right lane for all arrival sequences with exactly $k$ zeros. Proofs involve an adaptation of Andre's Reflection Principle. Other interesting connections also emerge, including to: Ballot numbers, the expected maximum number of heads or tails appearing in a sequence of $n$ coin flips, the largest domino snake that can be made using pieces up to $[n:n]$, and the longest trail on the complete graph $K_n$ with loops.


title: No More Than 6ft Apart: Robust K-Means via Radius Upper Bounds

abstract: Centroid based clustering methods such as k-means, k-medoids and k-centers are heavily applied as a go-to tool in exploratory data analysis. In many cases, those methods are used to obtain representative centroids of the data manifold for visualization or summarization of a dataset. Real world datasets often contain inherent abnormalities, e.g., repeated samples and sampling bias, that manifest imbalanced clustering. We propose to remedy such a scenario by introducing a maximal radius constraint $r$ on the clusters formed by the centroids, i.e., samples from the same cluster should not be more than $2r$ apart in terms of $\ell_2$ distance. We achieve this constraint by solving a semi-definite program, followed by a linear assignment problem with quadratic constraints. Through qualitative results, we show that our proposed method is robust towards dataset imbalances and sampling artifacts. To the best of our knowledge, ours is the first constrained k-means clustering method with hard radius constraints. Codes at https://bit.ly/kmeans-constrained


title: HyperTransformer: A Textural and Spectral Feature Fusion Transformer for Pansharpening

abstract: Pansharpening aims to fuse a registered high-resolution panchromatic image (PAN) with a low-resolution hyperspectral image (LR-HSI) to generate an enhanced HSI with high spectral and spatial resolution. Existing pansharpening approaches neglect using an attention mechanism to transfer HR texture features from PAN to LR-HSI features, resulting in spatial and spectral distortions. In this paper, we present a novel attention mechanism for pansharpening called HyperTransformer, in which features of LR-HSI and PAN are formulated as queries and keys in a transformer, respectively. HyperTransformer consists of three main modules, namely two separate feature extractors for PAN and HSI, a multi-head feature soft attention module, and a spatial-spectral feature fusion module. Such a network improves both spatial and spectral quality measures of the pansharpened HSI by learning cross-feature space dependencies and long-range details of PAN and LR-HSI. Furthermore, HyperTransformer can be utilized across multiple spatial scales at the backbone for obtaining improved performance. Extensive experiments conducted on three widely used datasets demonstrate that HyperTransformer achieves significant improvement over the state-of-the-art methods on both spatial and spectral quality measures. Implementation code and pre-trained weights can be accessed at https://github.com/wgcban/HyperTransformer.


title: The SEDIGISM survey: Molecular cloud morphology. I. Classification and star formation

abstract: We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspection also serves as a test for the J plots algorithm, as this is the first time it has been used on molecular gas. Generally, it has been found that the structure of molecular clouds is highly filamentary and our observations indeed verify that most of our molecular clouds are elongated structures. Based on our visual classification of the 10663 SEDIGISM clouds, 15% are ring-like, 57% are elongated, 15% are concentrated and 10% are clumpy clouds. The remaining clouds do not belong to any of these morphology classes and are termed unclassified. We compare the SEDIGISM molecular clouds with structures identified through other surveys, i.e. ATLASGAL elongated structures and the bubbles from Milky Way Project (MWP). We find that many of the ATLASGAL and MWP structures are velocity coherent. ATLASGAL elongated structures overlap with ~ 21% of the SEDIGISM elongated structures (elongated and clumpy clouds) and MWP bubbles overlap with ~ 25% of the SEDIGISM ring-like clouds. We also analyse the star-formation associated with different cloud morphologies using two different techniques. The first technique examines star formation efficiency (SFE) and the dense gas fraction (DGF), based on SEDIGISM clouds and ATLASGAL clumps data. The second technique uses the high-mass star formation (HMSF) threshold for molecular clouds. The results indicate that clouds with ring-like and clumpy morphologies show a higher degree of star formation.


title: ARM 4-BIT PQ: SIMD-based Acceleration for Approximate Nearest Neighbor Search on ARM

abstract: We accelerate the 4-bit product quantization (PQ) on the ARM architecture. Notably, the drastic performance of the conventional 4-bit PQ strongly relies on x64-specific SIMD register, such as AVX2; hence, we cannot yet achieve such good performance on ARM. To fill this gap, we first bundle two 128-bit registers as one 256-bit component. We then apply shuffle operations for each using the ARM-specific NEON instruction. By making this simple but critical modification, we achieve a dramatic speedup for the 4-bit PQ on an ARM architecture. Experiments show that the proposed method consistently achieves a 10x improvement over the naive PQ with the same accuracy.


title: Non-linear predictive vector quantization of speech

abstract: In this paper we propose a Non-Linear Predictive Vector quantizer (PVQ) for speech coding, based on Multi-Layer Perceptrons. We also propose a method to evaluate if a quantizer is well designed, and if it exploits the correlation between consecutive outputs. Although the results of the Non-linear PVQ do not improve the results of the non-linear scalar predictor, we check that there is some room for the PVQ improvement.


title: Parallel Fourier Ptychography reconstruction

abstract: Fourier ptychography has attracted a wide range of focus for its ability of large space-bandwidth-produce, and quantative phase measurement. It is a typical computational imaging technique which refers to optimizing both the imaging hardware and reconstruction algorithms simultaneously. The data redundancy and inverse problem algorithms are the sources of FPM's excellent performance. But at the same time, this large amount of data processing and complex algorithms also greatly reduce the imaging speed. In this article, we propose a parallel Fourier ptychography reconstruction framework consisting of three levels of parallel computing parts and implemented it with both central processing unit (CPU) and compute unified device architecture (CUDA) platform. In the conventional FPM reconstruction framework, the sample image is divided into multiple sub-regions for separately processing because the illumination angles for different subregions are varied for the same LED and different subregions contain different defocus distances due to the non-planar distribution or non-ideal posture of biological sample. We first build a parallel computing sub-framework in spatial domain based on the above-mentioned characteristics. And then, by utilizing the sequential characteristics of different spectrum regions to update, a parallel computing sub-framework in the spectrum domain is carried out in our scheme. The feasibility of the proposed parallel FPM reconstruction framework is verified with different experimental results acquired with the system we built.


title: Optimization of Traffic Control in MMAP[c]/PH[c]/S Catastrophic Queueing Model with PH Retrial Times and Controllable Preemptive Repeat Priority Policy

abstract: The presented study elaborates a multi-server catastrophic retrial queueing model considering preemptive repeat priority policy with phase-type (PH) distributed retrial times. For the sake of comprehension, the scenario of model operation prior and later to the occurrence of the disaster is referred to as the normal scenario and as the catastrophic scenario, respectively. In the normal scenario, the incoming heterogeneous calls are categorized as handoff calls and new calls. Handoff calls are provided controllable preemptive priority over new calls. In the catastrophic scenario, when a disaster causes the shut down of the entire system and failure of all functioning channels, a set of backup channels is quickly deployed to restore services. Due to the emergency situation in the concerned area, the incoming heterogeneous calls are divided into three categories: handoff, new call, and emergency calls. Emergency calls are provided controllable preemptive priority over new/handoff calls due to the pressing need to save lives in such situations. The Markov chain's ergodicity criteria are established by demonstrating that it belongs to the class of asymptotically quasi-Toeplitz Markov chains (AQTMC). Further, a multi-objective optimization problem to obtain optimal number of backup channels has been formulated and dealt by employing non-dominated sorting genetic algorithm-II (NSGA-II) approach.


title: Double Exponential method for Riemann Zeta, Lerch and Dirichlet L-functions

abstract: We provide efficient methods to evaluate the Riemann zeta, the Lerch zeta and the Dirichlet $L$-functions. The method uses the Riemann-Siegel (RS) type formulas and a modified double exponential (MDE) quadrature method near the saddle point of appropriate integrands. We provide simplified derivations of the RS formulas, containing finite series sums and residual integrals, for the Lerch and the Dirichlet $L$-functions. The MDE method allows us to remove the contribution of the singularities near the contour of integration for the residual integrals. The method allows us to evaluate the residual integrals to any prescribed accuracy. The numerical cost of evaluating them is minimal compared to the main series sums. Thus even highly oscillatory integrands for these functions can be evaluated with the same complexity as the RS formula. In particular, the numerical complexity to evaluate the $\zeta(s)$ and the Lerch zeta function with $s=\sigma+i t$ scales as $\sqrt{t}$, and for Dirichlet $L$-function it scales as $\max(q, \sqrt{q t})$. The method allows the automatic setting of integral cutoffs and finite discretization size to achieve prescribed accuracy. Furthermore, it ensures that every halving of the discretization interval almost leads to doubling the accuracy of the results. Thus, it allows for a further check on the accuracy of the results obtained.


title: Cellular Segmentation and Composition in Routine Histology Images using Deep Learning

abstract: Identification and quantification of nuclei in colorectal cancer haematoxylin \& eosin (H\&E) stained histology images is crucial to prognosis and patient management. In computational pathology these tasks are referred to as nuclear segmentation, classification and composition and are used to extract meaningful interpretable cytological and architectural features for downstream analysis. The CoNIC challenge poses the task of automated nuclei segmentation, classification and composition into six different types of nuclei from the largest publicly known nuclei dataset - Lizard. In this regard, we have developed pipelines for the prediction of nuclei segmentation using HoVer-Net and ALBRT for cellular composition. On testing on the preliminary test set, HoVer-Net achieved a PQ of 0.58, a PQ+ of 0.58 and finally a mPQ+ of 0.35. For the prediction of cellular composition with ALBRT on the preliminary test set, we achieved an overall $R^2$ score of 0.53, consisting of 0.84 for lymphocytes, 0.70 for epithelial cells, 0.70 for plasma and .060 for eosinophils.


title: Self-Supervised Learning for Joint Pushing and Grasping Policies in Highly Cluttered Environments

abstract: Robots often face situations where grasping a goal object is desirable but not feasible due to other present objects preventing the grasp action. We present a deep Reinforcement Learning approach to learn grasping and pushing policies for manipulating a goal object in highly cluttered environments to address this problem. In particular, a dual Reinforcement Learning model approach is proposed, which presents high resilience in handling complicated scenes, reaching an average of 98% task completion using primitive objects in a simulation environment. To evaluate the performance of the proposed approach, we performed two extensive sets of experiments in packed objects and a pile of object scenarios with a total of 1000 test runs in simulation. Experimental results showed that the proposed method worked very well in both scenarios and outperformed the recent state-of-the-art approaches. Demo video, trained models, and source code for the results reproducibility purpose are publicly available. https://sites.google.com/view/pushandgrasp/home


title: The Pristine survey -- XVII. The C-19 stream is dynamically hot and more extended than previously thought

abstract: The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($\sigma_{\rm [Fe/H]}$ $<$ 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream's properties. One of these new member stars is located 30 $\deg$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $\sigma_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$\deg\pm0.08\deg$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.


title: The Pristine survey XVIII: C-19: Tidal debris of a dark matter-dominated globular cluster?

abstract: The recently discovered C-19 stellar stream is a collection of kinematically associated metal-poor stars in the halo of the Milky Way lacking an obvious progenitor. The stream spans an arc of ~15 degrees in the sky, and orbit-fitting suggests an apocentric distance of ~20 kpc and a pericentre of ~10 kpc. The narrow metallicity dispersion of stars with available spectra, together with light element abundance variations, suggests a globular cluster (GC) origin. The observed metallicity ([Fe/H] ~ -3.4), however, is much lower than that of any known GC. In addition, the width and velocity dispersion of the stream are similar to those expected from disrupting dwarf galaxies, and substantially larger than the tidal debris of GCs able to disrupt on C-19's orbit. We propose here an unconventional model where the C-19 progenitor is a dark matter-dominated stellar system with GC-like abundance patterns. We use N-body simulations to show that the tidal disruption of a ~100 pc King-model stellar component embedded in a ~20 km/s cuspy cold dark matter halo yields debris consistent with C-19's observed width and velocity dispersion. The stellar component of the progenitor is fully disrupted, and is spread over two distinct streams; one corresponding to C-19 and another possibly hiding behind the Galactic plane. If such companion stream were found, it would suggest that dark matter-dominated dwarfs may also develop GC-like enrichment patterns, a finding that would inform our theoretical understanding of the formation of multiple populations in GCs and dwarf galaxies alike.


title: Preferential Accretion and Circumbinary Disk Precession in Eccentric Binary Systems

abstract: We present a suite of high resolution hydrodynamic simulations of binaries immersed in circumbinary accretion disks (CBDs). For the first time, we investigate the preferential accretion rate as a function of both eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$ in a densely sampled parameter space, finding that when compared with circular binaries, mass ratios grow more efficiently in binaries on moderately eccentric orbits ($0.0 \lesssim e_{\rm b} \lesssim 0.4$), and high eccentricities ($e_{\rm b} \gtrsim 0.6 $) suppress mass ratio growth. We suggest that this non-monotonic preferential accretion behaviour may produce an observable shift in the mass ratio distributions of stellar binaries and massive black hole binaries. We further find that the response of a CBD can be divided into three regimes, depending on eccentricity and mass ratio: (i) CBDs around circular binaries always precess freely, whereas CBDs around eccentric binaries either (ii) undergo forced precession or (iii) remain locked at an angle with respect to the binary periapsis. Forced precession in eccentric binaries is associated with strong modulation of individual accretion rates on the precession timescale, a potentially observable signature in accreting binaries with short orbital periods. We provide CBD locking angles and precession rates as a function of $e_{\rm b}$ and $q_{\rm b}$ for our simulation suite.


title: Suspicious Siblings: The Distribution of Mass and Spin Across Component Black Holes in Isolated Binary Evolution

abstract: The LIGO and Virgo gravitational-wave detectors have uncovered binary black hole systems with definitively nonzero spins, as well as systems with significant spin residing in the more massive black hole of the pair. We investigate the ability of isolated binary evolution in forming such highly spinning, asymmetric-mass systems through both accretion onto the first-born black hole and tidal spin-up of the second-born black hole using a rapid population synthesis approach with detailed considerations of spin-up through tidal interactions. Even with the most optimistic assumptions regarding the efficiency at which an accreting star receives material from a donor, we find that it is difficult to form systems with significant mass asymmetry and moderate or high spins in the primary black hole component. Assuming efficient angular momentum transport within massive stars and Eddington-limited accretion onto black holes, we find that $< 1.5\%$ of systems in the underlying binary black hole population have a primary black hole spin greater than 0.2 and a mass asymmetry of greater than 2:1 in our most optimistic models, with most models finding that this criteria is only met in $\sim 0.01\%$ of systems. The production of systems with significant mass asymmetries and spin in the primary black hole component is thus an unlikely byproduct of isolated evolution unless highly super-Eddington accretion is invoked or angular momentum transport in massive stars is less efficient than typically assumed.


title: Purely Virtual Particles in Quantum Gravity, Inflationary Cosmology and Collider Physics

abstract: We review the concept of purely virtual particle and its uses in quantum gravity, primordial cosmology and collider physics. The fake particle, or "fakeon", which mediates interactions without appearing among the incoming and outgoing states, can be introduced by means of a new diagrammatics. The renormalization coincides with the one of the parent Euclidean diagrammatics, while unitarity follows from spectral optical identities, which can be derived by means of algebraic operations. The classical limit of a theory of physical particles and fakeons is described by an ordinary Lagrangian plus Hermitian, micro acausal and micro nonlocal self-interactions. Quantum gravity propagates the graviton, a massive scalar field (the inflaton) and a massive spin-2 fakeon, and leads to a constrained primordial cosmology, which predicts the tensor-to-scalar ratio $r$ in the window $0.4\lesssim 1000r\lesssim 3.5$. The interpretation of inflation as a cosmic RG flow allows us to calculate the perturbation spectra to high orders in the presence of the Weyl squared term. In models of new physics beyond the standard model, fakeons evade various phenomenological bounds, because they are less constrained than normal particles. The resummation of self-energies reveals that it is impossible to get too close to the fakeon peak. The related peak uncertainty, equal to the fakeon width divided by 2, is expected to be observable.


title: Massive black hole formation in dense stellar environments: Enhanced X-ray detection rates in high velocity dispersion nuclear star clusters

abstract: We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass BH into a massive BH via tidal captures. Furthermore, there is a velocity dispersion threshold ($40\;\rm{km\;s^{-1}}$) above which NSCs should inevitably form a massive BH. To provide an observational test of these theories, we measure X-ray emission from NSCs and compare to the measured velocity dispersion and tidal capture runaway timescale. We find that NSCs above the $40\;\rm{km\;s^{-1}}$ threshold are X-ray detected at roughly twice the rate of those below (after accounting for contamination from X-ray binaries). These results are consistent with a scenario in which dense, high-velocity NSCs can form massive BHs, providing a formation pathway that does not rely on conditions found only at high redshift.


title: SU(2)-Symmetric Spin-Boson Model: Quantum Criticality, Fixed-Point Annihilation, and Duality

abstract: The annihilation of two intermediate-coupling renormalization-group (RG) fixed points is of interest in diverse fields from statistical mechanics to high-energy physics, but has so far only been studied using perturbative techniques. Here we present high-accuracy quantum Monte Carlo results for the SU(2)-symmetric $S=1/2$ spin-boson (or Bose-Kondo) model. We study the model with a power-law bath spectrum $\propto \omega^s$ where, in addition to a critical phase predicted by perturbative RG, a stable strong-coupling phase is present. Using a detailed scaling analysis, we provide direct numerical evidence for the collision and annihilation of two RG fixed points at $s^\ast = 0.6540(2)$, causing the critical phase to disappear for $s<s^\ast$. In particular, we uncover a surprising duality between the two fixed points, corresponding to a reflection symmetry of the RG beta function, which we utilize to make analytical predictions at strong coupling which are in excellent agreement with numerics. Our work makes phenomena of fixed-point annihilation accessible to large-scale simulations, and we comment on the consequences for impurity moments in critical magnets.

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


title: Robustness versus sensitivity in non-Hermitian topological lattices probed by pseudospectra

abstract: Non-Hermitian topological systems simultaneously posses two antagonistic features: ultra sensitivity due to exceptional points and robustness of topological zero energy modes, and it is unclear which one prevails under different perturbations. We study that question by applying the pseudospectra theory on the prototypical non-Hermitian SSH (NHSSH) lattice. Topological modes are robust with respect to chiral perturbations and sensitive to parity-time (PT) symmetry preserving perturbations. In fact, the chiral symmetry exactly at the exceptional point leads to the suppression of its sensitivity, corresponding to a lower order exceptional point. However, counterintuitively, they are most sensitive with respect to unstructured perturbations, leaving the fingerprint of the pertinent higher order exceptional (HEP) point.


title: Constraining 3-3-1 Models at the LHC and Future Hadron Colliders

abstract: In this work, we derive lower mass bounds on the Z' gauge boson based on the dilepton data from LHC with 13 TeV of center-of-mass energy, and forecast the sensitivity of the High-Luminosity-LHC with $L=3000 fb^{-1}$, the High-Energy LHC with $\sqrt{s}=27$ TeV, and also at the Future Circular Collider with $\sqrt{s}=100$ TeV. We take into account the presence of exotic and invisible decays of the Z' gauge boson to find a more conservative and robust limit, different from previous studies. We investigate the impact of these new decays channels for several benchmark models in the scope of two different 3-3-1 models. We found that in the most constraining cases, LHC with $139fb^{-1}$ can impose $m_{Z^{\prime}}>4$ TeV. Moreover, we forecast HL-LHC, HE-LHC, and FCC bounds that yield $m_{Z^{\prime}}>5.8$ TeV, $m_{Z^{\prime}}>9.9$ TeV, and $m_{Z^{\prime}}> 27$ TeV, respectively. Lastly, put our findings into perspective with dark matter searches to show the region of parameter space where a dark matter candidate with the right relic density is possible.


title: Quantum algorithms for grid-based variational time evolution

abstract: The simulation of quantum dynamics calls for quantum algorithms working in first quantized grid encodings. Here, we propose a variational quantum algorithm for performing quantum dynamics in first quantization. In addition to the usual reduction in circuit depth conferred by variational approaches, this algorithm also enjoys several advantages compared to previously proposed ones. For instance, variational approaches suffer from the need for a large number of measurements. However, the grid encoding of first quantized Hamiltonians only requires measuring in position and momentum bases, irrespective of the system size. Their combination with variational approaches is therefore particularly attractive. Moreover, heuristic variational forms can be employed to overcome the limitation of the hard decomposition of Trotterized first quantized Hamiltonians into quantum gates. We apply this quantum algorithm to the dynamics of several systems in one and two dimensions. Our simulations exhibit the previously observed numerical instabilities of variational time propagation approaches. We show how they can be significantly attenuated through subspace diagonalization at a cost of an additional $\mathcal{O}(MN^2)$ 2-qubit gates where $M$ is the number of dimensions and $N^M$ is the total number of grid points.


title: Metallicity, ionization parameter, and pressure variations of HII regions in the TYPHOON spiral galaxies

abstract: We present a spatially-resolved HII region study of the gas-phase metallicity, ionization parameter, and ISM pressure maps of 6 local star-forming and face-on spiral galaxies from the TYPHOON program. Self-consistent metallicity, ionization parameter, and pressure maps are calculated simultaneously through an iterative process to provide useful measures of the local chemical abundance and its relation to localized ISM properties. We constrain the presence of azimuthal variations in metallicity by measuring the residual metallicity offset $\Delta$(O/H) after subtracting the linear fits to the radial metallicity profiles. We however find weak evidence of azimuthal variations in most of the galaxies, with small (mean 0.03 dex) scatter. The galaxies instead reveal that HII regions with enhanced and reduced abundances are found distributed throughout the disk. While the spiral pattern plays a role in organizing the ISM, it alone does not establish the relatively uniform azimuthal variations we observe. Differences in the metal abundances are more likely driven by the strong correlations with the local physical conditions. We find a strong and positive correlation between the ionization parameter and the local abundances as measured by the relative metallicity offset $\Delta$(O/H), indicating a tight relationship between local physical conditions and their localized enrichment of the ISM. Additionally, we demonstrate the impact of unresolved observations on the measured ISM properties by rebinning the datacubes to simulate low-resolution (1 kpc) observations, typical of large IFU surveys. We find that the ionization parameter and ISM pressure diagnostics are impacted by the loss of resolution such that their measured values are larger relative to the measured values on sub-HII region scales.


title: Constraints on cosmologies inside black holes

abstract: We study the construction of holographic models with closed FRW cosmologies on the worldvolume of a constant-tension brane inside a Schwarzschild-AdS black hole. In dimensions d>2, having a smooth Euclidean solution where the brane does not self-intersect limits the brane tension to T<T_*, preventing us from realising a separation of scales between the brane and bulk curvature scales. We show that adding interface branes to this model does not relax the condition on the brane tension.


title: Non-Hermitian Quantum Quenches in Holography

abstract: The notion of non-Hermitian PT symmetric quantum theory has recently been generalized to the gauge/gravity duality. We study the evolution of such non-Hermitian holographic field theories when the couplings are varied with time with particular emphasis on the question non-unitary time vs. unitary time evolution. We show that a non-unitary time evolution in the dual quantum theory corresponds to a violation of the Null Energy Condition (NEC) in the bulk of the asymptotically AdS spacetime. We find that upon varying the non-Hermitian coupling the horizon of a bulk AdS black hole shrinks. On the other hand varying the Hermitian coupling in the presence of a constant non-Hermitian coupling still violates the NEC but results in a growing horizon. We also show that by introducing a non-Hermitian gauge field the time evolution can be made unitary, e.g. the NEC in the bulk is obeyed, and an exactly equivalent purely Hermitian description can be given.


title: Rounding near-optimal quantum strategies for nonlocal games to strategies using maximally entangled states

abstract: We establish approximate rigidity results for boolean constraint system (BCS) nonlocal games. In particular, we show that near-perfect quantum strategies are approximate representations of the corresponding BCS algebra in the little Frobenius norm. Likewise, for the class of XOR nonlocal games, we show that near-optimal quantum strategies are approximate representations of the corresponding *-algebra associated with the game. In both cases, the norm of the approximate representations is independent of the quantum state employed in the strategy. We also show that approximate representations of the BCS (resp. XOR-algebra) are close to near-perfect (resp. near-optimal) quantum strategies employing a maximally entangled state for the corresponding game. As a corollary, any near-perfect BCS (resp. near-optimal XOR) quantum strategy is close to a near-perfect (resp. near-optimal) quantum strategy using a maximally entangled state. Lastly, we show that every synchronous algebra is *-isomorphic to a certain BCS algebra, allowing us to apply our results to the class of synchronous nonlocal games as well.


title: Fermi and eRosita bubbles as relics of the past activity of the Galactic black hole

abstract: The newly launched X-ray satellite, eRosita, has recently revealed two gigantic bubbles extending to ~80 degrees above and below the Galactic center. The morphology of these "eRosita bubbles" bears a remarkable resemblance to the Fermi bubbles previously discovered by the Fermi Gamma-ray Space Telescope and its counterpart, the microwave haze. The physical origin of these striking structures has been intensely debated; however, because of their symmetry about the Galactic center, they likely originate from some energetic outbursts from the Galactic center in the past. Here we propose a theoretical model in which the eRosita bubbles, Fermi bubbles, and the microwave haze could be simultaneously explained by a single event of jet activity from the central supermassive black hole a few million years ago. Using numerical simulations, we show that this model could successfully reproduce the morphology and multi-wavelength spectra of the observed bubbles and haze, which allows us to derive critical constraints on the energetics and timescales of the outburst. This study serves as an important step forward in our understanding of the past Galactic center activity of our Milky Way Galaxy, and may bring valuable insights into the broader picture of supermassive black hole-galaxy co-evolution in the context of galaxy formation.


title: Linear Run Time of Persistent Homology Computation with GPU Parallelization

abstract: Persistent homology is a crucial invariant that is used in many areas to understand data. The $O(N^4)$ run time is a hindrance to its use on most large datasets. We give a parallelization method to utilize multi-core machines and clusters. We implement the computation of the $0^{th}$ persistent homology with OpenMP parallelization and observe a 1.75 fold performance increase by using 2 threads on a dual core machine. We also benchmark the computation using larger numbers of threads and show that the thread computational overhead decreases performance. With GPU parallelization, we analytically and empirically decrease the run time scaling from $O(N^4)$ to $O(N^3)$ and even $O(N^2)$ where $N$ is the number of data points, for a large enough GPU. Next, we analytically show run time scaling $O(N)$ for an even larger GPU.


title: Curvature in the Balance: The Weyl Functional and Scalar Curvature of 4-Manifolds

abstract: The infimum of the Weyl functional is shown to be surprisingly small on many compact 4-manifolds that admit positive-scalar-curvature metrics. Results are also proved that systematically compare the scalar and self-dual Weyl curvatures of certain almost-Kaehler 4-manifolds.


title: The Renovated Thacher Observatory and First Science Results

abstract: Located on the campus of the Thacher School in Southern California, the Thacher Observatory has a legacy of astronomy research and education that dates back to the late 1950's. In 2016, the observatory was fully renovated with upgrades including a new 0.7-m telescope, a research grade camera, and a slit dome with full automation capabilities. The low-elevation site is bordered by the Los Padres National Forest and therefore affords dark to very dark skies allowing for accurate and precise photometric observations. We present a characterization of the site including sky brightness, weather, and seeing, and we demonstrate the on-sky performance of the facility. Our primary research programs are based around our multi-band photometric capabilities and include photometric monitoring of variable sources, a nearby supernova search and followup program, a quick response transient followup effort, and exoplanet and eclipsing binary light curves. Select results from these programs are included in this work which highlight the broad range of science available to an automated observatory with a moderately sized telescope.


title: Machine Learning for CUDA+MPI Design Rules

abstract: We present a new strategy for automatically exploring the design space of key CUDA+MPI programs and providing design rules that discriminate slow from fast implementations. In such programs, the order of operations (e.g., GPU kernels, MPI communication) and assignment of operations to resources (e.g., GPU streams) makes the space of possible designs enormous. Systems experts have the task of redesigning and reoptimizing these programs to effectively utilize each new platform. This work provides a prototype tool to reduce that burden. In our approach, a directed acyclic graph of CUDA and MPI operations defines the design space for the program. Monte-Carlo tree search discovers regions of the design space that have large impact on the program's performance. A sequence-to-vector transformation defines features for each explored implementation, and each implementation is assigned a class label according to its relative performance. A decision tree is trained on the features and labels to produce design rules for each class; these rules can be used by systems experts to guide their implementations. We demonstrate our strategy using a key kernel from scientific computing -- sparse-matrix vector multiplication -- on a platform with multiple MPI ranks and GPU streams.


title: Global pointwise estimates of positive solutions to sublinear equations

abstract: We give bilateral pointwise estimates for positive solutions $u$ to the sublinear integral equation \[ u = \mathbf{G}(\sigma u^q) + f \quad \textrm{in} \,\, \Omega,\] for $0 < q < 1$, where $\sigma\ge 0$ is a measurable function, or a Radon measure, $f \ge 0$, and $\mathbf{G}$ is the integral operator associated with a positive kernel $G$ on $\Omega\times\Omega$. Our main results, which include the existence criteria and uniqueness of solutions, hold for quasi-metric, or quasi-metrically modifiable kernels $G$. As a consequence, we obtain bilateral estimates, along with the existence and uniqueness, for positive solutions $u$, possibly unbounded, to sublinear elliptic equations involving the fractional Laplacian, \[ (-\Delta)^{\frac{\alpha}{2}} u = \sigma u^q + \mu \quad \textrm{in} \,\, \Omega, \qquad u=0 \, \, \textrm{in} \,\, \Omega^c, \] where $0<q<1$, and $\mu, \sigma \ge 0$ are measurable functions, or Radon measures, on a bounded uniform domain $\Omega \subset \mathbf{R}^n$ for $0 < \alpha \le 2$, or on the entire space $\mathbf{R}^n$, a ball or half-space, for $0 < \alpha <n$.


title: Gas inflows in the polar ring of NGC 4111: the birth of an AGN

abstract: We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ \AA velocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$\gamma$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring.


title: BoostMIS: Boosting Medical Image Semi-supervised Learning with Adaptive Pseudo Labeling and Informative Active Annotation

abstract: In this paper, we propose a novel semi-supervised learning (SSL) framework named BoostMIS that combines adaptive pseudo labeling and informative active annotation to unleash the potential of medical image SSL models: (1) BoostMIS can adaptively leverage the cluster assumption and consistency regularization of the unlabeled data according to the current learning status. This strategy can adaptively generate one-hot "hard" labels converted from task model predictions for better task model training. (2) For the unselected unlabeled images with low confidence, we introduce an Active learning (AL) algorithm to find the informative samples as the annotation candidates by exploiting virtual adversarial perturbation and model's density-aware entropy. These informative candidates are subsequently fed into the next training cycle for better SSL label propagation. Notably, the adaptive pseudo-labeling and informative active annotation form a learning closed-loop that are mutually collaborative to boost medical image SSL. To verify the effectiveness of the proposed method, we collected a metastatic epidural spinal cord compression (MESCC) dataset that aims to optimize MESCC diagnosis and classification for improved specialist referral and treatment. We conducted an extensive experimental study of BoostMIS on MESCC and another public dataset COVIDx. The experimental results verify our framework's effectiveness and generalisability for different medical image datasets with a significant improvement over various state-of-the-art methods.


title: Discrete self-similar and ergodic Markov chains

abstract: The first aim of this paper is to introduce a class of Markov chains on $\mathbb{Z}_+$ which are discrete self-similar in the sense that their semigroups satisfy an invariance property expressed in terms of a discrete random dilation operator. After showing that this latter property requires the chains to be upward skip-free, we first establish a gateway relation, a concept introduced in [26], between the semigroup of such chains and the one of spectrally negative self-similar Markov processes on $\mathbb{R}_+$. As a by-product, we prove that each of these Markov chains, after an appropriate scaling, converge in the Skorohod metric, to the associated self-similar Markov process. By a linear perturbation of the generator of these Markov chains, we obtain a class of ergodic Markov chains, which are non-reversible. By means of intertwining and interweaving relations, where the latter was recently introduced in [27], we derive several deep analytical properties of such ergodic chains including the description of the spectrum, the spectral expansion of their semigroups, the study of their convergence to equilibrium in the $\Phi$-entropy sense as well as their hypercontractivity property.


title: A Semi-Monolithic Detector providing intrinsic DOI-encoding and sub-200 ps CRT TOF-Capabilities for Clinical PET Applications

abstract: Current clinical PET systems utilize detectors where the scintillator typically contains single elements of 3 - 6 mm width and about 20 mm height. While providing good TOF performance, this design limits the spatial resolution and causes radial astigmatism as the DOI remains unknown. We propose an alternative, aiming to combine the advantages of current detectors with the DOI capabilities shown for monolithic concepts, based on semi-monolithic scintillators (slabs). An array of 8 monolithic LYSO slabs of dimensions 3.9 x 32 x 19 mm3 was read out by a 64-channel photosensor containing digital SiPMs. The position estimation in the detector's monolithic and DOI direction was based on a calibration with a fan beam collimator and the machine learning technique gradient tree boosting (GTB). We achieved a positioning performance in terms of mean absolute error (MAE) of 1.44 mm and 2.12 mm for DOI considering a wide energy window of 300 - 700 keV. The energy resolution was determined to be 11.3%. We established both an analytical and machine-learning-based timing calibration approach and applied them for a first-photon trigger. The analytical timing calibration corrects for electronic and optical time skews leading to 240 ps coincidence resolving time (CRT) for a pair of slab-detectors. The CRT was significantly improved by utilizing GTB to predict the time difference based on specific training data and applied on top of the analytical calibration. We achieved 209 ps for the wide energy window and 198 ps for a narrow selection around the photopeak (411 - 561 keV). To maintain the detector's sensitivity, no filters were applied to the data during processing. Overall, the semi-monolithic detector provides attractive performance characteristics. Especially, a good CRT can be achieved while introducing DOI capabilities to the detector, making the concept suitable for clinical PET scanners.


title: Beyond Time Complexity: Data Movement Complexity Analysis for Matrix Multiplication

abstract: Data movement is becoming the dominant contributor to the time and energy costs of computation across a wide range of application domains. However, time complexity is inadequate to analyze data movement. This work expands upon Data Movement Distance, a recently proposed framework for memory-aware algorithm analysis, by 1) demonstrating that its assumptions conform with microarchitectural trends, 2) applying it to six variants of matrix multiplication, and 3) showing it to be capable of asymptotically differentiating algorithms with the same time complexity but different memory behavior, as well as locality optimized vs. non-optimized versions of the same algorithm. In doing so, we attempt to bridge theory and practice by combining the operation count analysis used by asymptotic time complexity with per-operation data movement cost resulting from hierarchical memory structure. Additionally, this paper derives the first fully precise, fully analytical form of recursive matrix multiplication's miss ratio curve on LRU caching systems. Our results indicate that the Data Movement Distance framework is a powerful tool going forward for engineers and algorithm designers to understand the algorithmic implications of hierarchical memory.


title: Groups having all elements off a normal subgroup with prime power order

abstract: We consider a finite group $G$ with a normal subgroup $N$ so that all elements of $G \setminus N$ have prime power order. We prove that if there is a prime $p$ so that all the elements in $G \setminus N$ have $p$-power order, then either $G$ is a $p$-group or $G = PN$ where $P$ is a Sylow $p$-subgroup and $(G,P,P \cap N)$ is a Frobenius-Wielandt triple. We also prove that if all the elements of $G \setminus N$ have prime power orders and the orders are divisible by two primes $p$ and $q$, then $G$ is a $\{ p, q \}$-group and $G/N$ is either a Frobenius group or a $2$-Frobenius group. If all the elements of $G \setminus N$ have prime power orders and the orders are divisible by at least three primes, then all elements of $G$ have prime power order and $G/N$ is nonsolvable.


title: Application placement study of environment adaptive software

abstract: To use heterogeneous hardware, programmers must have sufficient technical skills to utilize OpenMP, CUDA, and OpenCL. On the basis of this, I have proposed environment-adaptive software that enables automatic conversion, configuration, and high performance operation of once written code, in accordance with the hardware. However, although it has been considered to convert the code according to the offload devices, there has been no study where to place the offloaded applications to satisfy users' requirements of price and response time. In this paper, as a new element of environment-adapted software, I examine a method to calculate appropriate locations using linear programming method. I confirm that applications can be arranged appropriately through simulations.


title: Linear and nonlinear edge dynamics of trapped fractional quantum Hall droplets beyond the chiral Luttinger liquid paradigm

abstract: We report numerical studies of the linear and nonlinear edge dynamics of a non-harmonically confined macroscopic fractional quantum Hall fluid. In the long-wavelength and weak excitation limit, observable consequences of the fractional transverse conductivity are recovered. The first non-universal corrections to the chiral Luttinger liquid theory are then characterized: for a weak excitation in the linear response regime, cubic corrections to the linear wave dispersion and a broadening of the dynamical structure factor of the edge excitations are identified; for stronger excitations, sizable nonlinear effects are found in the dynamics. The numerically observed features are quantitatively captured by a nonlinear chiral Luttinger liquid quantum Hamiltonian that reduces to a driven Korteweg-de Vries equation in the semiclassical limit. Experimental observability of our predictions is finally discussed.


title: Evolving symbolic density functionals

abstract: Systematic development of accurate density functionals has been a decades-long challenge for scientists. Despite the emerging application of machine learning (ML) in approximating functionals, the resulting ML functionals usually contain more than tens of thousands parameters, which makes a huge gap in the formulation with the conventional human-designed symbolic functionals. We propose a new framework, Symbolic Functional Evolutionary Search (SyFES), that automatically constructs accurate functionals in the symbolic form, which is more explainable to humans, cheaper to evaluate, and easier to integrate to existing density functional theory codes than other ML functionals. We first show that without prior knowledge, SyFES reconstructed a known functional from scratch. We then demonstrate that evolving from an existing functional $\omega$B97M-V, SyFES found a new functional, GAS22 (Google Accelerated Science 22), that performs better for the majority of molecular types in the test set of Main Group Chemistry Database (MGCDB84). Our framework opens a new direction in leveraging computing power for the systematic development of symbolic density functionals.


title: Proof the non-existence of causal classical electrodynamics of point charged particles

abstract: Up until now, a consistent causal theory of point charged particles (for example electrons) interacting with electromagnetic field is not known. The well-known problem is that the standard Lorentz force alone (in the case of point particles) does not lead to a theory in which momentum and energy are conserved. The need of radiation reaction force (or self-force) thus arises. The well known candidate for such force, the Lorentz-Abraham-Dirac reaction force, gives non-causal particle behavior, i.e., the particle starts to move before the arrival of external electromagnetic fields. Alternative, causal proposals provide non-physical behavior of the particle -- the particle moves with non-zero acceleration long after any external forces acted on it. Below, we question the existence of a causal theory. We show that for certain electromagnetic pulse of radiation and point particle being initially at rest, there {\it does not exists} a causal particle trajectory, such that the particle ends up moving with constant velocity (when no external electromagnetic fields act on it for very long time). This shows that the proper, causal electrodynamics of point particles does not exist.


title: Resonance $X(4630)$

abstract: We investigate the structure $X(4630)$ discovered by the LHCb collaboration in the process $B^{+}\rightarrow J/\psi \phi K^{+}$ as a resonance in the $ J/\psi \phi $ mass distribution. We explore this resonance as a diquark-antidiquark state $X=[cs][\overline{c}\overline{s}]$ with spin-parities $J^{\mathrm{PC}}=1^{-+}$. Its mass and current coupling are calculated using the QCD two-point sum rule method by taking into account vacuum condensates up to dimension $10$. We also study decays of this tetraquark to mesons $J/\psi\phi $, $\eta _{c}\eta ^{(\prime )}$, and $\chi _{c1}\eta ^{(\prime )}$, and compute partial widths of these channels. To this end, we employ the light-cone sum rule approach and technical methods of soft-meson approximation to extract strong coupling at relevant tetraquark-meson-meson vertices. Our predictions for the mass $m=(4632\pm 60)~\mathrm{MeV}$ and width $\Gamma=(159\pm 31)~\mathrm{MeV}$ of $X$ are in a very nice agreement with recent measurements of the LHCb collaboration. These results allow us to interpret the resonance $X(4630)$ as the tetraquark $X$ with spin-parities $J^{\mathrm{PC}}=1^{-+}$.


title: Ridges, Neural Networks, and the Radon Transform

abstract: A ridge is a function that is characterized by a one-dimensional profile (activation) and a multidimensional direction vector. Ridges appear in the theory of neural networks as functional descriptors of the effect of a neuron, with the direction vector being encoded in the linear weights. In this paper, we investigate properties of the Radon transform in relation to ridges and to the characterization of neural networks. We introduce a broad category of hyper-spherical Banach subspaces (including the relevant subspace of measures) over which the back-projection operator is invertible. We also give conditions under which the back-projection operator is extendable to the full parent space with its null space being identifiable as a Banach complement. Starting from first principles, we then characterize the sampling functionals that are in the range of the filtered Radon transform. Next, we extend the definition of ridges for any distributional profile and determine their (filtered) Radon transform in full generality. Finally, we apply our formalism to clarify and simplify some of the results and proofs on the optimality of ReLU networks that have appeared in the literature.


title: Reinventing High Performance Computing: Challenges and Opportunities

abstract: The world of computing is in rapid transition, now dominated by a world of smartphones and cloud services, with profound implications for the future of advanced scientific computing. Simply put, high-performance computing (HPC) is at an important inflection point. For the last 60 years, the world's fastest supercomputers were almost exclusively produced in the United States on behalf of scientific research in the national laboratories. Change is now in the wind. While costs now stretch the limits of U.S. government funding for advanced computing, Japan and China are now leaders in the bespoke HPC systems funded by government mandates. Meanwhile, the global semiconductor shortage and political battles surrounding fabrication facilities affect everyone. However, another, perhaps even deeper, fundamental change has occurred. The major cloud vendors have invested in global networks of massive scale systems that dwarf today's HPC systems. Driven by the computing demands of AI, these cloud systems are increasingly built using custom semiconductors, reducing the financial leverage of traditional computing vendors. These cloud systems are now breaking barriers in game playing and computer vision, reshaping how we think about the nature of scientific computation. Building the next generation of leading edge HPC systems will require rethinking many fundamentals and historical approaches by embracing end-to-end co-design; custom hardware configurations and packaging; large-scale prototyping, as was common thirty years ago; and collaborative partnerships with the dominant computing ecosystem companies, smartphone, and cloud computing vendors.


title: Femtosecond Diffraction and Dynamic High Pressure Science

abstract: Solid state material at high pressure is prevalent throughout the Universe, and an understanding of the structure of matter under such extreme conditions, gleaned from x-ray diffraction, has been pursued for the best part of a century. The highest pressures that can be reached to date (2 TPa) in combination with x-ray diffraction diagnosis have been achieved by dynamic compression via laser ablation . The past decade has witnessed remarkable advances in x-ray technologies, with novel x-ray Free-Electron-Lasers (FELs) affording the capacity to produce high quality single-shot diffraction data on time scales below 100 fsec. We provide a brief history of the field of dynamic compression, spanning from when the x-ray sources were almost always laser-plasma based, to the current state-of-the art diffraction capabilities provided by FELs. We give an overview of the physics of dynamic compression, diagnostic techniques, and the importance of understanding how the rate of compression influences the final temperatures reached. We provide illustrative examples of experiments performed on FEL facilities that are starting to give insight into how materials deform at ultra-high strain rates, their phase diagrams, and the types of states that can be reached. We emphasize that there often appear to be differences in the crystalline phases observed between use of static and dynamic compression techniques. We give our perspective on both the current state of this rapidly evolving field, and some glimpses of how we see it developing in the near-to-medium term.


title: Nodal Precession and Tidal Evolution of Two Hot-Jupiters: WASP-33 b and KELT-9 b

abstract: Hot Jupiters orbiting rapidly rotating stars on inclined orbits undergo tidally induced nodal precession measurable over several years of observations. The Hot Jupiters WASP-33 b and KELT-9 b are particularly interesting targets as they are among the hottest planets found to date, orbiting relatively massive stars. Here, we analyze archival and new data that span 11 and 5 years for WASP-33 b and KELT-9 b, respectively, in order to to model and improve upon their tidal precession parameters. Our work confirms the nodal precession for WASP-33 b and presents the first clear detection of the precession of KELT-9 b. We determine that WASP-33 and KELT-9 have gravitational quadrupole moments $(6.3^{+1.2}_{-0.8})\times10^{-5}$ and $(3.26^{+0.93}_{-0.80})\times10^{-4}$, respectively. We estimate the planets' precession periods to be $1460^{+170}_{-130}$ years and $890^{+200}_{-140}$ years, respectively, and that they will cease to transit their host stars around the years $2090^{+17}_{-10}$~CE and $2074^{+12}_{-10}$~CE, respectively. Additionally, we investigate both planets' tidal and orbital evolution, suggesting that a high-eccentricity tidal migration scenario is possible to produce both system architectures and that they will most likely not be engulfed by their hosts before the end of their main sequence lifetimes.


title: Homomorphically Encrypted Computation using Stochastic Encodings

abstract: Homomorphic encryption (HE) is a privacy-preserving technique that enables computation directly over ciphertext. Unfortunately, a key challenge for HE is that implementations can be impractically slow and have limits on computation that can be efficiently implemented. For instance, in Boolean constructions of HE like TFHE, arithmetic operations need to be decomposed into constituent elementary logic gates to implement so performance depends on logical circuit depth. For even heavily quantized fixed-point arithmetic operations, these HE circuit implementations can be slow. This paper explores the merit of using stochastic computing (SC) encodings to reduce the logical depth required for HE computation to enable more efficient implementations. Contrary to computation in the plaintext space where many efficient hardware implementations are available, HE provides support for only a limited number of primitive operators and their performance may not directly correlate to their plaintext performance. Our results show that by layering SC encodings on top of TFHE, we observe similar challenges and limitations that SC faces in the plaintext space. Additional breakthroughs would require more support from the HE libraries to make SC with HE a viable solution.


title: Uncertainty Propagation for General Stochastic Hybrid Systems on Compact Lie Groups

abstract: This paper deals with uncertainty propagation of general stochastic hybrid systems (GSHS) where the continuous state space is a compact Lie group. A computational framework is proposed to solve the Fokker-Planck (FP) equation that describes the time evolution of the probability density function for the state of GSHS. The FP equation is split into two parts: the partial differential operator corresponding to the continuous dynamics, and the integral operator arising from the discrete dynamics. These two parts are solved alternatively using the operator splitting technique. Specifically, the partial differential equation is solved by the spectral method where the density function is decomposed into a linear combination of a complete orthonormal function basis brought forth by the Peter-Weyl theorem, thereby resulting an ordinary differential equation. Next, the integral equation is solved by approximating the integral by a finite summation using a quadrature rule. The proposed method is then applied to a three-dimensional rigid body pendulum colliding with a wall, evolving on the product of the three-dimensional special orthogonal group and the Euclidean space. It is illustrated that the proposed method exhibits numerical results consistent with a Monte Carlo simulation, while explicitly generating the density function that carries the complete stochastic information of the hybrid state.


title: Structured Pruning is All You Need for Pruning CNNs at Initialization

abstract: Pruning is a popular technique for reducing the model size and computational cost of convolutional neural networks (CNNs). However, a slow retraining or fine-tuning procedure is often required to recover the accuracy loss caused by pruning. Recently, a new research direction on weight pruning, pruning-at-initialization (PAI), is proposed to directly prune CNNs before training so that fine-tuning or retraining can be avoided. While PAI has shown promising results in reducing the model size, existing approaches rely on fine-grained weight pruning which requires unstructured sparse matrix computation, making it difficult to achieve real speedup in practice unless the sparsity is very high. This work is the first to show that fine-grained weight pruning is in fact not necessary for PAI. Instead, the layerwise compression ratio is the main critical factor to determine the accuracy of a CNN model pruned at initialization. Based on this key observation, we propose PreCropping, a structured hardware-efficient model compression scheme. PreCropping directly compresses the model at the channel level following the layerwise compression ratio. Compared to weight pruning, the proposed scheme is regular and dense in both storage and computation without sacrificing accuracy. In addition, since PreCropping compresses CNNs at initialization, the computational and memory costs of CNNs are reduced for both training and inference on commodity hardware. We empirically demonstrate our approaches on several modern CNN architectures, including ResNet, ShuffleNet, and MobileNet for both CIFAR-10 and ImageNet.


title: AgileWatts: An Energy-Efficient CPU Core Idle-State Architecture for Latency-Sensitive Server Applications

abstract: User-facing applications running in modern datacenters exhibit irregular request patterns and are implemented using a multitude of services with tight latency requirements. These characteristics render ineffective existing energy conserving techniques when processors are idle due to the long transition time from a deep idle power state (C-state). While prior works propose management techniques to mitigate this inefficiency, we tackle it at its root with AgileWatts (AW): a new deep C-state architecture optimized for datacenter server processors targeting latency-sensitive applications. AW is based on three key ideas. First, AW eliminates the latency overhead of saving/restoring the core context (i.e., micro-architectural state) when powering-off/-on the core in a deep idle power state by i) implementing medium-grained power-gates, carefully distributed across the CPU core, and ii) retaining context in the power-ungated domain. Second, AW eliminates the flush latency overhead (several tens of microseconds) of the L1/L2 caches when entering a deep idle power state by keeping L1/L2 cache content power-ungated. A minimal control logic also remains power-ungated to serve cache coherence traffic (i.e., snoops) seamlessly. AW implements sleep-mode in caches to reduce caches leakage power consumption and lowers a core voltage to the minimum operational voltage level to minimize the leakage power of the power-ungated domain. Third, using a state-of-the-art power efficient all-digital phase-locked loop (ADPLL) clock generator, AW keeps the PLL active and locked during the idle state, further cutting precious microseconds of wake-up latency at a negligible power cost. Our evaluation with an accurate simulator calibrated against an Intel Skylake server shows that AW reduces the energy consumption of Memcached by up to 71% (35% on average) with up to 1% performance degradation.


title: Proof Methods in Random Matrix Theory

abstract: In this survey article, we give an introduction to two methods of proof in random matrix theory: The method of moments and the Stieltjes transform method. We thoroughly develop these methods and apply them to show both the semicircle law and the Marchenko-Pastur law for random matrices with independent entries. The material is presented in a pedagogical manner and is suitable for anyone who has followed a course in measure-theoretic probability theory.


title: A streamable large-scale clinical EEG dataset for Deep Learning

abstract: Deep Learning has revolutionized various fields, including Computer Vision, Natural Language Processing, as well as Biomedical research. Within the field of neuroscience, specifically in electrophysiological neuroimaging, researchers are starting to explore leveraging deep learning to make predictions on their data without extensive feature engineering. The availability of large-scale datasets is a crucial aspect of allowing the experimentation of Deep Learning models. We are publishing the first large-scale clinical EEG dataset that simplifies data access and management for Deep Learning. This dataset contains eyes-closed EEG data prepared from a collection of 1,574 juvenile participants from the Healthy Brain Network. We demonstrate a use case integrating this framework, and discuss why providing such neuroinformatics infrastructure to the community is critical for future scientific discoveries.


title: Optimal Clock Synchronization with Signatures

abstract: Cryptographic signatures can be used to increase the resilience of distributed systems against adversarial attacks, by increasing the number of faulty parties that can be tolerated. While this is well-studied for consensus, it has been underexplored in the context of fault-tolerant clock synchronization, even in fully connected systems. Here, the honest parties of an $n$-node system are required to compute output clocks of small skew (i.e., maximum phase offset) despite local clock rates varying between $1$ and $\vartheta>1$, end-to-end communication delays varying between $d-u$ and $d$, and the interference from malicious parties. So far, it is only known that clock pulses of skew $d$ can be generated with (trivially optimal) resilience of $\lceil n/2\rceil-1$ (PODC `19), improving over the tight bound of $\lceil n/3\rceil-1$ holding without signatures for \emph{any} skew bound (STOC `84, PODC `85). Since typically $d\gg u$ and $\vartheta-1\ll 1$, this is far from the lower bound of $u+(\vartheta-1)d$ that applies even in the fault-free case (IPL `01). We prove matching upper and lower bounds of $\Theta(u+(\vartheta-1)d)$ on the skew for the resilience range from $\lceil n/3\rceil$ to $\lceil n/2\rceil-1$. The algorithm showing the upper bound is, under the assumption that the adversary cannot forge signatures, deterministic. The lower bound holds even if clocks are initially perfectly synchronized, message delays between honest nodes are known, $\vartheta$ is arbitrarily close to one, and the synchronization algorithm is randomized. This has crucial implications for network designers that seek to leverage signatures for providing more robust time. In contrast to the setting without signatures, they must ensure that an attacker cannot easily bypass the lower bound on the delay on links with a faulty endpoint.


title: Building 3D Generative Models from Minimal Data

abstract: We propose a method for constructing generative models of 3D objects from a single 3D mesh and improving them through unsupervised low-shot learning from 2D images. Our method produces a 3D morphable model that represents shape and albedo in terms of Gaussian processes. Whereas previous approaches have typically built 3D morphable models from multiple high-quality 3D scans through principal component analysis, we build 3D morphable models from a single scan or template. As we demonstrate in the face domain, these models can be used to infer 3D reconstructions from 2D data (inverse graphics) or 3D data (registration). Specifically, we show that our approach can be used to perform face recognition using only a single 3D template (one scan total, not one per person). We extend our model to a preliminary unsupervised learning framework that enables the learning of the distribution of 3D faces using one 3D template and a small number of 2D images. This approach could also provide a model for the origins of face perception in human infants, who appear to start with an innate face template and subsequently develop a flexible system for perceiving the 3D structure of any novel face from experience with only 2D images of a relatively small number of familiar faces.


title: Physics-informed data based neural networks for two-dimensional turbulence

abstract: Turbulence remains a problem that is yet to be fully understood, with experimental and numerical studies aiming to fully characterise the statistical properties of turbulent flows. Such studies require huge amount of resources to capture, simulate, store and analyse the data. In this work, we present physics-informed neural network (PINN) based methods to predict flow quantities and features of two-dimensional turbulence with the help of sparse data in a rectangular domain with periodic boundaries. While the PINN model can reproduce all the statistics at large scales, the small scale properties are not captured properly. We introduce a new PINN model that can effectively capture the energy distribution at small scales performing better than the standard PINN based approach. It relies on the training of the low and high wavenumber behaviour separately leading to a better estimate of the full turbulent flow. With 0.1 % training data, we observe that the new PINN model captures the turbulent field at inertial scales leading to a general agreement of the kinetic energy spectra upto eight to nine decades as compared with the solutions from direct numerical simulation (DNS). We further apply these techniques to successfully capture the statistical behaviour of large scale modes in the turbulent flow. We believe such methods to have significant applications in enhancing the retrieval of existing turbulent data sets at even shorter time intervals.


title: Improved Wavelets for Image Compression from Unitary Circuits

abstract: We benchmark the efficacy of several novel orthogonal, symmetric, dilation-3 wavelets, derived from a unitary circuit based construction, towards image compression. The performance of these wavelets is compared across several photo databases against the CDF-9/7 wavelets in terms of the minimum number of non-zero wavelet coefficients needed to obtain a specified image quality, as measured by the multi-scale structural similarity index (MS-SSIM). The new wavelets are found to consistently offer better compression efficiency than the CDF-9/7 wavelets across a broad range of image resolutions and quality requirements, averaging 7-8% improved compression efficiency on high-resolution photo images when high-quality (MS-SSIM = 0.99) is required.


title: UVCGAN: UNet Vision Transformer cycle-consistent GAN for unpaired image-to-image translation

abstract: Unpaired image-to-image translation has broad applications in art, design, and scientific simulations. One early breakthrough was CycleGAN that emphasizes one-to-one mappings between two unpaired image domains via generative-adversarial networks (GAN) coupled with the cycle-consistency constraint, while more recent works promote one-to-many mapping to boost diversity of the translated images. Motivated by scientific simulation and one-to-one needs, this work revisits the classic CycleGAN framework and boosts its performance to outperform more contemporary models without relaxing the cycle-consistency constraint. To achieve this, we equip the generator with a Vision Transformer (ViT) and employ necessary training and regularization techniques. Compared to previous best-performing models, our model performs better and retains a strong correlation between the original and translated image. An accompanying ablation study shows that both the gradient penalty and self-supervised pre-training are crucial to the improvement. To promote reproducibility and open science, the source code, hyperparameter configurations, and pre-trained model are available at https://github.com/LS4GAN/uvcgan.


title: Thermal drag effect in quantum Hall circuits

abstract: We study the thermal drag between two mesoscopic quantum Hall (QH) circuits. Each circuit consists of Ohmic contact perfectly coupled to quantum Hall edge states. The drag is caused by strong capacitive coupling between Ohmic contacts. The non-equilibirum conditions and the electron-electron interaction are taken into account by using the non-equilibrium bosonization technique. The thermal drag current in the passive circuit, the noise power of the corresponding heat current, and the Fano factor are calculated and analyzed.


title: Scale-invariance at the core of quantum black holes

abstract: We study spherically-symmetric solutions to a modified Einstein-Hilbert action with Renormalization Group scale-dependent couplings, inspired by Weinberg's Asymptotic Safety scenario for Quantum Gravity. The Renormalization Group scale is identified with the Tolman temperature for an isolated gravitational system in thermal equilibrium with Hawking radiation. As a result, the point of infinite local temperature is shifted from the classical black-hole horizon to the origin and coincides with a timelike curvature singularity. Close to the origin, the spacetime is determined by the scale-dependence of the cosmological constant in the vicinity of the Reuter fixed point: the free components of the metric can be derived analytically and are characterized by a radial power law with exponent $\alpha = \sqrt{3}-1$. Away from the fixed point, solutions for different masses are studied numerically and smoothly interpolate between the Schwarzschild exterior and the scale-invariant interior. Whereas the exterior of objects with astrophysical mass is described well by vacuum General Relativity, deviations become significant at a Planck distance away from the classical horizon and could lead to observational signatures. We further highlight potential caveats in this intriguing result with regard to our choice of scale-identification and identify future avenues to better understand quantum black holes in relation to the key feature of scale-invariance.


title: Improving sandwich variance estimation for marginal Cox analysis of cluster randomized trials

abstract: Cluster randomized trials (CRTs) frequently recruit a small number of clusters, therefore necessitating the application of small-sample corrections for valid inference. A recent systematic review indicated that CRTs reporting right-censored, time-to-event outcomes are not uncommon, and that the marginal Cox proportional hazards model is one of the common approaches used for primary analysis. While small-sample corrections have been studied under marginal models with continuous, binary and count outcomes, no prior research has been devoted to the development and evaluation of bias-corrected sandwich variance estimators when clustered time-to-event outcomes are analyzed by the marginal Cox model. To improve current practice, we propose 9 bias-corrected sandwich variance estimators for the analysis of CRTs using the marginal Cox model, and report on a simulation study to evaluate their small-sample properties. Our results indicate that the optimal choice of bias-corrected sandwich variance estimator for CRTs with survival outcomes can depend on the variability of cluster sizes, and can also slightly differ whether it is evaluated according to relative bias or type I error rate. Finally, we illustrate the new variance estimators in a real-world CRT where the conclusion about intervention effectiveness differs depending on the use of small-sample bias corrections. The proposed sandwich variance estimators are implemented in an R package CoxBcv.


title: Triplet-Pair Spin Signatures from Macroscopically Aligned Heteroacenes in an Oriented Single Crystal

abstract: The photo-driven process of singlet fission generates coupled triplet pairs (TT) with fundamentally intriguing and potentially useful properties. The quintet 5TT0 sublevel is particularly interesting for quantum information because it is highly entangled, addressable with microwave pulses, and could be detected using optical techniques. Previous theoretical work on a model Hamiltonian and nonadiabatic transition theory, called the JDE model, has determined that this sublevel can be selectively populated if certain conditions are met. Among the most challenging, the molecules within the dimer undergoing singlet fission must have their principal magnetic axes parallel to one another and to an applied Zeeman field B0. Here, we present time-resolved paramagnetic resonance (TR-EPR) spectroscopy of a single crystal sample of a novel tetracenethiophene compound featuring arrays of dimers aligned in this manner, mounted so that the orientation of the field relative to the molecular axes could be controlled. The observed spin sublevel populations in the paired TT and unpaired (T+T) triplets are consistent with predictions from the JDE model, including preferential 5TT0 formation at z||B0, with one caveat - two 5TT spin sublevels have little to no population. This may be due to crossings between the 5TT and 3TT manifolds in the field range investigated by TR-EPR, consistent with the inter-triplet exchange energy determined by monitoring photoluminescence at varying magnetic fields.


title: On Beltrami equations with the inverse conditions and hydrodynamic normalization

abstract: We consider problems concerning the existence of solutions of the Beltrami equations and their convergence in the entire complex plane. We are mainly interested in the case when these solutions satisfy the so-called hydrodynamic normalization condition in the neighborhood of infinity. Under conditions related to dilations of inverse mappings, we have established the existence of such solutions in the class of continuous Sobolev mappings. We have also obtained results on the locally uniform limit of a sequence of such solutions.


title: Heterogeneous Multi-Resource Allocation with Subset Demand Requests

abstract: We consider the problem of allocating multiple heterogeneous resources geographically and over time to meet demands that require some subset of the available resource types simultaneously at a specified time, location, and duration. The objective is to maximize the total reward accrued from meeting (a subset of) demands. We model this problem as an integer program, show that it is NP-hard, and analyze the complexity of various special cases. We introduce approximation algorithms and an extension to our problem that considers travel costs. Finally, we test the performance of the integer programming model in an extensive computational study.


title: The efficiency of Quantum Mechanical Carnot Engine using the Woods Saxon model

abstract: The quantum engine cycle serves as an analogous representation of classical heat engines for microscopic systems and the quantum regime of thermal devices is composed of a single element. In this work, the Quantum-Mechanical properties of a non-linear quantum oscillator described by the Woods-Saxon [WS] model are examined. The Quantum-Mechanical analogue of the Carnot cycle was constructed using changes in both the width L of the well and the quantum state of the potential well. The efficiency of the quantum engine, consisting of adiabatic and isothermal processes based on the Woods-Saxon [WS] potential is derived. The result is shown to be analogous to that of the classical engine and found to agree, within an appropriate limit, with existing results obtained from other potential models. This implies that the [WS] potential can be used as an alternative model in quantum engines.


title: SDSS IV MaNGA: Visual Morphological and Statistical Characterization of the DR15 sample

abstract: We present a detailed visual morphological classification for the 4614 MaNGA galaxies in SDSS Data Release 15, using image mosaics generated from a combination of r-band (SDSS and deeper DESI Legacy Surveys) images and their digital post-processing. We distinguish 13 Hubble types and identify the presence of bars and bright tidal debris. After correcting the MaNGA sample for volume completeness, we calculate the morphological fractions, the bi-variate distribution of type and stellar mass M*-where we recognise a morphological transition "valley" around S0a-Sa types- and the variations of the g-i colour and luminosity-weighted age over this distribution. We identified bars in 46.8% of galaxies, present in all Hubble types later than S0. This fraction amounts to a factor ~2 larger when compared with other works for samples in common. We detected 14% of galaxies with tidal features, with the fraction changing with M* and morphology. For 355 galaxies, the classification was uncertain; they are visually faint, mostly of low/intermediate masses, low concentrations, and disky in nature. Our morphological classification agrees well with other works for samples in common, though some particular differences emerge, showing that our image procedures allow us to identify a wealth of added value information as compared to SDSS-based previous estimates. Based on our classification, we also propose an alternative criteria for the E-S0 separation, in the structural semi-major to semi-minor axis versus bulge to total light ratio (b/a-B/T) and concentration versus semi-major to semi-minor axis (C-b/a) space.


title: Torus bundles over lens spaces

abstract: Let $p$ be an odd prime and let $\rho:\mathbb{Z}/p\rightarrow\operatorname{GL}_n(\mathbb{Z})$ be an action of $\mathbb{Z}/p$ on a lattice and let $\Gamma:=\mathbb{Z}^n\rtimes_{\rho}\mathbb{Z}/p$ be the corresponding semidirect product. The torus bundle $M:=T^n_{\rho}\times_{\mathbb{Z}/p}S^{\ell}$ over the lens space $S^{\ell}/\mathbb{Z}/p$ has fundamental group $\Gamma$. When $\mathbb{Z}/p$ fixes only the origin of $\mathbb{Z}^n$, Davis and L\"uck \cite{DavisLuckTorusBundles} compute the $L$-groups $L^{\langle j\rangle}_m(\mathbb{Z}[\Gamma])$ and the structure set $\mathcal{S}^{geo,s}(M)$. In this paper, we extend these computations to all actions of $\mathbb{Z}/p$ on $\mathbb{Z}^n$. In particular, we compute $L^{\langle j\rangle}_m(\mathbb{Z}[\Gamma])$ and $\mathcal{S}^{geo,s}(M)$ in a case where $\underline{E}\Gamma$ has a non-discrete singular set.


title: Exchange interactions and spin dynamics in the layered honeycomb ferromagnet CrI$_3$

abstract: We derive the microscopic spin Hamiltonian for rhombohedral CrI$_3$ using extensive first-principles density functional theory (DFT) calculations which incorporate spin-orbit coupling and Hubbard U. Our calculations indicate a dominant nearest-neighbor ferromagnetic Heisenberg exchange with weaker further-neighbor Heisenberg terms. In addition, we find a Dzyaloshinskii-Moriya interaction which primarily drives a topological gap in the spin-wave spectrum at the Dirac point, and uncover a non-negligible antiferromagnetic Kitaev coupling between the S=3/2 Cr moments. The out-of-plane magnetic moment is stabilized by weak symmetric bond-dependent terms and a local single-ion anisotropy. Using linear spin wave theory, we find that our exchange parameters are in reasonably good agreement with inelastic neutron scattering (INS) experiments. Employing classical Monte Carlo simulations, we study the magnetic phase transition temperature $T_c$ and its evolution with an applied in-plane magnetic field. We further demonstrate how future high-resolution INS experiments on the magnon dispersion of single crystals in an in-plane magnetic field may be used to quantitatively extract the strength of the antiferromagnetic Kitaev exchange coupling.


title: NEMAR: An open access data, tools, and compute resource operating on NeuroElectroMagnetic data

abstract: To take advantage of recent and ongoing advances in large-scale computational methods, and to preserve the scientific data created by publicly funded research projects, data archives must be created as well as standards for specifying, identifying, and annotating deposited data. The OpenNeuro.org archive, begun as a repository for magnetic resonance imaging (MRI) data, is such an archive. We present a gateway to OpenNeuro for human electrophysiology data (BIDS-formatted EEG and MEG, as well as intracranial data). The NEMAR gateway allows users to visualize electrophysiological data, including time-domain and frequency-domain dynamics time locked to sets of experimental events recorded using BIDS- and HED-formatted data annotation. In addition, NEMAR allows users to process archived EEG data on the XSEDE high-performance resources at SDSC in conjunction with the Neuroscience Gateway (nsgportal.org), a freely available and easy to use portal to leverage high-performance computing resources for neuroscience research.


title: Coverage Properties of Empirical Bayes Intervals

abstract: This note is an invited discussion of the article "Confidence Intervals for Nonparametric Empirical Bayes Analysis" by Ignatiadis and Wager. In this discussion, I review some goals of empirical Bayes data analysis and the contribution of Ignatiadis and Wager. Differences between across-group inference and group-specific inference are discussed. Standard empirical Bayes interval procedures focus on controlling the across-group average coverage rate. However, if group-specific inferences are of primary interest, confidence intervals with group-specific coverage control may be preferable.


title: Bayesian Optimization Meets Hybrid Zero Dynamics: Safe Parameter Learning for Bipedal Locomotion Control

abstract: In this paper, we propose a multi-domain control parameter learning framework that combines Bayesian Optimization (BO) and Hybrid Zero Dynamics (HZD) for locomotion control of bipedal robots. We leverage BO to learn the control parameters used in the HZD-based controller. The learning process is firstly deployed in simulation to optimize different control parameters for a large repertoire of gaits. Next, to tackle the discrepancy between the simulation and the real world, the learning process is applied on the physical robot to learn for corrections to the control parameters learned in simulation while also respecting a safety constraint for gait stability. This method empowers an efficient sim-to-real transition with a small number of samples in the real world, and does not require a valid controller to initialize the training in simulation. Our proposed learning framework is experimentally deployed and validated on a bipedal robot Cassie to perform versatile locomotion skills with improved performance on smoothness of walking gaits and reduction of steady-state tracking errors.


title: Improving the Energy Efficiency and Robustness of tinyML Computer Vision using Log-Gradient Input Images

abstract: This paper studies the merits of applying log-gradient input images to convolutional neural networks (CNNs) for tinyML computer vision (CV). We show that log gradients enable: (i) aggressive 1.5-bit quantization of first-layer inputs, (ii) potential CNN resource reductions, and (iii) inherent robustness to illumination changes (1.7% accuracy loss across 1/32...8 brightness variation vs. up to 10% for JPEG). We establish these results using the PASCAL RAW image data set and through a combination of experiments using neural architecture search and a fixed three-layer network. The latter reveal that training on log-gradient images leads to higher filter similarity, making the CNN more prunable. The combined benefits of aggressive first-layer quantization, CNN resource reductions, and operation without tight exposure control and image signal processing (ISP) are helpful for pushing tinyML CV toward its ultimate efficiency limits.


title: Inverse design of plasma metamaterial devices with realistic elements

abstract: In an expansion of a previous study [1], we apply inverse design methods to produce two-dimensional plasma metamaterial devices with realistic plasma elements which incorporate quartz envelopes, collisionality (loss), non-uniform density profiles, and resistance to experimental error/perturbation. Backpropagated finite difference frequency domain simulations are used to design waveguides and demultiplexers operating under the transverse magnetic polarization. Optimal devices with realistic elements are compared to previous devices with idealized elements, and several parameter initialization schemes for the optimization algorithm are explored. Demultiplexing and waveguiding are demonstrated for microwave-regime devices composed of plasma elements with reasonable space-averaged plasma frequencies ~10 GHz and a collision frequency ~1 GHz, allowing for future in-situ training and experimental realization of these designs.


title: Show Me What and Tell Me How: Video Synthesis via Multimodal Conditioning

abstract: Most methods for conditional video synthesis use a single modality as the condition. This comes with major limitations. For example, it is problematic for a model conditioned on an image to generate a specific motion trajectory desired by the user since there is no means to provide motion information. Conversely, language information can describe the desired motion, while not precisely defining the content of the video. This work presents a multimodal video generation framework that benefits from text and images provided jointly or separately. We leverage the recent progress in quantized representations for videos and apply a bidirectional transformer with multiple modalities as inputs to predict a discrete video representation. To improve video quality and consistency, we propose a new video token trained with self-learning and an improved mask-prediction algorithm for sampling video tokens. We introduce text augmentation to improve the robustness of the textual representation and diversity of generated videos. Our framework can incorporate various visual modalities, such as segmentation masks, drawings, and partially occluded images. It can generate much longer sequences than the one used for training. In addition, our model can extract visual information as suggested by the text prompt, e.g., "an object in image one is moving northeast", and generate corresponding videos. We run evaluations on three public datasets and a newly collected dataset labeled with facial attributes, achieving state-of-the-art generation results on all four.


title: Style-ERD: Responsive and Coherent Online Motion Style Transfer

abstract: Motion style transfer is a common method for enriching character animation. Motion style transfer algorithms are often designed for offline settings where motions are processed in segments. However, for online animation applications, such as realtime avatar animation from motion capture, motions need to be processed as a stream with minimal latency. In this work, we realize a flexible, high-quality motion style transfer method for this setting. We propose a novel style transfer model, Style-ERD, to stylize motions in an online manner with an Encoder-Recurrent-Decoder structure, along with a novel discriminator that combines feature attention and temporal attention. Our method stylizes motions into multiple target styles with a unified model. Although our method targets online settings, it outperforms previous offline methods in motion realism and style expressiveness and provides significant gains in runtime efficiency


title: Electron acceleration from the interaction of three crossed parallel Alfv\'en waves

abstract: We study the non-linear interaction of three parallel Alfv\'en wave packets (AWP) in an initially uniform plasma using 2.5D particle-in-a-cell (PIC) numerical simulations. We aim to help to explain the observation of suprathermal electrons by multiple Alfv\'en waves collision in regions where these waves are trapped like IAR (Ionospheric Alfv\'en Resonator), Earth radiation belts or coronal magnetic loops. In the context of APAWI process that have been described by Mottez (2012, 2015), the interaction of two parallel Alfv\'en waves (AW) generates longitudinal density modulations and parallel electric field at the APAWI crossing region that can accelerate particles effectively in the direction of the background magnetic field. Our simulations show that when a third parallel AWP of different initial position arrives at the APAWI crossing region, it gives rise to a strong parallel electron beams ($V \sim 5--7VTe$) at longitudinal cavity density gradients. We suggest that velocity drift from outgoing AW generates interface waves in the transverse direction, which allows propagating waves to develop parallel electric fields by phase mixing process when the transverse scale of the wavy density gradient (oblique gradient) is in the range of the electron inertial length.

categories: astro-ph.SR astro-ph.EP physics.plasm-ph physics.space-ph


title: Machine Learning Simulates Agent-Based Model Towards Policy

abstract: Public Policies are not intrinsically positive or negative. Rather, policies provide varying levels of effects across different recipients. Methodologically, computational modeling enables the application of multiple influences on empirical data, thus allowing for heterogeneous response to policies. We use a random forest machine learning algorithm to emulate an agent-based model (ABM) and evaluate competing policies across 46 Metropolitan Regions (MRs) in Brazil. In doing so, we use input parameters and output indicators of 11,076 actual simulation runs and one million emulated runs. As a result, we obtain the optimal (and non-optimal) performance of each region over the policies. Optimum is defined as a combination of GDP production and the Gini coefficient inequality indicator for the full ensemble of Metropolitan Regions. Results suggest that MRs already have embedded structures that favor optimal or non-optimal results, but they also illustrate which policy is more beneficial to each place. In addition to providing MR-specific policies' results, the use of machine learning to simulate an ABM reduces the computational burden, whereas allowing for a much larger variation among model parameters. The coherence of results within the context of larger uncertainty--vis-\`a-vis those of the original ABM--reinforces robustness of the model. At the same time the exercise indicates which parameters should policymakers intervene on, in order to work towards precise policy optimal instruments.


title: Confirming Brennan's conjecture numerically on a counterexample to Thurston's $K=2$ conjecture

abstract: It was shown by Bishop that if Thurston's $K = 2$ conjecture holds for some planar domain, then Brennan's conjecture holds for the Riemann map of that domain as well. In this paper we show numerically that the original counterexample to Thurston's $K=2$ conjecture given by Epstein, Marden and Markovi\'{c} is not a counterexample to Brennan's conjecture.


title: Harmonic projections in negative curvature

abstract: In this paper we construct harmonic maps that are at a bounded distance from nearest-point retractions to convex sets, in negatively curved manifolds. Specifically, given a quasidisk $Q$ in hyperbolic space, we construct a harmonic map to the hyperbolic plane that corresponds to the nearest-point retraction to the convex hull of $Q$. If $M$ is a pinched Hadamard manifold so that its isometry group acts with cobounded orbits, and if $S$ is a set in the boundary at infinity of $M$, with the property that all elements of its orbit under the isometry group of $M$ have dimension less than $\frac{n-1}{2}$, we show that the nearest-point retraction to the convex hull of $S$ is a bounded distance away from some harmonic map.


title: The Status and Future of Color Transparency and Nuclear Filtering

abstract: 40 years after its introduction, the phenomenon of color transparency remains a domain of controversial interpretations of experimental data. We review present evidence for or against its manifestation in various exclusive hard scattering reactions. The nuclear transparency experiments reveal whether short distance processes dominate a scattering amplitude at some given kinematical point. We plead for a new round of nuclear transparency measurements in a variety of experimental set-ups, including near-forward exclusive reactions related to generalized parton distribution (GPD) physics and near-backward exclusive reactions related to transition distribution amplitudes (TDA) physics.