title: Independent Action Models and Prediction of Combination Treatment Effects for Response Rate, Duration of Response and Tumor Size Change in Oncology Drug Development

abstract: An unprecedented number of new cancer targets are in development, and most are being developed in combination therapies. Early oncology development is strategically challenged in choosing the best combinations to move forward to late stage development. The most common early endpoints to be assessed in such decision-making include objective response rate, duration of response and tumor size change. In this paper, using independent-drug-action and Bliss-drug-independence concepts as a foundation, we introduce simple models to predict combination therapy efficacy for duration of response and tumor size change. These models complement previous publications using the independent action models (Palmer 2017, Schmidt 2020) to predict progression-free survival and objective response rate and serve as new predictive models to understand drug combinations for early endpoints. The models can be applied to predict the combination treatment effect for early endpoints given monotherapy data, or to estimate the possible effect of one monotherapy in the combination if data are available from the combination therapy and the other monotherapy. Such quantitative work facilitates efficient oncology drug development.


title: FLAME: Taming Backdoors in Federated Learning (Extended Version 1)

abstract: Federated Learning (FL) is a collaborative machine learning approach allowing participants to jointly train a model without having to share their private, potentially sensitive local datasets with others. Despite its benefits, FL is vulnerable to backdoor attacks, in which an adversary injects manipulated model updates into the model aggregation process so that the resulting model will provide targeted false predictions for specific adversary-chosen inputs. Proposed defenses against backdoor attacks based on detecting and filtering out malicious model updates consider only very specific and limited attacker models, whereas defenses based on differential privacy-inspired noise injection significantly deteriorate the benign performance of the aggregated model. To address these deficiencies, we introduce FLAME, a defense framework that estimates the sufficient amount of noise to be injected to ensure the elimination of backdoors while maintaining the model performance. To minimize the required amount of noise, FLAME uses a model clustering and weight clipping approach. Our evaluation of FLAME on several datasets stemming from application areas including image classification, word prediction, and IoT intrusion detection demonstrates that FLAME removes backdoors effectively with a negligible impact on the benign performance of the models. Furthermore, following the considerable attention that our research has received after its presentation at USENIX SEC 2022, FLAME has become the subject of numerous investigations proposing diverse attack methodologies in an attempt to circumvent it. As a response to these endeavors, we provide a comprehensive analysis of these attempts. Our findings show that these papers (e.g., 3DFed [36]) have not fully comprehended nor correctly employed the fundamental principles underlying FLAME, i.e., our defense mechanism effectively repels these attempted attacks.


title: Navigation Framework for a Hybrid Steel Bridge Inspection Robot

abstract: Autonomous navigation is essential for steel bridge inspection robot to monitor and maintain the working condition of steel bridges. Majority of existing robotic solutions requires human support to navigate the robot doing the inspection. In this paper, a navigation framework is proposed for ARA robot [1], [2] to run on mobile mode. In this mode, the robot needs to cross and inspect all the available steel bars. The most significant contributions of this research are four algorithms, which can process the depth data, segment it into clusters, estimate the boundaries, construct a graph to represent the structure, generate a shortest inspection path with any starting and ending points, and determine available robot configuration for path planning. Experiments on steel bridge structures setup highlight the effective performance of the algorithms, and the potential to apply to the ARA robot to run on real bridge structures. We released our source code in Github for the research community to use.


title: Selberg's Central Limit Theorem for $L$-functions of level aspect, II

abstract: In this paper, we have proved Selberg's Central Limit Theorem for $GL(3)$ $L$-functions associated with the Hecke-Maass cusp form $f$. Moreover, we have proved the independence of the automorphic $L$-functions.


title: Handling many conversions per click in modeling delayed feedback

abstract: Predicting the expected value or number of post-click conversions (purchases or other events) is a key task in performance-based digital advertising. In training a conversion optimizer model, one of the most crucial aspects is handling delayed feedback with respect to conversions, which can happen multiple times with varying delay. This task is difficult, as the delay distribution is different for each advertiser, is long-tailed, often does not follow any particular class of parametric distributions, and can change over time. We tackle these challenges using an unbiased estimation model based on three core ideas. The first idea is to split the label as a sum of labels with different delay buckets, each of which trains only on mature label, the second is to use thermometer encoding to increase accuracy and reduce inference cost, and the third is to use auxiliary information to increase the stability of the model and to handle drift in the distribution.


title: TryOnGAN: Body-Aware Try-On via Layered Interpolation

abstract: Given a pair of images-target person and garment on another person-we automatically generate the target person in the given garment. Previous methods mostly focused on texture transfer via paired data training, while overlooking body shape deformations, skin color, and seamless blending of garment with the person. This work focuses on those three components, while also not requiring paired data training. We designed a pose conditioned StyleGAN2 architecture with a clothing segmentation branch that is trained on images of people wearing garments. Once trained, we propose a new layered latent space interpolation method that allows us to preserve and synthesize skin color and target body shape while transferring the garment from a different person. We demonstrate results on high resolution 512x512 images, and extensively compare to state of the art in try-on on both latent space generated and real images.


title: Scalable Parallel Linear Solver for Compact Banded Systems on Heterogeneous Architectures

abstract: A scalable algorithm for solving compact banded linear systems on distributed memory architectures is presented. The proposed method factorizes the original system into two levels of memory hierarchies, and solves it using parallel cyclic reduction on both distributed and shared memory. This method has a lower communication footprint across distributed memory partitions compared to conventional algorithms involving data transpose or re-partitioning. The algorithm developed in this work is generalized to cyclic compact banded systems with flexible data decompositions. For cyclic compact banded systems, the method is a direct solver with a deterministic operation and communication counts depending on the matrix size, its bandwidth, and the partition strategy. The implementation and runtime configuration details are discussed for performance optimization. Scalability is demonstrated on the linear solver as well as on a representative fluid mechanics application problem, in which the dominant computational cost is solving the cyclic tridiagonal linear systems of compact numerical schemes on a 3D periodic domain. The algorithm is particularly useful for solving the linear systems arising from the application of compact finite difference operators to a wide range of partial differential equation problems, such as but not limited to the numerical simulations of compressible turbulent flows, aeroacoustics, elastic-plastic wave propagation, and electromagnetics. It alleviates obstacles to their use on modern high performance computing hardware, where memory and computational power are distributed across nodes with multi-threaded processing units.


title: COVID19-HPSMP: COVID-19 Adopted Hybrid and Parallel Deep Information Fusion Framework for Stock Price Movement Prediction

abstract: The novel of coronavirus (COVID-19) has suddenly and abruptly changed the world as we knew at the start of the 3rd decade of the 21st century. Particularly, COVID-19 pandemic has negatively affected financial econometrics and stock markets across the globe. Artificial Intelligence (AI) and Machine Learning (ML)-based prediction models, especially Deep Neural Network (DNN) architectures, have the potential to act as a key enabling factor to reduce the adverse effects of the COVID-19 pandemic and future possible ones on financial markets. In this regard, first, a unique COVID-19 related PRIce MOvement prediction (COVID19 PRIMO) dataset is introduced in this paper, which incorporates effects of social media trends related to COVID-19 on stock market price movements. Afterwards, a novel hybrid and parallel DNN-based framework is proposed that integrates different and diversified learning architectures. Referred to as the COVID-19 adopted Hybrid and Parallel deep fusion framework for Stock price Movement Prediction (COVID19-HPSMP), innovative fusion strategies are used to combine scattered social media news related to COVID-19 with historical mark data. The proposed COVID19-HPSMP consists of two parallel paths (hence hybrid), one based on Convolutional Neural Network (CNN) with Local/Global Attention modules, and one integrated CNN and Bi-directional Long Short term Memory (BLSTM) path. The two parallel paths are followed by a multilayer fusion layer acting as a fusion centre that combines localized features. Performance evaluations are performed based on the introduced COVID19 PRIMO dataset illustrating superior performance of the proposed framework.


title: Theory and Applications of Financial Chaos Index

abstract: We develop a new stock market index that captures the chaos existing in the market by measuring the mutual changes of asset prices. This new index relies on a tensor-based embedding of the stock market information, which in turn frees it from the restrictive value- or capitalization-weighting assumptions that commonly underlie other various popular indexes. We show that our index is a robust estimator of the market volatility which enables us to characterize the market by performing the task of segmentation with a high degree of reliability. In addition, we analyze the dynamics and kinematics of the realized market volatility as compared to the implied volatility by introducing a time-dependent dynamical system model. Our computational results which pertain to the time period from January 1990 to December 2019 imply that there exist a bidirectional causal relation between the processes underlying the realized and implied volatility of the stock market within the given time period, where it is shown that the later has a stronger causal effect on the former as compared to the opposite. This result connotes that the implied volatility of the market plays a key role in characterization of the market's realized volatility.


title: Hyperboost: Hyperparameter Optimization by Gradient Boosting surrogate models

abstract: Bayesian Optimization is a popular tool for tuning algorithms in automatic machine learning (AutoML) systems. Current state-of-the-art methods leverage Random Forests or Gaussian processes to build a surrogate model that predicts algorithm performance given a certain set of hyperparameter settings. In this paper, we propose a new surrogate model based on gradient boosting, where we use quantile regression to provide optimistic estimates of the performance of an unobserved hyperparameter setting, and combine this with a distance metric between unobserved and observed hyperparameter settings to help regulate exploration. We demonstrate empirically that the new method is able to outperform some state-of-the art techniques across a reasonable sized set of classification problems.


title: High-energy emission from tidal disruption events in active galactic nuclei

abstract: Tidal disruption events (TDEs) taking place in active galactic nuclei (AGNs) are different from ordinary TDEs. In these events, the returning tidal debris stream drills through the pre-existing AGN accretion disk near the stream pericenter, destroying the inner disk in the process, and then intersects with the disk a second time at radii ranging from a few to hundreds of times the pericenter distance. The debris dynamics of such TDEs, and hence their appearance, are distinct from ordinary TDEs. Here we explore the observational signatures of this "second impact" of the stream with the disk. Strong shocks form as the dilute stream is stopped by the denser disk. Compton cooling of the shocked material produces hard X-rays, even soft gamma-rays, with most of the energy emitted between ~10 keV and 1 MeV. The luminosity follows the mass-return rate, peaking between ~$10^{42}$ and $10^{44}$ erg/s. The X-ray hardness and the smoothness of the light curve provide possible means for distinguishing the second impact from ordinary AGN flares, which exhibit softer spectra and more irregular light curves.


title: The quantum paraelectric phase of SrTiO$_3$ from first principles

abstract: We demonstrate how the quantum paraelectric ground state of SrTiO$_3$ can be accessed via a microscopic $ab~initio$ approach based on density functional theory. At low temperature the quantum fluctuations are strong enough to stabilize the paraelectric phase even though a classical description would predict a ferroelectric phase. We find that accounting for quantum fluctuations of the lattice and for the strong coupling between the ferroelectric soft mode and lattice elongation is necessary to achieve quantitative agreement with experimental frequency of the ferroelectric soft mode. The temperature dependent properties in SrTiO$_3$ are also well captured by the present microscopic framework.


title: The Conundrum of Strongly Coupled Supercurrent Flow in Both Under- and Over-doped Bi-2212 Round Wires

abstract: Understanding what makes Bi$_2$Sr$_2$Ca$_1$Cu$_2$O$_x$ (Bi-2212) the only high critical current density ($J_c$), high temperature superconductor (HTS) capable of being made as a round wire (RW) is important intellectually because high $J_c$ RW Bi-2212 breaks the paradigm that forces biaxially textured REBCO and uniaxially textured (Bi,Pb)$_2$Sr$_2$Ca$_2$Cu$_3$O$_x$ (Bi-2223) into tape geometries that reproduce the strong anisotropy of the native crystal structure and force expensive fabrication routes to ensure the best possible texture with minimum density of high angle grain boundaries. The biaxial growth texture of Bi-2212 developed during a partial melt heat treatment should favor high $J_c$, even though its $\sim$15$^{\circ}$ full width at half maximum (FWHM) grain-to-grain misorientation is well beyond the commonly accepted strong-coupling range. Its ability to be strongly overdoped should be valuable too, since underdoped cuprate grain boundaries are widely believed to be weakly linked. Accordingly, we here study property changes after oxygen underdoping the optimized, overdoped wire. While $J_c$ and vortex pinning diminish significantly in underdoped wires, we were not able to develop the prominent weak-link signature (a hysteretic $J_c$(H) characteristic) evident in even the very best Bi-2223 tapes with a $\sim$ 15$^{\circ}$ FWHM uniaxial texture. We attribute the high $J_c$ and lack of weak link signature in our Bi-2212 round wires to the high-aspect ratio, large-grain, basal-plane-faced grain morphology produced by partial-melt processing of Bi-2212 which enables $c$-axis Brick-Wall current flow when $ab$-plane transport is blocked. We conclude that the presently optimized biaxial texture of Bi-2212 intrinsically constitutes a strongly coupled current path, regardless of its oxygen doping state.


title: Most Elliptic Curves over Global Function Fields are Torsion Free

abstract: Given an elliptic curve $E$ over a global function field $K$, the Galois action on the $n$-torsion points of $E$ gives rise to a mod-n Galois representation $\rho_{E,n}$. For $K$ satisfying some mild conditions, we show that the set of $E$ for which $\rho_{E,n}$ is as large as possible for all $n$, has density $1$.


title: Magnet design optimization of a 100 MeV separated sector medical cyclotron and its injection line

abstract: This paper presents the magnetic design of a 100MeV Separated Sector Medical Cyclotron (SSMC) as well as the optimization of the injection line magnets where will be used to transport the 14MeV proton beam from the Azimuthally Varying Field (AVF) cyclotron to the SSMC. The study demonstrates that the isochronous magnetic field with a tolerance around 10-4 T can be obtained all along prescribed path by placing three sets of correction trim coils in the pole tips. The result of Betatron oscillations is sufficient to verify that the focusing forces of the magnetic field can hold particles close to the median plane of the magnet. In order to inject the beam by minimum loss from a 14MeV Cyclotron to SSMC, it is essential for the beam to fit into the acceptance ellipse of the separate sector cyclotron. The characteristics of the beam injection system is calculated and optimized by Trace-3D beam dynamic code.


title: Magnetic evaluation of a developed permanent magnet undulator as a linac-based THz coherent radiation source

abstract: Due to the various applications of THz sources, we developed a pure permanent magnet (PPM) undulator with 8 periods and min/max values of gap equal to 15/25, to use at the end of our current 10 MeV linear accelerator (linac) in order to produce THz coherent radiation. We used RADIA for simulation and calculation of our magnet and B2E code for calculation of its radiation harmonics. After fabrication, we also measured the magnetic field via high precision Hall probes. Calculations and measurements are in good agreements. These calculations show that it is possible to reach photon energies from 3 meV to 23 meV in the THz region.


title: Predicting CEO Compensation in Non-Controlled Public Corporations with the Canonical Regression Quantile Method

abstract: The use of the Canonical Regression Quantiles Index proved that non-controlled companies that engage in long-term operational and financial goals post superior future performance. The Index indicates that current CEO compensation influences future performance. The Index provides a method for determining CEO pay for the next 1-2 year and is a useful method to distinguish over/underpaid CEOs as an unbiased alternative to the peer groups comparison used by most compensation consultants. This determination is statistically weak, but future research using the Canonical Regression Quantiles with a larger data set may lead to increased sensitivity and a powerful unbiased method for replacing compensation consultants who are responsible for the decoupling of CEO compensation and corporate performance.


title: Cosmological constant as a finite temperature effect

abstract: We reexamine the cosmological constant problem in a finite temperature setup and propose an intriguing possibility of carrying out perturbative analysis by employing a renormalization scheme in which the renormalized Higgs mass (or resummed mass, to be more precise) is taken to be on the order of the CMB temperature. Our proposal hinges on the fact that although the physical value of the cosmological constant does not depend on one's renormalization scheme, whether or not a fine tuning is involved does. The cosmological constant problem is avoided in the sense that the renormalization process no longer requires finetuning. This is achieved essentially by renormalization scheme-independence of a physical quantity, which in turn is assured by bare perturbation theory. The proposal shifts the cosmological constant problem to a peculiarity of the consequent perturbation series for the Higgs mass (and other massive sectors of the Standard Model); the peculiarity is interpreted as an indicator of new physics after the expected mathematical structure of the series is scrutinized. Finite-temperature-induced complexification of the effective potential is observed and its interpretation is given. A consistency check in the cosmology context is suggested.


title: Implementation of conformal digital metasurfaces for THz polarimetric sensing

abstract: Monitoring and controlling the state of polarization of electromagnetic waves is of significant interest for various basic and practical applications such as linear position sensing and medical imaging. Here, we propose the first conformal digital metamaterial absorber to detect the polarization state of THz incident waves. The proposed polarimeter is capable of characterizing four independent polarization states of (TE, TM, $\pm {45^\circ}$ linear, and RCP/LCP) by observing the reflectivity of the structure with respect to the x- and y-direction. Besides, the proposed structure displays a strong absorptivity above 90\% up to the incidence angle of $50^{\circ}$ for oblique incident waves with different polarizations. By mere changing the bias voltage of two orthogonal VO2 microwires via two independent computer-programmed multichannel DC network, distinct conditions for reflected waves occurs under excitations of different polarizations, whereby the polarization state of the incident wave may readily be estimated. We believe that the proposed metasurface-based polarimeter can pave the way for polarization detection applications on curved surfaces.


title: Enumerating Labeled Graphs that Realize a Fixed Degree Sequence

abstract: A finite non-increasing sequence of positive integers $d = (d_1\geq \cdots\geq d_n)$ is called a degree sequence if there is a graph $G = (V,E)$ with $V = \{v_1,\ldots,v_n\}$ and $deg(v_i)=d_i$ for $i=1,\ldots,n$. In that case we say that the graph $G$ realizes the degree sequence $d$. We show that the exact number of labeled graphs that realize a fixed degree sequence satisfies a simple recurrence relation. Using this relation, we then obtain a recursive algorithm for the exact count. We also show that in the case of regular graphs the complexity of our algorithm is better than the complexity of the same enumeration that uses generating functions.


title: Continuous-variable error correction for general Gaussian noises

abstract: Quantum error correction is essential for robust quantum information processing with noisy devices. As bosonic quantum systems play a crucial role in quantum sensing, communication, and computation, it is important to design error correction codes suitable for these systems against various different types of noises. While most efforts aim at protecting qubits encoded into the infinite dimensional Hilbert space of a bosonic mode, [Phys. Rev. Lett. 125, 080503 (2020)] proposed an error correction code to maintain the infinite-dimensional-Hilbert-space nature of bosonic systems by encoding a single bosonic mode into multiple bosonic modes. Enabled by Gottesman-Kitaev-Preskill states as ancilla, the code overcomes the no-go theorem of Gaussian error correction. In this work, we generalize the error correction code to the scenario with general correlated and heterogeneous Gaussian noises, including memory effects. We introduce Gaussian pre-processing and post-processing to convert the general noise model to an independent but heterogeneous collection of additive white Gaussian noise channels and then apply concatenated codes in an optimized manner. To evaluate the performance, we develop a theory framework to enable the efficient calculation of the noise standard deviation after the error correction, despite the non-Gaussian nature of the codes. Our code provides the optimal scaling of the residue noise standard deviation with the number of modes and can be widely applied to distributed sensor-networks, network communication and composite quantum memory systems.


title: Bounded generation for congruence subgroups of ${\rm Sp}_4(R)$

abstract: This paper describes a bounded generation result concerning the minimal natural number $K$ such that for $Q(C_2,2R):=\{A\varepsilon_{\phi}(2x)A^{-1}|x\in R,A\in{\rm Sp}_4(R),\phi\in C_2\}$, one has $N_{C_2,2R}=\{X_1\cdots X_K|\forall 1\leq i\leq K:X_i\in Q(C_2,2R)\}$ for rings of algebraic integers $R$ and the principal congruence subgroup $N_{C_2,2R}$ in ${\rm Sp}_4(R).$ This gives an explicit version of an abstract bounded generation result of a similar type as presented by Morris. Furthermore, the result presented does not depend on several number-theoretic quantities unlike Morris' result. Using this bounded generation result, we further give explicit bounds for the strong boundedness of ${\rm Sp}_4(R)$ for certain examples of rings $R,$ thereby giving explicit versions of results in an earlier paper. We further give a classification of normally generating subsets of ${\rm Sp}_4(R)$ for $R$ a ring of algebraic integers.


title: Lipschitz geometry and combinatorics of abnormal surface germs

abstract: We study outer Lipschitz geometry of real semialgebraic or, more general, definable in a polynomially bounded o-minimal structure over the reals, surface germs. In particular, any definable H\"older triangle is either Lipschitz normally embedded or contains some "abnormal" arcs. We show that abnormal arcs constitute finitely many "abnormal zones" in the space of all arcs, and investigate geometric and combinatorial properties of abnormal surface germs. We establish a strong relation between geometry and combinatorics of abnormal H\"older triangles.


title: Review of Core Process Representation in Power System Operational Models: Gaps, Challenges, and Opportunities for Multisector Dynamics Research

abstract: Power grid operations increasingly interact with environmental systems and human systems such as transportation, agriculture, the economy, and financial markets. Our objective is to discuss the modelling gaps and opportunities to advance the science for multisector adaptation and tradeoffs. We focus on power system operational models, which typically represent key physical and economic aspects of grid operations over days to a year and assume a fixed power grid infrastructure. Due to computational burden, models are typically customized to reflect regional resource opportunities, data availability, and applications of interest. While there are model intercomparison papers, there is however no model-agnostic characterization and systematic overview of the state-of-the-art process representations in operational power system models. To address our objective, we conceptualize power system operational models with four core processes: physical grid assets (generation, transmission, loads, and storage), model objectives and purpose, institutions and decision agents, and performance metrics. We taxonomize the representations of these core processes based on a review of 23 existing open-source and commercial models. As we acknowledge the computational burden of certain representations, we leverage this taxonomy to describe tradeoffs in process fidelity and tractability that have been adopted by the research community to address interactions between the power grid and hydrometeorological uncertainties, global change, and/or technological innovation. The core process taxonomy along with the existing computational tradeoffs are used to identify technical gaps and recommend future model development needs and research directions to better represent power grid operations as part of integrated multisector dynamics modeling and interdisciplinary research.


title: Statistical challenges in the analysis of sequence and structure data for the COVID-19 spike protein

abstract: As the major target of many vaccines and neutralizing antibodies against SARS-CoV-2, the spike (S) protein is observed to mutate over time. In this paper, we present statistical approaches to tackle some challenges associated with the analysis of S-protein data. We build a Bayesian hierarchical model to study the temporal and spatial evolution of S-protein sequences, after grouping the sequences into representative clusters. We then apply sampling methods to investigate possible changes to the S-protein's 3-D structure as a result of commonly observed mutations. While the increasing spread of D614G variants has been noted in other research, our results also show that the co-occurring mutations of D614G together with S477N or A222V may spread even more rapidly, as quantified by our model estimates.


title: Demand Forecasting for Platelet Usage: from Univariate Time Series to Multivariate Models

abstract: Platelet products are both expensive and have very short shelf lives. As usage rates for platelets are highly variable, the effective management of platelet demand and supply is very important yet challenging. The primary goal of this paper is to present an efficient forecasting model for platelet demand at Canadian Blood Services (CBS). To accomplish this goal, four different demand forecasting methods, ARIMA (Auto Regressive Moving Average), Prophet, lasso regression (least absolute shrinkage and selection operator) and LSTM (Long Short-Term Memory) networks are utilized and evaluated. We use a large clinical dataset for a centralized blood distribution centre for four hospitals in Hamilton, Ontario, spanning from 2010 to 2018 and consisting of daily platelet transfusions along with information such as the product specifications, the recipients' characteristics, and the recipients' laboratory test results. This study is the first to utilize different methods from statistical time series models to data-driven regression and a machine learning technique for platelet transfusion using clinical predictors and with different amounts of data. We find that the multivariate approaches have the highest accuracy in general, however, if sufficient data are available, a simpler time series approach such as ARIMA appears to be sufficient. We also comment on the approach to choose clinical indicators (inputs) for the multivariate models.


title: Interpolation by holomorphic maps from the disc to the tetrablock

abstract: The tetrablock is the set $$ \mathcal{E}=\{x \in \mathbb{C}^3: \quad 1-x_1z-x_2w+x_3z w \neq 0 \quad whenever \quad |z|\leq 1, |w|\leq 1\}. $$ The closure of $\mathcal{E}$ is denoted by $\overline{\mathcal{E}}$. A tetra-inner function is an analytic map $x$ from the unit disc $ \mathbb{D} $ to $\overline{\mathcal{E}}$ such that, for almost all points $\lambda$ of the unit circle $ \mathbb{T}$, \[ \lim_{r\uparrow 1} x(r \lambda) \mbox{ exists and lies in } b \overline{\mathcal{E}}, \] where $b \overline{\mathcal{E}}$ denotes the distinguished boundary of $\overline{\mathcal{E}}$. There is a natural notion of degree of a rational tetra-inner function $ x$; it is simply the topological degree of the continuous map $ x|_\mathbb{T} $ from $ \mathbb{T} $ to $ b \overline{\mathcal{E}} $. In this paper we give a prescription for the construction of a general rational tetra-inner function of degree $n$. The prescription exploits a known construction of the finite Blaschke products of given degree which satisfy some interpolation conditions with the aid of a Pick matrix formed from the interpolation data. It is known that if $x= (x_1, x_2, x_3)$ is a rational tetra-inner function of degree $n$, then $x_1 x_2 - x_3$ either is identically $0$ or has precisely $n$ zeros in the closed unit disc $\overline{\mathbb{D}}$, counted with multiplicity. It turns out that a natural choice of data for the construction of a rational tetra-inner function $x= (x_1, x_2, x_3)$ consists of the points in $\overline{\mathbb{D}}$ for which $x_1 x_2 - x_3=0$ and the values of $x$ at these points.


title: Directed mixed membership stochastic blockmodel

abstract: Mixed membership problem for undirected network has been well studied in network analysis recent years. However, the more general case of mixed membership for directed network in which nodes can belong to multiple communities remains a challenge. Here, we propose an interpretable and identifiable model: directed mixed membership stochastic blockmodel (DiMMSB) for directed mixed membership networks. DiMMSB allows that row nodes and column nodes of the adjacency matrix can be different and these nodes may have distinct community structure in a directed network. We also develop an efficient spectral algorithm called DiSP designed based on simplex structures inherent in the left and right singular vectors of the population adjacency matrix to estimate the mixed memberships for both row nodes and column nodes in a directed network. We show that DiSP is asymptotically consistent under mild conditions by providing error bounds for the inferred membership vectors of each row node and each column node using delicate spectral analysis. Numerical results on computer-generated directed mixed membership networks support our theoretical findings and show that our DiSP outperforms its competitor in both error rates and run-time. Applications of DiSP to real-world directed networks demonstrate the advantages of DiSP in studying the asymmetric structure of directed networks.


title: Coding for Distributed Multi-Agent Reinforcement Learning

abstract: This paper aims to mitigate straggler effects in synchronous distributed learning for multi-agent reinforcement learning (MARL) problems. Stragglers arise frequently in a distributed learning system, due to the existence of various system disturbances such as slow-downs or failures of compute nodes and communication bottlenecks. To resolve this issue, we propose a coded distributed learning framework, which speeds up the training of MARL algorithms in the presence of stragglers, while maintaining the same accuracy as the centralized approach. As an illustration, a coded distributed version of the multi-agent deep deterministic policy gradient(MADDPG) algorithm is developed and evaluated. Different coding schemes, including maximum distance separable (MDS)code, random sparse code, replication-based code, and regular low density parity check (LDPC) code are also investigated. Simulations in several multi-robot problems demonstrate the promising performance of the proposed framework.


title: Cellular Automata on Quantum Annealing Systems

abstract: We present herein an introduction to implementing 2-color cellular automata on quantum annealing systems, such as the D-Wave quantum computer. We show that implementing nearest-neighbor cellular automata is possible. We present an implementation of Wolfram's cellular automata Rule 110, which has previously been shown to be a universal Turing machine, as a QUBO suitable for use on quantum annealing systems. We demonstrate back-propagation of cellular automata rule sets to determine initial cell states for a desired later system state. We show 2-D 2-color cellular automata, such as Conway's Game of Life, can be expressed for quantum annealing systems.


title: Parallel Hyperedge Replacement Grammars

abstract: In 2018, it was shown that all finitely generated virtually Abelian groups have multiple context-free word problems, and it is still an open problem as to where to precisely place the word problems of hyperbolic groups in the formal language hierarchy. Motivated by this, we introduce a new language class, the parallel hyperedge replacement string languages, containing all multiple context-free and ET0L languages. We show that parallel hyperedge replacement grammars can be "synchronised", which allows us to establish many useful formal language closure results relating to both the hypergraph and string languages generated by various families of parallel hyperedge replacement grammars, laying the foundations for future work in this area.


title: A Deterministic Parallel APSP Algorithm and its Applications

abstract: In this paper we show a deterministic parallel all-pairs shortest paths algorithm for real-weighted directed graphs. The algorithm has $\tilde{O}(nm+(n/d)^3)$ work and $\tilde{O}(d)$ depth for any depth parameter $d\in [1,n]$. To the best of our knowledge, such a trade-off has only been previously described for the real-weighted single-source shortest paths problem using randomization [Bringmann et al., ICALP'17]. Moreover, our result improves upon the parallelism of the state-of-the-art randomized parallel algorithm for computing transitive closure, which has $\tilde{O}(nm+n^3/d^2)$ work and $\tilde{O}(d)$ depth [Ullman and Yannakakis, SIAM J. Comput. '91]. Our APSP algorithm turns out to be a powerful tool for designing efficient planar graph algorithms in both parallel and sequential regimes. One notable ingredient of our parallel APSP algorithm is a simple deterministic $\tilde{O}(nm)$-work $\tilde{O}(d)$-depth procedure for computing $\tilde{O}(n/d)$-size hitting sets of shortest $d$-hop paths between all pairs of vertices of a real-weighted digraph. Such hitting sets have also been called $d$-hub sets. Hub sets have previously proved especially useful in designing parallel or dynamic shortest paths algorithms and are typically obtained via random sampling. Our procedure implies, for example, an $\tilde{O}(nm)$-time deterministic algorithm for finding a shortest negative cycle of a real-weighted digraph. Such a near-optimal bound for this problem has been so far only achieved using a randomized algorithm [Orlin et al., Discret. Appl. Math. '18].


title: 4 vs 7 sparse undirected unweighted Diameter is SETH-hard at time $n^{4/3}$

abstract: We show, assuming the Strong Exponential Time Hypothesis, that for every $\varepsilon > 0$, approximating undirected unweighted Diameter on $n$-vertex $n^{1+o(1)}$-edge graphs within ratio $7/4 - \varepsilon$ requires $n^{4/3 - o(1)}$ time. This is the first result that conditionally rules out a near-linear time $5/3$-approximation for undirected Diameter.


title: A Note on Rough Set Algebra and Core Regular Double Stone Algebras

abstract: Rough Set Theory (RST), first introduced by Pawlak in 1982, is an approach for dealing with information systems where knowledge is uncertain or incomplete.\cite{Pawlak} It is of fundamental importance in many subfields of artificial intelligence and cognitive science.\cite{RSTppf} Given a universe $U$ with an equivalence relation $\theta$, the pair $\langle U,\theta\rangle$ is referred to as an information system and we denote its collection of rough sets $R_\theta$. In our main Theorem we show $R_\theta$ with $|\theta_u| > 1\ \forall\ u \in U$ to be isomorphic to core regular double Stone algebras, CRDSA, that are complete and atomic, and that the crisp, or definable, sets form a complete atomistic Boolean algebra. These guarantees of infimum/supremeum for arbitrary subsets and formulations in terms of fundamental elements are likely useful if dealing with equivalence relations with an infinite number of partitions, such as projective Hilbert spaces. We further derive that every CRDSA is isomorphic to a subalgebra of a principal rough set algebra, $R_\theta$, for some approximation space $\langle U,\theta \rangle$. In our main Corollary we show explicitly how to embed $R_\theta$ into the CRDSA and first demonstrate by extending the culminating finite example of \cite{RCRDSA}. As our capstone, we consider the projective Hilbert space of complex numbers, $\mathbb{C}$ and show, among other things, the power set of the set of pure states is a complete, atomistic Boolean algebra. In closing, we suggest other Quantum relevant applications that may be useful, such as Hilbert spaces of operators


title: Hilbert's 17th problem in free skew fields

abstract: This paper solves the rational noncommutative analog of Hilbert's 17th problem: if a noncommutative rational function is positive semidefinite on all tuples of hermitian matrices in its domain, then it is a sum of hermitian squares of noncommutative rational functions. This result is a generalization and culmination of earlier positivity certificates for noncommutative polynomials or rational functions without hermitian singularities. More generally, a rational Positivstellensatz for free spectrahedra is given: a noncommutative rational function is positive semidefinite or undefined at every matricial solution of a linear matrix inequality $L\succeq0$ if and only if it belongs to the rational quadratic module generated by $L$. The essential intermediate step towards this Positivstellensatz for functions with singularities is an extension theorem for invertible evaluations of linear matrix pencils.


title: (Non)local logistic equations with Neumann conditions

abstract: We consider here a problem of population dynamics modeled on a logistic equation with both classical and nonlocal diffusion, possibly in combination with a pollination term. The environment considered is a niche with zero-flux, according to a new type of Neumann condition. We discuss the situations that are more favorable for the survival of the species, in terms of the first positive eigenvalue. Quite surprisingly, the eigenvalue analysis for the one dimensional case is structurally different than the higher dimensional setting, and it sensibly depends on the nonlocal character of the dispersal. The mathematical framework of this problem takes into consideration the equation $$ -\alpha\Delta u +\beta(-\Delta)^su =(m-\mu u)u+\tau\;J\star u \qquad{\mbox{in }}\; \Omega,$$ where $m$ can change sign. This equation is endowed with a set of Neumann condition that combines the classical normal derivative prescription and the nonlocal condition introduced in [S. Dipierro, X. Ros-Oton, E. Valdinoci, Rev. Mat. Iberoam. (2017)]. We will establish the existence of a minimal solution for this problem and provide a throughout discussion on whether it is possible to obtain non-trivial solutions (corresponding to the survival of the population). The investigation will rely on a quantitative analysis of the first eigenvalue of the associated problem and on precise asymptotics for large lower and upper bounds of the resource. In this, we also analyze the role played by the optimization strategy in the distribution of the resources, showing concrete examples that are unfavorable for survival, in spite of the large resources that are available in the environment.


title: Analyzing the Stability of Non-coplanar Circumbinary Planets using Machine Learning

abstract: Exoplanet detection in the past decade by efforts including NASA's Kepler and TESS missions has discovered many worlds that differ substantially from planets in our own Solar system, including more than 400 exoplanets orbiting binary or multi-star systems. This not only broadens our understanding of the diversity of exoplanets, but also promotes our study of exoplanets in the complex binary and multi-star systems and provides motivation to explore their habitability. In this study, we analyze orbital stability of exoplanets in non-coplanar circumbinary systems using a numerical simulation method, with which a large number of circumbinary planet samples are generated in order to quantify the effects of various orbital parameters on orbital stability. We also train a machine learning model that can quickly determine the stability of the circumbinary planetary systems. Our results indicate that larger inclinations of the planet tend to increase the stability of its orbit, but change in the planet's mass range between Earth and Jupiter has little effect on the stability of the system. In addition, we find that Deep Neural Networks (DNNs) have higher accuracy and precision than other machine learning algorithms.


title: On a family of $C^*$-subalgebras of Cuntz-Krieger algebras

abstract: In this paper, we study a family of $C^*$-subalgebras defined by fixed points of generalized gauge actions of a Cuntz-Krieger algebra, by introducing a family of \'etale groupoids whose associated $C^*$-algebras are these $C^*$-subalgebras. We know that topological conjugacy classes of one-sided topological Markov shifts are characterized in terms of the isomorphism classes of these \'etale groupoids.


title: Braid Loops with infinite monodromy on the Legendrian contact DGA

abstract: We present the first examples of elements in the fundamental group of the space of Legendrian links in the standard contact 3-sphere whose action on the Legendrian contact DGA is of infinite order. This allows us to construct the first families of Legendrian links that can be shown to admit infinitely many Lagrangian fillings by Floer-theoretic techniques. These families include the first known Legendrian links with infinitely many fillings that are not rainbow closures of positive braids, and the smallest Legendrian link with infinitely many fillings known to date. We discuss how to use our examples to construct other links with infinitely many fillings, in particular giving the first Floer-theoretic proof that Legendrian (n,m) torus links have infinitely many Lagrangian fillings, if n is greater than 3 and m greater than 6, or (n,m)=(4,4),(4,5). In addition, for any given higher genus, we construct a Weinstein 4-manifold homotopic to the 2-sphere whose wrapped Fukaya category can distinguish infinitely many exact closed Lagrangian surfaces of that genus. A key technical ingredient behind our results is a new combinatorial formula for decomposable cobordism maps between Legendrian contact DGAs with integer (group ring) coefficients.


title: Integration of Renewable Generators in Synthetic Electric Grids for Dynamic Analysis

abstract: This paper presents a method to better integrate dynamic models for renewable resources into synthetic electric grids. An automated dynamic models assignment process is proposed for wind and solar generators. A realistic composition ratio for different types of wind turbine generators (WTG) is assigned to each wind generator. Statistics summarized from real electric grid data form the bases in assigning proper models with reasonable parameters to each WTG. A similar process is used to assign appropriate models and parameters to each photovoltaic (PV) generator. Multiple control strategies of the renewable resources are considered and tested in case studies. Two large-scale synthetic network test cases are used as examples of modeling the dynamics of renewable generators. Several transient stability metrics are adopted to assess the stability level after being subject to N-1 contingency event. Representative contingency events are given to demonstrate the performance of the synthetic renewable generator models.


title: Scale-free tree network with an ultra-large diameter

abstract: Scale-free networks are prevalently observed in a great variety of complex systems, which triggers various researches relevant to networked models of such type. In this work, we propose a family of growth tree networks $\mathcal{T}_{t}$, which turn out to be scale-free, in an iterative manner. As opposed to most of published tree models with scale-free feature, our tree networks have the power-law exponent $\gamma=1+\ln5/\ln2$ that is obviously larger than $3$. At the same time, "small-world" property can not be found particularly because models $\mathcal{T}_{t}$ have an ultra-large diameter $D_{t}$ (i.e., $D_{t}\sim|\mathcal{T}_{t}|^{\ln3/\ln5}$) and a greater average shortest path length $\langle\mathcal{W}_{t}\rangle$ (namely, $\langle\mathcal{W}_{t}\rangle\sim|\mathcal{T}_{t}|^{\ln3/\ln5}$) where $|\mathcal{T}_{t}|$ represents vertex number. Next, we determine Pearson correlation coefficient and verify that networks $\mathcal{T}_{t}$ display disassortative mixing structure. In addition, we study random walks on tree networks $\mathcal{T}_{t}$ and derive exact solution to mean hitting time $\langle\mathcal{H}_{t}\rangle$. The results suggest that the analytic formula for quantity $\langle\mathcal{H}_{t}\rangle$ as a function of vertex number $|\mathcal{T}_{t}|$ shows a power-law form, i.e., $\langle\mathcal{H}_{t}\rangle\sim|\mathcal{T}_{t}|^{1+\ln3/\ln5}$. Accordingly, we execute extensive experimental simulations, and demonstrate that empirical analysis is in strong agreement with theoretical results. Lastly, we provide a guide to extend the proposed iterative manner in order to generate more general scale-free tree networks with large diameter.


title: Deep Scattering Network with Max-pooling

abstract: Scattering network is a convolutional network, consisting of cascading convolutions using pre-defined wavelets followed by the modulus operator. Since its introduction in 2012, the scattering network is used as one of few mathematical tools explaining components of the convolutional neural networks (CNNs). However, a pooling operator, which is one of main components of conventional CNNs, is not considered in the original scattering network. In this paper, we propose a new network, called scattering-maxp network, integrating the scattering network with the max-pooling operator. We model continuous max-pooling, apply it to the scattering network, and get the scattering-maxp network. We show that the scattering-maxp network shares many useful properties of the scattering network including translation invariance, but with much smaller number of parameters. Numerical experiments show that the use of scattering-maxp does not degrade the performance too much and is much faster than the original one, in image classification tasks.


title: Mesh Total Generalized Variation for Denoising

abstract: Total Generalized Variation (TGV) has recently been proven certainly successful in image processing for preserving sharp features as well as smooth transition variations. However, none of the existing works aims at numerically calculating TGV over triangular meshes. In this paper, we develop a novel numerical framework to discretize the second-order TGV over triangular meshes. Further, we propose a TGV-based variational model to restore the face normal field for mesh denoising. The TGV regularization in the proposed model is represented by a combination of a first- and second-order term, which can be automatically balanced. This TGV regularization is able to locate sharp features and preserve them via the first-order term, while recognize smoothly curved regions and recover them via the second-order term. To solve the optimization problem, we introduce an efficient iterative algorithm based on variable-splitting and augmented Lagrangian method. Extensive results and comparisons on synthetic and real scanning data validate that the proposed method outperforms the state-of-the-art methods visually and numerically.


title: Diminishing Uncertainty within the Training Pool: Active Learning for Medical Image Segmentation

abstract: Active learning is a unique abstraction of machine learning techniques where the model/algorithm could guide users for annotation of a set of data points that would be beneficial to the model, unlike passive machine learning. The primary advantage being that active learning frameworks select data points that can accelerate the learning process of a model and can reduce the amount of data needed to achieve full accuracy as compared to a model trained on a randomly acquired data set. Multiple frameworks for active learning combined with deep learning have been proposed, and the majority of them are dedicated to classification tasks. Herein, we explore active learning for the task of segmentation of medical imaging data sets. We investigate our proposed framework using two datasets: 1.) MRI scans of the hippocampus, 2.) CT scans of pancreas and tumors. This work presents a query-by-committee approach for active learning where a joint optimizer is used for the committee. At the same time, we propose three new strategies for active learning: 1.) increasing frequency of uncertain data to bias the training data set; 2.) Using mutual information among the input images as a regularizer for acquisition to ensure diversity in the training dataset; 3.) adaptation of Dice log-likelihood for Stein variational gradient descent (SVGD). The results indicate an improvement in terms of data reduction by achieving full accuracy while only using 22.69 % and 48.85 % of the available data for each dataset, respectively.


title: Joint User Activity and Data Detection in Grant-Free NOMA using Generative Neural Networks

abstract: Grant-free non-orthogonal multiple access (NOMA) is considered as one of the supporting technology for massive connectivity for future networks. In the grant-free NOMA systems with a massive number of users, user activity detection is of great importance. Existing multi-user detection (MUD) techniques rely on complicated update steps which may cause latency in signal detection. In this paper, we propose a generative neural network-based MUD (GenMUD) framework to utilize low-complexity neural networks, which are trained to reconstruct signals in a small fixed number of steps. By exploiting the uncorrelated user behaviours, we design a network architecture to achieve higher recovery accuracy with a low computational cost. Experimental results show significant performance gains in detection accuracy compared to conventional solutions under different channel conditions and user sparsity levels. We also provide a sparsity estimator through extensive experiments. Simulation results of the sparsity estimator showed high estimation accuracy, strong robustness to channel variations and neglectable impact on support detection accuracy.


title: Understanding the Error in Evaluating Adversarial Robustness

abstract: Deep neural networks are easily misled by adversarial examples. Although lots of defense methods are proposed, many of them are demonstrated to lose effectiveness when against properly performed adaptive attacks. How to evaluate the adversarial robustness effectively is important for the realistic deployment of deep models, but yet still unclear. To provide a reasonable solution, one of the primary things is to understand the error (or gap) between the true adversarial robustness and the evaluated one, what is it and why it exists. Several works are done in this paper to make it clear. Firstly, we introduce an interesting phenomenon named gradient traps, which lead to incompetent adversaries and are demonstrated to be a manifestation of evaluation error. Then, we analyze the error and identify that there are three components. Each of them is caused by a specific compromise. Moreover, based on the above analysis, we present our evaluation suggestions. Experiments on adversarial training and its variations indicate that: (1) the error does exist empirically, and (2) these defenses are still vulnerable. We hope these analyses and results will help the community to develop more powerful defenses.


title: GraphHop: An Enhanced Label Propagation Method for Node Classification

abstract: A scalable semi-supervised node classification method on graph-structured data, called GraphHop, is proposed in this work. The graph contains attributes of all nodes but labels of a few nodes. The classical label propagation (LP) method and the emerging graph convolutional network (GCN) are two popular semi-supervised solutions to this problem. The LP method is not effective in modeling node attributes and labels jointly or facing a slow convergence rate on large-scale graphs. GraphHop is proposed to its shortcoming. With proper initial label vector embeddings, each iteration of GraphHop contains two steps: 1) label aggregation and 2) label update. In Step 1, each node aggregates its neighbors' label vectors obtained in the previous iteration. In Step 2, a new label vector is predicted for each node based on the label of the node itself and the aggregated label information obtained in Step 1. This iterative procedure exploits the neighborhood information and enables GraphHop to perform well in an extremely small label rate setting and scale well for very large graphs. Experimental results show that GraphHop outperforms state-of-the-art graph learning methods on a wide range of tasks (e.g., multi-label and multi-class classification on citation networks, social graphs, and commodity consumption graphs) in graphs of various sizes. Our codes are publicly available on GitHub (https://github.com/TianXieUSC/GraphHop).


title: A design of human-like robust AI machines in object identification

abstract: This is a perspective paper inspired from the study of Turing Test proposed by A.M. Turing (23 June 1912 - 7 June 1954) in 1950. Following one important implication of Turing Test for enabling a machine with a human-like behavior or performance, we define human-like robustness (HLR) for AI machines. The objective of the new definition aims to enforce AI machines with HLR, including to evaluate them in terms of HLR. A specific task is discussed only on object identification, because it is the most common task for every person in daily life. Similar to the perspective, or design, position by Turing, we provide a solution of how to achieve HLR AI machines without constructing them and conducting real experiments. The solution should consists of three important features in the machines. The first feature of HLR machines is to utilize common sense from humans for realizing a causal inference. The second feature is to make a decision from a semantic space for having interpretations to the decision. The third feature is to include a "human-in-the-loop" setting for advancing HLR machines. We show an "identification game" using proposed design of HLR machines. The present paper shows an attempt to learn and explore further from Turing Test towards the design of human-like AI machines.


title: Resilient quantum gates on periodically driven Rydberg atoms

abstract: Fault-tolerant implementation of quantum gates is one of preconditions for realizing quantum computation. The platform of Rydberg atoms is one of the most promising candidates for achieving quantum computation. We propose to implement a controlled-$Z$ gate on Rydberg atoms where an amplitude-modulated field is employed to induce Rydberg antiblockade. Gate robustness against the fluctuations in the Rydberg-Rydberg interaction can be largely enhanced by adjusting amplitude-modulated field. Furthermore, we introduce a Landau-Zener-St\"{u}ckelberg transition on the target atom so as to improve the gate resilience to the deviation in the gate time and the drift in the pulse amplitude. With feasible experimental parameters, one can achieve the gate with low fidelity errors caused by atomic decay, interatomic dipole-dipole force, and Doppler effects. Finally, we generalize the gate scheme into multiqubit cases, where resilient multiqubit phase gates can be obtained in one step with an unchanged gate time as the number of qubits increases.


title: Symmetry of Narayana numbers and rowvacuation of root posets

abstract: For a Weyl group $W$ of rank $r$, the $W$-Catalan number is the number of antichains of the poset of positive roots, and the $W$-Narayana numbers refine the $W$-Catalan number by keeping track of the cardinalities of these antichains. The $W$-Narayana numbers are symmetric, i.e., the number of antichains of cardinality $k$ is the same as the number of cardinality $r-k$. However, this symmetry is far from obvious. Panyushev posed the problem of defining an involution on root poset antichains that exhibits the symmetry of the $W$-Narayana numbers. Rowmotion and rowvacuation are two related operators, defined as compositions of "toggles," that give a dihedral action on the set of antichains of any ranked poset. Rowmotion acting on root posets has been the subject of a significant amount of research in the recent past. We prove that for the root posets of classical types, rowvacuation is Panyushev's desired involution.


title: Copula Quadrant Similarity for Anomaly Scores

abstract: Practical anomaly detection requires applying numerous approaches due to the inherent difficulty of unsupervised learning. Direct comparison between complex or opaque anomaly detection algorithms is intractable; we instead propose a framework for associating the scores of multiple methods. Our aim is to answer the question: how should one measure the similarity between anomaly scores generated by different methods? The scoring crux is the extremes, which identify the most anomalous observations. A pair of algorithms are defined here to be similar if they assign their highest scores to roughly the same small fraction of observations. To formalize this, we propose a measure based on extremal similarity in scoring distributions through a novel upper quadrant modeling approach, and contrast it with tail and other dependence measures. We illustrate our method with simulated and real experiments, applying spectral methods to cluster multiple anomaly detection methods and to contrast our similarity measure with others. We demonstrate that our method is able to detect the clusters of anomaly detection algorithms to achieve an accurate and robust ensemble algorithm.


title: Modeling and mitigation of cross-talk effects in readout noise with applications to the Quantum Approximate Optimization Algorithm

abstract: We introduce a correlated measurement noise model that can be efficiently described and characterized, and which admits effective noise-mitigation on the level of marginal probability distributions. Noise mitigation can be performed up to some error for which we derive upper bounds. Characterization of the model is done efficiently using Diagonal Detector Overlapping Tomography -- a generalization of the recently introduced Quantum Overlapping Tomography to the problem of reconstruction of readout noise with restricted locality. The procedure allows to characterize $k$-local measurement cross-talk on $N$-qubit device using $O(k2^klog(N))$ circuits containing random combinations of X and identity gates. We perform experiments on 15 (23) qubits using IBM's (Rigetti's) devices to test both the noise model and the error-mitigation scheme, and obtain an average reduction of errors by a factor $>22$ ($>5.5$) compared to no mitigation. Interestingly, we find that correlations in the measurement noise do not correspond to the physical layout of the device. Furthermore, we study numerically the effects of readout noise on the performance of the Quantum Approximate Optimization Algorithm (QAOA). We observe in simulations that for numerous objective Hamiltonians, including random MAX-2-SAT instances and the Sherrington-Kirkpatrick model, the noise-mitigation improves the quality of the optimization. Finally, we provide arguments why in the course of QAOA optimization the estimates of the local energy (or cost) terms often behave like uncorrelated variables, which greatly reduces sampling complexity of the energy estimation compared to the pessimistic error analysis. We also show that similar effects are expected for Haar-random quantum states and states generated by shallow-depth random circuits.


title: Identification of Latent Variables From Graphical Model Residuals

abstract: Graph-based causal discovery methods aim to capture conditional independencies consistent with the observed data and differentiate causal relationships from indirect or induced ones. Successful construction of graphical models of data depends on the assumption of causal sufficiency: that is, that all confounding variables are measured. When this assumption is not met, learned graphical structures may become arbitrarily incorrect and effects implied by such models may be wrongly attributed, carry the wrong magnitude, or mis-represent direction of correlation. Wide application of graphical models to increasingly less curated "big data" draws renewed attention to the unobserved confounder problem. We present a novel method that aims to control for the latent space when estimating a DAG by iteratively deriving proxies for the latent space from the residuals of the inferred model. Under mild assumptions, our method improves structural inference of Gaussian graphical models and enhances identifiability of the causal effect. In addition, when the model is being used to predict outcomes, it un-confounds the coefficients on the parents of the outcomes and leads to improved predictive performance when out-of-sample regime is very different from the training data. We show that any improvement of prediction of an outcome is intrinsically capped and cannot rise beyond a certain limit as compared to the confounded model. We extend our methodology beyond GGMs to ordinal variables and nonlinear cases. Our R package provides both PCA and autoencoder implementations of the methodology, suitable for GGMs with some guarantees and for better performance in general cases but without such guarantees.


title: Infinitely Wide Tensor Networks as Gaussian Process

abstract: Gaussian Process is a non-parametric prior which can be understood as a distribution on the function space intuitively. It is known that by introducing appropriate prior to the weights of the neural networks, Gaussian Process can be obtained by taking the infinite-width limit of the Bayesian neural networks from a Bayesian perspective. In this paper, we explore the infinitely wide Tensor Networks and show the equivalence of the infinitely wide Tensor Networks and the Gaussian Process. We study the pure Tensor Network and another two extended Tensor Network structures: Neural Kernel Tensor Network and Tensor Network hidden layer Neural Network and prove that each one will converge to the Gaussian Process as the width of each model goes to infinity. (We note here that Gaussian Process can also be obtained by taking the infinite limit of at least one of the bond dimensions $\alpha_{i}$ in the product of tensor nodes, and the proofs can be done with the same ideas in the proofs of the infinite-width cases.) We calculate the mean function (mean vector) and the covariance function (covariance matrix) of the finite dimensional distribution of the induced Gaussian Process by the infinite-width tensor network with a general set-up. We study the properties of the covariance function and derive the approximation of the covariance function when the integral in the expectation operator is intractable. In the numerical experiments, we implement the Gaussian Process corresponding to the infinite limit tensor networks and plot the sample paths of these models. We study the hyperparameters and plot the sample path families in the induced Gaussian Process by varying the standard deviations of the prior distributions. As expected, the parameters in the prior distribution namely the hyper-parameters in the induced Gaussian Process controls the characteristic lengthscales of the Gaussian Process.


title: Machine Learning on Cloud with Blockchain: A Secure, Verifiable and Fair Approach to Outsource the Linear Regression for Data Analysis

abstract: Linear Regression (LR) is a classical machine learning algorithm which has many applications in the cyber physical social systems (CPSS) to shape and simplify the way we live, work and communicate. This paper focuses on the data analysis for CPSS when the Linear Regression is applied. The training process of LR is time-consuming since it involves complex matrix operations, especially when it gets a large scale training dataset In the CPSS. Thus, how to enable devices to efficiently perform the training process of the Linear Regression is of significant importance. To address this issue, in this paper, we present a secure, verifiable and fair approach to outsource LR to an untrustworthy cloud-server. In the proposed scheme, computation inputs/outputs are obscured so that the privacy of sensitive information is protected against cloud-server. Meanwhile, computation result from cloud-server is verifiable. Also, fairness is guaranteed by the blockchain, which ensures that the cloud gets paid only if he correctly performed the outsourced workload. Based on the presented approach, we exploited the fair, secure outsourcing system on the Ethereum blockchain. We analysed our presented scheme on theoretical and experimental, all of which indicate that the presented scheme is valid, secure and efficient.


title: k-core percolation on interdependent and interconnected multiplex networks

abstract: Many real-world networks are coupled together to maintain their normal functions. Here we study the robustness of multiplex networks with interdependent and interconnected links under k-core percolation, where a node fails when it connects to a threshold of less than k neighbors. By deriving the self-consistency equations, we solve the key quantities of interests such as the critical threshold and size of the giant component analytically and validate the theoretical results with numerical simulations. We find a rich phase transition phenomenon as we tune the inter-layer coupling strength. Specifically speaking, in the ER-ER multiplex networks, with the increase of coupling strength, the size of the giant component in each layer first undergoes a first-order transition and then a second-order transition and finally a first-order transition. This is due to the nature of inter-layer links with both connectivity and dependency simultaneously. The system is more robust if the dependency on the initial robust network is strong and more vulnerable if the dependency on the initial attacked network is strong. These effects are even amplified in the cascading process. When applying our model to the SF-SF multiplex networks, the type of transition changes. The system undergoes a first-order phase transition first only when the two layers' mutually coupling is very strong and a second-order transition in other conditions.


title: Distributed Arithmetic Coding for Sources with Hidden Markov Correlation

abstract: Distributed arithmetic coding (DAC) has been shown to be effective for Slepian-Wolf coding, especially for short data blocks. In this letter, we propose to use the DAC to compress momery-correlated sources. More specifically, the correlation between sources is modeled as a hidden Markov process. Experimental results show that the performance is close to the theoretical Slepian-Wolf limit.


title: Learning Temporal Dynamics from Cycles in Narrated Video

abstract: Learning to model how the world changes as time elapses has proven a challenging problem for the computer vision community. We propose a self-supervised solution to this problem using temporal cycle consistency jointly in vision and language, training on narrated video. Our model learns modality-agnostic functions to predict forward and backward in time, which must undo each other when composed. This constraint leads to the discovery of high-level transitions between moments in time, since such transitions are easily inverted and shared across modalities. We justify the design of our model with an ablation study on different configurations of the cycle consistency problem. We then show qualitatively and quantitatively that our approach yields a meaningful, high-level model of the future and past. We apply the learned dynamics model without further training to various tasks, such as predicting future action and temporally ordering sets of images. Project page: https://dave.ml/mmcc


title: Max-Affine Spline Insights Into Deep Network Pruning

abstract: In this paper, we study the importance of pruning in Deep Networks (DNs) and the yin & yang relationship between (1) pruning highly overparametrized DNs that have been trained from random initialization and (2) training small DNs that have been "cleverly" initialized. As in most cases practitioners can only resort to random initialization, there is a strong need to develop a grounded understanding of DN pruning. Current literature remains largely empirical, lacking a theoretical understanding of how pruning affects DNs' decision boundary, how to interpret pruning, and how to design corresponding principled pruning techniques. To tackle those questions, we propose to employ recent advances in the theoretical analysis of Continuous Piecewise Affine (CPA) DNs. From this perspective, we will be able to detect the early-bird (EB) ticket phenomenon, provide interpretability into current pruning techniques, and develop a principled pruning strategy. In each step of our study, we conduct extensive experiments supporting our claims and results; while our main goal is to enhance the current understanding towards DN pruning instead of developing a new pruning method, our spline pruning criteria in terms of layerwise and global pruning is on par with or even outperforms state-of-the-art pruning methods.


title: Dyson's disordered linear chain from a random matrix theory viewpoint

abstract: The first work of Dyson relating to random matrix theory, "The dynamics of a disordered linear chain", is reviewed. Contained in this work is an exact solution of a so-called Type I chain in the case of the disorder variables being given by a gamma distribution. The exact solution exhibits a singularity in the density of states about the origin, which has since been shown to be universal for one-dimensional tight binding models with off diagonal disorder. We discuss this context and also point out some universal features of the weak disorder expansion of the exact solution near the band edge. Further, a link between the exact solution, and a tridiagonal formalism of anti symmetric Gaussian $\beta$-ensembles with $\beta$ proportional to $1/N$, is made.


title: Neutron Response of the EJ-254 Boron-Loaded Plastic Scintillator

abstract: Organic scintillators doped with capture agents provide a detectable signal for neutrons over a broad energy range. This work characterizes the fast and slow neutron response of EJ-254, an organic plastic scintillator with 5% natural boron loading by weight. For fast neutrons, the primary mechanism for light generation in organic scintillators is n-p elastic scattering. To study the fast neutron response, the proton light yield of EJ-254 was measured at the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory. Using a broad-spectrum neutron source and a double time-of-flight technique, the EJ-254 proton light yield was obtained over the energy range of approximately 270 keV to 4.5 MeV and determined to be in agreement with other plastic scintillators comprised of the same polymer base. To isolate the slow neutron response, an AmBe source with polyethylene moderator was made incident on the EJ-254 scintillator surrounded by an array of EJ-309 observation detectors. Events in the EJ-254 target coincident with the signature 477.6 keV $\gamma$ ray (resulting from deexcitation of the residual $^{7}$Li nucleus following boron neutron capture) were identified. Pulse shape discrimination was used to evaluate the temporal differences in the response of EJ-254 scintillation signals arising from $\gamma$-ray and fast/slow neutron interactions. Clear separation between $\gamma$-ray and fast neutrons signals was not achieved and the neutron capture feature was observed to overlap both the $\gamma$-ray and fast neutron bands. Taking into account the electron light nonproportionality, the neutron-capture light yield in EJ-254 was determined to be 89.4$\pm$1.1 keVee.


title: Blockchain Aided Privacy-Preserving Outsourcing Algorithms of Bilinear Pairings for Internet of Things Devices

abstract: Bilinear pairing is a fundamental operation that is widely used in cryptographic algorithms (e.g., identity-based cryptographic algorithms) to secure IoT applications. Nonetheless, the time complexity of bilinear pairing is $O(n^3)$, making it a very time-consuming operation, especially for resource-constrained IoT devices. Secure outsourcing of bilinear pairing has been studied in recent years to enable computationally weak devices to securely outsource the bilinear pairing to untrustworthy cloud servers. However, the state-of-art algorithms often require to pre-compute and store some values, which results in storage burden for devices. In the Internet of Things, devices are generally with very limited storage capacity. Thus, the existing algorithms do not fit the IoT well. In this paper, we propose a secure outsourcing algorithm of bilinear pairings, which does not require pre-computations. In the proposed algorithm, the outsourcer side's efficiency is significantly improved compared with executing the original bilinear pairing operation. At the same time, the privacy of the input and output is ensured. Also, we apply the Ethereum blockchain in our outsourcing algorithm to enable fair payments, which ensures that the cloud server gets paid only when he correctly accomplished the outsourced work. The theoretical analysis and experimental results show that the proposed algorithm is efficient and secure.


title: User Response Prediction in Online Advertising

abstract: Online advertising, as the vast market, has gained significant attention in various platforms ranging from search engines, third-party websites, social media, and mobile apps. The prosperity of online campaigns is a challenge in online marketing and is usually evaluated by user response through different metrics, such as clicks on advertisement (ad) creatives, subscriptions to products, purchases of items, or explicit user feedback through online surveys. Recent years have witnessed a significant increase in the number of studies using computational approaches, including machine learning methods, for user response prediction. However, existing literature mainly focuses on algorithmic-driven designs to solve specific challenges, and no comprehensive review exists to answer many important questions. What are the parties involved in the online digital advertising eco-systems? What type of data are available for user response prediction? How to predict user response in a reliable and/or transparent way? In this survey, we provide a comprehensive review of user response prediction in online advertising and related recommender applications. Our essential goal is to provide a thorough understanding of online advertising platforms, stakeholders, data availability, and typical ways of user response prediction. We propose a taxonomy to categorize state-of-the-art user response prediction methods, primarily focus on the current progress of machine learning methods used in different online platforms. In addition, we also review applications of user response prediction, benchmark datasets, and open-source codes in the field.


title: Distributing hash families with few rows

abstract: Column replacement techniques for creating covering arrays rely on the construction of perfect and distributing hash families with few rows, having as many columns as possible for a specified number of symbols. To construct distributing hash families in which the number of rows is less than the strength, we examine a method due to Blackburn and extend it in three ways. First, the method is generalized from homogeneous hash families (in which every row has the same number of symbols) to heterogeneous ones. Second, the extension treats distributing hash families, in which only separation into a prescribed number of parts is required, rather than perfect hash families, in which columns must be completely separated. Third, the requirements on one of the main ingredients are relaxed to permit the use of a large class of distributing hash families, which we call fractal. Constructions for fractal perfect and distributing hash families are given, and applications to the construction of perfect hash families of large strength are developed.


title: DICE: Deep Significance Clustering for Outcome-Aware Stratification

abstract: We present deep significance clustering (DICE), a framework for jointly performing representation learning and clustering for "outcome-aware" stratification. DICE is intended to generate cluster membership that may be used to categorize a population by individual risk level for a targeted outcome. Following the representation learning and clustering steps, we embed the objective function in DICE with a constraint which requires a statistically significant association between the outcome and cluster membership of learned representations. DICE further includes a neural architecture search step to maximize both the likelihood of representation learning and outcome classification accuracy with cluster membership as the predictor. To demonstrate its utility in medicine for patient risk-stratification, the performance of DICE was evaluated using two datasets with different outcome ratios extracted from real-world electronic health records. Outcomes are defined as acute kidney injury (30.4\%) among a cohort of COVID-19 patients, and discharge disposition (36.8\%) among a cohort of heart failure patients, respectively. Extensive results demonstrate that DICE has superior performance as measured by the difference in outcome distribution across clusters, Silhouette score, Calinski-Harabasz index, and Davies-Bouldin index for clustering, and Area under the ROC Curve (AUC) for outcome classification compared to several baseline approaches.


title: A generalization of the Von Neumann extractor

abstract: An iterative randomness extraction algorithm which generalized the Von Neumann's extraction algorithm is detailed, analyzed and implemented in standard C++. Given a sequence of independently and identically distributed biased Bernoulli random variables, to extract randomness from the aforementioned sequence pertains to produce a new sequence of independently and identically distributed unbiased Bernoulli random variables. The iterative construction here is inspired from the work of Stout and Warren 1984 who modified appropriately the tree of probabilities produced by recursively repeating the Von Neumann's extraction algorithm. The correctness of the iterative algorithm is proven. The number of biased Bernoulli random variables needed to produce one unbiased instance is the complexity of interest. The complexity depends on the bias of the source. The expected complexity converges toward 3.10220648... when the bias tends to 0 and diverges when the bias tends to 1/2. In addition to the expected complexity, some other results that concern the limiting asymptotic construction, and that seem unnoticed in the literature so far, are proven.


title: Multitask Learning for Emotion and Personality Detection

abstract: In recent years, deep learning-based automated personality trait detection has received a lot of attention, especially now, due to the massive digital footprints of an individual. Moreover, many researchers have demonstrated that there is a strong link between personality traits and emotions. In this paper, we build on the known correlation between personality traits and emotional behaviors, and propose a novel multitask learning framework, SoGMTL that simultaneously predicts both of them. We also empirically evaluate and discuss different information-sharing mechanisms between the two tasks. To ensure the high quality of the learning process, we adopt a MAML-like framework for model optimization. Our more computationally efficient CNN-based multitask model achieves the state-of-the-art performance across multiple famous personality and emotion datasets, even outperforming Language Model based models.


title: SDP Achieves Exact Minimax Optimality in Phase Synchronization

abstract: We study the phase synchronization problem with noisy measurements $Y=z^*z^{*H}+\sigma W\in\mathbb{C}^{n\times n}$, where $z^*$ is an $n$-dimensional complex unit-modulus vector and $W$ is a complex-valued Gaussian random matrix. It is assumed that each entry $Y_{jk}$ is observed with probability $p$. We prove that an SDP relaxation of the MLE achieves the error bound $(1+o(1))\frac{\sigma^2}{2np}$ under a normalized squared $\ell_2$ loss. This result matches the minimax lower bound of the problem, and even the leading constant is sharp. The analysis of the SDP is based on an equivalent non-convex programming whose solution can be characterized as a fixed point of the generalized power iteration lifted to a higher dimensional space. This viewpoint unifies the proofs of the statistical optimality of three different methods: MLE, SDP, and generalized power method. The technique is also applied to the analysis of the SDP for $\mathbb{Z}_2$ synchronization, and we achieve the minimax optimal error $\exp\left(-(1-o(1))\frac{np}{2\sigma^2}\right)$ with a sharp constant in the exponent.


title: Simulations of globular clusters within their parent galaxies: multiple stellar populations and internal kinematics

abstract: Using three-dimensional smoothed particle hydrodynamics simulations, we investigate the formation of multiple stellar populations (MSPs) in globular clusters (GCs) within the context of their parent galaxies. In our scenario, the second generation (2G) of stars originate from both asymptotic giant branch (AGB) polluters and pristine gas accreted from the host galaxy. Previous theoretical and numerical studies have demonstrated that this 'AGB with dilution' model has the potential to alleviate several problems faced by the classical AGB scenario. However, the accretion of pristine gas on to the GC has yet to be investigated within the context of the parent galaxy. This paper presents the preliminary results from our original simulation code which models GC formation from giant molecular clouds in a host galaxy, and subsequent gas accretion on to the GC. By simulating the genesis of the 2G over a 370 Myr time frame, we demonstrate that the fraction of 2G stars are inextricably linked to the GC's environment. Our simulations rationalize the wide variety of abundance patterns, kinematics, and 2G concentrations by altering the initial conditions of both the GC progenitor and the host galaxy itself. Most notably, we reproduce a positive correlation between the fraction of 2G stars and the initial mass of the cluster. We discuss the physical implications of our scenario and compare our simulations with observations of the Galactic GC 47 Tucanae (47 Tuc). Finally, we present scaling relations that encompass the wider GC population and serve as a reference for future observations.


title: Attention Actor-Critic algorithm for Multi-Agent Constrained Co-operative Reinforcement Learning

abstract: In this work, we consider the problem of computing optimal actions for Reinforcement Learning (RL) agents in a co-operative setting, where the objective is to optimize a common goal. However, in many real-life applications, in addition to optimizing the goal, the agents are required to satisfy certain constraints specified on their actions. Under this setting, the objective of the agents is to not only learn the actions that optimize the common objective but also meet the specified constraints. In recent times, the Actor-Critic algorithm with an attention mechanism has been successfully applied to obtain optimal actions for RL agents in multi-agent environments. In this work, we extend this algorithm to the constrained multi-agent RL setting. The idea here is that optimizing the common goal and satisfying the constraints may require different modes of attention. By incorporating different attention modes, the agents can select useful information required for optimizing the objective and satisfying the constraints separately, thereby yielding better actions. Through experiments on benchmark multi-agent environments, we show the effectiveness of our proposed algorithm.


title: Polynomial Duality-Symmetric Lagrangians for Free p-Forms

abstract: We explore the properties of polynomial Lagrangians for chiral $p$-forms previously proposed by the last named author, and in particular, provide a self-contained treatment of the symmetries and equations of motion that shows a great economy and simplicity of this formalism. We further use analogous techniques to construct polynomial democratic Lagrangians for general $p$-forms where electric and magnetic potentials appear on equal footing as explicit dynamical variables. Due to our reliance on the differential form notation, the construction is compact and universally valid for forms of all ranks, in any number of dimensions.


title: Applying Transfer Learning for Improving Domain-Specific Search Experience Using Query to Question Similarity

abstract: Search is one of the most common platforms used to seek information. However, users mostly get overloaded with results whenever they use such a platform to resolve their queries. Nowadays, direct answers to queries are being provided as a part of the search experience. The question-answer (QA) retrieval process plays a significant role in enriching the search experience. Most off-the-shelf Semantic Textual Similarity models work fine for well-formed search queries, but their performances degrade when applied to a domain-specific setting having incomplete or grammatically ill-formed search queries in prevalence. In this paper, we discuss a framework for calculating similarities between a given input query and a set of predefined questions to retrieve the question which matches to it the most. We have used it for the financial domain, but the framework is generalized for any domain-specific search engine and can be used in other domains as well. We use Siamese network [6] over Long Short-Term Memory (LSTM) [3] models to train a classifier which generates unnormalized and normalized similarity scores for a given pair of questions. Moreover, for each of these question pairs, we calculate three other similarity scores: cosine similarity between their average word2vec embeddings [15], cosine similarity between their sentence embeddings [7] generated using RoBERTa [17] and their customized fuzzy-match score. Finally, we develop a metaclassifier using Support Vector Machines [19] for combining these five scores to detect if a given pair of questions is similar. We benchmark our model's performance against existing State Of The Art (SOTA) models on Quora Question Pairs (QQP) dataset as well as a dataset specific to the financial domain.


title: M\"{o}biusE: Knowledge Graph Embedding on M\"{o}bius Ring

abstract: In this work, we propose a novel Knowledge Graph Embedding (KGE) strategy, called M\"{o}biusE, in which the entities and relations are embedded to the surface of a M\"{o}bius ring. The proposition of such a strategy is inspired by the classic TorusE, in which the addition of two arbitrary elements is subject to a modulus operation. In this sense, TorusE naturally guarantees the critical boundedness of embedding vectors in KGE. However, the nonlinear property of addition operation on Torus ring is uniquely derived by the modulus operation, which in some extent restricts the expressiveness of TorusE. As a further generalization of TorusE, M\"{o}biusE also uses modulus operation to preserve the closeness of addition operation on it, but the coordinates on M\"{o}bius ring interacts with each other in the following way: {\em \color{red} any vector on the surface of a M\"{o}bius ring moves along its parametric trace will goes to the right opposite direction after a cycle}. Hence, M\"{o}biusE assumes much more nonlinear representativeness than that of TorusE, and in turn it generates much more precise embedding results. In our experiments, M\"{o}biusE outperforms TorusE and other classic embedding strategies in several key indicators.


title: Low-cost and high-performance data augmentation for deep-learning-based skin lesion classification

abstract: Although deep convolutional neural networks (DCNNs) have achieved significant accuracy in skin lesion classification comparable or even superior to those of dermatologists, practical implementation of these models for skin cancer screening in low resource settings is hindered by their limitations in computational cost and training dataset. To overcome these limitations, we propose a low-cost and high-performance data augmentation strategy that includes two consecutive stages of augmentation search and network search. At the augmentation search stage, the augmentation strategy is optimized in the search space of Low-Cost-Augment (LCA) under the criteria of balanced accuracy (BACC) with 5-fold cross validation. At the network search stage, the DCNNs are fine-tuned with the full training set in order to select the model with the highest BACC. The efficiency of the proposed data augmentation strategy is verified on the HAM10000 dataset using EfficientNets as a baseline. With the proposed strategy, we are able to reduce the search space to 60 and achieve a high BACC of 0.853 by using a single DCNN model without external database, suitable to be implemented in mobile devices for DCNN-based skin lesion detection in low resource settings.


title: Kullback-Leibler-Based Discrete Failure Time Models for Integration of Published Prediction Models with New Time-To-Event Dataset

abstract: Prediction of time-to-event data often suffers from rare event rates, small sample sizes, high dimensionality and low signal-to-noise ratios. Incorporating published prediction models from large-scale studies is expected to improve the performance of prognosis prediction on internal individual-level time-to-event data. However, existing integration approaches typically assume that underlying distributions from the external and internal data sources are similar, which is often invalid. To account for challenges including heterogeneity, data sharing, and privacy constraints, we propose a discrete failure time modeling procedure, which utilizes a discrete hazard-based Kullback-Leibler discriminatory information measuring the discrepancy between the published models and the internal dataset. Simulations show the advantage of the proposed method compared with those solely based on the internal data or published models. We apply the proposed method to improve prediction performance on a kidney transplant dataset from a local hospital by integrating this small-scale dataset with published survival models obtained from the national transplant registry.


title: Fractional Buffer Layers: Absorbing Boundary Conditions for Wave Propagation

abstract: We develop fractional buffer layers (FBLs) to absorb propagating waves without reflection in bounded domains. Our formulation is based on variable-order spatial fractional derivatives. We select a proper variable-order function so that dissipation is induced to absorb the coming waves in the buffer layers attached to the domain. In particular, we first design proper FBsL for the one-dimensional one-way and two-way wave propagation. Then, we extend our formulation to two-dimensional problems, where we introduce a consistent variable-order fractional wave equation. In each case, we obtain the fully discretized equations by employing a spectral collocation method in space and Crank-Nicolson or Adams-Bashforth method in time. We compare our results with the perfectly matched layer (PML) method and show the effectiveness of FBL in accurately suppressing any erroneously reflected waves, including corner reflections in two-dimensional rectangular domains. FBLs can be used in conjunction with any discretization method appropriate for fractional operators describing wave propagation in bounded or truncated domains.


title: Hydrodynamic description of $D$ meson production in high-energy heavy-ion collisions

abstract: The large values and the constituent-quark-number (NCQ) scaling of the elliptic flow of low-$p_T$ $D$ mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through the strong interactions with the quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-$p_T$ charm quarks, we compare the $D^0$ meson spectra and elliptic flow from a hydrodynamic model to the experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in local comoving frame of the medium before freeze-out. $D^0$ mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both hydrodynamic and transport model can describe the elliptic flow of $D^0$ mesons at $p_T<3$ GeV/$c$ as measured in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV. Though the experimental data on $D^0$ spectra are consistent with the hydrodynamic result at small $p_T\sim 1$ GeV/$c$, they deviate from the hydrodynamic model at high transverse momentum $p_T>2$ GeV/$c$. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Our comparative study indicates that charm quarks only approach to local thermal equilibrium at small $p_T$ even though they acquire sizable elliptic flow that is comparable to light-quark hadrons at both small and intermediate $p_T$.


title: Rapid ortho-to-para nuclear spin conversion of H2 on a silicate dust surface

abstract: The H2 molecule has two nuclear spin isomers, the so-called ortho and para isomers. Nuclear spin conversion (NSC) between these states is forbidden in the gas phase. The energy difference between the lowest ortho and para states is as large as 14.7 meV, corresponding to ~170 K. Therefore, each state of H2 differently affects not only the chemistry but also the macroscopic gas dynamics in space, and thus, the ortho-to-para abundance ratio (OPR) of H2 has significant impacts on various astronomical phenomena. For a long time, the OPR of nascent H2 upon formation on dust grains has been assumed to have a statistical value of three and to gradually equilibrate in the gas phase at the temperature of the circumstances. Recently, NSC of H2 was experimentally revealed to occur on water ice at very low temperatures and thus incorporated into gas-dust chemical models. However, H2 molecules should form well before dust grains are coated by water ice. Information about how the OPR of H2 behaves on bare silicate dust before ice-mantle formation is lacking. Knowing the influence of the OPR of H2 if the OPR changes even on a bare silicate surface within an astronomically meaningful time scale is desirable. We report the first laboratory measurements of NSC of H2 physisorbed on amorphous silicate (Mg2SiO4) at temperatures up to 18 K. The conversion was found to occur very rapidly.