title: Boson stars and black holes with wavy scalar hair

abstract: In this paper, we follow up on the discovery of a new type of solution in the Einstein-Maxwell system coupled minimally to a self-interacting complex scalar field. For sufficiently large gravitational coupling and sufficiently small electromagnetic coupling we demonstrate that boson stars as well as black holes can carry scalar hair that shows a distinct new feature~: a number of spatial oscillations in the scalar field away from the core or horizon, respectively. These spatial oscillations appear also in the curvature invariants and hence should be a detectable feature of the space-time. As a first hint that this is true, we show that the effective potential for null geodesics in this space-time possesses a local minimum indicating that in the spatial region where oscillations occur a new stable photon sphere should be possible. We also study the interior of the black holes with scalar hair and show that the curvature singularity appears at a finite value of the radius and that black holes with wavy scalar hair have this singularity very close to the center.


title: Analysis of a diffuse interface method for the Stokes-Darcy coupled problem

abstract: We consider the interaction between a free flowing fluid and a porous medium flow, where the free flowing fluid is described using the time dependent Stokes equations, and the porous medium flow is described using Darcy's law in the primal formulation. To solve this problem numerically, we use a diffuse interface approach, where the weak form of the coupled problem is written on an extended domain which contains both Stokes and Darcy regions. This is achieved using a phase-field function which equals one in the Stokes region and zero in the Darcy region, and smoothly transitions between these two values on a diffuse region of width $\mathcal{O}(\epsilon)$ around the Stokes-Darcy interface. We prove convergence of the diffuse interface formulation to the standard, sharp interface formulation, and derive rates of convergence. This is performed by deriving a priori error estimates for discretizations of the diffuse interface method, and by analyzing the modeling error of the diffuse interface approach at the continuous level. The convergence rates are also shown computationally in a numerical example.


title: Correlative electron and light spectroscopy of strongly-coupled mid-infrared plasmon and phonon polaritons

abstract: The ability to control and modify infrared excitations in condensed matter is of both fundamental and application interests. Here we explore a system supporting low-energy excitations, in particular, mid-infrared localized plasmon modes and phonon polaritons that are tuned to be strongly coupled. We study the coupled modes by using far-field infrared spectroscopy, state-of-the-art monochromated electron energy-loss spectroscopy, numerical simulations and analytical modeling. We demonstrate that the electron probe facilitates a precise characterization of polaritons constituting the coupled system, and enables an active control over the coupling and the resulting sample response both in frequency and space. Our work establishes a rigorous description of the spectral features observed in light- and localized electron-based spectroscopies, which can be extended to analogous optical systems with applications in heat management, and electromagnetic field concentration or nanofocusing.


title: Dense Out-of-Distribution Detection by Robust Learning on Synthetic Negative Data

abstract: Standard machine learning is unable to accommodate inputs which do not belong to the training distribution. The resulting models often give rise to confident incorrect predictions which may lead to devastating consequences. This problem is especially demanding in the context of dense prediction since input images may be only partially anomalous. Previous work has addressed dense out-of-distribution detection by discriminative training with respect to off-the-shelf negative datasets. However, real negative data are unlikely to cover all modes of the entire visual world. To this end, we extend this approach by generating synthetic negative patches along the border of the inlier manifold. We leverage a jointly trained normalizing flow due to coverage-oriented learning objective and the capability to generate samples at different resolutions. We detect anomalies according to a principled information-theoretic criterion which can be consistently applied through training and inference. The resulting models set the new state of the art on benchmarks for out-of-distribution detection in road-driving scenes and remote sensing imagery, in spite of minimal computational overhead.


title: Software Engineering Education Knowledge Versus Industrial Needs

abstract: Contribution: Determine and analyze the gap between software practitioners' education outlined in the 2014IEEE/ACM Software Engineering Education Knowledge (SEEK) and industrial needs pointed by Wikipedia articles referenced in Stack Overflow (SO) posts. Background: Previous work has uncovered deficiencies in the coverage of computer fundamentals, people skills, software processes, and human-computer interaction, suggesting rebalancing. Research Questions: 1) To what extent are developers' needs, in terms of Wikipedia articles referenced in SO posts, covered by the SEEK knowledge units? 2) How does the popularity of Wikipedia articles relate to their SEEK coverage? 3) What areas of computing knowledge can be better covered by the SEEK knowledge units? 4) Why are Wikipedia articles covered by the SEEK knowledge units cited on SO? Methodology: Wikipedia articles were systematically collected from SO posts. The most cited were manually mapped to the SEEK knowledge units, assessed according to their degree of coverage. Articles insufficiently covered by the SEEK were classified by hand using the 2012 ACM Computing Classification System. A sample of posts referencing sufficiently covered articles was manually analyzed. A survey was conducted on software practitioners to validate the study findings. Findings: SEEK appears to cover sufficiently computer science fundamentals, software design and mathematical concepts, but less so areas like the World Wide Web, software engineering components, and computer graphics. Developers seek advice, best practices and explanations about software topics, and code review assistance. Future SEEK models and the computing education could dive deeper in information systems, design, testing, security, and soft skills.


title: Practical limits on Nanosatellite Telescope Pointing: The Impact of Disturbances and Photon Noise

abstract: Accurate and stable spacecraft pointing is a requirement of many astronomical observations. Pointing particularly challenges nanosatellites because of an unfavorable surface area to mass ratio and proportionally large volume required for even the smallest attitude control systems. This work explores the limitations on astrophysical attitude knowledge and control in a regime unrestricted by actuator precision or actuator-induced disturbances such as jitter. The external disturbances on an archetypal 6U CubeSat are modeled and the limiting sensing knowledge is calculated from the available stellar flux and grasp of a telescope within the available volume. These inputs are integrated using a model-predictive control scheme. For a simple test case at 1 Hz, with an 85 mm telescope and a single 11th magnitude star, the achievable body pointing is predicted to be 0.39 arcseconds. For a more general limit, integrating available star light, the achievable attitude sensing is approximately 1 milliarcsecond, which leads to a predicted body pointing accuracy of 20 milliarcseconds after application of the control model. These results show significant room for attitude sensing and control systems to improve before astrophysical and environmental limits are reached.


title: Towards Hardware Support for FPGA Resource Elasticity

abstract: FPGAs are increasingly being deployed in the cloud to accelerate diverse applications. They are to be shared among multiple tenants to improve the total cost of ownership. Partial reconfiguration technology enables multi-tenancy on FPGA by partitioning it into regions, each hosting a specific application's accelerator. However, the region's size can not be changed once they are defined, resulting in the underutilization of FPGA resources. This paper argues to divide the acceleration requirements of an application into multiple small computation modules. The devised FPGA shell can reconfigure the available PR regions with those modules and enable them to communicate with each other over Crossbar interconnect with the Wishbone bus interface. For each PR region being reconfigured, it updates the register file with the valid destination addresses and the bandwidth allocation of the interconnect. Any invalid communication request originating from the Wishbone master interface is masked in the corresponding master port of the crossbar. The allocated bandwidth for the PR region is ensured by the weighted round-robin arbiter in the slave port of the crossbar. Finally, the envisioned resource manager can increase or decrease the number of PR regions allocated to an application based on its acceleration requirements and PR regions' availability.


title: P-REx II. Off-line Performance on VLTI/GRAVITY

abstract: For sensitive optical interferometry, it is crucial to control the evolution of the optical path difference (OPD) of the wavefront between the individual telescopes of the array. The OPD between a pair of telescopes is induced by differential optical properties such as atmospheric refraction, telescope alignment, etc. This has classically been measured using a fringe tracker that provides corrections to a piston actuator to account for this difference. An auxiliary method, known as the Piston Reconstruction Experiment (P-REx) has been developed to measure the OPD, or differential 'piston' of the wavefront, induced by the atmosphere at each telescope. Previously, this method was outlined and results obtained from LBT adaptive optics (AO) data for a single telescope aperture were presented. P-REx has now been applied off-line to previously acquired VLT's GRAVITY CIAO wavefront sensing data to estimate the atmospheric OPD for the six VLTI baselines. Comparisons with the OPD obtained from the VLTI GRAVITY fringe tracker were made. The results indicate that the telescope and instrumental noise of the combined VLTI and GRAVITY systems dominate over the atmospheric turbulence contributions. However, good agreement between simulated and on-sky P-REx data indicates that if the telescope and instrumental noise were reduced to atmospheric piston noise levels, P-REx has the potential to reduce the OPD root mean square of piston turbulence by up to a factor of 10 for frequencies down to 1 Hz. In such conditions, P-REx will assist in pushing the sensitivity limits of optical fringe tracking with long baseline interferometers.


title: Performing Bayesian analyses with AZURE2 using BRICK: an application to the ${}^7$Be system

abstract: Phenomenological $R$-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical challenge has always been the estimation of the uncertainty on both the cross sections at the energies of interest and the fit parameters, which can take the form of standard level parameters. Frequentist ($\chi^2$-based) estimation has been the norm. In this work, a Markov Chain Monte Carlo sampler, \texttt{emcee}, has been implemented for the $R$-matrix code \texttt{AZURE2}, creating the Bayesian $R$-matrix Inference Code Kit (\texttt{BRICK}). Bayesian uncertainty estimation has then been carried out for a simultaneous $R$-matrix fit of the $^3$He$(\alpha,\gamma)^7$Be and $^3$He$(\alpha,\alpha)^3$He reactions in order to gain further insight into the fitting of capture and scattering data. Both data sets constrain the values of the bound state $\alpha$-particle asymptotic normalization coefficients in $^7$Be. The analysis highlights the need for low-energy scattering data with well-documented uncertainty information and shows how misleading results can be obtained in its absence.


title: Faster Deep Ensemble Averaging for Quantification of DNA Damage from Comet Assay Images With Uncertainty Estimates

abstract: Several neurodegenerative diseases involve the accumulation of cellular DNA damage. Comet assays are a popular way of estimating the extent of DNA damage. Current literature on the use of deep learning to quantify DNA damage presents an empirical approach to hyper-parameter optimization and does not include uncertainty estimates. Deep ensemble averaging is a standard approach to estimating uncertainty but it requires several iterations of network training, which makes it time-consuming. Here we present an approach to quantify the extent of DNA damage that combines deep learning with a rigorous and comprehensive method to optimize the hyper-parameters with the help of statistical tests. We also use an architecture that allows for a faster computation of deep ensemble averaging and performs statistical tests applicable to networks using transfer learning. We applied our approach to a comet assay dataset with more than 1300 images and achieved an $R^2$ of 0.84, where the output included the confidence interval for each prediction. The proposed architecture is an improvement over the current approaches since it speeds up the uncertainty estimation by 30X while being statistically more rigorous.


title: The behavior of essential dimension under specialization

abstract: Let $A$ be a discrete valuation ring with generic point $\eta$ and closed point $s$. We show that in a family of torsors over $\operatorname{Spec}(A)$, the essential dimension of the torsor above $s$ is less than or equal to the essential dimension of the torsor above $\eta$. We give two applications of this result, one in mixed characteristic, the other in equal characteristic.


title: ABC of the Future

abstract: Approximate Bayesian computation (ABC) has advanced in two decades from a seminal idea to a practically applicable inference tool for simulator-based statistical models, which are becoming increasingly popular in many research domains. The computational feasibility of ABC for practical applications has been recently boosted by adopting techniques from machine learning to build surrogate models for the approximate likelihood or posterior and by the introduction of a general-purpose software platform with several advanced features, including automated parallelization. Here we demonstrate the strengths of the advances in ABC by going beyond the typical benchmark examples and considering real applications in astronomy, infectious disease epidemiology, personalised cancer therapy and financial prediction. We anticipate that the emerging success of ABC in producing actual added value and quantitative insights in the real world will continue to inspire a plethora of further applications across different fields of science, social science and technology.


title: Recurrent Neural Networks (RNNs) with dimensionality reduction and break down in computational mechanics; application to multi-scale localization step

abstract: Artificial Neural Networks (NNWs) are appealing functions to substitute high dimensional and non-linear history-dependent problems in computational mechanics since they offer the possibility to drastically reduce the computational time. This feature has recently been exploited in the context of multi-scale simulations, in which the NNWs serve as surrogate model of micro-scale finite element resolutions. Nevertheless, in the literature, mainly the macro-stress-macro-strain response of the meso-scale boundary value problem was considered and the micro-structure information could not be recovered in a so-called localization step. In this work, we develop Recurrent Neural Networks (RNNs) as surrogates of the RVE response while being able to recover the evolution of the local micro-structure state variables for complex loading scenarios. The main difficulty is the high dimensionality of the RNNs output which consists in the internal state variable distribution in the micro-structure. We thus propose and compare several surrogate models based on a dimensionality reduction: i) direct RNN modeling with implicit NNW dimensionality reduction, ii) RNN with PCA dimensionality reduction, and iii) RNN with PCA dimensionality reduction and dimensionality break down, i.e. the use of several RNNs instead of a single one. Besides, we optimize the sequential training strategy of the latter surrogate for GPU usage in order to speed up the process. Finally, through RNN modeling of the principal components coefficients, the connection between the physical state variables and the hidden variables of the RNN is revealed, and exploited in order to select the hyper-parameters of the RNN-based surrogate models in their design stage.


title: Impact of class imbalance on chest x-ray classifiers: towards better evaluation practices for discrimination and calibration performance

abstract: This work aims to analyze standard evaluation practices adopted by the research community when assessing chest x-ray classifiers, particularly focusing on the impact of class imbalance in such appraisals. Our analysis considers a comprehensive definition of model performance, covering not only discriminative performance but also model calibration, a topic of research that has received increasing attention during the last years within the machine learning community. Firstly, we conducted a literature study to analyze common scientific practices and confirmed that: (1) even when dealing with highly imbalanced datasets, the community tends to use metrics that are dominated by the majority class; and (2) it is still uncommon to include calibration studies for chest x-ray classifiers, albeit its importance in the context of healthcare. Secondly, we perform a systematic experiment on two major chest x-ray datasets to explore the behavior of several performance metrics under different class ratios and show that widely adopted metrics can conceal the performance in the minority class. Finally, we recommend the inclusion of complementary metrics to better reflect the system's performance in such scenarios. Our study indicates that current evaluation practices adopted by the research community for chest x-ray computer-aided diagnosis systems may not reflect their performance in real clinical scenarios, and suggest alternatives to improve this situation.


title: Emerging Devices and Packaging Strategies for Electronic-Photonic AI Accelerators

abstract: The field of mimicking the structure of the brain on a chip is experiencing interest driven by the demand for machine intelligent applications. However, the power consumption and available performance of machine-learning (ML) accelerating hardware still leave much desire for improvement. In this letter, we share viewpoints, challenges, and prospects of electronic-photonic neural network (NN) accelerators. Combining electronics with photonics offers synergistic co-design strategies for high-performance AI Application-specific integrated circuits (ASICs) and systems. Taking advantages of photonic signal processing capabilities and combining them with electronic logic control and data storage is an emerging prospect. However, the optical component library leaves much to be desired and is challenged by the enormous size of photonic devices. Within this context, we will review the emerging electro-optic materials, functional devices, and systems packaging strategies that, when realized, provide significant performance gains and fuel the ongoing AI revolution, leading to a stand-alone photonics-inside AI ASIC 'black-box' for streamlined plug-and-play board integration in future AI processors.


title: Automatic Meta-Path Discovery for Effective Graph-Based Recommendation

abstract: Heterogeneous Information Networks (HINs) are labeled graphs that depict relationships among different types of entities (e.g., users, movies and directors). For HINs, meta-path-based recommenders (MPRs) utilize meta-paths (i.e., abstract paths consisting of node and link types) to predict user preference, and have attracted a lot of attention due to their explainability and performance. We observe that the performance of MPRs is highly sensitive to the meta-paths they use, but existing works manually select the meta-paths from many possible ones. Thus, to discover effective meta-paths automatically, we propose the Reinforcement learning-based Meta-path Selection (RMS) framework. Specifically, we define a vector encoding for meta-paths and design a policy network to extend meta-paths. The policy network is trained based on the results of downstream recommendation tasks and an early stopping approximation strategy is proposed to speed up training. RMS is a general model, and it can work with all existing MPRs. We also propose a new MPR called RMS-HRec, which uses an attention mechanism to aggregate information from the meta-paths. We conduct extensive experiments on real datasets. Compared with the manually selected meta-paths, the meta-paths identified by RMS consistently improve recommendation quality. Moreover, RMS-HRec outperforms state-of-the-art recommender systems by an average of 7% in hit ratio. The codes and datasets are available on https://github.com/Stevenn9981/RMS-HRec.


title: Relic galaxy analogues in Illustris TNG50: the formation pathways of surviving red nuggets in a cosmological simulation

abstract: Relic galaxies are massive compact quiescent galaxies that formed at high-redshift and remained almost unchanged since then. In this work, we search for analogues to relic galaxies in the TNG50 cosmological simulations to understand relic formation and test the ability of TNG50 to reproduce such rare objects. Using stellar mass, age, radius, quiescence and stellar assembly criteria, we find 5 subhalos in TNG50 that could be potential relic analogues. We compare their properties with other constraints imposed by a sample of 13 observed relic galaxies. We find one analogue in TNG50 that simultaneously satisfies most of the available observational constraints, such as metallicity and morphology. It also shows similarities to the confirmed relic NGC 1277, regarding environment and dark matter fraction. By taking into account a degree of relicness, a second relic analogue may also be considered. However, the central parts of the brightness and density profiles of the analogues are less steep than that of real relic galaxies, possibly due to limited numerical resolution. We identify two formation pathways of relic analogues in TNG50 depending on their environment: they either have their remaining gas stripped during the infall into a cluster at $z \lesssim 1.2$ or consume it before $z > 1.5$. They are then deprived of significant star formation, leaving their stellar populations almost unaltered during the last 9 Gyr. We also find that the analogue progenitors at $z \sim 4$ inhabit more massive halos than progenitors of quiescent galaxies with similar stellar mass at $z \sim 0$.


title: The Extended Local Supercluster

abstract: It has long been established but seldom noticed that we are in a region at least 170 Mpc across in which different types of galaxies show different degrees of alignment with the plane of the de Vaucouleurs Local Supercluster. While clusters of galaxies and radio galaxies at redshifts less than 0.02 are concentrated at low supergalactic latitudes, the most luminous galaxies in the infrared, LIRGs, show little correlation with this plane. The most luminous early-type galaxies are concentrated at low supergalactic latitudes, but similarly luminous spirals are not noticeably so. The cross-correlations of the positions of what might be termed field galaxies with positions of clusters and LIRGs offer a measure of the situation. The mean density at distance 0.5 Mpc from a LIRG is comparable to the mean density at that distance from a cluster of galaxies, but the mean density 5 Mpc from a LIRG is well below the mean density at that distance from a cluster and not much greater than the cosmic mean density. Discussion of issues arising is brief.


title: Neutrino Tagging, a new tool for accelerator based neutrino experiments

abstract: This article describes a new experimental method for accelerator based neutrino experiments called neutrino tagging. The method consists in exploiting the neutrino production mechanism, the $\pi^{\pm}\to\mu^{\pm}\nu_\mu$ decay, to kinematically reconstruct the neutrino properties from the decay incoming and outgoing charged particles. The reconstruction of these particles relies on the recent progress and on-going developments in silicon particle detector technology. A detailed description of the method and achievable key performances is presented, together with its potential benefits for short and long baseline experiments. Then, a novel configuration for long baseline experiments is discussed in which a tagged beam would be employed together with mega-ton scale natural deep water Cherenkov detectors. The coarseness of this type of detectors is overcome by the precision of the tagging and, conversely, the rate limitation imposed by the tagging is outweighed by the virtually unlimited size of the detector. These mutual benefits result in an affordable design for next generations of long based line experiments. The physics potential of such experiments is quantified using the Protvino to KM3NeT/ORCA setup as a case study for which an unprecedented sensitivity to the leptonic CP violation could be achieved.


title: A first-order condition for the independence on $p$ of weak gradients

abstract: It is well known that on arbitrary metric measure spaces, the notion of minimal $p$-weak upper gradient may depend on $p$. In this paper we investigate how a first-order condition of the metric-measure structure, that we call Bounded Interpolation Property, guarantees that in fact such dependence is not present. We also show that the Bounded Interpolation Property is stable for pointed measure Gromov Hausdorff convergence and holds on a large class of spaces satisfying curvature dimension conditions.


title: Skyrmion lattice in centrosymmetric magnets with local Dzyaloshinsky-Moriya interaction

abstract: It is common for the local inversion symmetry to break in crystals, even though the whole crystal has global inversion symmetry. This local inversion symmetry breaking allows for a local Dzyaloshinsky-Moriya interaction (DMI) in magnetic crystals. Here we show that the local DMI can stabilize a skyrmion as a metastable excitation or as a skyrmion crystal in equilibrium. We consider the crystal structure with layered structure as an example, where local inversion is violated in each layer but a global inversion center exists in the middle of the two layers. These skyrmions come in pairs that are related by inversion symmetry. The two skyrmions with opposite helicity in a pair form a bound state. We study the properties of a skyrmion pair in the ferromagnetic background and determine the equilibrium phase diagram, where a robust lattice of skyrmion pairs is stabilized. Our results point to a new direction to search for the skyrmion lattice in centrosymmetric magnets.


title: On The Limiting Distribution of Free Path Lengths for Flat Surfaces with Circular Obstacles

abstract: In this note, we prove the existence of a limiting distribution of the free path lengths on flat surfaces with circular obstacles as the radius of the obstacles goes to zero. Moreover, we relate this distribution to the distribution of the heights of zippered rectangle decompositions of flat surfaces.


title: Asymptotics of quantum invariants of surface diffeomorphisms I: conjecture and algebraic computations

abstract: The Kashaev-Murakami-Murakami Volume Conjecture connects the hyperbolic volume of a knot complement to the asymptotics of certain evaluations of the colored Jones polynomials of the knot. We introduce a closely related conjecture for diffeomorphisms of surfaces, backed up by numerical evidence. The conjecture involves isomorphisms between certain representations of the Kauffman bracket skein algebra of the surface, and the bulk of the article is devoted to the development of explicit methods to compute these isomorphisms. These combinatorial and algebraic techniques are exploited in two subsequent articles, which prove the conjecture for a large family of diffeomorphisms of the one-puncture torus and are much more analytic.


title: Systolic blood pressure estimation using ECG and PPG in patients undergoing surgery

abstract: Background and Objectives: In a significant portion of surgeries, blood pressure (BP) is often measured non-invasively in an intermittent manner. This practice has a risk of missing clinically relevant BP changes between two adjacent intermittent BP measurements. This study proposes a method to non-invasively estimate systolic blood pressure (SBP) with high accuracy in patients undergoing surgery. Methods: Continuous arterial BP, electrocardiography (ECG), and photoplethysmography (PPG) signals were acquired from 29 patients undergoing surgery. After extracting 9 features from the PPG and ECG signals, we dynamically selected features upon each intermittent measurement (every 10 min) of SBP based on feature robustness and the principle of correlation-based feature selection. Finally, multiple linear regression models were built to combine these features to estimate SBP every 30 s. Results: Compared to the reference SBP, the proposed method achieved a mean of difference at 0.08 mmHg, a standard deviation of difference at 7.97 mmHg, and a correlation coefficient at 0.89 (p < 0.001). Conclusions: This study demonstrates the feasibility of non-invasively estimating SBP every 30 s with high accuracy during surgery by using ECG, PPG, and intermittent SBP measurements every 10 min, which meets the standard of the Association for the Advancement of Medical Instrumentation. The proposed method has the potential to enhance BP monitoring in the operating room, improving patient outcomes and experiences.


title: A fast computational model for the electrophysiology of the whole human heart

abstract: In this study we present a novel computational model for unprecedented simulations of the whole cardiac electrophysiology. According to the heterogeneous electrophysiologic properties of the heart, the whole cardiac geometry is decomposed into a set of coupled conductive media having different topology and electrical conductivities: (i) a network of slender bundles comprising a fast conduction atrial network, the AV-node and the ventricular bundles; (ii) the Purkinje network; and (iii) the atrial and ventricular myocardium. The propagation of the action potential in these conductive media is governed by the bidomain/monodomain equations, which are discretized in space using an in-house finite volume method and coupled to three different cellular models, the Courtemanche model [1] for the atrial myocytes, the Stewart model [2] for the Purkinje Network and the ten Tusscher-Panfilov model [3] for the ventricular myocytes. The developed numerical model correctly reproduces the cardiac electrophysiology of the whole human heart in healthy and pathologic conditions and it can be tailored to study and optimize resynchronization therapies or invasive surgical procedures. Importantly, the whole solver is GPU-accelerated using CUDA Fortran providing an unprecedented speedup, thus opening the way for systematic parametric studies and uncertainty quantification analyses.


title: SoK: Privacy-preserving Deep Learning with Homomorphic Encryption

abstract: Outsourced computation for neural networks allows users access to state of the art models without needing to invest in specialized hardware and know-how. The problem is that the users lose control over potentially privacy sensitive data. With homomorphic encryption (HE) computation can be performed on encrypted data without revealing its content. In this systematization of knowledge, we take an in-depth look at approaches that combine neural networks with HE for privacy preservation. We categorize the changes to neural network models and architectures to make them computable over HE and how these changes impact performance. We find numerous challenges to HE based privacy-preserving deep learning such as computational overhead, usability, and limitations posed by the encryption schemes.


title: LFT Representation of a Class of Nonlinear Systems: A Data-Driven Approach

abstract: This paper focuses on developing a method to obtain an uncertain linear fractional transformation (LFT) system that adequately captures the dynamics of a nonlinear time-invariant system over some desired envelope. First, the nonlinear system is approximated as a polynomial nonlinear state-space (PNLSS) system, and a linear parameter-varying (LPV) representation of the PNLSS model is obtained. To reduce the potentially large number of scheduling parameters in the resulting LPV system, an approach based on the cascade feedforward neural network (CFNN) is proposed. We account for the approximation and reduction errors through the addition of a norm-bounded, causal, dynamic uncertainty in the LFT system. Falsification is used in a novel way to perform guided simulations for deriving a norm bound on the dynamic uncertainty and generating data for training the CFNN. Then, robustness analysis using integral quadratic constraint theory is carried out to choose an LFT representation that leads to a computationally tractable analysis problem and useful analysis results. Finally, the proposed approach is used to obtain an LFT representation for the nonlinear equations of motion of a fixed-wing unmanned aircraft system.


title: Variational quantum circuits to prepare low energy symmetry states

abstract: We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a Symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(~{\alpha}) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ hamiltonian using rotation and reflection symmetry and H2 hamiltonian within S_z = 0 subspace using S^2 symmetry. We have found the variationally trained unitary surprisingly giving very good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.


title: Chance and the Continuum Hypothesis

abstract: This paper presents and defends an argument that the continuum hypothesis is false, based on considerations about objective chance and an old theorem due to Banach and Kuratowski. More specifically, I argue that the probabilistic inductive methods standardly used in science presuppose that every proposition about the outcome of a chancy process has a certain chance between 0 and 1. I also argue in favour of the standard view that chances are countably additive. Since it is possible to randomly pick out a point on a continuum, for instance using a roulette wheel or by flipping a countable infinity of fair coins, it follows, given the axioms of ZFC, that there are many different cardinalities between countable infinity and the cardinality of the continuum.


title: Cardinality of Real Numbers and Set Axiomatization

abstract: The cardinality of real numbers depends on the chosen axiomatic system and specific models of set theory. Under the ZF framework, real numbers are nondenumerable, but their exact cardinality remains debated due to possible extensions of the ZF axioms. We explored this issue by analyzing sets within the ZFC model and Wang's $\Sigma$ model. Set $S'$, within the ZFC framework, includes all pairs of algebraic numbers with a transcendental number inserted between each pair by a function $f$, which inserts an arbitrarily chosen transcendental number into every interval bounded by algebraic endpoints. This approach enables the representation of any number using these intervals. By employing the AC, one can identify intervals from any Cauchy sequence within $S'$. As each interval contains at least one number due to the function $f$, the intersections of these intervals contain numbers defined by these sequences, making a complete and countable set $S'$. The second set, created by inserting specifically defined transcendental numbers between algebraic numbers, starts with the sequence of algebraic numbers as the 0-th layer in Wang's $\Sigma$ model. A diagonal procedure applied to that layer generates a set of transcendental numbers which, combined with the 0-th layer, forms the $\Sigma_1$ layer. Continuing this recursive process results in a set equivalent to $S'$. This suggests that all constructive real numbers are countable, challenging the traditional diagonal argument. The AC is crucial for selecting transcendental numbers and identifying intervals and sequences, making the set $S'$. This might resonate with constructivists who oppose the AC. However, further analysis revealed conflicts between proofs of non-denumerability and the Nested Interval Property, indicating a need to reevaluate the countability of real numbers.


abstract: We prove versions of Ekeland, Takahashi and Caristi principles in preordered quasi-metric spaces, the equivalence between these principles, as well as their equivalence to some completeness results for the underlying quasi-metric space. These extend the results proved in S.~Cobza\c{s}, Topology Appl. \textbf{265} (2019), 106831, 22, for quasi-metric spaces. The key tools are Picard sequences for some special set-valued mappings on a preordered quasi-metric space $X$, defined in terms of the preorder and of a function $\varphi$ on $X$. Key words: preordered quasi-metric space; completeness in quasi-metric spaces; variational principles; Ekeland variational principle; Takahashi minimization principle; fixed point; Caristi fixed point theorem.


title: Quantum geodesics on $\lambda$-Minkowski spacetime

abstract: We apply a recent formalism of quantum geodesics to the well-known bicrossproduct model $\lambda$-Minkowski quantum spacetime $[x^i,t]=\imath\lambda_p x^i$ with its flat quantum metric as a model of quantum gravity effects, with $\lambda_p$ the Planck scale. As examples, quantum geodesic flow of a plane wave gets an order $\lambda_p$ frequency dependent correction to the classical geodesic velocity. A quantum geodesic flow with classical velocity $v$ of a Gaussian with width $\sqrt{2\beta}$ initially centred at the origin changes its shape but its centre of mass moves with ${<x>\over<t>}=v(1+{\lambda_p^2\over 2\beta}+O(\lambda^3_p))$, an order $\lambda_p^2$ correction. This implies, at least within perturbation theory, that a `point particle' cannot be modelled as an infinitely sharp Gaussian due to quantum gravity corrections. For contrast, we also look at quantum geodesics on the noncommutative torus with a 2D curved weak quantum Levi-Civita connection.


title: On soliton solutions and soliton interactions of Kulish-Sklyanin and Hirota-Ohta systems

abstract: In this paper we consider a simplest two-dimensional reduction of the remarkable three-dimensional Hirota-Ohta system. The Lax pair of the Hirota-Ohta system was extended to a Lax triad by adding extra third linear equation, whose compatibility conditions with the Lax pair of the Hirota-Ohta imply another remarkable systems: the Kulish-Sklyanin system (KSS) together with its first higher commuting flow, which we can call as vector complex MKdV. This means that any common particular solution of these both two-dimensional integrable systems yields a corresponding particular solution of the three-dimensional Hirota-Ohta system. Using the dressing Zakharov-Shabat method we derive the $N$-soliton solutions of these systems and analyze their interactions, i.e. derive explicitly the shifts of the relative center-of-mass coordinates and the phases as functions of the discrete eigenvalues of the Lax operator. Next we relate to these NLEE system of Hirota--Ohta type and obtain its $N$-soliton solutions


title: Two-center basis generator method calculations for Li$^{3+}$, C$^{3+}$ and O$^{3+}$ ion impact on ground state hydrogen

abstract: The two-center basis generator method is used to obtain cross sections for excitation, capture, and ionization in Li$^{3+}$, C$^{3+}$, and O$^{3+}$ collisions with ground-state hydrogen at projectile energies from 1 to 100 keV/u. The interaction of the C$^{3+}$ and O$^{3+}$ projectiles with the active electron is represented by a model potential. Comparisons of cross sections with previously reported data show overall good agreement while discrepancies in capture for C$^{3+}$ collisions at low energies are noted. The present results show that excitation and ionization are similar across the three collision systems, which indicates that these cross sections are mostly dependent on the net charge of the projectile only. The situation is different for the capture channel.


title: Rank growth of elliptic curves over $N$-th root extensions

abstract: Fix an elliptic curve $E$ over a number field $F$ and an integer $n$ which is a power of $3$. We study the growth of the Mordell--Weil rank of $E$ after base change to the fields $K_d = F(\sqrt[2n]{d})$. If $E$ admits a $3$-isogeny, then we show that the average ``new rank'' of $E$ over $K_d$, appropriately defined, is bounded as the height of $d$ goes to infinity. When $n = 3$, we moreover show that for many elliptic curves $E/\mathbb{Q}$, there are no new points on $E$ over $\mathbb{Q}(\sqrt[6]d)$, for a positive proportion of integers $d$. This is a horizontal analogue of a well-known result of Cornut and Vatsal. As a corollary, we show that Hilbert's tenth problem has a negative solution over a positive proportion of pure sextic fields $\mathbb{Q}(\sqrt[6]{d})$. The proofs combine our recent work on ranks of abelian varieties in cyclotomic twist families with a technique we call the ``correlation trick'', which applies in a more general context where one is trying to show simultaneous vanishing of multiple Selmer groups. We also apply this technique to families of twists of Prym surfaces, which leads to bounds on the number of rational points in sextic twist families of bielliptic genus 3 curves.


title: Absorption by stringy black holes

abstract: We investigate the absorption of planar massless scalar waves by a charged rotating stringy black hole, namely a Kerr--Sen black hole. We compute numerically the absorption cross section and compare our results with those of the Kerr-Newman black hole, a classical general relativity solution. In order to better compare both charged black holes, we define the ratio of the black hole charge to the extreme charge as $Q$. We conclude that Kerr--Sen and Kerr-Newman black holes have a similar absorption cross section, with the difference increasing for higher values of $Q$.


title: Quantum nonequilibrium dynamics from Knizhnik-Zamolodchikov equations

abstract: We consider a set of non-stationary quantum models. We show that their dynamics can be studied using links to Knizhnik-Zamolodchikov (KZ) equations for correlation functions in conformal field theories. We specifically consider the boundary Wess-Zumino-Novikov-Witten model, where equations for correlators of primary fields are defined by an extension of KZ equations and explore the links to dynamical systems. As an example of the workability of the proposed method, we provide an exact solution to a dynamical system that is a specific multi-level generalization of the two-level Landau-Zenner system known in the literature as the Demkov-Osherov model. The method can be used to study the nonequilibrium dynamics in various multi-level systems from the solution of the corresponding KZ equations.

categories: quant-ph cond-mat.quant-gas hep-th math-ph math.MP


title: HSPACE: Synthetic Parametric Humans Animated in Complex Environments

abstract: Advances in the state of the art for 3d human sensing are currently limited by the lack of visual datasets with 3d ground truth, including multiple people, in motion, operating in real-world environments, with complex illumination or occlusion, and potentially observed by a moving camera. Sophisticated scene understanding would require estimating human pose and shape as well as gestures, towards representations that ultimately combine useful metric and behavioral signals with free-viewpoint photo-realistic visualisation capabilities. To sustain progress, we build a large-scale photo-realistic dataset, Human-SPACE (HSPACE), of animated humans placed in complex synthetic indoor and outdoor environments. We combine a hundred diverse individuals of varying ages, gender, proportions, and ethnicity, with hundreds of motions and scenes, as well as parametric variations in body shape (for a total of 1,600 different humans), in order to generate an initial dataset of over 1 million frames. Human animations are obtained by fitting an expressive human body model, GHUM, to single scans of people, followed by novel re-targeting and positioning procedures that support the realistic animation of dressed humans, statistical variation of body proportions, and jointly consistent scene placement of multiple moving people. Assets are generated automatically, at scale, and are compatible with existing real time rendering and game engines. The dataset with evaluation server will be made available for research. Our large-scale analysis of the impact of synthetic data, in connection with real data and weak supervision, underlines the considerable potential for continuing quality improvements and limiting the sim-to-real gap, in this practical setting, in connection with increased model capacity.


title: Discovery of periodic variables in the hot subdwarf catalogue through a systematic search in Zwicky Transient Facility data

abstract: We conduct a systematic search for periodic variables in the hot subdwarf catalogue using data from the Zwicky Transient Facility. We present the classification of 67 HW Vir binaries, 496 reflection effect, pulsation or rotation sinusoids, 11 eclipsing signals, and 4 ellipsoidally modulated binaries. Of these, 486 are new discoveries that have not been previously published including a new mass-transferring hot subdwarf binary candidate. These sources were determined by applying the Lomb-Scargle and Box Least Squares periodograms along with manual inspection. We calculated variability statistics on all periodic sources, and compared our results to traditional methods of determining astrophysical variability. We find that $\approx60$\% percent of variable targets, mostly sinusoidal variability, would have been missed using a traditional $varindex$ cut. Most HW Virs, eclipsing systems and all ellipsoidal variables were recovered with a $varindex\,>0.02$. We also find a significant reddening effect, with some variable hot subdwarfs meshing with the main sequence stripe in the Hertzprung-Russell Diagram. Examining the positions of the variable stars in Galactic coordinates, we discover a higher proportion of variable stars within $|b|<25^\circ$ of the Galactic Plane, suggesting that the Galactic Plane may be fertile grounds for future discoveries if photometric surveys can effectively process the clustered field.


title: A Lightweight Approach to Computing Message Races with an Application to Causal-Consistent Reversible Debugging

abstract: This paper presents a lightweight formalism (a trace) to model message-passing concurrent executions where some common common problems can be identified, like lost or delayed messages, some forms of deadlock, etc. In particular, we consider (potential) message races that can be useful to analyze alternative executions. We consider a particular application for our developments in the context of a causal-consistent reversible debugging framework for Erlang programs


title: Measuring Attribution in Natural Language Generation Models

abstract: With recent improvements in natural language generation (NLG) models for various applications, it has become imperative to have the means to identify and evaluate whether NLG output is only sharing verifiable information about the external world. In this work, we present a new evaluation framework entitled Attributable to Identified Sources (AIS) for assessing the output of natural language generation models, when such output pertains to the external world. We first define AIS and introduce a two-stage annotation pipeline for allowing annotators to appropriately evaluate model output according to AIS guidelines. We empirically validate this approach on generation datasets spanning three tasks (two conversational QA datasets, a summarization dataset, and a table-to-text dataset) via human evaluation studies that suggest that AIS could serve as a common framework for measuring whether model-generated statements are supported by underlying sources. We release guidelines for the human evaluation studies.


title: Fundamental Analytic Solutions for the Kulish-Sklyanin Model with Constant Boundary Conditions

abstract: In the present paper we analyze the construction of fundamental analytic solutions (FAS) for the generalized Kulish-Sklyanin models (KSM) for vanishing (VBC) and constant boundary conditions (CBC). Using FAS one can reduce the direct and inverse scattering problems for the Lax operator to a Riemann-Hilbert problem (RHP). For VBC we find two FAS $\chi^+(x,t,\lambda) $ and $\chi^-(x,t,\lambda) $ analytic in the upper/lower $\mathbb{C}_\pm$ complex $\lambda$-plane. The RHP consists in: given the sewing function $G(x,t,\lambda)$ to constructing both $\chi^\pm(x,t,\lambda) $ in their regions of analyticity. For CBC the problem becomes more complicated, because now the RHP must be formulated on a Riemannian surface of genus 1.


title: Sparsified Secure Aggregation for Privacy-Preserving Federated Learning

abstract: Secure aggregation is a popular protocol in privacy-preserving federated learning, which allows model aggregation without revealing the individual models in the clear. On the other hand, conventional secure aggregation protocols incur a significant communication overhead, which can become a major bottleneck in real-world bandwidth-limited applications. Towards addressing this challenge, in this work we propose a lightweight gradient sparsification framework for secure aggregation, in which the server learns the aggregate of the sparsified local model updates from a large number of users, but without learning the individual parameters. Our theoretical analysis demonstrates that the proposed framework can significantly reduce the communication overhead of secure aggregation while ensuring comparable computational complexity. We further identify a trade-off between privacy and communication efficiency due to sparsification. Our experiments demonstrate that our framework reduces the communication overhead by up to 7.8x, while also speeding up the wall clock training time by 1.13x, when compared to conventional secure aggregation benchmarks.


title: Network analysis of nanoscale energy conversion processes

abstract: Energy conversion in nanosized devices is studied in the framework of state-space models. We use a network representation of the underlying master equation to describe the dynamics by a graph. Particular segments of this network represent input and output processes that provide a way to introduce a coupling to several heat reservoirs and particle reservoirs. In addition, the network representation scheme allows one to decompose the stationary dynamics as cycles. The cycle analysis is a convenient tool for analyse models of machine operations, which are characterized by different nanoscale energy conversion processes. By introducing the cycle affinity, we are able to calculate the zero-current limit. The zero-current limit can be mapped to the zero-affinity limit in a network representation scheme. For example, for systems with competing external driving forces the open-circuit voltage can be determined by setting the cycle affinity zero. This framework is used to derive open-circuit voltage with respect to microscopic material energetics and different coupling to particle and temperature reservoirs.


title: ATP Synthase: A Moonlighting Enzyme with Unprecedented Functions

abstract: ATP synthase's intrinsic molecular electrostatic potential (MESP) adds constructively to, and hence reinforces, the chemiosmotic voltage. This ATP synthase voltage represents a new free energy term that appears to have been overlooked. This term is at least roughly equal in order of magnitude and opposite in sign to the energy needed to be dissipated as a Maxwell's demon (Landauer principle).


title: A novel Bayesian approach to the computation of the configurational density of states

abstract: In this work we develop and implement a novel Bayesian method for computing the DOS of a system. This method is based on the use of a test function with adjustable parameters and we use Bayes theorem to find the best parameters given a certain number of measurements done on the system. This measurements can be done in any ensemble defined by a distribution function. We found that the algorithm can find the DOS in a reasonable amount of time, and that if the test function is suitable enough, the DOS found by the algorithm is very close to the true DOS.


title: Learning to Walk with Dual Agents for Knowledge Graph Reasoning

abstract: Graph walking based on reinforcement learning (RL) has shown great success in navigating an agent to automatically complete various reasoning tasks over an incomplete knowledge graph (KG) by exploring multi-hop relational paths. However, existing multi-hop reasoning approaches only work well on short reasoning paths and tend to miss the target entity with the increasing path length. This is undesirable for many reason-ing tasks in real-world scenarios, where short paths connecting the source and target entities are not available in incomplete KGs, and thus the reasoning performances drop drastically unless the agent is able to seek out more clues from longer paths. To address the above challenge, in this paper, we propose a dual-agent reinforcement learning framework, which trains two agents (GIANT and DWARF) to walk over a KG jointly and search for the answer collaboratively. Our approach tackles the reasoning challenge in long paths by assigning one of the agents (GIANT) searching on cluster-level paths quickly and providing stage-wise hints for another agent (DWARF). Finally, experimental results on several KG reasoning benchmarks show that our approach can search answers more accurately and efficiently, and outperforms existing RL-based methods for long path queries by a large margin.


title: An orbital-based representation for accurate Quantum Machine Learning

abstract: We introduce an electronic structure based representation for quantum machine learning (QML) of electronic properties throughout chemical compound space. The representation is constructed using computationally inexpensive ab initio calculations and explicitly accounts for changes in the electronic structure. We demonstrate the accuracy and flexibility of resulting QML models when applied to property labels such as total potential energy, HOMO and LUMO energies, ionization potential, and electron affinity, using as data sets for training and testing entries from the QM7b, QM7b-T, QM9, and LIBE libraries. For the latter, we also demonstrate the ability of this approach to account for molecular species of different charge and spin multiplicity, resulting in QML models that infer total potential energies based on geometry, charge, and spin as input.


title: Dynamics of Li deposition on epitaxial graphene/Ru(0001) islands

abstract: Li metal has been deposited on the surface of a Ru(0001) single crystal containing patches of monolayer-thick epitaxial graphene islands. The use of low-energy electron microscopy and diffraction allowed us to {\em in situ} monitor the process by measuring the local work function as well as to study the system in real and reciprocal space, comparing the changes taking place on the graphene with those on the bare Ru(0001) surface. It is found that Li deposition decreases the work function of the graphene islands but to a much smaller degree than of the Ru(0001) surface, as corresponds to its intercalation below the graphene overlayer. Finally, the diffusion process of Li out of the graphene islands has been monitored by photoelectron microscopy using a visible-light laser.


title: Evolution of Transthyretin

abstract: Evolution of the amino acid sequences of transthyretin (TTR) can provide additional information about its dynamics that both complements and extends the already extensive static structural data. Protein dynamics is largely driven by interactions between the protein itself and the thin water film that covers it. Here those interactions are connected to profiles of water waves connecting domain pivots (hydrophobic extrema), as well as a broad hydrophilic central hinge region that has evolved to provide human TTR with greater flexibility and stability. This central region has a density of single mutations related to amyloid polyneuropathy that is three times higher than other regions.


title: Vegetation Impact on Atmospheric Moisture Transport under Increasing Land-Ocean Temperature Contrasts

abstract: Destabilization of the water cycle threatens human lives and livelihoods. Meanwhile our understanding of whether and how changes in vegetation cover could trigger abrupt transitions in moisture regimes remains incomplete. This challenge calls for better evidence as well as for the theoretical concepts to describe it. Here we briefly summarise the theoretical questions surrounding the role of vegetation cover in the dynamics of a moist atmosphere. We discuss the previously unrecognized sensitivity of local wind power to condensation rate as revealed by our analysis of the continuity equation for a gas mixture. Using the framework of condensation-induced atmospheric dynamics, we then show that with the temperature contrast between land and ocean increasing up to a critical threshold, ocean-to-land moisture transport reaches a tipping point where it can stop or even reverse. Land-ocean temperature contrasts are affected by both global and regional processes, in particular, by the surface fluxes of sensible and latent heat that are strongly influenced by vegetation. Our results clarify how a disturbance of natural vegetation cover, e.g., by deforestation, can disrupt large-scale atmospheric circulation and moisture transport. In view of the increasing pressure on natural ecosystems, successful strategies of mitigating climate change require taking into account the impact of vegetation on moist atmospheric dynamics. Our analysis provides a theoretical framework to assess this impact. The available data for Eurasia indicate that the observed climatological land-ocean temperature contrasts are close to the threshold. This can explain the increasing fluctuations in the continental water cycle including droughts and floods and signifies a yet greater potential importance for large-scale forest conservation.


title: Classifying invariants for $E_1$: A tail of a generic real

abstract: Let $E$ be an analytic equivalence relation on a Polish space. We introduce a framework for studying the possible "reasonable" complete classifications and the complexity of possible classifying invariants for $E$, such that: (1) the standard results and intuitions regarding classifications by countable structures are preserved in this framework; (2) this framework respects Borel reducibility; (3) this framework allows for a precise study of the possible invariants of certain equivalence relations which are not classifiable by countable structures, such as $E_1$. In this framework we show that $E_1$ can be classified, with classifying invariants which are $\kappa$-sequences of $E_0$-classes where $\kappa=\mathfrak{b}$, and it cannot be classified in such a manner if $\kappa<\mathbf{add}(\mathcal{B})$. These results depend on analyzing the following sub-model of a Cohen real extension, introduced by Kanovei-Sabok-Zapletal (2013) and Larson-Zapletal (2020). Let $\left<c_n:\,n<\omega\right>$ be a generic sequence of Cohen reals, and define the tail intersection model $$M=\bigcap_{n<\omega}V[\left<c_m:\,m\geq n\right>].$$ An analysis of reals in $M$ will provide lower bounds for the possible invariants for $E_1$. We also extend the characterization of turbulence from Larson-Zapletal (2020) in terms of intersection models.


title: Exploration of Spinodal Decomposition in Multi-Principal Element Alloys (MPEAs) using CALPHAD Modeling

abstract: Researchers attributed the orderly arranged nanoscale phases observed in many multi-principal element alloys (MPEAs) to spinodal/spinodal-mediated phase transformation pathways. However, spinodal decomposition is not well understood in multicomponent systems. Although the theoretical background is available, CALPHAD databases were not used to explore the miscibility gap in MPEAs. In this work, we develop a CALPHAD framework utilizing the Hessian of free energy to study the stability of solid solutions in MPEAs. In particular, we utilize the geometry of higher dimensional Gibbs simplex in conjunction with the Hessian to calculate concentration modulations in early stages of spinodal decomposition. We apply this framework to a diverse set of multi-phase MPEAs that have been studied in the literature including TiZrNbTa (BCC/BCC), Fe15Co15Ni20Mn20Cu30 (FCC/FCC), Al0.5NbTa0.8Ti1.5V0.2Zr (BCC/B2), AlCo0.4Cr0.6FeNi (BCC/B2) and Al0.5Cr0.9FeNi2.5V0.2 (FCC/L12). We show that the MPEA systems are unstable only to certain concentration modulations which could be further explored to design microstructurally engineered alloys.


title: Backhaul-aware Drone Base Station Placement and Resource Management for FSO-based Drone-assisted Mobile Networks

abstract: In drone-assisted mobile networks, Drone-mounted Base Stations (DBSs) are responsively and flexibly deployed over any Places of Interest (PoI), such as sporadic hotspots and disaster-struck areas, where the existing mobile network infrastructure is unable to provide wireless coverage. In this paper, a DBS is an aerial base station to relay traffic between a nearby Macro Base Station (MBS) and the users. In addition, Free Space Optics (FSO) is applied as the backhauling solution to significantly increase the capacity of the backhaul link between an MBS and a DBS. Most of the existing DBS placement solutions assume the FSO-based backhaul link provides sufficient link capacity, which may not be true, especially when a DBS is placed far away from an MBS (e.g., > 10 km in disaster-struck areas) or in a bad weather condition. In this paper, we formulate a problem to jointly optimize bandwidth allocation and DBS placement by considering the FSO-based backhaul link capacity constraint. A Backhaul awaRe bandwidth allOcAtion and DBS placement (BROAD) algorithm is designed to efficiently solve the problem, and the performance of the algorithm is demonstrated via extensive simulations.


title: Progressively Enlargement of Filtrations and Control Problems for Step Processes

abstract: In the present paper we address stochastic optimal control problems for a step process $(X,\mathbb{F})$ under a progressive enlargement of the filtration. The global information is obtained adding to the reference filtration $\mathbb{F}$ the point process $H=1_{[\tau,+\infty)}$. Here $\tau$ is a random time that can be regarded as the occurrence time of an external shock event. We study two classes of control problems, over $[0,T]$ and over the random horizon $[0,T \wedge \tau]$. We solve these control problems following a dynamical approach based on a class of BSDEs driven by the jump measure $\mu^ Z$ of the semimartingale $Z=(X,H)$, which is a step process with respect to the enlarged filtration $\mathbb G$. The BSDEs that we consider can be solved in $\mathbb{G}$ thanks to a martingale representation theorem which we also establish here. To solve the BSDEs and the control problems we need to ensure that $Z$ is quasi-left continuous in the enlarged filtration $\mathbb{G}$. Therefore, in addition to the $\mathbb{F}$-quasi left continuity of $X$, we assume some further conditions on $\tau$: the {\it avoidance} of $\mathbb{F}$-stopping times and the {\it immersion} property, or alternatively {\it Jacod's absolutely continuity} hypothesis.


title: Monotonicity of Steklov eigenvalues on graphs and applications

abstract: In this paper, we obtain monotonicity of Steklov eigenvalues on graphs which as a special case on trees extends the results of He-Hua [Calc. Var. Partial Differential Equations 61 (2022), no. 3, Paper No. 101, arXiv: 2103.07696] to higher Steklov eigenvalues and gives affirmative answers to two problems proposed in He-Hua [arXiv: 2103.07696]. As applications of the monotonicity of Steklov eigenvalues, we obtain some estimates for Steklov eigenvalues on trees generalizing the isodiametric estimate for the first positive Steklov eigenvalues on trees in He-Hua [arXiv:2011.11014].


title: Rediscovering Affordance: A Reinforcement Learning Perspective

abstract: Affordance refers to the perception of possible actions allowed by an object. Despite its relevance to human-computer interaction, no existing theory explains the mechanisms that underpin affordance-formation; that is, how affordances are discovered and adapted via interaction. We propose an integrative theory of affordance-formation based on the theory of reinforcement learning in cognitive sciences. The key assumption is that users learn to associate promising motor actions to percepts via experience when reinforcement signals (success/failure) are present. They also learn to categorize actions (e.g., "rotating" a dial), giving them the ability to name and reason about affordance. Upon encountering novel widgets, their ability to generalize these actions determines their ability to perceive affordances. We implement this theory in a virtual robot model, which demonstrates human-like adaptation of affordance in interactive widgets tasks. While its predictions align with trends in human data, humans are able to adapt affordances faster, suggesting the existence of additional mechanisms.


title: A formal approach to good practices in Pseudo-Labeling for Unsupervised Domain Adaptive Re-Identification

abstract: The use of pseudo-labels prevails in order to tackle Unsupervised Domain Adaptive (UDA) Re-Identification (re-ID) with the best performance. Indeed, this family of approaches has given rise to several UDA re-ID specific frameworks, which are effective. In these works, research directions to improve Pseudo-Labeling UDA re-ID performance are varied and mostly based on intuition and experiments: refining pseudo-labels, reducing the impact of errors in pseudo-labels... It can be hard to deduce from them general good practices, which can be implemented in any Pseudo-Labeling method, to consistently improve its performance. To address this key question, a new theoretical view on Pseudo-Labeling UDA re-ID is proposed. The contributions are threefold: (i) A novel theoretical framework for Pseudo-Labeling UDA re-ID, formalized through a new general learning upper-bound on the UDA re-ID performance. (ii) General good practices for Pseudo-Labeling, directly deduced from the interpretation of the proposed theoretical framework, in order to improve the target re-ID performance. (iii) Extensive experiments on challenging person and vehicle cross-dataset re-ID tasks, showing consistent performance improvements for various state-of-the-art methods and various proposed implementations of good practices.


title: A flexible and low-cost open-source IPMC mezzanine for ATCA boards based on OpenIPMC

abstract: This work presents the development of an Intelligent Platform Management Controller mezzanine in a Mini DIMM form factor for use in electronic boards compliant to the PICMG Advanced Telecommunication Computing Architecture (ATCA) standard. The module is based on an STMicroelectronics STM32H745 microcontroller running the OpenIPMC open-source software. The mezzanine has been successfully tested on a variety of ATCA boards being proposed for the upgrade of the experiments at the HL-LHC, with its design and firmware being distributed under open-source hardware license.


title: Spaces of extremal magnitude

abstract: Magnitude is a numerical invariant of compact metric spaces. Its theory is most mature for spaces satisfying the classical condition of being of negative type, and the magnitude of such a space lies in the interval $[1, \infty]$. Until now, no example with magnitude $\infty$ was known. We construct some, thus answering a question open since 2010. We also give a sufficient condition for the magnitude of a space to converge to 1 as it is scaled down to a point, unifying and generalizing previously known conditions.


title: Comparison of the material quality of AlxIn1-xN (x ~ 0-0.50) films deposited on Si (100) and (111) by reactive RF sputtering

abstract: Here, we compare the material quality of AlxIn1-xN layers deposited on Si with different crystallographic orientations, (100) and (111), via radio-frequency (RF) sputtering. To modulate their Al content, the Al RF power was varied from 0 to 225 W, whereas the In RF power and deposition temperature were fixed at 30 W and 300oC, respectively. X-ray diffraction measurements reveal a c-axis-oriented wurtzite structure with no phase separation regardless of the Al content (x = 0-0.50), which increases with the Al power supply. The surface morphology of the AlxIn1-xN layers improves with increasing Al content and it is similar for samples grown on both Si substrates (the root-mean-square roughness decreases from 12 nm to 2.5 nm). Furthermore, from TEM images we notice a similar grain-like columnar morphology and defect density on samples deposited on both Si substrates under the same conditions. Simultaneously grown AlxIn1-xN-on-sapphire samples point to a residual n-type carrier concentration in the 1020-1021 cm-3 range. The optical band gap energy of these layers evolves from 1.75 eV to 2.56 eV with increasing Al content, consistent with the blue shift of their low-temperature photoluminescence. In general, the material quality of the AlxIn1-xN films on Si is similar for both crystallographic orientations. Nonetheless, samples deposited on sapphire show an improved structural and morphological characteristic likely due to the lower difference in lattice constants between the nitride and the sapphire substrate.


title: Degenerations, fibrations and higher rank Landau-Ginzburg models

abstract: We study semi-stable degenerations of quasi-Fano varieties to unions of two pieces. We conjecture that the higher rank Landau-Ginzburg models mirror to these two pieces can be glued together to lower rank Landau-Ginzburg models which are mirror to the original quasi-Fano varieties. We prove this conjecture by relating their Euler characteristics, generalized functional invariants as well as periods. We also use it to conjecture a relation between the degenerations to the normal cones and the fibrewise compactifications of higher rank Landau-Ginzburg models. Furthermore, we use it to iterate the Doran-Harder-Thompson conjecture and obtain higher codimension Calabi-Yau fibrations.


title: Computing Viscous Flow Along a 3D Open Tube Using the Immerse Interface Method

abstract: In a companion study \cite{patterson2020computing2D}, we present a numerical method for simulating 2D viscous flow through an open compliant closed channel, drive by pressure gradient. We consider the highly viscous regime, where fluid dynamics is described by the Stokes equations, and the less viscous regime described by the Navier-Stokes equations. In this study, we extend the method to 3D tubular flow. The problem is formulated in axisymmetric cylindrical coordinates, an approach that is natural for tubular flow simulations and that substantially reduces computational cost. When the elastic tubular walls are stretched or compressed, they exert forces on the fluid. These singular forces introduce unsmoothness into the fluid solution. As in the companion 2D study \cite{patterson2020computing2D}, we extend the immersed interface method to an open tube, and we compute solution to the model equations using the resulting method. Numerical results indicate that this new method preserves sharp jumps in the solution and its derivatives, and converges with second-order accuracy in both space and time.


title: Revealing sign-reversal $s^{+-}$-wave pairing by quasiparticle interference in the heavy-fermion superconductor CeCu$_2$Si$_2$

abstract: Recent observations of two nodeless gaps in superconducting CeCu$_2$Si$_2$ have raised intensive debates as to its exact gap structure of either sign-reversal ($s^{+-}$) or sign-preserving ($s^{++}$) pairing. Here we investigate the quasiparticle interference (QPI) using realistic Fermi surface topology for both weak and strong interband impurity scatterings. Our calculations of the QPI and integrated antisymmetrized local density of states reveal qualitative distinctions between $s^{+-}$ and $s^{++}$ pairing states, which include the intragap impurity resonance and a significant energy-dependence difference between two gap energies. Our predictions provide a guide for phase-sensitive QPI measurements to uncover decisively the true pairing symmetry in the heavy-fermion superconductor CeCu$_2$Si$_2$.


title: The Role of the Gradient Term of the Bochner-Kodaira Formula in Coherent Sheaf Extension

abstract: In applying the Bochner-Kodaira formula with boundary term to solve the $\bar\partial$ equation with $L^2$ estimates, the gradient term is usually not used. Two potentially important applications of the use of the gradient term are the strong rigidity for holomorphic vector bundles and the very ampleness part of the Fujita conjecture. In this note we use the gradient term to construct holomorphic sections to prove the Thullen-type extension across codimension $1$ for holomorphic vector bundles with Hermitian metric whose curvature is $L^p$ for some $p>1$. This construction of sections points out a typical way of how the gradient term can be used.


title: Wavelet-based estimation of power densities of size-biased data

abstract: We propose a new wavelet-based method for density estimation when the data are size-biased. More specifically, we consider a power of the density of interest, where this power exceeds 1/2. Warped wavelet bases are employed, where warping is attained by some continuous cumulative distribution function. A special case is the conventional orthonormal wavelet estimation, where the warping distribution is the standard continuous uniform. We show that both linear and nonlinear wavelet estimators are consistent, with optimal and/or near-optimal rates. Monte Carlo simulations are performed to compare four special settings which are easy to interpret in practice. An application with a real dataset on fatal traffic accidents involving alcohol illustrates the method. We observe that warped bases provide more flexible and superior estimates for both simulated and real data. Moreover, we find that estimating the power of a density (for instance, its square root) further improves the results.


title: Dewetting Characteristics of Contact Lenses Coated with Wetting Agents

abstract: Hypothesis: Although wetting agents have been developed to limit tear film dewetting over contact lenses, systematic analyses correlating wetting agents properties to mechanisms of the tear film destabilization are not readily available. Clarifying destabilization characteristics across key physio-chemical variables will provide a rational basis for identifying optimal wetting agents. Experiments: We employ an in-house, in vitro platform to comprehensively evaluate drainage and dewetting dynamics of five wetting agents across seventeen different formulations and two model tear film solutions: phosphate-buffered saline (PBS) and artificial tear solution (ATS). We consider the film thickness evolution, film thickness at breakup, dewetted front propagation, and develop correlations to contact angle to compare the samples. Findings: Zwitterionic wetting agents effectively stabilize the tear film by reducing the film thickness at the onset of dewetting, and delaying dewetted region propagation across the lens. Furthermore, tuning wetting agent surface concentrations in binary mixtures can enhance wetting characteristics. Finally, despite disparities in wetting agent molecular properties, the time to dewet $50\%$ of the lens scales linearly with the product of the receding contact angle and contact angle hysteresis. Hence, we fundamentally establish the importance of minimizing both the absolute contact angle values and contact angle hysteresis for effective wetting performance.


title: Concave-Convex PDMP-based sampling

abstract: Recently non-reversible samplers based on simulating piecewise deterministic Markov processes (PDMPs) have shown potential for efficient sampling in Bayesian inference problems. However, there remains a lack of guidance on how to best implement these algorithms. If implemented poorly, the computational costs of simulating event times can out-weigh the statistical efficiency of the non-reversible dynamics. Drawing on the adaptive rejection literature, we propose the concave-convex adaptive thinning approach for simulating a piecewise deterministic Markov process (CC-PDMP). This approach provides a general guide for constructing bounds that may be used to facilitate PDMP-based sampling. A key advantage of this method is its additive structure - adding concave-convex decompositions yields a concave-convex decomposition. This facilitates swapping priors, simple implementation and computationally efficient thinning. In particular, our approach is well suited to local PDMP simulation where known conditional independence of the target can be exploited for potentially huge computational gains. We provide an R package for implementing the CC-PDMP approach and illustrate how our method outperforms existing approaches to simulating events in the PDMP literature.


title: Dynamics of compression-driven gas-liquid displacement in a capillary tube

abstract: We study two-phase displacement via the steady compression of an air reservoir connected to an oil-filled capillary tube. Our experiments and modelling reveal complex displacement dynamics depending on compression rate and reservoir volume that, for large reservoirs, depend on a single dimensionless compressibility number. We identify two distinct displacement regimes, separated by a critical value of the compressibility number. While the subcritical regime exhibits quasi-steady displacement after an initial transient, the supercritical regime exhibits burst-like expulsion.


title: Monitoring Deforestation Using Multivariate Bayesian Online Changepoint Detection with Outliers

abstract: Near real time change detection is important for a variety of Earth monitoring applications and remains a high priority for remote sensing science. Data sparsity, subtle changes, seasonal trends, and the presence of outliers make detecting actual landscape changes challenging. Adams and MacKay (2007) introduced Bayesian Online Changepoint Detection (BOCPD), a computationally efficient, exact Bayesian method for change detection. Incorporation of prior information allows for relaxed dependence on dense data and an extensive stable period, making this method applicable to relatively short time series and multiple changepoint detection. In this paper we conduct BOCPD with a multivariate linear regression framework that supports seasonal trends. We introduce a mechanism to make BOCPD robust against occasional outliers without compromising the computational efficiency of an exact posterior change distribution nor the detection latency. We show via simulations that the method effectively detects change in the presence of outliers. The method is then applied to monitor deforestation in Myanmar where we show superior performance compared to current online changepoint detection methods.


title: Global Existence and Large Time Behaviors of Strong Solutions to the Kinetic Cucker--Smale Model Coupled with the Three Dimensional Incompressible Navier--Stokes Equations

abstract: In this paper, we study global existence and large time behaviors of strong solutions to the kinetic Cucker--Smale model coupled with the three dimensional incompressible Navier--Stokes equations in the whole space. Using the maximal regularity estimate on the Stokes equations, global-in-time strong solutions to the Cauchy problem of the coupled system are obtained under small initial data regime. The optimal decay rate of the system, in the sense that the energy of the system decays at a rate coinciding with the one corresponding to the underling incompressible Navier--Stokes equations without the coupling term, is also analyzed. Compared with previous results set in spatial-periodic or bounded domains, we circumvent the difficulty caused by unboundedness of the domain using the Fourier splitting method.


title: A machine learning analysis of the relationship between some underlying medical conditions and COVID-19 susceptibility

abstract: For the past couple years, the Coronavirus, commonly known as COVID-19, has significantly affected the daily lives of all citizens residing in the United States by imposing several, fatal health risks that cannot go unnoticed. In response to the growing fear and danger COVID-19 inflicts upon societies in the USA, several vaccines and boosters have been created as a permanent remedy for individuals to take advantage of. In this paper, we investigate the relationship between the COVID-19 vaccines and boosters and the total case count for the Coronavirus across multiple states in the USA. Additionally, this paper discusses the relationship between several, selected underlying health conditions with COVID-19. To discuss these relationships effectively, this paper will utilize statistical tests and machine learning methods for analysis and discussion purposes. Furthermore, this paper reflects upon conclusions made about the relationship between educational attainment, race, and COVID-19 and the possible connections that can be established with underlying health conditions, vaccination rates, and COVID-19 total case and death counts.


title: Protonation-induced discrete superconducting phases in bulk FeSe single crystals

abstract: The superconducting transition temperature, $T_{\rm{c}}$, of FeSe can be significantly enhanced several-fold by applying pressure, electron doping, intercalating spacing layer, and reducing dimensionality. Various ordered electronic phases, such as nematicity and spin density waves, have also been observed accompanying high-$T_{\rm{c}}$ superconductivity. Investigation on the evolution of the electronic structure with $T_{\rm{c}}$ is essential to understanding electronic behavior and high-$T_{\rm{c}}$ superconductivity in FeSe and its derived superconductors. In this report, we have found a series of discrete superconducting phases, with a maximum $T_{\rm{c}}$ up to 44 K, in H$^+$-intercalated FeSe single crystals using an ionic liquid gating method. Accompanied with the increase of $T_{\rm{c}}$, suppression of the nematic phase and evolution from non-Fermi-liquid to Fermi-liquid behavior was observed. An abrupt change in the Fermi surface topology was proposed to explain the discrete superconducting phases. A band structure that favors the high-$T_{\rm{c}}$ superconducting phase was also revealed.


title: Visualizing Atomically-Layered Magnetism in CrSBr

abstract: Two-dimensional (2D) materials can host stable, long-range magnetic phases in the presence of underlying magnetic anisotropy. The ability to realize the full potential of 2D magnets necessitates systematic investigation of the role of individual atomic layers and nanoscale inhomogeneity ($\textit{i.e.}$, strain) on the emergence and stability of both intra- and interlayer magnetic phases. Here, we report multifaceted spatial-dependent magnetism in few-layer CrSBr using magnetic force microscopy (MFM) and Monte Carlo-based magnetic simulations. We perform nanoscale visualization of the magnetic sheet susceptibility from raw MFM data and force-distance curves, revealing a characteristic onset of both intra- and interlayer magnetic correlations as a function of temperature and layer-thickness. We demonstrate that the presence of a single uncompensated layer in odd-layer terraces significantly reduces the stability of the low-temperature antiferromagnetic (AFM) phase and gives rise to multiple coexisting magnetic ground states at temperatures close to the bulk N\'eel temperature ($\textit{T}$$_N$). Furthermore, the AFM phase can be reliably suppressed using modest fields (~300 Oe) from the MFM probe, behaving as a nanoscale magnetic switch. Our prototypical study of few-layer CrSBr demonstrates the critical role of layer parity on field-tunable 2D magnetism and provides vital design criteria for future nanoscale magnetic devices. Moreover, we provide a roadmap for using MFM for nano-magnetometry of 2D materials, despite the ubiquitous absence of bulk zero-field magnetism in magnetized sheets.


title: Quantile Regression with Multiple Proxy Variables

abstract: Data integration has become increasingly popular owing to the availability of multiple data sources. This study considered quantile regression estimation when a key covariate had multiple proxies across several datasets. In a unified estimation procedure, the proposed method incorporates multiple proxies that have various relationships with the unobserved covariates. The proposed approach allows the inference of both the quantile function and unobserved covariates. Moreover, it does not require the quantile function's linearity and, simultaneously, accommodates both the linear and nonlinear proxies. Simulation studies have demonstrated that this methodology successfully integrates multiple proxies and revealed quantile relationships for a wide range of nonlinear data. The proposed method is applied to administrative data obtained from the Survey of Household Finances and Living Conditions provided by Statistics Korea, to specify the relationship between assets and salary income in the presence of multiple income records.


title: Statistics and Topology of Fluctuating Ribbons

abstract: Ribbons are a class of slender structures whose length, width, and thickness are widely separated from each other. This scale separation gives a ribbon unusual mechanical properties in athermal macroscopic settings, e.g. it can bend without twisting, but cannot twist without bending. Given the ubiquity of ribbon-like biopolymers in biology and chemistry, here we study the statistical mechanics of microscopic inextensible, fluctuating ribbons loaded by forces and torques. We show that these ribbons exhibit a range of topologically and geometrically complex morphologies exemplified by three phases - a twist-dominated helical phase (HT), a writhe-dominated helical phase (HW), and an entangled phase - that arise as the applied torque and force is varied. Furthermore, the transition from HW to HT phases is characterized by the spontaneous breaking of parity symmetry and the disappearance of perversions that characterize chirality reversals. This leads to a universal response curve of a topological quantity, the link, as a function of the applied torque that is similar to magnetization curves in second-order phase transitions.


title: Efficient decision tree training with new data structure for secure multi-party computation

abstract: We propose a secure multi-party computation (MPC) protocol that constructs a secret-shared decision tree for a given secret-shared dataset. The previous MPC-based decision tree training protocol (Abspoel et al. 2021) requires $O(2^hmn\log n)$ comparisons, being exponential in the tree height $h$ and with $n$ and $m$ being the number of rows and that of attributes in the dataset, respectively. The cause of the exponential number of comparisons in $h$ is that the decision tree training algorithm is based on the divide-and-conquer paradigm, where dummy rows are added after each split in order to hide the number of rows in the dataset. We resolve this issue via secure data structure that enables us to compute an aggregate value for every group while hiding the grouping information. By using this data structure, we can train a decision tree without adding dummy rows while hiding the size of the intermediate data. We specifically describes a decision tree training protocol that requires only $O(hmn\log n)$ comparisons when the input attributes are continuous and the output attribute is binary. Note that the order is now \emph{linear} in the tree height $h$. To demonstrate the practicality of our protocol, we implement it in an MPC framework based on a three-party secret sharing scheme. Our implementation results show that our protocol trains a decision tree with a height of 5 in 33 seconds for a dataset of 100,000 rows and 10 attributes.


title: 3D Point Cloud Reconstruction and SLAM as an Input

abstract: To handle the different types of surface reconstruction tasks, we have replicated as well as modified a few of reconstruction methods and have made comparisons between the traditional method and data-driven method for reconstruction the surface of an object with dense point cloud as input. On top of that, we proposed a system using tightly-coupled SLAM as an input to generate deskewed point cloud and odometry and a Truncated Signed Distance Function based Surface Reconstruction Library. To get higher accuracy, IMU(Inertial Measurement Unit) pre-integration and pose graph optimization are conduct in the SLAM part. With the help of the Robot Operating System, we could build a system containing those two parts, which can conduct a real-time outdoor surface reconstruction.


title: Annealed Leap-Point Sampler for Multimodal Target Distributions

abstract: In Bayesian statistics, exploring multimodal posterior distribution poses major challenges for existing techniques such as Markov Chain Monte Carlo (MCMC). These problems are exacerbated in high-dimensional settings where MCMC methods typically rely upon localised proposal mechanisms. This paper introduces the Annealed Leap-Point Sampler (ALPS), which augments the target distribution state space with modified annealed (cooled) target distributions, in contrast to traditional approaches which have employed tempering. The temperature of the coldest state is chosen such that its corresponding annealed target density can be sufficiently well-approximated by a Laplace approximation. As a result, a Gaussian mixture independence Metropolis-Hastings sampler can perform mode-jumping proposals even in high-dimensional problems. The ability of this procedure to "mode hop" at this super-cold state is then filtered through to the target state using a sequence of tempered targets in a similar way to that in parallel tempering methods. ALPS also incorporates the best aspects of current gold-standard approaches to multimodal sampling in high-dimensional contexts. A theoretical analysis of the ALPS approach in high dimensions is given, providing practitioners with a gauge on the optimal setup as well as the scalability of the algorithm. For a $d$-dimensional problem the it is shown that the coldest inverse temperature level required for the ALPS only needs to be linear in the dimension, $\mathcal{O}(d)$, and this means that for a collection of multimodal problems the algorithmic cost is polynomial, $\mathcal{O}\left(d^{3}\right)$. ALPS is illustrated on a complex multimodal posterior distribution that arises from a seemingly-unrelated regression (SUR) model of longitudinal data from U.S. manufacturing firms.


title: Optimal Variable Clustering for High-Dimensional Matrix Valued Data

abstract: Matrix valued data has become increasingly prevalent in many applications. Most of the existing clustering methods for this type of data are tailored to the mean model and do not account for the dependence structure of the features, which can be very informative, especially in high-dimensional settings or when mean information is not available. To extract the information from the dependence structure for clustering, we propose a new latent variable model for the features arranged in matrix form, with some unknown membership matrices representing the clusters for the rows and columns. Under this model, we further propose a class of hierarchical clustering algorithms using the difference of a weighted covariance matrix as the dissimilarity measure. Theoretically, we show that under mild conditions, our algorithm attains clustering consistency in the high-dimensional setting. While this consistency result holds for our algorithm with a broad class of weighted covariance matrices, the conditions for this result depend on the choice of the weight. To investigate how the weight affects the theoretical performance of our algorithm, we establish the minimax lower bound for clustering under our latent variable model in terms of some cluster separation metric. Given these results, we identify the optimal weight in the sense that using this weight guarantees our algorithm to be minimax rate-optimal. The practical implementation of our algorithm with the optimal weight is also discussed. Simulation studies show that our algorithm performs better than existing methods in terms of the adjusted Rand index (ARI). The method is applied to a genomic dataset and yields meaningful interpretations.


title: Towards Understanding Human Functional Brain Development with Explainable Artificial Intelligence: Challenges and Perspectives

abstract: The last decades have seen significant advancements in non-invasive neuroimaging technologies that have been increasingly adopted to examine human brain development. However, these improvements have not necessarily been followed by more sophisticated data analysis measures that are able to explain the mechanisms underlying functional brain development. For example, the shift from univariate (single area in the brain) to multivariate (multiple areas in brain) analysis paradigms is of significance as it allows investigations into the interactions between different brain regions. However, despite the potential of multivariate analysis to shed light on the interactions between developing brain regions, artificial intelligence (AI) techniques applied render the analysis non-explainable. The purpose of this paper is to understand the extent to which current state-of-the-art AI techniques can inform functional brain development. In addition, a review of which AI techniques are more likely to explain their learning based on the processes of brain development as defined by developmental cognitive neuroscience (DCN) frameworks is also undertaken. This work also proposes that eXplainable AI (XAI) may provide viable methods to investigate functional brain development as hypothesised by DCN frameworks.


title: Cluster-guided Image Synthesis with Unconditional Models

abstract: Generative Adversarial Networks (GANs) are the driving force behind the state-of-the-art in image generation. Despite their ability to synthesize high-resolution photo-realistic images, generating content with on-demand conditioning of different granularity remains a challenge. This challenge is usually tackled by annotating massive datasets with the attributes of interest, a laborious task that is not always a viable option. Therefore, it is vital to introduce control into the generation process of unsupervised generative models. In this work, we focus on controllable image generation by leveraging GANs that are well-trained in an unsupervised fashion. To this end, we discover that the representation space of intermediate layers of the generator forms a number of clusters that separate the data according to semantically meaningful attributes (e.g., hair color and pose). By conditioning on the cluster assignments, the proposed method is able to control the semantic class of the generated image. Our approach enables sampling from each cluster by Implicit Maximum Likelihood Estimation (IMLE). We showcase the efficacy of our approach on faces (CelebA-HQ and FFHQ), animals (Imagenet) and objects (LSUN) using different pre-trained generative models. The results highlight the ability of our approach to condition image generation on attributes like gender, pose and hair style on faces, as well as a variety of features on different object classes.


title: TSAX is Trending

abstract: Time series mining is an important branch of data mining, as time series data is ubiquitous and has many applications in several domains. The main task in time series mining is classification. Time series representation methods play an important role in time series classification and other time series mining tasks. One of the most popular representation methods of time series data is the Symbolic Aggregate approXimation (SAX). The secret behind its popularity is its simplicity and efficiency. SAX has however one major drawback, which is its inability to represent trend information. Several methods have been proposed to enable SAX to capture trend information, but this comes at the expense of complex processing, preprocessing, or post-processing procedures. In this paper we present a new modification of SAX that we call Trending SAX (TSAX), which only adds minimal complexity to SAX, but substantially improves its performance in time series classification. This is validated experimentally on 50 datasets. The results show the superior performance of our method, as it gives a smaller classification error on 39 datasets compared with SAX.