title: A Probabilistic Domain-knowledge Framework for Nosocomial Infection Risk Estimation of Communicable Viral Diseases in Healthcare Personnel: A Case Study for COVID-19

abstract: Hospital-acquired infections of communicable viral diseases (CVDs) are posing a tremendous challenge to healthcare workers globally. Healthcare personnel (HCP) is facing a consistent risk of hospital-acquired infections, and subsequently higher rates of morbidity and mortality. We proposed a domain knowledge-driven infection risk model to quantify the individual HCP and the population-level healthcare facility risks. For individual-level risk estimation, a time-variant infection risk model is proposed to capture the transmission dynamics of CVDs. At the population-level, the infection risk is estimated using a Bayesian network model constructed from three feature sets including individual-level factors, engineering control factors, and administrative control factors. The sensitivity analyses indicated that the uncertainty in the individual infection risk can be attributed to two variables: the number of close contacts and the viral transmission probability. The model validation was implemented in the transmission probability model, individual level risk model, and population-level risk model using a Coronavirus disease case study. Regarding the first, multivariate logistic regression was applied for a cross-sectional data in the UK with an AIC value of 7317.70 and a 10-fold cross validation accuracy of 78.23%. For the second model, we collected laboratory-confirmed COVID-19 cases of HCP in different occupations. The occupation-specific risk evaluation suggested the highest-risk occupations were registered nurses, medical assistants, and respiratory therapists, with estimated risks of 0.0189, 0.0188, and 0.0176, respectively. To validate the population-level risk model, the infection risk in Texas and California was estimated. The proposed model will significantly influence the PPE allocation and safety plans for HCP


title: Non-dimensional Newton-Puiseux Expansions

abstract: Recent results in the theory and application of Newton-Puiseux expansions, i.e. fractional power series solutions of equations, suggest further developments within a more abstract algebraic-geometric framework, involving in particular the theory of toric varieties and ideals. Here, we present a number of such developments, especially in relation to the equations of van der Pol, Riccati, and Schr\"{o}dinger. Some pure mathematical concepts we are led to are Graver, Gr\"{o}bner, lattice and circuit bases, combinatorial geometry and differential algebra, and algebraic-differential equations. Two techniques are coordinated: classical dimensional analysis (DA) in applied mathematics and science, and a polynomial and differential-polynomial formulation of asymptotic expansions, referred to here as Newton-Puiseux (NP) expansions. The latter leads to power series with rational exponents, which depend on the choice of the dominant part of the equations, often referred to as the method of "dominant balance". The paper shows, using a new approach to DA based on toric ideals, how dimensionally homogeneous equations may best be non-dimensionalised, and then, with examples involving differential equations, uses recent work on NP expansions to show how non-dimensional parameters affect the results. Our approach finds a natural home within computational algebraic geometry, i.e. at the interface of abstract and algorithmic mathematics.


title: A proposal to integrate Data Envelopment Analysis and Le Chatelier principle

abstract: The article aims to estimate the short and long run efficient production frontiers for firm operating in an industry. The Le Chatelier principle represents the theoretical framework, while the tool used to find the frontiers is the non-parametric Data Envelopment Analysis technique. The proposal tries to overcome the main limitations that other efforts done in this context seem to have.


title: A framework for comprehensible multi-modal detection of cyber threats

abstract: Detection of malicious activities in corporate environments is a very complex task and much effort has been invested into research of its automation. However, vast majority of existing methods operate only in a narrow scope which limits them to capture only fragments of the evidence of malware's presence. Consequently, such approach is not aligned with the way how the cyber threats are studied and described by domain experts. In this work, we discuss these limitations and design a detection framework which combines observed events from different sources of data. Thanks to this, it provides full insight into the attack life cycle and enables detection of threats that require this coupling of observations from different telemetries to identify the full scope of the incident. We demonstrate applicability of the framework on a case study of a real malware infection observed in a corporate network.


title: Direct calculation of the ZZ-interaction rates in the multi-mode circuit-QED

abstract: Hamiltonians of the superconducting qubits of Transmon type involve non-zero ZZ-interaction terms due to their finite and small anharmonicities. These terms might lead to the unwanted accumulation of spurious phases during the execution of the two-qubit gates. Exact calculation of the ZZ-interaction rates requires the full diagonalization of the circuit Hamiltonians which very quickly becomes computationally demanding as the number of the modes in the coupler circuit increases. Here we propose a direct analytical method for the accurate estimation of the ZZ-interaction rates between low-anharmonicity qubits in the dispersive limit of the multi-mode circuit-QED. We observe very good agreement between the predictions of our method and the measurement data collected from the multi-qubit devices. Our method being an extension of our previous work in [1] is a new addition to the toolbox of the quantum microwave engineers as it relates the ZZ-interaction rates directly to the entries of the impedance matrix defined between the qubit ports.


title: On a Proof of the Strong Law of Large Numbers

abstract: Inspired by a recent preprint of N. Curien, we provided what may be a new and elementary proof of the Law of Large Numbers.


title: Dependence of magnetic domain patterns on plasma-induced differential oxidation of CoPd thin films

abstract: We demonstrate the evolution of the micro-patterned magnetic domains in CoPd thin films pretreated with e-beam lithography and O2 plasma. During the days-long oxidation, significantly different behaviors of the patterned magnetic domains under magnetization reversal are observed via magneto-optic Kerr effect microscopy on different days. The evolution of the magnetic behaviors indicate critical changes in the local magnetic anisotropy energies due to the Co oxides that evolve into different oxide forms, which are characterized by micro-area X-ray absorption spectroscopy and X-ray photoelectron spectroscopy. The coercive field of the area pre-exposed to plasma can decrease to a value 10 Oe smaller than that unexposed to plasma, whereas after a longer duration of oxidation the coercive field can instead become larger in the area pre-exposed to plasma than that unexposed, leading to an opposite magnetic pattern. Various forms of oxidation can therefore provide an additional dimension for magnetic-domain engineering to the current conventional lithographies.


title: Robust near-diagonal Green function estimates

abstract: We prove sharp near-diagonal pointwise bounds for the Green function $G_\Omega(x,y)$ for nonlocal operators of fractional order $\alpha \in (0,2)$. The novelty of our results is two-fold: the estimates are robust as $\alpha \to 2-$ and we prove the bounds without making use of the Dirichlet heat kernel $p_\Omega(t;x,y)$. In this way we can cover cases, in which the Green function satisfies isotropic bounds but the heat kernel does not.


title: C60 Cation as the Carrier of the 9577 and 9632 Angstrom Diffuse Interstellar Bands: Further Support from the VLT/X-Shooter Spectra

abstract: Ever since their first detection over 100 years ago, the mysterious diffuse interstellar bands (DIBs), a set of several hundred broad absorption features seen against distant stars in the optical and near infrared wavelength range, largely remain unidentified. The close match both in wavelengths and in relative strengths recently found between the experimental absorption spectra of gas-phase buckminsterfullerene ions (C60+) and four DIBs at 9632, 9577, 9428 and 9365 Angstrom (and, to a lesser degree, a weaker DIB at 9348 Angstrom) suggests C60+ as a promising carrier. However, arguments against the C60+ identification remain and are mostly concerned with the large variation in the intensity ratios of the 9632 and 9577 DIBs. In this work, we search for these DIBs in the ESO VLT/X-shooter archival data and identify the 9632, 9577, 9428 and 9365 Angstrom DIBs in a sample of 25 stars. While the 9428 and 9365 Angstrom DIBs are too noisy to allow any reliable analysis, the 9632 and 9577 Angstrom DIBs are unambiguously detected and, after correcting for telluric water vapor absorption, their correlation can be used to probe their origin. To this end, we select a sub-sample of nine hot, O- or B0-type stars of which the stellar Mg II contamination to the 9632 Angstrom DIB is negligibly small. We find their equivalent widths, after normalized by reddening to eliminate their common correlation with the density of interstellar clouds, exhibit a tight, positive correlation. This supports C60+ as the carrier of the 9632 and 9577 Angstrom DIBs.

categories: astro-ph.GA cond-mat.mtrl-sci physics.atm-clus physics.chem-ph


title: Symbolic Security Predicates: Hunt Program Weaknesses

abstract: Dynamic symbolic execution (DSE) is a powerful method for path exploration during hybrid fuzzing and automatic bug detection. We propose security predicates to effectively detect undefined behavior and memory access violation errors. Initially, we symbolically execute program on paths that don't trigger any errors (hybrid fuzzing may explore these paths). Then we construct a symbolic security predicate to verify some error condition. Thus, we may change the program data flow to entail null pointer dereference, division by zero, out-of-bounds access, or integer overflow weaknesses. Unlike static analysis, dynamic symbolic execution does not only report errors but also generates new input data to reproduce them. Furthermore, we introduce function semantics modeling for common C/C++ standard library functions. We aim to model the control flow inside a function with a single symbolic formula. This assists bug detection, speeds up path exploration, and overcomes overconstraints in path predicate. We implement the proposed techniques in our dynamic symbolic execution tool Sydr. Thus, we utilize powerful methods from Sydr such as path predicate slicing that eliminates irrelevant constraints. We present Juliet Dynamic to measure dynamic bug detection tools accuracy. The testing system also verifies that generated inputs trigger sanitizers. We evaluate Sydr accuracy for 11 CWEs from Juliet test suite. Sydr shows 95.59% overall accuracy. We make Sydr evaluation artifacts publicly available to facilitate results reproducibility.


title: A classification of nonexpansive Bratteli-Vershik systems

abstract: We study simple, properly ordered nonexpansive Bratteli-Vershik ($BV$) systems. Correcting a mistake in an earlier paper, we redefine the classes standard nonexpansive ($SNE$) and strong standard nonexpansive ($SSNE$). We define also the classes of very well timed and well timed systems, their opposing classes of untimed and very untimed systems (which feature, as subclasses of "Case (2)", in the work of Downarowicz and Maass as well as Hoynes on expansiveness of $BV$ systems of finite topological rank), and several related classes according to the existence of indistinguishable pairs (of some "depth") and their synchronization ("common cuts"). We establish some properties of these types of systems and some relations among them. We provide several relevant examples, including a problematic one that is conjugate to a well timed system while also (vacuously) in the classes "Case (2)". We prove that the class of all simple, properly ordered nonexpansive $BV$ systems is the disjoint union of the ones conjugate to well timed systems and those conjugate to untimed systems, thereby showing that nonexpansiveness in $BV$ systems arises in one of two mutually exclusive ways.


title: The double smeared null energy condition

abstract: The null energy condition (NEC), an important assumption of the Penrose singularity theorem, is violated by quantum fields. The natural generalization of the NEC in quantum field theory, the renormalized null energy averaged over a finite null segment, is known to be unbounded from below. Here, we propose an alternative, the double smeared null energy condition (DSNEC), stating that the null energy smeared over two null directions has a finite lower bound. We rigorously derive DSNEC from general worldvolume bounds for free quantum fields in Minkowski spacetime. Our method allows for future systematic inclusion of curvature corrections. As a further application of the techniques we develop, we prove additional lower bounds on the expectation values of various operators such as conserved higher spin currents. DSNEC provides a natural starting point for proving singularity theorems in semi-classical gravity.


title: Alexandrov-Fenchel inequality for convex hypersurfaces with capillary boundary in a ball

abstract: In this paper, we first introduce the quermassintegrals for convex hypersurfaces with capillary boundary in the unit Euclidean ball $\mathbb{B}^{n+1}$ and derive its first variational formula. Then by using a locally constrained nonlinear curvature flow, which preserves the $n$-th quermassintegral and non-decreases the $k$-th quermassintegral, we obtain the Alexandrov-Fenchel inequality for convex hypersurfaces with capillary boundary in $\mathbb{B}^{n+1}$. This generalizes the result in \cite{SWX} for convex hypersurfaces with free boundary in $\mathbb{B}^{n+1}$.


title: Optimal Discrete Morse Theory Simplification (Expository Survey)

abstract: A central problem in topological data analysis is that of computing the homology of a given simplicial complex. Said complexes can have arbitrary large number of simplices, as can happen, for example, if the space is the Rips-Vietoris or Cech complex of a large data cloud. Thus, pre-processing the simplicial complex to get a smaller complex with the same homology groups and then applying the homology algorithm to the smaller one, has been an active research topic in the last years. In this survey, we discuss some recent papers that examine the complexity of this simplification via Discrete Morse Theory. This survey was prepared as a final project for a course on Computational Topology at The Ohio State University.


title: Data Compression: Multi-Term Approach

abstract: In terms of signal samples, we propose and justify a new rank reduced multi-term transform, abbreviated as MTT, which, under certain conditions, may provide better-associated accuracy than that of known optimal rank reduced transforms. The basic idea is to construct the transform with more parameters to optimize than those in the known optimal transforms. This is realized by the extension of the known transform structures to the form that includes additional terms - the MTT has four matrices to minimize the cost. The MTT structure has also a special transformation that decreases the numerical load. As a result, the MTT performance is improved by the variation of the MTT components.


title: Spin Hurwitz theory and Miwa transform for the Schur Q-functions

abstract: Schur functions are the common eigenfunctions of generalized cut-and-join operators which form a closed algebra. They can be expressed as differential operators in time-variables and also through the eigenvalues of auxiliary $N\times N$ matrices $X$, known as Miwa variables. Relevant for the cubic Kontsevich model and also for spin Hurwitz theory is an alternative set of Schur Q-functions. They appear in representation theory of the Sergeev group, which is a substitute of the symmetric group, related to the queer Lie superalgebras $\mathfrak{q}(N)$.. The corresponding spin $\hat{\cal W}$-operators were recently found in terms of time-derivatives, but a substitute of the Miwa parametrization remained unknown, which is an essential complication for the matrix model technique and further developments. We demonstrate that the Miwa representation, in this case, involves a fermionic matrix $\Psi$ in addition to $X$, but its realization using supermatrices is {\it not} quite naive.


title: Power-of-two sampling in redundancy systems: the impact of assignment constraints

abstract: A classical sampling strategy for load balancing policies is power-of-two, where any server pair is sampled with equal probability. This does not cover practical settings with assignment constraints which force non-uniform sampling. While intuition suggests that non-uniform sampling adversely impacts performance, this was only supported through simulations, and rigorous statements have remained elusive. Building on product-form distributions for redundancy systems, we prove the stochastic dominance of uniform sampling for a four-server system as well as arbitrary-size systems in light traffic.


title: Theoretical and Empirical Analysis of a Fast Algorithm for Extracting Polygons from Signed Distance Bounds

abstract: Recently there has been renewed interest in signed distance bound representations due to their unique properties for 3D shape modelling. This is especially the case for deep learning-based bounds. However, it is beneficial to work with polygons in most computer-graphics applications. Thus, in this paper we introduce and investigate an asymptotically fast method for transforming signed distance bounds into polygon meshes. This is achieved by combining the principles of sphere tracing (or ray marching) with traditional polygonization techniques, such as Marching Cubes. We provide theoretical and experimental evidence that this approach is of the $O(N^2\log N)$ computational complexity for a polygonization grid with $N^3$ cells. The algorithm is tested on both a set of primitive shapes as well as signed distance bounds generated from point clouds by machine learning (and represented as neural networks). Given its speed, implementation simplicity and portability, we argue that it could prove useful during the modelling stage as well as in shape compression for storage. The code is available here: https://github.com/nenadmarkus/gridhopping


title: Universality of Hofstadter butterflies on hyperbolic lattices

abstract: Motivated by recent experimental breakthroughs in realizing hyperbolic lattices in superconducting waveguides and electric circuits, we compute the Hofstadter butterfly on regular hyperbolic tilings. By utilizing large hyperbolic lattices with periodic boundary conditions, we obtain the true hyperbolic bulk spectrum that is unaffected by contributions from boundary states. Our results reveal that the butterfly spectrum with large extended gapped regions prevails and that its shape is universally determined by the number of edges of the fundamental tile, while the fractal structure is lost in such a non-Euclidean case. We explain how these intriguing features are related to the nature of Landau levels in hyperbolic space, and how they could be verified experimentally.

categories: cond-mat.mes-hall cond-mat.other cond-mat.stat-mech math-ph math.MP


title: Approximating bottleneck spanning trees on partitioned tuples of points

abstract: We present approximation algorithms for the following NP-hard optimization problems related to bottleneck spanning trees in metric spaces. 1. The disjoint bottleneck spanning tree problem: Given $n$ pairs of points in a metric space, find two disjoint trees each containing exactly one point from each pair and minimize the largest edge length (over all edges of both trees). It is known that approximating this problem by a factor better than 2 is NP-hard. We present a 4-approximation algorithm for this problem. This improves upon the previous best known approximation ratio of $9$. Our algorithm extends to a $(3k-2)$-approximation for a more general case where points are partitioned into $k$-tuples and we seek $k$ disjoint trees. 2. The generalized bottleneck spanning tree problem: Given $n$ points in some metric space that are partitioned into clusters of size at most 2, find a tree that contains exactly one point from each cluster and minimizes the largest edge length. We show that it is NP-hard to approximate this problem by a factor better than 2, and present a 3-approximation algorithm. 3. The partitioned bottleneck spanning tree problem: Given $kn$ points in some metric space, find $k$ trees each containing exactly $n$ points and minimize the largest edge length (over all edges of the $k$ trees). We show that it is NP-hard to approximate this problem by a factor better than 2 for any $k\ge 2$. We present an $\alpha$-approximation algorithm for this problem where $\alpha=2$ for $k=2,3$ and $\alpha=3$ for $k\ge 4$. Towards obtaining these approximation ratios we present tight upper bounds on the edge lengths of $k$ equal-size disjoint trees that can be obtained from the nodes of a given tree. This result is of independent interest. Our hardness proofs imply that it is NP-hard to approximate the non-metric version of the above problems within any constant factor.


title: MAJORCA: Multi-Architecture JOP and ROP Chain Assembler

abstract: Nowadays, exploits often rely on a code-reuse approach. Short pieces of code called gadgets are chained together to execute some payload. Code-reuse attacks can exploit vulnerabilities in the presence of operating system protection that prohibits data memory execution. The ROP chain construction task is the code generation for the virtual machine defined by an exploited executable. It is crucial to understand how powerful ROP attacks can be. Such knowledge can be used to improve software security. We implement MAJORCA that generates ROP and JOP payloads in an architecture agnostic manner and thoroughly consider restricted symbols such as null bytes that terminate data copying via strcpy. The paper covers the whole code-reuse payloads construction pipeline: cataloging gadgets, chaining them in DAG, scheduling, linearizing to the ready-to-run payload. MAJORCA automatically generates both ROP and JOP payloads for x86 and MIPS. MAJORCA constructs payloads respecting restricted symbols both in gadget addresses and data. We evaluate MAJORCA performance and accuracy with rop-benchmark and compare it with open-source compilers. We show that MAJORCA outperforms open-source tools. We propose a ROP chaining metric and use it to estimate the probabilities of successful ROP chaining for different operating systems with MAJORCA as well as other ROP compilers to show that ROP chaining is still feasible. This metric can estimate the efficiency of OS defences.


title: Geometric nature of relations on plabic graphs and totally non-negative Grassmannians

abstract: The standard parametrization of totally non-negative Grassmannians was obtained by A. Postnikov [45] introducing the boundary measurement map in terms of discrete path integration on planar bicolored (plabic) graphs in the disk. An alternative parametrization was proposed by T. Lam [38] introducing systems of relations on vectors on such graphs, depending on some signatures defined on edges. The problem of characterizing the signatures corresponding to the totally non-negative cells, was left open in [38]. In our paper we provide an explicit construction of such signatures, satisfying both the full rank condition and the total non-negativity property on the full positroid cell. If the graph $\mathcal G$ satisfies the following natural constraint: each edge belongs to some oriented path from the boundary to the boundary, then such signature is unique up to a vertex gauge transformation. Such signature is uniquely identified by geometric indices (local winding and intersection number) ruled by the orientation $\mathcal O$ and gauge ray direction $\mathfrak l$ on $\mathcal G$. Moreover, we provide a combinatorial representation of geometric signatures by showing that the total signature of every finite face just depends on the number of white vertices on it. The latter characterization is a Kasteleyn-type property [7,1] and we conjecture a mechanical-statistical interpretation of such relations. An explicit connection between the solution of Lam system of relations and the value of Postnikov boundary measurement map is established using the generalization of Talaska formula [50] obtained in [6]. In particular, the components of the edge vectors are rational in the edge weights with subtraction-free denominators. Finally, we provide explicit formulas for transformations of signatures under Postnikov moves and reductions, and amalgamations of networks.


title: The Local Economic Impact of Mineral Mining in Africa: Evidence from Four Decades of Satellite Imagery

abstract: Using state-of-the-art techniques in computer vision, we analyze one million satellite images covering 12% of the African continent between 1984 and 2019 to track local development around 1,658 mineral deposits. We use stacked event studies and difference-in-difference models to estimate the impact of mine openings and closings. The magnitude of the effect of mine openings is considerable - after 15 years, urban areas within 20km of an opening mine almost double in size. We find strong evidence of a political resource curse at the local level. Although mining boosts the local economy in democratic countries, these gains are meager in autocracies and come at the expense of tripling the likelihood of conflict relative to prior the onset of mining. Furthermore, our results suggest that the growth acceleration in mining areas is only temporary and diminishes with the closure of the mine.


title: Maximum-distance Race Strategies for a Fully Electric Endurance Race Car

abstract: This paper presents a bi-level optimization framework to compute the maximum-distance stint and charging strategies for a fully electric endurance race car. Thereby, the lower level computes the minimum-stint-time Powertrain Operation (PO) for a given battery energy budget and stint length, whilst the upper level leverages that information to jointly optimize the stint length, charge time and number of pit stops, in order to maximize the driven distance in the course of a fixed-time endurance race. Specifically, we first extend a convex lap time optimization framework to capture multiple laps and force-based electric motor models, and use it to create a map linking the charge time and stint length to the achievable stint time. Second, we leverage the map to frame the maximum-race-distance problem as a mixed-integer second order conic program that can be efficiently solved to the global optimum with off-the-shelf optimization algorithms. Finally, we showcase our framework on a 6 h race around the Zandvoort circuit. Our results show that a flat-out strategy can be extremely detrimental, and that, compared to when the stints are optimized for a fixed number of pit stops, jointly optimizing the stints and number of pit stops can increase the driven distance of several laps.


title: New scenario of QCD axion clump formation I: Linear analysis

abstract: The QCD axion acquires the potential through the non-perturbative effect of the QCD matters around the QCD phase transition. During this period, the direct interaction between the axion and the QCD matters sets in. Focusing on the impact of this direct interaction, we propose two scenarios where the fluctuation of the axion can rapidly grow, potentially leading to the formation of axion miniclusters even if the Peccei-Quinn (PQ) symmetry was already broken during inflation. The first scenario assumes that the primordial curvature perturbation at the horizon scale during the QCD epoch was significantly enhanced and the second one assumes that the initial misalignment was tuned around the hilltop of the potential.


title: Electroosmosis in nanopores: Computational methods and technological applications

abstract: Electroosmosis is a fascinating effect where liquid motion is induced by an applied electric field. Counter ions accumulate in the vicinity of charged surfaces, triggering a coupling between liquid mass transport and external electric field. In nanofluidic technologies, where surfaces play an exacerbated role, electroosmosis is thus of primary importance. Its consequences on the transport properties in biological and synthetic nanopores are subtle and intricate. Thorough understanding is therefore challenging yet crucial to fully assess the mechanisms at play. Here, we review recent progress on computational techniques for the analysis of electroosmosis and discuss technological applications, in particular for sensing devices.


title: Bayesian Approach to Inverse Problems: an Application to NNPDF Closure Testing

abstract: We discuss the Bayesian approach to the solution of inverse problems and apply the formalism to analyse the closure tests performed by the NNPDF collaboration. Starting from a comparison with the approach that is currently used for the determination of parton distributions (PDFs) by the NNPDF collaboration, we discuss some analytical results that can be obtained for linear problems and use these results as a guidance for the more complicated non-linear problems. We show that, in the case of Gaussian distributions, the posterior probability density of the parametrized PDFs is fully determined by the results of the NNPDF fitting procedure. In the particular case that we consider, the fitting procedure and the Bayesian analysis yield exactly the same result. Building on the insight that we obtain from the analytical results, we introduce new estimators to assess the statistical faithfulness of the fit results in closure tests. These estimators are defined in data space, and can be studied analytically using the Bayesian formalism in a linear model in order to clarify their meaning. Finally we present numerical results from a number of closure tests performed with current NNPDF methodologies. These further tests allow us to validate the NNPDF4.0 methodology and provide a quantitative comparison of the NNPDF4.0 and NNPDF3.1 methodologies. As PDFs determinations move into precision territory, the need for a careful validation of the methodology becomes increasingly important: the error bar has become the focal point of contemporary PDFs determinations. In this perspective, theoretical assumptions and other sources of error are best formulated and analysed in the Bayesian framework, which provides an ideal language to address the precision and the accuracy of current fits.


title: Science Case for the new High-Intensity Muon Beams HIMB at PSI

abstract: In April 2021, scientists active in muon physics met to discuss and work out the physics case for the new High-Intensity Muon Beams (HIMB) project at PSI that could deliver of order $10^{10}$\,s$^{-1}$ surface muons to experiments. Ideas and concrete proposals were further substantiated over the following months and assembled in the present document. The high intensities will allow for completely new experiments with considerable discovery potential and unique sensitivities. The physics case is outstanding and extremely rich, ranging from fundamental particle physics via chemistry to condensed matter research and applications in energy research and elemental analysis. In all these fields, HIMB will ensure that the facilities S$\mu$S and CHRISP on PSI's High Intensity Proton Accelerator complex HIPA remain world-leading, despite the competition of muon facilities elsewhere.


title: Evaluation of Deep Learning Topcoders Method for Neuron Individualization in Histological Macaque Brain Section

abstract: Cell individualization has a vital role in digital pathology image analysis. Deep Learning is considered as an efficient tool for instance segmentation tasks, including cell individualization. However, the precision of the Deep Learning model relies on massive unbiased dataset and manual pixel-level annotations, which is labor intensive. Moreover, most applications of Deep Learning have been developed for processing oncological data. To overcome these challenges, i) we established a pipeline to synthesize pixel-level labels with only point annotations provided; ii) we tested an ensemble Deep Learning algorithm to perform cell individualization on neurological data. Results suggest that the proposed method successfully segments neuronal cells in both object-level and pixel-level, with an average detection accuracy of 0.93.


title: Early Myocardial Infarction Detection over Multi-view Echocardiography

abstract: Myocardial infarction (MI) is the leading cause of mortality in the world that occurs due to a blockage of the coronary arteries feeding the myocardium. An early diagnosis of MI and its localization can mitigate the extent of myocardial damage by facilitating early therapeutic interventions. Following the blockage of a coronary artery, the regional wall motion abnormality (RWMA) of the ischemic myocardial segments is the earliest change to set in. Echocardiography is the fundamental tool to assess any RWMA. Assessing the motion of the left ventricle (LV) wall only from a single echocardiography view may lead to missing the diagnosis of MI as the RWMA may not be visible on that specific view. Therefore, in this study, we propose to fuse apical 4-chamber (A4C) and apical 2-chamber (A2C) views in which a total of 12 myocardial segments can be analyzed for MI detection. The proposed method first estimates the motion of the LV wall by Active Polynomials (APs), which extract and track the endocardial boundary to compute myocardial segment displacements. The features are extracted from the A4C and A2C view displacements, which are concatenated and fed into the classifiers to detect MI. The main contributions of this study are 1) creation of a new benchmark dataset by including both A4C and A2C views in a total of 260 echocardiography recordings, which is publicly shared with the research community, 2) improving the performance of the prior work of threshold-based APs by a Machine Learning based approach, and 3) a pioneer MI detection approach via multi-view echocardiography by fusing the information of A4C and A2C views. Experimental results show that the proposed method achieves 90.91% sensitivity and 86.36% precision for MI detection over multi-view echocardiography. The software implementation is shared at https://github.com/degerliaysen/MultiEchoAI.

categories: eess.IV cs.AI cs.CV cs.LG physics.med-ph


title: Distribution-Invariant Differential Privacy

abstract: Differential privacy is becoming one gold standard for protecting the privacy of publicly shared data. It has been widely used in social science, data science, public health, information technology, and the U.S. decennial census. Nevertheless, to guarantee differential privacy, existing methods may unavoidably alter the conclusion of the original data analysis, as privatization often changes the sample distribution. This phenomenon is known as the trade-off between privacy protection and statistical accuracy. In this work, we mitigate this trade-off by developing a distribution-invariant privatization (DIP) method to reconcile both high statistical accuracy and strict differential privacy. As a result, any downstream statistical or machine learning task yields essentially the same conclusion as if one used the original data. Numerically, under the same strictness of privacy protection, DIP achieves superior statistical accuracy in a wide range of simulation studies and real-world benchmarks.


title: HARPO: Learning to Subvert Online Behavioral Advertising

abstract: Online behavioral advertising, and the associated tracking paraphernalia, poses a real privacy threat. Unfortunately, existing privacy-enhancing tools are not always effective against online advertising and tracking. We propose Harpo, a principled learning-based approach to subvert online behavioral advertising through obfuscation. Harpo uses reinforcement learning to adaptively interleave real page visits with fake pages to distort a tracker's view of a user's browsing profile. We evaluate Harpo against real-world user profiling and ad targeting models used for online behavioral advertising. The results show that Harpo improves privacy by triggering more than 40% incorrect interest segments and 6x higher bid values. Harpo outperforms existing obfuscation tools by as much as 16x for the same overhead. Harpo is also able to achieve better stealthiness to adversarial detection than existing obfuscation tools. Harpo meaningfully advances the state-of-the-art in leveraging obfuscation to subvert online behavioral advertising


title: Sterile Neutrino Searches with MicroBooNE: Electron Neutrino Disappearance

abstract: A sterile neutrino is a well motivated minimal new physics model that leaves an imprint in neutrino oscillations. Over the last two decades, a number of hints pointing to a sterile neutrino have emerged, many of which are pointing near $m_4\sim1$ eV. Here we show how MicroBooNE data can be used to search for electron neutrino disappearance using each of their four analysis channels. We find a hint for oscillations with the highest single channel significance of $2.4\sigma$ (using the Feldman-Cousins approach) coming from the Wire-Cell analysis and a simplified treatment of the experimental systematics. The preferred parameters are $\sin^2(2\theta_{14})=0.35^{+0.19}_{-0.16}$ and $\Delta m^2_{41}=1.25^{+0.74}_{-0.39}$ eV$^2$. This region of parameter space is in good agreement with existing hints from source experiments, is at a similar frequency but higher mixing than indicated by reactor anti-neutrinos, and is at the edge of the region allowed by solar neutrino data. Existing unanalyzed data from MicroBooNE could increase the sensitivity to the $>3\sigma$ level.


title: PIMIP: An Open Source Platform for Pathology Information Management and Integration

abstract: Digital pathology plays a crucial role in the development of artificial intelligence in the medical field. The digital pathology platform can make the pathological resources digital and networked, and realize the permanent storage of visual data and the synchronous browsing processing without the limitation of time and space. It has been widely used in various fields of pathology. However, there is still a lack of an open and universal digital pathology platform to assist doctors in the management and analysis of digital pathological sections, as well as the management and structured description of relevant patient information. Most platforms cannot integrate image viewing, annotation and analysis, and text information management. To solve the above problems, we propose a comprehensive and extensible platform PIMIP. Our PIMIP has developed the image annotation functions based on the visualization of digital pathological sections. Our annotation functions support multi-user collaborative annotation and multi-device annotation, and realize the automation of some annotation tasks. In the annotation task, we invited a professional pathologist for guidance. We introduce a machine learning module for image analysis. The data we collected included public data from local hospitals and clinical examples. Our platform is more clinical and suitable for clinical use. In addition to image data, we also structured the management and display of text information. So our platform is comprehensive. The platform framework is built in a modular way to support users to add machine learning modules independently, which makes our platform extensible.


title: The Weingarten map and curvatures of SL(n, R)

abstract: In this article, we shall derive by elementary calculations the Gauss map, spherical image, Weingarten map and the curvatures at identity of the special linear group, that is, the matrices of determinant 1. We could not find any reference where this simple, nice computation has been explicitly put down.


title: Human-Centric Decision Support Tools: Insights from Real-World Design and Implementation

abstract: Decision support tools enable improved decision-making for challenging decision problems by empowering stakeholders to process, analyze, visualize, and otherwise make sense of a variety of key factors. Their intentional design is a critical component of the value they create. All decision-support tools share in common that there is a complex decision problem to be solved for which decision-support is useful, and moreover, that appropriate analytics expertise is available to produce solutions to the problem setting at hand. When well-designed, decision support tools reduce friction and increase efficiency in providing support for the decision-making process, thereby improving the ability of decision-makers to make quality decisions. On the other hand, the presence of overwhelming, superfluous, insufficient, or ill-fitting information and software features can have an adverse effect on the decision-making process and, consequently, outcomes. We advocate for an innovative, and perhaps overlooked, approach to designing effective decision support tools: genuinely listening to the project stakeholders, to ascertain and appreciate their real needs and perspectives. By prioritizing stakeholder needs, a foundation of mutual trust and understanding is established with the design team. We maintain this trust is critical to eventual tool acceptance and adoption, and its absence jeopardizes the future use of the tool, which would leave its analytical insights for naught. We discuss examples across multiple contexts to underscore our collective experience, highlight lessons learned, and present recommended practices to improve the design and eventual adoption of decision dupport tools.


title: A systematic study of the valence electronic structure of cyclo(Gly-Phe), cyclo(Trp-Tyr) and cyclo(Trp-Trp) dipeptides in gas phase

abstract: The electronic energy levels of cyclo(Glycine-Phenylalanine), cyclo(Tryptophan-Tyrosine) and cyclo(Tryptophan-Tryptophan) dipeptides are investigated with a joint experimental and theoretical approach. Experimentally, valence photoelectron spectra in the gas phase are measured using VUV radiation. Theoretically, we first obtain low-energy conformers through an automated conformer-rotamer ensemble sampling scheme based on tight-binding simulations. Then, different first principles computational schemes are considered to simulate the spectra: Hartree-Fock (HF), density functional theory (DFT) within the B3LYP approximation, the quasi--particle GW correction, and the quantum-chemistry CCSD method. Theory allows to assign the main features of the spectra. A discussion on the role of electronic correlation is provided, by comparing computationally cheaper DFT scheme (and GW) results with the accurate CCSD method.


title: On the evaluation of the Appell $F_2$ double hypergeometric function

abstract: The transformation theory of the Appell $F_2(a,b_1,b_2;c_1,c_2;x,y)$ double hypergeometric function is used to obtain a set of series representations of $F_2$ which provide an efficient way to evaluate $F_2$ for real values of its arguments $x$ and $y$ and generic complex values of its parameters $a,b_1, b_2, c_1$ and $c_2$ (i.e. in the nonlogarithmic case). This study rests on a classical approach where the usual double series representation of $F_2$ and other double hypergeometric series that appear in the intermediate steps of the calculations are written as infinite sums of one variable hypergeometric series, such as the Gauss $_2F_1$ or the $_3F_2$, various linear transformations of the latter being then applied to derive known and new formulas. Using the three well-known Euler transformations of $F_2$ on these results allows us to obtain a total of 44 series which form the basis of the Mathematica package AppellF2, dedicated to the evaluation of $F_2$. A brief description of the package and of the numerical analysis that we have performed to test it are also presented.


title: Computing the Dimension of a Bipartition Matrix

abstract: This article presents a computer program that computes the dimension of a bipartition matrix. Its dimension has three independent components: row dimension, column dimension, and entry dimension. The program applies four routines of independent interest: (1) a routine that factors a bipartition as a formal product of indecomposables, (2) a routine that computes the inverse and recovers the bipartition, (3) a routine that factors a bipartition matrix as a formal product of indecomposables, and (4) a routine that calculates the "transpose-rotation" of a bipartition matrix (the column dimension is the row dimension of the transpose-rotation).


title: Arbitrary order principal directions and how to compute them

abstract: Curvature principal directions on geometric surfaces are a ubiquitous concept of Geometry Processing techniques. However they only account for order 2 differential quantities, oblivious of higher order differential behaviors. In this paper, we extend the concept of principal directions to higher orders for surfaces in R^3 by considering symmetric differential tensors. We further show how they can be explicitly approximated on point set surfaces and that they convey valuable geometric information, that can help the analysis of 3D surfaces.


title: The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

abstract: We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound.


title: A Stochastic Graph-based Model for the Simulation of SARS-CoV-2 Transmission

abstract: In this work we propose the design principles of a stochastic graph-based model for the simulation of SARS-CoV-2 transmission. The proposed approach incorporates three sub-models, namely, the spatial model, the mobility model, and the propagation model, in order to develop a realistic environment for the study of the properties exhibited by the spread of SARS-CoV-2. The spatial model converts images of real cities taken from Google Maps into undirected weighted graphs that capture the spatial arrangement of the streets utilized next for the mobility of individuals. The mobility model implements a stochastic agent-based approach, developed in order to assign specific routes to individuals moving in the city, through the use of stochastic processes, utilizing the weights of the underlying graph to deploy shortest path algorithms. The propagation model implements both the epidemiological model and the physical substance of the transmission of an airborne virus considering the transmission parameters of SARS-CoV-2. Finally, we integrate these sub-models in order to derive an integrated framework for the study of the epidemic dynamics exhibited through the transmission of SARS-CoV-2.


title: Gradients are Not All You Need

abstract: Differentiable programming techniques are widely used in the community and are responsible for the machine learning renaissance of the past several decades. While these methods are powerful, they have limits. In this short report, we discuss a common chaos based failure mode which appears in a variety of differentiable circumstances, ranging from recurrent neural networks and numerical physics simulation to training learned optimizers. We trace this failure to the spectrum of the Jacobian of the system under study, and provide criteria for when a practitioner might expect this failure to spoil their differentiation based optimization algorithms.


title: Secure and Reliable Transfer Learning Framework for 6G-enabled Internet of Vehicles

abstract: In the coming 6G era, Internet of Vehicles (IoV) has been evolving towards 6G-enabled IoV with super-high data rate, seamless networking coverage, and ubiquitous intelligence by Artificial Intelligence (AI). Transfer Learning (TL) has great potential to empower promising 6G-enabled IoV, such as smart driving assistance, with its outstanding features including enhancing quality and quantity of training data, speeding up learning processes and reducing computing demands. Although TL had been widely adopted in wireless applications (e.g., spectrum management and caching), its reliability and security in 6G-enabled IoV were still not well investigated. For instance, malicious vehicles in source domains may transfer and share untrustworthy models (i.e., knowledge) about connection availability to target domains, thus adversely affecting the performance of learning processes. Therefore, it is important to select and also incentivize trustworthy vehicles to participate in TL. In this article, we first introduce the integration of TL and 6G-enbaled IoV and provide TL applications for 6G-enabled IoV. We then design a secure and reliable transfer learning framework by using reputation to evaluate the reliability of pre-trained models and utilizing the consortium blockchain to achieve secure and efficient decentralized reputation management. Moreover, a deep learning-based auction scheme for the TL model market is designed to motivate high-reputation vehicles to participate in model sharing. Finally, the simulation results demonstrate that the proposed framework is secure and reliable with well-design incentives for TL in 6G-enabled IoV.


title: Cross-lingual Adaption Model-Agnostic Meta-Learning for Natural Language Understanding

abstract: Meta learning with auxiliary languages has demonstrated promising improvements for cross-lingual natural language processing. However, previous studies sample the meta-training and meta-testing data from the same language, which limits the ability of the model for cross-lingual transfer. In this paper, we propose XLA-MAML, which performs direct cross-lingual adaption in the meta-learning stage. We conduct zero-shot and few-shot experiments on Natural Language Inference and Question Answering. The experimental results demonstrate the effectiveness of our method across different languages, tasks, and pretrained models. We also give analysis on various cross-lingual specific settings for meta-learning including sampling strategy and parallelism.


title: On the efficiency of a general attack against the MOBS cryptosystem

abstract: All instances of the semidirect key exchange protocol, a generalisation of the famous Diffie-Hellman key exchange protocol, satisfy the so-called "telescoping equality"; in some cases, this equality has been used to construct an attack. In this report we present computational evidence suggesting that an instance of the scheme called `MOBS' is an example of a scheme where the telescoping equality has too many solutions to be a practically viable means to conduct an attack.


title: On the functional CLT for slowly mixing triangular arrays

abstract: In \cite{MPU} a functional CLT was obtained for triangular arrays satisfying the Lindeberg condition, that the sum of the individual variances is at most the same order as the variance of the underlying sum, and under the optimal mixing rats $\sum_{n}\rho(2^n)<\infty$, where $\rho(\cdot)$ are the $\rho$-mixing coefficients of the array. In this paper we will present alternative conditions which do not involve the assumption on the sum of variances, and instead we will assume certain maximal moment assumptions (which we can verify for $\phi$-mixing arrays) and mixing rates of the form $\sum_n\rho(e^{G(n)})<\infty$ where $G(n)$ grows sub-linearly fast in $n$ (e.g. $G(n)=n/\ln(\ln n)$). We will also discuss alternative conditions to the ones in the functional CLT for $\alpha$-mixing triangular arrays which was obtained in \cite{MP}.


title: BagBERT: BERT-based bagging-stacking for multi-topic classification

abstract: This paper describes our submission on the COVID-19 literature annotation task at Biocreative VII. We proposed an approach that exploits the knowledge of the globally non-optimal weights, usually rejected, to build a rich representation of each label. Our proposed approach consists of two stages: (1) A bagging of various initializations of the training data that features weakly trained weights, (2) A stacking of heterogeneous vocabulary models based on BERT and RoBERTa Embeddings. The aggregation of these weak insights performs better than a classical globally efficient model. The purpose is the distillation of the richness of knowledge to a simpler and lighter model. Our system obtains an Instance-based F1 of 92.96 and a Label-based micro-F1 of 91.35.


title: Observationally Constraining the Starspot Properties of Magnetically Active M67 Sub-Subgiant S1063

abstract: Our understanding of the impact of magnetic activity on stellar evolution continues to unfold. This impact is seen in sub-subgiant stars, defined to be stars that sit below the subgiant branch and red of the main sequence in a cluster color-magnitude diagram. Here we focus on S1063, a prototypical sub-subgiant in open cluster M67. We use a novel technique combining a two-temperature spectral decomposition and light curve analysis to constrain starspot properties over a multi-year time frame. Using a high-resolution near-infrared IGRINS spectrum and photometric data from K2 and ASAS-SN, we find a projected spot filling factor of 32 $\pm$ 7% with a spot temperature of 4000 $\pm$ 200 K. This value anchors the variability seen in the light curve, indicating the spot filling factor of S1063 ranged from 20% to 45% over a four-year time period with an average spot filling factor of 30%. These values are generally lower than those determined from photometric model comparisons but still indicate that S1063, and likely other sub-subgiants, are magnetically active spotted stars. We find observational and theoretical comparisons of spotted stars are nuanced due to the projected spot coverage impacting estimates of the surface-averaged effective temperature. The starspot properties found here are similar to those found in RS CVn systems, supporting classifying sub-subgiants as another type of active giant star binary system. This technique opens the possibility of characterizing the surface conditions of many more spotted stars than previous methods, allowing for larger future studies to test theoretical models of magnetically active stars.


title: Homological properties of pinched Veronese rings

abstract: Pinched Veronese rings are formed by removing an algebra generator from a Veronese subring of a polynomial ring. We study the homological properties of such rings, including the Cohen-Macaulay, Gorenstein, and complete intersection properties. Greco and Martino classified Cohen-Macaulayness of pinched Veronese rings by the maximum entry of the exponent vector of the pinched monomial; we re-prove their results with semigroup methods and correct an omission of a small class of examples of Cohen-Macaulay pinched Veronese rings. When the underlying field is of prime characteristic, we show that pinched Veronese rings exhibit a variety of F-singularities, including F-regular, F-injective, and F-nilpotent. We also compute upper bounds on the Frobenius test exponents of pinched Veronese rings, a computational invariant which controls the Frobenius closure of all parameter ideals simultaneously.


title: Internet of Things (IoT) Connectivity in 6G: An Interplay of Time, Space, Intelligence, and Value

abstract: Internet of Things (IoT) connectivity has a prominent presence in the 5G wireless communication systems. As these systems are being deployed, there is a surge of research efforts and visions towards 6G wireless systems. In order to position the evolution of IoT within the 6G systems, this paper first takes a critical view on the way IoT connectivity is supported within 5G. Following that, the wireless IoT evolution is discussed through multiple dimensions: time, space, intelligence, and value. We also conjecture that the focus will broaden from IoT devices and their connections towards the emergence of complex IoT environments, seen as building blocks of the overall IoT ecosystem.


title: A study of sunspot 3 minute oscillations using ALMA and GST

abstract: Waves and oscillations are important solar phenomena, not only because they can propagate and dissipate energy in the chromosphere, but also because they carry information about the structure of the atmosphere in which they propagate. The nature of the three-minute oscillations observed in the umbral region of sunspots is considered to be an effect of propagation of magnetohydrodynamic (MHD) waves upward from below the photosphere. We present a study of sunspot oscillations and wave propagation in NOAA AR 12470 using an approximately one-hour long data set acquired on 2015 December 17 by the Atacama Large Millimeter/submillimeter Array (ALMA), the Goode Solar Telescope (GST) operating at the Big Bear Solar Observatory (BBSO), the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), and the Interface Region Imaging Spectrograph (IRIS). The ALMA data are unique in providing a time-series of direct temperature measurements in the sunspot chromosphere. The two-second cadence of ALMA images allows us to well resolve the three-minute periods typical of sunspot oscillations in the chromosphere. Fourier analysis is applied to ALMA Band 3 ($\sim$100 GHz, $\sim$3 mm) and GST H$\alpha$ data sets to obtain power spectra as well as oscillation phase information. We analysed properties of the wave propagation by combining multiple wavelengths that probe physical parameters of solar atmosphere at different heights. We find that the ALMA temperature fluctuations are consistent with that expected for a propagating acoustic wave, with a slight asymmetry indicating non-linear steepening.


title: Congruences of elliptic curves arising from non-surjective mod $N$ Galois representations

abstract: We study $N$-congruences between quadratic twists of elliptic curves. If $N$ has exactly two distinct prime factors we show that these are parametrised by double covers of certain modular curves. In many, but not all cases, the modular curves in question correspond to the normaliser of a Cartan subgroup of $\mathrm{GL}_2(\mathbb{Z}/N\mathbb{Z})$. By computing explicit models for these double covers we find all pairs $(N, r)$ such that there exist infinitely many $j$-invariants of elliptic curves $E/\mathbb{Q}$ which are $N$-congruent with power $r$ to a quadratic twist of $E$. We also find an example of a $48$-congruence over $\mathbb{Q}$. We make a conjecture classifying nontrivial $(N,r)$-congruences between quadratic twists of elliptic curves over $\mathbb{Q}$. Finally, we give a more detailed analysis of the level $15$ case. We use elliptic Chabauty to determine the rational points on a modular curve of genus $2$ whose Jacobian has rank $2$ and which arises as a double cover of the modular curve $X(\mathrm{ns} 3^+, \mathrm{ns} 5^+)$. As a consequence we obtain a new proof of the class number $1$ problem.


title: SwAMP: Swapped Assignment of Multi-Modal Pairs for Cross-Modal Retrieval

abstract: We tackle the cross-modal retrieval problem, where learning is only supervised by relevant multi-modal pairs in the data. Although the contrastive learning is the most popular approach for this task, it makes potentially wrong assumption that the instances in different pairs are automatically irrelevant. To address the issue, we propose a novel loss function that is based on self-labeling of the unknown semantic classes. Specifically, we aim to predict class labels of the data instances in each modality, and assign those labels to the corresponding instances in the other modality (i.e., swapping the pseudo labels). With these swapped labels, we learn the data embedding for each modality using the supervised cross-entropy loss. This way, cross-modal instances from different pairs that are semantically related can be aligned to each other by the class predictor. We tested our approach on several real-world cross-modal retrieval problems, including text-based video retrieval, sketch-based image retrieval, and image-text retrieval. For all these tasks our method achieves significant performance improvement over the contrastive learning.


title: (Mis)perceptions and Engagement on Twitter: COVID-19 Vaccine Rumors on Efficacy and Mass Immunization Effort

abstract: This paper reports the findings of a 606-participant study where we analyzed the perception and engagement effects of COVID-19 vaccine rumours on Twitter pertaining to (a) vaccine efficacy; and (b) mass immunization efforts in the United States. Misperceptions regarding vaccine efficacy were successfully induced through simple content alterations and the addition of popular anti COVID-19 hashtags to otherwise valid Twitter content. Twitter's misinformation contextual tags caused a "backfire effect" for the skeptic, vaccine-hesitant reinforcing their opposition stance. While the majority of the participants staunchly refrain from engaging with the COVID-19 rumours, the skeptic, vaccine-hesitant ones were open to comment, re-tweet, like and share the vaccine efficacy rumors. We discuss the implications of our results in the context of broadening the effort for dispelling rumors about COVID-19 on social media.


title: Geometric Bounds on the Fastest Mixing Markov Chain

abstract: In the Fastest Mixing Markov Chain problem, we are given a graph $G = (V, E)$ and desire the discrete-time Markov chain with smallest mixing time $\tau$ subject to having equilibrium distribution uniform on $V$ and non-zero transition probabilities only across edges of the graph. It is well-known that the mixing time $\tau_\textsf{RW}$ of the lazy random walk on $G$ is characterised by the edge conductance $\Phi$ of $G$ via Cheeger's inequality: $\Phi^{-1} \lesssim \tau_\textsf{RW} \lesssim \Phi^{-2} \log |V|$. Analogously, we characterise the fastest mixing time $\tau^\star$ via a Cheeger-type inequality but for a different geometric quantity, namely the vertex conductance $\Psi$ of $G$: $\Psi^{-1} \lesssim \tau^\star \lesssim \Psi^{-2} (\log |V|)^2$. This characterisation forbids fast mixing for graphs with small vertex conductance. To bypass this fundamental barrier, we consider Markov chains on $G$ with equilibrium distribution which need not be uniform, but rather only $\varepsilon$-close to uniform in total variation. We show that it is always possible to construct such a chain with mixing time $\tau \lesssim \varepsilon^{-1} (\operatorname{diam} G)^2 \log |V|$. Finally, we discuss analogous questions for continuous-time and time-inhomogeneous chains.


title: Quaternary quartic forms and Gorenstein rings

abstract: A quaternary quartic form, a quartic form in four variables, is the dual socle generator of an Artinian Gorenstein ring of codimension and regularity 4. We present a classification of quartic forms in terms of rank and powersum decompositions which corresponds to the classification by the Betti tables of the corresponding Artinian Gorenstein rings. This gives a stratification of the space of quaternary quartic forms which we compare with the Noether-Lefschetz stratification. We discuss various phenomena related to this stratification. We study the geometry of powersum varieties for a general form in each stratum. In particular, we show that the powersum variety $VSP(F,9)$ of a general quartic with singular middle catalecticant is again a quartic surface, thus giving a rational map between two divisors in the space of quartics. Finally, we provide various explicit constructions of general Artinian Gorenstein rings corresponding to each stratum and discuss their lifting to higher dimension. These provide constructions of codimension four varieties, which include canonical surfaces, Calabi-Yau threefolds and Fano fourfolds. In the particular case of quaternary quartics, our results yield answers to questions posed by Geramita, Iarrobino-Kanev, and Reid.


title: Efficient Projection-Free Online Convex Optimization with Membership Oracle

abstract: In constrained convex optimization, existing methods based on the ellipsoid or cutting plane method do not scale well with the dimension of the ambient space. Alternative approaches such as Projected Gradient Descent only provide a computational benefit for simple convex sets such as Euclidean balls, where Euclidean projections can be performed efficiently. For other sets, the cost of the projections can be too high. To circumvent these issues, alternative methods based on the famous Frank-Wolfe algorithm have been studied and used. Such methods use a Linear Optimization Oracle at each iteration instead of Euclidean projections; the former can often be performed efficiently. Such methods have also been extended to the online and stochastic optimization settings. However, the Frank-Wolfe algorithm and its variants do not achieve the optimal performance, in terms of regret or rate, for general convex sets. What is more, the Linear Optimization Oracle they use can still be computationally expensive in some cases. In this paper, we move away from Frank-Wolfe style algorithms and present a new reduction that turns any algorithm A defined on a Euclidean ball (where projections are cheap) to an algorithm on a constrained set C contained within the ball, without sacrificing the performance of the original algorithm A by much. Our reduction requires O(T log T) calls to a Membership Oracle on C after T rounds, and no linear optimization on C is needed. Using our reduction, we recover optimal regret bounds [resp. rates], in terms of the number of iterations, in online [resp. stochastic] convex optimization. Our guarantees are also useful in the offline convex optimization setting when the dimension of the ambient space is large.


title: Look Before You Leap: Safe Model-Based Reinforcement Learning with Human Intervention

abstract: Safety has become one of the main challenges of applying deep reinforcement learning to real world systems. Currently, the incorporation of external knowledge such as human oversight is the only means to prevent the agent from visiting the catastrophic state. In this paper, we propose MBHI, a novel framework for safe model-based reinforcement learning, which ensures safety in the state-level and can effectively avoid both "local" and "non-local" catastrophes. An ensemble of supervised learners are trained in MBHI to imitate human blocking decisions. Similar to human decision-making process, MBHI will roll out an imagined trajectory in the dynamics model before executing actions to the environment, and estimate its safety. When the imagination encounters a catastrophe, MBHI will block the current action and use an efficient MPC method to output a safety policy. We evaluate our method on several safety tasks, and the results show that MBHI achieved better performance in terms of sample efficiency and number of catastrophes compared to the baselines.


title: Multi-Task Neural Processes

abstract: Neural processes have recently emerged as a class of powerful neural latent variable models that combine the strengths of neural networks and stochastic processes. As they can encode contextual data in the network's function space, they offer a new way to model task relatedness in multi-task learning. To study its potential, we develop multi-task neural processes, a new variant of neural processes for multi-task learning. In particular, we propose to explore transferable knowledge from related tasks in the function space to provide inductive bias for improving each individual task. To do so, we derive the function priors in a hierarchical Bayesian inference framework, which enables each task to incorporate the shared knowledge provided by related tasks into its context of the prediction function. Our multi-task neural processes methodologically expand the scope of vanilla neural processes and provide a new way of exploring task relatedness in function spaces for multi-task learning. The proposed multi-task neural processes are capable of learning multiple tasks with limited labeled data and in the presence of domain shift. We perform extensive experimental evaluations on several benchmarks for the multi-task regression and classification tasks. The results demonstrate the effectiveness of multi-task neural processes in transferring useful knowledge among tasks for multi-task learning and superior performance in multi-task classification and brain image segmentation.


title: Collision-Free Continuum Deformation Coordination of a Multi-Quadcopter System Using Cooperative Localization

abstract: This paper integrates cooperative localization with continuum deformation coordination of a multi-quadcopter system {(MQS) to assure safety and optimality of the quadcopter team coordination in the presence of position uncertainty. We first consider the MQS as a finite number of particles of a deformable triangle in a $3$-D motion space and define their continuum deformation coordination as a leader-follower problem in which leader quadcopters} can estimate (know) their positions but follower quadcopters rely on relative position measurements to localize themselves and estimate the leaders' positions. We then propose a navigation strategy for the MQS to plan and acquire the desired continuum deformation coordination, in the presence of measurement noise, disturbance, and position uncertainties, such that collision is avoided and rotor angular speeds of all quadcopters remain bounded. We show the efficacy of the proposed strategy by simulating the continuum deformation coordination of an MQS with eight quadcopters.


title: $L^p$-theory for the exterior Stokes problem with Navier's type slip-without-friction boundary conditions

abstract: In this paper, we consider the stationary Stokes equations in an exterior domain three-dimensional under a slip boundary condition without friction. We set the problem in weighted Sobolev spaces in order to control the behavior at infinity of the solutions. In this work, we try to investigate the existence and uniqueness of the weak solutions related to this problem in $L^P-$theory when $p>3$. Our proof are based on obtaining inf-sup conditions that play a fundamental role.


title: Understanding COVID-19 Vaccine Reaction through Comparative Analysis on Twitter

abstract: Although multiple COVID-19 vaccines have been available for several months now, vaccine hesitancy continues to be at high levels in the United States. In part, the issue has also become politicized, especially since the presidential election in November. Understanding vaccine hesitancy during this period in the context of social media, including Twitter, can provide valuable guidance both to computational social scientists and policy makers. Rather than studying a single Twitter corpus, this paper takes a novel view of the problem by comparatively studying two Twitter datasets collected between two different time periods (one before the election, and the other, a few months after) using the same, carefully controlled data collection and filtering methodology. Our results show that there was a significant shift in discussion from politics to COVID-19 vaccines from fall of 2020 to spring of 2021. By using clustering and machine learning-based methods in conjunction with sampling and qualitative analysis, we uncover several fine-grained reasons for vaccine hesitancy, some of which have become more (or less) important over time. Our results also underscore the intense polarization and politicization of this issue over the last year.


title: Very weak solution for the exterior stationary Stokes equations with Navier slip boundary condition

abstract: In some problems of fluid mechanics, it is possible to be confronted with data that are not regular, that is why we are interested here in the search for the so-called very weak solutions for the stationary Stokes problem with Navier-type boundary conditions in a three-dimensional exterior domain. The problem describes the flow of a viscous and incompressible fluid past an obstacle where we assume that the fluid may slip on the boundary of the obstacle. Because the flow domain is unbounded, we set the problem in weighted Sobolev spaces in order to control the behavior at infinity of the solutions. Our purpose is to prove the existence and the uniqueness of a very weak solution in a Hilbertian framework.


title: A Two-Stage Approach towards Generalization in Knowledge Base Question Answering

abstract: Most existing approaches for Knowledge Base Question Answering (KBQA) focus on a specific underlying knowledge base either because of inherent assumptions in the approach, or because evaluating it on a different knowledge base requires non-trivial changes. However, many popular knowledge bases share similarities in their underlying schemas that can be leveraged to facilitate generalization across knowledge bases. To achieve this generalization, we introduce a KBQA framework based on a 2-stage architecture that explicitly separates semantic parsing from the knowledge base interaction, facilitating transfer learning across datasets and knowledge graphs. We show that pretraining on datasets with a different underlying knowledge base can nevertheless provide significant performance gains and reduce sample complexity. Our approach achieves comparable or state-of-the-art performance for LC-QuAD (DBpedia), WebQSP (Freebase), SimpleQuestions (Wikidata) and MetaQA (Wikimovies-KG).


title: Palette: Image-to-Image Diffusion Models

abstract: This paper develops a unified framework for image-to-image translation based on conditional diffusion models and evaluates this framework on four challenging image-to-image translation tasks, namely colorization, inpainting, uncropping, and JPEG restoration. Our simple implementation of image-to-image diffusion models outperforms strong GAN and regression baselines on all tasks, without task-specific hyper-parameter tuning, architecture customization, or any auxiliary loss or sophisticated new techniques needed. We uncover the impact of an L2 vs. L1 loss in the denoising diffusion objective on sample diversity, and demonstrate the importance of self-attention in the neural architecture through empirical studies. Importantly, we advocate a unified evaluation protocol based on ImageNet, with human evaluation and sample quality scores (FID, Inception Score, Classification Accuracy of a pre-trained ResNet-50, and Perceptual Distance against original images). We expect this standardized evaluation protocol to play a role in advancing image-to-image translation research. Finally, we show that a generalist, multi-task diffusion model performs as well or better than task-specific specialist counterparts. Check out https://diffusion-palette.github.io for an overview of the results.


title: Data-Driven AI Model Signal-Awareness Enhancement and Introspection

abstract: AI modeling for source code understanding tasks has been making significant progress, and is being adopted in production development pipelines. However, reliability concerns, especially whether the models are actually learning task-related aspects of source code, are being raised. While recent model-probing approaches have observed a lack of signal awareness in many AI-for-code models, i.e. models not capturing task-relevant signals, they do not offer solutions to rectify this problem. In this paper, we explore data-driven approaches to enhance models' signal-awareness: 1) we combine the SE concept of code complexity with the AI technique of curriculum learning; 2) we incorporate SE assistance into AI models by customizing Delta Debugging to generate simplified signal-preserving programs, augmenting them to the training dataset. With our techniques, we achieve up to 4.8x improvement in model signal awareness. Using the notion of code complexity, we further present a novel model learning introspection approach from the perspective of the dataset.


title: Existence and decay of traveling waves for the nonlocal Gross-Pitaevskii equation

abstract: We consider the nonlocal Gross-Pitaevskii equation that models a Bose gas with general nonlocal interactions between particles in one spatial dimension, with constant density far away. We address the problem of the existence of traveling waves with nonvanishing conditions at infinity, i.e. dark solitons. Under general conditions on the interactions, we prove existence of dark solitons for almost every subsonic speed. Moreover, we show existence in the whole subsonic regime for a family of potentials. The proofs rely on a Mountain Pass argument combined with the so-called "monotonicity trick", as well as on a priori estimates for the Palais-Smale sequences. Finally, we establish properties of the solitons such us exponential decay at infinity and analyticity.


title: Two-tone Doppler cooling of radial two-dimensional crystals in a radiofrequency ion trap

abstract: We study the Doppler-cooling of radial two-dimensional (2D) Coulomb crystals of trapped barium ions in a radiofrequency trap. Ions in radial 2D crystals experience micromotion of an amplitude that increases linearly with the distance from the trap center, leading to a position-dependent frequency modulation of laser light in each ion's rest frame. We use two tones of Doppler-cooling laser light separated by approximately 100~MHz to efficiently cool distinct regions in the crystals with differing amplitudes of micromotion. This technique allows us to trap and cool more than 50 ions populating 4 shells in a radial two-dimensional crystal, where with a single tone of Doppler cooling light we are limited to 30 ions in 3 shells. We also individually characterize the micromotion of all ions within the crystals, and use this information to locate the center of the trap and to determine the Matthieu parameters $q_{x}$ and $q_{y}$.


title: Well-mixed Lotka-Volterra model with random strongly competitive interactions

abstract: The random Lotka-Volterra model is widely used to describe the dynamical and thermodynamic features of ecological communities. In this work, we consider random symmetric interactions between species and analyze the strongly competitive interaction case. We investigate different scalings for the distribution of the interactions with the number of species and try to bridge the gap with previous works. Our results show two different behaviors for the mean abundance at zero and finite temperature respectively, with a continuous crossover between the two. We confirm and extend previous results obtained for weak interactions: at zero temperature, even in the strong competitive interaction limit, the system is in a multiple-equilibria phase, whereas at finite temperature only a unique stable equilibrium can exist. Finally, we establish the qualitative phase diagrams in both cases and compare the two species abundance distributions.


title: On recovering quadratic pencils with singular coefficients and entire functions in the boundary conditions

abstract: In the paper, we study an inverse spectral problem for quadratic pencils of the Sturm--Liouville operators with singular coefficients and entire functions in the boundary conditions. We prove that a subspectrum is sufficient for recovering the pencil if this subspectrum generates a complete functional system. As well, we obtain an algorithm solving the inverse problem and alternative conditions on the subspectrum. Finally, these results are applied to studying a partial inverse problem.


title: Twisted Nakano-Positivity of Fields of Hilbert Spaces

abstract: In a previous article (arXiv:2111.03143), we generalized Berndtsson's Nakano-positivity by retaining the same consequences under weaker hypotheses. In this article, we propose to further generalize our "twisted" Nakano-positivity theorem to fields of Hilbert spaces associated to possibly unbounded Stein manifolds. We achieve this result using exhaustion arguments. As direct applications, we prove $\log$-plurisubharmonic variation results for a certain class of non-trivial families of Stein manifolds.


title: C$^{18}$O emission as an effective measure of gas masses of protoplanetary disks

abstract: Many astrochemical models of observed CO isotopologue line emission, earlier considered a good proxy measure of H$_2$ and hence disk gas mass, favor large deviations in the carbon and oxygen gas phase abundances and argue that severe gas phase CO depletion makes it a poor mass tracer. Here, we show that C$^{18}$O line emission is an effective measure of the gas mass, and despite its complex chemistry, a possibly better tracer than HD. Our models are able to reproduce C$^{18}$O emission from recent ALMA surveys and the TW Hya disk to within a factor of $\sim 2-3$ using carbon and oxygen abundances characteristic of the interstellar medium (C/H$=1.4 \times 10^{-4}$; O/H$=3.2\times 10^{-4}$) without having to invoke unusual chemical processing. Our gas and dust disk structure calculations considering hydrostatic pressure equilibrium and our treatment of the CO conversion on grains are primarily responsible for the very different conclusions on disk masses and CO depletion. As did previous studies, we find that a gas phase C/O of $\sim 1-2$ can explain observed hydrocarbon emission from the TW Hya disk; but significantly, we find that CO isotopologue emission is only marginally affected by the C/O ratio. We therefore conclude that C$^{18}$O emission provides estimates of disk masses that are uncertain only to within a factor of a few, and describe a simplified modeling procedure to obtain gas disk masses from C$^{18}$O emission lines.


title: Searching in the Forest for Local Bayesian Optimization

abstract: Because of its sample efficiency, Bayesian optimization (BO) has become a popular approach dealing with expensive black-box optimization problems, such as hyperparameter optimization (HPO). Recent empirical experiments showed that the loss landscapes of HPO problems tend to be more benign than previously assumed, i.e. in the best case uni-modal and convex, such that a BO framework could be more efficient if it can focus on those promising local regions. In this paper, we propose BOinG, a two-stage approach that is tailored toward mid-sized configuration spaces, as one encounters in many HPO problems. In the first stage, we build a scalable global surrogate model with a random forest to describe the overall landscape structure. Further, we choose a promising subregion via a bottom-up approach on the upper-level tree structure. In the second stage, a local model in this subregion is utilized to suggest the point to be evaluated next. Empirical experiments show that BOinG is able to exploit the structure of typical HPO problems and performs particularly well on mid-sized problems from synthetic functions and HPO.


title: What Makes Online Communities 'Better'? Measuring Values, Consensus, and Conflict across Thousands of Subreddits

abstract: Making online social communities 'better' is a challenging undertaking, as online communities are extraordinarily varied in their size, topical focus, and governance. As such, what is valued by one community may not be valued by another. However, community values are challenging to measure as they are rarely explicitly stated. In this work, we measure community values through the first large-scale survey of community values, including 2,769 reddit users in 2,151 unique subreddits. Through a combination of survey responses and a quantitative analysis of public reddit data, we characterize how these values vary within and across communities. Amongst other findings, we show that community members disagree about how safe their communities are, that longstanding communities place 30.1% more importance on trustworthiness than newer communities, and that community moderators want their communities to be 56.7% less democratic than non-moderator community members. These findings have important implications, including suggesting that care must be taken to protect vulnerable community members, and that participatory governance strategies may be difficult to implement. Accurate and scalable modeling of community values enables research and governance which is tuned to each community's different values. To this end, we demonstrate that a small number of automatically quantifiable features capture a significant yet limited amount of the variation in values between communities with a ROC AUC of 0.667 on a binary classification task. However, substantial variation remains, and modeling community values remains an important topic for future work. We make our models and data public to inform community design and governance.


title: Magnetic fields in star formation: a complete compilation of all the DCF estimations

abstract: The Davis-Chandrasekhar-Fermi (DCF) method provides an indirect way to estimate the magnetic field strength from statistics of magnetic field orientations. We compile all the previous DCF estimations from polarized dust emission observations and re-calculate the magnetic field strength of the selected samples with the new DCF correction factors in Liu et al. (2021). We find the magnetic field scales with the volume density as $B \propto n^{0.57}$. However, the estimated power-law index of the observed $B-n$ relation has large uncertainties and may not be comparable to the $B-n$ relation of theoretical models. A clear trend of decreasing magnetic viral parameter (i.e., increasing mass-to-flux ratio in units of critical value) with increasing column density is found in the sample, which suggests the magnetic field dominates the gravity at lower densities but cannot compete with the gravity at higher densities. This finding also indicates that the magnetic flux is dissipated at higher column densities due to ambipolar diffusion or magnetic recennection, and the accumulation of mass at higher densities may be by mass flows along the magnetic field lines. Both sub-Alfv\'{e}nic and super-Alfv\'{e}nic states are found in the sample, with the average state being approximately trans-Alfv\'{e}nic.


title: Electronic Conduction Through Monolayer Amorphous Carbon Nano-Junctions

abstract: In molecular electronic conduction, exotic lattice morphologies often give rise to exotic behaviors. Among 2D systems, graphene is a notable example. Recently, a stable amorphous version of graphene called Monolayer Amorphous Carbon (MAC) was synthesized. MAC poses a new set of questions regarding the effects of disorder on conduction. In this Letter, we perform ensemble-level computational analysis of the coherent electronic transmission through MAC nano-fragments in search of defining characteristics. Our analysis, relying on a semi-empirical Hamiltonian (Pariser-Parr-Pople) and Landauer theory, showed that states near the Fermi energy ($E_F$) in MAC inherit partial characteristics of analogous surface states in graphene nano-fragments. Away from $E_F$, current is carried by a set of delocalized states which transition into a subset of insulating interior states at the band edges. Finally, we also found that quantum interference between the frontier orbitals is a common feature among MAC nano-fragments.


title: Convergence of Stratified MCMC Sampling of Non-Reversible Dynamics

abstract: We present a form of stratified MCMC algorithm built with non-reversible stochastic dynamics in mind. It can also be viewed as a generalization of the exact milestoning method, or form of NEUS. We prove convergence of the method under certain assumptions, with expressions for the convergence rate in terms of the process's behavior within each stratum and large scale behavior between strata. We show that the algorithm has a unique fixed point which corresponds to the invariant measure of the process without stratification. We will show how the speeds of two versions of the new algorithm, one with an extra eigenvalue problem step and one without, relate to the mixing rate of a discrete process on the strata, and the mixing probability of the process being sampled within each stratum. The eigenvalue problem version also relates to local and global perturbation results of discrete Markov chains, such as those given by Van Koten, Weare et. al.


title: On the Equity-Efficiency Trade-off in Food-Bank Network Operations

abstract: In this paper, we present a novel modelling perspective to the food-bank donation allocation problem under equity and efficiency performance measures. Using a penalty factor in the objective function, our model explicitly accounts for both efficiency and equity, simultaneously. We give the tightest lower and upper bounds of the penalty factor, which can conveniently characterize closed-form optimal solutions for the perfect efficiency and perfect equity cases. Testing our model on the full spectrum of our penalty factor, using real data from Feeding America, we demonstrate that the solutions from our model dominate those of a benchmark from the literature in terms of equity, efficiency, and utilization equity (utiloquity). Our sensitivity analysis demonstrates that the society should put its priority on helping eliminate poverty before investing on capacity expansions in charity organizations like food-banks. This will ensure that adding more capacity to the network will always lead to a decrease in the price of equity for the food-banks. On the other hand, we observed that encouraging the society towards charitability is always beneficial for the food-banks, albeit with diminishing returns. Finally, our experiments demonstrate that reducing poverty, as the most important element in achieving higher equity, is dependent on reducing demand variability, as opposed to the average level of demand.


title: Calculation of Electric Quadrupole Linestrengths for Diatomic Molecules: Application to the H2, CO, HF and O2 Molecules

abstract: We present a unified variational treatment of the electric quadrupole (E2) matrix elements, Einstein coefficients, and line strengths for general open-shell diatomic molecules in the general purpose diatomic code \Duo. Transformation relations between the Cartesian representation (typically used in electronic structure calculations) to the tensorial representation (required for spectroscopic applications) of the electric quadrupole moment components are derived. The implementation has been validated against accurate theoretical calculations and experimental measurements of quadrupole intensities of $^1\text{H}_2$ available in the literature. We also present accurate electronic structure calculations of the electric quadrupole moment functions for the $X^1\Sigma^+$ electronic states of $\text{CO}$ and $\text{HF}$ at the CCSD(T) and MRCI levels of theory, respectively, as well for the $a^1\Delta_g$ -- $b^1\Sigma_g^+$ quadrupole transition moment of $\text{O}_2$ with MRCI level of theory. Accurate infrared E2 line lists for $^{12}\text{C}^{16}\text{O}$ and $^1\text{H}^{19}\text{F}$ are provided. A demonstration of spectroscopic applications is presented by simulating E2 spectra for $^{12}\text{C}^{16}\text{O}$, $\text{H}^{19}\text{F}$ and $^{16}\text{O}_2$ (Noxon $a^1\Delta_g$ -- $b^1\Sigma_g^+$ band).


title: Physics-enhanced deep surrogates for partial differential equations

abstract: Many physics and engineering applications demand Partial Differential Equations (PDE) property evaluations that are traditionally computed with resource-intensive high-fidelity numerical solvers. Data-driven surrogate models provide an efficient alternative but come with a significant cost of training. Emerging applications would benefit from surrogates with an improved accuracy-cost tradeoff, while studied at scale. Here we present a "physics-enhanced deep-surrogate" ("PEDS") approach towards developing fast surrogate models for complex physical systems, which is described by PDEs. Specifically, a combination of a low-fidelity, explainable physics simulator and a neural network generator is proposed, which is trained end-to-end to globally match the output of an expensive high-fidelity numerical solver. Experiments on three exemplar testcases, diffusion, reaction-diffusion, and electromagnetic scattering models, show that a PEDS surrogate can be up to 3$\times$ more accurate than an ensemble of feedforward neural networks with limited data ($\approx 10^3$ training points), and reduces the training data need by at least a factor of 100 to achieve a target error of 5%. Experiments reveal that PEDS provides a general, data-driven strategy to bridge the gap between a vast array of simplified physical models with corresponding brute-force numerical solvers modeling complex systems, offering accuracy, speed, data efficiency, as well as physical insights into the process.


title: Reply to Comment on "TVOR: Finding Discrete Total Variation Outliers among Histograms"

abstract: In this paper, we respond to a critique of one of our papers previously published in this journal, entitled "TVOR: Finding Discrete Total Variation Outliers among Histograms". Our paper proposes a method for smoothness outliers detection among histograms by using the relation between their discrete total variations (DTV) and their respective sample sizes. In this response, we demonstrate point by point that, contrary to its claims, the critique has not found any mistakes or problems in our paper, either in the used datasets, methodology, or in the obtained top outlier candidates. On the contrary, the critique's claims can easily be shown to be mathematically unfounded, to directly contradict the common statistical theorems, and to go against well established demographic terms. Exactly this is done in the reply here by providing both theoretical and experimental evidence. Additionally, due to critique's complaint, a more extensive research on top outlier candidate, i.e. the Jasenovac list is conducted and in order to clear any of the critique's doubts, new evidence of its problematic nature unseen in other lists are presented. This reply is accompanied by additional theoretical explanations, simulations, and experimental results that not only confirm the earlier findings, but also present new data. The source code is at https://github.com/DiscreteTotalVariation/TVOR.


title: An Integrated Vaccination Site Selection and Dose Allocation Problem with Fairness Concerns

abstract: Fairness in vaccination is not only important from a social justice point of view, but experience has shown that a fair distribution of vaccine proves more effective in public immunization by preventing highly-concentrated infected areas to form among the population. In this paper, we address fairness from two simultaneous points of view: equity and accessibility. Equity in our setting means that as far as possible, each demand zone should receive a fair-share of the total doses available. On the other hand, accessibility means that as far as possible, each demand zone should have equal travel distance to access their assigned vaccination site.


title: Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

abstract: The CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $\alpha_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $\alpha_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $\alpha_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $\alpha_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $\alpha_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $\alpha_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.