title: Hurricane-blackout-heatwave Compound Hazard Risk and Resilience in a Changing Climate

abstract: Hurricanes have caused power outages and blackouts, affecting millions of customers and inducing severe social and economic impacts. The impacts of hurricane-caused blackouts may worsen due to increased heat extremes and possibly increased hurricanes under climate change. We apply hurricane and heatwave projections with power outage and recovery process analysis to investigate how the emerging hurricane-blackout-heatwave compound hazard may vary in a changing climate, for Harris County in Texas (including major part of Houston City) as an example. We find that, under the high-emissions scenario RCP8.5, the expected percent of customers experiencing at least one longer-than-5-day hurricane-induced power outage in a 20-year period would increase significantly from 14% at the end of the 20th century to 44% at the end of the 21st century in Harris County. The expected percent of customers who may experience at least one longer-than-5-day heatwave without power (to provide air conditioning) would increase alarmingly, from 0.8% to 15.5%. These increases of risk may be largely avoided if the climate is well controlled under the stringent mitigation scenario RCP2.6. We also reveal that a moderate enhancement of critical sectors of the distribution network can significantly improve the resilience of the entire power grid and mitigate the risk of the future compound hazard. Together these findings suggest that, in addition to climate mitigation, climate adaptation actions are urgently needed to improve the resilience of coastal power systems.


title: Isolated singularities in the heat equation behaving like fractional Brownian motions

abstract: We consider solutions of the linear heat equation in $\mathbb{R}^N$ with isolated singularities. It is assumed that the position of a singular point depends on time and is H\"older continuous with the exponent $\alpha \in (0,1)$. We show that any isolated singularity is removable if it is weaker than a certain order depending on $\alpha$. We also show the optimality of the removability condition by showing the existence of a solution with a nonremovable singularity. These results are applied to the case where the singular point behaves like a fractional Brownian motion with the Hurst exponent $H \in (0,1/2] $. It turns out that $H=1/N$ is critical.


title: Adversarial Disentanglement of Speaker Representation for Attribute-Driven Privacy Preservation

abstract: In speech technologies, speaker's voice representation is used in many applications such as speech recognition, voice conversion, speech synthesis and, obviously, user authentication. Modern vocal representations of the speaker are based on neural embeddings. In addition to the targeted information, these representations usually contain sensitive information about the speaker, like the age, sex, physical state, education level or ethnicity. In order to allow the user to choose which information to protect, we introduce in this paper the concept of attribute-driven privacy preservation in speaker voice representation. It allows a person to hide one or more personal aspects to a potential malicious interceptor and to the application provider. As a first solution to this concept, we propose to use an adversarial autoencoding method that disentangles in the voice representation a given speaker attribute thus allowing its concealment. We focus here on the sex attribute for an Automatic Speaker Verification (ASV) task. Experiments carried out using the VoxCeleb datasets have shown that the proposed method enables the concealment of this attribute while preserving ASV ability.


title: Ratio of counts vs ratio of rates in Poisson processes

abstract: The often debated issue of `ratios of small numbers of events' is approached from a probabilistic perspective, making a clear distinction between the predictive problem (forecasting numbers of events we might count under well stated assumptions, and therefore of their ratios) and inferential problem (learning about the relevant parameters of the related probability distribution, in the light of the observed number of events). The quantities of interests and their relations are visualized in a graphical model (`Bayesian network'), very useful to understand how to approach the problem following the rules of probability theory. In this paper, written with didactic intent, we discuss in detail the basic ideas, however giving some hints of how real life complications, like (uncertain) efficiencies and possible background and systematics, can be included in the analysis, as well as the possibility that the ratio of rates might depend on some physical quantity. The simple models considered in this paper allow to obtain, under reasonable assumptions, closed expressions for the rates and their ratios. Monte Carlo methods are also used, both to cross check the exact results and to evaluate by sampling the ratios of counts in the cases in which large number approximation does not hold. In particular it is shown how to make approximate inferences using a Markov Chain Monte Carlo using JAGS/rjags. Some examples of R and JAGS code are provided.


title: Understanding How Dimension Reduction Tools Work: An Empirical Approach to Deciphering t-SNE, UMAP, TriMAP, and PaCMAP for Data Visualization

abstract: Dimension reduction (DR) techniques such as t-SNE, UMAP, and TriMAP have demonstrated impressive visualization performance on many real world datasets. One tension that has always faced these methods is the trade-off between preservation of global structure and preservation of local structure: these methods can either handle one or the other, but not both. In this work, our main goal is to understand what aspects of DR methods are important for preserving both local and global structure: it is difficult to design a better method without a true understanding of the choices we make in our algorithms and their empirical impact on the lower-dimensional embeddings they produce. Towards the goal of local structure preservation, we provide several useful design principles for DR loss functions based on our new understanding of the mechanisms behind successful DR methods. Towards the goal of global structure preservation, our analysis illuminates that the choice of which components to preserve is important. We leverage these insights to design a new algorithm for DR, called Pairwise Controlled Manifold Approximation Projection (PaCMAP), which preserves both local and global structure. Our work provides several unexpected insights into what design choices both to make and avoid when constructing DR algorithms.


title: Codimensional Incremental Potential Contact

abstract: We extend the incremental potential contact (IPC) model for contacting elastodynamics to resolve systems composed of codimensional DOFs in arbitrary combination. This enables a unified, interpenetration-free, robust, and stable simulation framework that couples codimension-0,1,2, and 3 geometries seamlessly with frictional contact. Extending IPC to thin structures poses new challenges in computing strain, modeling thickness and determining collisions. To address these challenges we propose three corresponding contributions. First, we introduce a C2 constitutive barrier model that directly enforces strain limiting as an energy potential while preserving rest state. This provides energetically-consistent strain limiting models (both isotropic and anisotropic) for cloth that enable strict satisfaction of strain-limit inequalities with direct coupling to both elastodynamics and contact via minimization of the incremental potential. Second, to capture the geometric thickness of codimensional domains we extend the IPC model to directly enforce distance offsets. Our treatment imposes a strict guarantee that mid-surfaces (resp. mid-lines) of shells (resp. rods) will not move closer than applied thickness values. This enables us to account for thickness in the contact behavior of codimensional structures and so robustly capture challenging contacting geometries; a number of which, to our knowledge, have not been simulated before. Third, codimensional models, especially with modeled thickness, mandate strict accuracy requirements that pose a severe challenge to all existing continuous collision detection (CCD) methods. To address these limitations we develop a new, efficient, simple-to-implement additive CCD (ACCD) method that applies conservative advancement to iteratively refine a lower bound for deforming primitives, converging to time of impact.


title: Entangled particle-states localized on curved space-time

abstract: In this work, based on a recently introduced localization scheme for scalar fields, we argue that the geometry of the space-time, where the particle states of a scalar field are localized, is intimately related to the quantum entanglement of these states. More specifically, we show that on curved space-time can only be localized entangled states, while separable states are located on flat space-time. Our result goes in parallel with recent theoretical developments in the context of AdS/CFT correspondence which uncovered connections between gravity and quantum entanglement.


title: Exposing the trion's fine structure by controlling the carrier concentration in hBN-encapsulated MoS$_2$

abstract: Atomically thin materials, like semiconducting transition metal dichalcogenides, are highly sensitive to the environment. This opens up an opportunity to externally control their properties by changing their surroundings. In this work, high-quality van der Waals heterostructures assembled from hBN-encapsulated monolayer MoS$_2$ are studied with the aid of photoluminescence, photoluminescence excitation, and reflectance contrast experiments. We demonstrate that carrier concentration in MoS$_2$ monolayers, arising from charge transfer from impurities in the substrate, can be significantly tuned within one order of magnitude by the modification of the bottom hBN flake thickness. The studied structures, characterized by spectral lines approaching the narrow homogeneously broadened limit enabled observations of subtle optical and spin-valley properties of excitonic complexes. Our results allowed us to resolve three optically-active negatively charged excitons in MoS$_2$ monolayers, which are assigned to the intravalley singlet, intervalley singlet, and intervalley triplet states.


title: Planet cartography with neural learned regularization

abstract: Finding potential life harboring exo-Earths is one of the aims of exoplanetary science. Detecting signatures of life in exoplanets will likely first be accomplished by determining the bulk composition of the planetary atmosphere via reflected/transmitted spectroscopy. However, a complete understanding of the habitability conditions will surely require mapping the presence of liquid water, continents and/or clouds. Spin-orbit tomography is a technique that allows us to obtain maps of the surface of exoplanets around other stars using the light scattered by the planetary surface. We leverage the potential of deep learning and propose a mapping technique for exo-Earths in which the regularization is learned from mock surfaces. The solution of the inverse mapping problem is posed as a deep neural network that can be trained end-to-end with suitable training data. We propose in this work to use methods based on the procedural generation of planets, inspired by what we found on Earth. We also consider mapping the recovery of surfaces and the presence of persistent cloud in cloudy planets. We show that the a reliable mapping can be carried out with our approach, producing very compact continents, even when using single passband observations. More importantly, if exoplanets are partially cloudy like the Earth is, we show that one can potentially map the distribution of persistent clouds that always occur on the same position on the surface (associated to orography and sea surface temperatures) together with non-persistent clouds that move across the surface. This will become the first test one can perform on an exoplanet for the detection of an active climate system. For small rocky planets in the habitable zone of their stars, this weather system will be driven by water, and the detection can be considered as a strong proxy for truly habitable conditions.


title: Combining Reinforcement Learning with Lin-Kernighan-Helsgaun Algorithm for the Traveling Salesman Problem

abstract: We address the Traveling Salesman Problem (TSP), a famous NP-hard combinatorial optimization problem. And we propose a variable strategy reinforced approach, denoted as VSR-LKH, which combines three reinforcement learning methods (Q-learning, Sarsa and Monte Carlo) with the well-known TSP algorithm, called Lin-Kernighan-Helsgaun (LKH). VSR-LKH replaces the inflexible traversal operation in LKH, and lets the program learn to make choice at each search step by reinforcement learning. Experimental results on 111 TSP benchmarks from the TSPLIB with up to 85,900 cities demonstrate the excellent performance of the proposed method.


title: Multi-Objective Interpolation Training for Robustness to Label Noise

abstract: Deep neural networks trained with standard cross-entropy loss memorize noisy labels, which degrades their performance. Most research to mitigate this memorization proposes new robust classification loss functions. Conversely, we propose a Multi-Objective Interpolation Training (MOIT) approach that jointly exploits contrastive learning and classification to mutually help each other and boost performance against label noise. We show that standard supervised contrastive learning degrades in the presence of label noise and propose an interpolation training strategy to mitigate this behavior. We further propose a novel label noise detection method that exploits the robust feature representations learned via contrastive learning to estimate per-sample soft-labels whose disagreements with the original labels accurately identify noisy samples. This detection allows treating noisy samples as unlabeled and training a classifier in a semi-supervised manner to prevent noise memorization and improve representation learning. We further propose MOIT+, a refinement of MOIT by fine-tuning on detected clean samples. Hyperparameter and ablation studies verify the key components of our method. Experiments on synthetic and real-world noise benchmarks demonstrate that MOIT/MOIT+ achieves state-of-the-art results. Code is available at https://git.io/JI40X.


title: Scalar Fully-heavy Tetraquark States $QQ^\prime \bar{Q} \bar{Q^\prime}$ in QCD Sum Rules

abstract: Very recently, the LHCb Collaboration observed distinct structures with the $cc\bar{c}\bar{c}$ in the $J/\Psi$-pair mass spectrum. In this work, we construct four scalar ($J^{PC} = 0^{++}$) $[8_c]_{Q\bar{Q^\prime}}\otimes [8_c]_{Q^\prime \bar{Q}}$ type currents to investigate the fully-heavy tetraquark state $Q Q^\prime \bar{Q} \bar{Q^\prime}$ in the framework of QCD sum rules, where $Q=c, b$ and $Q^\prime = c, b$. Our results suggest that the broad structure around 6.2-6.8 GeV can be interpreted as the $0^{++}$ octet-octet tetraquark states with masses $6.44\pm 0.11$ GeV and $6.52\pm 0.10$ GeV, and the narrow structure around $6.9$ GeV can be interpreted as the $0^{++}$ octet-octet tetraquark states with masses $6.87\pm 0.11$ GeV and $6.96\pm 0.11$ GeV, respectivley. Extending to the b-quark sector,the masses of their fully-bottom partners are found to be around 18.38-18.59 GeV. Additionally, we also analyze the spectra of the $[c\bar{c}][b\bar{b}]$ and $[c\bar{b}] [b \bar{c}]$ tetraquark states, which lie in the range of 12.51-12.74 GeV and 12.49-12.81 GeV, respectively.


title: Bridging Bayesian, frequentist and fiducial (BFF) inferences using confidence distribution

abstract: Bayesian, frequentist and fiducial (BFF) inferences are much more congruous than they have been perceived historically in the scientific community (cf., Reid and Cox 2015; Kass 2011; Efron 1998). Most practitioners are probably more familiar with the two dominant statistical inferential paradigms, Bayesian inference and frequentist inference. The third, lesser known fiducial inference paradigm was pioneered by R.A. Fisher in an attempt to define an inversion procedure for inference as an alternative to Bayes' theorem. Although each paradigm has its own strengths and limitations subject to their different philosophical underpinnings, this article intends to bridge these different inferential methodologies through the lenses of confidence distribution theory and Monte-Carlo simulation procedures. This article attempts to understand how these three distinct paradigms, Bayesian, frequentist, and fiducial inference, can be unified and compared on a foundational level, thereby increasing the range of possible techniques available to both statistical theorists and practitioners across all fields.


title: Surface Diffusion Control Enables Tailored Aspect Ratio Nanostructures in Area-Selective Atomic Layer Deposition

abstract: Area-selective atomic layer deposition is a key technology for modern microelectronics as it eliminates alignment errors inherent to conventional approaches by enabling material deposition only in specific areas. Typically, the selectivity originates from surface modifications of the substrate that allow or block precursor adsorption. The control of the deposition process currently remains a major challenge as the selectivity of the no-growth areas is lost quickly. Here, we show that surface modifications of the substrate strongly manipulate the surface diffusion. The selective deposition of TiO$_2$ on poly (methyl methacrylate) and SiO$_2$ yields localized nanostructures with tailored aspect ratios. Controlling the surface diffusion allows to tune such nanostructures as it boosts the growth rate at the interface of the growth and no-growth areas. Kinetic Monte-Carlo calculations reveal that species move from high to low diffusion areas. Further, we identify the catalytic activity of TiCl$_4$ during the formation of carboxylic acid on poly (methyl methacrylate) as the reaction mechanism responsible for the loss of selectivity, and show that process optimization leads to higher selectivity. Our work enables the precise control of area-selective atomic layer deposition on the nanoscale, and offers new strategies in area-selective deposition processes by exploiting surface diffusion effects.


title: Arbitrarily Accurate Analytical Approximations for the Error Function

abstract: In this paper a spline based integral approximation is utilized to propose a sequence of approximations to the error function that converge at a significantly faster manner than the default Taylor series. The approximations can be improved by utilizing the approximation erf(x) approximately equal to one for x>>1. Two generalizations are possible, the first is based on demarcating the integration interval into m equally spaced sub-intervals. The second, it based on utilizing a larger fixed sub-interval, with a known integral, and a smaller sub-interval whose integral is to be approximated. Both generalizations lead to significantly improved accuracy. Further, the initial approximations, and the approximations arising from the first generalization, can be utilized as the inputs to a custom dynamical system to establish approximations with better convergence properties. Indicative results include those of a fourth order approximation, based on four sub-intervals, which leads to a relative error bound of 1.43 x 10-7 over the positive real line. Various approximations, that achieve the set relative error bounds of 10-4, 10-6, 10-10 and 10-16, over the positive real, are specified. Applications include, first, the definition of functions that are upper and lower bounds, of arbitrary accuracy, for the error function. Second, new series for the error function. Third, new sequences of approximations for exp(-x2) which have significantly higher convergence properties that a Taylor series approximation. Fourth, the definition of a complementary demarcation function eC(x) which satisfies the constraint eC(x)^2 + erf(x)^2 = 1. Fifth, arbitrarily accurate approximations for the power and harmonic distortion for a sinusoidal signal subject to a error function nonlinearity. Sixth, approximate expressions for the linear filtering of a step signal that is modelled by the error function.


title: On the Optimal Deployment of Power Beacons for Massive Wireless Energy Transfer

abstract: Wireless energy transfer (WET) is emerging as an enabling green technology for Internet of Things (IoT) networks. WET allows the IoT devices to wirelessly recharge their batteries with energy from external sources such as dedicated radio frequency transmitters called power beacons (PBs). In this paper, we investigate the optimal deployment of PBs that guarantees a network-wide energy outage constraint. Optimal positions for the PBs are determined by maximizing the average incident power for the worst location in the service area since no information about the sensor deployment is provided. Such network planning guarantees the fairest harvesting performance for all the IoT devices. Numerical simulations evidence that our proposed optimization framework improves the energy supply reliability compared to benchmark schemes. Additionally, we show that although both, the number of deployed PBs and the number of antennas per PB, introduce performance improvements, the former has a dominant role. Finally, our proposal allows to extend the coverage area while keeping the total power budget fixed, which additionally reduces the level of electromagnetic radiation in the vicinity of PBs.


title: Active Learning Methods for Efficient Hybrid Biophysical Variable Retrieval

abstract: Kernel-based machine learning regression algorithms (MLRAs) are potentially powerful methods for being implemented into operational biophysical variable retrieval schemes. However, they face difficulties in coping with large training datasets. With the increasing amount of optical remote sensing data made available for analysis and the possibility of using a large amount of simulated data from radiative transfer models (RTMs) to train kernel MLRAs, efficient data reduction techniques will need to be implemented. Active learning (AL) methods enable to select the most informative samples in a dataset. This letter introduces six AL methods for achieving optimized biophysical variable estimation with a manageable training dataset, and their implementation into a Matlab-based MLRA toolbox for semi-automatic use. The AL methods were analyzed on their efficiency of improving the estimation accuracy of leaf area index and chlorophyll content based on PROSAIL simulations. Each of the implemented methods outperformed random sampling, improving retrieval accuracy with lower sampling rates. Practically, AL methods open opportunities to feed advanced MLRAs with RTM-generated training data for development of operational retrieval models.


title: Multi-temporal and multi-source remote sensing image classification by nonlinear relative normalization

abstract: Remote sensing image classification exploiting multiple sensors is a very challenging problem: data from different modalities are affected by spectral distortions and mis-alignments of all kinds, and this hampers re-using models built for one image to be used successfully in other scenes. In order to adapt and transfer models across image acquisitions, one must be able to cope with datasets that are not co-registered, acquired under different illumination and atmospheric conditions, by different sensors, and with scarce ground references. Traditionally, methods based on histogram matching have been used. However, they fail when densities have very different shapes or when there is no corresponding band to be matched between the images. An alternative builds upon \emph{manifold alignment}. Manifold alignment performs a multidimensional relative normalization of the data prior to product generation that can cope with data of different dimensionality (e.g. different number of bands) and possibly unpaired examples. Aligning data distributions is an appealing strategy, since it allows to provide data spaces that are more similar to each other, regardless of the subsequent use of the transformed data. In this paper, we study a methodology that aligns data from different domains in a nonlinear way through {\em kernelization}. We introduce the Kernel Manifold Alignment (KEMA) method, which provides a flexible and discriminative projection map, exploits only a few labeled samples (or semantic ties) in each domain, and reduces to solving a generalized eigenvalue problem. We successfully test KEMA in multi-temporal and multi-source very high resolution classification tasks, as well as on the task of making a model invariant to shadowing for hyperspectral imaging.


title: Graph Coloring with Quantum Annealing

abstract: We develop a heuristic graph coloring approximation algorithm that uses the D-Wave 2X as an independent set sampler and evaluate its performance against a fully classical implementation. A randomly generated set of small but hard graph instances serves as our test set. Our performance analysis suggests limited quantum advantage in the hybrid quantum-classical algorithm. The quantum edge holds over multiple metrics and suggests that graph problem applications are a good fit for quantum annealers.


title: Ultramassive black holes in the most massive galaxies: $M_{\rm BH}-\sigma$ versus $M_{\rm BH}-R_{\rm b}$

abstract: [Abridged] We investigate the nature of the relations between black hole (BH) mass ($M_{\rm BH}$) and the central velocity dispersion ($\sigma$) and, for core-S\'ersic galaxies, the size of the depleted core ($R_{\rm b}$). Our sample of 144 galaxies with dynamically determined $M_{\rm BH}$ encompasses 24 core-S\'ersic galaxies, thought to be products of gas-poor mergers, and reliably identified based on high-resolution HST imaging. For core-S\'ersic galaxies -- i.e., combining normal-core ($R_{\rm b} < 0.5 $ kpc) and large-core galaxies ($R_{\rm b} \gtrsim 0.5$ kpc), we find that $M_{\rm BH}$ correlates remarkably well with $R_{\rm b}$ such that $M_{\rm BH} \propto R_{\rm b}^{1.20 \pm 0.14}$ (rms scatter in log $M_{\rm BH}$ of $\Delta_{\rm rms} \sim 0.29$ dex), confirming previous works on the same galaxies except three new ones. Separating the sample into S\'ersic, normal-core and large-core galaxies, we find that S\'ersic and normal-core galaxies jointly define a single log-linear $M_{\rm BH}-\sigma$ relation $M_{\rm BH} \propto \sigma^{ 4.88 \pm 0.29}$ with $\Delta_{\rm rms} \sim 0.47$ dex, however, at the high-mass end large-core galaxies (four with measured $M_{\rm BH}$) are offset upward from this relation by ($2.5-4) \times \sigma_{\rm s}$, explaining the previously reported steepening of the $M_{\rm BH}-\sigma$ relation for massive galaxies. Large-core spheroids have magnitudes $M_{V} \le -23.50$ mag, half-light radii Re $>$ 10 kpc and are extremely massive $M_{*} \ge 10^{12}M_{\odot}$. Furthermore, these spheroids tend to host ultramassive BHs ($M_{\rm BH} \ge 10^{10}M_{\odot}$) tightly connected with their $R_{\rm b}$ rather than $\sigma$. The less popular $M_{\rm BH}-R_{\rm b}$ relation exhibits $\sim$ 62% less scatter in log $M_{\rm BH}$ than the $M_{\rm BH}- \sigma$ relations.


title: Weyl-Invariant Gravity and the Nature of Dark Matter

abstract: The apparent missing mass in galaxies and galaxy clusters, commonly viewed as evidence for dark matter, could possibly originate from gradients in the gravitational coupling parameter, $G$, and active gravitational mass, $M_{act}$, rather than hypothetical beyond-the-standard-model particles. We argue that in (the weak field limit of) a Weyl-invariant extension of General Relativity, one can simply affect the change $\Phi_{b}(x)\rightarrow\Phi_{b}(x) + \Phi_{DM}(x)$, where $\Phi_{b}$ is the baryon-sourced potential and $\Phi_{DM}$ is the `excess' potential. This is compensated by gradients of $GM_{act}$ and a fractional increase of $O(-4\Phi_{DM}(x))$ in the baryon density, well below current detection thresholds on all relevant scales.


title: Quantum Technology for Economists

abstract: Research on quantum technology spans multiple disciplines: physics, computer science, engineering, and mathematics. The objective of this manuscript is to provide an accessible introduction to this emerging field for economists that is centered around quantum computing and quantum money. We proceed in three steps. First, we discuss basic concepts in quantum computing and quantum communication, assuming knowledge of linear algebra and statistics, but not of computer science or physics. This covers fundamental topics, such as qubits, superposition, entanglement, quantum circuits, oracles, and the no-cloning theorem. Second, we provide an overview of quantum money, an early invention of the quantum communication literature that has recently been partially implemented in an experimental setting. One form of quantum money offers the privacy and anonymity of physical cash, the option to transact without the involvement of a third party, and the efficiency and convenience of a debit card payment. Such features cannot be achieved in combination with any other form of money. Finally, we review all existing quantum speedups that have been identified for algorithms used to solve and estimate economic models. This includes function approximation, linear systems analysis, Monte Carlo simulation, matrix inversion, principal component analysis, linear regression, interpolation, numerical differentiation, and true random number generation. We also discuss the difficulty of achieving quantum speedups and comment on common misconceptions about what is achievable with quantum computing.


title: Rotation-Invariant Autoencoders for Signals on Spheres

abstract: Omnidirectional images and spherical representations of $3D$ shapes cannot be processed with conventional 2D convolutional neural networks (CNNs) as the unwrapping leads to large distortion. Using fast implementations of spherical and $SO(3)$ convolutions, researchers have recently developed deep learning methods better suited for classifying spherical images. These newly proposed convolutional layers naturally extend the notion of convolution to functions on the unit sphere $S^2$ and the group of rotations $SO(3)$ and these layers are equivariant to 3D rotations. In this paper, we consider the problem of unsupervised learning of rotation-invariant representations for spherical images. In particular, we carefully design an autoencoder architecture consisting of $S^2$ and $SO(3)$ convolutional layers. As 3D rotations are often a nuisance factor, the latent space is constrained to be exactly invariant to these input transformations. As the rotation information is discarded in the latent space, we craft a novel rotation-invariant loss function for training the network. Extensive experiments on multiple datasets demonstrate the usefulness of the learned representations on clustering, retrieval and classification applications.


title: Privacy-Preserving Synthetic Smart Meters Data

abstract: Power consumption data is very useful as it allows to optimize power grids, detect anomalies and prevent failures, on top of being useful for diverse research purposes. However, the use of power consumption data raises significant privacy concerns, as this data usually belongs to clients of a power company. As a solution, we propose a method to generate synthetic power consumption samples that faithfully imitate the originals, but are detached from the clients and their identities. Our method is based on Generative Adversarial Networks (GANs). Our contribution is twofold. First, we focus on the quality of the generated data, which is not a trivial task as no standard evaluation methods are available. Then, we study the privacy guarantees provided to members of the training set of our neural network. As a minimum requirement for privacy, we demand our neural network to be robust to membership inference attacks, as these provide a gateway for further attacks in addition to presenting a privacy threat on their own. We find that there is a compromise to be made between the privacy and the performance provided by the algorithm.


title: Circular geodesics in Kerr-Newman-Kasuya black hole

abstract: This article explores the characteristics of ergoregion, horizons and circular geodesics around a Kerr-Newman-Kasuya black hole. We investigate the effect of spin and dyonic charge parameters on ergoregion, event horizon and static limit surface of the said black hole. We observed that both electric, as well as magnetic charge parameters, results in decreasing the radii of event horizon and static limit, whereas increasing the area of ergoregion. The obtained results are compared with that acquired from Kerr and Schwarzschild black holes. Moreover, we figured out the photons orbit of circular null geodesics and studied the angular velocity of a particle within ergoregion.


title: Analyzing Finite Neural Networks: Can We Trust Neural Tangent Kernel Theory?

abstract: Neural Tangent Kernel (NTK) theory is widely used to study the dynamics of infinitely-wide deep neural networks (DNNs) under gradient descent. But do the results for infinitely-wide networks give us hints about the behavior of real finite-width ones? In this paper, we study empirically when NTK theory is valid in practice for fully-connected ReLU and sigmoid DNNs. We find out that whether a network is in the NTK regime depends on the hyperparameters of random initialization and the network's depth. In particular, NTK theory does not explain the behavior of sufficiently deep networks initialized so that their gradients explode as they propagate through the network's layers: the kernel is random at initialization and changes significantly during training in this case, contrary to NTK theory. On the other hand, in the case of vanishing gradients, DNNs are in the the NTK regime but become untrainable rapidly with depth. We also describe a framework to study generalization properties of DNNs, in particular the variance of network's output function, by means of NTK theory and discuss its limits.


title: Escape dynamics of active particles in multistable potentials

abstract: Rare transitions between long-lived metastable states underlie a great variety of physical, chemical and biological processes. Our quantitative understanding of reactive mechanisms has been driven forward by the insights of transition state theory. In particular, the dynamic framework developed by Kramers marks an outstanding milestone for the field. Its predictions, however, do not apply to systems driven by non-conservative forces or correlated noise histories. An important class of such systems are active particles, prominent in both biology and nanotechnology. Here, we trap a silica nanoparticle in a bistable potential. To emulate an active particle, we subject the particle to an engineered external force that mimics self-propulsion. We investigate the active particle's transition rate between metastable states as a function of friction and correlation time of the active force. Our experiments reveal the existence of an optimal correlation time where the transition rate is maximized. This novel \emph{active turnover} is reminiscent of the much celebrated Kramers turnover despite its fundamentally different origin. Our observations are quantitatively supported by a theoretical analysis of a one-dimensional model. Besides providing a deeper understanding of the escape dynamics of active particles in multistable potentials, our work establishes a new, versatile experimental platform to study particle dynamics in non-equilibrium settings.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.stat-mech physics.optics


title: Transfer Learning for Human Activity Recognition using Representational Analysis of Neural Networks

abstract: Human activity recognition (HAR) research has increased in recent years due to its applications in mobile health monitoring, activity recognition, and patient rehabilitation. The typical approach is training a HAR classifier offline with known users and then using the same classifier for new users. However, the accuracy for new users can be low with this approach if their activity patterns are different than those in the training data. At the same time, training from scratch for new users is not feasible for mobile applications due to the high computational cost and training time. To address this issue, we propose a HAR transfer learning framework with two components. First, a representational analysis reveals common features that can transfer across users and user-specific features that need to be customized. Using this insight, we transfer the reusable portion of the offline classifier to new users and fine-tune only the rest. Our experiments with five datasets show up to 43% accuracy improvement and 66% training time reduction when compared to the baseline without using transfer learning. Furthermore, measurements on the Nvidia Jetson Xavier-NX hardware platform reveal that the power and energy consumption decrease by 43% and 68%, respectively, while achieving the same or higher accuracy as training from scratch.


title: The Role that Gaiters, Masks and Face Shields Can Play in Limiting the Transmission of Respiratory Droplets

abstract: The efficacy of face masks, neck gaiters and face shields are predicted and contrasted. In particular, a Lattice Spring Model of a neck gaiter serves as the input to a Lattice Boltzmann simulation. The Lattice Boltzmann method is used to capture the fluid dynamics both through and around various face coverings. The evaporation and transport of respiratory droplets is simulated in this flow field, and the efficacy of the different face coverings at stopping the spread of respiratory droplets are elucidated and contrasted. In agreement with recent experimental studies, we predict that neck gaiters are at least as effective as cloth masks, and both are far superior to a face shield.


title: Implementing strong interference in ultrathin film top absorbers for tandem solar cells

abstract: Strong interference in ultrathin film semiconductor absorbers on metallic back reflectors has been shown to enhance the light harvesting efficiency of solar cell materials. However, metallic back reflectors are not suitable for tandem cell configurations because photons cannot be transmitted through the device. Here, we introduce a method to implement strong interference in ultrathin film top absorbers in a tandem cell configuration through use of distributed Bragg reflectors (DBRs). We showcase this by designing and fabricating a photoelectrochemical-photovoltaic (PEC-PV) stacked tandem cell in a V-shaped configuration where short wavelength photons are reflected back to the photoanode material (hematite, Fe2O3), whereas long wavelength photons are transmitted to the bottom silicon PV cell. We employ optical simulations to determine the optimal thicknesses of the DBR layers and the V-shape angle to maximize light absorption in the ultrathin (10 nm thick) hematite film. The DBR spectral response can be tailored to allow for a more than threefold enhancement in absorbed photons compared to a layer of the same thickness on transparent current collectors. Using a DBR to couple a bottom silicon PV cell with an ultrathin hematite top PEC cell, we demonstrate unassisted solar water splitting and show that DBRs can be designed to enhance strong interference in ultrathin films while enabling stacked tandem cell configuration.


title: Kadomtsev-Petviashvili turning points and CKP hierarchy

abstract: A characterization of the Kadomtsev-Petviashvili hierarchy of type C (CKP) in terms of the KP tau-function is given. Namely, we prove that the CKP hierarchy can be identified with the restriction of odd times flows of the KP hierarchy on the locus of turning points of the second flow. The notion of CKP tau-function is clarified and connected with the KP tau function. Algebraic-geometrical solutions and in particular elliptic solutions are discussed in detail. The new identity for theta-functions of curves with holomorphic involution having fixed points is obtained.


title: A Novel Transformation Approach of Shared-link Coded Caching Schemes for Multiaccess Networks

abstract: This paper considers the multiaccess coded caching systems formulated by Hachem et al., including a central server containing $N$ files connected to $K$ cache-less users through an error-free shared link, and $K$ cache-nodes, each equipped with a cache memory size of $M$ files. Each user has access to $L$ neighbouring cache-nodes with a cyclic wrap-around topology. The coded caching scheme proposed by Hachem et al. suffers from the case that $L$ does not divide $K$, where the needed number of transmissions (a.k.a. load) is at most four times the load expression for the case where $L$ divides $K$. Our main contribution is to propose a novel {\it transformation} approach to smartly extend the schemes satisfying some conditions for the well known shared-link caching systems to the multiaccess caching systems. Then we can get many coded caching schemes with different subpacketizations for multiaccess coded caching system. These resulting schemes have the maximum local caching gain (i.e., the cached contents stored at any $L$ neighbouring cache-nodes are different such that the number of retrieval packets by each user from the connected cache-nodes is maximal) and the same coded caching gain as the original schemes. Applying the transformation approach to the well-known shared-link coded caching scheme proposed by Maddah-Ali and Niesen, we obtain a new multiaccess coded caching scheme that achieves the same load as the scheme of Hachem et al. but for any system parameters. Under the constraint of the cache placement used in this new multiaccess coded caching scheme, our delivery strategy is approximately optimal when $K$ is sufficiently large. Finally, we also show that the transmission load of the proposed scheme can be further reduced by compressing the multicast message.


title: On the observability of field-assisted birefringent Delbr\"uck scattering

abstract: We consider the scattering of an x-ray free-electron laser (XFEL) beam on the superposition of a strong magnetic field $\bf{B}_{\rm ext}$ with the Coulomb field $\bf{E}_{\rm ext}$ of a nucleus with charge number $Z$. In contrast to pure Delbr\"uck scattering (Coulomb field only), the magnetic field $\bf{B}_{\rm ext}$ introduces an asymmetry (i.e., polarization dependence) and renders the effective interaction volume quite large, while the nuclear Coulomb field facilitates a significant momentum transfer $\Delta\bf k$. For a field strength of $B_{\rm ext}=10^6~\rm T$ (corresponding to an intensity of order $10^{22}~\rm W/cm^2$) and an XFEL frequency of 24 keV, we find a differential cross section $d\sigma/d\Omega\sim10^{-25}~Z^2/(\Delta{\bf k})^2$ in forward direction for one nucleus. Thus, this effect might be observable in the near future at facilities such as the Helmholtz International Beamline for Extreme Fields (HIBEF) at the European XFEL.


title: Occupational segregation in a Roy model with composition preferences

abstract: We propose a model of labor market sector self-selection that combines comparative advantage, as in the Roy model, and sector composition preference. Two groups choose between two sectors based on heterogeneous potential incomes and group compositions in each sector. Potential incomes incorporate group specific human capital accumulation and wage discrimination. Composition preferences are interpreted as reflecting group specific amenity preferences as well as homophily and aversion to minority status. We show that occupational segregation is amplified by the composition preferences and we highlight a resulting tension between redistribution and diversity. The model also exhibits tipping from extreme compositions to more balanced ones. Tipping occurs when a small nudge, associated with affirmative action, pushes the system to a very different equilibrium, and when the set of equilibria changes abruptly when a parameter governing the relative importance of pecuniary and composition preferences crosses a threshold.


title: Testing the black-hole area law with GW150914

abstract: We present observational confirmation of Hawking's black-hole area theorem based on data from GW150914, finding agreement with the prediction with 97% (95%) probability when we model the ringdown including (excluding) overtones of the quadrupolar mode. We obtain this result from a new time-domain analysis of the pre- and postmerger data. We also confirm that the inspiral and ringdown portions of the signal are consistent with the same remnant mass and spin, in agreement with general relativity.


title: Finite Quantum Field Theory and Renormalization Group

abstract: Renormalization group methods are applied to a scalar field within a finite, nonlocal quantum field theory formulated perturbatively in Euclidean momentum space. It is demonstrated that the triviality problem in scalar field theory, the Higgs boson mass hierarchy problem and the stability of the vacuum do not arise as issues in the theory. The scalar Higgs field has no Landau pole.


title: A Concentration Inequality for the Facility Location Problem

abstract: We give a concentration inequality for a stochastic version of the facility location problem. We show the objective $C_n = \min_{F \subseteq [0,1]^2}|F|+\sum_{x\in X}\min_{f\in F}\|x-f\|$ is concentrated in an interval of length $O(n^{1/6})$ and $\E[C_n]=\Theta(n^{2/3})$ if the input $X$ consists of i.i.d. uniform points in the unit square. Our main tool is to use a geometric quantity, previously used in the design of approximation algorithms for the facility location problem, to analyze a martingale process. Many of our techniques generalize to other settings.


title: $\alpha$-$\beta$ phase transition of zirconium predicted by on-the-fly machine-learned force field

abstract: The accurate prediction of solid-solid structural phase transitions at finite temperature is a challenging task, since the dynamics is so slow that direct simulations of the phase transitions by first-principles (FP) methods are typically not possible. Here, we study the $\alpha$-$\beta$ phase transition of Zr at ambient pressure by means of on-the-fly machine-learned force fields. These are automatically generated during FP molecular dynamics (MD) simulations without the need of human intervention, while retaining almost FP accuracy. Our MD simulations successfully reproduce the first-order displacive nature of the phase transition, which is manifested by an abrupt jump of the volume and a cooperative displacement of atoms at the phase transition temperature. The phase transition is further identified by the simulated x-ray powder diffraction, and the predicted phase transition temperature is in reasonable agreement with experiment. Furthermore, we show that using a singular value decomposition and pseudo inversion of the design matrix generally improves the machine-learned force field compared to the usual inversion of the squared matrix in the regularized Bayesian regression.


title: Lorentz Symmetry Group, Retardation, Intergalactic Mass Depletion and Mechanisms Leading to Galactic Rotation Curves

abstract: The general theory of relativity (GR) is symmetric under smooth coordinate transformations, also known as diffeomorphisms. The general coordinate transformation group has a linear subgroup denoted as the Lorentz group of symmetry, which is also maintained in the weak field approximation to GR. The dominant operator in the weak field equation of GR is thus the d'Alembert (wave) operator, which has a retarded potential solution. Galaxies are huge physical systems with dimensions of many tens of thousands of light years. Thus, any change at the galactic center will be noticed at the rim only tens of thousands of years later. Those retardation effects are neglected in the present day galactic modelling used to calculate rotational velocities of matter in the rims of the galaxy and surrounding gas. The significant differences between the predictions of Newtonian instantaneous action at a distance and observed velocities are usually explained by either assuming dark matter or by modifying the laws of gravity (MOND). In this paper, we will show that, by taking general relativity seriously without neglecting retardation effects, one can explain the radial velocities of galactic matter in the M33 galaxy without postulating dark matter.


title: Quantum mechanics of stationary states of particles in a space-time of classical black holes

abstract: We consider interactions of scalar particles, photons, and fermions in Schwarzschild, Reissner-Nordstr\"om, Kerr, and Kerr-Newman gravitational and electromagnetic fields with a zero and nonzero cosmological constant. We also consider interactions of scalar particles, photons, and fermions with nonextremal rotating charged black holes in a minimal five-dimensional gauge supergravity. We analyze the behavior of effective potentials in second-order relativistic Schr\"odinger-type equations. In all cases, we establish the existence of the regime of particle "falling" on event horizons. An alternative can be collapsars with fermions in stationary bound states without a regime of particles "falling".


title: Reduced motion artifacts and speed improvements in enhanced line-scanning fiber bundle endomicroscopy

abstract: Significance: Confocal laser scanning enables optical sectioning in fiber bundle endomicroscopy but limits the frame rate. To be able to better explore tissue morphology it is useful to stitch sequentially acquired frames into a mosaic. However, low frame rates limit the maximum probe translation speed. Line-scanning confocal endomicroscopy provides higher frame rates, but residual out-of-focus light degrades images. Subtraction based approaches can suppress this residue at the expense of introducing motion artifacts. Aim: To generate high frame rate endomicroscopy images with improved optical sectioning, we develop a high-speed subtraction method that only requires the acquisition of a single camera frame. Approach: The rolling shutter of a CMOS camera acts as both the aligned and offset detector slits required for subtraction-based sectioning enhancement. Two images of the bundle are formed on different regions of the camera, allowing both images to be acquired simultaneously. Results: We confirm improved optical sectioning compared to conventional line-scanning, particularly far from focus, and show that motion artifacts are not introduced. We demonstrate high-speed mosaicing at frame rates of up to 240 Hz. Conclusion: High-speed acquisition of optically sectioned images using the new subtraction based approach leads to improved mosaicing at high frame rates.


title: Relativistic Landau levels for a fermion-antifermion pair interacting through Dirac oscillator interaction

abstract: We introduce a unique model for a fermion antifermion pair interacting through Dirac oscillator interaction in the presence of external uniform magnetic field. In order to acquire a non perturbative energy spectrum for such a system we solve the corresponding form of a fully covariant two body Dirac equation. The dynamic symmetry of the system allows to study in three dimensions and the corresponding equation leads a $4\times4$ dimensional matrix equation for such a static and spinless composite system. We can obtain an exact solution of the matrix equation and arrive at a spectrum (in closed form) in energy domain. The obtained energy spectrum shows that the composite system that under scrutiny behaves like a single relativistic quantum oscillator. As a result, we obtain relativistic Landau levels of a fermion antifermion pair interacting through Dirac oscillator interaction and determine the components of the corresponding bi-spinor. We think that our results can provide enlightening informations about the quark antiquark systems and thus mass formula for mesons.


title: Bayesian Learning of LF-MMI Trained Time Delay Neural Networks for Speech Recognition

abstract: Discriminative training techniques define state-of-the-art performance for automatic speech recognition systems. However, they are inherently prone to overfitting, leading to poor generalization performance when using limited training data. In order to address this issue, this paper presents a full Bayesian framework to account for model uncertainty in sequence discriminative training of factored TDNN acoustic models. Several Bayesian learning based TDNN variant systems are proposed to model the uncertainty over weight parameters and choices of hidden activation functions, or the hidden layer outputs. Efficient variational inference approaches using a few as one single parameter sample ensure their computational cost in both training and evaluation time comparable to that of the baseline TDNN systems. Statistically significant word error rate (WER) reductions of 0.4%-1.8% absolute (5%-11% relative) were obtained over a state-of-the-art 900 hour speed perturbed Switchboard corpus trained baseline LF-MMI factored TDNN system using multiple regularization methods including F-smoothing, L2 norm penalty, natural gradient, model averaging and dropout, in addition to i-Vector plus learning hidden unit contribution (LHUC) based speaker adaptation and RNNLM rescoring. Consistent performance improvements were also obtained on a 450 hour HKUST conversational Mandarin telephone speech recognition task. On a third cross domain adaptation task requiring rapidly porting a 1000 hour LibriSpeech data trained system to a small DementiaBank elderly speech corpus, the proposed Bayesian TDNN LF-MMI systems outperformed the baseline system using direct weight fine-tuning by up to 2.5\% absolute WER reduction.


title: Retrieval of Case 2 Water Quality Parameters with Machine Learning

abstract: Water quality parameters are derived applying several machine learning regression methods on the Case2eXtreme dataset (C2X). The used data are based on Hydrolight in-water radiative transfer simulations at Sentinel-3 OLCI wavebands, and the application is done exclusively for absorbing waters with high concentrations of coloured dissolved organic matter (CDOM). The regression approaches are: regularized linear, random forest, Kernel ridge, Gaussian process and support vector regressors. The validation is made with and an independent simulation dataset. A comparison with the OLCI Neural Network Swarm (ONSS) is made as well. The best approached is applied to a sample scene and compared with the standard OLCI product delivered by EUMETSAT/ESA


title: Invariant scalar-flat K\"ahler metrics on line bundles over generalized flag varieties

abstract: Let $G$ be a simply-connected semisimple compact Lie group, $X$ a compact K\"ahler manifold homogeneous under $G$, and $L$ a negative $G$-equivariant holomorphic line bundle over $X$. We prove that all $G$-invariant K\"ahler metrics on the total space of $L$ arise from the Calabi ansatz. Using this, we then show that there exists a unique $G$-invariant scalar-flat K\"ahler metric in each K\"ahler class of $L$.


title: Step Momentum Operator

abstract: In the present study, the concept of a quantum particle with step momentum is introduced. The energy eigenvalues and eigenfunctions of such particles are obtained in the context of the generalized momentum operator, proposed recently in [13, 14]. While, the number of bound states with real energy for the particles with Hermitian step momentum inside a square well is infinite, it is finite for a particle with \mathcal{PT} -symmetric momentum.


title: Experimental decoy state BB84 quantum key distribution through a turbulent channel

abstract: In free-space Quantum Key Distribution in turbulent conditions, scattering and beam wandering cause intensity fluctuations which increase the detected signal-to-noise ratio. This effect can be mitigated by rejecting received bits when the channel's transmittance is below a threshold. Thus, the overall error rate is reduced and the secure key rate increases despite the deletion of bits. In Decoy State BB84 QKD, several methods to find the ideal threshold have already been proposed. One promising method is the Prefixed-Threshold Real-time Selection (P-RTS) where a cutoff can be chosen prior to data collection and independently of the intensity distribution. In this work we perform finite-size Decoy State BB84 QKD in a laboratory setting where we simulate the atmospheric turbulence using an acousto-optical modulator. We show that P-RTS can yield considerably higher secure key rates, especially at higher losses. In addition, we verify that P-RTS gives a robust cutoff even when the predicted transmittance does not match the actual transmittance.


title: Toroidalization of Locally Toroidal Morphisms

abstract: The problem of toroidalization is to construct a toroidal lifting of a dominant morphism $\varphi:X\to Y$ of algebraic varieties by blowing up in the target and domain. This paper contains a solution to this problem when $\varphi$ is locally toroidal.


title: Portfolio Optimisation within a Wasserstein Ball

abstract: We study the problem of active portfolio management where an investor aims to outperform a benchmark strategy's risk profile while not deviating too far from it. Specifically, an investor considers alternative strategies whose terminal wealth lie within a Wasserstein ball surrounding a benchmark's -- being distributionally close -- and that have a specified dependence/copula -- tying state-by-state outcomes -- to it. The investor then chooses the alternative strategy that minimises a distortion risk measure of terminal wealth. In a general (complete) market model, we prove that an optimal dynamic strategy exists and provide its characterisation through the notion of isotonic projections. We further propose a simulation approach to calculate the optimal strategy's terminal wealth, making our approach applicable to a wide range of market models. Finally, we illustrate how investors with different copula and risk preferences invest and improve upon the benchmark using the Tail Value-at-Risk, inverse S-shaped, and lower- and upper-tail distortion risk measures as examples. We find that investors' optimal terminal wealth distribution has larger probability masses in regions that reduce their risk measure relative to the benchmark while preserving the benchmark's structure.


title: Sitnikov-type solution for the motion of infinitesimal mass in BiER4BP

abstract: In this paper, we present a new ansatz for approximated solving equations of motion of the infinitesimal mass m in case of bi-elliptic restricted problem of four bodies (BiER4BP) (where three primaries M1, M2, M3 are rotating around their common centre of mass on elliptic orbits with hierarchical configuration M3 < M2 << M1). A new type of the solving procedure is implemented here to obtain the coordinates of the infinitesimal mass m. Meanwhile, the system of equations of motion has been successfully explored with respect to the existence of semi-analytical (approximated) way for presentation of the solution. We obtain as follows: 1) the solution for coordinates {x, y} = {0, 0} is approximately satisfied both the first and second equations of motion if we take into consideration assumption {M3, M2} << M1, 2) the expression for coordinate z(f) is given by the equation of 2-nd order, which describes Sitnikov-type approximated solution. It means that test particle is moving along the z-axis, outward the common barycenter of the system (but perpendicular to the plane of the mutual rotation of all the primaries).


title: Dissipation-induced antiferromagnetic-like frustration in coupled photonic resonators

abstract: We propose a photonic quantum simulator for anti-ferromagnetic spin systems based on reservoir engineering. We consider a scheme where quadratically driven dissipative Kerr cavities are indirectly coupled via lossy ancillary cavities. We show that the ancillary cavities can produce an effective dissipative and Hamiltonian anti-ferromagnetic-like coupling between the cavities. By solving the master equation for a triangular cavity configuration, we demonstrate that the non-equilibrium steady state of the system bears full analogy with the ground state of an antiferromagnetic Ising model, exhibiting key signatures of frustration. We show that when the effective photon hopping amplitude is zero, the engineered non-local dissipation alone is capable of inducing antiferromagnetic interaction and frustration. This simple scheme can be generalised to arbitrary lattice geometries, providing a fully controllable recipe for simulating antiferromagnetism and frustration on a controlled quantum optical platform.


title: Reciprocity law of finite Galois extension fields using Jacobian Varitey

abstract: For finite Galois extension fields defined by odd degree irreducible polynomials over algebraic integer ring, we observe "Reciprocity Law" through Jacobian Variety by embedding all roots of the polynomials into 2-torsion points of Jacobian Variety. Furthermore, Galois group of the minimal splitting field of such a polynomial is a subgroup of general linear group with coefficient $\mathbb{F}_{2}$.


title: Revised planet brightness temperatures using the Planck/LFI 2018 data release

abstract: We present new estimates of the brightness temperatures of Jupiter, Saturn, Uranus, and Neptune based on the measurements carried in 2009--2013 by PLANCK/LFI at 30, 44, and 70 GHz and released to the public in 2018. This work extends the results presented in the 2013 and 2015 PLANCK/LFI Calibration Papers, based on the data acquired in 2009--2011. PLANCK observed each planet up to eight times during the nominal mission. We processed time-ordered data from the 22 LFI radiometers to derive planet antenna temperatures for each planet and transit. We accounted for the beam shape, radiometer bandpasses, and several systematic effects. We compared our results with the results from the ninth year of WMAP, PLANCK/HFI observations, and existing data and models for planetary microwave emissivity. For Jupiter, we obtain Tb = 144.9, 159.8, 170.5 K (+/- 0.2 K at 1 sigma, with temperatures expressed using the Rayleigh-Jeans scale) at 30, 44 and 70 GHz, respectively, or equivalently a band averaged Planck temperature TbBA=144.7$, 160.3, 171.2 K in good agreement with WMAP and existing models. A slight excess at 30 GHz with respect to models is interpreted as an effect of synchrotron emission. Our measures for Saturn agree with the results from WMAP for rings Tb = 9.2 +/- 1.4, 12.6 +/- 2.3, 16.2 +/- 0.8 K, while for the disc we obtain Tb = 140.0 +/- 1.4, 147.2 +/- 1.2, 150.2 +/- 0.4 K, or equivalently a TbBA=139.7, 147.8, 151.0 K. Our measures for Uranus (Tb = 152 +/- 6, 145 +/- 3, 132.0 +/- 2 K, or TbBA=152, 145, 133 K and Neptune Tb = 154 +/- 11, 148 +/- 9, 128 +/- 3 K, or TbBA=154 , 149, 128 K) agree closely with WMAP and previous data in literature.


title: Gibbs posterior concentration rates under sub-exponential type losses

abstract: Bayesian posterior distributions are widely used for inference, but their dependence on a statistical model creates some challenges. In particular, there may be lots of nuisance parameters that require prior distributions and posterior computations, plus a potentially serious risk of model misspecification bias. Gibbs posterior distributions, on the other hand, offer direct, principled, probabilistic inference on quantities of interest through a loss function, not a model-based likelihood. Here we provide simple sufficient conditions for establishing Gibbs posterior concentration rates when the loss function is of a sub-exponential type. We apply these general results in a range of practically relevant examples, including mean regression, quantile regression, and sparse high-dimensional classification. We also apply these techniques in an important problem in medical statistics, namely, estimation of a personalized minimum clinically important difference.


title: Business Cycles as Collective Risk Fluctuations

abstract: We suggest use continuous numerical risk grades [0,1] of R for a single risk or the unit cube in Rn for n risks as the economic domain. We consider risk ratings of economic agents as their coordinates in the economic domain. Economic activity of agents, economic or other factors change agents risk ratings and that cause motion of agents in the economic domain. Aggregations of variables and transactions of individual agents in small volume of economic domain establish the continuous economic media approximation that describes collective variables, transactions and their flows in the economic domain as functions of risk coordinates. Any economic variable A(t,x) defines mean risk XA(t) as risk weighted by economic variable A(t,x). Collective flows of economic variables in bounded economic domain fluctuate from secure to risky area and back. These fluctuations of flows cause time oscillations of macroeconomic variables A(t) and their mean risks XA(t) in economic domain and are the origin of any business and credit cycles. We derive equations that describe evolution of collective variables, transactions and their flows in the economic domain. As illustration we present simple self-consistent equations of supply-demand cycles that describe fluctuations of supply, demand and their mean risks.


title: Dynamical diffeomorphisms

abstract: We construct a general effective dynamics for diffeomorphisms of spacetime, in a fixed external metric. Though related to familiar models of scalar fields as coordinates, our models have subtly different properties, both at kinematical and dynamical level. The energy-momentum tensor consists of two independently conserved parts. The background solution is the identity diffeomorphism and the energy-momentum tensor of this solution gives rise to an effective cosmological constant.


title: Equidistribution of primitive lattices in $\mathbb{R}^n$

abstract: We count primitive lattices of rank $d$ inside $\mathbb{Z}^{n}$ as their covolume tends to infinity, with respect to certain parameters of such lattices. These parameters include, for example, the subsapce that a lattice spans, namely its projection to the Grassmannian; its homothety class; and its equivalence class modulo rescaling and rotation, often referred to as a shape. We add to a prior work of Schmidt by allowing sets in the spaces of parameters that are general enough to conclude joint equidistribution of these parameters. In addition to the primitive $d$-lattices themselves, we also consider their orthogonal complements in $\mathbb{Z}^{n}$, and show that the equidistribution occurs jointly for primitive lattices and their orthogonal complements. Finally, our asymptotic formulas for the number of primitive lattices include an explicit error term.


title: Halfspace type theorems for self-shrinkers in arbitrary codimension

abstract: In this paper, we generalize some halfspace type theorems for self-shrinkers of codimension 1 to the case of arbitrary codimension.


title: Graph-based open-ended survey on concerns related to COVID-19

abstract: The COVID-19 pandemic is an unprecedented public health crisis with broad social and economic consequences. We conducted four surveys between April and August 2020 using the graph-based open-ended survey (GOS) framework, and investigated the most pressing concerns and issues for the general public in Japan. The GOS framework is a hybrid of the two traditional survey frameworks that allows respondents to post their opinions in a free-format style, which can subsequently serve as one of the choice items for other respondents, just as in a multiple-choice survey. As a result, this framework generates an opinion graph that relates opinions and respondents. We can also construct annotated opinion graphs to achieve a higher resolution. By clustering the annotated opinion graphs, we revealed the characteristic evolution of the response patterns as well as the interconnectedness and multi-faceted nature of opinions. Substantively, our notable finding is that "social pressure," not "infection risk," was one of the major concerns of our respondents. Social pressure refers to criticism and discrimination that they anticipate receiving from others should they contract COVID-19. It is possible that the collectivist nature of Japanese culture coupled with the government's policy of relying on personal responsibility to combat COVID-19 explains some of the above findings, as the latter has led to the emergence of vigilantes. The presence of mutual surveillance can contribute to growing skepticism toward others as well as fear of ostracism, which may have negative consequences at both the societal and individual levels.


title: Emotive Response to a Hybrid-Face Robot and Translation to Consumer Social Robots

abstract: We introduce the conceptual formulation, design, fabrication, control and commercial translation with IoT connection of a hybrid-face social robot and validation of human emotional response to its affective interactions. The hybrid-face robot integrates a 3D printed faceplate and a digital display to simplify conveyance of complex facial movements while providing the impression of three-dimensional depth for natural interaction. We map the space of potential emotions of the robot to specific facial feature parameters and characterise the recognisability of the humanoid hybrid-face robot's archetypal facial expressions. We introduce pupil dilation as an additional degree of freedom for conveyance of emotive states. Human interaction experiments demonstrate the ability to effectively convey emotion from the hybrid-robot face to human observers by mapping their neurophysiological electroencephalography (EEG) response to perceived emotional information and through interviews. Results show main hybrid-face robotic expressions can be discriminated with recognition rates above 80% and invoke human emotive response similar to that of actual human faces as measured by the face-specific N170 event-related potentials in EEG. The hybrid-face robot concept has been modified, implemented, and released in the commercial IoT robotic platform Miko (My Companion), an affective robot with facial and conversational features currently in use for human-robot interaction in children by Emotix Inc. We demonstrate that human EEG responses to Miko emotions are comparative to neurophysiological responses for actual human facial recognition. Finally, interviews show above 90% expression recognition rates in our commercial robot. We conclude that simplified hybrid-face abstraction conveys emotions effectively and enhances human-robot interaction.


title: SSCNav: Confidence-Aware Semantic Scene Completion for Visual Semantic Navigation

abstract: This paper focuses on visual semantic navigation, the task of producing actions for an active agent to navigate to a specified target object category in an unknown environment. To complete this task, the algorithm should simultaneously locate and navigate to an instance of the category. In comparison to the traditional point goal navigation, this task requires the agent to have a stronger contextual prior to indoor environments. We introduce SSCNav, an algorithm that explicitly models scene priors using a confidence-aware semantic scene completion module to complete the scene and guide the agent's navigation planning. Given a partial observation of the environment, SSCNav first infers a complete scene representation with semantic labels for the unobserved scene together with a confidence map associated with its own prediction. Then, a policy network infers the action from the scene completion result and confidence map. Our experiments demonstrate that the proposed scene completion module improves the efficiency of the downstream navigation policies. Video, code, and data: https://sscnav.cs.columbia.edu/


title: Correlated electrons in a Crystalline Topological insulator

abstract: We study the effects of e-e interaction in a 3D Crystalline Topological Insulator by adding on-site repulsion to the single-particle Hamiltonian and solving the many-body problem within Cluster Perturbation Theory. The goal is to clarify how many body effects modify the topological phase that stems from the crystal symmetries. Tuning the strength of the on-site interaction we show that band inversion disappears in the bulk and surface states loose their metallic character.


title: Human Motion Tracking by Registering an Articulated Surface to 3-D Points and Normals

abstract: We address the problem of human motion tracking by registering a surface to 3-D data. We propose a method that iteratively computes two things: Maximum likelihood estimates for both the kinematic and free-motion parameters of a kinematic human-body representation, as well as probabilities that the data are assigned either to a body part, or to an outlier cluster. We introduce a new metric between observed points and normals on one side, and a parameterized surface on the other side, the latter being defined as a blending over a set of ellipsoids. We claim that this metric is well suited when one deals with either visual-hull or visual-shape observations. We illustrate the method by tracking human motions using sparse visual-shape data (3-D surface points and normals) gathered from imperfect silhouettes.


title: Digital Gimbal: End-to-end Deep Image Stabilization with Learnable Exposure Times

abstract: Mechanical image stabilization using actuated gimbals enables capturing long-exposure shots without suffering from blur due to camera motion. These devices, however, are often physically cumbersome and expensive, limiting their widespread use. In this work, we propose to digitally emulate a mechanically stabilized system from the input of a fast unstabilized camera. To exploit the trade-off between motion blur at long exposures and low SNR at short exposures, we train a CNN that estimates a sharp high-SNR image by aggregating a burst of noisy short-exposure frames, related by unknown motion. We further suggest learning the burst's exposure times in an end-to-end manner, thus balancing the noise and blur across the frames. We demonstrate this method's advantage over the traditional approach of deblurring a single image or denoising a fixed-exposure burst on both synthetic and real data.


title: Finding equilibrium in two-stage traffic assignment model

abstract: Authors describe a two-stage traffic assignment model. It contains of two blocks. The first block consists of model for calculating correspondence (demand) matrix, whereas the second block is a traffic assignment model. The first model calculates a matrix of correspondences using a matrix of transport costs. It characterizes the required volumes of movement from one area to another. The second model describes how exactly the needs for displacement, specified by the correspondence matrix, are distributed along the possible paths. It works on the basis of the Nash--Wardrop equilibrium (each driver chooses the shortest path). Knowing the ways of distribute flows along the paths, it is possible to calculate the cost matrix. Equilibrium in a two-stage model is a fixed point in the sequence of these two models. The article proposes a method of reducing the problem of finding the equilibrium to the problem of the convex non-smooth optimization. Also a numerical method for solving the obtained optimization problem is proposed. Numerical experiments were carried out for the small towns.


title: A Geometric Framework for Pitch Estimation on Acoustic Musical Signals

abstract: This paper presents a geometric approach to pitch estimation (PE)-an important problem in Music Information Retrieval (MIR), and a precursor to a variety of other problems in the field. Though there exist a number of highly-accurate methods, both mono-pitch estimation and multi-pitch estimation (particularly with unspecified polyphonic timbre) prove computationally and conceptually challenging. A number of current techniques, whilst incredibly effective, are not targeted towards eliciting the underlying mathematical structures that underpin the complex musical patterns exhibited by acoustic musical signals. Tackling the approach from both a theoretical and experimental perspective, we present a novel framework, a basis for further work in the area, and results that (whilst not state of the art) demonstrate relative efficacy. The framework presented in this paper opens up a completely new way to tackle PE problems, and may have uses both in traditional analytical approaches, as well as in the emerging machine learning (ML) methods that currently dominate the literature.


title: Accomplice Manipulation of the Deferred Acceptance Algorithm

abstract: The deferred acceptance algorithm is an elegant solution to the stable matching problem that guarantees optimality and truthfulness for one side of the market. Despite these desirable guarantees, it is susceptible to strategic misreporting of preferences by the agents on the other side. We study a novel model of strategic behavior under the deferred acceptance algorithm: manipulation through an accomplice. Here, an agent on the proposed-to side (say, a woman) partners with an agent on the proposing side -- an accomplice -- to manipulate on her behalf (possibly at the expense of worsening his match). We show that the optimal manipulation strategy for an accomplice comprises of promoting exactly one woman in his true list (i.e., an inconspicuous manipulation). This structural result immediately gives a polynomial-time algorithm for computing an optimal accomplice manipulation. We also study the conditions under which the manipulated matching is stable with respect to the true preferences. Our experimental results show that accomplice manipulation outperforms self manipulation both in terms of the frequency of occurrence as well as the quality of matched partners.


title: Coxeter factorizations with generalized Jucys-Murphy weights and Matrix Tree theorems for reflection groups

abstract: We prove universal (case-free) formulas for the weighted enumeration of factorizations of Coxeter elements into products of reflections valid in any well-generated reflection group $W$, in terms of the spectrum of an associated operator, the $W$-Laplacian. This covers in particular all finite Coxeter groups. The results of this paper include generalizations of the Matrix Tree and Matrix Forest theorems to reflection groups, and cover reduced (shortest length) as well as arbitrary length factorizations. Our formulas are relative to a choice of weighting system that consists of $n$ free scalar parameters and is defined in terms of a tower of parabolic subgroups. To study such systems we introduce (a class of) variants of the Jucys-Murphy elements for every group, from which we define a new notion of `tower equivalence' of virtual characters. A main technical point is to prove the tower equivalence between virtual characters naturally appearing in the problem, and exterior products of the reflection representation of $W$. Finally we study how this $W$-Laplacian matrix we introduce can be used in other problems in Coxeter combinatorics. We explain how it defines analogues of trees for $W$ and how it relates them to Coxeter factorizations, we give new numerological identities between the Coxeter number of $W$ and those of its parabolic subgroups, and finally, when $W$ is a Weyl group, we produce a new, explicit formula for the volume of the corresponding root zonotope.


title: Numerical analysis of a wave equation for lossy media obeying a frequency power law

abstract: We study a wave equation with a nonlocal time fractional damping term that models the effects of acoustic attenuation characterized by a frequency dependence power law. First we prove existence of a unique solution to this equation with particular attention paid to the handling of the fractional derivative. Then we derive an explicit time stepping scheme based on the finite element method in space and a combination of convolution quadrature and second order central differences in time. We conduct a full error analysis of the mixed time discretization and in turn the fully space time discretized scheme. Error estimates are given for both smooth solutions and solutions with a singularity at $t = 0$ of a type that is typical for equations involving fractional time-derivatives. A number of numerical results are presented to support the error analysis.


title: Minimizing Spectral Risk Measures Applied to Markov Decision Processes

abstract: We study the minimization of a spectral risk measure of the total discounted cost generated by a Markov Decision Process (MDP) over a finite or infinite planning horizon. The MDP is assumed to have Borel state and action spaces and the cost function may be unbounded above. The optimization problem is split into two minimization problems using an infimum representation for spectral risk measures. We show that the inner minimization problem can be solved as an ordinary MDP on an extended state space and give sufficient conditions under which an optimal policy exists. Regarding the infinite dimensional outer minimization problem, we prove the existence of a solution and derive an algorithm for its numerical approximation. Our results include the findings in B\"auerle and Ott (2011) in the special case that the risk measure is Expected Shortfall. As an application, we present a dynamic extension of the classical static optimal reinsurance problem, where an insurance company minimizes its cost of capital.


title: Your College Dorm and Dormmates: Fair Resource Sharing with Externalities

abstract: We study a fair resource sharing problem, where a set of resources are to be shared among a group of agents. Each agent demands one resource and each resource can serve a limited number of agents. An agent cares about what resource they get as well as the externalities imposed by their mates, who share the same resource with them. Clearly, the strong notion of envy-freeness, where no agent envies another for their resource or mates, cannot always be achieved and we show that even deciding the existence of such a strongly envy-free assignment is an intractable problem. Hence, a more interesting question is whether (and in what situations) a relaxed notion of envy-freeness, the Pareto envy-freeness, can be achieved. Under this relaxed notion, an agent envies another only when they envy both the resource and the mates of the other agent. In particular, we are interested in a dorm assignment problem, where students are to be assigned to dorms with the same capacity and they have dichotomous preference over their dormmates. We show that when the capacity of each dorm is 2, a Pareto envy-free assignment always exists and we present a polynomial-time algorithm to compute such an assignment. Nevertheless, the result breaks immediately when the capacity increases to 3, in which case even Pareto envy-freeness cannot be guaranteed. In addition to the existential results, we also investigate the utility guarantees of (Pareto) envy-free assignments in our model.


title: Congruences concerning generalized central trinomial coefficients

abstract: For any $n\in\mathbb{N}=\{0,1,2,\ldots\}$ and $b,c\in\mathbb{Z}$, the generalized central trinomial coefficient $T_n(b,c)$ denotes the coefficient of $x^n$ in the expansion of $(x^2+bx+c)^n$. Let $p$ be an odd prime. In this paper, we determine the summation $\sum_{k=0}^{p-1}T_k(b,c)^2/m^k$ modulo $p^2$ for integers $m$ with certain restrictions. As applications, we confirm some conjectural congruences of Sun [Sci. China Math. 57 (2014), 1375--1400].


title: Construction of optimal spectral methods in phase retrieval

abstract: We consider the phase retrieval problem, in which the observer wishes to recover a $n$-dimensional real or complex signal $\mathbf{X}^\star$ from the (possibly noisy) observation of $|\mathbf{\Phi} \mathbf{X}^\star|$, in which $\mathbf{\Phi}$ is a matrix of size $m \times n$. We consider a \emph{high-dimensional} setting where $n,m \to \infty$ with $m/n = \mathcal{O}(1)$, and a large class of (possibly correlated) random matrices $\mathbf{\Phi}$ and observation channels. Spectral methods are a powerful tool to obtain approximate observations of the signal $\mathbf{X}^\star$ which can be then used as initialization for a subsequent algorithm, at a low computational cost. In this paper, we extend and unify previous results and approaches on spectral methods for the phase retrieval problem. More precisely, we combine the linearization of message-passing algorithms and the analysis of the \emph{Bethe Hessian}, a classical tool of statistical physics. Using this toolbox, we show how to derive optimal spectral methods for arbitrary channel noise and right-unitarily invariant matrix $\mathbf{\Phi}$, in an automated manner (i.e. with no optimization over any hyperparameter or preprocessing function).


title: GMM-Based Generative Adversarial Encoder Learning

abstract: While GAN is a powerful model for generating images, its inability to infer a latent space directly limits its use in applications requiring an encoder. Our paper presents a simple architectural setup that combines the generative capabilities of GAN with an encoder. We accomplish this by combining the encoder with the discriminator using shared weights, then training them simultaneously using a new loss term. We model the output of the encoder latent space via a GMM, which leads to both good clustering using this latent space and improved image generation by the GAN. Our framework is generic and can be easily plugged into any GAN strategy. In particular, we demonstrate it both with Vanilla GAN and Wasserstein GAN, where in both it leads to an improvement in the generated images in terms of both the IS and FID scores. Moreover, we show that our encoder learns a meaningful representation as its clustering results are competitive with the current GAN-based state-of-the-art in clustering.


title: Two dimensional non-Hermitian harmonic oscillator: coherent states

abstract: In this study, we introduce a two dimensional complex harmonic oscillator potential with space and time reflection symmetries. The corresponding time independent Schr\"odinger equation yields real eigenvalues with complex eigenfunctions. We also construct the coherent state of the system by using a superposition of 12 eigenfunctions. Using the complex correspondence principle for the probability density we investigate the possible modifications in the probability densities due to the non-Hermitian aspect of the Hamiltonian.


title: Intersymbol and Intercarrier Interference in OFDM Systems: Unified Formulation and Analysis

abstract: A new unified formulation for seven different orthogonal frequency-division multiplexing (OFDM) systems is presented. The proposed formulation relies on six parameters and encompasses conventional OFDM systems, with windowing in the transmitter and/or in the receiver, and also with a cyclic prefix (CP) or both a CP and cyclic suffix (CS). A new equivalent channel matrix that is useful for calculating both the received signal and the intersymbol and intercarrier interference power is defined and characterized. Unlike previous works, this new channel matrix formulates the channel convolution with no restrictions on the length of the channel impulse response. Moreover, it includes the overlap-and-add procedure performed in the transmitter of windowed OFDM systems. Furthermore, new theoretical expressions for the intersymbol and intercarrier interference and also for the signal-to-interference-plus-noise ratio are derived.


title: Single-Trajectory Characterization of Active Swimmers in a Flow

abstract: We develop a maximum likelihood method to infer relevant physical properties of elongated active particles. Using individual trajectories of advected swimmers as input, we are able to accurately determine their rotational diffusion coefficients and an effective measure of their aspect ratio, also providing reliable estimators for the uncertainties of such quantities. We validate our theoretical construction using numerically generated active trajectories upon no-flow, simple shear, and Poiseuille flow, with excellent results. Being designed to rely on single-particle data, our method eases applications in experimental conditions where swimmers exhibit a strong morphological diversity. We briefly discuss some of such ongoing experimental applications, specifically, in the characterization of swimming E.coli in a flow.


title: Cross-Modal Collaborative Representation Learning and a Large-Scale RGBT Benchmark for Crowd Counting

abstract: Crowd counting is a fundamental yet challenging task, which desires rich information to generate pixel-wise crowd density maps. However, most previous methods only used the limited information of RGB images and cannot well discover potential pedestrians in unconstrained scenarios. In this work, we find that incorporating optical and thermal information can greatly help to recognize pedestrians. To promote future researches in this field, we introduce a large-scale RGBT Crowd Counting (RGBT-CC) benchmark, which contains 2,030 pairs of RGB-thermal images with 138,389 annotated people. Furthermore, to facilitate the multimodal crowd counting, we propose a cross-modal collaborative representation learning framework, which consists of multiple modality-specific branches, a modality-shared branch, and an Information Aggregation-Distribution Module (IADM) to capture the complementary information of different modalities fully. Specifically, our IADM incorporates two collaborative information transfers to dynamically enhance the modality-shared and modality-specific representations with a dual information propagation mechanism. Extensive experiments conducted on the RGBT-CC benchmark demonstrate the effectiveness of our framework for RGBT crowd counting. Moreover, the proposed approach is universal for multimodal crowd counting and is also capable to achieve superior performance on the ShanghaiTechRGBD dataset. Finally, our source code and benchmark are released at {\url{http://lingboliu.com/RGBT_Crowd_Counting.html}}.


title: Mutually polar retractions on convex cones

abstract: Two retractions Q and R on closed convex cones M and respectively N of a Banach space are called mutually polar if Q+R=I and QR=RQ=0. This note investigates the existence of a pair of mutually polar retractions for given cones M and N. It is shown that if dim N=1 (or dim M=1) then the retractions are subadditive with respect to the order relation their cone range endow.


title: Spaces of Lorentzian and real stable polynomials are Euclidean balls

abstract: We prove that projective spaces of Lorentzian and real stable polynomials are homeomorphic to closed Euclidean balls. This solves a conjecture of June Huh and the author. The proof utilizes and refines a connection between the symmetric exclusion process in Interacting Particle Systems and the geometry of polynomials.


title: The Simons Observatory: Magnetic Sensitivity Measurements of Microwave SQUID Multiplexers

abstract: The Simons Observatory (SO) will be a cosmic microwave background (CMB) survey experiment with three small-aperture telescopes and one large-aperture telescope, which will observe from the Atacama Desert in Chile. In total, SO will field $\sim$70,000 transition-edge sensor (TES) bolometers in six spectral bands centered between 27 and 280 GHz in order to achieve the sensitivity necessary to measure or constrain numerous cosmological quantities. The SO Universal Focal Plane Modules (UFMs) each contain a 150 mm diameter TES detector array, horn or lenslet optical coupling, cold readout components, and magnetic shielding. SO will use a microwave SQUID multiplexing ($\mu$MUX) readout at an initial multiplexing factor of $\sim$1000; the cold (100 mK) readout components are packaged in a $\mu$MUX readout module, which is part of the UFM, and can also be characterized independently. The 100 mK stage TES bolometer arrays and microwave SQUIDs are sensitive to magnetic fields, and their measured response will vary with the degree to which they are magnetically shielded. We present measurements of the magnetic pickup of test microwave SQUID multiplexers as a study of various shielding configurations for the Simons Observatory. We discuss how these measurements motivated the material choice and design of the UFM magnetic shielding.


title: Beyond 4D Tracking: Using Cluster Shapes for Track Seeding

abstract: Tracking is one of the most time consuming aspects of event reconstruction at the Large Hadron Collider (LHC) and its high-luminosity upgrade (HL-LHC). Innovative detector technologies extend tracking to four-dimensions by including timing in the pattern recognition and parameter estimation. However, present and future hardware already have additional information that is largely unused by existing track seeding algorithms. The shape of clusters provides an additional dimension for track seeding that can significantly reduce the combinatorial challenge of track finding. We use neural networks to show that cluster shapes can reduce significantly the rate of fake combinatorical backgrounds while preserving a high efficiency. We demonstrate this using the information in cluster singlets, doublets and triplets. Numerical results are presented with simulations from the TrackML challenge.

categories: physics.ins-det hep-ex hep-ph physics.data-an stat.ML


title: Frequency Sub-Sampling of Ultrasound Non-Destructive Measurements: Acquisition, Reconstruction and Performance

abstract: In ultrasound nondestructive testing, a widespread approach is to take synthetic aperture measurements from the surface of a specimen to detect and locate defects within it. Based on these measurements, imaging is usually performed using the Synthetic Aperture Focusing Technique (SAFT). However, SAFT is sub-optimal in terms of resolution and requires oversampling in time domain to obtain a fine grid for the Delay-and-Sum (DAS). On the other hand, parametric reconstruction algorithms give better resolution, but their usage for imaging becomes computationally expensive due to the size of the parameter space and the large amount of measurement data in realistic 3-D scenarios. In the literature, the remedies to this are twofold: First, the amount of measurement data can be reduced using state of the art sub-Nyquist sampling approaches to measure Fourier coefficients instead of time domain samples. Second, parametric reconstruction algorithms mostly rely on matrix-vector operations that can be implemented efficiently by exploiting the underlying model structure. In this paper, we propose and compare different strategies to choose the Fourier coefficients to be measured. Their asymptotic performance is compared by numerically evaluating the Cram\'er-Rao-Bound for the localizability of the defect coordinates. These subsampling strategies are then combined with an $\ell_1$-minimization scheme to compute 3-D reconstructions from the low-rate measurements. Compared to conventional DAS, this allows us to formulate a fully physically motivated forward model. To enable this, the projection operations of the forward model matrix are implemented matrix-free by exploiting the underlying 2-level Toeplitz structure. Finally, we show that high resolution reconstructions from as low as a single Fourier coefficient per scan are possible based on simulated data as well as on measurements.


title: Shear Viscosity in the Strong Interaction Regime of a p-wave Superfluid Fermi Gas

abstract: The $p$-wave superfluid state is a promising spin-triplet and non $s$-wave pairing state in an ultracold Fermi gas. In this work we study the low-temperature shear viscosity of a one-component $p$-wave superfluid Fermi gas, by means of Kubo formalism. Our study is done in the strong-coupling limit where Fermi superfluid reduces into a system of composite bosons. Taking into account ${{p}_{x}}$-wave Cooper channel in the self-energy, the viscous relaxation rates are determined. The relaxation rates related to these interactions are calculated as a function of temperature. Their temperature dependence is different from the $s$-wave superfluid Fermi gas, and this is due to the anisotropic pairing interaction in the $p$-wave superfluid. Our results contribute to understand how this anisotropy affects transport properties of this unconventional superfluid Fermi gas in low temperature limit.


title: The Astrophysical Distance Scale III: Distance to the Local Group Galaxy WLM using Multi-Wavelength Observations of the Tip of the Red Giant Branch, Cepheids, and JAGB Stars

abstract: The local determination of the Hubble Constant sits at a crossroad. Current estimates of the local expansion rate of the Universe differ by about 1.7-sigma, derived from the Cepheid and TRGB based calibrations, applied to type Ia supernovae. To help elucidate possible sources of systematic error causing the tension, we show in this study the recently developed distance indicator, the J-region Asymptotic Giant Branch (JAGB) method (arXiv:2005.10792), can serve as an independent cross-check and comparison with other local distance indicators. Furthermore, we make the case that the JAGB method has substantial potential as an independent, precise and accurate calibrator of type Ia supernovae for the determination of H0. Using the Local Group galaxy, WLM we present distance comparisons between the JAGB method, a TRGB measurement at near-infrared (JHK) wavelengths, a TRGB measurement in the optical I band, and a multi-wavelength Cepheid period-luminosity relation determination. We find: $\mu_0$ (JAGB) = 24.97 +/- 0.02 (stat) +/- 0.04 (sys) mag, $\mu_0$ (TRGB NIR) = 24.98 +/- 0.04 stat) +/- 0.07 (sys) mag, $\mu_0$ (TRGB F814W) = 24.93 +/- 0.02 (stat) +/- 0.06 (sys) mag, $\mu_0$ (Cepheids) = 24.98 +/- 0.03 (stat) +/- 0.04 (sys) mag. All four methods are in good agreement, confirming the local self-consistency of the four distance scales at the 3% level, and adding confidence that the JAGB method is as accurate and as precise a distance indicator as either of the other three astrophysically-based methods.