title: Global Pointer: Novel Efficient Span-based Approach for Named Entity Recognition

abstract: Named entity recognition (NER) task aims at identifying entities from a piece of text that belong to predefined semantic types such as person, location, organization, etc. The state-of-the-art solutions for flat entities NER commonly suffer from capturing the fine-grained semantic information in underlying texts. The existing span-based approaches overcome this limitation, but the computation time is still a concern. In this work, we propose a novel span-based NER framework, namely Global Pointer (GP), that leverages the relative positions through a multiplicative attention mechanism. The ultimate goal is to enable a global view that considers the beginning and the end positions to predict the entity. To this end, we design two modules to identify the head and the tail of a given entity to enable the inconsistency between the training and inference processes. Moreover, we introduce a novel classification loss function to address the imbalance label problem. In terms of parameters, we introduce a simple but effective approximate method to reduce the training parameters. We extensively evaluate GP on various benchmark datasets. Our extensive experiments demonstrate that GP can outperform the existing solution. Moreover, the experimental results show the efficacy of the introduced loss function compared to softmax and entropy alternatives.


title: Joint Beamforming Design in DFRC Systems for Wideband Sensing and OFDM Communications

abstract: Dual-function radar-communication (DFRC) systems, which can efficiently utilize the congested spectrum and costly hardware resources by employing one common waveform for both sensing and communication (S&C), have attracted increasing attention. While the orthogonal frequency division multiplexing (OFDM) technique has been widely adopted to support high-quality communications, it also has great potentials of improving radar sensing performance and providing flexible S&C. In this paper, we propose to jointly design the dual-functional transmit signals occupying several subcarriers to realize multi-user OFDM communications and detect one moving target in the presence of clutter. Meanwhile, the signals in other frequency subcarriers can be optimized in a similar way to perform other tasks. The transmit beamforming and receive filter are jointly optimized to maximize the radar output signal-to-interference-plus-noise ratio (SINR), while satisfying the communication SINR requirement and the power budget. An majorization minimization (MM) method based algorithm is developed to solve the resulting non-convex optimization problem. Numerical results reveal the significant wideband sensing gain brought by jointly designing the transmit signals in different subcarriers, and demonstrate the advantages of our proposed scheme and the effectiveness of the developed algorithm.


title: Derivation of a macroscopic model for Brownian hard needles

abstract: We study the role of anisotropic steric interactions in a system of hard Brownian needles. Despite having no volume, non-overlapping needles exclude a volume in configuration space that influences the macroscopic evolution of the system. Starting from the stochastic particle system, we use the method of matched asymptotic expansions and conformal mapping to systematically derive a nonlinear nonlocal partial differential equation for the evolution of the population density in position and orientation. We consider the regime of high rotational diffusion, resulting in an equation for the spatial density that allows us to compare the effective excluded volume of a hard-needles system with that of a hard-spheres system. We further consider spatially homogeneous solutions and find an isotropic to nematic transition as density increases, consistent with Onsager's theory.


title: Dark path holonomic qudit computation

abstract: Non-adiabatic holonomic quantum computation is a method used to implement high-speed quantum gates with non-Abelian geometric phases associated with paths in state space. Due to their noise tolerance, these phases can be used to construct error resilient quantum gates. We extend the holonomic dark path qubit scheme in [M.-Z. Ai {\t et al.}, Fundam. Res. {\bf 2}, 661 (2022)] to qudits. Specifically, we demonstrate one- and two-qudit universality by using the dark path technique. Explicit qutrit ($d=3$) gates are demonstrated and the scaling of the number of loops with the dimension $d$ is addressed. This scaling is linear and we show how any diagonal qudit gate can be implemented efficiently in any dimension.


title: Multi-Axis Control of a Qubit in the Presence of Unknown Non-Markovian Quantum Noise

abstract: In this paper, we consider the problem of open-loop control of a qubit that is coupled to an unknown fully quantum non-Markovian noise (either bosonic or fermionic). A graybox model that is empirically obtained from measurement data is employed to approximately represent the unknown quantum noise. The estimated model is then used to calculate the open-loop control pulses under constraints on the pulse amplitude and timing. For the control pulse optimization, we explore the use of gradient descent and genetic optimization methods. We consider the effect of finite sampling on estimating expectation values of observables and show results for single- and multi-axis control of a qubit.


title: Sinusoidal Sensitivity Calculation for Line Segment Geometries

abstract: Purpose: Provide a closed-form solution to the sinusoidal coil sensitivity model proposed by Kern et al. Solution allows for the precise computations of varied, simulated bias fields which can be directly applied onto raw intensity datasets. Methods: Fourier distribution theory and standard integration techniques were used to calculate the Fourier transform of measured magnetic field produced from line segment sources. Results: A $L^1_{\rm loc}(\mathbb{R}^3)$ function is derived in full generality for arbitrary line segment geometries. Sampling criteria and equivalence to the original sinusoidal model are discussed. Lastly a CUDA accelerated implementation $\texttt{biasgen}$ is provided for on-demand sensitivity and bias generation. Conclusion: Given the modeling flexibility of the simulated procedure, practitioners will now have access to a more diverse ecosystem of simulated datasets which may be used to quantitatively compare prospective debiasing methods.


title: Probing critical behavior of long-range transverse-field Ising model through quantum Kibble-Zurek mechanism

abstract: The trapped ion quantum simulator has demonstrated qualitative properties of different physical models for up to tens of ions. In particular, a linear ion chain naturally hosts long-range Ising interactions under the laser driving, which has been used for various phenomena such as quantum phase transition, localization, thermalization and information propagation. For near-term practical usage, a central task is to find more quantitative applications of the noisy quantum simulators that are robust to small errors in the parameters. Here we report the quantum simulation of a long-range transverse-field Ising model using up to 61 ions and probe the critical behavior of its quantum phase transition through the Kibble-Zurek mechanism. By calibrating and verifying the coupling coefficients, we realize the same model for increasing ion numbers, so as to extract a critical exponent free of the finite size effect. For ferromagnetic interaction, our experimental result agrees well with the previous numerical predictions. As for the anti-ferromagnetic case, signals are too weak to fit a critical exponent due to the frustration in the interaction, but still consistent with the theory.


title: The reaction $\pi N \to \omega N$ in a dynamical coupled-channel approach

abstract: A refined investigation on light flavor meson-baryon scatterings is performed using a dynamical coupled-channel approach, the J\"ulich-Bonn model, that respects unitartiy and analyticity constraints. The channel space of $\pi N$, $\pi \Delta$, $\sigma N$, $\rho N$, $\eta N$, $K \Lambda$ and $K \Sigma$ is extended by adding the $\omega N$ final state. The spectra of $N^*$ and $\Delta$ resonances are extracted in terms of complex poles of the scattering amplitudes, based on the result of a global fit to a worldwide collection of data, in the energy region from the $\pi N$ threshold to center-of-mass energy $z=2.3$ GeV. A negative value of the $\omega N$ elastic spin-averaged scattering length is extracted, questioning the existence of bound states of the $\omega$ meson in the nuclear matter.


title: Exploring Resolution and Degradation Clues as Self-supervised Signal for Low Quality Object Detection

abstract: Image restoration algorithms such as super resolution (SR) are indispensable pre-processing modules for object detection in low quality images. Most of these algorithms assume the degradation is fixed and known a priori. However, in practical, either the real degradation or optimal up-sampling ratio rate is unknown or differs from assumption, leading to a deteriorating performance for both the pre-processing module and the consequent high-level task such as object detection. Here, we propose a novel self-supervised framework to detect objects in degraded low resolution images. We utilizes the downsampling degradation as a kind of transformation for self-supervised signals to explore the equivariant representation against various resolutions and other degradation conditions. The Auto Encoding Resolution in Self-supervision (AERIS) framework could further take the advantage of advanced SR architectures with an arbitrary resolution restoring decoder to reconstruct the original correspondence from the degraded input image. Both the representation learning and object detection are optimized jointly in an end-to-end training fashion. The generic AERIS framework could be implemented on various mainstream object detection architectures with different backbones. The extensive experiments show that our methods has achieved superior performance compared with existing methods when facing variant degradation situations. Code would be released at https://github.com/cuiziteng/ECCV_AERIS.


title: Enhancing the Robustness via Adversarial Learning and Joint Spatial-Temporal Embeddings in Traffic Forecasting

abstract: Traffic forecasting is an essential problem in urban planning and computing. The complex dynamic spatial-temporal dependencies among traffic objects (e.g., sensors and road segments) have been calling for highly flexible models; unfortunately, sophisticated models may suffer from poor robustness especially in capturing the trend of the time series (1st-order derivatives with time), leading to unrealistic forecasts. To address the challenge of balancing dynamics and robustness, we propose TrendGCN, a new scheme that extends the flexibility of GCNs and the distribution-preserving capacity of generative and adversarial loss for handling sequential data with inherent statistical correlations. On the one hand, our model simultaneously incorporates spatial (node-wise) embeddings and temporal (time-wise) embeddings to account for heterogeneous space-and-time convolutions; on the other hand, it uses GAN structure to systematically evaluate statistical consistencies between the real and the predicted time series in terms of both the temporal trending and the complex spatial-temporal dependencies. Compared with traditional approaches that handle step-wise predictive errors independently, our approach can produce more realistic and robust forecasts. Experiments on six benchmark traffic forecasting datasets and theoretical analysis both demonstrate the superiority and the state-of-the-art performance of TrendGCN. Source code is available at https://github.com/juyongjiang/TrendGCN.


title: Immersions of punctured 4-manifolds

abstract: Motivated by applications to pulling back quantum cellular automata from one manifold to another, we study the existence of immersions between certain smooth 4-manifolds. We show that they lead to a very interesting partial order on closed 4-manifolds.


title: Compact spaces associated to separable Banach lattices

abstract: We study the class of compact spaces that appear as structure spaces of separable Banach lattices. In other words, we analyze what $C(K)$ spaces appear as principal ideals of separable Banach lattices. Among other things, it is shown that every such compactum $K$ admits a strictly positive regular Borel measure of countable type that is analytic, and in the nonmetrizable case these compacta are saturated with copies of $\beta\mathbb N$. Some natural questions about this class are left open.


title: Tailoring to the Tails: Risk Measures for Fine-Grained Tail Sensitivity

abstract: Expected risk minimization (ERM) is at the core of many machine learning systems. This means that the risk inherent in a loss distribution is summarized using a single number - its average. In this paper, we propose a general approach to construct risk measures which exhibit a desired tail sensitivity and may replace the expectation operator in ERM. Our method relies on the specification of a reference distribution with a desired tail behaviour, which is in a one-to-one correspondence to a coherent upper probability. Any risk measure, which is compatible with this upper probability, displays a tail sensitivity which is finely tuned to the reference distribution. As a concrete example, we focus on divergence risk measures based on f-divergence ambiguity sets, which are a widespread tool used to foster distributional robustness of machine learning systems. For instance, we show how ambiguity sets based on the Kullback-Leibler divergence are intricately tied to the class of subexponential random variables. We elaborate the connection of divergence risk measures and rearrangement invariant Banach norms.


title: Large vocabulary speech recognition for languages of Africa: multilingual modeling and self-supervised learning

abstract: Almost none of the 2,000+ languages spoken in Africa have widely available automatic speech recognition systems, and the required data is also only available for a few languages. We have experimented with two techniques which may provide pathways to large vocabulary speech recognition for African languages: multilingual modeling and self-supervised learning. We gathered available open source data and collected data for 15 languages, and trained experimental models using these techniques. Our results show that pooling the small amounts of data available in multilingual end-to-end models, and pre-training on unsupervised data can help improve speech recognition quality for many African languages.


title: On the Diophantine equation $U_n-b^m = c$

abstract: Let $(U_n)_{n\in \mathbb{N}}$ be a fixed linear recurrence sequence defined over the integers (with some technical restrictions). We prove that there exist effectively computable constants $B$ and $N_0$ such that for any $b,c\in \mathbb{Z}$ with $b> B$ the equation $U_n - b^m = c$ has at most two distinct solutions $(n,m)\in \mathbb{N}^2$ with $n\geq N_0$ and $m\geq 1$. Moreover, we apply our result to the special case of Tribonacci numbers given by $T_1= T_2=1$, $T_3=2$ and $T_{n}=T_{n-1}+T_{n-2}+T_{n-3}$ for $n\geq 4$. By means of the LLL-algorithm and continued fraction reduction we are able to prove $N_0=1.1\cdot 10^{37}$ and $B=e^{438}$. The corresponding reduction algorithm is implemented in Sage.


title: The first minutes of a binary-driven hypernova

abstract: We simulate the first minutes of the evolution of a binary-driven hypernova (BdHN) event, with a special focus on the associated accretion processes of supernova (SN) ejecta onto the newborn neutron star ($\nu$NS) and the NS companion. We calculate the rotational evolution of the $\nu$NS and the NS under the torques exerted by the accreted matter and the magnetic field. We take into account general relativistic effects and use realistic hypercritical accretion rates obtained from three-dimensional smoothed-particle-hydrodynamics (SPH) numerical simulations of the BdHN for a variety of orbital periods. We show that the rotation power of the $\nu$NS has a unique double-peak structure while that of the NS has a single peak. These peaks are of comparable intensity and can occur very close in time or even simultaneously depending on the orbital period and the initial angular momentum of the stars. We outline the consequences of the above features in the early emission and their consequent observation in long gamma-ray bursts (GRBs).


title: Activity Detection in Distributed MIMO: Distributed AMP via Likelihood Ratio Fusion

abstract: We develop a new algorithm for activity detection for grant-free multiple access in distributed multiple-input multiple-output (MIMO). The algorithm is a distributed version of the approximate message passing (AMP) based on a soft combination of likelihood ratios computed independently at multiple access points. The underpinning theoretical basis of our algorithm is a new observation that we made about the state evolution in the AMP. Specifically, with a minimum mean-square error denoiser, the state maintains a block-diagonal structure whenever the covariance matrices of the signals have such a structure. We show by numerical examples that the algorithm outperforms competing schemes from the literature.


title: On Hermitian manifolds with Bismut-Strominger parallel torsion

abstract: In this article, we study Hermitian manifolds whose Bismut-Strominger connection has parallel torsion tensor, which will be called {\em Bismut torsion parallel manifolds,} or {\em BTP} manifolds for short. We obtain a necessary and sufficient condition characterizing this class in terms of the curvature tensor. In particular, Bismut flat or Bismut K\"ahler-like manifolds are {\em BTP} manifolds, known by our earlier results. We also obtain structural results for non-balanced {\em BTP} manifolds, and classification theorems for non-balanced {\em BTP} threefolds and balanced {\em BTP} threefolds.


title: Pushing Least-Squares Deconvolution to the next level: application to binary stars

abstract: Eclipsing, spectroscopic double-lined (SB2) binaries remain to be the prime source of precise and accurate fundamental properties of stars. Furthermore, high-cadence spectroscopic observations of the eclipse phases allow us to resolve the Rossiter-McLaughlin effect whose modelling offers the means to probe spin-orbit misalignment in binaries. In this study, we develop the LSDBinary algorithm that is capable of working with both in-eclipse and out-of-eclipse spectra of SB2 binaries as input and delivers the LSD profiles, LSD-based model spectra, and precise RVs of both binary components as output. We offer an option to account for the Rossiter-McLaughlin effect in the calculation of the initial guess LSD profiles and components' flux ratio such that the effect can be modelled within the algorithm itself. We provide an extensive test of the LSDBinary software package on simulated spectra of artificial binaries. We study the effects of signal-to-noise-ratio of input spectra, resolving power of the instrument, uncertain atmospheric parameters of stars, and orbital properties of the binary system on the resulting LSD profiles and RVs measured from them. We find that atmospheric parameters have negligible effect on the shape of the computed LSD profiles while affecting mostly their global scaling. Our results are barely sensitive to signal-to-noise ratio of the input spectra provided they contain sufficient number of spectral lines, such as in A-type stars and later. Finally, the orbital inclination angle and components' radii ratio are found to have the largest effect on the shapes of the LSD profiles and RV curves extracted from them. The LSDBinary algorithm is specifically developed to perform detailed spectroscopic studies of eclipsing SB2 systems whose orbital configuration and components' atmospheric parameters are estimated by other means.


title: On the Strict Picard Spectrum of Commutative Ring Spectra

abstract: We compute the connective spectra of maps from $\mathbb{Z}$ to the Picard spectra of the spherical Witt vectors associated with perfect rings of characteristic $p$. As an application, we determine the connective spectrum of maps from $\mathbb{Z}$ to the Picard spectrum of the sphere spectrum.


title: Analytical disk-cylinder interaction potential laws for the computational modeling of adhesive, deformable (nano)fibers

abstract: The analysis of complex fibrous systems or materials on the micro- and nanoscale, which have a high practical relevance for many technical or biological systems, requires accurate analytical descriptions of the adhesive and repulsive forces acting on the fiber surfaces. While such analytical expressions are generally needed both for theoretical studies and for computer-based simulations, the latter motivates us here to derive disk-cylinder interaction potential laws that are valid for arbitrary mutual orientations in the decisive regime of small surface separations. The chosen type of fundamental point-pair interaction follows the simple Lennard-Jones model with inverse power laws for both the adhesive van der Waals part and the steric, repulsive part. We present three different solutions, ranging from highest accuracy to the best trade-off between simplicity of the expression and sufficient accuracy for our intended use. The validity of simplifying approximations and the accuracy of the derived potential laws is thoroughly analyzed, using both numerical and analytical reference solutions for specific interaction cases. Most importantly, the correct asymptotic scaling behavior in the decisive regime of small separations is achieved, and also the theoretically predicted $(1\!/\!\sin\!\alpha)$-angle dependence (for non-parallel cylinders) is obtained by the proposed analytical solutions. As we show in the outlook to our current research, the derived analytical disk-cylinder interaction potential laws may be used to formulate highly efficient computational models for the interaction of arbitrarily curved fibers, such that the disk represents the cross-section of the first and the cylinder a local approximation to the shape of the second fiber.


title: PGX: A Multi-level GNN Explanation Framework Based on Separate Knowledge Distillation Processes

abstract: Graph Neural Networks (GNNs) are widely adopted in advanced AI systems due to their capability of representation learning on graph data. Even though GNN explanation is crucial to increase user trust in the systems, it is challenging due to the complexity of GNN execution. Lately, many works have been proposed to address some of the issues in GNN explanation. However, they lack generalization capability or suffer from computational burden when the size of graphs is enormous. To address these challenges, we propose a multi-level GNN explanation framework based on an observation that GNN is a multimodal learning process of multiple components in graph data. The complexity of the original problem is relaxed by breaking into multiple sub-parts represented as a hierarchical structure. The top-level explanation aims at specifying the contribution of each component to the model execution and predictions, while fine-grained levels focus on feature attribution and graph structure attribution analysis based on knowledge distillation. Student models are trained in standalone modes and are responsible for capturing different teacher behaviors, later used for particular component interpretation. Besides, we also aim for personalized explanations as the framework can generate different results based on user preferences. Finally, extensive experiments demonstrate the effectiveness and fidelity of our proposed approach.


title: On the weak second-order optimality condition for nonlinear semidefinite and second-order cone programming

abstract: Second-order necessary optimality conditions for nonlinear conic programming problems that depend on a single Lagrange multiplier are usually built under nondegeneracy and strict complementarity. In this paper we establish a condition of such type for two classes of nonlinear conic problems, namely semidefinite and second-order cone programming, assuming Robinson's constraint qualification and a weak constant rank-type property which are, together, strictly weaker than nondegeneracy. Our approach is done via a penalty-based strategy, which is aimed at providing strong global convergence results for first- and second-order algorithms. Since we are not assuming strict complementarity, the critical cone does not reduce to a subspace, thus, the second-order condition we arrive at is defined in terms of the lineality space of the critical cone. In the case of nonlinear programming, this condition reduces to the standard second-order condition widely used as second-order stationarity measure in the algorithmic practice.


title: Universal exploration dynamics of random walks

abstract: The territory explored by a random walk is a key property that may be quantified by the number of distinct sites that the random walk visits up to a given time. The extent of this spatial exploration characterizes many important physical, chemical, and ecological phenomena. In spite of its fundamental interest and wide utility, the number of visited sites gives only an incomplete picture of this exploration. In this work, we introduce a more fundamental quantity, the elapsed time $\tau_n$ between visits to the $n^{\rm th}$ and the $(n+1)^{\rm st}$ distinct sites, from which the full dynamics about the visitation statistics can be obtained. To determine the distribution of these inter-visit times $\tau_n$, we develop a theoretical approach that relies on a mapping with a trapping problem, in which, in contrast to previously studied situations, the spatial distribution of traps is continuously updated by the random walk itself. Despite the geometrical complexity of the territory explored by a random walk (typically aspherical, as well as containing holes and islands at all scales), we find that the distribution of the $\tau_n$ can be accounted for by simple analytical expressions. Processes as varied as regular diffusion, anomalous diffusion, and diffusion in disordered media and fractals, fall into the same universality classes for the temporal history of distinct sites visited. We confirm our theoretical predictions by Monte Carlo and exact enumeration methods. We also determine additional basic exploration observables, such as the perimeter of the visited domain or the number of islands of unvisited sites enclosed within this domain, thereby illustrating the generality of our approach. Because of their fundamental character and their universality, these inter-visit times represent a promising tool to unravel many more aspects of the exploration dynamics of random walks.


title: Cohort comfort models -- Using occupants' similarity to predict personal thermal preference with less data

abstract: We introduce Cohort Comfort Models, a new framework for predicting how new occupants would perceive their thermal environment. Cohort Comfort Models leverage historical data collected from a sample population, who have some underlying preference similarity, to predict thermal preference responses of new occupants. Our framework is capable of exploiting available background information such as physical characteristics and one-time on-boarding surveys (satisfaction with life scale, highly sensitive person scale, the Big Five personality traits) from the new occupant as well as physiological and environmental sensor measurements paired with thermal preference responses. We implemented our framework in two publicly available datasets containing longitudinal data from 55 people, comprising more than 6,000 individual thermal comfort surveys. We observed that, a Cohort Comfort Model that uses background information provided very little change in thermal preference prediction performance but uses none historical data. On the other hand, for half and one third of each dataset occupant population, using Cohort Comfort Models, with less historical data from target occupants, Cohort Comfort Models increased their thermal preference prediction by 8~\% and 5~\% on average, and up to 36~\% and 46~\% for some occupants, when compared to general-purpose models trained on the whole population of occupants. The framework is presented in a data and site agnostic manner, with its different components easily tailored to the data availability of the occupants and the buildings. Cohort Comfort Models can be an important step towards personalization without the need of developing a personalized model for each new occupant.


title: Instance As Identity: A Generic Online Paradigm for Video Instance Segmentation

abstract: Modeling temporal information for both detection and tracking in a unified framework has been proved a promising solution to video instance segmentation (VIS). However, how to effectively incorporate the temporal information into an online model remains an open problem. In this work, we propose a new online VIS paradigm named Instance As Identity (IAI), which models temporal information for both detection and tracking in an efficient way. In detail, IAI employs a novel identification module to predict identification number for tracking instances explicitly. For passing temporal information cross frame, IAI utilizes an association module which combines current features and past embeddings. Notably, IAI can be integrated with different image models. We conduct extensive experiments on three VIS benchmarks. IAI outperforms all the online competitors on YouTube-VIS-2019 (ResNet-101 43.7 mAP) and YouTube-VIS-2021 (ResNet-50 38.0 mAP). Surprisingly, on the more challenging OVIS, IAI achieves SOTA performance (20.6 mAP). Code is available at https://github.com/zfonemore/IAI


title: Double Mode Cepheids from the Zwicky Transient Facility Survey

abstract: Multi-mode Cepheids pulsate simultaneously in more than one mode of oscillation. They provide an independent means to test stellar models and pulsation theories. They can also be used to derive metallicities. In recent years, the number of known multi-mode Cepheids has increased dramatically with the discovery of a large number of Galactic double-mode Cepheids. To date, 209 double-mode Cepheids have been detected in the Galactic bulge and disk, mostly based on the Optical Gravitational Lensing Experiment's (OGLE) catalog. In this paper, we conduct a comprehensive search for double-mode Cepheids in the northern sky based on Zwicky Transient Facility Data Release 5. We found 72 such objects in the Milky Way. The periods of the 30 sample objects already included in the OGLE catalog show excellent agreement with the OGLE periods. The period ratios of our new Cepheids are consistent with those of known double-mode Cepheids, as evidenced by their loci in the so-called `Petersen diagram'. Compared with OGLE, the completeness of our double-mode Cepheid sample is around 71\%. The much improved temporal sampling of the Zwicky Transient Facility offers significant scope to find more double-mode Cepheids, especially at the distribution's short-period end.


title: Validation of EEG forward modeling approaches in the presence of anisotropy in the source space

abstract: The quality of the inverse approach in electroencephalography (EEG) source analysis is - among other things - depending on the accuracy of the forward modeling approach, i.e., the simulation of the electric potential for a known dipole source in the brain. Here, we use multilayer sphere modeling scenarios to investigate the performance of three different finite element method (FEM) based EEG forward approaches - subtraction, Venant and partial integration - in the presence of tissue conductivity anisotropy in the source space. In our studies, the effect of anisotropy on the potential is related to model errors when ignoring anisotropy and to numerical errors, convergence behavior and computational speed of the different FEM approaches. Three different source space anisotropy models that best represent adult, child and premature baby volume conduction scenarios, are used. Major findings of the study include (1) source space conductivity anisotropy has a significant effect on electric potential computation: The effect increases with increasing anisotropy ratio; (2) with numerical errors far below anisotropy effects, all three FEM approaches are able to model source space anisotropy accordingly, with the Venant approach offering the best compromise between accuracy and computational speed; (3) FE meshes have to be fine enough in the subdomain between the source and the sensors that capture its main activity. We conclude that, especially for the analysis of cortical development, but also for more general applications using EEG source analysis techniques, source space conductivity anisotropy should be modeled and the FEM Venant approach is an appropriate method.


title: Conflict-free joint sampling for preference satisfaction through quantum interference

abstract: Collective decision-making is vital for recent information and communications technologies. In our previous research, we mathematically derived conflict-free joint decision-making that optimally satisfies players' probabilistic preference profiles. However, two problems exist regarding the optimal joint decision-making method. First, as the number of choices increases, the computational cost of calculating the optimal joint selection probability matrix explodes. Second, to derive the optimal joint selection probability matrix, all players must disclose their probabilistic preferences. Now, it is noteworthy that explicit calculation of the joint probability distribution is not necessarily needed; what is necessary for collective decisions is sampling. This study examines several sampling methods that converge to heuristic joint selection probability matrices that satisfy players' preferences. We show that they can significantly reduce the above problems of computational cost and confidentiality. We analyze the probability distribution each of the sampling methods converges to, as well as the computational cost required and the confidentiality secured. In particular, we introduce two conflict-free joint sampling methods through quantum interference of photons. The first system allows the players to hide their choices while satisfying the players' preferences almost perfectly when they have the same preferences. The second system, where the physical nature of light replaces the expensive computational cost, also conceals their choices under the assumption that they have a trusted third party. This paper has been published in Phys. Rev. Applied 18, 064018 (2022) (DOI: 10.1103/PhysRevApplied.18.064018).

categories: quant-ph cs.AI physics.app-ph physics.optics


title: Neural Network Verification using Residual Reasoning

abstract: With the increasing integration of neural networks as components in mission-critical systems, there is an increasing need to ensure that they satisfy various safety and liveness requirements. In recent years, numerous sound and complete verification methods have been proposed towards that end, but these typically suffer from severe scalability limitations. Recent work has proposed enhancing such verification techniques with abstraction-refinement capabilities, which have been shown to boost scalability: instead of verifying a large and complex network, the verifier constructs and then verifies a much smaller network, whose correctness implies the correctness of the original network. A shortcoming of such a scheme is that if verifying the smaller network fails, the verifier needs to perform a refinement step that increases the size of the network being verified, and then start verifying the new network from scratch - effectively "wasting" its earlier work on verifying the smaller network. In this paper, we present an enhancement to abstraction-based verification of neural networks, by using residual reasoning: the process of utilizing information acquired when verifying an abstract network, in order to expedite the verification of a refined network. In essence, the method allows the verifier to store information about parts of the search space in which the refined network is guaranteed to behave correctly, and allows it to focus on areas where bugs might be discovered. We implemented our approach as an extension to the Marabou verifier, and obtained promising results.


title: Deep Feature Learning for Medical Acoustics

abstract: The purpose of this paper is to compare different learnable frontends in medical acoustics tasks. A framework has been implemented to classify human respiratory sounds and heartbeats in two categories, i.e. healthy or affected by pathologies. After obtaining two suitable datasets, we proceeded to classify the sounds using two learnable state-of-art frontends -- LEAF and nnAudio -- plus a non-learnable baseline frontend, i.e. Mel-filterbanks. The computed features are then fed into two different CNN models, namely VGG16 and EfficientNet. The frontends are carefully benchmarked in terms of the number of parameters, computational resources, and effectiveness. This work demonstrates how the integration of learnable frontends in neural audio classification systems may improve performance, especially in the field of medical acoustics. However, the usage of such frameworks makes the needed amount of data even larger. Consequently, they are useful if the amount of data available for training is adequately large to assist the feature learning process.


title: Optimistic Gradient Descent Ascent in Zero-Sum and General-Sum Bilinear Games

abstract: We study the convergence of Optimistic Gradient Descent Ascent in unconstrained bilinear games. In a first part, we consider the zero-sum case and extend previous results by Daskalakis et al. in 2018, Liang and Stokes in 2019, and others: we prove, for any payoff matrix, the exponential convergence of OGDA to a saddle point and also provide a new, optimal, geometric ratio for the convergence. We also characterize the step sizes inducing convergence, and are able to deduce the optimal step size for the speed of convergence. In a second part, we introduce OGDA for general-sum bilinear games: we show that in an interesting class of games, either OGDA converges exponentially fast to a Nash equilibrium, or the payoffs for both players converge exponentially fast to $+\infty$ (which might be interpreted as endogenous emergence of coordination, or cooperation, among players). We also give sufficient conditions for convergence of OGDA to a Nash equilibrium. These conditions are used to increase the speed of convergence of a min-max problem involving a matrix $A$, by introducing a general-sum game using the Moore-Penrose inverse matrix of $A$. This shows for the first time, at our knowledge, that general-sum games can be used to optimally improve algorithms designed for min-max problems. We finally illustrate our results on simple examples of Generative Adversarial Networks.


title: QCD in the cores of neutron stars

abstract: I discuss why state-of-the art perturbative QCD calculations of the equation of state at large chemical potential that are reliable at asymptotically high densities constrain the same equation of state at neutron-star densities. I describe how these theoretical calculations affect the EOS at lower density. I argue that the ab-initio calculations in QCD offer significant information about the equation of state of the neutron-star matter, which is complementary to the current astrophysical observations.


title: A Fixpoint Characterization of Three-Valued Disjunctive Hybrid MKNF Knowledge Bases

abstract: The logic of hybrid MKNF (minimal knowledge and negation as failure) is a powerful knowledge representation language that elegantly pairs ASP (answer set programming) with ontologies. Disjunctive rules are a desirable extension to normal rule-based reasoning and typically semantic frameworks designed for normal knowledge bases need substantial restructuring to support disjunctive rules. Alternatively, one may lift characterizations of normal rules to support disjunctive rules by inducing a collection of normal knowledge bases, each with the same body and a single atom in its head. In this work, we refer to a set of such normal knowledge bases as a head-cut of a disjunctive knowledge base. The question arises as to whether the semantics of disjunctive hybrid MKNF knowledge bases can be characterized using fixpoint constructions with head-cuts. Earlier, we have shown that head-cuts can be paired with fixpoint operators to capture the two-valued MKNF models of disjunctive hybrid MKNF knowledge bases. Three-valued semantics extends two-valued semantics with the ability to express partial information. In this work, we present a fixpoint construction that leverages head-cuts using an operator that iteratively captures three-valued models of hybrid MKNF knowledge bases with disjunctive rules. This characterization also captures partial stable models of disjunctive logic programs since a program can be expressed as a disjunctive hybrid MKNF knowledge base with an empty ontology. We elaborate on a relationship between this characterization and approximators in AFT (approximation fixpoint theory) for normal hybrid MKNF knowledge bases.


title: Jumping Evaluation of Nested Regular Path Queries

abstract: Nested regular path queries are used for querying graph databases and RDF triple stores. We propose a new algorithm for evaluating nested regular path queries on a graph from a set of start nodes in combined linear time. We show that this complexity upper bound can be reduced by making it dependent on the size of the query's top-down needed subgraph, a notion that we introduce. For many queries in practice, the top-down needed subgraph is way smaller than the whole graph. Our algorithm is based on a novel compilation schema from nested regular path queries to monadic datalog queries. Its complexity upper bound follows from known properties of top-down datalog evaluation. As an application, we show that our algorithm permits to reformulate in simple terms a variant of a very efficient automata-based algorithm proposed by Maneth and Nguyen that evaluates navigational path queries in datatrees based on indexes and jumping. Moreover, it overcomes some limitations of Maneth and Nguyen's: it is not bound to trees and applies to graphs; it is not limited to forward navigational XPath but can treat any nested regular path query and it can be implemented efficiently without any dedicated techniques, by using any efficient datalog evaluator such as LogicBlox.


title: Tree-Like Justification Systems are Consistent

abstract: Justification theory is an abstract unifying formalism that captures semantics of various non-monotonic logics. One intriguing problem that has received significant attention is the consistency problem: under which conditions are justifications for a fact and justifications for its negation suitably related. Two variants of justification theory exist: one in which justifications are trees and one in which they are graphs. In this work we resolve the consistency problem once and for all for the tree-like setting by showing that all reasonable tree-like justification systems are consistent.


title: A Preliminary Data-driven Analysis of Common Errors Encountered by Novice SPARC Programmers

abstract: Answer Set Programming (ASP), a modern development of Logic Programming, enables a natural integration of Computing with STEM subjects. This integration addresses a widely acknowledged challenge in K-12 education, and early empirical results on ASP-based integration are promising. Although ASP is considered a simple language when compared with imperative programming languages, programming errors can still be a significant barrier for students. This is particularly true for K-12 students who are novice users of ASP. Categorizing errors and measuring their difficulty has yielded insights into imperative languages like Java. However, little is known about the types and difficulty of errors encountered by K-12 students using ASP. To address this, we collected high school student programs submitted during a 4-session seminar teaching an ASP language known as SPARC. From error messages in this dataset, we identify a collection of error classes, and measure how frequently each class occurs and how difficult it is to resolve.


title: On Model Reconciliation: How to Reconcile When Robot Does not Know Human's Model?

abstract: The Model Reconciliation Problem (MRP) was introduced to address issues in explainable AI planning. A solution to a MRP is an explanation for the differences between the models of the human and the planning agent (robot). Most approaches to solving MRPs assume that the robot, who needs to provide explanations, knows the human model. This assumption is not always realistic in several situations (e.g., the human might decide to update her model and the robot is unaware of the updates). In this paper, we propose a dialog-based approach for computing explanations of MRPs under the assumptions that (i) the robot does not know the human model; (ii) the human and the robot share the set of predicates of the planning domain and their exchanges are about action descriptions and fluents' values; (iii) communication between the parties is perfect; and (iv) the parties are truthful. A solution of a MRP is computed through a dialog, defined as a sequence of rounds of exchanges, between the robot and the human. In each round, the robot sends a potential explanation, called proposal, to the human who replies with her evaluation of the proposal, called response. We develop algorithms for computing proposals by the robot and responses by the human and implement these algorithms in a system that combines imperative means with answer set programming using the multi-shot feature of clingo.


title: An Iterative Fixpoint Semantics for MKNF Hybrid Knowledge Bases with Function Symbols

abstract: Hybrid Knowledge Bases based on Lifschitz's logic of Minimal Knowledge with Negation as Failure are a successful approach to combine the expressivity of Description Logics and Logic Programming in a single language. Their syntax, defined by Motik and Rosati, disallows function symbols. In order to define a well-founded semantics for MKNF HKBs, Knorr et al. define a partition of the modal atoms occurring in it, called the alternating fixpoint partition. In this paper, we propose an iterated fixpoint semantics for HKBs with function symbols. We prove that our semantics extends Knorr et al.'s, in that, for a function-free HKBs, it coincides with its alternating fixpoint partition. The proposed semantics lends itself well to a probabilistic extension with a distribution semantic approach, which is the subject of future work.


title: A Gaze into the Internal Logic of Graph Neural Networks, with Logic

abstract: Graph Neural Networks share with Logic Programming several key relational inference mechanisms. The datasets on which they are trained and evaluated can be seen as database facts containing ground terms. This makes possible modeling their inference mechanisms with equivalent logic programs, to better understand not just how they propagate information between the entities involved in the machine learning process but also to infer limits on what can be learned from a given dataset and how well that might generalize to unseen test data. This leads us to the key idea of this paper: modeling with the help of a logic program the information flows involved in learning to infer from the link structure of a graph and the information content of its nodes properties of new nodes, given their known connections to nodes with possibly similar properties. The problem is known as graph node property prediction and our approach will consist in emulating with help of a Prolog program the key information propagation steps of a Graph Neural Network's training and inference stages. We test our a approach on the ogbn-arxiv node property inference benchmark. To infer class labels for nodes representing papers in a citation network, we distill the dependency trees of the text associated to each node into directed acyclic graphs that we encode as ground Prolog terms. Together with the set of their references to other papers, they become facts in a database on which we reason with help of a Prolog program that mimics the information propagation in graph neural networks predicting node properties. In the process, we invent ground term similarity relations that help infer labels in the test set by propagating node properties from similar nodes in the training set and we evaluate their effectiveness in comparison with that of the graph's link structure. Finally, we implement explanation generators that unveil performance upper bounds inherent to the dataset. As a practical outcome, we obtain a logic program, that, when seen as machine learning algorithm, performs close to the state of the art on the node property prediction benchmark.


title: Knowledge Authoring with Factual English

abstract: Knowledge representation and reasoning (KRR) systems represent knowledge as collections of facts and rules. Like databases, KRR systems contain information about domains of human activities like industrial enterprises, science, and business. KRRs can represent complex concepts and relations, and they can query and manipulate information in sophisticated ways. Unfortunately, the KRR technology has been hindered by the fact that specifying the requisite knowledge requires skills that most domain experts do not have, and professional knowledge engineers are hard to find. One solution could be to extract knowledge from English text, and a number of works have attempted to do so (OpenSesame, Google's Sling, etc.). Unfortunately, at present, extraction of logical facts from unrestricted natural language is still too inaccurate to be used for reasoning, while restricting the grammar of the language (so-called controlled natural language, or CNL) is hard for the users to learn and use. Nevertheless, some recent CNL-based approaches, such as the Knowledge Authoring Logic Machine (KALM), have shown to have very high accuracy compared to others, and a natural question is to what extent the CNL restrictions can be lifted. In this paper, we address this issue by transplanting the KALM framework to a neural natural language parser, mStanza. Here we limit our attention to authoring facts and queries and therefore our focus is what we call factual English statements. Authoring other types of knowledge, such as rules, will be considered in our followup work. As it turns out, neural network based parsers have problems of their own and the mistakes they make range from part-of-speech tagging to lemmatization to dependency errors. We present a number of techniques for combating these problems and test the new system, KALMFL (i.e., KALM for factual language), on a number of benchmarks, which show KALMFL achieves correctness in excess of 95%.


title: A Model-Oriented Approach for Lifting Symmetries in Answer Set Programming

abstract: When solving combinatorial problems, pruning symmetric solution candidates from the search space is essential. Most of the existing approaches are instance-specific and focus on the automatic computation of Symmetry Breaking Constraints (SBCs) for each given problem instance. However, the application of such approaches to large-scale instances or advanced problem encodings might be problematic since the computed SBCs are propositional and, therefore, can neither be meaningfully interpreted nor transferred to other instances. As a result, a time-consuming recomputation of SBCs must be done before every invocation of a solver. To overcome these limitations, we introduce a new model-oriented approach for Answer Set Programming that lifts the SBCs of small problem instances into a set of interpretable first-order constraints using a form of machine learning called Inductive Logic Programming. After targeting simple combinatorial problems, we aim to extend our method to be applied also for advanced decision and optimization problems.


title: Tools and Methodologies for Verifying Answer Set Programs

abstract: Answer Set Programming (ASP) is a powerful declarative programming paradigm commonly used for solving challenging search and optimization problems. The modeling languages of ASP are supported by sophisticated solving algorithms (solvers) that make the solution search efficient while enabling the programmer to model the problem at a high level of abstraction. As an approach to Knowledge Representation and Reasoning, ASP benefits from its simplicity, conciseness and rigorously defined semantics. These characteristics make ASP a straightforward way to develop formally verifiable programs. In the context of artificial intelligence (AI), the clarity of ASP programs lends itself to the construction of explainable, trustworthy AI. In support of these goals, my research is concerned with extending the theory and tools supporting the verification of ASP progams.


title: An ASP Framework for Efficient Urban Traffic Optimization

abstract: Avoiding congestion and controlling traffic in urban scenarios is becoming nowadays of paramount importance due to the rapid growth of our cities' population and vehicles. The effective control of urban traffic as a means to mitigate congestion can be beneficial in an economic, environmental and health way. In this paper, a framework which allows to efficiently simulate and optimize traffic flow in a large roads' network with hundreds of vehicles is presented. The framework leverages on an Answer Set Programming (ASP) encoding to formally describe the movements of vehicles inside a network. Taking advantage of the ability to specify optimization constraints in ASP and the off-the-shelf solver Clingo, it is then possible to optimize the routes of vehicles inside the network to reduce a range of relevant metrics (e.g., travel times or emissions). Finally, an analysis on real-world traffic data is performed, utilizing the state-of-the-art Urban Mobility Simulator (SUMO) to keep track of the state of the network, test the correctness of the solution and to prove the efficiency and capabilities of the presented solution.


title: Viscous flow through a finite-width slit: Boundary conditions and dissipation

abstract: We study the hydrodynamic viscous electronic transport in a two-dimensional sample separated into two semi-infinite planes by a one-dimensional infinite barrier. The semi-infinite planes are electrically connected via the finite-size slit in the barrier. We calculate the current through the slit assuming finite voltage drop between the planes and neglecting disorder-induced Ohmic resistance, so that dissipation and resistance are purely viscosity-induced. We find that the only solution to the Stokes equation in this geometry, which yields a finite dissipation at finite resistance (and, hence, is not self-contradictory), is the one that fulfills both the no-stress and no-slip boundary conditions simultaneously. As a remarkable consequence, the obtained velocity profile satisfies the so-called "partial-slip" (Maxwell) boundary condition for any value of the slip length, which drops out from all final results. We also calculate the electronic temperature profile for the small and large heat conductivity, and find asymmetric (with respect to the barrier) temperature patterns in the former case.


title: Planning and Scheduling in Digital Health with Answer Set Programming

abstract: In the hospital world there are several complex combinatory problems, and solving these problems is important to increase the degree of patients' satisfaction and the quality of care offered. The problems in the healthcare are complex since to solve them several constraints and different type of resources should be taken into account. Moreover, the solutions must be evaluated in a small amount of time to ensure the usability in real scenarios. We plan to propose solutions to these kind of problems both expanding already tested solutions and by modelling solutions for new problems, taking into account the literature and by using real data when available. Solving these kind of problems is important but, since the European Commission established with the General Data Protection Regulation that each person has the right to ask for explanation of the decision taken by an AI, without developing Explainability methodologies the usage of AI based solvers e.g. those based on Answer Set programming will be limited. Thus, another part of the research will be devoted to study and propose new methodologies for explaining the solutions obtained.


title: Decomposition Strategies for Solving Scheduling Problems in Industrial Applications

abstract: This article presents an overview of a research study of a crucial optimization problem in the Computer Science/Operations research field: The Job-shop Scheduling Problem (JSP). The JSP is a challenging task in which a set of operations must be processed using a set of scarce machines to optimize a particular objective. The main purpose of the JSP is to determine the execution order of the processes assigned to each machine to optimize an objective. Our main interest in this study is to investigate developing decomposition strategies using logic programming to solve the JSP. We split our goal into two main phases. The first phase is to apply the decomposition approach and evaluate the proposed model by solving a set of known benchmark instances. The second phase is to apply the successful decomposition methods obtained from the first phase to solve a scheduling problem in the real-life application. In the current state, we finished the first phase and started the second one aiming to have a model that can provide a schedule of a factory for a short-time period.


title: Extended GUP formulation and the role of momentum cut-off

abstract: We analyze the extension of the GUP theory deriving from the modified uncertainty principle in agreement with the string low energy limit, which represents one of the most general formulations satisfying the Jacobi identity, in the context of the associative algebras. After providing some physical insights on the nature of the considered approaches exploiting the cosmological arena, first, we show how a natural formulation of the theory in an infinite momentum space does not lead to the emergence of a nonzero minimal uncertainty in position, then we construct a truncated formulation of the theory in momentum space, proving that only in this case we can recover the desired feature of the presence of a nonzero minimal uncertainty in position, which - as usual in these theories - can be interpreted as a phenomenological manifestation of cut-off physics effects. Both quantization schemes are completely characterized and finally applied to study wave packets' behavior and their evolution in time. The obtained results can shed light on which generalizations of the GUP theory are more coherent with the string low energy limit, in view of the existence of a minimum length in the form of a minimal uncertainty in position.


title: Automorphic Integrals with Log-polynomial Period Functions and Arithmetical Identities

abstract: Building on the works of S. Bochner on equivalence of modular relation with functional equation associated to the Dirichlet series, K. Chandrasekharan and R. Narasimhan obtained new equivalences between the functional equation and some arithmetical identities. Sister Ann M. Heath considered the functional equation in the Hawkins and Knopp context and showed its equivalence to two arithmetical identities associated with entire modular cusp integrals involving rational period functions for the full modular group. In this paper we use techniques of Chandrasekharan and Narasimhan and extend the results of Sister Ann M. Heath to entire automorphic integrals involving rational period functions on discrete Hecke group. Moreover, we establish equivalence of two arithmetical identities with a functional equation associated with automorphic integrals involving log-polynomial-period functions on the Hecke groups.


title: Growth principles of natural hypergraphs

abstract: Several systems can be represented by hypergraphs, an extension of graphs with associations between any number of vertices. These natural hypergraphs doe not appear at once. They are generated by some dynamical process of hypergraph evolution. Here I investigate what are the minimal growth principles of natural hypergraphs. I postulate edge duplication and vertex addition at edge duplications as the key principles of hypergraph growth. The implementation of these two principles induce the emergence of preferential attachment, power law degree distribution, the small-world property, high clustering coefficient and the founder effect. This work clarifies the distinction between principles, emergent properties and context specific details in the context of hypergraph growth dynamics.


title: Dynamical quantum phase transitions in Stark quantum spin chains

abstract: We investigate the nonequilibrium dynamics of one-dimension spin models in the presence of a uniform force. The linear potential induces delocalization-localization transition in the free particles model which is known as the Wannier-Stark effect. We study dynamical quantum phase transition (DQPT) due to sudden global quenches across a quantum critical point when the system undergoes a localization-delocalization transition. In this regard, we consider the XX and XXZ spin chains and explore two types of quenches with and without ramping through the delocalization-localization point. The XX model was mapped to the free fermion particles, so both analytical and numerical results were provided. Results unveil that the dynamical signature of localization-delocalization transition can be characterized by the nonanalyticities in dynamical free energy (corresponds to the zero points in the Loschmit echo). We also explore the interaction effects considering XXZ spin chains, using the time-dependent extension of the numerical DMRG technique. Our results show that depending on the anisotropic parameter $\Delta\lessgtr1.0$, if both the initial and post-quench Hamiltonian are in the same phase or not, DQPTs may happen. Moreover, the interrelation between DQPTs with different correlation measures such as the equilibrium order parameters or entanglement entropy production of the system remains unclear. We provide more analyses on the feature of DQPTs, in both types of quenches, by connecting them to the average local magnetization, entanglement entropy production, and the Schmidt gap.


title: The Value and Use of Data in Chemical Engineering Practice

abstract: The ability to generate, organize, analyze, understand and leverage data for sound decision making is a central activity of chemical engineers. Chemical engineers are responsible for the safe, profitable and environmentally friendly operation of chemical facilities; thus, they are expected to be good at designing and operating chemical processes. To this end, they make use of models which involves planning and conducting experiments in the laboratory or a pilot plant, analyzing the generated data and making use of it for designing large scale industrial systems. In an operating plant, they need to analyze data so as to achieve maximum efficiencies, reduce the use of precious natural resources, minimize environmental degradation, keep the plant safe as well as to help generate value for the customers and stakeholders. This chapter provides a non-technical view of how chemical engineers translate data to the best possible decisions that result in reliable, safe and profitable process design and operations.


title: Scattering theory with both regular and singular perturbations

abstract: We provide an asymptotic completeness criterion and a representation formula for the scattering matrix of the scattering couple $(A_B,A)$, where both $A$ and $A_B$ are self-adjoint operator and $A_B$ formally corresponds to adding to $A$ two terms, one regular and the other singular. In particular, our abstract results apply to the couple $(\Delta_B,\Delta)$, where $\Delta$ is the free self-adjoint Laplacian in $L^2(\mathbb{R}^3)$ and $\Delta_B$ is a self-adjoint operator in a class of Laplacians with both a regular perturbation, given by a short-range potential, and a singular one describing boundary conditions (like Dirichlet, Neumann and semi-transparent $\delta$ and $\delta'$ ones) at the boundary of a open, bounded Lipschitz domain. The results hinge upon a limiting absorption principle for $A_B$ and a Krein-like formula for the resolvent difference $(-A_B+z)^{-1}-(-A+z)^{-1}$ which puts on an equal footing the regular (here, in the case of the Laplacian, a Kato-Rellich potential suffices) and the singular perturbations.


title: Fixed-Point Automatic Differentiation of Forward--Backward Splitting Algorithms for Partly Smooth Functions

abstract: A large class of non-smooth practical optimization problems can be written as minimization of a sum of smooth and partly smooth functions. We consider such structured problems which also depend on a parameter vector and study the problem of differentiating its solution mapping with respect to the parameter which has far reaching applications in sensitivity analysis and parameter learning optmization problems. We show that under partial smoothness and other mild assumptions, Automatic Differentiation (AD) of the sequence generated by proximal splitting algorithms converges to the derivative of the solution mapping. For a variant of automatic differentiation, which we call Fixed-Point Automatic Differentiation (FPAD), we remedy the memory overhead problem of the Reverse Mode AD and moreover provide faster convergence theoretically. We numerically illustrate the convergence and convergence rates of AD and FPAD on Lasso and Group Lasso problems and demonstrate the working of FPAD on prototypical practical image denoising problem by learning the regularization term.


title: Complete classification of global solutions to the obstacle problem

abstract: The characterization of global solutions to the obstacle problems in $\mathbb{R}^N$, or equivalently of null quadrature domains, has been studied over more than 90 years. In this paper we give a conclusive answer to this problem by proving the following long-standing conjecture: The coincidence set of a global solution to the obstacle problem is either a half-space, an ellipsoid, a paraboloid, or a cylinder with an ellipsoid or a paraboloid as base.


title: Mean flow data assimilation based on physics-informed neural networks

abstract: Physics-informed neural networks (PINNs) can be used to solve partial differential equations (PDEs) and identify hidden variables by incorporating the governing equations into neural network training. In this study, we apply PINNs to the assimilation of turbulent mean flow data and investigate the method's ability to identify inaccessible variables and closure terms from sparse data. Using high-fidelity large-eddy simulation (LES) data and particle image velocimetry (PIV) measured mean fields, we show that PINNs are suitable for simultaneously identifying multiple missing quantities in turbulent flows and providing continuous and differentiable mean fields consistent with the provided PDEs. In this way, consistent and complete mean states can be provided, which are essential for linearized mean field methods. The presented method does not require a grid or discretization scheme, is easy to implement, and can be used for a wide range of applications, making it a very promising tool for mean field-based methods in fluid mechanics.


title: MorDeephy: Face Morphing Detection Via Fused Classification

abstract: Face morphing attack detection (MAD) is one of the most challenging tasks in the field of face recognition nowadays. In this work, we introduce a novel deep learning strategy for a single image face morphing detection, which implies the discrimination of morphed face images along with a sophisticated face recognition task in a complex classification scheme. It is directed onto learning the deep facial features, which carry information about the authenticity of these features. Our work also introduces several additional contributions: the public and easy-to-use face morphing detection benchmark and the results of our wild datasets filtering strategy. Our method, which we call MorDeephy, achieved the state of the art performance and demonstrated a prominent ability for generalising the task of morphing detection to unseen scenarios.


title: Data-free Backdoor Removal based on Channel Lipschitzness

abstract: Recent studies have shown that Deep Neural Networks (DNNs) are vulnerable to the backdoor attacks, which leads to malicious behaviors of DNNs when specific triggers are attached to the input images. It was further demonstrated that the infected DNNs possess a collection of channels, which are more sensitive to the backdoor triggers compared with normal channels. Pruning these channels was then shown to be effective in mitigating the backdoor behaviors. To locate those channels, it is natural to consider their Lipschitzness, which measures their sensitivity against worst-case perturbations on the inputs. In this work, we introduce a novel concept called Channel Lipschitz Constant (CLC), which is defined as the Lipschitz constant of the mapping from the input images to the output of each channel. Then we provide empirical evidences to show the strong correlation between an Upper bound of the CLC (UCLC) and the trigger-activated change on the channel activation. Since UCLC can be directly calculated from the weight matrices, we can detect the potential backdoor channels in a data-free manner, and do simple pruning on the infected DNN to repair the model. The proposed Channel Lipschitzness based Pruning (CLP) method is super fast, simple, data-free and robust to the choice of the pruning threshold. Extensive experiments are conducted to evaluate the efficiency and effectiveness of CLP, which achieves state-of-the-art results among the mainstream defense methods even without any data. Source codes are available at https://github.com/rkteddy/channel-Lipschitzness-based-pruning.


title: Explanation of Machine Learning Models of Colon Cancer Using SHAP Considering Interaction Effects

abstract: When using machine learning techniques in decision-making processes, the interpretability of the models is important. Shapley additive explanation (SHAP) is one of the most promising interpretation methods for machine learning models. Interaction effects occur when the effect of one variable depends on the value of another variable. Even if each variable has little effect on the outcome, its combination can have an unexpectedly large impact on the outcome. Understanding interactions is important for understanding machine learning models; however, naive SHAP analysis cannot distinguish between the main effect and interaction effects. In this paper, we introduce the Shapley-Taylor index as an interpretation method for machine learning models using SHAP considering interaction effects. We apply the method to the cancer cohort data of Kyushu University Hospital (N=29,080) to analyze what combination of factors contributes to the risk of colon cancer.


title: On the non-universality of deep learning: quantifying the cost of symmetry

abstract: We prove limitations on what neural networks trained by noisy gradient descent (GD) can efficiently learn. Our results apply whenever GD training is equivariant, which holds for many standard architectures and initializations. As applications, (i) we characterize the functions that fully-connected networks can weak-learn on the binary hypercube and unit sphere, demonstrating that depth-2 is as powerful as any other depth for this task; (ii) we extend the merged-staircase necessity result for learning with latent low-dimensional structure [ABM22] to beyond the mean-field regime. Under cryptographic assumptions, we also show hardness results for learning with fully-connected networks trained by stochastic gradient descent (SGD).


title: Effectiveness of the Krotov method in controlling open quantum systems

abstract: We apply the Krotov method for open and closed quantum systems with the objective of finding optimized controls to manipulate qubit/qutrit systems in the presence of the external environment. In the case of unitary optimization, the Krotov method is first applied to a quantum system neglecting its interaction with the environment. The resulting controls from the unitary optimization are then used to drive the system along with the environmental noise. In the case of non-unitary optimization, the Krotov method already takes into account the noise during the optimization process. We consider two distinct computational task: target-state preparation and quantum gate implementation. These tasks are carried out in simple qubit/qutrit systems and also in systems presenting leakage states. For the state-preparation cases, the controls from the non-unitary optimization outperform the controls from the unitary optimization. However, as we show here, this is not always true for the implementation of quantum gates. There are some situations where the unitary optimization performs equally well compared to the non-unitary optimization. We verify that these situations corresponds to either the absence of leakage states or to the effects of dissipation being spread uniformly over the system, including non-computational levels. For such cases, the quantum gate implementation must cover the entire Hilbert space and there is no way to dodge dissipation. On the other hand, if the subspace containing the computational levels and its complement are differently affected by dissipation, the non-unitary optimization becomes effective.


title: Multi-fidelity surrogate modeling using long short-term memory networks

abstract: When evaluating quantities of interest that depend on the solutions to differential equations, we inevitably face the trade-off between accuracy and efficiency. Especially for parametrized, time dependent problems in engineering computations, it is often the case that acceptable computational budgets limit the availability of high-fidelity, accurate simulation data. Multi-fidelity surrogate modeling has emerged as an effective strategy to overcome this difficulty. Its key idea is to leverage many low-fidelity simulation data, less accurate but much faster to compute, to improve the approximations with limited high-fidelity data. In this work, we introduce a novel data-driven framework of multi-fidelity surrogate modeling for parametrized, time-dependent problems using long short-term memory (LSTM) networks, to enhance output predictions both for unseen parameter values and forward in time simultaneously - a task known to be particularly challenging for data-driven models. We demonstrate the wide applicability of the proposed approaches in a variety of engineering problems with high- and low-fidelity data generated through fine versus coarse meshes, small versus large time steps, or finite element full-order versus deep learning reduced-order models. Numerical results show that the proposed multi-fidelity LSTM networks not only improve single-fidelity regression significantly, but also outperform the multi-fidelity models based on feed-forward neural networks.


title: Influence of errors on the transport of quantum information through distant quantum dot spin qubits

abstract: The ability to connect distant qubits plays a fundamental role in quantum computing. Therefore, quantum systems candidates for quantum computation must be able to interact all their constituent qubits. Here, we model the quantum dot spin qubits by a spin chain with nearest-neighbors interaction. Within this model, we can perform the interaction of distant qubits by the action of consecutive SWAP gates. The SWAP gate exchange the information of two different qubits and it is obtained by a time-dependent interaction of nearest-neighbors qubits that is switched on and off as the quantum information is propagated through the system. By using this scheme, we also are able to implement the CNOT gate, which is a fundamental gate to obtain universal quantum computation. These gates are probed in a system free from decoherence, which provides a very efficient connection between distant qubits. Furthermore, we analyze the situation when the dissipation is present. To perform such a task, we consider dephasing and amplitude-damping types of errors in each site of the spin chain. We found that the order of the SWAP and CNOT gates is important and it can lead to a relevant difference in fidelity when the number of qubits is large.


title: A reformulation of collision avoidance algorithm based on artificial potential fields for fixed-wing UAVs in a dynamic environment

abstract: As mini UAVs become increasingly useful in the civilian work domain, the need for a method for them to operate safely in a cluttered environment is growing, especially for fixed-wing UAVs as they are incapable of following the stop-decide-execute methodology. This paper presents preliminary research to design a reactive collision avoidance algorithm based on the improved definition of the repulsive forces used in the Artificial potential field algorithms to allow feasible and safe navigation of fixed-wing UAVs in cluttered, dynamic environments. We present simulation results of the improved definition in multiple scenarios, and we have also discussed possible future studies to improve upon these results.


title: A Design of Low-Projection SCMA Codebooks for Ultra-Low Decoding Complexity in Downlink IoT Networks

abstract: This paper conceives a novel sparse code multiple access (SCMA) codebook design which is motivated by the strong need for providing ultra-low decoding complexity and good error performance in downlink Internet-of-things (IoT) networks, in which a massive number of low-end and low-cost IoT communication devices are served. By focusing on the typical Rician fading channels, we analyze the pair-wise probability of superimposed SCMA codewords and then deduce the design metrics for multi-dimensional constellation construction and sparse codebook optimization. For significant reduction of the decoding complexity, we advocate the key idea of projecting the multi-dimensional constellation elements to a few overlapped complex numbers in each dimension, called low projection (LP). An emerging modulation scheme, called golden angle modulation (GAM), is considered for multi-stage LP optimization, where the resultant multi-dimensional constellation is called LP-GAM. Our analysis and simulation results show the superiority of the proposed LP codebooks (LPCBs) including one-shot decoding convergence and excellent error rate performance. In particular, the proposed LPCBs lead to decoding complexity reduction by at least $97\%$ compared to that of the conventional codebooks, whilst owning large minimum Euclidean distance. Some examples of the proposed LPCBs are available at \url{https://github.com/ethanlq/SCMA-codebook}.


title: Model reduction for constrained mechanical systems via spectral submanifolds

abstract: Dynamical systems are often subject to algebraic constraints in conjunction to their governing ordinary differential equations. In particular, multibody systems are commonly subject to configuration constraints that define kinematic compatibility between the motion of different bodies. A full-scale numerical simulation of such constrained problems is challenging, making reduced-order models (ROMs) of paramount importance. In this work, we show how to use spectral submanifolds (SSMs) to construct rigorous ROMs for mechanical systems with configuration constraints. These SSM-based ROMs enable the direct extraction of backbone curves and forced response curves, and facilitate efficient bifurcation analysis. We demonstrate the effectiveness of this SSM-based reduction procedure on several examples of varying complexity, including nonlinear finite-element models of multi-body systems. We also provide an open-source implementation of the proposed method that also contains all details of our numerical examples.


title: Accelerating the Sinkhorn algorithm for sparse multi-marginal optimal transport by fast Fourier transforms

abstract: We consider the numerical solution of the discrete multi-marginal optimal transport (MOT) by means of the Sinkhorn algorithm. In general, the Sinkhorn algorithm suffers from the curse of dimensionality with respect to the number of marginals. If the MOT cost function decouples according to a tree or circle, its complexity is linear in the number of marginal measures. In this case, we speed up the convolution with the radial kernel required in the Sinkhorn algorithm by non-uniform fast Fourier methods. Each step of the proposed accelerated Sinkhorn algorithm with a tree-structured cost function has a complexity of $\mathcal O(K N)$ instead of the classical $\mathcal O(K N^2)$ for straightforward matrix-vector operations, where $K$ is the number of marginals and each marginal measure is supported on at most $N$ points. In case of a circle-structured cost function, the complexity improves from $\mathcal O(K N^3)$ to $\mathcal O(K N^2)$. This is confirmed by numerical experiments.


title: Motivating explanations in Bayesian networks using MAP-independence

abstract: In decision support systems the motivation and justification of the system's diagnosis or classification is crucial for the acceptance of the system by the human user. In Bayesian networks a diagnosis or classification is typically formalized as the computation of the most probable joint value assignment to the hypothesis variables, given the observed values of the evidence variables (generally known as the MAP problem). While solving the MAP problem gives the most probable explanation of the evidence, the computation is a black box as far as the human user is concerned and it does not give additional insights that allow the user to appreciate and accept the decision. For example, a user might want to know to whether an unobserved variable could potentially (upon observation) impact the explanation, or whether it is irrelevant in this aspect. In this paper we introduce a new concept, MAP- independence, which tries to capture this notion of relevance, and explore its role towards a potential justification of an inference to the best explanation. We formalize several computational problems based on this concept and assess their computational complexity.


title: Searches for exclusive Higgs and $Z$ boson decays into a vector quarkonium state and a photon using $139$ fb$^{-1}$ of ATLAS $\sqrt{s}=13$ TeV proton$-$proton collision data

abstract: Searches for the exclusive decays of Higgs and $Z$ bosons into a vector quarkonium state and a photon are performed in the $\mu^+\mu^- \gamma$ final state with a proton$-$proton collision data sample corresponding to an integrated luminosity of $139$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN Large Hadron Collider. The observed data are compatible with the expected backgrounds. The 95% confidence-level upper limits on the branching fractions of the Higgs boson decays into $J/\psi \gamma$, $\psi(2S) \gamma$, and $\Upsilon(1S,2S,3S) \gamma$ are found to be $2.0\times10^{-4}$, $10.5\times10^{-4}$, and $(2.5,4.2,3.4)\times10^{-4}$, respectively, assuming Standard Model production of the Higgs boson. The corresponding 95% CL upper limits on the branching fractions of the $Z$ boson decays are $1.2\times10^{-6}$, $2.4\times10^{-6}$, and $(1.1,1.3,2.4)\times10^{-6}$. An observed 95% CL interval of $(-133,175)$ is obtained for the $\kappa_c/\kappa_\gamma$ ratio of Higgs boson coupling modifiers, and a 95% CL interval of $(-37,40)$ is obtained for $\kappa_b/\kappa_\gamma$.


title: Watson-Crick conjugates of words and languages

abstract: This paper explores the concept of Watson-Crick conjugates, also known as $\theta$-conjugates, of words and languages. This concept extends the classical idea of conjugates by incorporating the Watson-Crick complementarity of DNA sequences, from the perspective of DNA computing. Our investigation initially focuses on the properties of $\theta$-conjugates of words. We then define $\theta$-conjugates of a language and study closure properties of certain families of languages under the $\theta$-conjugate operation. Furthermore, we analyze the iterated $\theta$-conjugate of both words and languages. Finally, we delve into the idea of $\theta$-conjugate-free languages and examine the decidability problems surrounding $\theta$-conjugate-freeness for different classes of languages


title: Measurement of the branching fraction and search for $CP$ violation in $D^0\rightarrow K^0_S\,K^0_S\,\pi^+\pi^-$ decays at Belle

abstract: We measure the branching fraction for the Singly Cabibbo-suppressed decay $D^0\rightarrow K^0_S\,K^0_S\,\pi^+\pi^-$, and we search for $CP$ violation via a measurement of the $CP$ asymmetry $A^{}_{CP}$ and also the $T$-odd triple-product asymmetry $a_{CP}^{T}$. The later two measurements are complementary. We use 922~fb$^{-1}$ of data recorded by the Belle experiment, which ran at the KEKB asymmetric-energy $e^+ e^-$ collider. The branching fraction is measured relative to the Cabibbo-favored normalization channel $D^0\rightarrow K^0_S\,\pi^+\pi^-$. Singly Cabibbo-suppressed charm decays are expected to have especially good sensitivity to new physics effects.


title: A Strengthened SDP Relaxation for Quadratic Optimization Over the Stiefel Manifold

abstract: We study semidefinite programming (SDP) relaxations for the NP-hard problem of globally optimizing a quadratic function over the Stiefel manifold. We introduce a strengthened relaxation based on two recent ideas in the literature: (i) a tailored SDP for objectives with a block-diagonal Hessian; (ii) and the use of the Kronecker matrix product to construct SDP relaxations. Using synthetic instances on four problem classes, we show that, in general, our relaxation significantly strengthens existing relaxations, although at the expense of longer solution times.


title: ULX pulsar Swift J0243.6+6124 observations with NuSTAR -- dominance of reflected emission in the super-Eddington state

abstract: We report the discovery of the bright reflected emission component in the super-Eddington state of the ULX pulsar Swift J0243.6+6124, based on the NuSTAR observations of the source during its 2017 outburst. The flux of the reflected emission is weakly variable over the pulsar phase while the direct emission shows significantly larger pulsation amplitude. We propose that in this system the neutron star finds itself in the centre of the well formed by the inner edge of the geometrically thick super-Eddington accretion disc truncated by the magnetic field of the pulsar. The aspect ratio of the well is H/R \sim 1. The inner edge of the truncated disc is continuously illuminated by the emission of the accretion column giving rise to the weakly variable reflected emission. As the neutron star rotates, its emission sweeps through the line of sight, giving rise to the pulsating direct emission. From Doppler broadening of the iron line, we measure the truncation radius of the accretion disc \sim 50 R_g. The inferred dipole component of the magnetic field is consistent with previous estimates favouring a not very strong field. The uniqueness of this system is determined by its moderately super-Eddington accretion rate and the moderate magnetic field so that the inner edge of the truncated geometrically thick accretion disc is seen from the neutron star at a large solid angle.


title: A PTAS for Minimizing Weighted Flow Time on a Single Machine

abstract: An important objective in scheduling literature is to minimize the sum of weighted flow times. We are given a set of jobs where each job is characterized by a release time, a processing time, and a weight. Our goal is to find a preemptive schedule on a single machine that minimizes the sum of the weighted flow times of the jobs, where the flow time of a job is the time between its completion time and its release time. The currently best known polynomial time algorithm for the problem is a (2+eps)-approximation by Rohwedder and Wiese [STOC 2021] which builds on the prior break-through result by Batra, Garg, and Kumar [FOCS 2018] who found the first pseudo-polynomial time constant factor approximation algorithm for the problem, and on the result by Feige, Kulkarni, and Li [SODA 2019] who turned the latter into a polynomial time algorithm. However, it remains open whether the problem admits a PTAS. We answer this question in the affirmative and present a polynomial time (1+eps)-approximation algorithm for weighted flow time on a single machine. We rely on a reduction of the problem to a geometric covering problem, which was introduced in the mentioned (2+eps)-approximation algorithm, losing a factor 1+eps in the approximation ratio. However, unlike that algorithm, we solve the resulting instances of this problem exactly, rather than losing a factor 2+eps. Key for this is to identify and exploit structural properties of instances of the geometric covering problem which arise in the reduction from weighted flow time.


title: Time-Frequency Distributions of Heart Sound Signals: A Comparative Study using Convolutional Neural Networks

abstract: Time-Frequency Distributions (TFDs) support the heart sound characterisation and classification in early cardiac screening. However, despite the frequent use of TFDs in signal analysis, no study comprehensively compared their performances on deep learning for automatic diagnosis. Furthermore, the combination of signal processing methods as inputs for Convolutional Neural Networks (CNNs) has been proved as a practical approach to increasing signal classification performance. Therefore, this study aimed to investigate the optimal use of TFD/ combined TFDs as input for CNNs. The presented results revealed that: 1) The transformation of the heart sound signal into the TF domain achieves higher classification performance than using of raw signals. Among the TFDs, the difference in the performance was slight for all the CNN models (within $1.3\%$ in average accuracy). However, Continuous wavelet transform (CWT) and Chirplet transform (CT) outperformed the rest. 2) The appropriate increase of the CNN capacity and architecture optimisation can improve the performance, while the network architecture should not be overly complicated. Based on the ResNet or SEResNet family results, the increase in the number of parameters and the depth of the structure do not improve the performance apparently. 3) Combining TFDs as CNN inputs did not significantly improve the classification results. The findings of this study provided the knowledge for selecting TFDs as CNN input and designing CNN architecture for heart sound classification.


title: The spins of stripped B stars support magnetic internal angular momentum transport

abstract: In order to predict the spins of stellar remnants we need to understand the evolution of the internal rotation of stars, and to identify at which stage the rotation of the contracting cores of evolved stars decouples from their expanding envelopes. The donor stars of mass transferring binaries lose almost their entire envelope and may thus offer a direct view on their core rotation. After the mass transfer event they contract and fade rapidly, although they are well observable when caught in the short-lived B-star phase. The B-type primary of the galactic binary system LB-1, which was originally suggested to contain a massive black hole, is nicely explained as a stripped star accompanied by a fainter Be star. The narrow absorption lines in the primary's spectrum signify extremely slow rotation, atypical of B-type main-sequence stars. Here we investigate the evolution of mass donors in generic grids of detailed binary evolution models, where both stars include differential rotation, internal angular momentum transport, and spin-orbit coupling. Whereas the mass gainers are typically spun-up during the mass transfer, we find that the spins of the stripped donor models depend sensitively on the employed mechanism for internal angular momentum transport. Purely hydrodynamic transport cannot explain the observed slow rotation, while models including magnetic angular momentum transport are able to reproduce the observed rotation of LB-1 and similar stars, independent of the initial rotation rate. In such models the spin of the white dwarfs that emerge at the end of the evolution is independent of the mass stripping. We find evidence that the mass transfer in LB-1 was moderately non-conservative.


title: A Lightweight Machine Learning Pipeline for LiDAR-simulation

abstract: Virtual testing is a crucial task to ensure safety in autonomous driving, and sensor simulation is an important task in this domain. Most current LiDAR simulations are very simplistic and are mainly used to perform initial tests, while the majority of insights are gathered on the road. In this paper, we propose a lightweight approach for more realistic LiDAR simulation that learns a real sensor's behavior from test drive data and transforms this to the virtual domain. The central idea is to cast the simulation into an image-to-image translation problem. We train our pix2pix based architecture on two real world data sets, namely the popular KITTI data set and the Audi Autonomous Driving Dataset which provide both, RGB and LiDAR images. We apply this network on synthetic renderings and show that it generalizes sufficiently from real images to simulated images. This strategy enables to skip the sensor-specific, expensive and complex LiDAR physics simulation in our synthetic world and avoids oversimplification and a large domain-gap through the clean synthetic environment.


title: Thermal effects on nuclear matter properties

abstract: A quantitative description of the properties of hot nuclear matter will be needed for the interpretation of the available and forthcoming astrophysical data, providing information on the post merger phase of a neutron star coalescence. We have employed a recently developed theoretical model, based on a phenomenological nuclear Hamiltonian including two- and three-nucleon potentials, to study the temperature dependence of average and single-particle properties of nuclear matter relevant to astrophysical applications. The possibility to represent the results of microscopic calculations using simple and yet physically motivated parametrisations of thermal effects, suitable for use in numerical simulations of astrophysical processes, is also discussed.


title: Flow Equation and Fermion Gap in the Holographic Superconductors

abstract: We reconsider the fermion spectral function in the presence of the Cooper pair condensation and identified the interaction type of complex scalar and fermion, which gives consistent results with the expected s-wave superconductor for the first time. We derive the matrix Riccati equation, which allows the precise calculation of the fermion spectral function. Apart from the gap structure, we studied the effect of the chemical potential and the density and compared it with the BCS theory. We found that two theories give similar results in small chemical potential but very different ones in the high-density case, which we attribute to the correlation effect.


title: Out of the BLEU: how should we assess quality of the Code Generation models?

abstract: In recent years, researchers have created and introduced a significant number of various code generation models. As human evaluation of every new model version is unfeasible, the community adopted automatic evaluation metrics such as BLEU to approximate the results of human judgement. These metrics originate from the machine translation domain and it is unclear whether they are applicable for the code generation tasks and how well they agree with the human evaluation on this task. There are also other metrics, CodeBLEU and RUBY, developed to estimate the similarity of code, that take into account the properties of source code. However, for these metrics there are hardly any studies on their agreement with the human evaluation. Despite all that, minimal differences in the metric scores have been used in recent papers to claim superiority of some code generation models over the others. In this paper, we present a study on the applicability of six metrics -- BLEU, ROUGE-L, METEOR, ChrF, CodeBLEU, and RUBY -- for evaluation of code generation models. We conduct a study on two different code generation datasets and use human annotators to assess the quality of all models run on these datasets. The results indicate that for the CoNaLa dataset of Python one-liners, none of the metrics can correctly emulate human judgement on which model is better with >95% certainty if the difference in model scores is less than 5 points. For the HearthStone dataset, which consists of classes of a particular structure, a difference in model scores of at least 2 points is enough to claim the superiority of one model over the other. Our findings suggest that the ChrF metric is a better fit for the evaluation of code generation models than the commonly used BLEU and CodeBLEU. Yet, finding a metric for code generation that closely agrees with humans requires additional work.