title: Internet-Augmented Dialogue Generation

abstract: The largest store of continually updating knowledge on our planet can be accessed via internet search. In this work we study giving access to this information to conversational agents. Large language models, even though they store an impressive amount of knowledge within their weights, are known to hallucinate facts when generating dialogue (Shuster et al., 2021); moreover, those facts are frozen in time at the point of model training. In contrast, we propose an approach that learns to generate an internet search query based on the context, and then conditions on the search results to finally generate a response, a method that can employ up-to-the-minute relevant information. We train and evaluate such models on a newly collected dataset of human-human conversations whereby one of the speakers is given access to internet search during knowledgedriven discussions in order to ground their responses. We find that search-query based access of the internet in conversation provides superior performance compared to existing approaches that either use no augmentation or FAISS-based retrieval (Lewis et al., 2020).


title: Beyond Goldfish Memory: Long-Term Open-Domain Conversation

abstract: Despite recent improvements in open-domain dialogue models, state of the art models are trained and evaluated on short conversations with little context. In contrast, the long-term conversation setting has hardly been studied. In this work we collect and release a human-human dataset consisting of multiple chat sessions whereby the speaking partners learn about each other's interests and discuss the things they have learnt from past sessions. We show how existing models trained on existing datasets perform poorly in this long-term conversation setting in both automatic and human evaluations, and we study long-context models that can perform much better. In particular, we find retrieval-augmented methods and methods with an ability to summarize and recall previous conversations outperform the standard encoder-decoder architectures currently considered state of the art.


title: Pulsational instability of pre-main sequence models from accreting protostars. I. Constraining the input physics for accretion with spectroscopic parameters and stellar pulsations

abstract: Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion on protostellar seeds with observed spectroscopic parameters and stellar pulsations of young stellar objects and pre-main sequence stars. Methods: We conducted a literature search for spectroscopic samples of young stellar objects and pre-main sequence stars including all previously known pulsators. The sample size of pulsating pre-main sequence stars is increased by analysing TESS observations and presenting discoveries in CoRoT data. We employ MESA and GYRE to calculate evolutionary tracks of accreting protostellar seeds in a constant accretion scenario, the subsequent pre-main sequence evolution, and their pulsation properties. The results are then compared with observations to constrain the input physics. Results: We discuss 16 formerly unknown pulsating pre-main sequence stars and candidates that are either of SPB, $\delta$ Scuti,$\gamma$ Doradus or $\delta$ Scuti - $\gamma$ Doradus hybrid type. We find that evolutionary tracks with a mass accretion rate of $5\times10^{-6} M_\odot/{\rm yr}$ and fraction of injected accretion energy of $\beta=0.1$ provide the best results in enveloping the spectroscopic parameters of pre-main sequence stars in a constant accretion scenario. The calculated instability regions constrain the atmospheric boundary conditions to Eddington Gray atmospheres; we discuss the future potential for additional constraints by instability regions that are dependent on radial order. We present a possible candidate for pulsations in M-type young stellar objects.


title: Effect of optically thin cooling curves on condensation formation: Case study using thermal instability

abstract: Non-gravitationally induced condensations are observed in many astrophysical environments. Such structures are formed due to energy loss by optically thin radiative emission. Instead of solving the full radiative transfer equations, precomputed cooling curves are typically used in numerical simulations. In the literature, there exists a wide variety of cooling curves and they are quite often used as unquestionable ingredients. We determine the effect of the optically thin cooling curves on the formation and evolution of condensations. We perform a case study using thermal instability as a mechanism to form in situ condensations. We compare 2D numerical simulations with different cooling curves using interacting slow magnetohydrodynamic (MHD) waves as trigger for the thermal instability. Furthermore, we discuss a bootstrap measure to investigate the far non-linear regime of thermal instability. In the appendix, we include the details of all cooling curves implemented in MPI-AMRVAC and briefly discuss a hydrodynamic variant of the slow MHD waves setup for thermal instability. For all tested cooling curves, condensations are formed. However, the growth rates of the thermal instability are different. Also, the morphology of the formed condensation widely varies. We find fragmentation that is affected by the low-temperature treatment of the cooling curves. Condensations formed using cooling curves that vanish for temperatures lower than 20 000 K seem to be more stable against dynamical instabilities. The non-linear regime and fragmentation in the hydrodynamic case differ greatly from the MHD case. We advocate the use of modern cooling curves, based on accurate computations and up-to-date atomic parameters and solar abundances. Our bootstrap procedure can be used in future multi-dimensional simulations, to study fine-structure dynamics in solar prominences.


title: Revealing asymmetrical dust distribution in the inner regions of HD 141569

abstract: We obtained polarimetric differential imaging of a gas-rich debris disk around HD 141569A with SPHERE in the H-band to compare the scattering properties of the innermost ring at 44 au with former observations in total intensity with the same instrument. In polarimetric imaging, we observed that the intensity of the ring peaks in the south-east, mostly in the forward direction, whereas in total intensity imaging, the ring is detected only at the south. This noticeable characteristic suggests a non-uniform dust density in the ring. We implemented a density function varying azimuthally along the ring and generated synthetic images both in polarimetry and in total intensity, which are then compared to the actual data. We find that the dust density peaks in the south-west at an azimuthal angle of $220^{\circ} \sim 238^{\circ}$ with a rather broad width of $61^{\circ} \sim 127^{\circ}$. Although there are still uncertainties that remain in the determination of the anisotropic scattering factor, the implementation of an azimuthal density variation to fit the data proved to be robust. Upon elaborating on the origin of this dust density distribution, we conclude that it could be the result of a massive collision when we account for the effect of the high gas mass that is present in the system on the dynamics of grains. Using the outcome of this modelization, we further measured the polarized scattering phase function for the observed scattering angle between 33$^{\circ}$ and 147$^{\circ}$ as well as the spectral reflectance of the southern part of the ring between 0.98 $\mu$m and 2.1 $\mu$m. We tentatively derived the grain properties by comparing these quantities with MCFOST models and assuming Mie scattering. Our preliminary interpretation indicates a mixture of porous sub-micron sized astro-silicate and carbonaceous grains.


title: Unexplained Spectral Phenomena in the Interstellar Medium: an introduction

abstract: There exists a number of astronomical spectral phenomena that have remained unidentified after decades of extensive observations. The diffuse interstellar bands, the 220 nm feature, unidentified infrared emission bands, extended red emissions, and 21 and 30 $\mu$m emission features are seen in a wide variety of astrophysical environments. The strengths of these features suggest that they originate from chemical compounds made of common elements, possibly organic in nature. The quest to understand how such organic materials are synthesized and distributed across the Galaxy represents a major challenge to our understanding of the chemical content of the Universe.


title: Globally Convergent Multilevel Training of Deep Residual Networks

abstract: We propose a globally convergent multilevel training method for deep residual networks (ResNets). The devised method can be seen as a novel variant of the recursive multilevel trust-region (RMTR) method, which operates in hybrid (stochastic-deterministic) settings by adaptively adjusting mini-batch sizes during the training. The multilevel hierarchy and the transfer operators are constructed by exploiting a dynamical system's viewpoint, which interprets forward propagation through the ResNet as a forward Euler discretization of an initial value problem. In contrast to traditional training approaches, our novel RMTR method also incorporates curvature information on all levels of the multilevel hierarchy by means of the limited-memory SR1 method. The overall performance and the convergence properties of our multilevel training method are numerically investigated using examples from the field of classification and regression.


title: Discovery and characterization of five new eclipsing AM CVn systems

abstract: AM CVn systems are ultra-compact, helium-rich, accreting binaries with degenerate or semi-degenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 minutes. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra of two long-period systems (61.5 and 53.3 minutes) show many emission and absorption lines, indicating the presence of N, O, Na, Mg, Si, and Ca, and also the K and Zn, elements which have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is $\approx0.8$$\mathrm{M_{\odot}}$, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of $\approx$ 2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with He-star progenitors, although the observed spectral features seem to contradict this. The discovery of 5 new eclipsing AM~CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency. Based on this estimate, we expect to find another 1--4 eclipsing AM CVn systems as ZTF continues to obtain data. This will further increase our understanding of the population, but will require high precision data to better characterize these 5 systems and any new discoveries.


title: Eco-evolutionary dynamics of cooperation in the presence of policing

abstract: Ecology and evolution are inherently linked, and studying a mathematical model that considers both holds promise of insightful discoveries related to the dynamics of cooperation. In the present article, we use the prisoner's dilemma (PD) game as a basis for long-term apprehension of the essential social dilemma related to cooperation among unrelated individuals. We upgrade the contemporary PD game with an inclusion of evolution-induced act of punishment as a third competing strategy in addition to the traditional cooperators and defectors. In a population structure, the abundance of ecologically-viable free space often regulates the reproductive opportunities of the constituents. Hence, additionally, we consider the availability of free space as an ecological footprint, thus arriving at a simple eco-evolutionary model, which displays fascinating complex dynamics. As possible outcomes, we report the individual dominance of cooperators and defectors as well as a plethora of mixed states, where different strategies coexist followed by maintaining the diversity in a socio-ecological framework. These states can either be steady or oscillating, whereby oscillations are sustained by cyclic dominance among different combinations of cooperators, defectors, and punishers. We also observe a novel route to cyclic dominance where cooperators, punishers, and defectors enter a coexistence via an inverse Hopf bifurcation that is followed by an inverse period doubling route.


title: Optimal tests of the composite null hypothesis arising in mediation analysis

abstract: The indirect effect of an exposure on an outcome through an intermediate variable can be identified by a product of regression coefficients under certain causal and regression modeling assumptions. Thus, the null hypothesis of no indirect effect is a composite null hypothesis, as the null holds if either regression coefficient is zero. A consequence is that existing hypothesis tests are either severely underpowered near the origin (i.e., when both coefficients are small with respect to standard errors) or do not preserve type 1 error uniformly over the null hypothesis space. We propose hypothesis tests that (i) preserve level alpha type 1 error, (ii) meaningfully improve power when both true underlying effects are small relative to sample size, and (iii) preserve power when at least one is not. One approach gives a closed-form test that is minimax optimal with respect to local power over the alternative parameter space. Another uses sparse linear programming to produce an approximately optimal test for a Bayes risk criterion. We provide an R package that implements the minimax optimal test.


title: Real-Time Face Recognition System for Remote Employee Tracking

abstract: During the COVID-19 pandemic, most of the human-to-human interactions have been stopped. To mitigate the spread of deadly coronavirus, many offices took the initiative so that the employees can work from home. But, tracking the employees and finding out if they are really performing what they were supposed to turn out to be a serious challenge for all the companies and organizations who are facilitating "Work From Home". To deal with the challenge effectively, we came up with a solution to track the employees with face recognition. We have been testing this system experimentally for our office. To train the face recognition module, we used FaceNet with KNN using the Labeled Faces in the Wild (LFW) dataset and achieved 97.8\% accuracy. We integrated the trained model into our central system, where the employees log their time. In this paper, we discuss in brief the system we have been experimenting with and the pros and cons of the system.


title: Algebraic hyperbolicity for surfaces in smooth projective toric threefolds with Picard rank 2 and 3

abstract: Algebraic hyperbolicity serves as a bridge between differential geometry and algebraic geometry. Generally, it is difficult to show that a given projective variety is algebraically hyperbolic. However, it was established recently that a very general surface of degree at least five in projective space is algebraically hyperbolic. We are interested in generalizing the study of surfaces in projective space to surfaces in smooth projective toric threefolds with Picard rank 2 or 3. Following Kleinschmidt and Batyrev, we explore the combinatorial description of smooth projective toric threefolds with Picard rank 2 and 3. We then use Haase and Ilten's method of finding algebraically hyperbolic surfaces in toric threefolds. As a result, we determine many algebraically hyperbolic surfaces in each of these varieties.


title: Real-Time Violence Detection Using CNN-LSTM

abstract: Violence rates however have been brought down about 57% during the span of the past 4 decades yet it doesn't change the way that the demonstration of violence actually happens, unseen by the law. Violence can be mass controlled sometimes by higher authorities, however, to hold everything in line one must "Microgovern" over each movement occurring in every road of each square. To address the butterfly effects impact in our setting, I made a unique model and a theorized system to handle the issue utilizing deep learning. The model takes the input of the CCTV video feeds and after drawing inference, recognizes if a violent movement is going on. And hypothesized architecture aims towards probability-driven computation of video feeds and reduces overhead from naively computing for every CCTV video feeds.


title: A Channel Coding Benchmark for Meta-Learning

abstract: Meta-learning provides a popular and effective family of methods for data-efficient learning of new tasks. However, several important issues in meta-learning have proven hard to study thus far. For example, performance degrades in real-world settings where meta-learners must learn from a wide and potentially multi-modal distribution of training tasks; and when distribution shift exists between meta-train and meta-test task distributions. These issues are typically hard to study since the shape of task distributions, and shift between them are not straightforward to measure or control in standard benchmarks. We propose the channel coding problem as a benchmark for meta-learning. Channel coding is an important practical application where task distributions naturally arise, and fast adaptation to new tasks is practically valuable. We use our MetaCC benchmark to study several aspects of meta-learning, including the impact of task distribution breadth and shift, which can be controlled in the coding problem. Going forward, MetaCC provides a tool for the community to study the capabilities and limitations of meta-learning, and to drive research on practically robust and effective meta-learners.


title: Flow onset for a single bubble in yield-stress fluids

abstract: We use computational methods to determine the minimal yield-stress required in order to hold static a buoyant bubble in a yield-stress liquid. The static limit is governed by the bubble shape, the dimensionless surface tension ($\gamma$) and the ratio of the yield-stress to the buoyancy stress ($Y$). For a given geometry, bubbles are static for $Y > Y_c$, which we determine for a range of shapes. Given that surface tension is negligible, long prolate bubbles require larger yield-stress to hold static compared to oblate bubbles. Non-zero $\gamma$ increases $Y_c$ and for large $\gamma$ the yield-capillary number ($Y/\gamma$) determines the static boundary. In this limit, although bubble shape is important, bubble orientation is not. 2D planar and axisymmetric bubbles are studied.


title: A multiple criteria approach for ship risk classification: An alternative to the Paris MoU Ship Risk Profile

abstract: The Paris Memorandum of Understanding on Port State Control (Paris MoU) is responsible for controlling substandard shipping in European waters and, consequently, increasing the standards of safety, pollution prevention, and onboard living and working conditions. Since 2011, the Memorandum adopted a system of points, named the "Ship Risk Profile" (SRP), under which each ship is assigned a risk profile according to its score on a set of criteria. Being a multiple criteria decision aiding (MCDA) tool at its core, comprising criteria, weights, and risk categories, limited research has been performed on the SRP from an MCDA perspective. The purpose of this paper is to propose an MCDA approach for ship risk classification through the Deck of Cards Method (DCM). The DCM is particularly suitable within this context as it allows, intuitively for the decision-maker, to model preference among different criteria and among different levels on criteria scales. First, a framework is built, based on the criteria established in the SRP. Second, the DCM is used to build the MCDA model, including the definition of the criteria value functions and criteria weights. Finally, the proposed MCDA model is applied to a dataset of ships and the results are discussed. Robust results have been obtained with the proposed approach, making it a potential alternative to the current SRP.


title: Prediction of Blood Lactate Values in Critically Ill Patients: A Retrospective Multi-center Cohort Study

abstract: Purpose. Elevations in initially obtained serum lactate levels are strong predictors of mortality in critically ill patients. Identifying patients whose serum lactate levels are more likely to increase can alert physicians to intensify care and guide them in the frequency of tending the blood test. We investigate whether machine learning models can predict subsequent serum lactate changes. Methods. We investigated serum lactate change prediction using the MIMIC-III and eICU-CRD datasets in internal as well as external validation of the eICU cohort on the MIMIC-III cohort. Three subgroups were defined based on the initial lactate levels: i) normal group (<2 mmol/L), ii) mild group (2-4 mmol/L), and iii) severe group (>4 mmol/L). Outcomes were defined based on increase or decrease of serum lactate levels between the groups. We also performed sensitivity analysis by defining the outcome as lactate change of >10% and furthermore investigated the influence of the time interval between subsequent lactate measurements on predictive performance. Results. The LSTM models were able to predict deterioration of serum lactate values of MIMIC-III patients with an AUC of 0.77 (95% CI 0.762-0.771) for the normal group, 0.77 (95% CI 0.768-0.772) for the mild group, and 0.85 (95% CI 0.840-0.851) for the severe group, with a slightly lower performance in the external validation. Conclusion. The LSTM demonstrated good discrimination of patients who had deterioration in serum lactate levels. Clinical studies are needed to evaluate whether utilization of a clinical decision support tool based on these results could positively impact decision-making and patient outcomes.


title: Isometries of lattices and automorphisms of K3 surfaces

abstract: The aim of this paper is to give necessary and sufficient conditions for an integral polynomial to be the characteristic polynomial of a semi-simple isometry of some even unimodular lattice of given signature. This result has applications applications to automorphisms of K3 surfaces; in particular, we show that every Salem number of degree 4, 6, 8,1 2, 14 or 16 is the dynamical degree of an automorphism of a non-projective K3 surface.


title: Searching for Dark Matter Particles with Compton Scattering Approach

abstract: The dark matter puzzle is one of the most important fundamental physics questions in 21 century. There is no doubt that solving the puzzle will be a new milestone for human beings in the way of deeper understanding the mother nature. Here we propose to use the Shanghai laser electron gamma source (SLEGS) to search for dark matter candidates particles, including dark pseudo scalar particles, dark scalar particles, and dark photons. Our simulations show that electron facilities like SLEGS with some upgrading could be competitive platforms in searching for light dark matter particles with mass under tens of keV.


title: Updated physics performance of the ESSnuSB experiment

abstract: In this paper, we present the physics performance of the ESSnuSB experiment in the standard three flavor scenario using the updated neutrino flux calculated specifically for the ESSnuSB configuration and updated migration matrices for the far detector. Taking conservative systematic uncertainties corresponding to a normalization error of $5\%$ for signal and $10\%$ for background, we find that there is $10\sigma$ $(13\sigma)$ CP violation discovery sensitivity for the baseline option of 540 km (360 km) at $\delta_{\rm CP} = \pm 90^\circ$. The corresponding fraction of $\delta_{\rm CP}$ for which CP violation can be discovered at more than $5 \sigma$ is $70\%$. Regarding CP precision measurements, the $1\sigma$ error associated with $\delta_{\rm CP} = 0^\circ$ is around $5^\circ$ and with $\delta_{\rm CP} = -90^\circ$ is around $14^\circ$ $(7^\circ)$ for the baseline option of 540 km (360 km). For hierarchy sensitivity, one can have $3\sigma$ sensitivity for 540 km baseline except $\delta_{\rm CP} = \pm 90^\circ$ and $5\sigma$ sensitivity for 360 km baseline for all values of $\delta_{\rm CP}$. The octant of $\theta_{23}$ can be determined at $3 \sigma$ for the values of: $\theta_{23} > 51^\circ$ ($\theta_{23} < 42^\circ$ and $\theta_{23} > 49^\circ$) for baseline of 540 km (360 km). Regarding measurement precision of the atmospheric mixing parameters, the allowed values at $3 \sigma$ are: $40^\circ < \theta_{23} < 52^\circ$ ($42^\circ < \theta_{23} < 51.5^\circ$) and $2.485 \times 10^{-3}$ eV$^2 < \Delta m^2_{31} < 2.545 \times 10^{-3}$ eV$^2$ ($2.49 \times 10^{-3}$ eV$^2 < \Delta m^2_{31} < 2.54 \times 10^{-3}$ eV$^2$) for the baseline of 540 km (360 km).


title: Photon radiation effects in kinematic reconstruction of top quarks

abstract: Kinematic reconstruction of top quarks allows to define a set of kinematic observables relevant to various physics processes that involve top quarks and provides an additional handle for the suppression of background events. Radiation of photons in association with the top quarks alters the kinematics and the topology of the event, leading to visible systematic effects in measurable observables. The present study introduces an improved reconstruction of the top quark kinematics in the presence of photon radiation. The results are presented for processes with top quark pair production, as well as for singly-produced top quarks.


title: Regular Ideals of Locally-Convex Higher-Rank Graph Algebras

abstract: We give a vertex set description for basic, graded, regular ideals of locally-convex Kumjian-Pask Algebras. We also show that Condition (B) is preserved when taking the quotient by a basic, graded, regular ideal. We further show that when a locally-convex, row-finite k-graph satisfies Condition (B), all regular ideals are graded. We then show the same things hold for gauge-invariant, regular ideals in locally-convex k-graph C*-algebras.


title: Turbulence-induced clustering in compressible active fluids

abstract: We study a novel phase of active polar fluids, which is characterized by the continuous creation and destruction of dense clusters due to self-sustained turbulence. This state arises due to the interplay of the self-advection of the aligned swimmers and their defect topology. The typical cluster size is determined by the characteristic vortex size. Our results are obtained by investigating a continuum model of compressible polar active fluids, which incorporates typical experimental observations in bacterial suspensions by assuming a non-monotone dependence of speed on density.


title: A Comparison of Modern General-Purpose Visual SLAM Approaches

abstract: Advancing maturity in mobile and legged robotics technologies is changing the landscapes where robots are being deployed and found. This innovation calls for a transformation in simultaneous localization and mapping (SLAM) systems to support this new generation of service and consumer robots. No longer can traditionally robust 2D lidar systems dominate while robots are being deployed in multi-story indoor, outdoor unstructured, and urban domains with increasingly inexpensive stereo and RGB-D cameras. Visual SLAM (VSLAM) systems have been a topic of study for decades and a small number of openly available implementations have stood out: ORB-SLAM3, OpenVSLAM and RTABMap. This paper presents a comparison of these 3 modern, feature rich, and uniquely robust VSLAM techniques that have yet to be benchmarked against each other, using several different datasets spanning multiple domains negotiated by service robots. ORB-SLAM3 and OpenVSLAM each were not compared against at least one of these datasets previously in literature and we provide insight through this lens. This analysis is motivated to find general purpose, feature complete, and multi-domain VSLAM options to support a broad class of robot applications for integration into the new and improved ROS 2 Nav2 System as suitable alternatives to traditional 2D lidar solutions.


title: Coarse grid corrections in Krylov subspace evaluations of the matrix exponential

abstract: A coarse grid correction (CGC) approach is proposed to enhance the efficiency of the matrix exponential and $\varphi$ matrix function evaluations. The approach is intended for iterative methods computing the matrix-vector products with these functions. It is based on splitting the vector by which the matrix function is multiplied into a smooth part and a remaining part. The smooth part is then handled on a coarser grid, whereas the computations on the original grid are carried out with a relaxed stopping criterion tolerance. Estimates on the error are derived for the two-grid and multigrid variants of the proposed CGC algorithm. Numerical experiments demonstrate the efficiency of the algorithm, when employed in combination with Krylov subspace and Chebyshev polynomial expansion methods.


title: Efficient quantum programming using EASE gates on a trapped-ion quantum computer

abstract: Parallel operations in conventional computing have proven to be an essential tool for efficient and practical computation, and the story is not different for quantum computing. Indeed, there exists a large body of works that study advantages of parallel implementations of quantum gates for efficient quantum circuit implementations. Here, we focus on the recently invented efficient, arbitrary, simultaneously entangling (EASE) gates, available on a trapped-ion quantum computer. Leveraging its flexibility in selecting arbitrary pairs of qubits to be coupled with any degrees of entanglement, all in parallel, we show an $n$-qubit Clifford circuit can be implemented using $6\log(n)$ EASE gates, an $n$-qubit multiply-controlled NOT gate can be implemented using $3n/2$ EASE gates, and an $n$-qubit permutation can be implemented using six EASE gates. We discuss their implications to near-term quantum chemistry simulations and the state of the art pattern matching algorithm. Given Clifford + multiply-controlled NOT gates form a universal gate set for quantum computing, our results imply efficient quantum computation by EASE gates, in general.


title: Hyperuniform scalar random fields for lensless, multispectral imaging systems

abstract: We propose a novel framework for the systematic design of lensless imaging systems based on the hyperuniform random field solutions of nonlinear reaction-diffusion equations from pattern formation theory. Specifically, we introduce a new class of imaging point-spread-functions (PSFs) with enhanced isotropic behavior and controllable sparsity. We investigate the PSFs and the modulated transfer functions (MTFs) for a number of nonlinear models and demonstrate that two-phase isotropic random fields with hyperuniform disorder are ideally suited to construct imaging PSFs with improved performances compared to PSFs based on the Perlin noise. Additionally, we introduce a phase retrieval algorithm based on the non-paraxial Rayleigh-Sommerfeld diffraction theory and introduce diffractive phase plates with PSFs designed from hyperuniform random fields, called hyperuniform phase plates (HPPs). Finally, using high-fidelity object reconstruction, we demonstrate improved image quality using engineered HPPs across the visible range. The proposed framework is suitable for high-performance lensless imaging systems for on-chip microscopy and spectroscopy applications.


title: On Bayesian data assimilation for PDEs with ill-posed forward problems

abstract: We study Bayesian data assimilation (filtering) for time-evolution PDEs, for which the underlying forward problem may be very unstable or ill-posed. Such PDEs, which include the Navier-Stokes equations of fluid dynamics, are characterized by a high sensitivity of solutions to perturbations of the initial data, a lack of rigorous global well-posedness results as well as possible non-convergence of numerical approximations. Under very mild and readily verifiable general hypotheses on the forward solution operator of such PDEs, we prove that the posterior measure expressing the solution of the Bayesian filtering problem is stable with respect to perturbations of the noisy measurements, and we provide quantitative estimates on the convergence of approximate Bayesian filtering distributions computed from numerical approximations. For the Navier-Stokes equations, our results imply uniform stability of the filtering problem even at arbitrarily small viscosity, when the underlying forward problem may become ill-posed, as well as the compactness of numerical approximants in a suitable metric on time-parametrized probability measures.


title: Ranked Sparsity: A Cogent Regularization Framework for Selecting and Estimating Feature Interactions and Polynomials

abstract: We explore and illustrate the concept of ranked sparsity, a phenomenon that often occurs naturally in modeling applications when an expected disparity exists in the quality of information between different feature sets. Its presence can cause traditional and modern model selection methods to fail because such procedures commonly presume that each potential parameter is equally worthy of entering into the final model - we call this presumption "covariate equipoise". However, this presumption does not always hold, especially in the presence of derived variables. For instance, when all possible interactions are considered as candidate predictors, the premise of covariate equipoise will often produce over-specified and opaque models. The sheer number of additional candidate variables grossly inflates the number of false discoveries in the interactions, resulting in unnecessarily complex and difficult-to-interpret models with many (truly spurious) interactions. We suggest a modeling strategy that requires a stronger level of evidence in order to allow certain variables (e.g. interactions) to be selected in the final model. This ranked sparsity paradigm can be implemented with the sparsity-ranked lasso (SRL). We compare the performance of SRL relative to competing methods in a series of simulation studies, showing that the SRL is a very attractive method because it is fast, accurate, and produces more transparent models (with fewer false interactions). We illustrate its utility in an application to predict the survival of lung cancer patients using a set of gene expression measurements and clinical covariates, searching in particular for gene-environment interactions.


title: Resource Allocation in a Quantum Key Distribution Network with LEO and GEO trusted-repeaters

abstract: Quantum Key Distribution~(QKD) is a technology that enables the exchange of private encryption keys between two legitimate parties, using protocols that involve quantum mechanics principles. The rate at which secret keys can be exchanged depends on the attenuation that is experienced. Therefore, it is more convenient to replace many terrestrial fiber segments (and repeaters) by just few optical satellite links that would enable flexible global coverage. Then, the satellite nodes can take the role of trusted-relays, forwarding the secret keys from source to destination. However, since the rate at which secret keys can be generated in each quantum link is limited, it is very important to select the intermediate satellite nodes to inter-connect ground stations efficiently. This paper studies the most convenient allocation of resources in a QKD network that combines complementary connectivity services of GEO and LEO satellites. The aim of the centralized routing algorithm is to select the most convenient trusted-relays to forward the secret keys between pairs of ground stations, verifying the constraints that satellite-to-ground and inter-satellite quantum channels have.


title: Depth Estimation from Monocular Images and Sparse radar using Deep Ordinal Regression Network

abstract: We integrate sparse radar data into a monocular depth estimation model and introduce a novel preprocessing method for reducing the sparseness and limited field of view provided by radar. We explore the intrinsic error of different radar modalities and show our proposed method results in more data points with reduced error. We further propose a novel method for estimating dense depth maps from monocular 2D images and sparse radar measurements using deep learning based on the deep ordinal regression network by Fu et al. Radar data are integrated by first converting the sparse 2D points to a height-extended 3D measurement and then including it into the network using a late fusion approach. Experiments are conducted on the nuScenes dataset. Our experiments demonstrate state-of-the-art performance in both day and night scenes.


title: Coarse-grained Stochastic Model of Myosin-Driven Vesicles into Dendritic Spines

abstract: We study the dynamics of membrane vesicle motor transport into dendritic spines, which are bulbous intracellular compartments in neurons that play a key role in transmitting signals between neurons. We consider the stochastic analog of the vesicle transport model in [Park and Fai, The Dynamics of Vesicles Driven Into Closed Constrictions by Molecular Motors. Bull. Math. Biol. 82, 141 (2020)]. The stochastic version, which may be considered as an agent-based model, relies mostly on the action of individual myosin motors to produce vesicle motion. To aid in our analysis, we coarse-grain this agent-based model using a master equation combined with a partial differential equation describing the probability of local motor positions. We confirm through convergence studies that the coarse-graining captures the essential features of bistability in velocity (observed in experiments) and waiting-time distributions to switch between steady-state velocities. Interestingly, these results allow us to reformulate the translocation problem in terms of the mean first passage time for a run-and-tumble particle moving on a finite domain with absorbing boundaries at the two ends. We conclude by presenting numerical and analytical calculations of vesicle translocation.


title: Speeding up reactive transport simulations in cement systems by surrogate geochemical modeling: deep neural networks and k-nearest neighbors

abstract: We accelerate reactive transport (RT) simulation by replacing the geochemical solver in the RT code by a surrogate model or emulator, considering either a trained deep neural network (DNN) or a k-nearest neighbor (kNN) regressor. We focus on 2D leaching of hardened cement paste under diffusive or advective-dispersive transport conditions, a solid solution representation of the calcium silicate hydrates and either 4 or 7 chemical components, and use the HPx (coupled Hydrus-PHREEQC model) reactive transport code as baseline. We find that after training, both our DNN-based and kNN-based codes, HPx-DNN and HPx-kNN, can make satisfactorily to very accurate predictions while providing either a 3 to 9 speedup factor compared to HPx with parallelized geochemical calculations over 4 cores. Benchmarking against single-threaded HPx, these speedup factors become 8 to 33. Overall, HPx-DNN are HPx-kNN are found to achieve a close to optimal speedup when DNN regression and kNN search are performed on a GPU. Importantly, for the more complex 7-components cement system, no emulator that is globally accurate over the full space of possible geochemical conditions could be devised. Instead we therefore build "local" emulators that are only valid over a relevant fraction of the input parameter space. This space is identified by running a coarse and thus computationally cheap full RT simulation, and subsequently explored by kernel density sampling. Future work will focus on improving accuracy for this type of cement systems.


title: Partially Observable Markov Decision Processes (POMDPs) and Robotics

abstract: Planning under uncertainty is critical to robotics. The Partially Observable Markov Decision Process (POMDP) is a mathematical framework for such planning problems. It is powerful due to its careful quantification of the non-deterministic effects of actions and partial observability of the states. But precisely because of this, POMDP is notorious for its high computational complexity and deemed impractical for robotics. However, since early 2000, POMDPs solving capabilities have advanced tremendously, thanks to sampling-based approximate solvers. Although these solvers do not generate the optimal solution, they can compute good POMDP solutions that significantly improve the robustness of robotics systems within reasonable computational resources, thereby making POMDPs practical for many realistic robotics problems. This paper presents a review of POMDPs, emphasizing computational issues that have hindered its practicality in robotics and ideas in sampling-based solvers that have alleviated such difficulties, together with lessons learned from applying POMDPs to physical robots.


title: Do quasar X-ray and UV flux measurements provide a useful test of cosmological models?

abstract: The recent compilation of quasar (QSO) X-ray and UV flux measurements include QSOs that appear to not be standardizable via the X-ray luminosity and UV luminosity ($L_X-L_{UV}$) relation and so should not be used to constrain cosmological model parameters. Here we show that the largest of seven sub-samples in this compilation, the SDSS-4XMM QSOs that contribute about 2/3 of the total QSOs, have $L_X-L_{UV}$ relations that depend on the cosmological model assumed and also on redshift, and is the main cause of the similar problem discovered earlier for the full QSO compilation. The second and third biggest sub-samples, the SDSS-Chandra and XXL QSOs that together contribute about 30% of the total QSOs, appear standardizable, but provide only weak constraints on cosmological parameters that are not inconsistent with the standard spatially-flat $\Lambda$CDM model or with constraints from better-established cosmological probes.


title: "Fast" and Furious focal-plane wavefront sensing at W. M. Keck Observatory

abstract: High quality, repeatable point-spread functions are important for science cases like direct exoplanet imaging, high-precision astrometry, and high-resolution spectroscopy of exoplanets. For such demanding applications, the initial on-sky point-spread function delivered by the adaptive optics system can require further optimization to correct unsensed static aberrations and calibration biases. We investigated using the Fast and Furious focal-plane wavefront sensing algorithm as a potential solution. This algorithm uses a simple model of the optical system and focal plane information to measure and correct the point-spread function phase, without using defocused images, meaning it can run concurrently with science. On-sky testing demonstrated significantly improved PSF quality in only a few iterations, with both narrow and broadband filters. These results suggest this algorithm is a useful path forward for creating and maintaining high-quality, repeatable on-sky adaptive optics point-spread functions.


title: Optimal-Design Domain-Adaptation for Exposure Prediction in Two-Stage Epidemiological Studies

abstract: In the first stage of a two-stage study, the researcher uses a statistical model to impute the unobserved exposures. In the second stage, imputed exposures serve as covariates in epidemiological models. Imputation error in the first stage operate as measurement errors in the second stage, and thus bias exposure effect estimates. This study aims to improve the estimation of exposure effects by sharing information between the first and second stage. At the heart of our estimator is the observation that not all second-stage observations are equally important to impute. We thus borrow ideas from the optimal-experimental-design theory, to identify individuals of higher importance. We then improve the imputation of these individuals using ideas from the machine-learning literature of domain-adaptation. Our simulations confirm that the exposure effect estimates are more accurate than the current best practice. An empirical demonstration yields smaller estimates of PM effect on hyperglycemia risk, with tighter confidence bands. Sharing information between environmental scientist and epidemiologist improves health effect estimates. Our estimator is a principled approach for harnessing this information exchange, and may be applied to any two stage study.


title: Measuring inter-cluster similarities with Alpha Shape TRIangulation in loCal Subspaces (ASTRICS) facilitates visualization and clustering of high-dimensional data

abstract: Clustering and visualizing high-dimensional (HD) data are important tasks in a variety of fields. For example, in bioinformatics, they are crucial for analyses of single-cell data such as mass cytometry (CyTOF) data. Some of the most effective algorithms for clustering HD data are based on representing the data by nodes in a graph, with edges connecting neighbouring nodes according to some measure of similarity or distance. However, users of graph-based algorithms are typically faced with the critical but challenging task of choosing the value of an input parameter that sets the size of neighbourhoods in the graph, e.g. the number of nearest neighbours to which to connect each node or a threshold distance for connecting nodes. The burden on the user could be alleviated by a measure of inter-node similarity that can have value 0 for dissimilar nodes without requiring any user-defined parameters or thresholds. This would determine the neighbourhoods automatically while still yielding a sparse graph. To this end, I propose a new method called ASTRICS to measure similarity between clusters of HD data points based on local dimensionality reduction and triangulation of critical alpha shapes. I show that my ASTRICS similarity measure can facilitate both clustering and visualization of HD data by using it in Stage 2 of a three-stage pipeline: Stage 1 = perform an initial clustering of the data by any method; Stage 2 = let graph nodes represent initial clusters instead of individual data points and use ASTRICS to automatically define edges between nodes; Stage 3 = use the graph for further clustering and visualization. This trades the critical task of choosing a graph neighbourhood size for the easier task of essentially choosing a resolution at which to view the data. The graph and consequently downstream clustering and visualization are then automatically adapted to the chosen resolution.


title: CLEDGE: A Hybrid Cloud-Edge Computing Framework over Information Centric Networking

abstract: In today's era of Internet of Things (IoT), where massive amounts of data are produced by IoT and other devices, edge computing has emerged as a prominent paradigm for low-latency data processing. However, applications may have diverse latency requirements: certain latency-sensitive processing operations may need to be performed at the edge, while delay-tolerant operations can be performed on the cloud, without occupying the potentially limited edge computing resources. To achieve that, we envision an environment where computing resources are distributed across edge and cloud offerings. In this paper, we present the design of CLEDGE (CLoud + EDGE), an information-centric hybrid cloud-edge framework, aiming to maximize the on-time completion of computational tasks offloaded by applications with diverse latency requirements. The design of CLEDGE is motivated by the networking challenges that mixed reality researchers face. Our evaluation demonstrates that CLEDGE can complete on-time more than 90% of offloaded tasks with modest overheads.


title: Quantifying the economic response to COVID-19 mitigations and death rates via forecasting Purchasing Managers' Indices using Generalised Network Autoregressive models with exogenous variables

abstract: Knowledge of the current state of economies, how they respond to COVID-19 mitigations and indicators, and what the future might hold for them is important. We use recently-developed generalised network autoregressive (GNAR) models, using trade-determined networks, to model and forecast the Purchasing Managers' Indices for a number of countries. We use networks that link countries where the links themselves, or their weights, are determined by the degree of export trade between the countries. We extend these models to include node-specific time series exogenous variables (GNARX models), using this to incorporate COVID-19 mitigation stringency indices and COVID-19 death rates into our analysis. The highly parsimonious GNAR models considerably outperform vector autoregressive models in terms of mean-squared forecasting error and our GNARX models themselves outperform GNAR ones. Further mixed frequency modelling predicts the extent to which that the UK economy will be affected by harsher, weaker or no interventions.


title: CMU-GPR Dataset: Ground Penetrating Radar Dataset for Robot Localization and Mapping

abstract: There has been exciting recent progress in using radar as a sensor for robot navigation due to its increased robustness to varying environmental conditions. However, within these different radar perception systems, ground penetrating radar (GPR) remains under-explored. By measuring structures beneath the ground, GPR can provide stable features that are less variant to ambient weather, scene, and lighting changes, making it a compelling choice for long-term spatio-temporal mapping. In this work, we present the CMU-GPR dataset--an open-source ground penetrating radar dataset for research in subsurface-aided perception for robot navigation. In total, the dataset contains 15 distinct trajectory sequences in 3 GPS-denied, indoor environments. Measurements from a GPR, wheel encoder, RGB camera, and inertial measurement unit were collected with ground truth positions from a robotic total station. In addition to the dataset, we also provide utility code to convert raw GPR data into processed images. This paper describes our recording platform, the data format, utility scripts, and proposed methods for using this data.


title: An epsilon-near-zero-based Dallenbach absorber

abstract: We report a theoretical analysis of total absorption conditions in a structure consisting of a lossy coating layer on top of a specular metal substrate, which is known as a Dallenbach absorber. All possible combinations of complex refractive indices of the coating material and thicknesses of the coating layer providing complete absorption are described for the case of normal incidence of a plane wave. The work of coatings based on epsilon-near-zero (ENZ) materials is analyzed in detail. Simple analytical relations are obtained for the conditions of total absorption. It is shown that the range of angles where strong absorption occurs is limited by the phenomenon of total external reflection. The characteristics of coatings based on indium-tin oxide, which is an example of an ENZ material, are analyzed.


title: Multi-Level Contrastive Learning for Few-Shot Problems

abstract: Contrastive learning is a discriminative approach that aims at grouping similar samples closer and diverse samples far from each other. It it an efficient technique to train an encoder generating distinguishable and informative representations, and it may even increase the encoder's transferability. Most current applications of contrastive learning benefit only a single representation from the last layer of an encoder.In this paper, we propose a multi-level contrasitive learning approach which applies contrastive losses at different layers of an encoder to learn multiple representations from the encoder. Afterward, an ensemble can be constructed to take advantage of the multiple representations for the downstream tasks. We evaluated the proposed method on few-shot learning problems and conducted experiments using the mini-ImageNet and the tiered-ImageNet datasets. Our model achieved the new state-of-the-art results for both datasets, comparing to previous regular, ensemble, and contrastive learing (single-level) based approaches.


title: Two photon exchange contribution to small angle $\mu p$ and $ep$ elastic scattering

abstract: We consider the role of the two-photon corrections to the small angle $\mu p$ and $ep$ elastic scattering and the expected $\sigma(\mu^- p)/\sigma(\mu^+ p)$ ratio.


title: Self-Supervised Contrastive Learning with Adversarial Perturbations for Defending Word Substitution-based Attacks

abstract: In this paper, we present an approach to improve the robustness of BERT language models against word substitution-based adversarial attacks by leveraging adversarial perturbations for self-supervised contrastive learning. We create a word-level adversarial attack generating hard positives on-the-fly as adversarial examples during contrastive learning. In contrast to previous works, our method improves model robustness without using any labeled data. Experimental results show that our method improves robustness of BERT against four different word substitution-based adversarial attacks, and combining our method with adversarial training gives higher robustness than adversarial training alone. As our method improves the robustness of BERT purely with unlabeled data, it opens up the possibility of using large text datasets to train robust language models against word substitution-based adversarial attacks.


title: Photonic analogue of Mollow triplet with on-chip photon pair generation in dressed modes

abstract: Making analogy with atomic physics is a powerful tool for photonic technology, witnessed by the recent development in topological photonics and non-Hermitian photonics based on parity-time symmetry. Mollow triplet is a prominent atomic effect with both fundamental and technological importance. Here we demonstrate the analogue of Mollow triplet with quantum photonic systems. Photonic entanglement is generated with spontaneous nonlinear processes in dressed photonic modes, which are introduced through coherent multimode coupling. We further demonstrate the flexibility of the photonic system to realize different configurations of dressed states, leading to the controlled modification of Mollow triplet. Our work would enable the investigation of complex atomic processes and the realization of unique quantum functionalities based on photonic systems.


title: Understanding Heating in Active Region Cores through Machine Learning II. Classifying Observations

abstract: Constraining the frequency of energy deposition in magnetically-closed active region cores requires sophisticated hydrodynamic simulations of the coronal plasma and detailed forward modeling of the optically-thin line-of-sight integrated emission. However, understanding which set of model inputs best matches a set of observations is complicated by the need for any proposed heating model to simultaneously satisfy multiple observable constraints. In this paper, we train a random forest classification model on a set of forward-modeled observable quantities, namely the emission measure slope, the peak temperature of the emission measure distribution, and the time lag and maximum cross-correlation between multiple pairs of AIA channels. We then use our trained model to classify the heating frequency in every pixel of active region NOAA 1158 using the observed emission measure slopes, peak temperatures, time lags, and maximum cross-correlations and are able to map the heating frequency across the entire active region. We find that high-frequency heating dominates in the inner core of the active region while intermediate frequency dominates closer to the periphery of the active region. Additionally, we assess the importance of each observed quantity in our trained classification model and find that the emission measure slope is the dominant feature in deciding with which heating frequency a given pixel is most consistent. The technique presented here offers a very promising and widely applicable method for assessing observations in terms of detailed forward models given an arbitrary number of observable constraints.


title: Nuclear Data Sheets for A=201

abstract: Evaluated nuclear structure and decay data for all nuclei with mass number A=201 (201Os, 201Ir, 201Pt, 201Au, 201Hg, 201Tl, 201Pb, 201Bi, 201Po, 201At, 201Rn, 201Fr, 201Ra) are presented. All available experimental data are compiled and evaluated, and best values for level and gamma-ray energies, quantum numbers, lifetimes, gamma-ray intensities and transition probabilities, as well as other nuclear properties, are recommended. Inconsistencies and discrepancies that exist in the literature are discussed. A number of computer codes (https://wwwnds. iaea.org/public/ensdf pgm/index.htm) developed by members of the NSDD network were used during the evaluation process. For example, the reported absolute gamma-ray emission probabilities and their uncertainties in various decay data sets were determined using the GABS code. The gamma-ray transition probabilities were determined using the RULER code and the corresponding uncertainties were determined using a Monte-Carlo approach. This work supersedes the earlier evaluation by F.G. Kondev (2007Ko06), published in Nuclear Data Sheets 108, 365 (2007).


title: Extremal Steering Assemblages

abstract: Non-local correlations between a fully characterised quantum system and an untrusted black box device are described by an assemblage of conditional quantum states. These assemblages form a convex set, whose extremal points are relevant in many operational contexts. We give necessary and sufficient conditions for an assemblage to be extremal using linear independence conditions, and an algorithm to decompose a generic assemblage into extremal points. A Matlab implementation of this algorithm is provided in the supplementary material.


title: Streaming Velocity Effects on the Post-reionization 21 cm Baryon Acoustic Oscillation Signal

abstract: The relative velocity between baryons and dark matter in the early Universe can suppress the formation of small-scale baryonic structure and leave an imprint on the baryon acoustic oscillation (BAO) scale at low redshifts after reionization. This "streaming velocity" affects the post-reionization gas distribution by directly reducing the abundance of pre-existing mini-halos ($\lesssim 10^7 M_{\bigodot}$) that could be destroyed by reionization and indirectly modulating reionization history via photoionization within these mini-halos. In this work, we investigate the effect of streaming velocity on the BAO feature in HI 21 cm intensity mapping after reionization, with a focus on redshifts $3.5\lesssim z\lesssim5.5$. We build a spatially modulated halo model that includes the dependence of the filtering mass on the local reionization redshift and thermal history of the intergalactic gas. In our fiducial model, we find isotropic streaming velocity bias coefficients $b_v$ ranging from $-0.0043$ at $z=3.5$ to $-0.0273$ at $z=5.5$, which indicates that the BAO scale is stretched (i.e., the peaks shift to lower $k$). In particular, streaming velocity shifts the transverse BAO scale between 0.121% ($z=3.5$) and 0.35% ($z=5.5$) and shifts the radial BAO scale between 0.167% ($z=3.5$) and 0.505% ($z=5.5$). These shifts exceed the projected error bars from the more ambitious proposed hemispherical-scale surveys in HI (0.13% at $1\sigma$ per $\Delta z = 0.5$ bin).


title: Magnetic Control of Soft Chiral Phonons in PbTe

abstract: PbTe crystals have a soft transverse optical phonon mode in the terahertz frequency range, which is known to efficiently decay into heat-carrying acoustic phonons, resulting in anomalously low thermal conductivity. Here, we studied this phonon via polarization-dependent terahertz spectroscopy. We observed softening of this mode with decreasing temperature, indicative of incipient ferroelectricity, which we explain through a model including strong anharmonicity with a quartic displacement term. In magnetic fields up to 25T, the phonon mode split into two modes with opposite handedness, exhibiting circular dichroism. Their frequencies displayed Zeeman splitting together with an overall diamagnetic shift with increasing magnetic field. Using a group-theoretical approach, we demonstrate that these observations are results of magnetic field-induced morphic changes in the crystal symmetries through the Lorentz force exerted on the lattice ions. This study thus reveals a novel process of controlling phonon properties in a soft ionic lattice by a strong magnetic field.


title: Algorithmic insights on continual learning from fruit flies

abstract: Continual learning in computational systems is challenging due to catastrophic forgetting. We discovered a two layer neural circuit in the fruit fly olfactory system that addresses this challenge by uniquely combining sparse coding and associative learning. In the first layer, odors are encoded using sparse, high dimensional representations, which reduces memory interference by activating non overlapping populations of neurons for different odors. In the second layer, only the synapses between odor activated neurons and the output neuron associated with the odor are modified during learning; the rest of the weights are frozen to prevent unrelated memories from being overwritten. We show empirically and analytically that this simple and lightweight algorithm significantly boosts continual learning performance. The fly associative learning algorithm is strikingly similar to the classic perceptron learning algorithm, albeit two modifications, which we show are critical for reducing catastrophic forgetting. Overall, fruit flies evolved an efficient lifelong learning algorithm, and circuit mechanisms from neuroscience can be translated to improve machine computation.


title: Adversarial Attack for Uncertainty Estimation: Identifying Critical Regions in Neural Networks

abstract: We propose a novel method to capture data points near decision boundary in neural network that are often referred to a specific type of uncertainty. In our approach, we sought to perform uncertainty estimation based on the idea of adversarial attack method. In this paper, uncertainty estimates are derived from the input perturbations, unlike previous studies that provide perturbations on the model's parameters as in Bayesian approach. We are able to produce uncertainty with couple of perturbations on the inputs. Interestingly, we apply the proposed method to datasets derived from blockchain. We compare the performance of model uncertainty with the most recent uncertainty methods. We show that the proposed method has revealed a significant outperformance over other methods and provided less risk to capture model uncertainty in machine learning.


title: Adiabatic soliton management: Controlling solitary wave motion while keeping the wave envelope unchanged

abstract: We describe how to control the motion -- both speed and propagation direction -- of a nonlinear traveling wave in real time via soliton management with time-varying dispersion/diffusion and loss/gain terms. When carried out subject to certain parameter constraints we derive, the approach allows for real-time manipulation of wave motion without compromising the structure or stability of the wave, maintaining strict adiabaticity in time. This makes our approach a promising tool for physical systems which admit a degree of control through adjustable dispersion and loss/gain terms, including problems arising in nonlinear optics, atomic physics, Bose-Einstein condensates, and thermally driven chemical reactions, where one is interested in modifying the position or speed of a wave while still maintaining adiabaticity. We illustrate both the utility and simplicity of the method through several canonical examples, demonstrating how the motion of KdV solitons, Fisher-KPP solitary wavefronts, BBM solitons, nonlinear Klein-Gordon equations, and interacting nonlinear waves in a vector NLS system can be controlled in real time.


title: Mechanism of primordial black holes production and secondary gravitational waves in $\alpha$-attractor Galileon inflationary scenario

abstract: We study the process of the Primordial Black Holes (PBHs) production in the novel framework, namely $\alpha$-attractor Galileon inflation (G-inflation) model. In our framework, we take the Galileon function as $G(\phi)=G_{I}(\phi)\left(1+G_{II}(\phi)\right)$, where the part $G_{I}(\phi)$ is motivated from the $\alpha$-attractor inflationary scenario in its original non-canonical frame, and it ensures for the model to be consistent with the Planck 2018 observations at the CMB scales. The part $G_{II}(\phi)$ is invoked to enhance the curvature perturbations at some smaller scales which in turn gives rise to PBHs formation. By fine-tuning of the model parameters, we find three parameter sets which successfully produce a sufficiently large peak in the curvature power spectrum. We show that these parameter sets produce PBHs with masses ${\cal O}(10)M_\odot$, ${\cal O}(10^{-5})M_\odot$, and ${\cal O}(10^{-13})M_\odot$ which can explain the LIGO events, the ultrashort-timescale microlensing events in OGLE data, and around $0.98\%$ of the current Dark Matter (DM) content of the universe, respectively. Additionally, we study the secondary Gravitational Waves (GWs) in our setup and show that our model anticipates the peak of their present fractional energy density as $\Omega_{GW0} \sim 10^{-8}$ for all the three parameter sets, but at different frequencies. These predictions can be located well inside the sensitivity region of some GWs detectors, and therefore the compatibility of our model can be assessed in light of the future data. We further estimate the tilts of the included GWs spectrum in the different ranges of frequency, and confirm that spectrum follows the power-law relation $\Omega_{GW0}\sim f^{n}$ in those frequency bands.


title: Double Categories of Relations

abstract: A double category of relations is essentially a cartesian equipment with strong, discrete and functorial tabulators and for which certain local products satisfy a Frobenius Law. A double category of relations is equivalent to a double category whose proarrows are relations on some ordinary category admitting a proper and stable factorization system. This characterization is based closely on the recent characterization of double categories of spans due to Aleiferi. The overall development can be viewed as a double-categorical version of that of the notion of a "tabular allegory" or that of a "functionally complete bicategory of relations."


title: Hybrid Analog and Digital Beamforming Design for Channel Estimation in Correlated Massive MIMO Systems

abstract: In this paper, we study the channel estimation problem in correlated massive multiple-input-multiple-output (MIMO) systems with a reduced number of radio-frequency (RF) chains. Importantly, other than the knowledge of channel correlation matrices, we make no assumption as to the structure of the channel. To address the limitation in the number of RF chains, we employ hybrid beamforming, comprising a low dimensional digital beamformer followed by an analog beamformer implemented using phase shifters. Since there is no dedicated RF chain per transmitter/receiver antenna, the conventional channel estimation techniques for fully-digital systems are impractical. By exploiting the fact that the channel entries are uncorrelated in its eigen-domain, we seek to estimate the channel entries in this domain. Due to the limited number of RF chains, channel estimation is typically performed in multiple time slots. Under a total energy budget, we aim to design the hybrid transmit beamformer (precoder) and the receive beamformer (combiner) in each training time slot, in order to estimate the channel using the minimum mean squared error criterion. To this end, we choose the precoder and combiner in each time slot such that they are aligned to transmitter and receiver eigen-directions, respectively. Further, we derive a water-filling-type expression for the optimal energy allocation at each time slot. This expression illustrates that, with a low training energy budget, only significant components of the channel need to be estimated. In contrast, with a large training energy budget, the energy is almost equally distributed among all eigen-directions. Simulation results show that the proposed channel estimation scheme can efficiently estimate correlated massive MIMO channels within a few training time slots.


title: Correlation detection in trees for planted graph alignment

abstract: Motivated by alignment of correlated sparse random graphs, we introduce a hypothesis testing problem of deciding whether or not two random trees are correlated. We obtain sufficient conditions under which this testing is impossible or feasible. We propose MPAlign, a message-passing algorithm for graph alignment inspired by the tree correlation detection problem. We prove MPAlign to succeed in polynomial time at partial alignment whenever tree detection is feasible. As a result our analysis of tree detection reveals new ranges of parameters for which partial alignment of sparse random graphs is feasible in polynomial time. We then conjecture that graph alignment is not feasible in polynomial time when the associated tree detection problem is impossible. If true, this conjecture together with our sufficient conditions on tree detection impossibility would imply the existence of a hard phase for graph alignment, i.e. a parameter range where alignment cannot be done in polynomial time even though it is known to be feasible in non-polynomial time.

categories: cs.DS cs.LG math.PR math.ST stat.ML stat.TH


title: A DPG method for shallow shells

abstract: We develop and analyze a discontinuous Petrov--Galerkin method with optimal test functions (DPG method) for a shallow shell model of Koiter type. It is based on a uniformly stable ultraweak formulation and thus converges robustly quasi-uniformly. Numerical experiments for various cases, including the Scordelis--Lo cylindrical roof, elliptic and hyperbolic geometries, illustrate its performance. The built-in DPG error estimator gives rise to adaptive mesh refinements that are capable to resolve boundary and interior layers. The membrane locking is dealt with by raising the polynomial degree only of the tangential displacement trace variable.


title: Deterministic and Las Vegas Algorithms for Sparse Nonnegative Convolution

abstract: Computing the convolution $A\star B$ of two length-$n$ integer vectors $A,B$ is a core problem in several disciplines. It frequently comes up in algorithms for Knapsack, $k$-SUM, All-Pairs Shortest Paths, and string pattern matching problems. For these applications it typically suffices to compute convolutions of nonnegative vectors. This problem can be classically solved in time $O(n\log n)$ using the Fast Fourier Transform. However, often the involved vectors are sparse and hence one could hope for output-sensitive algorithms to compute nonnegative convolutions. This question was raised by Muthukrishnan and solved by Cole and Hariharan (STOC '02) by a randomized algorithm running in near-linear time in the (unknown) output-size $t$. Chan and Lewenstein (STOC '15) presented a deterministic algorithm with a $2^{O(\sqrt{\log t\cdot\log\log n})}$ overhead in running time and the additional assumption that a small superset of the output is given; this assumption was later removed by Bringmann and Nakos (ICALP '21). In this paper we present the first deterministic near-linear-time algorithm for computing sparse nonnegative convolutions. This immediately gives improved deterministic algorithms for the state-of-the-art of output-sensitive Subset Sum, block-mass pattern matching, $N$-fold Boolean convolution, and others, matching up to log-factors the fastest known randomized algorithms for these problems. Our algorithm is a blend of algebraic and combinatorial ideas and techniques. Additionally, we provide two fast Las Vegas algorithms for computing sparse nonnegative convolutions. In particular, we present a simple $O(t\log^2t)$ time algorithm, which is an accessible alternative to Cole and Hariharan's algorithm. We further refine this new algorithm to run in Las Vegas time $O(t\log t\cdot\log\log t)$, matching the running time of the dense case apart from the $\log\log t$ factor.


title: Multiple recurrence and popular differences for polynomial patterns in rings of integers

abstract: We demonstrate that the phenomenon of popular differences (aka the phenomenon of large intersections) holds for natural families of polynomial patterns in rings of integers of number fields. If $K$ is a number field with ring of integers $\mathcal{O}_K$ and $E \subseteq \mathcal{O}_K$ has positive upper Banach density $d^*(E) = \delta > 0$, we show, inter alia: 1. If $p(x) \in K[x]$ is an intersective $\mathcal{O}_K$-valued polynomial and $r, s \in \mathcal{O}_K$ are distinct and nonzero, then for any $\varepsilon > 0$, the set of $n \in \mathcal{O}_K$ such that \[ d^* \left( \{ x \in \mathcal{O}_K : \{x, x + rp(n), x + sp(n)\} \subseteq E \} \right) > \delta^3 - \varepsilon. \] is syndetic. Moreover, if $\frac{s}{r} \in \mathbb{Q}$, then there are syndetically many $n \in \mathcal{O}_K$ such that \[ d^* \left( \{ x \in \mathcal{O}_K : \{x, x + rp(n), x + sp(n), x + (r+s)p(n)\} \subseteq E \} \right) > \delta^4 - \varepsilon. \] 2. If $\{p_1, \dots, p_k\} \subseteq K[x]$ is a jointly intersective family of linearly independent $\mathcal{O}_K$-valued polynomials, then the set of $n \in \mathcal{O}_K$ such that \[ d^* \left( \{ x \in \mathcal{O}_K : \{x, x + p_1(n), \dots, x + p_k(n)\} \subseteq E \} \right)> \delta^{k+1} - \varepsilon \] is syndetic. These two results generalize and extend previous work of Frantzikinakis and Kra on polynomial configurations in $\mathbb{Z}$ and build upon recent work of the authors and Best on linear patterns in general abelian groups. The above combinatorial results follow from multiple recurrence results in ergodic theory, which require a sharpening of existing tools for handling polynomial multiple ergodic averages. A key advancement made in this paper is a new result on the equidistribution of polynomial orbits in nilmanifolds, which can be seen as a far-reaching generalization of Weyl's equidistribution theorem.


title: Dirac equation on a catenoid bridge: a supersymmetric approach

abstract: In this paper, we study the Dirac equation for an electron constrained to move on a catenoid surface. We decoupled the two components of the spinor and obtained two Klein-Gordon-like equations. Analytical solutions were obtained using supersymmetric quantum mechanics for two cases, namely, the constant Fermi velocity and the position-dependent Fermi velocity cases.


title: Bardeen-Cooper-Schrieffer pairing of composite fermions

abstract: Topological pairing of composite fermions has led to remarkable ideas, such as excitations obeying non-Abelian braid statistics and topological quantum computation. We construct a $p$-wave paired Bardeen-Cooper-Schrieffer (BCS) wave function for composite fermions in the torus geometry, which is a convenient geometry for formulating momentum space pairing as well as for revealing the underlying composite-fermion Fermi sea. Following the standard BCS approach, we minimize the Coulomb interaction energy at half filling in the lowest and the second Landau levels, which correspond to filling factors $\nu=1/2$ and $\nu=5/2$ in GaAs quantum wells, by optimizing two variational parameters that are analogous to the gap and the Debye cut-off energy of the BCS theory. Our results show no evidence for pairing at $\nu=1/2$ but a clear evidence for pairing at $\nu=5/2$. To a good approximation, the highest overlap between the exact Coulomb ground state at $\nu=5/2$ and the BCS state is obtained for parameters that minimize the energy of the latter, thereby providing support for the physics of composite-fermion pairing as the mechanism for the $5/2$ fractional quantum Hall effect. We discuss the issue of modular covariance of the composite-fermion BCS wave function, and calculate its Hall viscosity and pair correlation function. By similar methods, we look for but do not find an instability to $s$-wave pairing for a spin-singlet composite-fermion Fermi sea at half-filled lowest Landau level in a system where the Zeeman splitting has been set to zero.


title: Triple top quark production in Standard Model

abstract: The triple top quark production processes have been investigated and calculated in the scope of the Standard Model. The cross sections for different channels are provided for the proton-proton collision energies of $\sqrt{s}=$ 14 and 100 TeV. The importance of the electroweak contribution has been demonstrated. For the main channels the interference between the gluon and the weak bosons mediated contributions is negative and significant, therefore the complete set of diagrams has to be taken into account. Estimated total rates are about 1.9 fb for $\sqrt{s}=14$ TeV and 530 fb for $\sqrt{s}=100$ TeV. A simple estimation of the uncertainty of the calculated cross sections gives about 20\%. The integrated luminosity of 3 ab$^{-1}$ at HL-LHC allows expecting about 5700 events, giving a chance to detect this rare process.


title: Evaluation of Human-AI Teams for Learned and Rule-Based Agents in Hanabi

abstract: Deep reinforcement learning has generated superhuman AI in competitive games such as Go and StarCraft. Can similar learning techniques create a superior AI teammate for human-machine collaborative games? Will humans prefer AI teammates that improve objective team performance or those that improve subjective metrics of trust? In this study, we perform a single-blind evaluation of teams of humans and AI agents in the cooperative card game Hanabi, with both rule-based and learning-based agents. In addition to the game score, used as an objective metric of the human-AI team performance, we also quantify subjective measures of the human's perceived performance, teamwork, interpretability, trust, and overall preference of AI teammate. We find that humans have a clear preference toward a rule-based AI teammate (SmartBot) over a state-of-the-art learning-based AI teammate (Other-Play) across nearly all subjective metrics, and generally view the learning-based agent negatively, despite no statistical difference in the game score. This result has implications for future AI design and reinforcement learning benchmarking, highlighting the need to incorporate subjective metrics of human-AI teaming rather than a singular focus on objective task performance.


title: Velocity Dispersions of Clusters in the Dark Energy Survey Y3 redMaPPer Catalog

abstract: We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $\lambda$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $\lambda < 70$ systems in our sample, but more than half (55\%) of $\lambda < 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature.


title: Northern Hemisphere UHECR data further supports CMBR photons with weak U(1) component

abstract: I further test the theory that the Cosmic Background Radiation (CBR) is made of photons with a weak U(1) component. Tipler (2005) has previously argued that the consistency of the Standard Model (SM) with the Second Law of Thermodynamics requires the early universe CBR to be composed mainly of an SU(2) gauge field. The U(1) field would be suppressed. One of the consequences of this approach is that the Ultra High Energy Cosmic Rays (UHECR) would be able to propagate much further than conventionally accepted, as an SU(2) dominant CBR would be largely unable to couple with right handed fermions. I test if this novel theory solves the problem of UHECR origin by finding suitable candidates up to a redshift z = 0.1 within three degrees of the arrival direction. Utilizing the Fly's Eye Northern Hemisphere UHECR data, I identified candidates with 80\% success for the Northern Sky UHECR (98.7\% if certain celestial objects which are likely to be Active Galactic Nuclei are absolutely identified as such). This is parsimonious with the CBR theory, which has other important implications for the Standard Model, early universe cosmology, and the origin of matter and anti-matter. This builds off of Tipler and Piasecki (2018), where we have reviewed UHECR data from the Southern Hemisphere to get a 86\% success identification rate. We predict that the remaining UHECR not paired with a potential source will have sources identified upon closer telescopic investigations of these regions. Other recent experiments further suggest an SU(2) dominant composition of the CBR photons. The problem of UHECR origin may be solved.


title: Statistical Learning for Fluid Flows: Sparse Fourier divergence-free approximations

abstract: We reconstruct the velocity field of incompressible flows given a finite set of measurements. For the spatial approximation, we introduce the Sparse Fourier divergence-free (SFdf) approximation based on a discrete $L^2$ projection. Within this physics-informed type of statistical learning framework, we adaptively build a sparse set of Fourier basis functions with corresponding coefficients by solving a sequence of minimization problems where the set of basis functions is augmented greedily at each optimization problem. We regularize our minimization problems with the seminorm of the fractional Sobolev space in a Tikhonov fashion. In the Fourier setting, the incompressibility (divergence-free) constraint becomes a finite set of linear algebraic equations. We couple our spatial approximation with the truncated Singular Value Decomposition (SVD) of the flow measurements for temporal compression. Our computational framework thus combines supervised and unsupervised learning techniques. We assess the capabilities of our method in various numerical examples arising in fluid mechanics.


title: Multi-task Learning with Cross Attention for Keyword Spotting

abstract: Keyword spotting (KWS) is an important technique for speech applications, which enables users to activate devices by speaking a keyword phrase. Although a phoneme classifier can be used for KWS, exploiting a large amount of transcribed data for automatic speech recognition (ASR), there is a mismatch between the training criterion (phoneme recognition) and the target task (KWS). Recently, multi-task learning has been applied to KWS to exploit both ASR and KWS training data. In this approach, an output of an acoustic model is split into two branches for the two tasks, one for phoneme transcription trained with the ASR data and one for keyword classification trained with the KWS data. In this paper, we introduce a cross attention decoder in the multi-task learning framework. Unlike the conventional multi-task learning approach with the simple split of the output layer, the cross attention decoder summarizes information from a phonetic encoder by performing cross attention between the encoder outputs and a trainable query sequence to predict a confidence score for the KWS task. Experimental results on KWS tasks show that the proposed approach achieves a 12% relative reduction in the false reject ratios compared to the conventional multi-task learning with split branches and a bi-directional long short-team memory decoder.


title: AGNIFS survey of local AGN: GMOS-IFU data and outflows in 30 sources

abstract: We analyse optical datacubes of the inner kiloparsec of 30 local ($z\le0.02$) active galactic nuclei (AGN) hosts that our research group, AGNIFS, has collected over the past decade via observations with the integral field units of the Gemini Multi-Object Spectrographs. Spatial resolutions range between $50~{\rm pc}$ and $300~{\rm pc}$ and spectral coverage is from $4800~\mathring{A}$ or $5600~\mathring{A}$ to $7000~\mathring{A}$, at velocity resolutions of $\approx 50~{\rm km~s^{-1}}$. We derive maps of the gas excitation and kinematics, determine the AGN ionisation axis -- which has random orientation relative to the galaxy, and the kinematic major axes of the emitting gas. We find that rotation dominates the gas kinematics in most cases, but is disturbed by the presence of inflows and outflows. Outflows have been found in 21 nuclei, usually along the ionisation axis. The gas velocity dispersion is traced by $W_{80}$ (velocity width encompassing 80 per cent of the line flux), adopted as a tracer of outflows. In 7 sources $W_{80}$ is enhanced perpendicularly to the ionisation axis, indicating lateral expansion of the outflow. We have estimated mass-outflow rates $\dot{M}$ and powers $\dot{E}$, finding median values of $\log\,[\dot{M}/({\rm\,M_\odot\,yr^{-1}})]=-2.1_{-1.0}^{+1.6}$ and $\log\,[\dot{E}/({\rm\,erg\,s^{-1}})]=38.5_{-0.9}^{+1.8}$, respectively. Both quantities show a mild correlation with the AGN luminosity ($L_{\rm AGN}$). $\dot{E}$ is of the order of 0.01 $L_{\rm AGN}$ for 4 sources, but much lower for the majority (9) of the sources, with a median value of $\log\,[\dot{E}/L_{\rm AGN}]=-5.34_{-0.9}^{+3.2}$ indicating that typical outflows in the local Universe are unlikely to significantly impact their host galaxy evolution.


title: An undecidable extension of Morley's theorem on the number of countable models

abstract: We show that Morley's theorem on the number of countable models of a countable first-order theory becomes an undecidable statement when extended to second-order logic. More generally, we calculate the number of equivalence classes of $\sigma$-projective equivalence relations in several models of set theory. Our methods include random and Cohen forcing, Woodin cardinals and Inner Model Theory.


title: A Proof of Merca's Conjectures on Sums of Odd Divisor Functions

abstract: In a recent paper, Merca posed three conjectures on congruences for specific convolutions of a sum of odd divisor functions with a generating function for generalized $m$-gonal numbers. Extending Merca's work, we complete the proof of these conjectures.


title: Quasi Differential Quotients

abstract: We explore basic properties and some applications of Quasi Differential Quotients ($QDQ$s) and the related $QDQ$-approximating multi-cones. A $QDQ$, which is a special kind of H.Sussmann's Approximate Generalized Differential Quotient ($AGDQ$), consists in a notion of generalized differentiation for set-valued maps. $QDQ$s have the advantage over $AGDQ$s of allowing a genuine, non-punctured, Open Mapping result, so implying stronger set-separation theorems. They have already proved quite useful in the investigation of some connections occurring between infimum gap phenomena and the normality of minima. Moreover, $QDQ$-approximating multi-cones are fit in optimal control to deduce Maximum Principles that involve (set-valued) Lie brackets of nonsmooth vector fields.


title: On abelian $\ell$-towers of multigraphs III

abstract: Let $\ell$ be a rational prime. Previously, abelian $\ell$-towers of multigraphs were introduced which are analogous to $\Z_{\ell}$-extensions of number fields. It was shown that for towers of bouquets, the growth of the $\ell$-part of the number of spanning trees behaves in a predictable manner (analogous to a well-known theorem of Iwasawa for $\Z_{\ell}$-extensions of number fields). In this paper, we extend this result to abelian $\ell$-towers over an arbitrary connected multigraph (not necessarily simple and not necessarily regular). In order to carry this out, we employ integer-valued polynomials to construct power series with coefficients in $\Z_\ell$ arising from cyclotomic number fields, different than the power series appearing in the prequel. This allows us to study the special value at $u=1$ of the Artin--Ihara $L$-function, when the base multigraph is not necessarily a bouquet.


title: Obtaining Causal Information by Merging Datasets with MAXENT

abstract: The investigation of the question "which treatment has a causal effect on a target variable?" is of particular relevance in a large number of scientific disciplines. This challenging task becomes even more difficult if not all treatment variables were or even cannot be observed jointly with the target variable. Another similarly important and challenging task is to quantify the causal influence of a treatment on a target in the presence of confounders. In this paper, we discuss how causal knowledge can be obtained without having observed all variables jointly, but by merging the statistical information from different datasets. We show how the maximum entropy principle can be used to identify edges among random variables when assuming causal sufficiency and an extended version of faithfulness, and when only subsets of the variables have been observed jointly.


title: The Extinction and Distance of the MBM Molecular Clouds at High Galactic Latitude

abstract: Based on the accurate color excess $E_{\rm G_{BP},G_{RP}}$ of more than 4 million stars and $E_{\rm NUV,G_{BP}}$ of more than 1 million stars from \citet{2021ApJS..254...38S}, the distance and the extinction of the molecular clouds in the MBM catalog at $|b|>20^{\circ}$ are studied in combination with the distance measurement of \emph{Gaia}/EDR3. The distance as well as the color excess is determined for 66 molecular clouds. The color excess ratio $E_{\rm G_{BP},G_{RP}}/E_{\rm NUV,G_{BP}}$ is derived for 39 of them, which is obviously larger and implies more small particles at smaller extinction. In addition, the scale height of the dust disk is found to be about 100 pc and becomes large at the anticenter direction due to the disk flaring.


title: Improving application performance with biased distributions of quantum states

abstract: We consider the properties of a specific distribution of mixed quantum states of arbitrary dimension that can be biased towards a specific mean purity. In particular, we analyze mixtures of Haar-random pure states with Dirichlet-distributed coefficients. We analytically derive the concentration parameters required to match the mean purity of the Bures and Hilbert--Schmidt distributions in any dimension. Numerical simulations suggest that this value recovers the Hilbert--Schmidt distribution exactly, offering an alternative and intuitive physical interpretation for ensembles of Hilbert--Schmidt-distributed random quantum states. We then demonstrate how substituting these Dirichlet-weighted Haar mixtures in place of the Bures and Hilbert--Schmidt distributions results in measurable performance advantages in machine-learning-based quantum state tomography systems and Bayesian quantum state reconstruction. Finally, we experimentally characterize the distribution of quantum states generated by both a cloud-accessed IBM quantum computer and an in-house source of polarization-entangled photons. In each case, our method can more closely match the underlying distribution than either Bures or Hilbert--Schmidt distributed states for various experimental conditions.


title: The principal Erd\H{o}s--Gallai differences of a degree sequence

abstract: The Erd\H{o}s--Gallai criteria for recognizing degree sequences of simple graphs involve a system of inequalities. Given a fixed degree sequence, we consider the list of differences of the two sides of these inequalities. These differences have appeared in varying contexts, including characterizations of the split and threshold graphs, and we survey their uses here. Then, enlarging upon properties of these graph families, we show that both the last term and the maximum term of the principal Erd\H{o}s--Gallai differences of a degree sequence are preserved under graph complementation and are monotonic under the majorization order and Rao's order on degree sequences.


title: TNQMetro: Tensor-network based package for efficient quantum metrology computations

abstract: TNQMetro is a numerical package written in Python for calculations of fundamental quantum bounds on measurement precision. Thanks to the usage of the tensor-network formalism it can beat the curse of dimensionality and provides an efficient framework to calculate bounds for finite size system as well as determine the asymptotic scaling of precision in systems where quantum enhancement amounts to a constant factor improvement over the Standard Quantum Limit. It is written in a user-friendly way so that the basic functions do not require any knowledge of tensor networks.


title: Next generation design and prospects for CANNEX

abstract: The Casimir And Non-Newtonian force EXperiment (CANNEX) implements the unique geometry of macroscopic plane parallel plates that guarantees an optimum sensitivity with respect to interfacial forces and their gradients. Based on experience from the recently completed proof-of-principle phase, we have started a re-design of the setup aiming to reduce systematic effects and maximize the achievable sensitivity. Several propositions have been made to measure Casimir forces in and out of thermal equilibrium, hypothetical axion and axion-like dark matter interactions, and forces originating from chameleon or symmetron dark energy interactions. In the present article, we give details on the design for the next implementation stage of CANNEX and discuss the experimental opportunities as well as limitations expected for this new setup.


title: Ecohydrological land reanalysis

abstract: The accurate estimation of terrestrial water and vegetation is a grand challenge in hydrometeorology. Many previous studies developed land data assimilation systems (LDASs) and provided global-scale land surface datasets by integrating numerical simulation and satellite data. However, vegetation dynamics has not been explicitly solved in these land reanalysis datasets. Here we present the newly developed land reanalysis dataset, ECoHydrological Land reAnalysis (ECHLA). ECHLA is generated by sequentially assimilating C- and X- band microwave brightness temperature satellite observations into a land surface model which can explicitly simulate the dynamic evolution of vegetation biomass. The ECHLA dataset provides semi-global soil moisture from surface to 1.95m depth, Leaf Area Index (LAI), and vegetation water content and is available from 2003 to 2010 and from 2013 to 2019. We assess the performance of ECHLA to estimate soil moisture and vegetation dynamics by comparing the ECHLA dataset with independent satellite and in-situ observation data. We found that our sequential update by data assimilation substantially improves the skill to reproduce the seasonal cycle of vegetation. Data assimilation also contributes to improving the skill to simulate soil moisture mainly in the shallow soil layers (0-0.15m depth). The ECHLA dataset will be publicly available and expected to contribute to understanding terrestrial ecohydrological cycles and water-related natural disasters such as drought.


title: An Energy-Efficient Edge Computing Paradigm for Convolution-based Image Upsampling

abstract: A novel energy-efficient edge computing paradigm is proposed for real-time deep learning-based image upsampling applications. State-of-the-art deep learning solutions for image upsampling are currently trained using either resize or sub-pixel convolution to learn kernels that generate high fidelity images with minimal artifacts. However, performing inference with these learned convolution kernels requires memory-intensive feature map transformations that dominate time and energy costs in real-time applications. To alleviate this pressure on memory bandwidth, we confine the use of resize or sub-pixel convolution to training in the cloud by transforming learned convolution kernels to deconvolution kernels before deploying them for inference as a functionally equivalent deconvolution. These kernel transformations, intended as a one-time cost when shifting from training to inference, enable a systems designer to use each algorithm in their optimal context by preserving the image fidelity learned when training in the cloud while minimizing data transfer penalties during inference at the edge. We also explore existing variants of deconvolution inference algorithms and introduce a novel variant for consideration. We analyze and compare the inference properties of convolution-based upsampling algorithms using a quantitative model of incurred time and energy costs and show that using deconvolution for inference at the edge improves both system latency and energy efficiency when compared to their sub-pixel or resize convolution counterparts.