title: Towards the Combination of Model Checking and Runtime Verification on Multi-Agent Systems

abstract: Multi-Agent Systems (MAS) are notoriously complex and hard to verify. In fact, it is not trivial to model a MAS, and even when a model is built, it is not always possible to verify, in a formal way, that it is actually behaving as we expect. Usually, it is relevant to know whether an agent is capable of fulfilling its own goals. One possible way to check this is through Model Checking. Specifically, by verifying Alternating-time Temporal Logic (ATL) properties, where the notion of strategies for achieving goals can be described. Unfortunately, the resulting model checking problem is not decidable in general. In this paper, we present a verification procedure based on combining Model Checking and Runtime Verification, where sub-models of the MAS model belonging to decidable fragments are verified by a model checker, and runtime monitors are used to verify the rest. Furthermore, we implement our technique and show experimental results.


title: First Law of Mechanics for Spinning Compact Binaries: Dipolar Order

abstract: Building upon the Noether charge formalism of Iyer and Wald, we derive a variational formula for spacetimes admitting a Killing vector field, for a generic energy-momentum distribution with compact support. Applying this general result to the particular case of a binary system of spinning compact objects moving along an exactly circular orbit, modelled using the multipolar gravitational skeleton formalism, we derive a first law of compact binary mechanics at dipolar order. We prove the equivalence of this new result with the canonical Hamiltonian first law previously derived for binary systems of spinning compact objects, for spins colinear with the orbital angular momentum. This paper paves the way to an extension of the first law of binary mechanics to the next quadrupolar order, thereby accounting for the spin-induced and tidally-induced deformability of the compact bodies.


title: Improving Molecular Contrastive Learning via Faulty Negative Mitigation and Decomposed Fragment Contrast

abstract: Deep learning has been a prevalence in computational chemistry and widely implemented in molecule property predictions. Recently, self-supervised learning (SSL), especially contrastive learning (CL), gathers growing attention for the potential to learn molecular representations that generalize to the gigantic chemical space. Unlike supervised learning, SSL can directly leverage large unlabeled data, which greatly reduces the effort to acquire molecular property labels through costly and time-consuming simulations or experiments. However, most molecular SSL methods borrow the insights from the machine learning community but neglect the unique cheminformatics (e.g., molecular fingerprints) and multi-level graphical structures (e.g., functional groups) of molecules. In this work, we propose iMolCLR: improvement of Molecular Contrastive Learning of Representations with graph neural networks (GNNs) in two aspects, (1) mitigating faulty negative contrastive instances via considering cheminformatics similarities between molecule pairs; (2) fragment-level contrasting between intra- and inter-molecule substructures decomposed from molecules. Experiments have shown that the proposed strategies significantly improve the performance of GNN models on various challenging molecular property predictions. In comparison to the previous CL framework, iMolCLR demonstrates an averaged 1.3% improvement of ROC-AUC on 7 classification benchmarks and an averaged 4.8% decrease of the error on 5 regression benchmarks. On most benchmarks, the generic GNN pre-trained by iMolCLR rivals or even surpasses supervised learning models with sophisticated architecture designs and engineered features. Further investigations demonstrate that representations learned through iMolCLR intrinsically embed scaffolds and functional groups that can reason molecule similarities.


title: Radiation of Closed Strings between the Parallel Dynamical-Dressed Unstable Dp-Branes

abstract: We introduce a boundary state which is corresponding to a D$p$-brane with tangential dynamics in the presence of the Kalb-Ramond field, a tachyonic field and a $U(1)$ gauge potential in a special gauge. From the interaction of such branes radiation amplitude of a general massless closed string will be computed. The effects of the large distances of the branes on this radiation will be studied. In the large distances of the branes, the possibility of axion radiation will be investigated. Our calculations will be in the context of the bosonic string theory.


title: A Machine Learning Paradigm for Studying Pictorial Realism: Are Constable's Clouds More Real than His Contemporaries?

abstract: The British landscape painter John Constable is considered foundational for the Realist movement in 19th-century European painting. Constable's painted skies, in particular, were seen as remarkably accurate by his contemporaries, an impression shared by many viewers today. Yet, assessing the accuracy of realist paintings like Constable's is subjective or intuitive, even for professional art historians, making it difficult to say with certainty what set Constable's skies apart from those of his contemporaries. Our goal is to contribute to a more objective understanding of Constable's realism. We propose a new machine-learning-based paradigm for studying pictorial realism in an explainable way. Our framework assesses realism by measuring the similarity between clouds painted by artists noted for their skies, like Constable, and photographs of clouds. The experimental results of cloud classification show that Constable approximates more consistently than his contemporaries the formal features of actual clouds in his paintings. The study, as a novel interdisciplinary approach that combines computer vision and machine learning, meteorology, and art history, is a springboard for broader and deeper analyses of pictorial realism.


title: Max Noether Theorem for Singular Curves

abstract: Max Noether's Theorem asserts that if $\omega$ is the dualizing sheaf of a nonsingular nonhyperelliptic projective curve, then the natural morphisms $\text{Sym}^nH^0(\omega)\to H^0(\omega^n)$ are surjective for all $n\geq 1$. The result was extended for Gorenstein curves by many different authors in distinct ways. More recently, it was proved for curves with projectively normal canonical models, and curves whose non-Gorenstein points are bibranch at worse. Based on those works, we address the combinatorics of the general case and extend the result for any integral curve.


title: Complexity of warped conformal field theory

abstract: Warped conformal field theories in two dimensions are exotic nonlocal, Lorentz violating field theories characterized by Virasoro-Kac-Moody symmetries and have attracted a lot of attention as candidate boundary duals to warped AdS$_3$ spacetimes, thereby expanding the scope of holography beyond asymptotically AdS spacetimes. Here we investigate WCFT$_2$\,s using \emph{circuit complexity} as a tool. First we compute the holographic volume complexity (CV) which displays a linear UV divergence structure, more akin to that of a local CFT$_2$ and has a very complicated dependence on the Virasoro central charge $c$ and the $U(1)$ Kac-Moody level parameter $k$. Next we consider circuit complexity based on Virasoro-Kac-Moody symmetry gates where the complexity functional is the geometric (group) action on coadjoint orbits of the Virasoro-Kac-Moody group. We consider a special solution to extremization equations for which complexity scales linearly with ``time''. In the semiclassical limit (large $c,k$, while $c/k$ remains finite and small) both the holographic volume complexity and circuit complexity scales linearly with $k$.


title: Inner ideals of real simple Lie algebras

abstract: A classification up to automorphism of the inner ideals of the real finite-dimensional simple Lie algebras is given, jointly with precise descriptions in the case of the exceptional Lie algebras.


title: Assessment of Cyber-Physical Intrusion Detection and Classification for Industrial Control Systems

abstract: The increasing interaction of industrial control systems (ICSs) with public networks and digital devices introduces new cyber threats to power systems and other critical infrastructure. Recent cyber-physical attacks such as Stuxnet and Irongate revealed unexpected ICS vulnerabilities and a need for improved security measures. Intrusion detection systems constitute a key security technology, which typically monitors cyber network data for detecting malicious activities. However, a central characteristic of modern ICSs is the increasing interdependency of physical and cyber network processes. Thus, the integration of network and physical process data is seen as a promising approach to improve predictability in real-time intrusion detection for ICSs by accounting for physical constraints and underlying process patterns. This work systematically assesses machine learning-based cyber-physical intrusion detection and multi-class classification through a comparison to its purely network data-based counterpart and evaluation of misclassifications and detection delay. Multiple supervised detection and classification pipelines are applied on a recent cyber-physical dataset, which describes various cyber attacks and physical faults on a generic ICS. A key finding is that the integration of physical process data improves detection and classification of all considered attack types. In addition, it enables simultaneous processing of attacks and faults, paving the way for holistic cross-domain root cause identification.


title: Model Calibration of the Liquid Mercury Spallation Target using Evolutionary Neural Networks and Sparse Polynomial Expansions

abstract: The mercury constitutive model predicting the strain and stress in the target vessel plays a central role in improving the lifetime prediction and future target designs of the mercury targets at the Spallation Neutron Source (SNS). We leverage the experiment strain data collected over multiple years to improve the mercury constitutive model through a combination of large-scale simulations of the target behavior and the use of machine learning tools for parameter estimation. We present two interdisciplinary approaches for surrogate-based model calibration of expensive simulations using evolutionary neural networks and sparse polynomial expansions. The experiments and results of the two methods show a very good agreement for the solid mechanics simulation of the mercury spallation target. The proposed methods are used to calibrate the tensile cutoff threshold, mercury density, and mercury speed of sound during intense proton pulse experiments. Using strain experimental data from the mercury target sensors, the newly calibrated simulations achieve 7\% average improvement on the signal prediction accuracy and 8\% reduction in mean absolute error compared to previously reported reference parameters, with some sensors experiencing up to 30\% improvement. The proposed calibrated simulations can significantly aid in fatigue analysis to estimate the mercury target lifetime and integrity, which reduces abrupt target failure and saves a tremendous amount of costs. However, an important conclusion from this work points out to a deficiency in the current constitutive model based on the equation of state in capturing the full physics of the spallation reaction. Given that some of the calibrated parameters that show a good agreement with the experimental data can be nonphysical mercury properties, we need a more advanced two-phase flow model to capture bubble dynamics and mercury cavitation.

categories: physics.comp-ph cs.LG physics.acc-ph stat.AP


title: Smoothness of Directed Chain Stochastic Differential Equations

abstract: We study the smoothness of the solution of the directed chain stochastic differential equations, where each process is affected by its neighborhood process in an infinite directed chain graph, introduced by Detering et al. (2020). Because of the auxiliary process in the chain-like structure, classic methods of Malliavin derivatives are not directly applicable. Namely, we cannot make a connection between the Malliavin derivative and the first order derivative of the state process. It turns out that the partial Malliavin derivatives can be used here to fix this problem.


title: Association schemes on triples over few vertices

abstract: In this paper, we obtain classification results for higher-dimensional analogues of classical association schemes called association schemes on triples (ASTs). We present an algorithm that enumerates all ASTs on a fixed number of vertices whose nontrivial relations are invariant under the action of some group. Applying this algorithm to three, four, and five vertices along with appropriate group actions yields the unique AST over three vertices, the unique symmetric ASTs over four or five vertices, the unique AST over four vertices with two nontrivial relations, and the unique nontrivial circulant AST over five vertices.


title: An exploration of an early gravity transition in light of cosmological tensions

abstract: We study a step-like transition in the value of the effective Planck mass (or effective gravitational constant) on cosmological scales prior to recombination. We employ CMB, BAO, and SNIa data and find they are sufficient to strongly constrain our implementation of the Effective Field Theory of Dark Energy and Modified Gravity, used to model the transition, to a limited parameter space. The data prefer a $\sim 5\%$ shift in the value of the effective Planck mass ($<10 \%$ at $2 \sigma$) prior to recombination. This Transitional Planck Mass (TPM) model is free to undergo its transition at any point over multiple decades of scale factor prior to recombination, $\log_{10}(a) = -5.32^{+0.96}_{-0.72}$ (68\% CL). This lowers the sound horizon at last scattering, which increases the Hubble constant to $71.09 \pm 0.75$ km $\textrm{s}^{-1}\textrm{Mpc}^{-1}$ with a combination of local measurements as prior and to $69.22^{+0.67}_{-0.86}$ km $\textrm{s}^{-1}\textrm{Mpc}^{-1}$ when the prior is excluded. The TPM model improves $\chi^2$ with respect to $\Lambda$CDM by $\Delta \chi^2 = -23.72$ with the $H_0$ prior and $\Delta \chi^2 = -4.8$ without the prior. The model allows for both $H_0 > 70$ km$\textrm{s}^{-1}\textrm{Mpc}^{-1}$ and $S_8 < 0.80$ simultaneously with lower values of $S_8$ due to a reduction in the matter density $\Omega_m$ to offset the increase in $H_0$ relative to $\Lambda$CDM. While this is a particular modified gravity model, studying other variants of modified gravity may be a productive path for potentially resolving cosmological tensions, while avoiding the need for a cosmological constant.


title: ProxSkip: Yes! Local Gradient Steps Provably Lead to Communication Acceleration! Finally!

abstract: We introduce ProxSkip -- a surprisingly simple and provably efficient method for minimizing the sum of a smooth ($f$) and an expensive nonsmooth proximable ($\psi$) function. The canonical approach to solving such problems is via the proximal gradient descent (ProxGD) algorithm, which is based on the evaluation of the gradient of $f$ and the prox operator of $\psi$ in each iteration. In this work we are specifically interested in the regime in which the evaluation of prox is costly relative to the evaluation of the gradient, which is the case in many applications. ProxSkip allows for the expensive prox operator to be skipped in most iterations: while its iteration complexity is $\mathcal{O}\left(\kappa \log \frac{1}{\varepsilon}\right)$, where $\kappa$ is the condition number of $f$, the number of prox evaluations is $\mathcal{O}\left(\sqrt{\kappa} \log \frac{1}{\varepsilon}\right)$ only. Our main motivation comes from federated learning, where evaluation of the gradient operator corresponds to taking a local GD step independently on all devices, and evaluation of prox corresponds to (expensive) communication in the form of gradient averaging. In this context, ProxSkip offers an effective acceleration of communication complexity. Unlike other local gradient-type methods, such as FedAvg, SCAFFOLD, S-Local-GD and FedLin, whose theoretical communication complexity is worse than, or at best matching, that of vanilla GD in the heterogeneous data regime, we obtain a provable and large improvement without any heterogeneity-bounding assumptions.


title: The Anthropic (and Mis-) Principle revisited, Steven Weinberg in Memoriam

abstract: I review the role and meaning of the Anthropic Principle, particularly in its relevance to particle physics.


title: Machine Learning Models in Stock Market Prediction

abstract: The paper focuses on predicting the Nifty 50 Index by using 8 Supervised Machine Learning Models. The techniques used for empirical study are Adaptive Boost (AdaBoost), k-Nearest Neighbors (kNN), Linear Regression (LR), Artificial Neural Network (ANN), Random Forest (RF), Stochastic Gradient Descent (SGD), Support Vector Machine (SVM) and Decision Trees (DT). Experiments are based on historical data of Nifty 50 Index of Indian Stock Market from 22nd April, 1996 to 16th April, 2021, which is time series data of around 25 years. During the period there were 6220 trading days excluding all the non trading days. The entire trading dataset was divided into 4 subsets of different size-25% of entire data, 50% of entire data, 75% of entire data and entire data. Each subset was further divided into 2 parts-training data and testing data. After applying 3 tests- Test on Training Data, Test on Testing Data and Cross Validation Test on each subset, the prediction performance of the used models were compared and after comparison, very interesting results were found. The evaluation results indicate that Adaptive Boost, k- Nearest Neighbors, Random Forest and Decision Trees under performed with increase in the size of data set. Linear Regression and Artificial Neural Network shown almost similar prediction results among all the models but Artificial Neural Network took more time in training and validating the model. Thereafter Support Vector Machine performed better among rest of the models but with increase in the size of data set, Stochastic Gradient Descent performed better than Support Vector Machine.


title: Analytic Method for Estimating Aircraft Fix Displacement from Gyroscope's Allan-Deviation Parameters

abstract: The noise and drift requirements for a navigation-grade gyroscope are widely known, yet there is no simple analytic model of how the noise and drift of a gyroscope influence the fix displacement error (FDE) of an inertial navigation system (INS). This work derives simple analytical expressions for the cross-track and along-track errors of an aircraft whose INS consists solely of a three-axis gyroscope system with perfect knowledge of the vertical direction. The error signal of each gyroscope is Gaussian white noise and drift modeled as a first-order Markov random walk. These expressions provide a straightforward mean of calculating the FDE of an aircraft as a function of the flight duration, velocity, noise amplitude, drift amplitude, and drift's time constant. These expressions are validated with Monte-Carlo simulations of long flights. This model quantifies the noise-versus-drift trade-off for a gyroscope in an inertial navigation system. It can save time when estimating the noise and drift that a gyroscope must exhibit to satisfy a given position-error requirement, or vice versa. They are used in particular to confirm the values, often cited without demonstration, of the noise and drift required to meet the maximum position error of an aircraft imposed by the Federal Aviation Administration's required navigation performance 10 specification. Finally, it demonstrates that using the minimum in the measured Allan deviation of a gyroscope as a metric of the drift is incorrect, because it fails to capture the drift's time constant. The proper metric is the maximum in the Allan deviation.


title: Parameter Identification of a PN-Guided Incoming Missile Using an Improved Multiple-Model Mechanism

abstract: An active defense against an incoming missile requires information of it, including a guidance law parameter and a first-order lateral time constant. To this end, assuming that a missile with a proportional navigation (PN) guidance law attempts to attack an aerial target with bang-bang evasive maneuvers, a parameter identification model based on the gated recurrent unit (GRU) neural network is built in this paper. The analytic identification solutions for the guidance law parameter and the first-order lateral time constant are derived. The inputs of the identification model are available kinematic information between the aircraft and the missile, while the outputs contain the regression results of missile parameters. To increase the training speed and the identification accuracy of the Model, an output processing method called improved multiplemodel mechanism (IMMM) is proposed in this paper. The effectiveness of IMMM and the performance of the established model are demonstrated through numerical simulations under various engagement scenarios.


title: Optimal Redundancy Allocation in Coherent Systems with Heterogeneous Dependent Components

abstract: This paper is concerned with the optimal number of redundant allocation to $n$-component coherent systems consist of heterogeneous dependent components. We assume that the system is built of $L$ groups of different components, $L\geq 1$, where there are $n_i$ components in group $i$, and $\sum_{i=1}^{L}n_i=n$. The problem of interest is to allocate $v_i$ active redundant components to each component of type $i$, $i=1,\dots, L$. To get the optimal values of $v_i$, we propose two cost-based criteria. One of them is introduced based on the costs of renewing the failed components and the costs of refreshing the alive ones at the system failure time. The other criterion is proposed based on the costs of replacing the system at its failure time or at a predetermined time $\tau$, whichever occurs first. The expressions for the proposed functions are derived using the mixture representation of the system reliability function based on the notion of survival signature. We assume that a given copula function models the dependency structure between the components. In the particular case that the system is a series-parallel structure, we provide the formulas for the proposed cost-based functions. The results are discussed numerically for some specific coherent systems.


title: Towards Digital Twin Oriented Modelling of Complex Networked Systems and Their Dynamics: A Comprehensive Survey

abstract: This paper aims to provide a comprehensive critical overview on how entities and their interactions in Complex Networked Systems (CNS) are modelled across disciplines as they approach their ultimate goal of creating a Digital Twin (DT) that perfectly matches the reality. We propose a new framework to conceptually compare diverse existing modelling paradigms from different perspectives and create unified assessment criteria to assess their respective capabilities of reaching such an ultimate goal. Using the proposed criteria, we also appraise how far the reviewed current state-of-the-art approaches are from the idealised DTs. We also identify and propose potential directions and ways of building a DT-orientated CNS based on the convergence and integration of CNS and DT utilising a variety of cross-disciplinary techniques.


title: Playing against no-regret players

abstract: In increasingly different contexts, it happens that a human player has to interact with artificial players who make decisions following decision-making algorithms. How should the human player play against these algorithms to maximize his utility? Does anything change if he faces one or more artificial players? The main goal of the paper is to answer these two questions. Consider n-player games in normal form repeated over time, where we call the human player optimizer, and the (n -- 1) artificial players, learners. We assume that learners play no-regret algorithms, a class of algorithms widely used in online learning and decision-making. In these games, we consider the concept of Stackelberg equilibrium. In a recent paper, Deng, Schneider, and Sivan have shown that in a 2-player game the optimizer can always guarantee an expected cumulative utility of at least the Stackelberg value per round. In our first result, we show, with counterexamples, that this result is no longer true if the optimizer has to face more than one player. Therefore, we generalize the definition of Stackelberg equilibrium introducing the concept of correlated Stackelberg equilibrium. Finally, in the main result, we prove that the optimizer can guarantee at least the correlated Stackelberg value per round. Moreover, using a version of the strong law of large numbers, we show that our result is also true almost surely for the optimizer utility instead of the optimizer's expected utility.


title: Migration-Based Synchronization

abstract: A fundamental challenge in multi- and many-core systems is the correct execution of concurrent access to shared data. A common drawback from existing synchronization mechanisms is the loss of data locality as the shared data is transferred between the accessing cores. In real-time systems, this is especially important as knowledge about data access times is crucial to establish bounds on execution times and guarantee the meeting of deadlines.We propose in this paper a refinement of our previously sketched approach of Migration-Based Synchronization (MBS) as well as its first practical implementation. The core concept of MBS is the replacement of data migration with control-flow migration to achieve synchronized memory accesses with guaranteed data locality. This leads to both shorter and more predictable execution times for critical sections. As MBS can be used as a substitute for classical locks, it can be employed in legacy applications without code alterations.We further examine how the gained data locality improves the results of worst-case timing analyses and results in tighter bounds on execution and response time. We reason about the similarity of MBS to existing synchronization approaches and how it enables us to reuse existing analysis techniques.Finally, we evaluate our prototype implementation, showing that MBS can exploit data locality with similar overheads as traditional locking mechanisms.


title: k^{th} order Slant Hankel Operators on the Polydisk

abstract: In this paper, we initiate the notion of k^{th} order slant Hankel operators on L^2(T^n) for k greater than or equal to 2 and n greater than or equal to 1 where T^n denotes the n-torus. We give the necessary and sufficient condition for a bounded operator on L^2(T^n) to be a k^{th} order slant Hankel and discuss their commutative, compactness, hyponormal and isometric property.


title: Snowflake Point Deconvolution for Point Cloud Completion and Generation with Skip-Transformer

abstract: Most existing point cloud completion methods suffer from the discrete nature of point clouds and the unstructured prediction of points in local regions, which makes it difficult to reveal fine local geometric details. To resolve this issue, we propose SnowflakeNet with snowflake point deconvolution (SPD) to generate complete point clouds. SPD models the generation of point clouds as the snowflake-like growth of points, where child points are generated progressively by splitting their parent points after each SPD. Our insight into the detailed geometry is to introduce a skip-transformer in the SPD to learn the point splitting patterns that can best fit the local regions. The skip-transformer leverages attention mechanism to summarize the splitting patterns used in the previous SPD layer to produce the splitting in the current layer. The locally compact and structured point clouds generated by SPD precisely reveal the structural characteristics of the 3D shape in local patches, which enables us to predict highly detailed geometries. Moreover, since SPD is a general operation that is not limited to completion, we explore its applications in other generative tasks, including point cloud auto-encoding, generation, single image reconstruction, and upsampling. Our experimental results outperform state-of-the-art methods under widely used benchmarks.


title: Mixture-of-Experts with Expert Choice Routing

abstract: Sparsely-activated Mixture-of-experts (MoE) models allow the number of parameters to greatly increase while keeping the amount of computation for a given token or a given sample unchanged. However, a poor expert routing strategy (e.g. one resulting in load imbalance) can cause certain experts to be under-trained, leading to an expert being under or over-specialized. Prior work allocates a fixed number of experts to each token using a top-k function regardless of the relative importance of different tokens. To address this, we propose a heterogeneous mixture-of-experts employing an expert choice method. Instead of letting tokens select the top-k experts, we have experts selecting the top-k tokens. As a result, each token can be routed to a variable number of experts and each expert can have a fixed bucket size. We systematically study pre-training speedups using the same computational resources of the Switch Transformer top-1 and GShard top-2 gating of prior work and find that our method improves training convergence time by more than 2x. For the same computational cost, our method demonstrates higher performance in fine-tuning 11 selected tasks in the GLUE and SuperGLUE benchmarks. For a smaller activation cost, our method outperforms the T5 dense model in 7 out of the 11 tasks.


title: Exact bistability and time pseudo-crystallization of driven-dissipative fermionic lattices

abstract: The existence of bistability in quantum optical systems remains a intensely debated open question beyond the mean-field approximation. Quantum fluctuations are finite-size corrections to the mean-field approximation used because the full exact solution is unobtainable. Usually, quantum fluctuations destroy the bistability present on the mean-field level. Here, by identifying and using exact modulated semi-local dynamical symmetries in a certain quantum optical models of driven-dissipative fermionic chains we exactly prove bistability in precisely the quantum fluctuations. Surprisingly, rather than destroying bistability, the quantum fluctuations themselves exhibit bistability, even though it is absent on the mean-field level for our systems. Moreover, the models studied acquire additional thermodynamic dynamical symmetries that imply persistent periodic oscillations in the quantum fluctuations, constituting pseudo-variants of boundary time crystals. Physically, these emergent operators correspond to finite-frequency and finite-momentum semi-local Goldstone modes. Our work therefore provides to the best of our knowledge the first example of a provably bistable quantum optical system.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph


title: Comparison of Electron Capture Rates in the N=50 Region using 1D Simulations of Core-collapse Supernovae

abstract: Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parameterized N=50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.


title: The ALPs from the Top: Searching for long lived axion-like particles from exotic top decays

abstract: We propose a search for long lived axion-like particles (ALPs) in exotic top decays. Flavour-violating ALPs appear as low energy effective theories for various new physics scenarios such as t-channel dark sectors or Froggatt-Nielsen models. In this case the top quark may decay to an ALP and an up- or charm-quark. For masses in the few GeV range, the ALP is long lived across most of the viable parameter space, suggesting a dedicated search. We propose to search for these long lived ALPs in $t\bar{t}$ events, using one top quark as a trigger. We focus on ALPs decaying in the hadronic calorimeter, and show that the ratio of energy deposits in the electromagnetic and hadronic calorimeters as well as track vetoes can efficiently suppress Standard Model backgrounds. Our proposed search can probe exotic top branching ratios smaller than $10^{-4}$ with a conservative strategy at the upcoming LHC run, and potentially below the $10^{-7}$ level with more advanced methods. Finally we also show that measurements of single top production probe these branching ratios in the very short and very long lifetime limit at the $10^{-3}$ level.


title: Quantum Optimization of Maximum Independent Set using Rydberg Atom Arrays

abstract: Realizing quantum speedup for practically relevant, computationally hard problems is a central challenge in quantum information science. Using Rydberg atom arrays with up to 289 qubits in two spatial dimensions, we experimentally investigate quantum algorithms for solving the Maximum Independent Set problem. We use a hardware-efficient encoding associated with Rydberg blockade, realize closed-loop optimization to test several variational algorithms, and subsequently apply them to systematically explore a class of graphs with programmable connectivity. We find the problem hardness is controlled by the solution degeneracy and number of local minima, and experimentally benchmark the quantum algorithm's performance against classical simulated annealing. On the hardest graphs, we observe a superlinear quantum speedup in finding exact solutions in the deep circuit regime and analyze its origins.

categories: quant-ph cond-mat.dis-nn cond-mat.quant-gas physics.atom-ph


title: Requirements on common solutions to the LSND and MiniBooNE excesses: a post-MicroBooNE study

abstract: The strong statistical significance of an observed electron-like event excess in the MiniBooNE (MB) experiment, along with an earlier similar excess seen in the Liquid Scintillator Neutrino Detector (LSND), when interpreted in conjunction with recent MicroBooNE results may have brought us to the cusp of new physics discoveries. This has led to many attempts to understand these observations, both for each experiment individually and in conjunction, via physics beyond the Standard Model (SM). We provide an overview of the current situation, and discuss three major categories under which the many proposals for new physics fall. The possibility that the same new, non-oscillation physics explains both anomalies leads to new restrictions and requirements. An important class of such common solutions, which we focus on in this work, consists of a heavy ${\cal O}$(MeV$-$sub-GeV) sterile neutral fermion produced in the detectors, (via up-scattering of the incoming muon neutrinos), and subsequently decaying to photons or $e^+e^-$ pairs which mimic the observed signals. Such solutions are subject to strong demands from a) cross section requirements which would yield a sufficient number of total events in both LSND and MB, b) requirements imposed by the measured energy and angular distributions in both experiments and finally, c) consistency and compatibility of the new physics model and its particle content with other bounds from a diverse swathe of particle physics experiments. We find that these criteria often pull proposed solutions in different directions, and stringently limit the viable set of proposals which could resolve both anomalies. Our conclusions are relevant for both the general search for new physics and for the ongoing observations and analyses of the MicroBooNE experiment.


title: Explaining, Evaluating and Enhancing Neural Networks' Learned Representations

abstract: Most efforts in interpretability in deep learning have focused on (1) extracting explanations of a specific downstream task in relation to the input features and (2) imposing constraints on the model, often at the expense of predictive performance. New advances in (unsupervised) representation learning and transfer learning, however, raise the need for an explanatory framework for networks that are trained without a specific downstream task. We address these challenges by showing how explainability can be an aid, rather than an obstacle, towards better and more efficient representations. Specifically, we propose a natural aggregation method generalizing attribution maps between any two (convolutional) layers of a neural network. Additionally, we employ such attributions to define two novel scores for evaluating the informativeness and the disentanglement of latent embeddings. Extensive experiments show that the proposed scores do correlate with the desired properties. We also confirm and extend previously known results concerning the independence of some common saliency strategies from the model parameters. Finally, we show that adopting our proposed scores as constraints during the training of a representation learning task improves the downstream performance of the model.


title: Ephemeral Learning -- Augmenting Triggers with Online-Trained Normalizing Flows

abstract: The large data rates at the LHC require an online trigger system to select relevant collisions. Rather than compressing individual events, we propose to compress an entire data set at once. We use a normalizing flow as a deep generative model to learn the probability density of the data online. The events are then represented by the generative neural network and can be inspected offline for anomalies or used for other analysis purposes. We demonstrate our new approach for a toy model and a correlation-enhanced bump hunt.


title: Using Artificial Intelligence and real galaxy images to constrain parameters in galaxy formation simulations

abstract: Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner.


title: Exploring the physical properties of lensed star-forming clumps at $2\lesssim z \lesssim6$

abstract: We study the physical properties (size, stellar mass, luminosity, star formation rate) and scaling relations for a sample of 166 star-forming clumps with redshift $z \sim 2-6.2$. They are magnified by the Hubble Frontier Field galaxy cluster MACS~J0416 and have robust lensing magnification ($2\lesssim \mu \lesssim 82$) computed by using our high-precision lens model, based on 182 multiple images. Our sample extends by $\sim 3$ times the number of spectroscopically-confirmed lensed clumps at $z \gtrsim 2$. We identify clumps in ultraviolet continuum images and find that, whenever the effective spatial resolution (enhanced by gravitational lensing) increases, they fragment into smaller entities, likely reflecting the hierarchically-organized nature of star formation. Kpc-scale clumps, most commonly observed in field, are not found in our sample. The physical properties of our sample extend the parameter space typically probed by $z \gtrsim 1$ field observations and simulations, by populating the low mass (M$_\star \lesssim 10^7$ M$_\odot$), low star formation rate (SFR $\lesssim 0.5$ M$_\odot$ yr$^{-1}$), and small size (R$_\mathrm{eff} \lesssim 100$ pc) regime. The new domain probed by our study approaches the regime of compact stellar complexes and star clusters. In the mass-size plane, our sample spans the region between galaxies and globular clusters, with a few clumps in the region populated by young star clusters and globular-clusters. For the bulk of our sample, we measure star-formation rates which are higher than those observed locally in compact stellar systems, indicating different conditions for star formation at high redshift than in the local Universe.


title: The Very Knotty Lenser: exploring the role of regularization in source and potential reconstructions using Gaussian Process Regression

abstract: Reconstructing lens potentials and lensed sources can easily become an underconstrained problem, even when the degrees of freedom are low, due to degeneracies, particularly when potential perturbations superimposed on a smooth lens are included. Regularization has traditionally been used to constrain the solutions where the data failed to do so, e.g. in unlensed parts of the source. In this exploratory work, we go beyond the usual choices of regularization and adopt observationally motivated priors for the source brightness. We also perform a similar comparison when reconstructing lens potential perturbations, which are assumed to be stationary, i.e. permeate the entire field of view. We find that physically motivated priors lead to lower residuals, avoid overfitting, and are decisively preferred within a Bayesian quantitative framework in all the examples considered. For the perturbations, choosing the wrong regularization can have a detrimental effect that even high-quality data cannot correct for, while using a purely smooth lens model can absorb them to a very high degree and lead to biased solutions. Finally, our new implementation of the semi-linear inversion technique provides the first quantitative framework for measuring degeneracies between the source and the potential perturbations.


title: Hints of Early Dark Energy in Planck, SPT, and ACT data: new physics or systematics?

abstract: We investigate constraints on early dark energy (EDE) using ACT DR4, SPT-3G 2018, Planck polarization, and restricted Planck temperature data (at $\ell<650$), finding a $3.3\sigma$ preference ($\Delta\chi^2=-16.2$ for 3 additional degrees of freedom) for EDE over $\Lambda$CDM. The EDE contributes a maximum fractional energy density of $f_{\rm EDE}(z_c)=0.163^{+0.047}_{-0.04}$ at a redshift $z_c=3357\pm200$ and leads to a CMB inferred value of the Hubble constant $H_0=74.2^{+1.9}_{-2.1}$ km/s/Mpc. We find that Planck and ACT DR4 data provide the majority of the improvement in $\chi^2$, and that the inclusion of SPT-3G pulls the posterior of $f_{\rm EDE}(z_c)$ away from $\Lambda$CDM. This is the first time that a moderate preference for EDE has been reported for these three combined CMB data sets. We find that including measurements of supernovae luminosity distances and the baryon acoustic oscillation standard ruler only minimally affects the preference ($3.0\sigma$), while measurements that probe the clustering of matter at late times - the lensing potential power spectrum from Planck and $f \sigma_8$ from BOSS - decrease the significance of the preference to 2.6$\sigma$. Conversely, adding a prior on the $H_0$ value as reported by the SH0ES collaboration increases the preference to the $4-5\sigma$ level. In the absence of this prior, the inclusion of Planck TT data at $\ell>1300$ reduces the preference from $3.0\sigma$ to $2.3\sigma$ and the constraint on $f_{\rm EDE}(z_c)$ becomes compatible with $\Lambda$CDM at $1\sigma$. We explore whether systematic errors in the Planck polarization data may affect our conclusions and find that changing the TE polarization efficiencies significantly reduces the Planck preference for EDE. More work will be necessary to establish whether these hints for EDE within CMB data alone are the sole results of systematic errors or an opening to new physics.


title: Efficient Gravitational Wave Template Bank Generation with Differentiable Waveforms

abstract: The most sensitive search pipelines for gravitational waves from compact binary mergers use matched filters to extract signals from the noisy data stream coming from gravitational wave detectors. Matched-filter searches require banks of template waveforms covering the physical parameter space of the binary system. Unfortunately, template bank construction can be a time-consuming task. Here we present a new method for efficiently generating template banks that utilizes automatic differentiation to calculate the parameter space metric. Principally, we demonstrate that automatic differentiation enables accurate computation of the metric for waveforms currently used in search pipelines, whilst being computationally cheap. Additionally, by combining random template placement and a Monte Carlo method for evaluating the fraction of the parameter space that is currently covered, we show that search-ready template banks for frequency-domain waveforms can be rapidly generated. Finally, we argue that differentiable waveforms offer a pathway to accelerating stochastic placement algorithms. We implement all our methods into an easy-to-use Python package based on the jax framework, diffbank, to allow the community to easily take advantage of differentiable waveforms for future searches.


title: Synthetic Disinformation Attacks on Automated Fact Verification Systems

abstract: Automated fact-checking is a needed technology to curtail the spread of online misinformation. One current framework for such solutions proposes to verify claims by retrieving supporting or refuting evidence from related textual sources. However, the realistic use cases for fact-checkers will require verifying claims against evidence sources that could be affected by the same misinformation. Furthermore, the development of modern NLP tools that can produce coherent, fabricated content would allow malicious actors to systematically generate adversarial disinformation for fact-checkers. In this work, we explore the sensitivity of automated fact-checkers to synthetic adversarial evidence in two simulated settings: AdversarialAddition, where we fabricate documents and add them to the evidence repository available to the fact-checking system, and AdversarialModification, where existing evidence source documents in the repository are automatically altered. Our study across multiple models on three benchmarks demonstrates that these systems suffer significant performance drops against these attacks. Finally, we discuss the growing threat of modern NLG systems as generators of disinformation in the context of the challenges they pose to automated fact-checkers.


title: Probing neutron stars with the full pre-merger and post-merger gravitational wave signal from binary coalescences

abstract: The gravitational wave signal emitted during the coalescence of two neutron stars carries information about the stars' internal structure. During the long inspiral phase the main matter observable is the tidal interaction between the binary components, an effect that can be parametrically modeled with compact-binary solutions to General Relativity. After the binary merger the main observable is frequency modes of the remnant, most commonly giving rise to a short-duration signal accessible only through numerical simulations. The complicated morphology and the decreasing detector sensitivity in the relevant frequencies currently hinder detection of the post-merger signal and motivate separate analyses for the pre-merger and post-merger data. However, planned and ongoing detector improvements could soon put the post-merger signal within reach. In this study we target the whole pre-merger and post-merger signal without an artificial separation at the binary merger. We construct a hybrid analysis that models the inspiral with templates based on analytical calculations and calibrated to numerical relativity and the post-merger signal with a flexible morphology-independent analysis. Applying this analysis to GW170817 we find, as expected, that the post-merger signal remains undetected. We further study simulated signals and find that we can reconstruct the full signal and simultaneously estimate both the pre-merger tidal deformation and the post-merger signal frequency content. Our analysis allows us to study neutron star physics using all the data available and directly test the pre-merger and post-merger signal for consistency thus probing effects such as the onset of the hadron-quark phase transition.


title: A Bayesian Hierarchical Time Series Model for Reconstructing Hydroclimate from Multiple Proxies

abstract: We propose a Bayesian hierarchical model which produces probabilistic reconstructions of hydroclimatic variability in Queensland Australia. The model provides a standardised approach to hydroclimate reconstruction using multiple palaeoclimate proxy records derived from natural archives such as speleothems, ice cores and tree rings. The method combines time-series modelling with inverse prediction to quantify the relationships between a given hydroclimate index and relevant proxies over an instrumental period and subsequently reconstruct the hydroclimate back through time. We present case studies for Brisbane and Fitzroy catchments focusing on two hydroclimate indices, the Rainfall Index (RFI) and the Standardised Precipitation-Evapotranspiration Index (SPEI). The probabilistic nature of the reconstructions allows us to estimate the probability that a hydroclimate index in any reconstruction year was lower (higher) than the minimum (maximum) value observed over the instrumental period. In Brisbane, the RFI is unlikely (probabilities < 20%) to have exhibited extremes beyond the minimum/maximum values observed between 1889 and 2017. However, in Fitzroy there are several years during the reconstruction period where the RFI is likely (> 50% probability) to have exhibited behaviour beyond the minimum/maximum of what has been observed. For SPEI, the probability of observing such extremes since the end of the instrumental period in 1889 doesn't exceed 50% in any reconstruction year in Brisbane or Fitzroy.


title: Orbits of actions of group superschemes

abstract: Working over an algebraically closed field $\Bbbk$, we prove that all orbits of a left action of an algebraic group superscheme $G$ on a superscheme $X$ of finite type are locally closed. Moreover, such an orbit $Gx$, where $x$ is a $\Bbbk$-point of $X$, is closed if and only if $G_{ev}x$ is closed in $X_{ev}$, or equivalently, if and only if $G_{res}x$ is closed in $X_{res}$. Here $G_{ev}$ is the largest purely even group super-subscheme of $G$ and $G_{res}$ is $G_{ev}$ regarded as a group scheme. Similarly, $X_{ev}$ is the largest purely even super-subscheme of $X$ and $X_{res}$ is $X_{ev}$ regarded as a scheme. We also prove that $\mathrm{sdim}(Gx)=\mathrm{sdim}(G)-\mathrm{sdim}(G_x)$, where $G_x$ is the stabilizer of $x$.


title: SMEFT interpretation of $\Delta$F = 2 transitions

abstract: A model-independent anatomy of $\Delta F= 2$ transitions in the context of the Weak Effective Theory (WET) below the electroweak scale (EW) and the Standard Model Effective Field Theory (SMEFT) above the EW scale is discussed. Two master formulae for the BSM contribution of the mixing amplitude $M_{12}$, in terms of Wilson coefficients are presented. The coefficients entering these formulae contain all the information below the EW scale and the NP scale $\Lambda$, respectively. The renormalization group evolution from the top-quark Yukawa coupling has the largest impact on the result. The obtained expressions depend on whether the down-basis or the up-basis for SMEFT operators is considered.


title: Mid and Far-Infrared Color-Color Relations within Local Galaxies

abstract: We present an extensive archival analysis of a sample of local galaxies, combining multi-wavelength data from GALEX, Spitzer and Herschel to investigate "blue-side" mid-infrared (MIR) and "red-side" far-infrared (FIR) color-color correlations within the observed infrared spectral energy distributions (IR SEDs). Our sample largely consists of the KINGFISH galaxies, with the important addition of a select few including NGC5236 (M83) and NGC4449. With data from the far-ultraviolet FUV (${\sim}0.15$ $\mu$m) through 500 $\mu$m convolved to common angular resolution, we measure photometry of $kpc$-scale star-forming regions 36$"\times$36$"$ in size. Star formation rates (SFRs), stellar masses and metallicity distributions are derived throughout our sample. Focusing on the $f_{70}/f_{500}$ "FIR" and $f_{8}/f_{24}$ "MIR" flux density ratios (colors), we find that a sub-sample of galaxies demonstrate a strong IR color-color correlation within their star-forming regions, while others demonstrate uncorrelated colors. This division is driven by two main effects: 1) the local strength of star formation (SF) and 2) the metal content of the interstellar medium (ISM). Galaxies uniformly dominated by high surface densities of SF (e.g. NGC5236) demonstrate strong IR color-color correlations, while galaxies that exhibit lower levels of SF and mixed environments (e.g. NGC5457) demonstrate weaker or no correlation--explained by the increasing effect of varying ISM heating and metal content on the IR colors, specifically in the MIR. We find large dispersion in the SFR-$L_{8}$ (8 $\mu$m luminosity) relation that is traced by the metallicity distributions, consistent with extant studies, highlighting its problematic use as a SFR indicator across diverse systems/samples.


title: A mixed-integer programming model for identifying intuitive ambulance dispatching policies

abstract: Markov decision process models and algorithms can be used to identify optimal policies for dispatching ambulances to spatially distributed customers, where the optimal policies indicate the ambulance to dispatch to each customer type in each state. Since the optimal solutions are dependent on Markov state variables, they may not always correspond to a simple set of rules when implementing the policies in practice. Restricted policies that conform to a priority list for each type of customer may be desirable for use in practice, since such policies are transparent, explainable, and easy to implement. A priority list policy is an ordered list of ambulances that indicates the preferred order to dispatch the ambulances to a customer type subject to ambulance availability. This paper proposes a constrained Markov decision process model for identifying optimal priority list policies that is formulated as a mixed integer programming model, does not extend the Markov state space, and can be solved using standard algorithms. A series of computational examples illustrate the benefit of intuitive policies. The optimal mixed integer programming solutions to the computational examples have objective function values that are close to those of the unrestricted model and are superior to those of heuristics.


title: A Molecular Prior Distribution for Bayesian Inference Based on Wilson Statistics

abstract: Background and Objective: Wilson statistics describe well the power spectrum of proteins at high frequencies. Therefore, it has found several applications in structural biology, e.g., it is the basis for sharpening steps used in cryogenic electron microscopy (cryo-EM). A recent paper gave the first rigorous proof of Wilson statistics based on a formalism of Wilson's original argument. This new analysis also leads to statistical estimates of the scattering potential of proteins that reveal a correlation between neighboring Fourier coefficients. Here we exploit these estimates to craft a novel prior that can be used for Bayesian inference of molecular structures. Methods: We describe the properties of the prior and the computation of its hyperparameters. We then evaluate the prior on two synthetic linear inverse problems, and compare against a popular prior in cryo-EM reconstruction at a range of SNRs. Results: We show that the new prior effectively suppresses noise and fills-in low SNR regions in the spectral domain. Furthermore, it improves the resolution of estimates on the problems considered for a wide range of SNR and produces Fourier Shell Correlation curves that are insensitive to masking effects. Conclusions: We analyze the assumptions in the model, discuss relations to other regularization strategies, and postulate on potential implications for structure determination in cryo-EM.


title: Black-box Node Injection Attack for Graph Neural Networks

abstract: Graph Neural Networks (GNNs) have drawn significant attentions over the years and been broadly applied to vital fields that require high security standard such as product recommendation and traffic forecasting. Under such scenarios, exploiting GNN's vulnerabilities and further downgrade its classification performance become highly incentive for adversaries. Previous attackers mainly focus on structural perturbations of existing graphs. Although they deliver promising results, the actual implementation needs capability of manipulating the graph connectivity, which is impractical in some circumstances. In this work, we study the possibility of injecting nodes to evade the victim GNN model, and unlike previous related works with white-box setting, we significantly restrict the amount of accessible knowledge and explore the black-box setting. Specifically, we model the node injection attack as a Markov decision process and propose GA2C, a graph reinforcement learning framework in the fashion of advantage actor critic, to generate realistic features for injected nodes and seamlessly merge them into the original graph following the same topology characteristics. Through our extensive experiments on multiple acknowledged benchmark datasets, we demonstrate the superior performance of our proposed GA2C over existing state-of-the-art methods. The data and source code are publicly accessible at: https://github.com/jumxglhf/GA2C.


title: Higher dimensional CFTs as 2D conformally-equivariant topological field theories

abstract: Two and three-point functions of primary fields in four dimensional CFT have a simple space-time dependences factored out from the combinatoric structure which enumerates the fields and gives their couplings. This has led to the formulation of two dimensional topological field theories with $so(4,2)$ equivariance which are conjectured to be equivalent to higher dimensional conformal field theories. We review this CFT4/TFT2 construction in the simplest possible setting of a free scalar field, which gives an algebraic construction of the correlators in terms of an infinite dimensional $so(4,2)$ equivariant algebra with finite dimensional subspaces at fixed scaling dimension. Crossing symmetry of the CFT4 is related to associativity of the algebra. This construction is then extended to describe perturbative CFT4, by making use of deformed co-products. Motivated by the Wilson-Fisher CFT we outline the construction of U(so($d$,2)) equivariant TFT2 for non-integer $d$, in terms of diagram algebras and their representations.


title: Counterfactual Analysis of the Impact of the IMF Program on Child Poverty in the Global-South Region using Causal-Graphical Normalizing Flows

abstract: This work demonstrates the application of a particular branch of causal inference and deep learning models: \emph{causal-Graphical Normalizing Flows (c-GNFs)}. In a recent contribution, scholars showed that normalizing flows carry certain properties, making them particularly suitable for causal and counterfactual analysis. However, c-GNFs have only been tested in a simulated data setting and no contribution to date have evaluated the application of c-GNFs on large-scale real-world data. Focusing on the \emph{AI for social good}, our study provides a counterfactual analysis of the impact of the International Monetary Fund (IMF) program on child poverty using c-GNFs. The analysis relies on a large-scale real-world observational data: 1,941,734 children under the age of 18, cared for by 567,344 families residing in the 67 countries from the Global-South. While the primary objective of the IMF is to support governments in achieving economic stability, our results find that an IMF program reduces child poverty as a positive side-effect by about 1.2$\pm$0.24 degree (`0' equals no poverty and `7' is maximum poverty). Thus, our article shows how c-GNFs further the use of deep learning and causal inference in AI for social good. It shows how learning algorithms can be used for addressing the untapped potential for a significant social impact through counterfactual inference at population level (ACE), sub-population level (CACE), and individual level (ICE). In contrast to most works that model ACE or CACE but not ICE, c-GNFs enable personalization using \emph{`The First Law of Causal Inference'}.


title: Smooth time optimal trajectory generation for drones

abstract: In this paper, we address a minimum-time steering problem for a drone modeled as point mass with bounded acceleration, across a set of desired waypoints in the presence of gravity. We first provide a method to solve for the minimum-time control input that will steer the point mass between two waypoints based on a continuous-time problem formulation which we address by using Pontryagin's Minimum Principle. Subsequently, we solve for the time-optimal trajectory across the given set of waypoints by discretizing in the time domain and formulating the minimum-time problem as a nonlinear program (NLP). The velocities at each waypoint obtained from solving the NLP in the discretized domain are then used as boundary conditions to extend our two-point solution across those multiple waypoints. We apply this planning methodology to execute a surveying task that minimizes the time taken to completely explore a target area or volume. Numerical simulations and theoretical analyses of this new planning methodology are presented. The results from our approach are also compared to traditional polynomial trajectories like minimum snap planning.


title: Information Decomposition Diagrams Applied beyond Shannon Entropy: A Generalization of Hu's Theorem

abstract: In information theory, one major goal is to find useful functions that summarize the amount of information contained in the interaction of several random variables. Specifically, one can ask how the classical Shannon entropy, mutual information, and higher interaction information relate to each other. This is answered by Hu's theorem, which is widely known in the form of information diagrams: it relates shapes in a Venn diagram to information functions, thus establishing a bridge from set theory to information theory. In this work, we view random variables together with the joint operation as a monoid that acts by conditioning on information functions, and entropy as a function satisfying the chain rule of information. This abstract viewpoint allows to prove a generalization of Hu's theorem. It applies to Shannon and Tsallis entropy, (Tsallis) Kullback-Leibler Divergence, cross-entropy, Kolmogorov complexity, submodular information functions, and the generalization error in machine learning. Our result implies for Chaitin's Kolmogorov complexity that the interaction complexities of all degrees are in expectation close to Shannon interaction information. For well-behaved probability distributions on increasing sequence lengths, this shows that the per-bit expected interaction complexity and information asymptotically coincide, thus showing a strong bridge between algorithmic and classical information theory.


title: Algebraic cycles and functorial lifts from $G_2$ to $\mathrm{PGSp}_6$

abstract: We study instances of Beilinson-Tate conjectures for automorphic representations of $\mathrm{PGSp}_6$ whose Spin $L$-function has a pole at $s=1$. We construct algebraic cycles of codimension three in the Siegel-Shimura variety of dimension six and we relate its regulator to the residue at $s=1$ of the $L$-function of certain cuspidal forms of $\mathrm{PGSp}_6$. Using the exceptional theta correspondence between the split group of type $G_2$ and $\mathrm{PGSp}_6$ and assuming the non-vanishing of a certain archimedean integral, this allows us to confirm a conjecture of Gross and Savin on rank $7$ motives of type $G_2$.


title: A new visit to the variable stars in M56 and its colour magnitude diagram structure

abstract: We present a $VI$ CCD photometric study and a membership analysis of the globular cluster M56 (NGC 6779). This produced a decontaminated CMD from field stars, which enabled a better confrontation with theoretical isochrones; Zero-Age Horizontal Branches (ZAHB) and post-ZAHB evolutionary tracks. Post He-flash evolutionary models with a He-core mass of 0.5 M$_\odot$ and envelopes of 0.04 - 0.18 M$_\odot$, cover the complete Horizontal Branch. Models with total mass $\sim$ 0.68 M$_\odot$ explain the RR Lyrae, while those with a mass $\sim$ 0.56 M$_\odot$ and a very subtle envelope, explain the Pop II cepheids with a progenitor in the blue tail of the HB. %Using the mean $(V-I)$ at minimum color curve for the RRab stars, we derived a value for the reddening of $E(B-V)$ = 0.26 $\pm$ 0.02. Based on the Fourier decomposition of the $V$ light curve of a single cluster member RRc star, we determined a metallicity of $\rm [Fe/H]_{ ZW}$ = -1.96 $\pm$ 0.09. Several independent distance determination approaches lead to a mean distance to M56 of $\big< d \big>$ = 9.4 $\pm$ 0.4 kpc. Finally, we report 5 new variables: 1 SX Phe, 3 EB, and 1 RRc.


title: Real space representation of the quasiparticle self-consistent $GW$ self-energy and its application to defect calculations

abstract: The quasiparticle self-consistent QS$GW$ approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. We show that this representation is relatively short-ranged. This offers new opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-eneregy can be constructed from those of the host and a smaller defect containing cell. The self-energy of the periodic host can be constructed simply from a $GW$ calculation for the primitive cell. We show for the case of the As$_\mathrm{Ga}$ in GaAs that the defect part can already be well represented by a minimal 8 atom cell and allows us to construct the self-energy for a 64 cell in good agreement with direct QS$GW$ calculations for the large cell. Using this approach to an even larger 216 atom cell shows the defect band approaches an isolated defect level. The calculations also allow to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.


title: The theta invariants and the volume function on arithmetic varieties

abstract: We introduce a new arithmetic invariant for hermitian line bundles on an arithmetic variety. We use this invariant to measure the variation of the volume function with respect to the metric. The main result of this paper is a generalized Hodge index theorem on arithmetic toric varieties.


title: Provably Private Distributed Averaging Consensus: An Information-Theoretic Approach

abstract: In this work, we focus on solving a decentralized consensus problem in a private manner. Specifically, we consider a setting in which a group of nodes, connected through a network, aim at computing the mean of their local values without revealing those values to each other. The distributed consensus problem is a classic problem that has been extensively studied and its convergence characteristics are well-known. Alas, state-of-the-art consensus methods build on the idea of exchanging local information with neighboring nodes which leaks information about the users' local values. We propose an algorithmic framework that is capable of achieving the convergence limit and rate of classic consensus algorithms while keeping the users' local values private. The key idea of our proposed method is to carefully design noisy messages that are passed from each node to its neighbors such that the consensus algorithm still converges precisely to the average of local values, while a minimum amount of information about local values is leaked. We formalize this by precisely characterizing the mutual information between the private message of a node and all the messages that another adversary collects over time. We prove that our method is capable of preserving users' privacy for any network without a so-called "generalized leaf", and formalize the trade-off between privacy and convergence time. Unlike many private algorithms, any desired accuracy is achievable by our method, and the required level of privacy only affects the convergence time.


title: States of finite effect algebras

abstract: In this paper we characterize finite effect algebras which have a state. We construct two matrices $A$ and $B$ assigned to a finite effect algebra $E$ and show that if $E$ has a state then rank$A=$ rank$B$.


title: Equivariant Transporter Network

abstract: Transporter Net is a recently proposed framework for pick and place that is able to learn good manipulation policies from a very few expert demonstrations. A key reason why Transporter Net is so sample efficient is that the model incorporates rotational equivariance into the pick module, i.e. the model immediately generalizes learned pick knowledge to objects presented in different orientations. This paper proposes a novel version of Transporter Net that is equivariant to both pick and place orientation. As a result, our model immediately generalizes place knowledge to different place orientations in addition to generalizing pick knowledge as before. Ultimately, our new model is more sample efficient and achieves better pick and place success rates than the baseline Transporter Net model.


title: Density Matrices of Seniority-Zero Geminal Wavefunctions

abstract: Scalar products and density matrix elements of closed-shell pair geminal wavefunctions are evaluated directly in terms of the pair amplitudes, resulting in an analogue of Wick's theorem for fermions or bosons. This expression is in general intractable, but it is shown how it becomes feasible in three distinct ways for Richardson-Gaudin (RG) states, the antisymmetrized geminal power, and the antisymmetrized product of strongly-orthogonal geminals. Dissociation curves for hydrogen chains are computed with off-shell RG states and the antisymmetrized product of interacting geminals. Both are near exact suggesting that the incorrect results observed with ground state RG states are fixable using a different RG state.


title: Controlling two-dimensional chaotic transients with the safety function

abstract: In this work we deal with the H\'enon and the Lozi map for a choice of parameters where they show transient chaos. Orbits close to the chaotic saddle behave chaotically for a while to eventually escape to an external attractor. Traditionally, to prevent such an escape, the partial control technique has been applied. This method stands out for considering disturbances (noise) affecting the map and for finding a special region of the phase space, called the safe set, where the control required to sustain the orbits is small. However, in this work we will apply a new approach of the partial control method that has been recently developed. This new approach is based on finding a special function called the safety function, which allows to automatically find the minimum control necessary to avoid the escape of the orbits. Furthermore, we will show the strong connection between the safety function and the classical safe set. To illustrate that, we will compute for the first time, safety functions for the two-dimensional H\'enon and Lozi maps, where we also show the strong dependence of this function with the magnitude of disturbances affecting the map, and how this change drastically impacts the controlled orbits.


title: Multiple Ancillary Services Provision by Distributed Energy Resources in Active Distribution Networks

abstract: The electric power system is currently experiencing radical changes stemming from the increasing share of renewable energy resources and the consequent decommissioning of conventional power plants based on synchronous generators. Since the principal providers of ancillary services are being phased out, new flexibility and reserve providers are needed. The proliferation of Distributed Energy Resources (DERs) in modern distribution networks has opened new possibilities for distribution system operators, enabling them to fill the market gap by harnessing the DER flexibility. This paper introduces a novel centralized MPC-based controller that enables the concurrent provision of voltage support, primary and secondary frequency control by adjusting the setpoints of a heterogeneous group of DERs in active distribution grids. The input-multirate control framework is used to accommodate the distinct timescales and provision requirements of each ancillary service and to ensure that the available resources are properly allocated. Furthermore, an efficient way for incorporating network constraints in the formulation is proposed, where network decomposition is applied to a linear power flow formulation together with network reduction. In addition, different timescale dynamics of the employed DERs and their capability curves are included. The performance of the proposed controller is evaluated on several case studies via dynamic simulations of the IEEE 33-bus system.


title: A minimizing problem of a polyharmonic operator with Critical Exponent

abstract: In this work, we study the two following minimization problems for $r \in \mathbb{N}^{*}$, \begin{equation*} \begin{array}{ccc} S_{0,r}(\varphi)=\displaystyle\inf_{u\in H_{0}^{r}(\Omega)\,|u+\varphi\|_{L^{2^{*r}}}=1}\|u\|_{r}^{2}& \textrm{and}& S_{\theta,r}(\varphi)=\displaystyle\inf_{u\in H_{\theta}^{r}(\Omega)\, \|u+\varphi\|_{L^{2^{*r}}}=1}\|u\|_{r}^{2}, \end{array} \end{equation*} where $\Omega \subset \mathbb{R}^{N}, $ $N > 2r$, is a smooth bounded domain, $2^{*r}=\frac{2N}{N-2 r}$, $\varphi\in L^{2^{*r}} (\Omega) \cap C(\Omega)$ and the norm $\|. \|_{r}=\displaystyle{ \int_{\Omega} |(-\Delta)^{\alpha} .|^{2}dx}$ where $ \alpha=\frac{r}{2} $ if $r$ is even and $\|. \|_{r}=\displaystyle{ \int_{\Omega} |\nabla(-\Delta)^{\alpha} . |^{2}dx }$ where $\alpha = \frac{r-1}{2}$ if $r$ is odd. Firstly, we prove that, when $\varphi \not\equiv 0, $ the infimum in $S_{0,r}(\varphi)$ and $S_{\theta,r}(\varphi)$ are attained. Secondly, we show that $ S_{\theta,r}(\varphi)< S_{0,r}(\varphi) $ for a large class of $ \varphi$.


title: A two-phase rank-based algorithm for low-rank matrix completion

abstract: Matrix completion aims to recover an unknown low-rank matrix from a small subset of its entries. In many applications, the rank of the unknown target matrix is known in advance. In this paper, first we revisit a recently proposed rank-based heuristic for "known-rank" matrix completion and establish a condition under which the generated sequence is quasi-Fej\'er convergent to the solution set. Then, by including an acceleration mechanism similar to Nesterov's acceleration, we obtain a new heuristic. Even though the convergence of such heuristic cannot be granted in general, it turns out that it can be very useful as a warm-start phase, providing a suitable estimate for the regularization parameter and a good starting-point, to an accelerated Soft-Impute algorithm. Numerical experiments with both synthetic and real data show that the resulting two-phase rank-based algorithm can recover low-rank matrices, with relatively high precision, faster than other well-established matrix completion algorithms.


title: Optical quantum super-resolution imaging and hypothesis testing

abstract: Estimating the angular separation between two incoherent thermal sources is a challenging task for direct imaging, especially when it is smaller than or comparable to the Rayleigh length. In addition, the task of discriminating whether there are one or two sources followed by detecting the faint emission of a secondary source in the proximity of a much brighter one is in itself a severe challenge for direct imaging. Here, we experimentally demonstrate two tasks for superresolution imaging based on quantum state discrimination and quantum imaging techniques. We show that one can significantly reduce the probability of error for detecting the presence of a weak secondary source, especially when the two sources have small angular separations. In this work, we reduce the experimental complexity down to a single two-mode interferometer: we show that (1) this simple set-up is sufficient for the state discrimination task, and (2) if the two sources are of equal brightness, then this measurement can super-resolve their angular separation, saturating the quantum Cram\'er-Rao bound. By using a collection baseline of 5.3~mm, we resolve the angular separation of two sources that are placed 15~$\mu$m apart at a distance of 1.0~m with an accuracy of $1.7\%$--this is between 2 to 3 orders of magnitudes more accurate than shot-noise limited direct imaging.


title: Blockchain Driven Privacy Preserving Contact Tracing Framework in Pandemics

abstract: Contact tracing has been proven an effective approach to control the virus spread in pandemics like COVID-19 pandemic. As an emerging powerful decentralized technique, blockchain has been explored to ensure data privacy and security in contact tracing processes. However, existing works are mostly high-level designs with no sufficient demonstration and treat blockchain as separate storage system assisting third-party central servers, ignoring the importance and capability of consensus mechanism and incentive mechanism. In this paper, we propose a light-weight and fully third-party free Blockchain-Driven Contact Tracing framework (BDCT) to bridge the gap. In the BDCT framework, RSA encryption based transaction verification method (RSA-TVM) is proposed to ensure contact tracing correctness, which can achieve more than 96\% contact cases recording accuracy even each person has 60\% probability of failing to verify the contact information. Reputation Corrected Delegated Proof of Stake (RC-DPoS) consensus mechanism is proposed together with the incentive mechanism, which can ensure timeliness of reporting contact cases and keep blockchain decentralized. A novel contact tracing simulation environment is created, which considers three different contact scenarios based on population density. The simulation results demonstrate the effectiveness, robustness and attack resistance of RSA-TVM and RC-DPoS in the proposed BDCT.


title: Unsupervised strategies for identifying optimal parameters in Quantum Approximate Optimization Algorithm

abstract: As combinatorial optimization is one of the main quantum computing applications, many methods based on parameterized quantum circuits are being developed. In general, a set of parameters are being tweaked to optimize a cost function out of the quantum circuit output. One of these algorithms, the Quantum Approximate Optimization Algorithm stands out as a promising approach to tackling combinatorial problems. However, finding the appropriate parameters is a difficult task. Although QAOA exhibits concentration properties, they can depend on instances characteristics that may not be easy to identify, but may nonetheless offer useful information to find good parameters. In this work, we study unsupervised Machine Learning approaches for setting these parameters without optimization. We perform clustering with the angle values but also instances encodings (using instance features or the output of a variational graph autoencoder), and compare different approaches. These angle-finding strategies can be used to reduce calls to quantum circuits when leveraging QAOA as a subroutine. We showcase them within Recursive-QAOA up to depth $3$ where the number of QAOA parameters used per iteration is limited to $3$, achieving a median approximation ratio of $0.94$ for MaxCut over $200$ Erd\H{o}s-R\'{e}nyi graphs. We obtain similar performances to the case where we extensively optimize the angles, hence saving numerous circuit calls.


title: Differentially Private Federated Learning via Inexact ADMM with Multiple Local Updates

abstract: Differential privacy (DP) techniques can be applied to the federated learning model to statistically guarantee data privacy against inference attacks to communication among the learning agents. While ensuring strong data privacy, however, the DP techniques hinder achieving a greater learning performance. In this paper we develop a DP inexact alternating direction method of multipliers algorithm with multiple local updates for federated learning, where a sequence of convex subproblems is solved with the objective perturbation by random noises generated from a Laplace distribution. We show that our algorithm provides $\bar{\epsilon}$-DP for every iteration, where $\bar{\epsilon}$ is a privacy budget controlled by the user. We also present convergence analyses of the proposed algorithm. Using MNIST and FEMNIST datasets for the image classification, we demonstrate that our algorithm reduces the testing error by at most $31\%$ compared with the existing DP algorithm, while achieving the same level of data privacy. The numerical experiment also shows that our algorithm converges faster than the existing algorithm.


title: Power-efficient soliton microcombs

abstract: Laser frequency combs are enabling some of the most exciting scientific endeavours in the 21st century, ranging from the development of optical clocks to the calibration of the astronomical spectrographs used for searching Earth-like exoplanets. Today, dissipative Kerr solitons generated in microresonators offer the prospect of attaining frequency combs in miniaturized systems by capitalizing on advances in photonic integration. Most of the applications based on soliton microcombs rely on tuning a continuous-wave laser into a longitudinal mode of a microresonator whose dimensions are engineered to display anomalous dispersion at the pump laser frequency. In this configuration, however, nonlinear physics precludes from attaining dissipative Kerr solitons with high power conversion efficiency, with typical comb powers amounting to ~1% of the available laser power. Here, we demonstrate that this fundamental limitation can be overcome by inducing a controllable frequency shift to a selected cavity resonance. Experimentally, we realize this shift using two linearly coupled anomalous-dispersion microresonators (a photonic molecule), resulting in a coherent dissipative Kerr soliton with a conversion efficiency exceeding 50% and excellent line spacing stability. We describe the physical soliton dynamics in this configuration, and discover the system displays unusual characteristics, such as the possibility to backwards initiate solitons and stable operation with a blue detuned pump laser. By optimizing the microcomb power available on chip, these results facilitate the practical implementation of a scalable integrated photonic architecture for energy-efficient applications.


title: Minimizing travelling waves for the Gross-Pitaevskii equation on $\mathbb{R} \times \mathbb{T}$

abstract: We study the Gross-Pitaevskii equation in dimension two with periodic conditions in one direction, or equivalently on the product space $\mathbb{R} \times \mathbb{T}_L$ where $L > 0$ and $\mathbb{T}_L = \mathbb{R} / L \mathbb{Z}.$ We focus on the variational problem consisting in minimizing the Ginzburg-Landau energy under a fixed momentum constraint. We prove that there exists a threshold value for $L$ below which minimizers are the one-dimensional dark solitons, and above which no minimizer can be one-dimensional.


title: SN 2020jfo: A short plateau Type II supernova from a low mass progenitor

abstract: We present spectroscopic and photometric observations of the Type IIP supernova, SN 2020jfo, in ultraviolet and optical wavelengths. SN 2020jfo occurred in the spiral galaxy M61 (NGC 4303), with eight observed supernovae in the past 100 years. SN 2020jfo exhibited a short plateau lasting < 65 d, and achieved a maximum brightness in V-band of $M_V$ = -17.4 $\pm$ 0.4 mag at about 8.0 $\pm$ 0.5 d since explosion. From the bolometric light curve, we have estimated the mass of $^{56}$Ni synthesised in the explosion to be 0.033 $\pm$ 0.006 $M_\odot$. The observed spectral features are typical for a type IIP supernova except for shallow H$\alpha$ absorption throughout the evolution and the presence of stable $^{58}$Ni feature at 7378 \r{A}, in the nebular phase. Using hydrodynamical modelling in the MESA + STELLA framework, an ejecta mass of ~ 5 $M_\odot$ is estimated. Models also indicate SN 2020jfo could be the result of a Red Super Giant progenitor with $M_{ZAMS}$ ~ 12 $M_\odot$. Bolometric light curve modelling revealed the presence of a secondary radiation source for initial ~ 20 d, which has been attributed to interaction with a circumstellar material of mass ~ 0.2 $M_\odot$, which most likely was ejected due to enhanced mass loss about 20 years prior to the supernova explosion.


title: Milligram mass metrology for quantitative deposition of liquid samples

abstract: Quantitative dispensing of liquids by mass becomes increasingly difficult as sample volume decreases. Whereas ml quantities of liquid can be directly weighed as dispensed on an analytical balance or using pycnometer methods, {\mu}l samples are more challenging due to their milligram-scale mass and rapid evaporation. Here, a method for quantitative determination of the mass of 1 {\mu}l droplets of liquid is proposed. The weighing of a microcapillary pycnometer combined with traceable calibration masses and corrections for sample evaporation allows determination of the droplet mass on the order of 1 mg with relative combined expanded uncertainty of 0.3 %.


title: Functional Optical Coherence Tomography for Intrinsic Signal Optoretinography: Recent Developments and Deployment Challenges

abstract: Intrinsic optical signal (IOS) imaging of the retina, also termed as optoretinography (ORG), promises a noninvasive method for objective assessment of retinal function. By providing unparalleled capability to differentiate individual layers of the retina, functional optical coherence tomography (OCT) has been actively investigated for intrinsic signal ORG measurements. However, clinical deployment of functional OCT for quantitative ORG is still challenging due to the lack of a standardized imaging protocol and the complication of IOS sources and mechanisms. This article aims to summarize recent developments of functional OCT for ORG measurement, OCT intensity- and phase-based IOS processing. Technical challenges and perspectives of quantitative IOS analysis and ORG interpretations are discussed.

categories: physics.med-ph eess.IV physics.bio-ph physics.optics


title: Navier-Stokes equations in the half space with non compatible data

abstract: In this paper we shall consider the Navier-Stokes equations in the half plane with Euler-type initial conditions, i.e. initial conditions which have a non-zero tangential component at the boundary. Under analyticity assumptions for the data, we shall prove that the solution exists for a small time independent of the viscosity. The solution is constructed through a composite asymptotic expansion involving the solutions of the Euler and Prandtl equations, plus an error term. The norm of the error goes to zero with the square root of the viscosity. The Prandtl solution contains a singular term, which influences the regularity of the error. The error term is written as the sum of a first order Euler correction, a first order Prandtl correction, and a further error term. The use of an analytic setting is mainly due to the Prandtl equation.


title: Harmonic Conjugation in Harmonic Matroids

abstract: We study a generalization of the concept of harmonic conjugation from projective geometry and full algebraic matroids to a larger class of matroids called \emph{harmonic matroids}. We use harmonic conjugation to construct a projective plane of prime order in harmonic matroids without using the axioms of projective geometry. As a particular case we have a combinatorial construction of a projective plane of prime order in full algebraic matroids.


title: Deep Vision-Inspired Bubble Dynamics on Hybrid Nanowires with Dual Wettability

abstract: The boiling efficacy is intrinsically tethered to trade-offs between the desire for bubble nucleation and necessity of vapor removal. The solution to these competing demands requires the separation of bubble activity and liquid delivery, often achieved through surface engineering. In this study, we independently engineer bubble nucleation and departure mechanisms through the design of heterogeneous and segmented nanowires with dual wettability with the aim of pushing the limit of structure-enhanced boiling heat transfer performances. The demonstration of separating liquid and vapor pathways outperforms state-of-the-art hierarchical nanowires, in particular, at low heat flux regimes while maintaining equal performances at high heat fluxes. A machine vision-based framework realizes the autonomous curation and extraction of hidden big data along with bubble dynamics. The combined efforts of materials design, deep learning techniques, and data-driven approach shed light on the mechanistic relationship between vapor/liquid pathways, bubble statistics, and phase change performance.


title: Uniting Control and Data Parallelism: Towards Scalable Memory-Driven Dynamic Graph Processing

abstract: Control parallelism and data parallelism is mostly reasoned and optimized as separate functions. Because of this, workloads that are irregular, fine-grain and dynamic such as dynamic graph processing become very hard to scale. An experimental research approach to computer architecture that synthesizes prior techniques of parallel computing along with new innovations is proposed in this paper. We establish the background and motivation of the research undertaking and provide a detailed description of the proposed omputing system that is highly parallel non-von Neumann, memory-centric and memory-driven. We also present a message-driven (or even-driven) programming model called "diffusive computation" and provide insights into its properties using SSSP and Triangle Counting problems as examples.


title: Optical and Near-Infrared Excesses are Correlated in T Tauri Stars

abstract: Accretion is one of the defining characteristics of classical T Tauri stars, fueled by the presence of a circumstellar disk comprised of dust and gas. Accretion produces a UV and optical excess, while re-radiated emission at the inner edge of the dust component of the disk produces a near-infrared (NIR) excess. The interplay between stars and their disks helps regulate protoplanetary disk evolution and dispersal, which is key to a full understanding of planet formation. To investigate the relations between NIR excess and optical excess in both single and binary stars, we used an archival sample of spectroscopically characterized members of the Taurus star-forming region ($\tau \sim$ 1-2 Myr) with measured luminosities, spectral types, and optical veiling. We combined the archival sample with 2MASS and WISE NIR photometry and high-resolution imaging surveys. We found that NIR and optical excesses are correlated in multiple NIR photometric bands, suggesting that they are closely related, likely because more massive disks have higher inner dust disk walls and are also associated with higher accretion rates. We also found that multiplicity has no impact on accretion or inner disk properties in a sample with a wide range of separations, but the sample was too small to specifically investigate close binaries, where the effects of multiplicity on disk properties should be most significant.


title: Partitioning Dense Graphs with Hardware Accelerators

abstract: Graph partitioning is a fundamental combinatorial optimization problem that attracts a lot of attention from theoreticians and practitioners due to its broad applications. From multilevel graph partitioning to more general-purpose optimization solvers such as Gurobi and CPLEX, a wide range of approaches have been developed. Limitations of these approaches are important to study in order to break the computational optimization barriers of this problem. As we approach the limits of Moore's law, there is now a need to explore ways of solving such problems with special-purpose hardware such as quantum computers or quantum-inspired accelerators. In this work, we experiment with solving the graph partitioning on the Fujitsu Digital Annealer (a special-purpose hardware designed for solving combinatorial optimization problems) and compare it with the existing top solvers. We demonstrate limitations of existing solvers on many dense graphs as well as those of the Digital Annealer on sparse graphs which opens an avenue to hybridize these approaches.


title: Parallel inverse-problem solver for time-domain optical tomography with perfect parallel scaling

abstract: This paper presents an efficient parallel radiative transfer-based inverse-problem solver for time-domain optical tomography. The radiative transfer equation provides a physically accurate model for the transport of photons in biological tissue, but the high computational cost associated with its solution has hindered its use in time-domain optical-tomography and other areas. In this paper this problem is tackled by means of a number of computational and modeling innovations, including 1) A spatial parallel-decomposition strategy with perfect parallel scaling for the forward and inverse problems of optical tomography on parallel computer systems; and, 2) A Multiple Staggered Source method (MSS) that solves the inverse transport problem at a computational cost that is independent of the number of sources employed, and which significantly accelerates the reconstruction of the optical parameters: a six-fold MSS acceleration factor is demonstrated in this paper. Finally, this contribution presents 3) An intuitive derivation of the adjoint-based formulation for evaluation of functional gradients, including the highly-relevant general Fresnel boundary conditions -- thus, in particular, generalizing results previously available for vacuum boundary conditions. Solutions of large and realistic 2D inverse problems are presented in this paper, which were produced on a 256-core computer system. The combined parallel/MSS acceleration approach reduced the required computing times by several orders of magnitude, from months to a few hours.


title: Communication-Efficient Actor-Critic Methods for Homogeneous Markov Games

abstract: Recent success in cooperative multi-agent reinforcement learning (MARL) relies on centralized training and policy sharing. Centralized training eliminates the issue of non-stationarity MARL yet induces large communication costs, and policy sharing is empirically crucial to efficient learning in certain tasks yet lacks theoretical justification. In this paper, we formally characterize a subclass of cooperative Markov games where agents exhibit a certain form of homogeneity such that policy sharing provably incurs no suboptimality. This enables us to develop the first consensus-based decentralized actor-critic method where the consensus update is applied to both the actors and the critics while ensuring convergence. We also develop practical algorithms based on our decentralized actor-critic method to reduce the communication cost during training, while still yielding policies comparable with centralized training.


title: Route Discovery and Capacity of Ad hoc Networks

abstract: Throughput capacity of large ad hoc networks has been shown to scale adversely with the size of network $n$. However the need for the nodes to find or repair routes has not been analyzed in this context. In this paper, we explicitly take route discovery into account and obtain the scaling law for the throughput capacity under general assumptions on the network environment, node behavior, and the quality of route discovery algorithms. We also discuss a number of possible scenarios and show that the need for route discovery may change the scaling for the throughput capacity dramatically.