title: Dark Energy Survey Year 3 Results: Weak Lensing Shape Catalogue

abstract: We present and characterise the galaxy shape catalogue from the first 3 years of Dark Energy Survey (DES) observations, over an effective area of ~4143 deg$^2$ of the southern sky. We describe our data analysis process and our self-calibrating shear measurement pipeline METACALIBRATION, which builds and improves upon the pipeline used in the DES Year 1 analysis in several aspects. The DES Year 3 weak-lensing shape catalogue consists of 100,204,026 galaxies, measured in the $riz$ bands, resulting in a weighted source number density of $n_{\rm eff} = 5.59$ gal/arcmin$ ^{2}$ and corresponding shape noise $\sigma_e = 0.261$. We perform a battery of internal null tests on the catalogue, including tests on systematics related to the point-spread function (PSF) modelling, spurious catalogue B-mode signals, catalogue contamination, and galaxy properties.


title: Dark Energy Survey Year 3 Results: Point-Spread Function Modeling

abstract: We introduce a new software package for modeling the point-spread function (PSF) of astronomical images, called Piff (PSFs In the Full FOV), which we apply to the first three years (known as Y3) of the Dark Energy Survey (DES) data. We describe the relevant details about the algorithms used by Piff to model the PSF, including how the PSF model varies across the field of view (FOV). Diagnostic results show that the systematic errors from the PSF modeling are very small over the range of scales that are important for the DES Y3 weak lensing analysis. In particular, the systematic errors from the PSF modeling are significantly smaller than the corresponding results from the DES year one (Y1) analysis. We also briefly describe some planned improvements to Piff that we expect to further reduce the modeling errors in future analyses.


title: Dark Energy Survey internal consistency tests of the joint cosmological probes analysis with posterior predictive distributions

abstract: Beyond-$\Lambda$CDM physics or systematic errors may cause subsets of a cosmological data set to appear inconsistent when analyzed assuming $\Lambda$CDM. We present an application of internal consistency tests to measurements from the Dark Energy Survey Year 1 (DES Y1) joint probes analysis. Our analysis relies on computing the posterior predictive distribution (PPD) for these data under the assumption of $\Lambda$CDM. We find that the DES Y1 data have an acceptable goodness of fit to $\Lambda$CDM, with a probability of finding a worse fit by random chance of ${p = 0.046}$. Using numerical PPD tests, supplemented by graphical checks, we show that most of the data vector appears completely consistent with expectations, although we observe a small tension between large- and small-scale measurements. A small part (roughly 1.5%) of the data vector shows an unusually large departure from expectations; excluding this part of the data has negligible impact on cosmological constraints, but does significantly improve the $p$-value to 0.10. The methodology developed here will be applied to test the consistency of DES Year 3 joint probes data sets.


title: Dark Energy Survey Year 3 Results: Optimizing the Lens Sample in Combined Galaxy Clustering and Galaxy-Galaxy Lensing Analysis

abstract: We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 Metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in $i$-band as a function of (photometric) redshift, $z_{\rm phot}$, and benchmark the potential gains against those using the well established redMaGiC sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies $i < 4 \, z_{\rm phot} + 18$ and has $\sim 30\%$ wider redshift distributions but $\sim 3.5$ times more galaxies than redMaGiC. Assuming a wCDM model and equivalent scale cuts to mitigate nonlinear effects, this leads to $40\%$ increase in the figure of merit for the pair combinations of $\Omega_m$, $w$, and $\sigma_8$, and gains of $16\%$ in $\sigma_8$, $10\%$ in $\Omega_m$, and $12\%$ in $w$. Similarly, in LCDM we find an improvement of $19\%$ and $27\%$ on $\sigma_8$ and $\Omega_m$, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, $i < 22.2$, has $\sim 7$ times more galaxies and $\sim 20\%$ wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.


title: The Effects of Robin Boundary Condition on Thermal Convection in a Rotating Spherical Shell

abstract: Convection in a spherical shell is widely used to model fluid layers of planets and stars. The choice of thermal boundary conditions in such models is not always straightforward. To understand the implications of this choice, we report on the effects of the thermal boundary condition on thermal convection, in terms of instability onset, fully developed transport properties and flow structure. We use the Boussinesq approximation, and impose a Robin boundary condition at the top. This enforces the temperature anomaly and its radial derivative to be linearly coupled with a proportionality factor $\beta$. Using the height H of the fluid layer, we introduce the non-dimensional Biot number Bi* = $\beta$H. Varying Bi* allows us to transition from fixed temperature for Bi* = +$\infty$, to fixed thermal flux for Bi* = 0. The bottom boundary of the shell is kept isothermal. We find that the onset of convection is only affected by Bi* in the non-rotating case. Far from onset, considering an effective Rayleigh number and a generalized Nusselt number, we show that the Nusselt and P{\'e}clet numbers follow standard universal scaling laws, independent of Bi* in all cases considered. However, the large-scale flow structure keeps the signature of the boundary condition with more vigorous large scales for smaller Bi* , even though the global heat transfer and kinetic energy are the same. Finally, for all practical purposes, the Robin condition can be safely replaced by a fixed flux when Bi* < 0.03 and by a fixed temperature for Bi* > 30.


title: Digital Nudging with Recommender Systems: Survey and Future Directions

abstract: Recommender systems are nowadays a pervasive part of our online user experience, where they either serve as information filters or provide us with suggestions for additionally relevant content. These systems thereby influence which information is easily accessible to us and thus affect our decision-making processes though the automated selection and ranking of the presented content. Automated recommendations can therefore be seen as digital nudges, because they determine different aspects of the choice architecture for users. In this work, we examine the relationship between digital nudging and recommender systems, topics that so far were mostly investigated in isolation. Through a systematic literature search, we first identified 87 nudging mechanisms, which we categorize in a novel taxonomy. A subsequent analysis then shows that only a small part of these nudging mechanisms was previously investigated in the context of recommender systems. This indicates that there is a huge potential to develop future recommender systems that leverage the power of digital nudging in order to influence the decision-making of users. In this work, we therefore outline potential ways of integrating nudging mechanisms into recommender systems.


title: Robust ENF Estimation Based on Harmonic Enhancement and Maximum Weight Clique

abstract: We present a framework for robust electric network frequency (ENF) extraction from real-world audio recordings, featuring multi-tone ENF harmonic enhancement and graph-based optimal harmonic selection. Specifically, We first extend the recently developed single-tone ENF signal enhancement method to the multi-tone scenario and propose a harmonic robust filtering algorithm (HRFA). It can respectively enhance each harmonic component without cross-component interference, thus further alleviating the effects of unwanted noise and audio content on the much weaker ENF signal. In addition, considering the fact that some harmonic components could be severely corrupted even after enhancement, disturbing rather than facilitating ENF estimation, we propose a graph-based harmonic selection algorithm (GHSA), which finds the optimal combination of harmonic components for more accurate ENF estimation. Noticeably, the harmonic selection problem is equivalently formulated as a maximum weight clique (MWC) problem in graph theory, and the Bron-Kerbosch algorithm (BKA) is adopted in the GHSA. With the enhanced and optimally selected harmonic components, both the existing maximum likelihood estimator (MLE) and weighted MLE (WMLE) are incorporated to yield the final ENF estimation results. The proposed framework is extensively evaluated using both synthetic signals and our ENF-WHU dataset consisting of $130$ real-world audio recordings, demonstrating substantially improved capability of extracting the ENF from realistically noisy observations over the existing single- and multi-tone competitors. This work further improves the applicability of the ENF as a forensic criterion in real-world situations.


title: Searching beyond the usual papillomavirus suspects in squamous carcinomas of the vulva, penis and head and neck

abstract: Human Papillomaviruses (HPVs) are involved in the etiology of anogenital and head and neck cancers. The HPV DNA prevalence greatly differs by anatomical site. Indeed, the high rates of viral DNA prevalence in anal and cer-vical carcinomas contrast with the lower fraction of cancer cases attributable to HPVs in other anatomical sites, chiefly the vulva, the penis and head and neck. Here we analyzed 2635 Formalin Fixed Paraffin Embedded surgical samples that had previously tested negative for the presence of HPVs DNA using the SPF10/DEIA procedure, in order to identify the presence of other PVs not explicitly targeted by standard molecular epidemiologic approaches. All samples were reanalyzed using five broad-PV PCR primer sets (CP1/2, FAP6064/FAP64, SKF/SKR, MY9/MY11, MFI/MFII) targeting the main PV main clades. In head and neck carcinoma samples (n = 1141), we recovered DNA from two BetaHPVs, namely HPV20 and HPV21, and from three cutaneous AlphaPVs, namely HPV2, HPV57 and HPV61. In vulvar squamous cell carcinoma samples (n = 902), we found one of the samples containing DNA of one cutaneous HPV, namely HPV2, and 29 samples contained DNA from essentially mucosal HPVs. In penile squamous cell carcinoma samples (n = 592), we retrieved the DNA of HPV16 in 16 samples. Our results show first that the SPF10/DEIA is very sensitive, as we recovered only 2.1% (55/2635) false negative results; second, that although the DNA of cutaneous HPVs can be detected in cancer samples, their relative contribution remains anyway minor (0.23%; 6/2635) and may be neglected for screening and vaccination purposes; and third, their contribution to malignancy is not necessarily warranted and needs to be elucidated.


title: Eulerian ideals

abstract: Let $G$ be a simple graph and $I(X_G)=\varphi^{-1}(x_i^2-x_j^2 : i,j\in V_G)$, where $\varphi \colon K[E_G]\to K[V_G]$ is the homomorphism that sends an edge to the product of its vertices. The ideal $I(X_G)$ is Cohen--Macaulay, one-dimensional and binomial. If $G$ is bipartite, it is known that the Castelnuovo--Mumford regularity of $I(X_G)$ is equal to the maximum cardinality of a set of edges having no more than half of the edges of any Eulerian subgraph of $G$. Here, with respect to the grevlex order associated to an ordering of the edge set of $G$, we describe a Gr\"obner basis for $I(X_G)$, and we characterize the standard monomials of the ideal $(I(X_G),t_e)$ in terms of even sets of vertices marked with a parity. Using these results, we give a combinatorial interpretation of the degree of $I(X_G)$, via the set of even sets of vertices of $G$; and we show that the Castelnuovo--Mumford regularity of $I(X_G)$, for any graph, is the maximum cardinality of a set of edges having no more than half of the edges of any \emph{even} Eulerian subgraph of $G$ or, equivalently, the maximum cardinality of a minimum fixed parity $T$-join.


title: Physical Layer Security of Intelligent Reflective Surface Aided NOMA Networks

abstract: Intelligent reflective surface (IRS) technology is emerging as a promising performance enhancement technique for next-generation wireless networks. Hence, we investigate the physical layer security of the downlink in IRS-aided non-orthogonal multiple access networks in the presence of an eavesdropper, where an IRS is deployed for enhancing the quality by assisting the cell-edge user to communicate with the base station. To characterize the network's performance, the expected value of the new channel statistics is derived for the reflected links in the case of Nakagami-m fading. Furthermore, the performance of the proposed network is evaluated both in terms of the secrecy outage probability (SOP) and the average secrecy capacity (ASC). The closed-form expressions of the SOP and the ASC are derived. We also study the impact of various network parameters on the overall performance of the network considered. To obtain further insights, the secrecy diversity orders and the high signal-to-noise ratio slopes are obtained. We finally show that: 1) the expectation of the channel gain in the reflected links is determined both by the number of IRSs and by the Nakagami-m fading parameters; 2) The SOP of both receiver 1 and receiver 2 becomes unity, when the number of IRSs is sufficiently high; 3) The secrecy diversity orders are affected both by the number of IRSs and by the Nakagami-m fading parameters, whereas the high-SNR slopes are not affected by these parameters. Our Monte-Carlo simulations perfectly demonstrate the analytical results.


title: Local Two-Sample Testing over Graphs and Point-Clouds by Random-Walk Distributions

abstract: Rejecting the null hypothesis in two-sample testing is a fundamental tool for scientific discovery. Yet, aside from concluding that two samples do not come from the same probability distribution, it is often of interest to characterize how the two distributions differ. Given samples from two densities $f_1$ and $f_0$, we consider the task of localizing occurrences of the inequality $f_1 > f_0$. To avoid the challenges associated with high-dimensional space, we propose a general hypothesis testing framework where hypotheses are formulated adaptively to the data by conditioning on the combined sample from the two densities. We then investigate a special case of this framework where the notion of locality is captured by a random walk on a weighted graph constructed over this combined sample. We derive a tractable testing procedure for this case employing a type of scan statistic, and provide non-asymptotic lower bounds on the power and accuracy of our test to detect whether $f_1>f_0$ in a local sense. Furthermore, we characterize the test's consistency according to a certain problem-hardness parameter, and show that our test achieves the minimax detection rate for this parameter. We conduct numerical experiments to validate our method, and demonstrate our approach on two real-world applications: detecting and localizing arsenic well contamination across the United States, and analyzing two-sample single-cell RNA sequencing data from melanoma patients.


title: Design Of Drug-Like Protein-Protein Interaction Stabilizers Guided By Chelation-Controlled Bioactive Conformation Stabilization

abstract: The protein-protein interactions (PPIs) of 14-3-3 proteins are a model system for studying PPI stabilization. The complex natural product Fusicoccin A stabilizes many 14-3-3 PPIs but is not amenable for use in SAR studies, motivating the search for more drug-like chemical matter. However, drug-like 14-3-3 PPI stabilizers enabling such study have remained elusive. An X-ray crystal structure of a PPI in complex with an extremely low potency stabilizer uncovered an unexpected non-protein interacting, ligand-chelated Mg 2+ leading to the discovery of metal ion-dependent 14-3-3 PPI stabilization potency. This originates from a novel chelation-controlled bioactive conformation stabilization effect. Metal chelation has been associated with pan-assay interference compounds (PAINS) and frequent hitter behavior, but chelation can evidently also lead to true potency gains and find use as a medicinal chemistry strategy to guide compound optimization. To demonstrate this, we exploited the effect to design the first potent, selective and drug-like 14-3-3 PPI stabilizers.


title: On characterizations of amenable C*-dynamical systems and new examples

abstract: We show the equivalence of several amenability type conditions for C*-dynamical systems. As an application, based on the Pimsner--Meyer construction, we give a new method to produce amenable actions on simple C*-algebras.


title: Polariton response in the presence of Brownian dissipation from molecular vibrations

abstract: We study the elastic response of a stationarily driven system of a cavity field strongly coupled with molecular excitons, taking into account the main dissipation channels due to the finite cavity linewidth and molecular vibrations. We show that the frequently used coupled oscillator model fails in describing this response especially due to the non-Lorentzian dissipation of the molecules to their vibrations. Signatures of this failure are the temperature dependent minimum point of the polariton peak splitting, uneven polariton peak height at the minimum splitting, and the asymmetric shape of the polariton peaks even at the experimentally accessed "zero-detuning" point. Using a rather generic yet representative model of molecular vibrations, we predict the polariton response in various conditions, depending on the temperature, molecular Stokes shift and vibration frequencies, and the size of the Rabi splitting. Our results can be used as a sanity check of the experiments trying to "prove" results originating from strong coupling, such as vacuum-enhanced chemical reaction rate.


title: Development and validation of HERWIG 7 tunes from CMS underlying-event measurements

abstract: This paper presents a new set of parameters ("tunes") for the underlying-event model of the HERWIG 7 event generator. These parameters control the description of multiple-parton interactions (MPI) and colour reconnection in HERWIG 7, and are obtained from a fit to minimum-bias data collected by the CMS experiment at $\sqrt{s} =$ 0.9, 7, and 13 TeV. The tunes are based on the NNPDF 3.1 next-to-next-to-leading-order parton distribution function (PDF) set for the parton shower, and either a leading-order or next-to-next-to-leading-order PDF set for the simulation of MPI and the beam remnants. Predictions utilizing the tunes are produced for event shape observables in electron-positron collisions, and for minimum-bias, inclusive jet, top quark pair, and Z and W boson events in proton-proton collisions, and are compared with data. Each of the new tunes describes the data at a reasonable level, and the tunes using a leading-order PDF for the simulation of MPI provide the best description of the data.


title: Analyzing the Productivity of GitHub Teams based on Formation Phase Activity

abstract: Our goal is to understand the characteristics of high-performing teams on GitHub. Towards this end, we collect data from software repositories and evaluate teams by examining differences in productivity. Our study focuses on the team formation phase, the first six months after repository creation. To better understand team activity, we clustered repositories based on the proportion of their work activities and discovered three work styles in teams: toilers, communicators, and collaborators. Based on our results, we contend that early activities in software development repositories on GitHub establish coordination processes that enable effective collaborations over time.


title: Session-aware Recommendation: A Surprising Quest for the State-of-the-art

abstract: Recommender systems are designed to help users in situations of information overload. In recent years, we observed increased interest in session-based recommendation scenarios, where the problem is to make item suggestions to users based only on interactions observed in an ongoing session. However, in cases where interactions from previous user sessions are available, the recommendations can be personalized according to the users' long-term preferences, a process called session-aware recommendation. Today, research in this area is scattered and many existing works only compare session-aware with session-based models. This makes it challenging to understand what represents the state-of-the-art. To close this research gap, we benchmarked recent session-aware algorithms against each other and against a number of session-based recommendation algorithms and trivial extensions thereof. Our comparison, to some surprise, revealed that (i) item simple techniques based on nearest neighbors consistently outperform recent neural techniques and that (ii) session-aware models were mostly not better than approaches that do not use long-term preference information. Our work therefore not only points to potential methodological issues where new methods are compared to weak baselines, but also indicates that there remains a huge potential for more sophisticated session-aware recommendation algorithms.


title: C-ITS bundling for integrated traffic management

abstract: Cooperative Intelligent Transportation Systems (C-ITS) enable vehicles communication with each other (Vehicle-to-Vehicle, V2V) and with roadside infrastructure (Vehicle-to-Infrastructure, V2I). In the context of traffic efficiency, C-ITS technologies could assist in road network status visualization and monitoring, through data exchange, improving this way traffic control organization and traffic management implementation. Bundling is the provision of several C-ITS services as one combined service. The purpose of bundling is to harvest the usability of C-ITS services by developing a strategy for the operation and exploitation of services in real-time and within varying geographical areas. Two different dimensions of bundling have been recognized covering: 1) end-users, and 2) operators-managers. The objective of the operators-managers dimension is the integration of C-ITS services in operational traffic management. This work spotlights the operators-managers bundling dimension, presenting a framework based on a step-by-step approach for integrating C-ITS services in traditional traffic management.


title: Self-Supervised Learning from Contrastive Mixtures for Personalized Speech Enhancement

abstract: This work explores how self-supervised learning can be universally used to discover speaker-specific features towards enabling personalized speech enhancement models. We specifically address the few-shot learning scenario where access to cleaning recordings of a test-time speaker is limited to a few seconds, but noisy recordings of the speaker are abundant. We develop a simple contrastive learning procedure which treats the abundant noisy data as makeshift training targets through pairwise noise injection: the model is pretrained to maximize agreement between pairs of differently deformed identical utterances and to minimize agreement between pairs of similarly deformed nonidentical utterances. Our experiments compare the proposed pretraining approach with two baseline alternatives: speaker-agnostic fully-supervised pretraining, and speaker-specific self-supervised pretraining without contrastive loss terms. Of all three approaches, the proposed method using contrastive mixtures is found to be most robust to model compression (using 85% fewer parameters) and reduced clean speech (requiring only 3 seconds).


title: Hyperoctahedral Homology for Involutive Algebras

abstract: Hyperoctahedral homology is the homology theory associated to the hyperoctahedral crossed simplicial group. It is defined for involutive algebras over a commutative ring using functor homology and the hyperoctahedral bar construction of Fiedorowicz. The main result of the paper proves that hyperoctahedral homology is related to equivariant stable homotopy theory: for a discrete group of odd order, the hyperoctahedral homology of the group algebra is isomorphic to the homology of the fixed points under the involution of an equivariant infinite loop space built from the classifying space of the group.


title: Illumination Normalization by Partially Impossible Encoder-Decoder Cost Function

abstract: Images recorded during the lifetime of computer vision based systems undergo a wide range of illumination and environmental conditions affecting the reliability of previously trained machine learning models. Image normalization is hence a valuable preprocessing component to enhance the models' robustness. To this end, we introduce a new strategy for the cost function formulation of encoder-decoder networks to average out all the unimportant information in the input images (e.g. environmental features and illumination changes) to focus on the reconstruction of the salient features (e.g. class instances). Our method exploits the availability of identical sceneries under different illumination and environmental conditions for which we formulate a partially impossible reconstruction target: the input image will not convey enough information to reconstruct the target in its entirety. Its applicability is assessed on three publicly available datasets. We combine the triplet loss as a regularizer in the latent space representation and a nearest neighbour search to improve the generalization to unseen illuminations and class instances. The importance of the aforementioned post-processing is highlighted on an automotive application. To this end, we release a synthetic dataset of sceneries from three different passenger compartments where each scenery is rendered under ten different illumination and environmental conditions: see https://sviro.kl.dfki.de


title: Nonlinear Quantum Neuron: A Fundamental Building Block for Quantum Neural Networks

abstract: Quantum computing enables quantum neural networks (QNNs) to have great potentials to surpass artificial neural networks (ANNs). The powerful generalization of neural networks is attributed to nonlinear activation functions. Although various models related to QNNs have been developed, they are facing the challenge of merging the nonlinear, dissipative dynamics of neural computing into the linear, unitary quantum system. In this paper, we establish different quantum circuits to approximate nonlinear functions and then propose a generalizable framework to realize any nonlinear quantum neuron. We present two quantum neuron examples based on the proposed framework. The quantum resources required to construct a single quantum neuron are the polynomial, in function of the input size. Finally, both IBM Quantum Experience results and numerical simulations illustrate the effectiveness of the proposed framework.


title: The solar eclipse of August 30, 1905 at Alcala de Chisvert

abstract: Solar eclipses play an important role in the history of astronomy because their observation allowed astrophysics to expand. The eclipse of August 30th, 1905 is particularly interesting since the line of totality passes through Spain and Algeria, and that it will last 3 minutes and 45 seconds. Consequently, all the French observatories send astronomers to observe it. About a hundred amateur astronomers of the French Astronomical Society also go to various sites, as well as quite a few individuals who are simply curious, like Camille Saint-Sa\"ens. Count Aymar de la Baume Pluvinel and his two young assistants settle in at Alcala de Chisvert, between Tarragona and Valencia.


title: Continuous surrogate-based optimization algorithms are well-suited for expensive discrete problems

abstract: One method to solve expensive black-box optimization problems is to use a surrogate model that approximates the objective based on previous observed evaluations. The surrogate, which is cheaper to evaluate, is optimized instead to find an approximate solution to the original problem. In the case of discrete problems, recent research has revolved around surrogate models that are specifically constructed to deal with discrete structures. A main motivation is that literature considers continuous methods, such as Bayesian optimization with Gaussian processes as the surrogate, to be sub-optimal (especially in higher dimensions) because they ignore the discrete structure by, e.g., rounding off real-valued solutions to integers. However, we claim that this is not true. In fact, we present empirical evidence showing that the use of continuous surrogate models displays competitive performance on a set of high-dimensional discrete benchmark problems, including a real-life application, against state-of-the-art discrete surrogate-based methods. Our experiments on different discrete structures and time constraints also give more insight into which algorithms work well on which type of problem.


title: Misalignment Recognition in Acoustic Sensor Networks using a Semi-supervised Source Estimation Method and Markov Random Fields

abstract: In this paper, we consider the problem of acoustic source localization by acoustic sensor networks (ASNs) using a promising, learning-based technique that adapts to the acoustic environment. In particular, we look at the scenario when a node in the ASN is displaced from its position during training. As the mismatch between the ASN used for learning the localization model and the one after a node displacement leads to erroneous position estimates, a displacement has to be detected and the displaced nodes need to be identified. We propose a method that considers the disparity in position estimates made by leave-one-node-out (LONO) sub-networks and uses a Markov random field (MRF) framework to infer the probability of each LONO position estimate being aligned, misaligned or unreliable while accounting for the noise inherent to the estimator. This probabilistic approach is advantageous over naive detection methods, as it outputs a normalized value that encapsulates conditional information provided by each LONO sub-network on whether the reading is in misalignment with the overall network. Experimental results confirm that the performance of the proposed method is consistent in identifying compromised nodes in various acoustic conditions.


title: Detecting and Counting Small Subgraphs, and Evaluating a Parameterized Tutte Polynomial: Lower Bounds via Toroidal Grids and Cayley Graph Expanders

abstract: Given a graph property $\Phi$, we consider the problem $\mathtt{EdgeSub}(\Phi)$, where the input is a pair of a graph $G$ and a positive integer $k$, and the task is to decide whether $G$ contains a $k$-edge subgraph that satisfies $\Phi$. Specifically, we study the parameterized complexity of $\mathtt{EdgeSub}(\Phi)$ and of its counting problem $\#\mathtt{EdgeSub}(\Phi)$ with respect to both approximate and exact counting. We obtain a complete picture for minor-closed properties $\Phi$: the decision problem $\mathtt{EdgeSub}(\Phi)$ always admits an FPT algorithm and the counting problem $\#\mathtt{EdgeSub}(\Phi)$ always admits an FPTRAS. For exact counting, we present an exhaustive and explicit criterion on the property $\Phi$ which, if satisfied, yields fixed-parameter tractability and otherwise $\#\mathsf{W[1]}$-hardness. Additionally, most of our hardness results come with an almost tight conditional lower bound under the so-called Exponential Time Hypothesis, ruling out algorithms for $\#\mathtt{EdgeSub}(\Phi)$ that run in time $f(k)\cdot|G|^{o(k/\log k)}$ for any computable function $f$. As a main technical result, we gain a complete understanding of the coefficients of toroidal grids and selected Cayley graph expanders in the homomorphism basis of $\#\mathtt{EdgeSub}(\Phi)$. This allows us to establish hardness of exact counting using the Complexity Monotonicity framework due to Curticapean, Dell and Marx (STOC'17). Our methods can also be applied to a parameterized variant of the Tutte Polynomial $T^k_G$ of a graph $G$, to which many known combinatorial interpretations of values of the (classical) Tutte Polynomial can be extended. As an example, $T^k_G(2,1)$ corresponds to the number of $k$-forests in the graph $G$. Our techniques allow us to completely understand the parametrized complexity of computing the evaluation of $T^k_G$ at every pair of rational coordinates $(x,y)$.


title: Maximum Matchings and Popularity

abstract: Let $G$ be a bipartite graph where every node has a strict ranking of its neighbors. For every node, its preferences over neighbors extend naturally to preferences over matchings. Matching $N$ is more popular than matching $M$ if the number of nodes that prefer $N$ to $M$ is more than the number that prefer $M$ to $N$. A maximum matching $M$ in $G$ is a "popular max-matching" if there is no maximum matching in $G$ that is more popular than $M$. Such matchings are relevant in applications where the set of admissible solutions is the set of maximum matchings and we wish to find a best maximum matching as per node preferences. It is known that a popular max-matching always exists in $G$. Here we show a compact extended formulation for the popular max-matching polytope. So when there are edge costs, a min-cost popular max-matching in $G$ can be computed in polynomial time. This is in contrast to the min-cost popular matching problem which is known to be NP-hard. We also consider Pareto-optimality, which is a relaxation of popularity, and show that computing a min-cost Pareto-optimal matching/max-matching is NP-hard.


title: Answer Span Correction in Machine Reading Comprehension

abstract: Answer validation in machine reading comprehension (MRC) consists of verifying an extracted answer against an input context and question pair. Previous work has looked at re-assessing the "answerability" of the question given the extracted answer. Here we address a different problem: the tendency of existing MRC systems to produce partially correct answers when presented with answerable questions. We explore the nature of such errors and propose a post-processing correction method that yields statistically significant performance improvements over state-of-the-art MRC systems in both monolingual and multilingual evaluation.


title: Homothetic packings of centrally symmetric convex bodies

abstract: A centrally symmetric convex body is a convex compact set with non-empty interior that is symmetric about the origin. Of particular interest are those that are both smooth and strictly convex -- known here as regular symmetric bodies -- since they retain many of the useful properties of the $d$-dimensional Euclidean ball. We prove that for any given regular symmetric body $C$, a homothetic packing of copies of $C$ with randomly chosen radii will have a $(2,2)$-sparse planar contact graph. We further prove that there exists a comeagre set of centrally symmetric convex bodies $C$ where any $(2,2)$-sparse planar graph can be realised as the contact graph of a stress-free homothetic packing of $C$.


title: A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE). VIII. Modeling ram pressure stripping of diffuse gas in the Virgo cluster spiral galaxy NGC 4330

abstract: NGC 4330 is one of the Virgo galaxies whose UV distribution shows a tail structure. An associated tail structure is also observed in the HI and H$\alpha$ emission distributions. Previous dynamical modeling showed that the galaxy is approaching the cluster center and is therefore undergoing increasing ram pressure stripping. Recent stellar population fitting of deep optical spectra together with multiband photometry lead to the determination of the time when star formation was quenched in the galactic disk. We introduce a new version of the dynamical model that includes the diffuse ionized gas and aim to reproduce the HI, H$\alpha$, UV distributions together with the star formation histories of the outer gas-free parts of the galactic disk. The results of 50 simulations with five different Lorentzian temporal ram-pressure profiles and five different delays between the simulation onset and peak ram pressure are presented. The inclusion of diffuse gas stripping changes significantly the HI, UV, and H$\alpha$ emission distributions. The simulations with diffuse gas stripping naturally lead to vertical low surface density filaments in the downwind region of the galactic disk. These filaments occur less frequently in the simulations without diffuse gas stripping. The simulations with diffuse gas stripping lead to better joint fits to the SEDs and optical spectra. The HI, NUV, and H$\alpha$ morphologies of the model snapshots which best reproduce the SEDs and optical spectra are sufficiently different to permit a selection of best-fit models. We conclude that the inclusion of diffuse gas stripping significantly improves the resemblance between the model and observations. Our preferred model yields a time to peak ram pressure of 140 Myr in the future. The spatial coincidence of the radio continuum and diffuse H$\alpha$ tails suggests that both gas phases are stripped together.


title: Stochastic optimal control formalism for an open quantum system

abstract: A stochastic procedure is developed which allows one to express Pontryagin's maximum principle for dissipative quantum system solely in terms of stochastic wave functions. Time-optimal controls can be efficiently computed without computing the density matrix. Specifically, the proper dynamical update rules are presented for the stochastic costate variables introduced by Pontryagin's maximum principle and restrictions on the form of the terminal cost function are discussed. The proposed procedure is confirmed by comparing the results to those obtained from optimal control on Lindbladian dynamics. Numerically, the proposed formalism becomes time and memory efficient for large systems, and it can be generalized to describe non-Markovian dynamics.


title: Ackermann and Goodstein go functorial

abstract: We present variants of Goodstein's theorem that are equivalent to arithmetical comprehension and to arithmetical transfinite recursion, respectively, over a weak base theory. These variants differ from the usual Goodstein theorem in that they (necessarily) entail the existence of complex infinite objects. As part of our proof, we show that the Veblen hierarchy of normal functions on the ordinals is closely related to an extension of the Ackermann function by direct limits.


title: Wave topology brought to the coast

abstract: Since the pioneering work of Kelvin on Laplace tidal equations, a zoology of trapped waves have been found in the context of coastal dynamics. Among them, the one originally computed by Kelvin plays a particular role, as it is an unidirectional mode filling a frequency gap between different wave bands. The existence of such Kelvin waves is robust to changes in the boundary shape and in changes of the underlying model for the coast. This suggests a topological interpretation that has yet up to now remained elusive. Here we rectify the situation, by taking advantage of a reformulation of the shallow water dynamics that highlights an analogy with the celebrated Haldane model in condensed matter physics. For any profile of bottom topography, the number of modes that transit from one wave band to another in the dispersion relation is predicted by computing a first Chern number describing the topology of complex eigenmodes in a dual, simpler wave problem.


title: Nighttime monitoring of the aerosol content of the lower atmosphere by differential photometry of the anthropogenic skyglow

abstract: Nighttime monitoring of the aerosol content of the lower atmosphere is a challenging task, because appropriate reference natural light sources are lacking. Here we show that the anthropogenic night sky brightness due to city lights can be successfully used for estimating the aerosol optical depth of arbitrarily thick atmospheric layers. This method requires measuring the zenith night sky brightness with two detectors located at the limiting layer altitudes. Combined with an estimate of the overall atmospheric optical depth (available from ground-based measurements or specific satellite products), the ratio of these radiances provides a direct estimate of the differential aerosol optical depth of the air column between these two altitudes. These measurements can be made with single-channel low-cost radiance detectors widely used by the light pollution research community.


title: DeepFoldit -- A Deep Reinforcement Learning Neural Network Folding Proteins

abstract: Despite considerable progress, ab initio protein structure prediction remains suboptimal. A crowdsourcing approach is the online puzzle video game Foldit, that provided several useful results that matched or even outperformed algorithmically computed solutions. Using Foldit, the WeFold crowd had several successful participations in the Critical Assessment of Techniques for Protein Structure Prediction. Based on the recent Foldit standalone version, we trained a deep reinforcement neural network called DeepFoldit to improve the score assigned to an unfolded protein, using the Q-learning method with experience replay. This paper is focused on model improvement through hyperparameter tuning. We examined various implementations by examining different model architectures and changing hyperparameter values to improve the accuracy of the model. The new model hyper-parameters also improved its ability to generalize. Initial results, from the latest implementation, show that given a set of small unfolded training proteins, DeepFoldit learns action sequences that improve the score both on the training set and on novel test proteins. Our approach combines the intuitive user interface of Foldit with the efficiency of deep reinforcement learning.


title: Is Transfer Learning Necessary for Protein Landscape Prediction?

abstract: Recently, there has been great interest in learning how to best represent proteins, specifically with fixed-length embeddings. Deep learning has become a popular tool for protein representation learning as a model's hidden layers produce potentially useful vector embeddings. TAPE introduced a number of benchmark tasks and showed that semi-supervised learning, via pretraining language models on a large protein corpus, improved performance on downstream tasks. Two of the tasks (fluorescence prediction and stability prediction) involve learning fitness landscapes. In this paper, we show that CNN models trained solely using supervised learning both compete with and sometimes outperform the best models from TAPE that leverage expensive pretraining on large protein datasets. These CNN models are sufficiently simple and small that they can be trained using a Google Colab notebook. We also find for the fluorescence task that linear regression outperforms our models and the TAPE models. The benchmarking tasks proposed by TAPE are excellent measures of a model's ability to predict protein function and should be used going forward. However, we believe it is important to add baselines from simple models to put the performance of the semi-supervised models that have been reported so far into perspective.


title: Optomechanical crystals for spatial sensing of submicron sized analytes

abstract: Optomechanical crystal cavities have rich perspectives for detecting and indirectly analysing biological particles, such as proteins, bacteria and viruses. In this work we demonstrate the working principle of an optomechanical crystal cavity operating under ambient conditions as a sensor of submicrometer analytes by optically monitoring the frequency shift of thermally activated mechanical modes. The resonator has been specifically designed so that the cavity region supports a particular family of low modal-volume mechanical modes, commonly known as -pinch modes-. These involve the oscillation of only a couple of adjacent cavity cells that are relatively insensitive to perturbations in other parts of the resonator. The eigenfrequency of these modes decreases as the deformation is localized closer to the centre of the resonator. Thus, by identifying specific modes that undergo a frequency shift that amply exceeds the mechanical linewidth, it is possible to infer if there are particles deposited on the resonator, how many are there and their approximate position within the cavity region.


title: Olbertian partition function in scalar field theory

abstract: The Olbertian partition function is reformulated in terms of continuous (Abelian) fields described by the Landau-Ginzburg action, respectively Hamiltonian. In order do make some progress, the Gaussian approximation to the partition function is transformed into the Olbertian prior to adding the quartic Landau-Ginzburg term in the Hamiltonian. The final result is provided in the form of an expansion suitable for application of diagrammatic techniques once the nature of the field is given, i.e. once the field equations are written down such that the interactions can be formulated.


title: Light-scanning hand-held photoacoustic probe design

abstract: Significance: We proposed a new design of hand-held linear-array photoacoustic (PA) probe which can acquire multi images via motor moving. Moreover, images from different locations are utilized via imaging fusion for SNR enhancement. Aim: We devised an adjustable hand-held for the purpose of realizing different images at diverse location for further image fusion. For realizing the light spot which is more matched with the Ultrasonic transducer detection area, we specially design a light adjust unit. Moreover, due to no displacement among the images, there is no need to execute image register process. The program execution time be reduced, greatly. Approach: mechanical design; Montel carol simulation; no-registration image fusion; Spot compression. Results: Multiple PA images with different optical illumination areas were acquired. After image fusion, we obtained fused PA images with higher signal-to-noise-ratio (SNR) and image fidelity than each single PA image. A quantitative comparison shows that the SNR of fused image is improved by 36.06% in agar-milk phantom, and 44.69% in chicken breast phantom, respectively. Conclusions: In this paper, the light scanning adjustable hand-held PA imaging probe is proposed, which can realize PA imaging with different illumination positions via adjusting the optical unit.


title: Convergence results for an averaged LQR problem with applications to reinforcement learning

abstract: In this paper, we will deal with a Linear Quadratic Optimal Control problem with unknown dynamics. As a modeling assumption, we will suppose that the knowledge that an agent has on the current system is represented by a probability distribution $\pi$ on the space of matrices. Furthermore, we will assume that such a probability measure is opportunely updated to take into account the increased experience that the agent obtains while exploring the environment, approximating with increasing accuracy the underlying dynamics. Under these assumptions, we will show that the optimal control obtained by solving the "average" Linear Quadratic Optimal Control problem with respect to a certain $\pi$ converges to the optimal control driven related to the Linear Quadratic Optimal Control problem governed by the actual, underlying dynamics. This approach is closely related to model-based Reinforcement Learning algorithms where prior and posterior probability distributions describing the knowledge on the uncertain system are recursively updated. In the last section, we will show a numerical test that confirms the theoretical results.


title: Partial decay widths of $\phi(2170)$ to kaonic resonances

abstract: In this work we study the strong decays of $\phi(2170)$ to final states involving the kaonic resonances $K(1460)$, $K_1(1270)$ and $K_1(1400)$, on which experimental data have recently been extracted by the BESIII Collaboration. The formalism developed here is based on interpreting $\phi(2170)$ and $K(1460)$ as states arising from three-hadron dynamics, which is inspired by our earlier works. For $K_1(1270)$ and $K_1(1400)$ we investigate different descriptions, such as a mixture of states belonging to the nonet of axial resonances, or the former one as a state originating from the vector-pseudoscalar dynamics. The ratios among the partial widths of $K^+(1460)K^-$, $K^+_1(1400)K^-$ and $K^+_1(1270)K^-$ obtained are compatible with the experimental results, reinforcing the three-body nature of $\phi(2170)$. Within our formalism, we can also explain the suppressed decay of $\phi(2170)$ to $K^*(892) \bar K^*(892)$, as found by the BESIII Collaboration. Furthermore, our results can be useful in clarifying the properties of $K(1460)$, $K_1(1270)$ and $K_1(1400)$ when higher statistics data would be available.


title: DRAST -- A Deep Learning and AST Based Approach for Bug Localization

abstract: Context: Given a bug report and source code of the project, bug localization can help developers to focus on fixing probable buggy files rather than searching the entire source code repository. While existing research uses information retrieval (IR) and/or combination of machine learning (ML) or deep learning (DL) approaches, they focus primarily on benchmark Java projects, and also motivate the need for multi-language bug localization approach. Objective: To create a novel bug localization approach that leverages the syntactic structure of source code, bug report information and which can support multi-language projects along with a new dataset of C projects. Method: The proposed DRAST approach represents source code as code vectors by using its high-level AST and combines rVSM, an IR technique with ML/DL models such as Random Forest and Deep Neural Network regressor to rank the list of buggy files. We also use features such as textual similarity using IR techniques, lexical mismatch using DNNs, and history of the project using the metadata of BugC dataset. Results: We tested DRAST on seven projects from the BugC dataset, which consists of 2462 bug reports from 21 open-source C projects. The results show that DRAST can locate correct buggy files 90% of the time from top 1, 5, and 10 suggested files with MAP and MRR scores of above 90% for the randomly selected seven projects. We also tested DRAST on Tomcat and AspectJ, projects from benchmark dataset with better results at accuracy@1, MAP and MRR when compared with state-of-the-art. Conclusions: This paper presents a novel bug localization approach that works on C and Java projects and a bug localization C dataset along with a novel source code representation. The results for C projects using DRAST are promising and could motivate researchers/practitioners to focus on developing and creating multi-language bug localization approaches.


title: Compactness of Sobolev embeddings and decay of norms

abstract: We investigate the relationship between the compactness of embeddings of Sobolev spaces built upon rearrangement-invariant spaces into rearrangement-invariant spaces endowed with $d$-Ahlfors measures under certain restriction on the speed of its decay on balls. We show that the gateway to compactness of such embeddings, while formally describable by means of optimal embeddings and almost-compact embeddings, is quite elusive. It is known that such a Sobolev embedding is not compact when its target space has the optimal fundamental function. We show that, quite surprisingly, such a target space can actually be "fundamentally enlarged", and yet the resulting embedding remains noncompact. In order to do that, we develop two different approaches. One is based on enlarging the optimal target space itself, and the other is based on enlarging the Marcinkiewicz space corresponding to the optimal fundamental function.


title: Deep Cross-modal Hashing via Margin-dynamic-softmax Loss

abstract: Due to their high retrieval efficiency and low storage cost for cross-modal search task, cross-modal hashing methods have attracted considerable attention. For the supervised cross-modal hashing methods, how to make the learned hash codes preserve semantic information sufficiently contained in the label of datapoints is the key to further enhance the retrieval performance. Hence, almost all supervised cross-modal hashing methods usually depends on defining a similarity between datapoints with the label information to guide the hashing model learning fully or partly. However, the defined similarity between datapoints can only capture the label information of datapoints partially and misses abundant semantic information, then hinders the further improvement of retrieval performance. Thus, in this paper, different from previous works, we propose a novel cross-modal hashing method without defining the similarity between datapoints, called Deep Cross-modal Hashing via \textit{Margin-dynamic-softmax Loss} (DCHML). Specifically, DCHML first trains a proxy hashing network to transform each category information of a dataset into a semantic discriminative hash code, called proxy hash code. Each proxy hash code can preserve the semantic information of its corresponding category well. Next, without defining the similarity between datapoints to supervise the training process of the modality-specific hashing networks , we propose a novel \textit{margin-dynamic-softmax loss} to directly utilize the proxy hashing codes as supervised information. Finally, by minimizing the novel \textit{margin-dynamic-softmax loss}, the modality-specific hashing networks can be trained to generate hash codes which can simultaneously preserve the cross-modal similarity and abundant semantic information well.


title: Improving Sales Forecasting Accuracy: A Tensor Factorization Approach with Demand Awareness

abstract: Due to accessible big data collections from consumers, products, and stores, advanced sales forecasting capabilities have drawn great attention from many companies especially in the retail business because of its importance in decision making. Improvement of the forecasting accuracy, even by a small percentage, may have a substantial impact on companies' production and financial planning, marketing strategies, inventory controls, supply chain management, and eventually stock prices. Specifically, our research goal is to forecast the sales of each product in each store in the near future. Motivated by tensor factorization methodologies for personalized context-aware recommender systems, we propose a novel approach called the Advanced Temporal Latent-factor Approach to Sales forecasting (ATLAS), which achieves accurate and individualized prediction for sales by building a single tensor-factorization model across multiple stores and products. Our contribution is a combination of: tensor framework (to leverage information across stores and products), a new regularization function (to incorporate demand dynamics), and extrapolation of tensor into future time periods using state-of-the-art statistical (seasonal auto-regressive integrated moving-average models) and machine-learning (recurrent neural networks) models. The advantages of ATLAS are demonstrated on eight product category datasets collected by the Information Resource, Inc., where a total of 165 million weekly sales transactions from more than 1,500 grocery stores over 15,560 products are analyzed.


title: Relative relevance of mobility and driving force on edge dislocation climb by the vacancy mechanism

abstract: In this work we examine the driving force for edge dislocations to climb in $\alpha$-Fe from atomistic and mesoscale viewpoints. We study the bias for the climb process depending on the dislocation orientation and the applied stress as coming from both the gradient of the chemical potential and the transport coefficients. Both terms are modified by the applied stress and therefore contribute to climb. Surprisingly, even though the vacancy migration barrier distribution is modified by the external stress as obtained by nudged-elastic band calculations, the mobilities resulting from a kinetic Monte Carlo model applied on the obtained energy landscape are indistinguishable, independently of the stress. Moreover, an object kinetic Monte Carlo (OKMC) model including the effect of the dislocation strain field to first order shows indeed a slight anisotropic component in the diffusion in more complex dislocation configurations. However, the OKMC results highlight the fact that the thermodynamic component is the dominant driving force. We conclude that in $\alpha$-Fe under thermal conditions, the main source of bias is given by the difference in vacancy chemical potentials, which is small enough to hinder the process for dynamic atomistic simulations.


title: Fixed Parameter Approximation Scheme for Min-max $k$-cut

abstract: We consider the graph $k$-partitioning problem under the min-max objective, termed as Minmax $k$-cut. The input here is a graph $G=(V,E)$ with non-negative edge weights $w:E\rightarrow \mathbb{R}_+$ and an integer $k\geq 2$ and the goal is to partition the vertices into $k$ non-empty parts $V_1, \ldots, V_k$ so as to minimize $\max_{i=1}^k w(\delta(V_i))$. Although minimizing the sum objective $\sum_{i=1}^k w(\delta(V_i))$, termed as Minsum $k$-cut, has been studied extensively in the literature, very little is known about minimizing the max objective. We initiate the study of Minmax $k$-cut by showing that it is NP-hard and W[1]-hard when parameterized by $k$, and design a parameterized approximation scheme when parameterized by $k$. The main ingredient of our parameterized approximation scheme is an exact algorithm for Minmax $k$-cut that runs in time $(\lambda k)^{O(k^2)}n^{O(1)}$, where $\lambda$ is value of the optimum and $n$ is the number of vertices. Our algorithmic technique builds on the technique of Lokshtanov, Saurabh, and Surianarayanan (FOCS, 2020) who showed a similar result for Minsum $k$-cut. Our algorithmic techniques are more general and can be used to obtain parameterized approximation schemes for minimizing $\ell_p$-norm measures of $k$-partitioning for every $p\geq 1$.


title: Systoles and diameters of hyperbolic surfaces

abstract: In this article we explore the relationship between the systole and the diameter of closed hyperbolic orientable surfaces. We show that they satisfy a certain inequality, which can be used to deduce that their ratio has a (genus dependent) upper bound.


title: Radio counterparts of gamma-ray sources in the Cygnus region

abstract: The view of the gamma-ray universe is being continuously expanded by space high energy (HE) and ground based very-high energy (VHE) observatories. Yet, the angular resolution limitation still precludes a straightforward identification of these gamma-ray emitting sources. Radio observations are an effective tool for searching their possible counterparts at lower energies because the same population of relativistic electrons responsible for radio emission can also produce HE/VHE emission via inverse-Compton scattering. The Cygnus region is crowded by many gamma-ray sources, most of them remaining unidentified. In order to find possible counterparts to unidentified gamma-ray sources, we carried out a deep survey of the Cygnus region using the Giant Metrewave Radio Telescope at 610 MHz and 325 MHz. We did a detailed search for counterparts in the error circle of HE/VHE sources. We report 36 radio sources found in the error ellipse of 15 HE sources, and 11 in those of VHE sources. Eight sources have very steep radio spectral index alpha <-1.5, which are most likely to be pulsars and will be followed up for periodicity search. Such a significant number of pulsar candidates within the error circle of HE/VHE sources prompts fresh look at the energetics and efficacy of pulsars and pulsar wind nebulae in this context.


title: On the maximum order complexity of Thue-Morse and Rudin-Shapiro sequences along polynomial values

abstract: Both the Thue-Morse and Rudin-Shapiro sequences are not suitable sequences for cryptography since their expansion complexity is small and their correlation measure of order 2 is large. These facts imply that these sequences are highly predictable despite the fact that they have a large maximum order complexity. Sun and Winterhof (2019) showed that the Thue-Morse sequence along squares keeps a large maximum order complexity. Since, by Christol's theorem, the expansion complexity of this rarefied sequence is no longer bounded, this provides a potentially better candidate for cryptographic applications. Similar results are known for the Rudin-Shapiro sequence and more general pattern sequences. In this paper we generalize these results to any polynomial subsequence (instead of squares) and thereby answer an open problem of Sun and Winterhof. We conclude this paper by some open problems.


title: Semi-invariants of Binary Forms and Sylvester's Theorem

abstract: We obtain a combinatorial formula related to the shear transformation for semi-invariants of binary forms, which implies the classical characterization of semi-invariants in terms of a differential operator. Then, we present a combinatorial proof of an identity of Hilbert, which leads to a relation of Cayley on semi-invariants. This identity plays a crucial role in the original proof of Sylvester's theorem on semi-invariants in connection with the Gaussian coefficients. Moreover, we show that the additivity lemma of Pak and Panova which yields the strict unimodality of the Gaussian coefficients for $n,k \geq 8$ can be deduced from the ring property of semi-invariants.


title: Complex Query Answering with Neural Link Predictors

abstract: Neural link predictors are immensely useful for identifying missing edges in large scale Knowledge Graphs. However, it is still not clear how to use these models for answering more complex queries that arise in a number of domains, such as queries using logical conjunctions ($\land$), disjunctions ($\lor$) and existential quantifiers ($\exists$), while accounting for missing edges. In this work, we propose a framework for efficiently answering complex queries on incomplete Knowledge Graphs. We translate each query into an end-to-end differentiable objective, where the truth value of each atom is computed by a pre-trained neural link predictor. We then analyse two solutions to the optimisation problem, including gradient-based and combinatorial search. In our experiments, the proposed approach produces more accurate results than state-of-the-art methods -- black-box neural models trained on millions of generated queries -- without the need of training on a large and diverse set of complex queries. Using orders of magnitude less training data, we obtain relative improvements ranging from 8% up to 40% in Hits@3 across different knowledge graphs containing factual information. Finally, we demonstrate that it is possible to explain the outcome of our model in terms of the intermediate solutions identified for each of the complex query atoms. All our source code and datasets are available online, at https://github.com/uclnlp/cqd.


title: Threats and Opportunities: Blockchain Meets Quantum Computation

abstract: This article considered deficiencies of the flourishing blockchain technology manifested by the development of quantum computation. We show that the future blockchain technology would under constant threats from the following aspects: 1) Speed up the generation of nonces; 2) Faster searching for hash collisions; 3) Break the security of the classical encryption. We also demonstrate that incorporating some quantum properties into blockchain makes it more robust and more efficient. For example people can establish a quantum-security blockchain system that utilizes quantum key distribution (QKD), and quantum synchronization and detectable Byzantine agreement (DBA) can help the blockchain systems achieve faster consensus even if there exist a number of malicious nodes.


title: Defining rough sets as core-support pairs of three-valued functions

abstract: We answer the question what properties a collection $\mathcal{F}$ of three-valued functions on a set $U$ must fulfill so that there exists a quasiorder $\leq$ on $U$ such that the rough sets determined by $\leq$ coincide with the core--support pairs of the functions in $\mathcal{F}$. Applying this characterization, we give a new representation of rough sets determined by equivalences in terms of three-valued {\L}ukasiewicz algebras of three-valued functions.


title: A Comprehensive Comparison of Multi-Dimensional Image Denoising Methods

abstract: Filtering multi-dimensional images such as color images, color videos, multispectral images and magnetic resonance images is challenging in terms of both effectiveness and efficiency. Leveraging the nonlocal self-similarity (NLSS) characteristic of images and sparse representation in the transform domain, the block-matching and 3D filtering (BM3D) based methods show powerful denoising performance. Recently, numerous new approaches with different regularization terms, transforms and advanced deep neural network (DNN) architectures are proposed to improve denoising quality. In this paper, we extensively compare over 60 methods on both synthetic and real-world datasets. We also introduce a new color image and video dataset for benchmarking, and our evaluations are performed from four different perspectives including quantitative metrics, visual effects, human ratings and computational cost. Comprehensive experiments demonstrate: (i) the effectiveness and efficiency of the BM3D family for various denoising tasks, (ii) a simple matrix-based algorithm could produce similar results compared with its tensor counterparts, and (iii) several DNN models trained with synthetic Gaussian noise show state-of-the-art performance on real-world color image and video datasets. Despite the progress in recent years, we discuss shortcomings and possible extensions of existing techniques. Datasets and codes for evaluation are made publicly available at https://github.com/ZhaomingKong/Denoising-Comparison.


title: Extending Equational Monadic Reasoning with Monad Transformers

abstract: There is a recent interest for the verification of monadic programs using proof assistants. This line of research raises the question of the integration of monad transformers, a standard technique to combine monads. In this paper, we extend Monae, a Coq library for monadic equational reasoning, with monad transformers and we explain the benefits of this extension. Our starting point is the existing theory of modular monad transformers, which provides a uniform treatment of operations. Using this theory, we simplify the formalization of models in Monae and we propose an approach to support monadic equational reasoning in the presence of monad transformers. We also use Monae to revisit the lifting theorems of modular monad transformers by providing equational proofs and explaining how to patch a known bug using a non-standard use of Coq that combines impredicative polymorphism and parametricity.


title: HAVEN: A Unity-based Virtual Robot Environment to Showcase HRI-based Augmented Reality

abstract: Due to the COVID-19 pandemic, conducting Human-Robot Interaction (HRI) studies in person is not permissible due to social distancing practices to limit the spread of the virus. Therefore, a virtual reality (VR) simulation with a virtual robot may offer an alternative to real-life HRI studies. Like a real intelligent robot, a virtual robot can utilize the same advanced algorithms to behave autonomously. This paper introduces HAVEN (HRI-based Augmentation in a Virtual robot Environment using uNity), a VR simulation that enables users to interact with a virtual robot. The goal of this system design is to enable researchers to conduct HRI Augmented Reality studies using a virtual robot without being in a real environment. This framework also introduces two common HRI experiment designs: a hallway passing scenario and human-robot team object retrieval scenario. Both reflect HAVEN's potential as a tool for future AR-based HRI studies.


title: A search for variable subdwarf B stars in TESS Full Frame Images -- I. Variable objects in the southern ecliptic hemisphere

abstract: We report the results of our search for pulsating subdwarf B stars in Full Frame Images, sampled at 30 min cadence and collected during Year 1 of the TESS mission. Year 1 covers most of the southern ecliptic hemisphere. The sample of objects we checked for pulsations was selected from a subdwarf B stars database available to public. Only two positive detections have been achieved, however, as a by-product of our search we found 1807 variable objects, most of them not classified, hence their specific variability class cannot be confirmed at this stage. Our preliminary discoveries include: two new subdwarf B (sdB) pulsators, 26 variables with known sdB spectra, 83 non-classified pulsating stars, 83 eclipsing binaries (detached and semi-detached), a mix of 1535 pulsators and non-eclipsing binaries, two novae, and 77 variables with known (non-sdB) spectral classification. Among eclipsing binaries we identified two known HW Vir systems and four new candidates. The amplitude spectra of the two sdB pulsators are not rich in modes, but we derive estimates of the modal degree for one of them. In addition, we selected five sdBV candidates for mode identification among 83 pulsators and describe our results based on this preliminary analysis. Further progress will require spectral classification of the newly discovered variable stars, which hopefully include more subdwarf B stars.


title: Real-time noise cancellation with Deep Learning

abstract: Biological measurements are often contaminated with large amounts of non-stationary noise which require effective noise reduction techniques. We present a new real-time deep learning algorithm which produces adaptively a signal opposing the noise so that destructive interference occurs. As a proof of concept, we demonstrate the algorithm's performance by reducing electromyogram noise in electroencephalograms with the usage of a custom, flexible, 3D-printed, compound electrode. With this setup, an average of 4dB and a maximum of 10dB improvement of the signal-to-noise ratio of the EEG was achieved by removing wide band muscle noise. This concept has the potential to not only adaptively improve the signal-to-noise ratio of EEG but can be applied to a wide range of biological, industrial and consumer applications such as industrial sensing or noise cancelling headphones.


title: Nodal set of monochromatic waves satisfying the Random Wave Model

abstract: We construct deterministic solutions to the Helmholtz equation in $\mathbb{R}^m$ which behave accordingly to the Random Wave Model. We then find the number of their nodal domains, their nodal volume and the topologies and nesting trees of their nodal set in growing balls around the origin. The proof of the pseudo-random behaviour of the functions under consideration hinges on a de-randomisation technique pioneered by Bourgain and proceeds via computing their $L^p$-norms. The study of their nodal set relies on its stability properties and on the evaluation of their doubling index, in an average sense.


title: Numerical Dissipation Based Error Estimators and Grid Adaptation for Large Eddy Simulation

abstract: Grid adaptation for implicit Large Eddy Simulation (LES) is a non-trivial challenge due to the inherent coupling of the modeling and numerical errors. An attempt to address the challenge first requires a comprehensive assessment and then the development of error estimators to highlight regions that require refinement. Following the work of Schranner et al., a novel approach to estimate the numerical dissipation of the turbulent kinetic energy (TKE) equations is proposed. The presented approach allows the computation of the local numerical dissipation for arbitrary curvilinear grids through a post-processing procedure. This method, as well as empirical and kinetic-energy-based approaches, are employed to estimate the inherent numerical TKE. We incorporate the numerical TKE to evaluate an effective eddy viscosity, an effective Kolmogorov length scale, and an effective TKE to build a family of Index Quality (IQ) based error estimators. The proposed IQ based estimators are then assessed and utilized to show their effectiveness through an application of grid adaptation for the periodic hill test case and transitional flow over the SD 7003 airfoil. Numerical results are validated through a comparison against reference LES and experimental data. Flow over the adapted grids appear better abled to capture pertinent flow features and integrated functions, such as the lift and drag coefficients.


title: Understanding Pure Character-Based Neural Machine Translation: The Case of Translating Finnish into English

abstract: Recent work has shown that deeper character-based neural machine translation (NMT) models can outperform subword-based models. However, it is still unclear what makes deeper character-based models successful. In this paper, we conduct an investigation into pure character-based models in the case of translating Finnish into English, including exploring the ability to learn word senses and morphological inflections and the attention mechanism. We demonstrate that word-level information is distributed over the entire character sequence rather than over a single character, and characters at different positions play different roles in learning linguistic knowledge. In addition, character-based models need more layers to encode word senses which explains why only deeper models outperform subword-based models. The attention distribution pattern shows that separators attract a lot of attention and we explore a sparse word-level attention to enforce character hidden states to capture the full word-level information. Experimental results show that the word-level attention with a single head results in 1.2 BLEU points drop.


title: Transmission spectroscopy for the warm sub-Neptune HD3167c: evidence for molecular absorption and a possible high metallicity atmosphere

abstract: We present a transmission spectrum for the warm (500-600K) sub-Neptune HD3167c obtained using the Hubble Space Telescope Wide Field Camera 3 infrared spectrograph. We combine these data, which span the 1.125-1.643 micron wavelength range, with broadband transit measurements made using Kepler/K2 (0.6-0.9 micron) and Spitzer/IRAC (4-5 micron). We find evidence for absorption by at least one of H2O, HCN, CO2, and CH4 (Bayes factor 7.4; 2.5-sigma significance), although the data precision does not allow us to unambiguously discriminate between these molecules. The transmission spectrum rules out cloud-free hydrogen-dominated atmospheres with metallicities <100x solar at >5.8-sigma confidence. In contrast, good agreement with the data is obtained for cloud-free models assuming metallicities >700x solar. However, for retrieval analyses that include the effect of clouds, a much broader range of metallicities (including subsolar) is consistent with the data, due to the degeneracy with cloud-top pressure. Self-consistent chemistry models that account for photochemistry and vertical mixing are presented for the atmosphere of HD3167c. The predictions of these models are broadly consistent with our abundance constraints, although this is primarily due to the large uncertainties on the latter. Interior structure models suggest the core mass fraction is >40%, independent of a rock or water core composition, and independent of atmospheric envelope metallicity up to 1000x solar. We also report abundance measurements for fifteen elements in the host star, showing that it has a very nearly solar composition.


title: Accelerating combinatorial filter reduction through constraints

abstract: Reduction of combinatorial filters involves compressing state representations that robots use. Such optimization arises in automating the construction of minimalist robots. But exact combinatorial filter reduction is an NP-complete problem and all current techniques are either inexact or formalized with exponentially many constraints. This paper proposes a new formalization needing only a polynomial number of constraints, and characterizes these constraints in three different forms: nonlinear, linear, and conjunctive normal form. Empirical results show that constraints in conjunctive normal form capture the problem most effectively, leading to a method that outperforms the others. Further examination indicates that a substantial proportion of constraints remain inactive during iterative filter reduction. To leverage this observation, we introduce just-in-time generation of such constraints, which yields improvements in efficiency and has the potential to minimize large filters.


title: A note on the Jacobian Conjecture

abstract: Let $F:\Bbb C^n\to\Bbb C^n$ be a polynomial mapping with a non vanishing Jacobian. If the set $S_F$ of non-properness of $F$ is smooth, then $F$ is a surjective mapping. Moreover, the set $S_F$ can not be connected (this is the Nollet-Xavier Conjecture). Additionally, if $n=2$, then the set $S_F$ of non-properness of $F$ cannot be a curve without self-intersections.


title: Geometric conditions for saturating the data processing inequality

abstract: The data processing inequality (DPI) is a scalar inequality satisfied by distinguishability measures on density matrices. For some distinguishability measures, saturation of the scalar DPI implies an operator equation relating the arguments of the measure. These results are typically derived using functional analytic techniques. In a complementary approach, we use geometric techniques to derive a formula that gives an operator equation from DPI saturation for any distinguishability measure; moreover, for a broad class of distinguishability measures, the derived operator equation is sufficient to imply saturation as well. Our operator equation coincides with known results for the sandwiched R\'{e}nyi relative entropies, and gives new results for $\alpha$-$z$ R\'{e}nyi relative entropies and a family of of quantum $f$-divergences, which we compute explicitly.


title: Correlator Convolutional Neural Networks: An Interpretable Architecture for Image-like Quantum Matter Data

abstract: Machine learning models are a powerful theoretical tool for analyzing data from quantum simulators, in which results of experiments are sets of snapshots of many-body states. Recently, they have been successfully applied to distinguish between snapshots that can not be identified using traditional one and two point correlation functions. Thus far, the complexity of these models has inhibited new physical insights from this approach. Here, using a novel set of nonlinearities we develop a network architecture that discovers features in the data which are directly interpretable in terms of physical observables. In particular, our network can be understood as uncovering high-order correlators which significantly differ between the data studied. We demonstrate this new architecture on sets of simulated snapshots produced by two candidate theories approximating the doped Fermi-Hubbard model, which is realized in state-of-the art quantum gas microscopy experiments. From the trained networks, we uncover that the key distinguishing features are fourth-order spin-charge correlators, providing a means to compare experimental data to theoretical predictions. Our approach lends itself well to the construction of simple, end-to-end interpretable architectures and is applicable to arbitrary lattice data, thus paving the way for new physical insights from machine learning studies of experimental as well as numerical data.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.quant-gas cs.LG physics.comp-ph


title: Theoretical Characterization of Photoactive Molecular Systems Based on BODIPY-Derivatives for the Design of Organic Solar Cells

abstract: To search for high-efficiency narrow-band donor materials to improve the short-circuit current density ($J_{sc}$) of organic solar cells, a series of small molecules based in Bodipy-Triphenylamine were characterized using density functional theory (DFT) and time-dependent (TD-DFT) calculations. According to the energy of the exciton driving force they have the appropriate energy levels to match \textbf{\textit{PC$_{61}$BM}}. The properties affecting the open circuit voltage ($V_{oc}$), $J_{sc}$ and the fill factor ($FF$) were investigated by calculating the geometric structures, the boundary molecular orbital energy levels, absorption spectra, light collection efficiencies, chare transfer rates, and exciton binding energies. The results show that the \textbf{\textit{BTPA~III}} system has a lower LUMO level, high absorption efficiency, and exction dissociation than other molecular systems, facilitating the improvement of $V_{oc}$, $J_{sc}$ and $FF$. Finally, \textbf{\textit{BTPA~III}} would be the most promising of this series of donors and further increase the efficiency of the device.


title: Towards Obfuscated Malware Detection for Low Powered IoT Devices

abstract: With the increased deployment of IoT and edge devices into commercial and user networks, these devices have become a new threat vector for malware authors. It is imperative to protect these devices as they become more prevalent in commercial and personal networks. However, due to their limited computational power and storage space, especially in the case of battery-powered devices, it is infeasible to deploy state-of-the-art malware detectors onto these systems. In this work, we propose using and extracting features from Markov matrices constructed from opcode traces as a low cost feature for unobfuscated and obfuscated malware detection. We empirically show that our approach maintains a high detection rate while consuming less power than similar work.


title: Diffusion-redistanciation schemes for 2D and 3D constrained Willmore flow: application to the equilibrium shapes of vesicles

abstract: In this paper we present a novel algorithm for simulating geometrical flows, and in particular the Willmore flow, with conservation of volume and area. The idea is to adapt the class of diffusion-redistanciation algorithms to the Willmore flow in both two and three dimensions. These algorithms rely on alternating diffusions of the signed distance function to the interface and a redistanciation step, and with careful choice of the applied diffusions, end up moving the zero level-set of the distance function by some geometrical quantity without resorting to any explicit transport equation. The constraints are enforced between the diffusion and redistanciation steps via a simple rescaling method. The energy globally decreases at the end of each global step. The algorithms feature the computational efficiency of thresholding methods without requiring any adaptive remeshing thanks to the use of a signed distance function to describe the interface. This opens their application to dynamic fluid-structure simulations for large and realistic cases. The methodology is validated by computing the equilibrium shapes of two- and three-dimensional vesicles, as well as the Clifford torus.


title: The future of astronomy with small satellites

abstract: The number of small satellites has grown dramatically in the past decade from tens of satellites per year in the mid-2010s to a projection of tens of thousands in orbit by the mid-2020s. This presents both problems and opportunities for observational astronomy. Small satellites offer complementary cost-effective capabilities to both ground-based astronomy and larger space missions. Compared to ground-based astronomy, these advantages are not just in the accessibility of wavelength ranges where the Earth's atmosphere is opaque, but also in stable, high precision photometry, long-term monitoring and improved areal coverage. Astronomy has a long history of new observational parameter spaces leading to major discoveries. Here we discuss the potential for small satellites to explore new parameter spaces in astrophysics, drawing on examples from current and proposed missions, and spanning a wide range of science goals from binary stars, exoplanets and solar system science to the early Universe and fundamental physics.

categories: astro-ph.IM astro-ph.EP astro-ph.GA astro-ph.HE astro-ph.SR


title: Massively Parallel Graph Drawing and Representation Learning

abstract: To fully exploit the performance potential of modern multi-core processors, machine learning and data mining algorithms for big data must be parallelized in multiple ways. Today's CPUs consist of multiple cores, each following an independent thread of control, and each equipped with multiple arithmetic units which can perform the same operation on a vector of multiple data objects. Graph embedding, i.e. converting the vertices of a graph into numerical vectors is a data mining task of high importance and is useful for graph drawing (low-dimensional vectors) and graph representation learning (high-dimensional vectors). In this paper, we propose MulticoreGEMPE (Graph Embedding by Minimizing the Predictive Entropy), an information-theoretic method which can generate low and high-dimensional vectors. MulticoreGEMPE applies MIMD (Multiple Instructions Multiple Data, using OpenMP) and SIMD (Single Instructions Multiple Data, using AVX-512) parallelism. We propose general ideas applicable in other graph-based algorithms like \emph{vectorized hashing} and \emph{vectorized reduction}. Our experimental evaluation demonstrates the superiority of our approach.


title: The non-orientable 4-genus for knots with 10 crossings

abstract: Given a knot in the 3-sphere, the non-orientable 4-genus or 4-dimensional crosscap number of a knot is the minimal first Betti number of non-orientable surfaces, smoothly and properly embedded in the 4-ball, with boundary the knot. In this paper, we calculate the non-orientable 4-genus of knots with crossing number 10.


title: Sublinearly Morse Boundary II: Proper geodesic spaces

abstract: We build an analogue of the Gromov boundary for any proper geodesic metric space, hence for any finitely generated group. More precisely, for any proper geodesic metric space $X$ and any sublinear function $\kappa$, we construct a boundary for $X$, denoted $\mathcal{\partial}_{\kappa} X$, that is quasi-isometrically invariant and metrizable. As an application, we show that when $G$ is the mapping class group of a finite type surface, or a relatively hyperbolic group, then with minimal assumptions the Poisson boundary of $G$ can be realized on the $\kappa$-Morse boundary of $G$ equipped the word metric associated to any finite generating set.


title: Detecting spatial clusters in functional data: new scan statistic approaches

abstract: We have developed two scan statistics for detecting clusters of functional data indexed in space. The first method is based on an adaptation of a functional analysis of variance and the second one is based on a distribution-free spatial scan statistic for univariate data. In a simulation study, the distribution-free method always performed better than a nonparametric functional scan statistic, and the adaptation of the anova also performed better for data with a normal or a quasi-normal distribution. Our methods can detect smaller spatial clusters than the nonparametric method. Lastly, we used our scan statistics for functional data to search for spatial clusters of abnormal unemployment rates in France over the period 1998-2013 (divided into quarters).


title: BICEP / Keck XII: Constraints on axion-like polarization oscillations in the cosmic microwave background

abstract: We present a search for axion-like polarization oscillations in the cosmic microwave background (CMB) with observations from the Keck Array. A local axion field induces an all-sky, temporally sinusoidal rotation of CMB polarization. A CMB polarimeter can thus function as a direct-detection experiment for axion-like dark matter. We develop techniques to extract an oscillation signal. Many elements of the method are generic to CMB polarimetry experiments and can be adapted for other datasets. As a first demonstration, we process data from the 2012 observing season to set upper limits on the axion-photon coupling constant in the mass range $10^{-21}$-$10^{-18}~\mathrm{eV}$, which corresponds to oscillation periods on the order of hours to months. We find no statistically significant deviations from the background model. For periods larger than $24~\mathrm{hr}$ (mass $m < 4.8 \times 10^{-20}~\mathrm{eV}$), the median 95%-confidence upper limit is equivalent to a rotation amplitude of $0.68^\circ$, which constrains the axion-photon coupling constant to $g_{\phi\gamma} < \left ( 1.1 \times 10^{-11}~\mathrm{GeV}^{-1} \right ) m/\left (10^{-21}~\mathrm{eV} \right )$, if axion-like particles constitute all of the dark matter. The constraints can be improved substantially with data already collected by the BICEP series of experiments. Current and future CMB polarimetry experiments are expected to achieve sufficient sensitivity to rule out unexplored regions of the axion parameter space.


title: GLISSE: A GPU-optimized planetary system integrator with application to orbital stability calculations

abstract: We present a GPU-accelerated numerical integrator specifically optimized for stability calculations of small bodies in planetary systems. Specifically, the integrator is designed for cases when large numbers of test particles (tens or hundreds of thousands) need to be followed for long durations (millions of orbits) to assess the orbital stability of their initially "close-encounter free" orbits. The GLISSE (Gpu Long-term Integrator for Solar System Evolution) code implements several optimizations to achieve a roughly factor of 100 speed increase over running the same code on a CPU. We explain how various hardware speed bottlenecks can be avoided by the careful code design, although some of the choices restrict the usage to specific types of application. As a first application, we study the long-term stability of small bodies initially on orbits between Uranus and Neptune. We map out in detail the small portion of the phase space in which small bodies can survive for 4.5 billion years of evolution; the ability to integrate large numbers of particles allow us to identify for the first time how instability-inducing mean-motion resonance overlaps sharply define the stable regions. As a second application, we map the boundaries of 4 Gyr stability or transneptunian objects in the 5:2 and 3:1 mean-motion resonances, demonstrating that long-term perturbations remove the initially stable Neptune-crossing members.


title: Enhanced decoherence for a neutral particle sliding on a metallic surface in vacuum

abstract: Bodies in relative motion, spatially separated in vacuum, experience a tiny friction force known as quantum friction. This force has eluded experimental detection so far due to its small magnitude and short range. Herein, we give quantitative details so as to track traces of the quantum friction by measuring coherences in the atom. We notice that the environmentally induced decoherence can be decomposed into contributions of different signature: corrections induced by the electromagnetic vacuum in presence of the dielectric sheet and those induced by the motion of the particle. In this direction, we show that non-contact friction enhances the decoherence of the moving atom. Further, its effect can be enlarged by a thorough selection of the two-level particle and the Drude-Lorentz parameters of the material. In this context, we suggest that measuring decoherence times through velocity dependence of coherences could indirectly demonstrate the existence of quantum friction.


title: Hawking radiation of massive fields in 2D

abstract: We extend the classic results of the paper P. C. W. Davies, S. A. Fulling, and W. G. Unruh "Energy-momentum tensor near an evaporating black hole" by considering a massive scalar field in a two dimensions in the presence of a thin shell collapse. We show that outside the shell the WKB approximation is valid for any value of $r$ if $mr_{g} \gg 1$, where $m$ is the mass of the field, and $r_{g}$ is the Schwarzschild radius. Thus, we use semiclassical modes to calculate the flux in the vicinity of the shell, and at spatial infinity, $r \rightarrow +\infty$ at the final stage of the collapse, $t \rightarrow +\infty$ with the use of the covariant point-splitting regularization. We get that near the shell and at the spatial infinity the radiation is thermal with Hawking temperature. We obtain the negative flux $T_{vv}$ in the vicinity of the shell, which is similar to the classic result in the massless case.


title: Dynamics of primordial fields in quantum cosmological spacetimes

abstract: Quantum cosmological models are commonly described by means of semiclassical approximations in which a smooth evolution of the expectation values of elementary geometry operators replaces the classical and singular dynamics. The advantage of such descriptions is that they are relatively simple and display the classical behavior for large universes. However, they may smooth out an important inner structure and to include it a more detailed treatment is needed. The purpose of the present work is to investigate quantum uncertainty in the basic background variables and its influence on primordial gravitational waves. To this end we quantize a model of the Friedmann-Lemaitre-Robertson-Walker universe filled with a linear barotropic cosmological fluid and with gravitational waves. We carefully derive the dynamical equations for the perturbations in quantum spacetime. The quantization yields an equation of motion for the Fourier modes of gravitational radiation, which is a quantum extension to the usual parametric oscillator equation for gravitational waves propagating in an expanding universe. The two quantum effects from the cosmological background that enter the enhanced equation of motion are (i) a repulsive potential resolving the big bang singularity and replacing it with a big bounce and (ii) uncertainties in the numerical values for the background spacetime dynamical variables. First we study the former effect and its consequences for the primordial amplitude spectrum and carefully discuss the relation between the bounce scale and the physical predictions of the model. Next we investigate the latter effect, in particular the extent to which it may affect the primordial amplitude of gravitational waves. Making use of the WKB approximation we find an analytical formula for the amplitude spectrum as a function of the quantum dispersion of the background spacetime.


title: Learning with Molecules beyond Graph Neural Networks

abstract: We demonstrate a deep learning framework which is inherently based in the highly expressive language of relational logic, enabling to, among other things, capture arbitrarily complex graph structures. We show how Graph Neural Networks and similar models can be easily covered in the framework by specifying the underlying propagation rules in the relational logic. The declarative nature of the used language then allows to easily modify and extend the propagation schemes into complex structures, such as the molecular rings which we choose for a short demonstration in this paper.


title: A Dust Trap in the Young Multiple System HD 34700

abstract: Millimeter observations of disks around young stars reveal substructures indicative of gas pressure traps that may aid grain growth and planet formation. We present Submillimeter Array observations of HD 34700- two Herbig Ae stars in a close binary system (Aa/Ab, $\sim$0.25 AU), surrounded by a disk presenting a large cavity and spiral arms seen in scattered light, and two distant, lower mass companions. These observations include 1.3 mm continuum emission and the $^{12}$CO 2-1 line at $\sim0.5$" (178 AU) resolution. They resolve a prominent azimuthal asymmetry in the continuum, and Keplerian rotation of a circumbinary disk in the $^{12}$CO line. The asymmetry is located at a radius of $155^{+11}_{-7}$ AU, consistent with the edge of the scattered light cavity, being resolved in both radius ($72 ^{+14}_{-15}$ AU) and azimuth (FWHM = $64 ^{\circ +8}_{-7}$). The strong asymmetry in millimeter continuum emission could be evidence for a dust trap, together with the more symmetric morphology of $^{12}$CO emission and small grains. We hypothesize an unseen circumbinary companion, responsible for the cavity in scattered light and creating a vortex at the cavity edge that manifests in dust trapping. The disk mass has limitations imposed by the detection of $^{12}$CO and non-detection of $^{13}$CO. We discuss its consequences for the potential past gravitational instability of this system, likely accounting for the rapid formation of a circumbinary companion. We also report the discovery of resolved continuum emission associated with HD 34700B (projected separation $\sim1850$AU), which we explain through a circumstellar disk.


title: Physics-informed Neural-Network Software for Molecular Dynamics Applications

abstract: We have developed a novel differential equation solver software called PND based on the physics-informed neural network for molecular dynamics simulators. Based on automatic differentiation technique provided by Pytorch, our software allows users to flexibly implement equation of atom motions, initial and boundary conditions, and conservation laws as loss function to train the network. PND comes with a parallel molecular dynamics (MD) engine in order for users to examine and optimize loss function design, and different conservation laws and boundary conditions, and hyperparameters, thereby accelerate the PINN-based development for molecular applications.