title: Message Passing for Hyper-Relational Knowledge Graphs

abstract: Hyper-relational knowledge graphs (KGs) (e.g., Wikidata) enable associating additional key-value pairs along with the main triple to disambiguate, or restrict the validity of a fact. In this work, we propose a message passing based graph encoder - StarE capable of modeling such hyper-relational KGs. Unlike existing approaches, StarE can encode an arbitrary number of additional information (qualifiers) along with the main triple while keeping the semantic roles of qualifiers and triples intact. We also demonstrate that existing benchmarks for evaluating link prediction (LP) performance on hyper-relational KGs suffer from fundamental flaws and thus develop a new Wikidata-based dataset - WD50K. Our experiments demonstrate that StarE based LP model outperforms existing approaches across multiple benchmarks. We also confirm that leveraging qualifiers is vital for link prediction with gains up to 25 MRR points compared to triple-based representations.


title: Repulsive properties of hadrons in lattice QCD data and neutron stars

abstract: Second-order susceptibilities $\chi^{11}_{ij}$ of baryon, electric, and strangeness, $B$, $Q$, and $S$, charges, are calculated in the Chiral Mean Field (CMF) model and compared to available lattice QCD data. The susceptibilities are sensitive to the short range repulsive interactions between different hadron species, especially to the hardcore repulsion of hyperons. Decreasing the hyperons size, as compared to the size of the non-strange baryons, does improve significantly the agreement of the CMF model results with the Lattice QCD data. The electric charge-dependent susceptibilities are sensitive to the short range repulsive volume of mesons. The comparison with lattice QCD data suggests that strange baryons, non-strange mesons and strange mesons have significantly smaller excluded volumes than non-strange baryons. The CMF model with these modified hadron volumes allows for a mainly hadronic description of the QCD susceptibilities significantly above the chiral pseudo-critical temperature. This improved CMF model which is based on the lattice QCD data, has been used to study the properties of both cold QCD matter and neutron star matter. The phase structure in both cases is essentially unchanged, i.e. a chiral first-order phase transition occurs at low temperatures ($T_{\rm CP}\approx 17$ MeV), and hyperons survive deconfinement to higher densities than non-strange hadrons. The neutron star maximal mass remains close to 2.1$M_\odot$ and the mass-radius diagram is only modified slightly due to the appearance of hyperons and is in agreement with astrophysical observations.


title: Variational and parquet-diagram calculations for neutron matter. II. Twisted Chain Diagrams

abstract: We develop a manifestly microscopic method to deal with strongly interacting nuclear systems that have different interactions in spin-singlet and spin-triplet states. In a first step we analyze variational wave functions that have been suggested to describe such systems, and demonstrate that the so-called commutator contributions can have important effects whenever the interactions in the spin-singlet and the spin-triplet states are very different. We then identify these contributions as terms that correspond, in the language of perturbation theory, to non-parquet diagrams. We include these diagrams in a way that is suggested by the Jastrow-Feenberg approach and show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects.


title: An analogue of $k$-marked Durfee symbols for strongly unimodal sequences

abstract: In a seminal 2007 paper, Andrews introduced a class of combinatorial objects that generalize partitions called $k$-marked Durfee symbols. Multivariate rank generating functions for these objects have been shown by many to have interesting modularity properties at certain vectors of roots of unity. Motivated by recent studies of rank generating functions for strongly unimodal sequences, we apply methods of Andrews to define an analogous class of combinatorial objects called $k$-marked strongly unimodal symbols that generalize strongly unimodal sequences. We establish a multivariate rank generating function for these objects, which we study combinatorially. We conclude by discussing potential quantum modularity properties for this rank generating function at certain vectors of roots of unity.


title: Permutation binomials of the form $x^r(x^{q-1}+a)$ over $\mathbb F_{q^e}$

abstract: We present several existence and nonexistence results for permutation binomials of the form $x^r(x^{q-1}+a)$, where $e\geq 2$ and $a\in \mathbb{F}_{q^e}^*$. As a consequence, we obtain a complete characterization of such permutation binomials over $\mathbb{F}_{q^2}$, $\mathbb{F}_{q^3}$, $\mathbb{F}_{q^4}$, $\mathbb{F}_{p^5}$, and $\mathbb{F}_{p^6}$, where $p$ is an odd prime.


title: Efficient Portfolios

abstract: Given two random realized returns on an investment, which is to be preferred? This is a fundamental problem in finance that has no definitive solution except in the case one investment always returns more than the other. In 1952 Markowitz and Roy introduced the following criterion for risk vs. return in portfolio selection: if two portfolios have the same expected realized return then prefer the one with smaller variance. An efficient portfolio has the least variance among all portfolios having the same expected realized return. The primary contribution of this short note is observation that the CAPM formula holds for realized returns as random variables, not just their expectations. This follows directly from writing down a mathematical model for one period investments.


title: Long Distance Continuous Methane Emissions Monitoring with Dual Frequency Comb Spectroscopy: deployment and blind testing in complex emissions scenarios

abstract: Continuous monitoring of oil and gas infrastructure is of interest for improving emissions and safety by enabling rapid identification and repair of emission sources, especially large sources that are responsible for the bulk of total emissions. We have previously demonstrated dual frequency comb spectroscopy (DCS) coupled with atmospheric modeling and inversion techniques (the DCS Observing System) as a viable and accurate approach for detection, attribution and quantification of methane emissions at distances of more than 1 km under controlled, steady emissions scenarios. Here, we present the results of validation testing designed to mimic the complexity of operational well pad emissions from oil and gas production, and the first field measurements at an active oil and gas facility. The validation tests are performed single-blind (the measurement and data analysis team are not given information about the emissions) at the Methane Emissions Technology Evaluation Center (METEC) test facility. They consist of a series of scenarios ranging from a single, steady-rate emission point to multiple emission points that include intermittent releases (the METEC "R2" tests). Additionally, we present field measurements at an active natural gas storage facility demonstrating that the system can remotely and autonomously monitor methane emissions in a true industrial setting. This field verification is in a configuration designed for continuous and long-term characterization of operational and fugitive emissions. These demonstrations confirm that the DCS Observing System can provide high-confidence continuous monitoring of emissions from complex, operational facilities among natural gas infrastructure.


title: Linear Systematics Mitigation in Galaxy Clustering in the Dark Energy Survey Year 1 Data

abstract: We implement a linear model for mitigating the effect of observing conditions and other sources of contamination in galaxy clustering analyses. Our treatment improves upon the fiducial systematics treatment of the Dark Energy Survey (DES) Year 1 (Y1) cosmology analysis in four crucial ways. Specifically, our treatment: 1) does not require decisions as to which observable systematics are significant and which are not, allowing for the possibility of multiple maps adding coherently to give rise to significant bias even if no single map leads to a significant bias by itself; 2) characterizes both the statistical and systematic uncertainty in our mitigation procedure, allowing us to propagate said uncertainties into the reported cosmological constraints; 3) explicitly exploits the full spatial structure of the galaxy density field to differentiate between cosmology-sourced and systematics-sourced fluctuations within the galaxy density field; 4) is fully automated, and can therefore be trivially applied to any data set. The updated correlation function for the DES Y1 redMaGiC catalog minimally impacts the cosmological posteriors from that analysis. Encouragingly, our analysis does improve the goodness of fit statistic of the DES Y1 3$\times$2pt data set ($\Delta \chi^2 = -6.5$ with no additional parameters). This improvement is due in nearly equal parts to both the change in the correlation function and the added statistical and systematic uncertainties associated with our method. We expect the difference in mitigation techniques to become more important in future work as the size of cosmological data sets grows.


title: Controlling Style in Generated Dialogue

abstract: Open-domain conversation models have become good at generating natural-sounding dialogue, using very large architectures with billions of trainable parameters. The vast training data required to train these architectures aggregates many different styles, tones, and qualities. Using that data to train a single model makes it difficult to use the model as a consistent conversational agent, e.g. with a stable set of persona traits and a typical style of expression. Several architectures affording control mechanisms over generation architectures have been proposed, each with different trade-offs. However, it remains unclear whether their use in dialogue is viable, and what the trade-offs look like with the most recent state-of-the-art conversational architectures. In this work, we adapt three previously proposed controllable generation architectures to open-domain dialogue generation, controlling the style of the generation to match one among about 200 possible styles. We compare their respective performance and tradeoffs, and show how they can be used to provide insights into existing conversational datasets, and generate a varied set of styled conversation replies.


title: Effects of Non-locality in Gravity and Quantum Theory

abstract: Spacetime---the union of space and time---is both the actor and the stage during physical processes in our fascinating Universe. In Lorentz invariant local theories, the existence of a maximum signalling speed (the "speed of light") determines a notion of causality in spacetime, distinguishing the past from the future, and the cause from the effect. This thesis is dedicated to the study of \emph{deviations} from locality. Focussing on a particular class of \emph{non-local} theories that is both Lorentz invariant and free of ghosts, we aim to understand the effects of such non-local physics in both gravity and quantum theory. Non-local ghost-free theories are accompanied by a parameter $\ell$ of dimension length that parametrizes the scale of non-locality, and for that reason we strive to express all effects of non-locality in terms of this symbol. In the limiting case of $\ell=0$ one recovers the local theory, and the effects of non-locality vanish. In order to address these questions we develop the notion of non-local Green functions [...]. The results presented in this thesis establish several effects of a Lorentz invariant, ghost-free non-locality in the areas of both gravitational and quantum physics. (Full abstract in document.)


title: A bound for the exterior product of $S$-units

abstract: We generalize an inequality for the determinant of a real matrix proved by A. Schinzel, to more general exterior products of vectors in Euclidean space. We apply this inequality to the logarithmic embedding of $S$-units contained in a number field $k$. This leads to a bound for the exterior product of $S$-units expressed as a product of heights. Using a volume formula of P. McMullen we show that our inequality is sharp up to a constant that depends only on the rank of the $S$-unit group but not on the field $k$. Our inequality is related to a conjecture of F. Rodriguez Villegas.


title: Improving Medical Annotation Quality to Decrease Labeling Burden Using Stratified Noisy Cross-Validation

abstract: As machine learning has become increasingly applied to medical imaging data, noise in training labels has emerged as an important challenge. Variability in diagnosis of medical images is well established; in addition, variability in training and attention to task among medical labelers may exacerbate this issue. Methods for identifying and mitigating the impact of low quality labels have been studied, but are not well characterized in medical imaging tasks. For instance, Noisy Cross-Validation splits the training data into halves, and has been shown to identify low-quality labels in computer vision tasks; but it has not been applied to medical imaging tasks specifically. In this work we introduce Stratified Noisy Cross-Validation (SNCV), an extension of noisy cross validation. SNCV can provide estimates of confidence in model predictions by assigning a quality score to each example; stratify labels to handle class imbalance; and identify likely low-quality labels to analyze the causes. We assess performance of SNCV on diagnosis of glaucoma suspect risk from retinal fundus photographs, a clinically important yet nuanced labeling task. Using training data from a previously-published deep learning model, we compute a continuous quality score (QS) for each training example. We relabel 1,277 low-QS examples using a trained glaucoma specialist; the new labels agree with the SNCV prediction over the initial label >85% of the time, indicating that low-QS examples mostly reflect labeler errors. We then quantify the impact of training with only high-QS labels, showing that strong model performance may be obtained with many fewer examples. By applying the method to randomly sub-sampled training dataset, we show that our method can reduce labelling burden by approximately 50% while achieving model performance non-inferior to using the full dataset on multiple held-out test sets.


title: A Macaulay2 Package for Stanley Simplicial Poset Ideals

abstract: We give a description of a new Macaulay2 package called SimplicialPosets. This package provides functions for working with simplicial posets and calculating their generalized Stanley-Reisner ideals. For practical purposes, we also introduce of a new random model for a class of simplicial posets which generalizes existing models for random simplicial complexes such as the Kahle model.


title: Bifurcation Analysis of Systems with Delays: Methods and Their Use in Applications

abstract: This chapter presents a dynamical systems point of view of the study of systems with delays. The focus is on how advanced tools from bifurcation theory, as implemented for example in the package DDE-BIFTOOL, can be applied to the study of delay differential equations (DDEs) arising in applications, including those that feature state-dependent delays. We discuss the present capabilities of the most recent release of DDE-BIFTOOL. They include the numerical continuation of steady states, periodic orbits and their bifurcations of codimension one, as well as the detection of certain bifurcations of codimension two and the calculation of their normal forms. Two longer case studies, of a conceptual DDE model for the El Nino phenomenon and of a prototypical scalar DDE with two state-dependent feedback terms, demonstrate what kind of insights can be obtained in this way.


title: Distributed Differentially Private Mutual Information Ranking and Its Applications

abstract: Computation of Mutual Information (MI) helps understand the amount of information shared between a pair of random variables. Automated feature selection techniques based on MI ranking are regularly used to extract information from sensitive datasets exceeding petabytes in size, over millions of features and classes. Series of one-vs-all MI computations can be cascaded to produce n-fold MI results, rapidly pinpointing informative relationships. This ability to quickly pinpoint the most informative relationships from datasets of billions of users creates privacy concerns. In this paper, we present Distributed Differentially Private Mutual Information (DDP-MI), a privacy-safe fast batch MI, across various scenarios such as feature selection, segmentation, ranking, and query expansion. This distributed implementation is protected with global model differential privacy to provide strong assurances against a wide range of privacy attacks. We also show that our DDP-MI can substantially improve the efficiency of MI calculations compared to standard implementations on a large-scale public dataset.


title: An Intuitive Tutorial to Gaussian Process Regression

abstract: This tutorial aims to provide an intuitive introduction to Gaussian process regression (GPR). GPR models have been widely used in machine learning applications due to their representation flexibility and inherent capability to quantify uncertainty over predictions. The tutorial starts with explaining the basic concepts that a Gaussian process is built on, including multivariate normal distribution, kernels, non-parametric models, and joint and conditional probability. It then provides a concise description of GPR and an implementation of a standard GPR algorithm. In addition, the tutorial reviews packages for implementing state-of-the-art Gaussian process algorithms. This tutorial is accessible to a broad audience, including those new to machine learning, ensuring a clear understanding of GPR fundamentals.


title: Initial Guesses for Sequences of Linear Systems in a GPU-Accelerated Incompressible Flow Solver

abstract: We consider several methods for generating initial guesses when iteratively solving sequences of linear systems, showing that they can be implemented efficiently in GPU-accelerated PDE solvers, specifically solvers for incompressible flow. We propose new initial guess methods based on stabilized polynomial extrapolation and compare them to the projection method of Fischer [15], showing that they are generally competitive with projection schemes despite requiring only half the storage and performing considerably less data movement and communication. Our implementations of these algorithms are freely available as part of the libParanumal collection of GPU-accelerated flow solvers.


title: Behavioral Repertoires for Soft Tensegrity Robots

abstract: Mobile soft robots offer compelling applications in fields ranging from urban search and rescue to planetary exploration. A critical challenge of soft robotic control is that the nonlinear dynamics imposed by soft materials often result in complex behaviors that are counterintuitive and hard to model or predict. As a consequence, most behaviors for mobile soft robots are discovered through empirical trial and error and hand-tuning. A second challenge is that soft materials are difficult to simulate with high fidelity -- leading to a significant reality gap when trying to discover or optimize new behaviors. In this work we employ a Quality Diversity Algorithm running model-free on a physical soft tensegrity robot that autonomously generates a behavioral repertoire with no a priori knowledge of the robot dynamics, and minimal human intervention. The resulting behavior repertoire displays a diversity of unique locomotive gaits useful for a variety of tasks. These results help provide a road map for increasing the behavioral capabilities of mobile soft robots through real-world automation.


title: The Physics Inventory of Quantitative Literacy: A tool for assessing mathematical reasoning in introductory physics

abstract: One desired outcome of introductory physics instruction is that students will develop facility with reasoning quantitatively about physical phenomena. Little research has been done regarding how students develop the algebraic concepts and skills involved in reasoning productively about physics quantities, which is different from either understanding of physics concepts or problem-solving abilities. We introduce the Physics Inventory of Quantitative Literacy (PIQL) as a tool for measuring quantitative literacy, a foundation of mathematical reasoning, in the context of introductory physics. We present the development of the PIQL and evidence of its validity for use in calculus-based introductory physics courses. Unlike concept inventories, the PIQL is a reasoning inventory, and can be used to assess reasoning over the span of students' instruction in introductory physics. Although mathematical reasoning associated with the PIQL is taught in prior mathematics courses, pre/post test scores reveal that this reasoning isn't readily used by most students in physics, nor does it develop as part of physics instruction--even in courses that use high-quality, research-based curricular materials. As has been the case with many inventories in physics education, we expect use of the PIQL to support the development of instructional strategies and materials--in this case, designed to meet the course objective that all students become quantitatively literate in introductory physics.


title: Topological Strings on Toric geometries in the presence of Lagrangian branes

abstract: We propose expressions for refined open topological string partition function on certain non-compact Calabi Yau 3-folds with topological branes wrapped on the special lagrangian submanifolds. The corresponding web diagrams are partially compact and a lagrangian brane is inserted on one of the external legs. Partial compactification introduces a mass deformation in the corresponding gauge theory. We propose conjectures that equate these open topological string partition functions with the generating function of equivaraint indices on certain quiver moduli spaces. To obtain these conjectures we use the identification of topological string partition functions with equivariant indices on the instanton moduli spaces.


title: Online AUC Optimization for Sparse High-Dimensional Datasets

abstract: The Area Under the ROC Curve (AUC) is a widely used performance measure for imbalanced classification arising from many application domains where high-dimensional sparse data is abundant. In such cases, each $d$ dimensional sample has only $k$ non-zero features with $k \ll d$, and data arrives sequentially in a streaming form. Current online AUC optimization algorithms have high per-iteration cost $\mathcal{O}(d)$ and usually produce non-sparse solutions in general, and hence are not suitable for handling the data challenge mentioned above. In this paper, we aim to directly optimize the AUC score for high-dimensional sparse datasets under online learning setting and propose a new algorithm, \textsc{FTRL-AUC}. Our proposed algorithm can process data in an online fashion with a much cheaper per-iteration cost $\mathcal{O}(k)$, making it amenable for high-dimensional sparse streaming data analysis. Our new algorithmic design critically depends on a novel reformulation of the U-statistics AUC objective function as the empirical saddle point reformulation, and the innovative introduction of the "lazy update" rule so that the per-iteration complexity is dramatically reduced from $\mathcal{O}(d)$ to $\mathcal{O}(k)$. Furthermore, \textsc{FTRL-AUC} can inherently capture sparsity more effectively by applying a generalized Follow-The-Regularized-Leader (FTRL) framework. Experiments on real-world datasets demonstrate that \textsc{FTRL-AUC} significantly improves both run time and model sparsity while achieving competitive AUC scores compared with the state-of-the-art methods. Comparison with the online learning method for logistic loss demonstrates that \textsc{FTRL-AUC} achieves higher AUC scores especially when datasets are imbalanced.


title: A Real-Time Predictive Pedestrian Collision Warning Service for Cooperative Intelligent Transportation Systems Using 3D Pose Estimation

abstract: Minimizing traffic accidents between vehicles and pedestrians is one of the primary research goals in intelligent transportation systems. To achieve the goal, pedestrian orientation recognition and prediction of pedestrian's crossing or not-crossing intention play a central role. Contemporary approaches do not guarantee satisfactory performance due to limited field-of-view, lack of generalization, and high computational complexity. To overcome these limitations, we propose a real-time predictive pedestrian collision warning service (P2CWS) for two tasks: pedestrian orientation recognition (100.53 FPS) and intention prediction (35.76 FPS). Our framework obtains satisfying generalization over multiple sites because of the proposed site-independent features. At the center of the feature extraction lies 3D pose estimation. The 3D pose analysis enables robust and accurate recognition of pedestrian orientations and prediction of intentions over multiple sites. The proposed vision framework realizes 89.3% accuracy in the behavior recognition task on the TUD dataset without any training process and 91.28% accuracy in intention prediction on our dataset achieving new state-of-the-art performance. To contribute to the corresponding research community, we make our source codes public which are available at https://github.com/Uehwan/VisionForPedestrian


title: Stabilization of NbTe3, VTe3, and TiTe3 via Nanotube Encapsulation

abstract: The structure of MX3 transition metal trichalcogenides (TMTs, with M a transition metal and X a chalcogen) is typified by one-dimensional (1D) chains weakly bound together via van der Waals interactions. This structural motif is common across a range of M and X atoms (e.g. NbSe3, HfTe3, TaS3), but not all M and X combinations are stable. We report here that three new members of the MX3 family which are not stable in bulk, specifically NbTe3, VTe3, and TiTe3, can be synthesized in the few- to single-chain limit via nano-confined growth within the stabilizing cavity of multi-walled carbon nanotubes. Transmission electron microscopy (TEM) and atomic-resolution scanning transmission electron microscopy (STEM) reveal the chain-like nature and the detailed atomic structure. The synthesized materials exhibit behavior unique to few-chain quasi-1D structures, such as multi-chain spiraling and a trigonal anti-prismatic rocking distortion in the single-chain limit. Density functional theory (DFT) calculations provide insight into the crystal structure and stability of the materials, as well as their electronic structure.


title: Third order M{\o}ller-Plesset theory made more useful? The role of density functional theory orbitals

abstract: The practical utility of M{\o}ller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error (J. Phys. Chem. Lett. 10, 4170, 2019). In this work we examine whether density functional theory (DFT) optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, non-covalent interactions, and dipole moments relative to standard HF based MP theory. Indeed MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled cluster singles and doubles (CCSD) for several datasets. We also found that the results are not particularly functional sensitive in most cases, (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, non-iterative $O(N^6)$ scaling wave function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.


title: Analytical Inverse For The Symmetric Circulant Tridiagonal Matrix

abstract: Finding the inverse of a matrix is an open problem especially when it comes to engineering problems due to their complexity and running time (cost) of matrix inversion algorithms. An optimum strategy to invert a matrix is, first, to reduce the matrix to a simple form, only then beginning a mathematical procedure. For symmetric matrices, the preferred simple form is tridiagonal. This makes tridiagonal matrices of high interest in applied mathematics and engineering problems. This study presents a time efficient, exact analytical approach for finding the inverse, decomposition, and solving linear systems of equations where symmetric circulant matrix appears. This matrix appears in many researches and it is different from ordinary tridiagonal matrices as there are two corner elements. For finding the inverse matrix, a set of matrices are introduced that any symmetric circulant matrix could be decomposed into them. After that, the exact analytical inverse of this set is found which gives the inverse of circulant matrix. Moreover, solving linear equations can be carried out using implemented decomposition where this matrix appears as a coefficient matrix. The methods principal strength is that it is as stable as any other direct methods (i.e., execute in a predictable number of operations). It is straightforward, understandable, solid as a rock and an exceptionally good psychological backup for those times when something is going wrong and you think it might be your linear equation solver. The downside of present method, and every other direct method, is the accumulation of round off errors.


title: Three-Dimensional Reconstructions of Coronal Wave Surfaces Using a New Mask-Fitting Method

abstract: Coronal waves are large-scale disturbances often driven by coronal mass ejections (CMEs). We investigate a spectacular wave event on 7 March 2012, which is associated with an X5.4 flare (SOL2012-03-07). By using a running center-median (RCM) filtering method for the detection of temporal variations in extreme ultraviolet (EUV) images, we enhance the EUV disturbance observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) onboard the PRoject for Onboard Autonomy 2 (PROBA2). In coronagraph images, a halo front is observed to be the upper counterpart of the EUV disturbance. Based on the EUV and coronagraph images observed from three different perspectives, we have made three-dimensional (3D) reconstructions of the wave surfaces using a new mask-fitting method. The reconstructions are compared with those obtained from forward-fitting methods. We show that the mask fitting method can reflect the inhomogeneous coronal medium by capturing the concave shape of the shock wave front. Subsequently, we trace the developing concave structure and derive the deprojected wave kinematics. The speed of the 3D-wave nose increases from a low value below a few hundred $\mathrm{km\,s^{-1}}$ to a maximum value of about 3800 $\mathrm{km\,s^{-1}}$, and then slowly decreases afterwards. The concave structure starts to decelerate earlier and has significantly lower speeds than those of the wave nose. We also find that the 3D-wave in the extended corona has a much higher speed than the speed of EUV disturbances across the solar disk.


title: Dyson's Crank and the Mex of Integer Partitions

abstract: Andrews and Newman have recently introduced the notion of the mex of a partition, the smallest positive integer that is not a part. The concept has been used since at least 2011, though, with connections to Frobenius symbols. Recently the parity of the mex has been associated to the crank statistic named by Dyson in 1944. In this note, we extend and strengthen the connection between the crank and mex (along with a new generaliztion of the mex) by proving a number of properties that naturally relate these partition statistics.


title: Hamming OCR: A Locality Sensitive Hashing Neural Network for Scene Text Recognition

abstract: Recently, inspired by Transformer, self-attention-based scene text recognition approaches have achieved outstanding performance. However, we find that the size of model expands rapidly with the lexicon increasing. Specifically, the number of parameters for softmax classification layer and output embedding layer are proportional to the vocabulary size. It hinders the development of a lightweight text recognition model especially applied for Chinese and multiple languages. Thus, we propose a lightweight scene text recognition model named Hamming OCR. In this model, a novel Hamming classifier, which adopts locality sensitive hashing (LSH) algorithm to encode each character, is proposed to replace the softmax regression and the generated LSH code is directly employed to replace the output embedding. We also present a simplified transformer decoder to reduce the number of parameters by removing the feed-forward network and using cross-layer parameter sharing technique. Compared with traditional methods, the number of parameters in both classification and embedding layers is independent on the size of vocabulary, which significantly reduces the storage requirement without loss of accuracy. Experimental results on several datasets, including four public benchmaks and a Chinese text dataset synthesized by SynthText with more than 20,000 characters, shows that Hamming OCR achieves competitive results.


title: LTLf Synthesis under Partial Observability: From Theory to Practice

abstract: LTL synthesis is the problem of synthesizing a reactive system from a formal specification in Linear Temporal Logic. The extension of allowing for partial observability, where the system does not have direct access to all relevant information about the environment, allows generalizing this problem to a wider set of real-world applications, but the difficulty of implementing such an extension in practice means that it has remained in the realm of theory. Recently, it has been demonstrated that restricting LTL synthesis to systems with finite executions by using LTL with finite-horizon semantics (LTLf) allows for significantly simpler implementations in practice. With the conceptual simplicity of LTLf, it becomes possible to explore extensions such as partial observability in practice for the first time. Previous work has analyzed the problem of LTLf synthesis under partial observability theoretically and suggested two possible algorithms, one with 3EXPTIME and another with 2EXPTIME complexity. In this work, we first prove a complexity lower bound conjectured in earlier work. Then, we complement the theoretical analysis by showing how the two algorithms can be integrated in practice into an established framework for LTLf synthesis. We furthermore identify a third, MSO-based, approach enabled by this framework. Our experimental evaluation reveals very different results from what the theory seems to suggest, with the 3EXPTIME algorithm often outperforming the 2EXPTIME approach. Furthermore, as long as it is able to overcome an initial memory bottleneck, the MSO-based approach can often outperforms the others.


title: Symbolic Parity Game Solvers that Yield Winning Strategies

abstract: Parity games play an important role for LTL synthesis as evidenced by recent breakthroughs on LTL synthesis, which rely in part on parity game solving. Yet state space explosion remains a major issue if we want to scale to larger systems or specifications. In order to combat this problem, we need to investigate symbolic methods such as BDDs, which have been successful in the past to tackle exponentially large systems. It is therefore essential to have symbolic parity game solving algorithms, operating using BDDs, that are fast and that can produce the winning strategies used to synthesize the controller in LTL synthesis. Current symbolic parity game solving algorithms do not yield winning strategies. We now propose two symbolic algorithms that yield winning strategies, based on two recently proposed fixpoint algorithms. We implement the algorithms and empirically evaluate them using benchmarks obtained from SYNTCOMP 2020. Our conclusion is that the algorithms are competitive with or faster than an earlier symbolic implementation of Zielonka's recursive algorithm, while also providing the winning strategies.


title: Symbolic Execution + Model Counting + Entropy Maximization = Automatic Search Synthesis

abstract: We present a method of automatically synthesizing steps to solve search problems. Given a specification of a search problem, our approach uses symbolic execution to analyze the specification in order to extract a set of constraints which model the problem. These constraints are used in a process called model counting, which is leveraged to compute probability distributions relating search steps to predicates about an unknown target. The probability distribution functions determine an information gain objective function based on Shannon entropy, which, when maximized, yields the next optimal step of the search. We prove that our algorithm converges to a correct solution, and discuss computational complexity issues. We implemented a domain specific language in which to write search problem specifications, enabling our static analysis phase. Our experiments demonstrate the effectiveness of our approach on a set of search problem case studies inspired by the domains of software security, computational geometry, AI for games, and user preference ranking.


title: A Game Theoretical Semantics for Logics of Nonsense

abstract: Logics of non-sense allow a third truth value to express propositions that are \emph{nonsense}. These logics are ideal formalisms to understand how errors are handled in programs and how they propagate throughout the programs once they appear. In this paper, we give a Hintikkan game semantics for logics of non-sense and prove its correctness. We also discuss how a known solution method in game theory, the iterated elimination of strictly dominated strategies, relates to semantic games for logics of nonsense. Finally, we extend the logics of nonsense only by means of semantic games, developing a new logic of nonsense, and propose a new game semantics for Priest's Logic of Paradox.


title: Creating and using large grids of precalculated model atmospheres for a rapid analysis of stellar spectra

abstract: $Aims.$ We present a database of 43,340 atmospheric models ($\sim$80,000 models at the conclusion of the project) for stars with stellar masses between 9 and 120 $M_{\odot}$, covering the region of the OB main-sequence and Wolf-Rayet (W-R) stars in the Hertzsprung--Russell (H--R) diagram. $Methods.$ The models were calculated using the ABACUS I supercomputer and the stellar atmosphere code CMFGEN. $Results.$ The parameter space has six dimensions: the effective temperature $T_{eff}$, the luminosity $L$, the metallicity $Z$, and three stellar wind parameters: the exponent $\beta$, the terminal velocity $V_{\infty}$, and the volume filling factor $F_{cl}$. For each model, we also calculate synthetic spectra in the UV (900-2000 A), optical (3500-7000 A), and near-IR (10000-40000 A) regions. To facilitate comparison with observations, the synthetic spectra can be rotationally broadened using ROTIN3, by covering vsin(i) velocities between 10 and 350 km s$^{-1}$ with steps of 10 km s$^{-1}$. $Conclusions.$ We also present the results of the reanalysis of $\epsilon$ Ori using our grid to demonstrate the benefits of databases of precalculated models. Our analysis succeeded in reproducing the best-fit parameter ranges of the original study, although our results favor the higher end of the mass-loss range and a lower level of clumping. Our results indirectly suggest that the resonance lines in the UV range are strongly affected by the velocity-space porosity, as has been suggested by recent theoretical calculations and numerical simulations.


title: Bounded Game-Theoretic Semantics for Modal Mu-Calculus and Some Variants

abstract: We introduce a new game-theoretic semantics (GTS) for the modal mu-calculus. Our so-called bounded GTS replaces parity games with alternative evaluation games where only finite paths arise; infinite paths are not needed even when the considered transition system is infinite. The novel games offer alternative approaches to various constructions in the framework of the mu-calculus. For example, they have already been successfully used as a basis for an approach leading to a natural formula size game for the logic. While our main focus is introducing the new GTS, we also consider some applications to demonstrate its uses. For example, we consider a natural model transformation procedure that reduces model checking games to checking a single, fixed formula in the constructed models, and we also use the GTS to identify new alternative variants of the mu-calculus with PTime model checking.


title: Local Higher-Order Fixpoint Iteration

abstract: Local fixpoint iteration describes a technique that restricts fixpoint iteration in function spaces to needed arguments only. It has been studied well for first-order functions in abstract interpretation and also in model checking. Here we consider the problem for least and greatest fixpoints of arbitrary type order. We define an abstract algebra of simply typed higher-order functions with fixpoints that can express fixpoint evaluation problems as they occur routinely in various applications, including program verification. We present an algorithm that realises local fixpoint iteration for such higher-order fixpoints, prove its correctness and study its optimisation potential in the context of several applications.


title: Comparison of Algorithms for Simple Stochastic Games

abstract: Simple stochastic games are turn-based 2.5-player zero-sum graph games with a reachability objective. The problem is to compute the winning probability as well as the optimal strategies of both players. In this paper, we compare the three known classes of algorithms -- value iteration, strategy iteration and quadratic programming -- both theoretically and practically. Further, we suggest several improvements for all algorithms, including the first approach based on quadratic programming that avoids transforming the stochastic game to a stopping one. Our extensive experiments show that these improvements can lead to significant speed-ups. We implemented all algorithms in PRISM-games 3.0, thereby providing the first implementation of quadratic programming for solving simple stochastic games.


title: LTLf Synthesis on Probabilistic Systems

abstract: Many systems are naturally modeled as Markov Decision Processes (MDPs), combining probabilities and strategic actions. Given a model of a system as an MDP and some logical specification of system behavior, the goal of synthesis is to find a policy that maximizes the probability of achieving this behavior. A popular choice for defining behaviors is Linear Temporal Logic (LTL). Policy synthesis on MDPs for properties specified in LTL has been well studied. LTL, however, is defined over infinite traces, while many properties of interest are inherently finite. Linear Temporal Logic over finite traces (LTLf) has been used to express such properties, but no tools exist to solve policy synthesis for MDP behaviors given finite-trace properties. We present two algorithms for solving this synthesis problem: the first via reduction of LTLf to LTL and the second using native tools for LTLf. We compare the scalability of these two approaches for synthesis and show that the native approach offers better scalability compared to existing automaton generation tools for LTL.


title: Topology optimization enhances the distinguishability and reconstructability of electrical resistance tomography based sensors

abstract: In the majority of applications of electrical resistance tomography (ERT) the estimation problem consists of either the estimation of spatial conductivity change over an existing background or the estimation of spatial distribution of conductivity of the entire target, including the background. In some instances however, it is possible to design the background conductivity; an example of such application is the design of ERT-based sensors where the background conductivity can be engineered. In such applications the natural question is whether the background conductivity can be engineered in such a way to increase the distinguishability and further reconstructability of the sensor. The present paper, uses topology optimization to design the background conductivity to achieve optimal distinguishability. Then, ERT reconstructions suggest the enhancements of reconstructability using topology optimized sensor.


title: Canonicity in GFG and Transition-Based Automata

abstract: Minimization of deterministic automata on finite words results in a {\em canonical\/} automaton. For deterministic automata on infinite words, no canonical minimal automaton exists, and a language may have different minimal deterministic B\"uchi (DBW) or co-B\"uchi (DCW) automata. In recent years, researchers have studied {\em good-for-games\/} (GFG) automata -- nondeterministic automata that can resolve their nondeterministic choices in a way that only depends on the past. Several applications of automata in formal methods, most notably synthesis, that are traditionally based on deterministic automata, can instead be based on GFG automata. The {\em minimization\/} problem for DBW and DCW is NP-complete, and it stays NP-complete for GFG B\"uchi and co-B\"uchi automata. On the other hand, minimization of GFG co-B\"uchi automata with {\em transition-based\/} acceptance (GFG-tNCWs) can be solved in polynomial time. In these automata, acceptance is defined by a set $\alpha$ of transitions, and a run is accepting if it traverses transitions in $\alpha$ only finitely often. This raises the question of canonicity of minimal deterministic and GFG automata with transition-based acceptance. In this paper we study this problem. We start with GFG-tNCWs and show that the safe components (that is, these obtained by restricting the transitions to these not in $\alpha$) of all minimal GFG-tNCWs are isomorphic, and that by saturating the automaton with transitions in $\alpha$ we get isomorphism among all minimal GFG-tNCWs. Thus, a canonical form for minimal GFG-tNCWs can be obtained in polynomial time. We continue to DCWs with transition-based acceptance (tDCWs), and their dual tDBWs. We show that here, while no canonical form for minimal automata exists, restricting attention to the safe components is useful, and implies that the only minimal tDCWs that have no canonical form are these for which the transition to the GFG model results in strictly smaller automaton, which do have a canonical minimal form.


title: The Quotient in Preorder Theories

abstract: Seeking the largest solution to an expression of the form A x <= B is a common task in several domains of engineering and computer science. This largest solution is commonly called quotient. Across domains, the meanings of the binary operation and the preorder are quite different, yet the syntax for computing the largest solution is remarkably similar. This paper is about finding a common framework to reason about quotients. We only assume we operate on a preorder endowed with an abstract monotonic multiplication and an involution. We provide a condition, called admissibility, which guarantees the existence of the quotient, and which yields its closed form. We call preordered heaps those structures satisfying the admissibility condition. We show that many existing theories in computer science are preordered heaps, and we are thus able to derive a quotient for them, subsuming existing solutions when available in the literature. We introduce the concept of sieved heaps to deal with structures which are given over multiple domains of definition. We show that sieved heaps also have well-defined quotients.


title: Schizophrenia-mimicking layers outperform conventional neural network layers

abstract: We have reported nanometer-scale three-dimensional studies of brain networks of schizophrenia cases and found that their neurites are thin and tortuous compared to healthy controls. This suggests that connections between distal neurons are suppressed in microcircuits of schizophrenia cases. In this study, we applied these biological findings to the design of schizophrenia-mimicking artificial neural network to simulate the observed connection alteration in the disorder. Neural networks having a "schizophrenia connection layer" in place of a fully connected layer were subjected to image classification tasks using the MNIST and CIFAR-10 datasets. The results revealed that the schizophrenia connection layer is tolerant to overfitting and outperforms a fully connected layer. The outperformance was observed only for networks using band matrices as weight windows, indicating that the shape of the weight matrix is relevant to the network performance. A schizophrenia convolution layer was also tested using the VGG configuration, showing that 60% of the kernel weights of the last three convolution layers can be eliminated without loss of accuracy. The schizophrenia layers can be used instead of conventional layers without any change in the network configuration and training procedures; hence, neural networks can easily take advantage of these layers. The results of this study suggest that the connection alteration found in schizophrenia is not a burden to the brain, but has functional roles in brain performance.


title: Distributing Graph States Across Quantum Networks

abstract: Graph states are an important class of multipartite entangled quantum states. We propose a new approach for distributing graph states across a quantum network. We consider a quantum network consisting of nodes-quantum computers within which local operations are free-and EPR pairs shared between nodes that can continually be generated. We prove upper bounds for our approach on the number of EPR pairs consumed, number of timesteps taken, and amount of classical communication required, all of which are equal to or better than that of prior work. We also reduce the problem of minimizing the number of timesteps taken to distribute a graph state using our approach to a network flow problem having polynomial time complexity.


title: The $\ell$-modular representation of reductive groups over finite local rings of length two

abstract: Let $\mathcal{O}_2$ and $\mathcal{O}'_2$ be two distinct finite local rings of length two with residue field of characteristic $p$. Let $\mathbb{G}(\mathcal{O}_2)$ and $\mathbb{G}(\mathcal{O}'_2)$, be the group of points of any reductive group scheme $\mathbb{G}$ over $\mathbb{Z}$ such that $p$ is very good for $\mathbb{G} \times \mathbb{F}_q$. We prove that there exists an isomorphism of group algebra $K[\mathbb{G}(\mathcal{O}_2)] \cong K[\mathbb{G}(\mathcal{O}'_2)]$, where $K$ is a sufficiently large field of characteristic different from $p$.


title: Demand Responsive Dynamic Pricing Framework for Prosumer Dominated Microgrids using Multiagent Reinforcement Learning

abstract: Demand Response (DR) has a widely recognized potential for improving grid stability and reliability while reducing customers energy bills. However, the conventional DR techniques come with several shortcomings, such as inability to handle operational uncertainties and incurring customer disutility, impeding their wide spread adoption in real-world applications. This paper proposes a new multiagent Reinforcement Learning (RL) based decision-making environment for implementing a Real-Time Pricing (RTP) DR technique in a prosumer dominated microgrid. The proposed technique addresses several shortcomings common to traditional DR methods and provides significant economic benefits to the grid operator and prosumers. To show its better efficacy, the proposed DR method is compared to a baseline traditional operation scenario in a small-scale microgrid system. Finally, investigations on the use of prosumers energy storage capacity in this microgrid highlight the advantages of the proposed method in establishing a balanced market setup.


title: Leveraging Local and Global Descriptors in Parallel to Search Correspondences for Visual Localization

abstract: Visual localization to compute 6DoF camera pose from a given image has wide applications such as in robotics, virtual reality, augmented reality, etc. Two kinds of descriptors are important for the visual localization. One is global descriptors that extract the whole feature from each image. The other is local descriptors that extract the local feature from each image patch usually enclosing a key point. More and more methods of the visual localization have two stages: at first to perform image retrieval by global descriptors and then from the retrieval feedback to make 2D-3D point correspondences by local descriptors. The two stages are in serial for most of the methods. This simple combination has not achieved superiority of fusing local and global descriptors. The 3D points obtained from the retrieval feedback are as the nearest neighbor candidates of the 2D image points only by global descriptors. Each of the 2D image points is also called a query local feature when performing the 2D-3D point correspondences. In this paper, we propose a novel parallel search framework, which leverages advantages of both local and global descriptors to get nearest neighbor candidates of a query local feature. Specifically, besides using deep learning based global descriptors, we also utilize local descriptors to construct random tree structures for obtaining nearest neighbor candidates of the query local feature. We propose a new probabilistic model and a new deep learning based local descriptor when constructing the random trees. A weighted Hamming regularization term to keep discriminativeness after binarization is given in the loss function for the proposed local descriptor. The loss function co-trains both real and binary descriptors of which the results are integrated into the random trees.


title: HiCOMEX: Facial Action Unit Recognition Based on Hierarchy Intensity Distribution and COMEX Relation Learning

abstract: The detection of facial action units (AUs) has been studied as it has the competition due to the wide-ranging applications thereof. In this paper, we propose a novel framework for the AU detection from a single input image by grasping the \textbf{c}o-\textbf{o}ccurrence and \textbf{m}utual \textbf{ex}clusion (COMEX) as well as the intensity distribution among AUs. Our algorithm uses facial landmarks to detect the features of local AUs. The features are input to a bidirectional long short-term memory (BiLSTM) layer for learning the intensity distribution. Afterwards, the new AU feature continuously passed through a self-attention encoding layer and a continuous-state modern Hopfield layer for learning the COMEX relationships. Our experiments on the challenging BP4D and DISFA benchmarks without any external data or pre-trained models yield F1-scores of 63.7\% and 61.8\% respectively, which shows our proposed networks can lead to performance improvement in the AU detection task.


title: Pruning Convolutional Filters using Batch Bridgeout

abstract: State-of-the-art computer vision models are rapidly increasing in capacity, where the number of parameters far exceeds the number required to fit the training set. This results in better optimization and generalization performance. However, the huge size of contemporary models results in large inference costs and limits their use on resource-limited devices. In order to reduce inference costs, convolutional filters in trained neural networks could be pruned to reduce the run-time memory and computational requirements during inference. However, severe post-training pruning results in degraded performance if the training algorithm results in dense weight vectors. We propose the use of Batch Bridgeout, a sparsity inducing stochastic regularization scheme, to train neural networks so that they could be pruned efficiently with minimal degradation in performance. We evaluate the proposed method on common computer vision models VGGNet, ResNet, and Wide-ResNet on the CIFAR image classification task. For all the networks, experimental results show that Batch Bridgeout trained networks achieve higher accuracy across a wide range of pruning intensities compared to Dropout and weight decay regularization.


title: Distributionally Robust Optimization Approaches for a Stochastic Mobile Facility Fleet Sizing, Routing, and Scheduling Problem

abstract: We propose two distributionally robust optimization (DRO) models for a mobile facility (MF) fleet sizing, routing, and scheduling problem (MFRSP) with time-dependent and random demand, as well as methodologies for solving these models. Specifically, given a set of MFs, a planning horizon, and a service region, our models aim to find the number of MFs to use (i.e., fleet size) within the planning horizon and a route and schedule for each MF in the fleet. The objective is to minimize the fixed cost of establishing the MF fleet plus a risk measure (expectation or mean conditional value-at-risk) of the operational cost over all demand distributions defined by an ambiguity set. In the first model, we use an ambiguity set based on the demand's mean, support, and mean absolute deviation. In the second model, we use an ambiguity set that incorporates all distributions within a 1-Wasserstein distance from a reference distribution. To solve the proposed DRO models, we propose a decomposition-based algorithm. In addition, we derive valid lower bound inequalities that efficiently strengthen the master problem in the decomposition algorithm, thus improving convergence. We also derive two families of symmetry-breaking constraints that improve the solvability of the proposed models. Finally, we present extensive computational experiments comparing the operational and computational performance of the proposed models and a stochastic programming model, demonstrating where significant performance improvements could be gained and derive insights into the MFRSP.


title: On the control of interference and diffraction of a 3-level atom in a double-slit scheme with cavity fields

abstract: A double cavity with a quantum mechanical and a classical field is located immediately behind of a double-slit in order to analyse the wave-particle duality. Both fields have common nodes and antinodes through which a three-level atom passes after crossing the double-slit. The atom-field interaction is maximum when the atom crosses a common antinode and path-information can be recorded on the phase of the quantum field. On other hand, if the atom crosses a common node, the interaction is null and no path-information is stored. A quadrature measurement on the quantum field can reveal the path followed by the atom, depending on its initial amplitude $\alpha$ and the classical amplitude $\varepsilon$. In this report we show that the classical radiation acts like a focusing element of the interference and diffraction patterns and how it alters the visibility and distinguishabilily. Furthermore, in our double-slit scheme the two possible paths are correlated with the internal atomic states, which allows us to study the relationship between concurrence and wave-particle duality considering different cases.


title: Observing and modeling the sequential pairwise reactions that drive solid-state ceramic synthesis

abstract: Solid-state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing ceramic synthesis routes is usually a laborious, trial-and-error process, as heterogeneous mixtures of powder precursors often evolve through a complicated series of reaction intermediates. Here, we show that phase evolution from multiple precursors can be modeled as a sequence of pairwise interfacial reactions, with thermodynamic driving forces that can be efficiently calculated using ab initio methods. Using the synthesis of the classic high-temperature superconductor YBa$_2$Cu$_3$O$_{6+x}$ (YBCO) as a representative system, we rationalize how replacing the common BaCO$_3$ precursor with BaO$_2$ redirects phase evolution through a kinetically-facile pathway. Our model is validated from in situ X-ray diffraction and in situ microscopy observations, which show rapid YBCO formation from BaO$_2$ in only 30 minutes. By combining thermodynamic modeling with in situ characterization, we introduce a new computable framework to interpret and ultimately design synthesis pathways to complex ceramic materials.


title: Revisiting Design Choices in Proximal Policy Optimization

abstract: Proximal Policy Optimization (PPO) is a popular deep policy gradient algorithm. In standard implementations, PPO regularizes policy updates with clipped probability ratios, and parameterizes policies with either continuous Gaussian distributions or discrete Softmax distributions. These design choices are widely accepted, and motivated by empirical performance comparisons on MuJoCo and Atari benchmarks. We revisit these practices outside the regime of current benchmarks, and expose three failure modes of standard PPO. We explain why standard design choices are problematic in these cases, and show that alternative choices of surrogate objectives and policy parameterizations can prevent the failure modes. We hope that our work serves as a reminder that many algorithmic design choices in reinforcement learning are tied to specific simulation environments. We should not implicitly accept these choices as a standard part of a more general algorithm.


title: On IPW-based estimation of conditional average treatment effect

abstract: The research in this paper gives a systematic investigation on the asymptotic behaviours of four inverse probability weighting (IPW)-based estimators for conditional average treatment effect, with nonparametrically, semiparametrically, parametrically estimated and true propensity score, respectively. To this end, we first pay a particular attention to semiparametric dimension reduction structure such that we can well study the semiparametric-based estimator that can well alleviate the curse of dimensionality and greatly avoid model misspecification. We also derive some further properties of existing estimator with nonparametrically estimated propensity score. According to their asymptotic variance functions, the studies reveal the general ranking of their asymptotic efficiencies; in which scenarios the asymptotic equivalence can hold; the critical roles of the affiliation of the given covariates in the set of arguments of propensity score, the bandwidth and kernel selections. The results show an essential difference from the IPW-based (unconditional) average treatment effect(ATE). The numerical studies indicate that for high-dimensional paradigms, the semiparametric-based estimator performs well in general {whereas nonparametric-based estimator, even sometimes, parametric-based estimator, is more affected by dimensionality. Some numerical studies are carried out to examine their performances. A real data example is analysed for illustration.


title: Periodicities in an active region correlated with Type III radio bursts observed by Parker Solar Probe

abstract: Context. Periodicities have frequently been reported across many wavelengths in the solar corona. Correlated periods of ~5 minutes, comparable to solar p-modes, are suggestive of coupling between the photosphere and the corona. Aims. Our study investigates whether there are correlations in the periodic behavior of Type III radio bursts, indicative of non-thermal electron acceleration processes, and coronal EUV emission, assessing heating and cooling, in an active region when there are no large flares. Methods. We use coordinated observations of Type III radio bursts from the FIELDS instrument on Parker Solar Probe (PSP), of extreme ultraviolet emissions by the Solar Dynamics Observatory (SDO)/AIA and white light observations by SDO/HMI, and of solar flare x-rays by Nuclear Spectroscopic Telescope Array (NuSTAR) on April 12, 2019. Several methods for assessing periodicities are utilized and compared to validate periods obtained. Results. Periodicities of about 5 minutes in the EUV in several areas of an active region are well correlated with the repetition rate of the Type III radio bursts observed on both PSP and Wind. Detrended 211A and 171A light curves show periodic profiles in multiple locations, with 171A peaks lagging those seen in 211A. This is suggestive of impulsive events that result in heating and then cooling in the lower corona. NuSTAR x-rays provide evidence for at least one microflare during the interval of Type III bursts, but there is not a one-to-one correspondence between the x-rays and the Type-III bursts. Our study provides evidence for periodic acceleration of non-thermal electrons (required to generate Type III radio bursts) when there were no observable flares either in the x-ray data or the EUV. The acceleration process, therefore, must be associated with small impulsive events, perhaps nanoflares.


title: A projection-based model checking for heterogeneous treatment effect

abstract: In this paper, we investigate the hypothesis testing problem that checks whether part of covariates / confounders significantly affect the heterogeneous treatment effect given all covariates. This model checking is particularly useful in the case where there are many collected covariates such that we can possibly alleviate the typical curse of dimensionality. In the test construction procedure, we use a projection-based idea and a nonparametric estimation-based test procedure to construct an aggregated version over all projection directions. The resulting test statistic is then interestingly with no effect from slow convergence rate the nonparametric estimation usually suffers from. This feature makes the test behave like a global smoothing test to have ability to detect a broad class of local alternatives converging to the null at the fastest possible rate in hypothesis testing. Also, the test can inherit the merit of lobal smoothing tests to be sensitive to oscillating alternative models. The performance of the test is examined by numerical studies and the analysis for a real data example for illustration.


title: Perturbative Removal of a Sign Problem

abstract: This paper presents a method for alleviating sign problems in lattice path integrals, including those associated with finite fermion density in relativistic systems. The method makes use of information gained from some systematic expansion -- such as perturbation theory -- in order to accelerate the Monte Carlo. The method is exact, in the sense that no approximation to the lattice path integral is introduced. Thanks to the underlying systematic expansion, the method is systematically improvable, so that an arbitrary reduction in the sign problem can in principle be obtained. The Thirring model (in 0 + 1 and 1 + 1 dimensions) is used to demonstrate the ability of this method to reduce the finite-density sign problem.


title: Permanental Graphs

abstract: The two components for infinite exchangeability of a sequence of distributions $(P_n)$ are (i) consistency, and (ii) finite exchangeability for each $n$. A consequence of the Aldous-Hoover theorem is that any node-exchangeable, subselection-consistent sequence of distributions that describes a randomly evolving network yields a sequence of random graphs whose expected number of edges grows quadratically in the number of nodes. In this note, another notion of consistency is considered, namely, delete-and-repair consistency; it is motivated by the sense in which infinitely exchangeable permutations defined by the Chinese restaurant process (CRP) are consistent. A goal is to exploit delete-and-repair consistency to obtain a nontrivial sequence of distributions on graphs $(P_n)$ that is sparse, exchangeable, and consistent with respect to delete-and-repair, a well known example being the Ewens permutations \cite{tavare}. A generalization of the CRP$(\alpha)$ as a distribution on a directed graph using the $\alpha$-weighted permanent is presented along with the corresponding normalization constant and degree distribution; it is dubbed the Permanental Graph Model (PGM). A negative result is obtained: no setting of parameters in the PGM allows for a consistent sequence $(P_n)$ in the sense of either subselection or delete-and-repair.


title: Betti numbers of weighted oriented graphs

abstract: Let $\mathcal{D}$ be a weighted oriented graph and $I(\mathcal{D})$ be its edge ideal. In this paper, we investigate the Betti numbers of $I(\mathcal{D})$ via upper-Koszul simplicial complexes, Betti splittings and the mapping cone construction. In particular, we provide recursive formulas for the Betti numbers of edge ideals of several classes of weighted oriented graphs. We also identify classes of weighted oriented graphs whose edge ideals have a unique extremal Betti number which allows us to compute the regularity and projective dimension for the identified classes. Furthermore, we characterize the structure of a weighted oriented graph $\mathcal{D}$ on $n$ vertices such that $\textrm{pdim } (R/I(\mathcal{D}))=n$ where $R=k[x_1,\ldots, x_n]$.


title: Optimal energy conversion through anti-adiabatic driving breaking time-reversal symmetry

abstract: Starting with Carnot engine, the ideal efficiency of a heat engine has been associated with quasi-static transformations and vanishingly small output power. Here, we exactly calculate the thermodynamic properties of a isothermal heat engine, in which the working medium is a periodically driven underdamped harmonic oscillator, focusing instead on the opposite, anti-adiabatic limit, where the period of a cycle is the fastest time scale in the problem. We show that in that limit it is possible to approach the ideal energy conversion efficiency $\eta=1$, with finite output power and vanishingly small relative power fluctuations. The simultaneous realization of all the three desiderata of a heat engine is possible thanks to the breaking of time-reversal symmetry. We also show that non-Markovian dynamics can further improve the power-efficiency trade-off.


title: A Multi-Agent Deep Reinforcement Learning Approach for a Distributed Energy Marketplace in Smart Grids

abstract: This paper presents a Reinforcement Learning (RL) based energy market for a prosumer dominated microgrid. The proposed market model facilitates a real-time and demanddependent dynamic pricing environment, which reduces grid costs and improves the economic benefits for prosumers. Furthermore, this market model enables the grid operator to leverage prosumers storage capacity as a dispatchable asset for grid support applications. Simulation results based on the Deep QNetwork (DQN) framework demonstrate significant improvements of the 24-hour accumulative profit for both prosumers and the grid operator, as well as major reductions in grid reserve power utilization.


title: Modeling impurity concentrations in liquid argon detectors

abstract: Impurities in noble liquid detectors used for neutrino and dark matter experiments can significantly impact the quality of data. We present an experimentally verified model for describing the dynamics of impurity distributions in liquid argon (LAr) detectors. The model considers sources, sinks, and transport of impurities within and between the gas and liquid argon phases. Measurements of oxygen concentrations in a 20-L LAr multi-purpose test stand are compared to calculations made with this model to show that an accurate description of the concentrations under various operational conditions can be obtained. A result of this analysis is a determination of Henry's coefficient for oxygen in LAr. These calculations also show that some processes have small effects on the impurity dynamics and excluding them yields a solution as a sum of two exponential terms. This solution provides a simple way to extract Henry's coefficient with negligible approximation error. It is applied to the data and the Henry's coefficient for oxygen in LAr is obtained as 0.84$^{+0.09}_{-0.05}$, consistent with literature results. Based on the analysis of the data with the model, we further suggest that, for a large liquid argon detector, barriers to flow ("baffles") installed in the gas phase to restrict flow can help reduce the ultimate impurity concentration in the LAr.


title: Joint routing and pricing control in congested mixed autonomy networks

abstract: Routing controllability of connected and autonomous vehicles (CAVs) has been shown to reduce the adverse effects of selfish routing on the network efficiency. However, the assumption that CAV owners would readily allow themselves to be controlled externally by a central agency for the good of the system is unrealistic. In this paper, we propose a joint routing and pricing control scheme that aims to incentivize CAVs to seek centrally controlled system-optimal (SO) routing by saving on tolls while user equilibrium (UE) seeking human-driven vehicles (HVs) are subject to a congestion charge. The problem is formulated as a bi-level optimization program where the upper level optimizes the dynamic toll rates using the network fundamental diagram (NFD) and the lower level is a mixed equilibrium simulation-based dynamic traffic assignment model (SBDTA) considering different combinations of SO-seeking CAVs. We apply a feedback-based controller to solve for the optimal spatially differentiated distance-based congestion charge from which SO-seeking CAVs are exempt; but UE-seeking HVs are subject to the charge for entering the city center. To capture the distinct microscopic behavior of CAVs in the mixed autonomy traffic, we also implement an adaptive link fundamental diagram (FD) within the SBDTA model. The proposed joint control scheme encourages CAV owners to seek SO routing resulting in less total system travel time. It also discourages UE-seeking HVs from congesting the city center. We demonstrate the performance of the proposed scheme in both a small network and a large-scale network of Melbourne, Australia.


title: On the Speed of Sound in Hyperonic Stars

abstract: We build upon the remarkable, model independent constraints on the equation of state of dense baryonic matter established recently by Annala et al. [1]. Using the quark-meson coupling model, an approach to nuclear structure based upon the self-consistent adjustment of hadron structure to the local meson fields, we show that, once hyperons are allowed to appear in dense matter in \b{eta}-equilibrium, the equation of state is consistent with those constraints. As a result, while one cannot rule out the occurence of quark matter in the cores of massive neutron stars, the available constraints are also compatible with the presence of hyperons.


title: Counting perverse coherent systems on Calabi-Yau 4-folds

abstract: Nagao-Nakajima introduced counting invariants of stable perverse coherent systems on small resolutions of Calabi-Yau 3-folds and determined them on the resolved conifold. Their invariants recover DT/PT invariants and Szendr\"oi's non-commutative invariants in some chambers of stability conditions. In this paper, we study an analogue of their work on Calabi-Yau 4-folds. We define counting invariants for stable perverse coherent systems using primary insertions and compute them in all chambers of stability conditions. We also study counting invariants of local resolved conifold $\mathcal{O}_{\mathbb{P}^1}(-1,-1,0)$ defined using torus localization and tautological insertions. We conjecture a wall-crossing formula for them, which upon dimensional reduction recovers Nagao-Nakajima's wall-crossing formula on resolved conifold.


title: The impact of porous walls on the rheology of suspensions

abstract: We study the effect of isotropic porous walls on a plane Couette flow laden with spherical and rigid particles. We perform a parametric study varying the volume fraction between $0$ and $30\%$, the porosity between $0.3$ and $0.9$ and the non-dimensional permeability between $0$ and $7.9 \times 10^{-3}$ We find that the porous walls induce a progressive decrease in the suspension effective viscosity as the wall permeability increases. This behavior is explained by the weakening of the wall-blocking effect and by the appearance of a slip velocity at the interface of the porous medium, which reduces the shear rate in the channel. Therefore, particle rotation and the consequent velocity fluctuations in the two phases are dampened, leading to reduced particle interactions and particle stresses. Based on our numerical evidence, we provide a closed set of equations for the suspension viscosity, which can be used to estimate the suspension rheology in the presence of porous walls.


title: Interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal. IV: The $\gamma$ model and its phase diagram at $1<\gamma <2$

abstract: In this paper we continue our analysis of the interplay between the pairing and the non-Fermi liquid behavior in a metal for a set of quantum-critical (QC) systems with an effective dynamical electron-electron interaction $V(\Omega_m) \propto 1/|\Omega_m|^\gamma$ (the $\gamma$-model). In previous papers we studied the cases $0<\gamma <1$ and $\gamma \approx 1$. We argued that the pairing by a gapless boson is fundamentally different from BCS/Eliashberg pairing by a massive boson as for the former there exists an infinite number of topologically distinct solutions for the gap function $\Delta_n (\omega_m)$ at $T=0$ ($n=0,1,2...$), each with its own condensation energy $E_{c,n}$. Here we extend the analysis to larger $1< \gamma <2$. We argue that the discrete set of solutions survives, and the spectrum of $E_{c,n}$ gets progressively denser as $\gamma$ approaches $2$ and eventually becomes continuous at $\gamma \to 2$. This increases the strength of "longitudinal" gap fluctuations, which tend to reduce the actual superconducting $T_c$ and give rise to a pseudogap region of preformed pairs. We also detect two features on the real axis for $\gamma >1$ which become critical at $\gamma\to 2$. First, the density of states evolves towards a set of discrete $\delta-$functions. Second, an array of dynamical vortices emerges in the upper frequency half-plane. These two features come about because on a real axis, the real part of the interaction, $V' (\Omega) \propto \cos(\pi \gamma/2)/|\Omega|^\gamma$, becomes repulsive for $\gamma >1$, and the imaginary $V^{''} (\Omega) \propto \sin(\pi \gamma/2)/|\Omega|^\gamma$, gets progressively smaller at $\gamma \to 2$. The features on the real axis are consistent with the development of a continuum spectrum of $E_{c,n}$ obtained using $\Delta_n (\omega_m)$ on the Matsubara axis. We consider the case $\gamma =2$ separately in the next paper.


title: SPARC-LDPC Coding for MIMO Massive Unsourced Random Access

abstract: A joint sparse-regression-code (SPARC) and low-density-parity-check (LDPC) coding scheme for multiple-input multiple-output (MIMO) massive unsourced random access (URA) is proposed in this paper. Different from the state-of-the-art covariance-based maximum likelihood (CB-ML) detection scheme, we first split users' messages into two parts. The former part is encoded by SPARCs and tasked to recover part of the messages, the corresponding channel coefficients as well as the interleaving patterns by compressed sensing. The latter part is coded by LDPC codes and then interleaved by the interleave-division multiple access (IDMA) scheme. The decoding of the latter part is based on belief propagation (BP) joint with successive interference cancellation (SIC). Numerical results show our scheme outperforms the CB-ML scheme when the number of antennas at the base station is smaller than that of active users. The complexity of our scheme is with the order $\mathcal{O}\left(2^{B_p}ML+\widehat{K}ML\right)$ and lower than the CB-ML scheme. Moreover, our scheme has higher spectral efficiency (nearly $15$ times larger) than CB-ML as we only split messages into two parts.


title: Isomorphisms Between Bestsvina--Brady Groups and Right-angled Artin Groups

abstract: We give a combinatorial criterion for determining which Bestvina--Brady group is isomorphic to a right-angled Artin group.


title: Coherent X-rays with Tunable Time-Dependent Polarization

abstract: We describe a method for producing high power, coherent x-ray pulses from a free electron laser with femtosecond scale periodic temporal modulation of the polarization vector. This approach relies on the generation of a temporal intensity modulation after self seeding either by modulating the seed intensity or the beam current. After generating a coherent temporally modulated $s$-polarization pulse, the electron beam is delayed by half a modulation period and sent into a short orthogonally oriented undulator, serving as a $p$-polarization afterburner. We provide simulations of three configurations for realizing this polarization switching, namely, enhanced self seeding with an intensity modulation generated by 2 color self seeding, enhanced self seeding of a current modulated bunch, and regular self seeding of a current modulated bunch. Start to end simulations for the Linac Coherent Light Source-II are provided for the latter.


title: The radiative width of the Hoyle state from $\gamma$-ray spectroscopy

abstract: The cascading 3.21 MeV and 4.44 MeV electric quadrupole transitions have been observed from the Hoyle state at 7.65 MeV excitation energy in $^{12}$C, excited by the $^{12}$C(p,p$^{\prime}$) reaction at 10.7 MeV proton energy. From the proton-$\gamma$-$\gamma$ triple coincidence data, a value of ${\Gamma_{\rm rad}}/{\Gamma}=6.2(6) \times 10^{-4}$ was obtained for the radiative branching ratio. Using our results, together with ${\Gamma_{\pi}^{E0}}/{\Gamma}$ from Eriksen et al., Phys. Rev. C 102, 024320 and the currently adopted $\Gamma_{\pi}(E0)$ values, the radiative width of the Hoyle state is determined as $\Gamma_{\rm rad}=5.1(6) \times 10^{-3}$ eV. This value is about 34% higher than the currently adopted value and will impact on models of stellar evolution and nucleosynthesis.


title: CLASS: Cross-Level Attention and Supervision for Salient Objects Detection

abstract: Salient object detection (SOD) is a fundamental computer vision task. Recently, with the revival of deep neural networks, SOD has made great progresses. However, there still exist two thorny issues that cannot be well addressed by existing methods, indistinguishable regions and complex structures. To address these two issues, in this paper we propose a novel deep network for accurate SOD, named CLASS. First, in order to leverage the different advantages of low-level and high-level features, we propose a novel non-local cross-level attention (CLA), which can capture the long-range feature dependencies to enhance the distinction of complete salient object. Second, a novel cross-level supervision (CLS) is designed to learn complementary context for complex structures through pixel-level, region-level and object-level. Then the fine structures and boundaries of salient objects can be well restored. In experiments, with the proposed CLA and CLS, our CLASS net. consistently outperforms 13 state-of-the-art methods on five datasets.


title: Portable high-order finite element kernels I: Streaming Operations

abstract: This paper is devoted to the development of highly efficient kernels performing vector operations relevant in linear system solvers. In particular, we focus on the low arithmetic intensity operations (i.e., streaming operations) performed within the conjugate gradient iterative method, using the parameters specified in the CEED benchmark problems for high-order hexahedral finite elements. We propose a suite of new Benchmark Streaming tests to focus on the distinct streaming operations which must be performed. We implemented these new tests using the OCCA abstraction framework to demonstrate portability of these streaming operations on different GPU architectures, and propose a simple performance model for such kernels which can accurately capture data movement rates as well as kernel launch costs.


title: Pocket Diagnosis: Secure Federated Learning against Poisoning Attack in the Cloud

abstract: Federated learning has become prevalent in medical diagnosis due to its effectiveness in training a federated model among multiple health institutions (i.e. Data Islands (DIs)). However, increasingly massive DI-level poisoning attacks have shed light on a vulnerability in federated learning, which inject poisoned data into certain DIs to corrupt the availability of the federated model. Previous works on federated learning have been inadequate in ensuring the privacy of DIs and the availability of the final federated model. In this paper, we design a secure federated learning mechanism with multiple keys to prevent DI-level poisoning attacks for medical diagnosis, called SFPA. Concretely, SFPA provides privacy-preserving random forest-based federated learning by using the multi-key secure computation, which guarantees the confidentiality of DI-related information. Meanwhile, a secure defense strategy over encrypted locally-submitted models is proposed to defense DI-level poisoning attacks. Finally, our formal security analysis and empirical tests on a public cloud platform demonstrate the security and efficiency of SFPA as well as its capability of resisting DI-level poisoning attacks.


title: Finding high-order Hadamard matrices by using quantum computers

abstract: Solving hard problems is one of the most important issues in computing to be addressed by a quantum computer. Previously, we have shown that the H-SEARCH; which is the problem of finding a Hadamard matrix (H-matrix) among all possible binary matrices of corresponding order, is a hard problem that can be solved by a quantum computer. However, due to the limitation on the number of qubits and connections in present day quantum processors, only low orders H-SEARCH are implementable. In this paper, we show that by adopting classical construction/search techniques of the H-matrix, we can develop new quantum computing methods to find higher order H-matrices. Especially, the Turyn-based quantum computing method can be further developed to find an arbitrarily high order H-matrix by balancing the classical and quantum resources. This method is potentially capable to find some unknown H-matrices of practical and scientific interests, where a classical computer alone cannot do because of the exponential grow of the complexity. We present some results of finding H-matrix of order more than one hundred and a prototypical experiment to find even higher order matrix by using the classical-quantum resource balancing method. Although heuristic optimizations generally only achieve approximate solutions, whereas the exact one should be determined by exhaustive listing; which is difficult to perform, in the H-SEARCH we can assure such exactness in polynomial time by checking the orthogonality of the solution. Since quantum advantage over the classical computing should have been measured by comparing the performance in solving a problem up to a definitive solution, the proposed method may lead to an alternate route for demonstrating practical quantum supremacy in the near future.


title: Educational and Outreach Resource for Astroparticle Physics

abstract: The modern astrophysics is moving towards the enlarging of experiments and combining the channels for detecting the highest energy processes in the Universe. To obtain reliable data, the experiments should operate within several decades, which means that the data will be obtained and analyzed by several generations of physicists. Thus, for the stability of the experiments, it is necessary to properly maintain not only the data life cycle, but also the human aspects, for example, attracting, learning and continuity. To this end, an educational and outreach resource has been deployed in the framework of German-Russian Astroparticle Data Life Cycle Initiative (GRADLCI).


title: On Quantum Simulation Of Cosmic Inflation

abstract: In this paper, we generalize Jordan-Lee-Preskill, an algorithm for simulating flat-space quantum field theories, to 3+1 dimensional inflationary spacetime. The generalized algorithm contains the encoding treatment, the initial state preparation, the inflation process, and the quantum measurement of cosmological observables at late time. The algorithm is helpful for obtaining predictions of cosmic non-Gaussianities, serving as useful benchmark problems for quantum devices, and checking assumptions made about interacting vacuum in the inflationary perturbation theory. Components of our work also include a detailed discussion about the lattice regularization of the cosmic perturbation theory, a detailed discussion about the in-in formalism, a discussion about encoding using the HKLL-type formula that might apply for both dS and AdS spacetimes, a discussion about bounding curvature perturbations, a description of the three-party Trotter simulation algorithm for time-dependent Hamiltonians, a ground state projection algorithm for simulating gapless theories, a discussion about the quantum-extended Church-Turing Thesis, and a discussion about simulating cosmic reheating in quantum devices.


title: Stochastic Generalized Lotka-Volterra Model with An Application to Learning Microbial Community Structures

abstract: Inferring microbial community structure based on temporal metagenomics data is an important goal in microbiome studies. The deterministic generalized Lotka-Volterra differential (GLV) equations have been used to model the dynamics of microbial data. However, these approaches fail to take random environmental fluctuations into account, which may negatively impact the estimates. We propose a new stochastic GLV (SGLV) differential equation model, where the random perturbations of Brownian motion in the model can naturally account for the external environmental effects on the microbial community. We establish new conditions and show various mathematical properties of the solutions including general existence and uniqueness, stationary distribution, and ergodicity. We further develop approximate maximum likelihood estimators based on discrete observations and systematically investigate the consistency and asymptotic normality of the proposed estimators. Our method is demonstrated through simulation studies and an application to the well-known "moving picture" temporal microbial dataset.


title: A Coded Caching Scheme with Linear Sub-packetization and its Application to Multi-Access Coded Caching

abstract: This paper addresses the problem of exponentially increasing sub-packetization with the number of users in a centralized coded caching system by introducing a new coded caching scheme inspired by the symmetric neighboring consecutive side information index coding problem. The scheme has a placement policy where the number of sub-packets required grows only linearly with the number of users, with no restriction on file size, and a delivery policy which is instantaneously decodable. Further, an application of the new delivery scheme in a multi-access coded caching set-up is studied and a few results in that direction are presented. In particular, in the multi-access set-up, for cases where optimality rate-memory trade-off characterizations are available, it is shown that the new delivery scheme achieves optimal or near-optimal rates.


title: FusionStitching: Boosting Memory Intensive Computations for Deep Learning Workloads

abstract: We show in this work that memory intensive computations can result in severe performance problems due to off-chip memory access and CPU-GPU context switch overheads in a wide range of deep learning models. For this problem, current just-in-time (JIT) kernel fusion and code generation techniques have limitations, such as rough fusion plan exploration strategies and limited code generation ability. We propose FusionStitching, a deep learning compiler capable of fusing memory intensive operators, with varied data dependencies and non-homogeneous parallelism, into large GPU kernels to reduce global memory access and context switch overhead automatically. FusionStitching widens the range of operation combinations that fusion can target beyond previous JIT works by introducing data reuse of intermediate values. It explores large fusion spaces to decide optimal fusion plans with considerations of memory access costs, kernel calls and resource usage constraints. FusionStitching tunes the optimal stitching scheme with a domain-specific cost model efficiently. Experimental results show that FusionStitching can reach up to 2.21x speedup compared to state-of-the-art, with 1.45x on average. Besides these experimental results, we integrated our approach into a compiler product and deployed it onto a production cluster for AI workloads with thousands of GPUs. The system has been in operation for more than 4 months and saves 7,000 GPU hours on average for approximately 30,000 tasks per month.


title: Screening Effects on Electron Capture Rates and Type Ia Supernova Nucleosynthesis

abstract: Type Ia supernovae (SNe Ia) are believed to be a thermonuclear explosion of a white dwarf, but the mass of their progenitors is still an open problem. In near-Chandrasekhar-mass (near-M_Ch) models of SNe Ia, the central density reaches >10^9 g cm^{-3}. The electron chemical potential becomes higher than the Q-values of electron capture (EC) transitions between fp-shell nuclei, so a portion of the available electrons is captured by iron group elements and thus neutron-rich isotopes are formed. Since EC reaction rates are sensitive to the density, the degree of neutronization is a key to distinguish near- and sub-M_Ch models. In order to compare observations and theoretical models, an accurate treatment of EC reactions is necessary. In previous theoretical works, however, effects of electron screening on ECs are ignored. Screening lowers EC rates and thus leads to a higher electron fraction. We implement electron screening on ECs to calculate explosive SN Ia nucleosynthesis in a near-M_Ch single degenerate model. It is found that some of neutron-rich nuclear abundances, namely those of 46,48Ca, 50Ti, 54Cr, 58Fe, 64Ni and 67,70Zn, decrease when screening effects on ECs are considered. Of these, 50Ti, 54Cr and 58Fe are particularly interesting because a significant portion of the solar abundance of these nuclei is presumed to originate from SNe Ia. We conclude that implementing the screening effect on ECs in modern SN Ia models is desirable to precisely calculate abundances of neutron-rich nuclides.


title: Time-Optimal Self-Stabilizing Leader Election on Rings in Population Protocols

abstract: We propose a self-stabilizing leader election protocol on directed rings in the model of population protocols. Given an upper bound $N$ on the population size $n$, the proposed protocol elects a unique leader within $O(nN)$ expected steps starting from any configuration and uses $O(N)$ states. This convergence time is optimal if a given upper bound $N$ is asymptotically tight, i.e., $N=O(n)$.


title: Invariance principle for random walks on dynamically averaging random conductances

abstract: We prove an invariance principle for continuous-time random walks in a dynamically averaging environment on $\mathbb Z$. In the beginning, the conductances may fluctuate substantially, but we assume that as time proceeds, the fluctuations decrease according to a typical diffusive scaling and eventually approach constant unit conductances. The proof relies on a coupling with the standard continuous time simple random walk.


title: Gamow vectors formalism applied to the Loschmidt echo

abstract: Gamow vectors have been developed in order to give a mathematical description for quantum decay phenomena. Mainly, they have been applied to radioactive phenomena, scattering and to some decoherence models. They play a crucial role in the description of quantum irreversible processes, and in the formulation of time asymmetry in quantum mechanics. In this paper, we use this formalism to describe a well-known phenomenon of irreversibility: the Loschmidt echo. The standard approach considers that the irreversibility of this phenomenon is the result of an additional term in the backward Hamiltonian. Here, we use the non-Hermitian formalism, where the time evolution is non-unitary. Additionally, we compare the characteristic decay times of this phenomenon with the decoherence ones. We conclude that the Loschmidt echo and the decoherence can be considered as two aspects of the same phenomenon, and that there is a mathematical relationship between their corresponding characteristic times.


title: Semantics of a Relational {\lambda}-Calculus (Extended Version)

abstract: We extend the {\lambda}-calculus with constructs suitable for relational and functional-logic programming: non-deterministic choice, fresh variable introduction, and unification of expressions. In order to be able to unify {\lambda}-expressions and still obtain a confluent theory, we depart from related approaches, such as {\lambda}Prolog, in that we do not attempt to solve higher-order unification. Instead, abstractions are decorated with a location, which intuitively may be understood as its memory address, and we impose a simple coherence invariant: abstractions in the same location must be equal. This allows us to formulate a confluent small-step operational semantics which only performs first-order unification and does not require strong evaluation (below lambdas). We study a simply typed version of the system. Moreover, a denotational semantics for the calculus is proposed and reduction is shown to be sound with respect to the denotational semantics.