title: Lower Bounds on Adversarial Robustness from Optimal Transport

abstract: While progress has been made in understanding the robustness of machine learning classifiers to test-time adversaries (evasion attacks), fundamental questions remain unresolved. In this paper, we use optimal transport to characterize the minimum possible loss in an adversarial classification scenario. In this setting, an adversary receives a random labeled example from one of two classes, perturbs the example subject to a neighborhood constraint, and presents the modified example to the classifier. We define an appropriate cost function such that the minimum transportation cost between the distributions of the two classes determines the minimum $0-1$ loss for any classifier. When the classifier comes from a restricted hypothesis class, the optimal transportation cost provides a lower bound. We apply our framework to the case of Gaussian data with norm-bounded adversaries and explicitly show matching bounds for the classification and transport problems as well as the optimality of linear classifiers. We also characterize the sample complexity of learning in this setting, deriving and extending previously known results as a special case. Finally, we use our framework to study the gap between the optimal classification performance possible and that currently achieved by state-of-the-art robustly trained neural networks for datasets of interest, namely, MNIST, Fashion MNIST and CIFAR-10.


title: Search for the doubly charmed baryon $\Xi_{cc}^{+}$

abstract: A search for the doubly charmed baryon $\Xi_{cc}^{+}$ is performed through its decay to the $\Lambda_c^+ K^- \pi^+$ final state, using proton-proton collision data collected with the LHCb detector at centre-of-mass energies of 7, 8 and 13$\mathrm{\,Te\kern -0.1em V}$. The data correspond to a total integrated luminosity of $9\,\mathrm{fb}^{-1}$. No significant signal is observed in the mass range from 3.4 to 3.8$\mathrm{\,Ge\kern -0.1em V}/c^2$. Upper limits are set at $95\%$ credibility level on the ratio of the $\Xi_{cc}^{+}$ production cross-section times the branching fraction to that of $\Lambda_c^+$ and $\Xi_{cc}^{++}$ baryons. The limits are determined as functions of the $\Xi_{cc}^{+}$ mass for different lifetime hypotheses, in the rapidity range from 2.0 to 4.5 and the transverse momentum range from 4 to 15$\mathrm{\,Ge\kern -0.1em V}/c$.


title: Persistence of Excitation in Reproducing Kernel Hilbert Spaces, Positive Limit Sets, and Smooth Manifolds

abstract: This paper studies the relationship between the positive limit sets of continuous semiflows and the newly introduced definition of persistently excited (PE) sets and associated subspaces of reproducing kernel Hilbert (RKH) spaces. It is shown that if the RKH space contains a rich collection of cut-off functions, persistently excited sets are contained as subsets of the positive limit set of the semiflow. The paper demonstrates how the new PE condition can be used to guarantee convergence of function estimates in the RKH space embedding method for adaptive estimation. In particular, the paper is applied to uncertain ODE systems with positive limit sets given by certain types of smooth manifolds, and it establishes convergence of adaptive function estimates over the manifolds.


title: On the Optimality of Treating Inter-Cell Interference as Noise: Downlink Cellular Networks and Uplink-Downlink Duality

abstract: We consider the information-theoretic optimality of treating inter-cell interference as noise (multi-cell TIN) in downlink cellular networks. We focus on scenarios modeled by the Gaussian interfering broadcast channel (IBC), comprising $K$ mutually interfering Gaussian broadcast channels (BCs), each formed by a base station communicating independent messages to an arbitrary number of users. We establish a new power allocation duality between the IBC and its dual interfering multiple access channel (IMAC), which entails that the corresponding generalized degrees-of-freedom regions achieved through multi-cell TIN and power control (TINA regions) for both networks are identical. As by-products of this duality, we obtain an explicit characterization of the IBC TINA region from a previously established characterization of the IMAC TINA region; and identify a multi-cell convex-TIN regime in which the IBC TINA region is a polyhedron (hence convex) without the need for time-sharing. We then identify a smaller multi-cell TIN regime in which the IBC TINA region is optimal and multi-cell TIN achieves the entire capacity region of the IBC, up to a constant gap. This is accomplished by deriving a new genie-aided outer bound for the IBC, that reveals a novel BC-type order that holds amongst users in each constituent BC (or cell) under inter-cell interference, which in turn is not implied by previously known BC-type orders (i.e. degraded, less noisy and more capable orders). The multi-cell TIN regime that we identify for the IBC coincides with a corresponding multi-cell TIN regime previously identified for the IMAC, hence establishing a comprehensive uplink-downlink duality of multi-cell TIN in the GDoF (and approximate capacity) sense.


title: Quantum and semiclassical exceptional points of a linear system of coupled cavities with losses and gain within the Scully-Lamb laser theory

abstract: In the past few decades, many works have been devoted to the study of exceptional points (EPs), i.e., exotic degeneracies of non-Hermitian systems. The usual approach in those studies involves the introduction of a phenomenological effective non-Hermitian Hamiltonian (NHH), where the gain and losses are incorporated as the imaginary frequencies of fields and from which the Hamiltonian EPs (HEPs) are derived. Although this approach can provide valid equations of motion for the fields in the classical limit, its application in the derivation of EPs in the quantum regime is questionable. Recently, a framework [Minganti {\it et al.}, \href{https://doi.org/10.1103/PhysRevA.100.062131}{Phys. Rev. A {\bf 100}, 062131 (2019)}], which allows one to determine quantum EPs from a Liouvillian EP (LEP), rather than from an NHH, has been proposed. Compared to the NHHs, a Liouvillian naturally includes quantum noise effects via quantum-jump terms, thus allowing one to consistently determine its EPs purely in the quantum regime. In this work we study a non-Hermitian system consisting of coupled cavities with unbalanced gain and losses, where the gain is far from saturation, i.e, the system is assumed to be linear. We apply both formalisms, based on an NHH and a Liouvillian within the Scully-Lamb laser theory, to determine and compare the corresponding HEPs and LEPs in the semiclassical and quantum regimes. {Our results indicate that, although the overall spectral properties of the NHH and the corresponding Liouvillian for a given system can differ substantially, their LEPs and HEPs occur for the same combination of system parameters.}


title: Equivariant Cohomotopy implies orientifold tadpole cancellation

abstract: There are fundamental open problems in the precise global nature of RR-field tadpole cancellation conditions in string theory. Moreover, the non-perturbative lift as M5/MO5-anomaly cancellation in M-theory had been based on indirect plausibility arguments,lacking a microscopic underpinning in M-brane charge quantization. We provide a framework for answering these questions, crucial not only for mathematical consistency but also for phenomenological accuracy of string theory, by formulating the M-theory C-field on flat M-orientifolds in the generalized cohomology theory called Equivariant Cohomotopy. This builds on our previous results for smooth but curved spacetimes, showing in that setting that charge quantization in twisted Cohomotopy rigorously implies a list of expected anomaly cancellation conditions. Here we further expand this list by proving that brane charge quantization in unstable equivariant Cohomotopy implies the anomaly cancellation conditions for M-branes and D-branes on flat orbi-orientifolds. For this we (a) use an unstable refinement of the equivariant Hopf-tom Dieck theorem to derive local/twisted tadpole cancellation, and (b) the lift to super-differential cohomology to establish global/untwisted tadpole cancellation. Throughout, we use (c) the unstable Pontrjagin-Thom theorem to identify the brane/O-plane configurations encoded in equivariant Cohomotopy and (d) the Boardman homomorphism to equivariant K-theory to identify Chan-Paton representations of D-brane charge. We find that unstable equivariant Cohomotopy, but not its image in K-theory, distinguishes D-brane charge from the finite set of types of O-plane charges.


title: Box splines, tensor product multiplicities and the volume function

abstract: We study the relationship between the tensor product multiplicities of a compact semisimple Lie algebra $\mathfrak{g}$ and a special function $\mathcal{J}$ associated to $\mathfrak{g}$, called the volume function. The volume function arises in connection with the randomized Horn's problem in random matrix theory and has a related significance in symplectic geometry. Building on box spline deconvolution formulae of Dahmen-Micchelli and De Concini-Procesi-Vergne, we develop new techniques for computing the multiplicities from $\mathcal{J}$, answering a question posed by Coquereaux and Zuber. In particular, we derive an explicit algebraic formula for a large class of Littlewood-Richardson coefficients in terms of $\mathcal{J}$. We also give analogous results for weight multiplicities, and we show a number of further identities relating the tensor product multiplicities, the volume function and the box spline. To illustrate these ideas, we give new proofs of some known theorems.


title: Fine-Grained, Language-Based Access Control for Database-Backed Applications

abstract: Context: Database-backed applications often run queries with more authority than necessary. Since programs can access more data than they legitimately need, flaws in security checks at the application level can enable malicious or buggy code to view or modify data in violation of intended access control policies. Inquiry: Although database management systems provide tools to control access to data, these tools are not well-suited for modern applications which often have many users and consist of many different software components. First, databases are unaware of application users, and creating a new database user for each application user is impractical for applications with many users. Second, different components of the same application may require different authority, which would require creating different database users for different software components. Thus, it is difficult to use existing tools to properly limit the authority an application has when executing queries. For this reason, we consider a new, language-based approach to application-specific database security. Approach: Prior work has addressed the difficulty of running applications with least privilege using capability-based security and software contracts, which we adapt to the setting of database-backed applications. Knowledge: This paper's main contribution is the design and implementation of ShillDB, a language for writing secure database-backed applications. ShillDB enables reasoning about database access at the language level through capabilities, which limit which database tables a program can access, and contracts, which limit what operations a program can perform on those tables. ShillDB contracts are expressed as part of function interfaces, making it easy to specify different access control policies for different components. Contracts act as executable security documentation for ShillDB programs and are enforced by the language runtime. Further, ShillDB provides database access control guarantees independent of (and in addition to) the security mechanisms of the underlying database management system. Grounding: We have implemented a prototype of ShillDB and have used it to implement the backend for a lending library reservation system with contracts for each endpoint to evaluate the performance and usability of ShillDB. Further, we benchmark individual database operations in ShillDB to better understand the language's performance overhead. Importance: Our experience indicates that ShillDB is a practical language for enforcing database access control policies in realistic, multi-user applications and has reasonable performance overhead. ShillDB allows developers to reason about security at the component level, safely compose components, and reuse third-party components with their own application-specific database security policies.


title: Multiplicative functions in short arithmetic progressions

abstract: We study for bounded multiplicative functions $f$ sums of the form \begin{align*} \sum_{\substack{n\leq x \atop n\equiv a\pmod q}}f(n), \end{align*} establishing that their variance over residue classes $a \pmod q$ is small as soon as $q=o(x)$, for almost all moduli $q$, with a nearly power-saving exceptional set of $q$. This improves and generalizes previous results of Hooley on Barban-Davenport-Halberstam-type theorems for such $f$, and moreover our exceptional set is essentially optimal unless one is able to make progress on certain well-known conjectures. We are nevertheless able to prove stronger bounds for the number of the exceptional moduli $q$ in the cases where $q$ is restricted to be either smooth or prime, and conditionally on GRH we show that our variance estimate is valid for every $q$. These results are special cases of a "hybrid result" that works for sums of $f$ over almost all short intervals and arithmetic progressions simultaneously, thus generalizing the Matom\"aki-Radziwill theorem on multiplicative functions in short intervals. We also consider the maximal deviation of $f$ over all residue classes $a\pmod q$ for $q\leq x^{1/2-\varepsilon}$, and show that it is small for "smooth-supported" $f$, again apart from a nearly power-saving set of exceptional $q$, thus providing a smaller exceptional set than what follows from Bombieri-Vinogradov-type theorems. As an application of our methods, we consider Linnik-type problems for products of exactly three primes, and in particular prove results relating to a ternary version of a conjecture of Erd\H{o}s on representing every element of the multiplicative group $\mathbb{Z}_p^{\times}$ as the product of two primes less than $p$.


title: Human-Centric Program Synthesis

abstract: Program synthesis techniques offer significant new capabilities in searching for programs that satisfy high-level specifications. While synthesis has been thoroughly explored for input/output pair specifications (programming-by-example), this paper asks: what does program synthesis look like beyond examples? What actual issues in day-to-day development would stand to benefit the most from synthesis? How can a human-centric perspective inform the exploration of alternative specification languages for synthesis? I sketch a human-centric vision for program synthesis where programmers explore and learn languages and APIs aided by a synthesis tool.


title: CapExec: Towards Transparently-Sandboxed Services (Extended Version)

abstract: Network services are among the riskiest programs executed by production systems. Such services execute large quantities of complex code and process data from arbitrary and untrusted network sources, often with high levels of system privilege. It is desirable to confine system services to a least-privileged environment so that the potential damage from a malicious attacker can be limited, but existing mechanisms for sandboxing services require invasive and system-specific code changes and are insufficient to confine broad classes of network services. Rather than sandboxing one service at a time, we propose that the best place to add sandboxing to network services is in the service manager that starts those services. As a first step towards this vision, we propose CapExec, a process supervisor that can execute a single service within a sandbox based on a service declaration file in which, required resources whose limited access to are supported by Caper services, are specified. Using the Capsicum compartmentalization framework and its Casper service framework, CapExec provides robust application sandboxing without requiring any modifications to the application itself. We believe that this is a first step towards ubiquitous sandboxing of network services without the costs of virtualization. Keywords: application security, sandboxing, service manager, Capsicum, compartmentalization


title: On the minimal diameter of closed hyperbolic surfaces

abstract: We prove that the minimal diameter of a hyperbolic compact orientable surface of genus $g$ is asymptotic to $\log g$ as $g \to \infty$. The proof relies on a random construction, which we analyse using lattice point counting theory and the exploration of random trivalent graphs.


title: Unitary vs pseudo-unitary time evolution and statistical effects in the dynamical Sauter-Schwinger process

abstract: Dynamical Sauter-Schwinger mechanism of pair creation by a time-dependent electric field comprising of $N_{\rm rep}$ identical pulses is analyzed within the framework of the spinor and scalar quantum electrodynamics. For linearly polarized pulses, both theories predict that a single eigenmode of the matter wave follows the dynamics of a two-level system. This dynamics, however, is either governed by a Hermitian (for spin 1/2 particles) or pseudo-Hermitian (for spin 0 particles) Hamiltonian. Essentially, both theories lead to a Fraunhofer-type enhancement of the momentum distributions of created pairs. While in the fermionic case the enhancement is never perfect and it deteriorates with increasing the number of pulses in a train $N_{\rm rep}$, in the bosonic case we observe the opposite. More specifically, it is at exceptional points where the spectra of bosonic pairs scale exactly as $N_{\rm rep}^2$, and this scaling is even enhanced with increasing the number of pulses in a train.


title: Interaction networks for the identification of boosted $H\to b\overline{b}$ decays

abstract: We develop an algorithm based on an interaction network to identify high-transverse-momentum Higgs bosons decaying to bottom quark-antiquark pairs and distinguish them from ordinary jets that reflect the configurations of quarks and gluons at short distances. The algorithm's inputs are features of the reconstructed charged particles in a jet and the secondary vertices associated with them. Describing the jet shower as a combination of particle-to-particle and particle-to-vertex interactions, the model is trained to learn a jet representation on which the classification problem is optimized. The algorithm is trained on simulated samples of realistic LHC collisions, released by the CMS Collaboration on the CERN Open Data Portal. The interaction network achieves a drastic improvement in the identification performance with respect to state-of-the-art algorithms.


title: Non-Invasive Fuhrman Grading of Clear Cell Renal Cell Carcinoma Using Computed Tomography Radiomics Features and Machine Learning

abstract: Purpose: To identify optimal classification methods for computed tomography (CT) radiomics-based preoperative prediction of clear cells renal cell carcinoma (ccRCC) grade. Methods and material: Seventy one ccRCC patients were included in the study. Three image preprocessing techniques (Laplacian of Gaussian, wavelet filter, and discretization of the intensity values) were applied on tumor volumes. In total, 2530 radiomics features (tumor shape and size, intensity statistics, and texture) were extracted from each segmented tumor volume. Univariate analysis was performed to assess the association of each feature with the histological condition. In the case of multivariate analysis, the following was implemented: three feature selection including the least absolute shrinkage and selection operator (LASSO), students t-test and minimum Redundancy Maximum Relevance (mRMR) algorithms. These selected features were then used to construct three classification models (SVM, random forest, and logistic regression) to discriminate the high from low-grade ccRCC at nephrectomy. Lastly, multivariate model performance was evaluated on the bootstrapped validation cohort using the area under receiver operating characteristic curve (AUC). Results: Univariate analysis demonstrated that among different image sets, 128 bin discretized images have statistically significant different (q-value < 0.05) texture parameters with a mean of AUC 0.74 (q-value < 0.05). The three ML-based classifier shows proficient discrimination of the high from low-grade ccRCC. The AUC was 0.78 in logistic regression, 0.62 in random forest, and 0.83 in SVM model, respectively. Conclusion: Radiomics features can be a useful and promising non-invasive method for preoperative evaluation of ccRCC Fuhrman grades. Key words: RCC, Radiomics, Machine Learning, Computed Tomography


title: Lawson schemes for highly oscillatory stochastic differential equations and conservation of invariants

abstract: In this paper, we consider a class of stochastic midpoint and trapezoidal Lawson schemes for the numerical discretization of highly oscillatory stochastic differential equations. These Lawson schemes incorporate both the linear drift and diffusion terms in the exponential operator. We prove that the midpoint Lawson schemes preserve quadratic invariants and discuss this property as well for the trapezoidal Lawson scheme. Numerical experiments demonstrate that the integration error for highly oscillatory problems is smaller than that of some standard methods.


title: Communication-Constrained Routing and Traffic Control for Autonomous Vehicles

abstract: Autonomous vehicles (AV) is an advanced technology that can bring convenience, improve the road-network throughput, and reduce traffic accidents. To enable higher levels of automation (LoA), massive amounts of sensory data need to be uploaded to the network for processing, and then, maneuvering decisions must be returned to the AV. Furthermore, passengers might have a higher transmission rate demands for various data-hungry and delay-sensitive applications.


title: Attention Forcing for Sequence-to-sequence Model Training

abstract: Auto-regressive sequence-to-sequence models with attention mechanism have achieved state-of-the-art performance in many tasks such as machine translation and speech synthesis. These models can be difficult to train. The standard approach, teacher forcing, guides a model with reference output history during training. The problem is that the model is unlikely to recover from its mistakes during inference, where the reference output is replaced by generated output. Several approaches deal with this problem, largely by guiding the model with generated output history. To make training stable, these approaches often require a heuristic schedule or an auxiliary classifier. This paper introduces attention forcing, which guides the model with generated output history and reference attention. This approach can train the model to recover from its mistakes, in a stable fashion, without the need for a schedule or a classifier. In addition, it allows the model to generate output sequences aligned with the references, which can be important for cascaded systems like many speech synthesis systems. Experiments on speech synthesis show that attention forcing yields significant performance gain. Experiments on machine translation show that for tasks where various re-orderings of the output are valid, guiding the model with generated output history is challenging, while guiding the model with reference attention is beneficial.


title: Knowledge of Process-Structure-Property Relationships to Engineer Better Heat Treatments for Laser Powder Bed Fusion Additive Manufactured Inconel 718

abstract: Dislocation structures, chemical segregation, {\gamma ^{\prime}, {\gamma ^{\prime \prime}}, {\delta} precipitates and Laves phase were quantified within the microstructures of Inconel 718 (IN718) produced by laser powder bed fusion additive manufacturing (AM) and subjected to standard, direct aging, and modified multi-step heat treatments. Additionally, heat-treated samples still attached to the build plates vs. those removed were also documented for a standard heat treatment. The effects of the different resulting microstructures on room temperature strengths and elongations to failure is revealed. Knowledge derived from these process structure property relationships was used to engineer a super solvus solution anneal at 1020 degC for 15 minutes, followed by aging at 720 degC for 24 hours heat treatment for AM-IN718 that eliminates Laves and {\delta} phases, preserves AM specific dislocation cells that are shown to be stabilized by MC carbide particles, and precipitates dense {\gamma ^{\prime} and {\gamma ^{\prime \prime}} nanoparticle populations. This 'optimized for AM-IN718 heat treatment' results in superior properties relative to wrought/additively manufactured, then industry standard heat treated IN718: relative increases of 7/10 percent in yield strength, 2/7 percent in ultimate strength, and 23/57 percent in elongation to failure are realized, respectively, regardless of as-built vs. machined surface finishes.


title: Exascale Deep Learning to Accelerate Cancer Research

abstract: Deep learning, through the use of neural networks, has demonstrated remarkable ability to automate many routine tasks when presented with sufficient data for training. The neural network architecture (e.g. number of layers, types of layers, connections between layers, etc.) plays a critical role in determining what, if anything, the neural network is able to learn from the training data. The trend for neural network architectures, especially those trained on ImageNet, has been to grow ever deeper and more complex. The result has been ever increasing accuracy on benchmark datasets with the cost of increased computational demands. In this paper we demonstrate that neural network architectures can be automatically generated, tailored for a specific application, with dual objectives: accuracy of prediction and speed of prediction. Using MENNDL--an HPC-enabled software stack for neural architecture search--we generate a neural network with comparable accuracy to state-of-the-art networks on a cancer pathology dataset that is also $16\times$ faster at inference. The speedup in inference is necessary because of the volume and velocity of cancer pathology data; specifically, the previous state-of-the-art networks are too slow for individual researchers without access to HPC systems to keep pace with the rate of data generation. Our new model enables researchers with modest computational resources to analyze newly generated data faster than it is collected.


title: Polylogarithmic width suffices for gradient descent to achieve arbitrarily small test error with shallow ReLU networks

abstract: Recent theoretical work has guaranteed that overparameterized networks trained by gradient descent achieve arbitrarily low training error, and sometimes even low test error. The required width, however, is always polynomial in at least one of the sample size $n$, the (inverse) target error $1/\epsilon$, and the (inverse) failure probability $1/\delta$. This work shows that $\widetilde{\Theta}(1/\epsilon)$ iterations of gradient descent with $\widetilde{\Omega}(1/\epsilon^2)$ training examples on two-layer ReLU networks of any width exceeding $\mathrm{polylog}(n,1/\epsilon,1/\delta)$ suffice to achieve a test misclassification error of $\epsilon$. We also prove that stochastic gradient descent can achieve $\epsilon$ test error with polylogarithmic width and $\widetilde{\Theta}(1/\epsilon)$ samples. The analysis relies upon the separation margin of the limiting kernel, which is guaranteed positive, can distinguish between true labels and random labels, and can give a tight sample-complexity analysis in the infinite-width setting


title: Dirac energy spectrum and inverted band gap in metamorphic InAsSb/InSb superlattices

abstract: A Dirac-type energy spectrum was demonstrated in gapless ultra-short-period metamorphic InAsSb/InSb superlattices by angle-resolved photoemission spectroscopy (ARPES_ measurements. The Fermi velocity value 7.4x10^5 m/s in a gapless superlattice with a period of 6.2nm is in a good agreement with the results of magneto-absorption experiments. An "inverted" bandgap opens in the center of the Brillouin zone at higher temperatures and in the SL with a larger period. The ARPES data indicate the presence of a surface electron accumulation layer


title: Non-equilibrium photons from the bottom-up thermalization scenario

abstract: In this work, I calculate the $p_\perp$ resolved spectra for the three stages of the \textit{bottom-up} scenario, which are comparable to the thermal contribution, particularly at higher values of the saturation scale $Q_S^2$. Analytical solutions are obtained by including a parametrization of scaling solutions from far-from-equilibrium classical statistical lattice simulations into a small angle kinetic rate. Furthermore, a theoretically motivated ansatz is used to account for near-collinear enhancement of the low-$p_\perp$ radiation. The system is phenomenologically constrained using the charge hadron multiplicities from LHC and RHIC as in previous parametric estimates and fair agreement with the data available for photons was found. I find that for this realistic set of parameters, the contribution from the thermalizing glasma dominates the excess photons.


title: Quantum Microwave Radiometry with a Superconducting Qubit

abstract: The interaction of photons and coherent quantum systems can be employed to detect electromagnetic radiation with remarkable sensitivity. We introduce a quantum radiometer based on the photon-induced-dephasing process of a superconducting qubit for sensing microwave radiation at the sub-unit-photon level. Using this radiometer, we demonstrated the radiative cooling of a 1-K microwave resonator and measured its mode temperature with an uncertainty ~0.01 K. We have thus developed a precise tool for studying the thermodynamics of quantum microwave circuits, which provides new solutions for calibrating hybrid quantum systems and detecting candidate particles for dark matter.


title: Eigenvalues and dynamical degrees of self-maps on abelian varieties

abstract: Let $X$ be a smooth projective variety over an algebraically closed field, and $f\colon X\to X$ a surjective self-morphism of $X$. The $i$-th cohomological dynamical degree $\chi_i(f)$ is defined as the spectral radius of the pullback $f^{*}$ on the \'etale cohomology group $H^i_{\textrm{\'et}}(X, \mathbf{Q}_\ell)$ and the $k$-th numerical dynamical degree $\lambda_k(f)$ as the spectral radius of the pullback $f^{*}$ on the vector space $\mathsf{N}^k(X)_{\mathbf{R}}$ of real algebraic cycles of codimension $k$ on $X$ modulo numerical equivalence. Truong conjectured that $\chi_{2k}(f) = \lambda_k(f)$ for all $0 \le k \le \dim X$ as a generalization of Weil's Riemann hypothesis. We prove this conjecture in the case of abelian varieties. In the course of the proof we also obtain a new parity result on the eigenvalues of self-maps of abelian varieties in prime characteristic, which is of independent interest.


title: Energy-Based Models for Deep Probabilistic Regression

abstract: While deep learning-based classification is generally tackled using standardized approaches, a wide variety of techniques are employed for regression. In computer vision, one particularly popular such technique is that of confidence-based regression, which entails predicting a confidence value for each input-target pair (x,y). While this approach has demonstrated impressive results, it requires important task-dependent design choices, and the predicted confidences lack a natural probabilistic meaning. We address these issues by proposing a general and conceptually simple regression method with a clear probabilistic interpretation. In our proposed approach, we create an energy-based model of the conditional target density p(y|x), using a deep neural network to predict the un-normalized density from (x,y). This model of p(y|x) is trained by directly minimizing the associated negative log-likelihood, approximated using Monte Carlo sampling. We perform comprehensive experiments on four computer vision regression tasks. Our approach outperforms direct regression, as well as other probabilistic and confidence-based methods. Notably, our model achieves a 2.2% AP improvement over Faster-RCNN for object detection on the COCO dataset, and sets a new state-of-the-art on visual tracking when applied for bounding box estimation. In contrast to confidence-based methods, our approach is also shown to be directly applicable to more general tasks such as age and head-pose estimation. Code is available at https://github.com/fregu856/ebms_regression.


title: Multilepton signatures for scalar dark matter searches in coannihilation scenario

abstract: We revisit the scalar singlet dark matter (DM) scenario with a pair of dark lepton partners which form a vectorlike Dirac fermionic doublet. The extra doublet couples with the Standard Model (SM) leptonic doublet and the scalar singlet via a non-SM-like Yukawa structure. As a result, (i) since the extra fermionic states interact with other dark sector particles as well as the SM via gauge and Yukawa interactions, it gives rise to new DM annihilation processes including pair annihilation as well as coannihilation channels, and (ii) such a Yukawa structure opens up new production channels for leptonic final states giving much enhancement in cross sections to search for dark matter in the LHC. Using suitable kinematic observables, we train a {\em boosted decision tree} (BDT) classifier to separate enhanced but still feeble light leptonic signals from the background in an effective manner. On the other hand, the same technique is applied to study $\tau$-tagged jets in the search for DM signals.


title: Expert2Coder: Capturing Divergent Brain Regions Using Mixture of Regression Experts

abstract: fMRI semantic category understanding using linguistic encoding models attempts to learn a forward mapping that relates stimuli to the corresponding brain activation. State-of-the-art encoding models use a single global model (linear or non-linear) to predict brain activation given the stimulus. However, the critical assumption in these methods is that a priori different brain regions respond the same way to all the stimuli, that is, there is no modularity or specialization assumed for any region. This goes against the modularity theory, supported by many cognitive neuroscience investigations suggesting that there are functionally specialized regions in the brain. In this paper, we achieve this by clustering similar regions together and for every cluster we learn a different linear regression model using a mixture of linear experts model. The key idea here is that each linear expert captures the behaviour of similar brain regions. Given a new stimulus, the utility of the proposed model is twofold (i) predicts the brain activation as a weighted linear combination of the activations of multiple linear experts and (ii) to learn multiple experts corresponding to different brain regions. We argue that each expert captures activity patterns related to a particular region of interest (ROI) in the human brain. This study helps in understanding the brain regions that are activated together given different kinds of stimuli. Importantly, we suggest that the mixture of regression experts (MoRE) framework successfully combines the two principles of organization of function in the brain, namely that of specialization and integration. Experiments on fMRI data from paradigm 1 [1]where participants view linguistic stimuli show that the proposed MoRE model has better prediction accuracy compared to that of conventional models.


title: Affleck-Dine inflation

abstract: The Affleck-Dine mechanism in its simplest form provides baryogenesis from the out-of-equilibrium evolution of a complex scalar field with a simple renormalizable potential. We show that such a model, supplemented by nonminimal coupling to gravity, can also provide inflation, consistent with Planck constraints, simultaneously with the generation of the baryon asymmetry. The predictions of the model include significant tensor-to-scalar ratio and possibly baryon isocurvature fluctuations. The reheating temperature is calculable, making the model fully predictive. We require color triplet scalars for reheating and transfering the primordial baryon asymmetry to quarks; these could be observable at colliders. They can also be probed at higher scales by searches for quark compositeness in dijet angular distributions, and flavor-changing neutral current effects.


title: DBRec: Dual-Bridging Recommendation via Discovering Latent Groups

abstract: In recommender systems, the user-item interaction data is usually sparse and not sufficient for learning comprehensive user/item representations for recommendation. To address this problem, we propose a novel dual-bridging recommendation model (DBRec). DBRec performs latent user/item group discovery simultaneously with collaborative filtering, and interacts group information with users/items for bridging similar users/items. Therefore, a user's preference over an unobserved item, in DBRec, can be bridged by the users within the same group who have rated the item, or the user-rated items that share the same group with the unobserved item. In addition, we propose to jointly learn user-user group (item-item group) hierarchies, so that we can effectively discover latent groups and learn compact user/item representations. We jointly integrate collaborative filtering, latent group discovering and hierarchical modelling into a unified framework, so that all the model parameters can be learned toward the optimization of the objective function. We validate the effectiveness of the proposed model with two real datasets, and demonstrate its advantage over the state-of-the-art recommendation models with extensive experiments.


title: New Attacks and Defenses for Randomized Caches

abstract: The last level cache is vulnerable to timing based side channel attacks because it is shared by the attacker and the victim processes even if they are located on different cores. These timing attacks evict the victim cache lines using small conflict groups(SCG), and monitor the cache to observe when the victim uses these cache lines again. A conflict group is a collection of cache lines which will evict the target cache line. Randomization is often used by defenses to prevent creation of SCGs. We introduce new attacks to demonstrate that the current randomization schemes require an extremely high refresh rate to be secure, on average a 15\% performance overhead, and upto 50\% in the worst case. Next, we propose a new randomization strategy using an indirection table, which mitigates this issue. Addresses of cache lines are encrypted and used to lookup the indirection table entry. Each indirection table entry stores a mapping to a randomly chosen cache set. The cache line is placed into this randomly chosen set. The encryption key changes upto 50x faster than CEASER's default rate, by using evictions to trigger the re-randomization. Instead of moving cache lines, this mechanism re-randomizes one iTable entry at a time, whenever the cache lines corresponding to the iTable entry are naturally evicted. Thus, the miss rate is not much worse than the baseline. We quantitatively show that our scheme does almost as well as a fully associative cache to defend against these attacks. We also demonstrate new attacks that target the iTable by oversubscribing its entries, and quantitatively show that our scheme is resilient against new attacks for trillions of years. We estimate low area ( < 7\%) and power overhead compared to a baseline inclusive last-level cache. Lastly, we evaluate a low performance overhead (<4%) using the SPECrate 2017 and PARSEC 3.0 benchmarks.


title: Thermal Regulation and the Star-Forming Main Sequence

abstract: We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming "main sequence". To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.


title: Bifurcating entanglement-renormalization group flows of fracton stabilizer models

abstract: We investigate the entanglement-renormalization group flows of translation-invariant topological stabilizer models in three dimensions. Fracton models are observed to bifurcate under entanglement renormalization, generically returning at least one copy of the original model. Based on this behavior we formulate the notion of bifurcated equivalence for fracton phases, generalizing foliated fracton equivalence. The notion of quotient superselection sectors is also generalized accordingly. We calculate bifurcating entanglement-renormalization group flows for a wide range of examples and, based on those results, propose conjectures regarding the classification of translation-invariant topological stabilizer models in three dimensions.


title: Chemical Abundances in a Turbulent Medium -- H$_2$, OH$^+$, H$_2$O$^+$, ArH$^+$

abstract: Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters.


title: Fractons from confinement in one dimension

abstract: Recent work has shown that two seemingly different physical mechanisms, namely fracton behavior and confinement, can give rise to non-ergodicity in one-dimensional quantum many-body systems. In this work, we demonstrate an intrinsic link between these two mechanisms by studying the dynamics of one-dimensional confining theories, such as a U(1) gauge theory and a quantum Ising model. We show that, within certain parameter regimes, these models exhibit effective fracton dynamics, characterized by immobility of stable single-particle excitations and free motion of dipolar bound states. By perturbatively integrating out the linearly confining field, we obtain an effective fracton Hamiltonian for the confined charges which exhibits conservation of dipole moment. We discuss an intuitive understanding of these results in terms of the motion of the confining strings, leading to potential extensions to higher dimensions. We thereby interpret recent observations of nonthermal eigenstates and glassy dynamics in confining theories in terms of corresponding results in the fracton literature.


title: Modeling Turbulence in Galactic Centers

abstract: Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way's Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a Smoothed Particle Hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.


title: Josephson junctions of Weyl and multi-Weyl semimetals

abstract: We study a Josephson junction involving a Weyl and a multi-Weyl semimetal separated by a barrier region of width $d$ created by putting a gate voltage $U_0$ over the Weyl semimetal. The topological winding number of such a junction changes across the barrier. We show that $I_c R_N$ for such junctions, where $I_c$ is the critical current and $R_N$ the normal state resistance, in the thin barrier limit, has a universal value independent of the barrier potential. We provide an analytical expression of the Andreev bound states and use it to demonstrate that the universal value of $I_c R_N$ is a consequence of change in topological winding number across the junction. We also study AC Josephson effect in such a junction in the presence of an external microwave radiation, chart out its current-voltage characteristics, and show that the change in the winding number across the junction shapes the properties of its Shapiro steps. We discuss the effect of increasing barrier thickness $d$ on the above-mentioned properties and chart out experiments which may test our theory.


title: Quasars as standard candles II: The non linear relation between UV and X-ray emission at high redshifts

abstract: A tight non-linear relation between the X-ray and the optical-ultraviolet (UV) emission has been observed in Active Galactic Nuclei (AGN) over a wide range of redshift and several orders of magnitude in luminosity, suggesting the existence of an ubiquitous physical mechanism regulating the energy transfer between the accretion disc and the X-ray emitting corona. Recently, our group developed a method to use this relation in the observational cosmology, turning quasars into standardizable candles. This work has the main aim to investigate the potential evolution of this correction at high redshifts. We thus studied the $L_{\rm X}-L_{\rm UV}$ relation for a sample of quasars in the redshift range 4<$z$<7, adopting the selection criteria proposed in our previous work regarding their spectral properties. The resulting sample consists of 53 Type 1 (unobscured) quasars, observed either with Chandra or XMM-Newton, for which we performed a full spectral analysis, determining the rest-frame 2 keV flux density, as well as more general X-ray properties such as the estimate of photon index, and the soft (0.5-2 keV) and hard (2-10 keV) unabsorbed luminosities. We find that the relation shows no evidence for evolution with redshift. The intrinsic dispersion of the L$_X$-L$_{UV}$ for a sample free of systematics/contaminants is of the order of 0.22 dex, which is consistent with previous estimates from our group on quasars at lower redshift.


title: Spectral fingerprints of the non-linear dynamics of driven superconductors with dissipation

abstract: We study the out-of-equilibrium dynamics of a BCS superconductor in the presence of a periodic drive with a frequency $\omega_d$ larger than the equilibrium superconducting gap and an external bath providing dissipation. Similar to the dissipationless case, a subset of quasiparticles, resonant with the drive, synchronizes and produces interesting non-linear phenomena. For small dissipation Rabi-Higgs oscillations [Collado et al. Phys. Rev. B 98, 214519 (2018)] can be observed as a transient effect. At long-times in contrast, dissipation leads to a steady state with a population imbalance that increases quadratically for small drives and saturates for large drives as in non-linear quantum optics. We show also that second harmonic generation is allowed for a drive which acts on the BCS coupling constant. We also compute the intensity of time- and angle-resolved photo-emission spectroscopy (tr-ARPES) and time-resolved tunneling spectra. The excited quasiparticle population appears as a decrease (increase) in the photo-emission intensity at energy $-\omega_d/2$ ($+\omega_d/2$) measured from the chemical potential. The tunneling intensity shows a time-dependent depression at $\pm\omega_d/2$ due to the lacking of spectral density and population unbalance causing Pauli blockade. We predict that at short times, compared to the energy relaxation time, both experiments will show oscillations with the Rabi-Higgs frequency.


title: Non-equilibrium criticality in quench dynamics of infinite-range spin models

abstract: Long-range interacting spin systems are ubiquitous in physics and exhibit a variety of ground state disorder-to-order phase transitions. We consider a prototype of infinite-range interacting models known as the Lipkin-Meshkov-Glick (LMG) model describing the collective interaction of $N$ spins, and investigate the dynamical properties of fluctuations and correlations after a sudden quench of the Hamiltonian. Specifically, we focus on critical quenches, where the initial state and/or the quench Hamiltonian are critical. Depending on the type of quench, we identify three distinct behaviors where both the short-time dynamics and the stationary state at long times are effectively thermal, quantum, and genuinely non-equilibrium, characterized by distinct universality classes and static and dynamical critical exponents. These behaviors can be identified by an infrared effective temperature that is finite, zero, and infinite (the latter scaling with the system size as $N^{1/3}$), respectively. The quench dynamics is studied through a combination of exact numerics and analytical calculations utilizing the non-equilibrium Keldysh field theory. Our results are amenable to realization in experiments with trapped-ion experiments where long-range interactions naturally arise.

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


title: Extrasolar Kuiper Belts

abstract: Extrasolar debris disks are the dust disks found around nearby main sequence stars arising from the break-up of asteroids and comets orbiting the stars. Far-IR surveys (e.g., with Herschel) showed that ~20% of stars host detectable dust levels. While dust temperatures suggest a location at 10s of au comparable with our Kuiper belt, orders of magnitude more dust is required implying a planetesimal population more comparable with the primordial Kuiper belt. High resolution imaging (e.g., with ALMA) has mapped the nearest and brightest disks, providing evidence for structures shaped by an underlying planetary system. Some of these are analogous to structures in our own Kuiper belt (e.g., the hot and cold classical, resonant, scattered disk and cometary populations), while others have no Solar System counterpart. CO gas is seen in some debris disks, and inferred to originate in the destruction of planetesimals with a similar volatile-rich composition to Solar System comets. This chapter reviews our understanding of extrasolar Kuiper belts and of how our own Kuiper belt compares with those of neighbouring stars.


title: A Conceptual Introduction to Markov Chain Monte Carlo Methods

abstract: Markov Chain Monte Carlo (MCMC) methods have become a cornerstone of many modern scientific analyses by providing a straightforward approach to numerically estimate uncertainties in the parameters of a model using a sequence of random samples. This article provides a basic introduction to MCMC methods by establishing a strong conceptual understanding of what problems MCMC methods are trying to solve, why we want to use them, and how they work in theory and in practice. To develop these concepts, I outline the foundations of Bayesian inference, discuss how posterior distributions are used in practice, explore basic approaches to estimate posterior-based quantities, and derive their link to Monte Carlo sampling and MCMC. Using a simple toy problem, I then demonstrate how these concepts can be used to understand the benefits and drawbacks of various MCMC approaches. Exercises designed to highlight various concepts are also included throughout the article.


title: Laparoscopy Surgery CO2 Removal via Generative Adversary Network and Dark Channel Prior

abstract: Laparoscopic surgery uses a thin tube with a camera called a laparoscope, which is inserted into the abdomen through a small incision in the skin during surgery. This allows to a surgeon to see inside of the body without causing significant injury to the patient. These characteristics make laparoscopy a widely used technique. In laparoscopic surgery, image quality can be severely degraded by surgical smoke caused by the use of tissue dissection tools which reduces the visibility of the observed organs and tissues. This lack of visibility increases the possibility of errors and surgery time with the consequences that this may have on the patient's health. In this paper, we introduce a novel hybrid approach for computational smoke removal which is based on the combination of a widely dehazing method used: the dark channel prior (DCP) and a pixel-to-pixel neural network approach: Generative Adversary Network (GAN). The experimental results have proven that the proposed method achieves a better performance than the individual results of the DCP and GAN in terms of restoration quality, obtaining a PSNR value of 25 and SSIM index of 0.88 over a test set of synthetic images.


title: Pseudo-rotations and Steenrod squares revisited

abstract: In this note we prove that if a closed monotone symplectic manifold admits a Hamiltonian pseudo-rotation, which may be degenerate, then the quantum Steenrod square of the cohomology class Poincar\'{e} dual to the point must be deformed. This result gives restrictions on the existence of pseudo-rotations, implying a form of uniruledness by pseudo-holomorphic spheres, and generalizes a recent result of the author. The new component in the proof consists in an elementary calculation with capped periodic orbits.


title: Preference-Based Learning for Exoskeleton Gait Optimization

abstract: This paper presents a personalized gait optimization framework for lower-body exoskeletons. Rather than optimizing numerical objectives such as the mechanical cost of transport, our approach directly learns from user preferences, e.g., for comfort. Building upon work in preference-based interactive learning, we present the CoSpar algorithm. CoSpar prompts the user to give pairwise preferences between trials and suggest improvements; as exoskeleton walking is a non-intuitive behavior, users can provide preferences more easily and reliably than numerical feedback. We show that CoSpar performs competitively in simulation and demonstrate a prototype implementation of CoSpar on a lower-body exoskeleton to optimize human walking trajectory features. In the experiments, CoSpar consistently found user-preferred parameters of the exoskeleton's walking gait, which suggests that it is a promising starting point for adapting and personalizing exoskeletons (or other assistive devices) to individual users.


title: Bayesian approach to constraining the properties of ionized bubbles during reionization

abstract: A possible way to study the reionization of cosmic hydrogen is by observing the large ionized regions (bubbles) around bright individual sources, e.g., quasars, using the redshifted 21 cm signal. It has already been shown that matched filter-based methods are not only able to detect the weak 21 cm signal from these bubbles but also aid in constraining their properties. In this work, we extend the previous studies to develop a rigorous Bayesian framework to explore the possibility of constraining the parameters that characterize the bubbles. To check the accuracy with which we can recover the bubble parameters, we apply our method on mock observations appropriate for the upcoming SKA1-low. For a region of size $\gtrsim 50$ cMpc around a typical quasar at redshift 7, we find that $\approx 20$ h of integration with SKA1-low will be able to constrain the size and location of the bubbles, as well as the difference in the neutral hydrogen fraction inside and outside the bubble, with $\lesssim 10\%$ precision. The recovery of the parameters are more precise and the SNR of the detected signal is higher when the bubble sizes are larger and their shapes are close to spherical. Our method can be useful in identifying regions in the observed field which contain large ionized regions and hence are interesting for following up with deeper integration times.


title: J0453+1559: a neutron star-white dwarf binary from a thermonuclear electron-capture supernova?

abstract: The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictive concerning the possibility to produce such a low-mass neutron star remnant. Here we challenge the neutron star interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen-neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several $0.1\,M_\odot$ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR~J0453+1559 by a neutron star--white dwarf system. More work on tECSNe is needed to assess the viability of this scenario.


title: Accurate many-body electronic structure near the basis set limit: application to the chromium dimer

abstract: We describe a method for computing near-exact energies for correlated systems with large Hilbert spaces. The method efficiently identifies the most important basis states (Slater determinants) and performs a variational calculation in the subspace spanned by these determinants. A semistochastic approach is then used to add a perturbative correction to the variational energy to compute the total energy. The size of the variational space is progressively increased until the total energy converges to within the desired tolerance. We demonstrate the power of the method by computing a near-exact potential energy curve (PEC) for a very challenging molecule -- the chromium dimer.


title: Diverse outcomes of planet formation and composition around low-mass stars and brown dwarfs

abstract: The detection of Earth-size exoplanets around low-mass stars -- in stars such as Proxima Centauri and TRAPPIST-1 -- provide an exceptional chance to improve our understanding of the formation of planets around M stars and brown dwarfs. We explore the formation of such planets with a population synthesis code based on a planetesimal-driven model previously used to study the formation of the Jovian satellites. Because the discs have low mass and the stars are cool, the formation is an inefficient process that happens at short periods, generating compact planetary systems. Planets can be trapped in resonances and we follow the evolution of the planets after the gas has dissipated and they undergo orbit crossings and possible mergers. We find that formation of planets above Mars mass and in the planetesimal accretion scenario, is only possible around stars with masses $M_{\star} \ge 0.07 M_{sun}$ and discs of $M_{disc} \ge 10^{-2}~M_{sun} $. We find that planets above Earth-mass form around stars with masses larger than $0.15 M_{\oplus}$ while planets larger than $5 M_{\oplus}$do not form in our model, even not under the most optimal conditions (massive disc), showing that planets such as GJ 3512b form with another, more efficient mechanism. Our results show that the majority of planets form with a significant water fraction; that most of our synthetic planetary systems have 1, 2, or 3 planets, but those with 4, 5, 6, and 7 planets are also common, confirming that compact planetary systems with many planets should be a relatively common outcome of planet formation around small stars.


title: Variational point-obstacle avoidance on Riemannian manifolds

abstract: In this letter we study variational obstacle avoidance problems on complete Riemannian manifolds. The problem consists of minimizing an energy functional depending on the velocity, covariant acceleration and a repulsive potential function used to avoid a static obstacle on the manifold, among a set of admissible curves. We derive the dynamical equations for extrema of the variational problem, in particular on compact connected Lie groups and Riemannian symmetric spaces. Numerical examples are presented to illustrate the proposed method.


title: How long is the convex minorant of a one-dimensional random walk?

abstract: We prove distributional limit theorems for the length of the largest convex minorant of a one-dimensional random walk with independent identically distributed increments. Depending on the increment law, there are several regimes with different limit distributions for this length. Among other tools, a representation of the convex minorant of a random walk in terms of uniform random permutations is utilized.


title: Cataclysmic variable evolution and the white dwarf mass problem: A Review

abstract: Although the theory of cataclysmic variable (CV) evolution is able to explain several observational aspects, strong discrepancies have existed for decades between observations and theoretical predictions of the orbital period distribution, the location of the minimum period, and the space density of CVs. Moreover, it has been shown in the last decade that the average white dwarf (WD) mass observed in CVs is significantly higher than the average mass in single WDs or in detached progenitors of CVs, and that there is an absence of helium-core WDs in CVs which is not observed in their immediate detached progenitors. This highly motivated us to revise the theory of CV formation and evolution. A new empirical model for angular momentum loss in CVs was developed in order to explain the high average WD mass observed and the absence of systems with helium-core WDs. This model seems to help, at the same time, with all of the above mentioned disagreements between theory and observations. Moreover, it also provides us with a very likely explanation for the existence of low-mass WDs without a companion. Here we will review the standard model for CV evolution and the disagreements that have existed for decades between simulations and observations with their possible solutions and/or improvements. We will also summarize the recently confirmed disagreement related to the average WD mass and the fraction of helium-core WDs among CVs, as well as the development of an empirical model that allows us to solve all the disagreements, discussing the physics that could be involved.


title: Learning Generalizable Locomotion Skills with Hierarchical Reinforcement Learning

abstract: Learning to locomote to arbitrary goals on hardware remains a challenging problem for reinforcement learning. In this paper, we present a hierarchical learning framework that improves sample-efficiency and generalizability of locomotion skills on real-world robots. Our approach divides the problem of goal-oriented locomotion into two sub-problems: learning diverse primitives skills, and using model-based planning to sequence these skills. We parametrize our primitives as cyclic movements, improving sample-efficiency of learning on a 18 degrees of freedom robot. Then, we learn coarse dynamics models over primitive cycles and use them in a model predictive control framework. This allows us to learn to walk to arbitrary goals up to 12m away, after about two hours of training from scratch on hardware. Our results on a Daisy hexapod hardware and simulation demonstrate the efficacy of our approach at reaching distant targets, in different environments and with sensory noise.


title: Crowdsourcing via Pairwise Co-occurrences: Identifiability and Algorithms

abstract: The data deluge comes with high demands for data labeling. Crowdsourcing (or, more generally, ensemble learning) techniques aim to produce accurate labels via integrating noisy, non-expert labeling from annotators. The classic Dawid-Skene estimator and its accompanying expectation maximization (EM) algorithm have been widely used, but the theoretical properties are not fully understood. Tensor methods were proposed to guarantee identification of the Dawid-Skene model, but the sample complexity is a hurdle for applying such approaches---since the tensor methods hinge on the availability of third-order statistics that are hard to reliably estimate given limited data. In this paper, we propose a framework using pairwise co-occurrences of the annotator responses, which naturally admits lower sample complexity. We show that the approach can identify the Dawid-Skene model under realistic conditions. We propose an algebraic algorithm reminiscent of convex geometry-based structured matrix factorization to solve the model identification problem efficiently, and an identifiability-enhanced algorithm for handling more challenging and critical scenarios. Experiments show that the proposed algorithms outperform the state-of-art algorithms under a variety of scenarios.


title: Model Pruning Enables Efficient Federated Learning on Edge Devices

abstract: Federated learning (FL) allows model training from local data collected by edge/mobile devices while preserving data privacy, which has wide applicability to image and vision applications. A challenge is that client devices in FL usually have much more limited computation and communication resources compared to servers in a datacenter. To overcome this challenge, we propose PruneFL -- a novel FL approach with adaptive and distributed parameter pruning, which adapts the model size during FL to reduce both communication and computation overhead and minimize the overall training time, while maintaining a similar accuracy as the original model. PruneFL includes initial pruning at a selected client and further pruning as part of the FL process. The model size is adapted during this process, which includes maximizing the approximate empirical risk reduction divided by the time of one FL round. Our experiments with various datasets on edge devices (e.g., Raspberry Pi) show that: (i) we significantly reduce the training time compared to conventional FL and various other pruning-based methods; (ii) the pruned model with automatically determined size converges to an accuracy that is very similar to the original model, and it is also a lottery ticket of the original model.


title: Electrostatic dust-acoustic rogue waves in an electron depleted dusty plasma

abstract: The formation of the gigantic dust-acoustic rouge waves (DARWs) in an electron depleted unmagnetized opposite polarity dusty plasma system is theoretically observed for the first time. The nonlinear Schr\"{o}dinger equation (derived by utilizing the reductive perturbation method) has been analytically as well as numerically analyzed to identify the basic features (viz., height, thickness, and modulational instability, etc.) of DARWs. The results obtained form this investigation should be useful in understanding the basic properties of these rouge waves which can predict to be formed in electron depleted unmagnetized opposite polarity dusty plasma systems like mesosphere, F-rings of Saturn, and cometary atmosphere, etc.


title: Approximation Algorithms for Process Systems Engineering

abstract: Designing and analyzing algorithms with provable performance guarantees enables efficient optimization problem solving in different application domains, e.g.\ communication networks, transportation, economics, and manufacturing. Despite the significant contributions of approximation algorithms in engineering, only limited and isolated works contribute from this perspective in process systems engineering. The current paper discusses three representative, NP-hard problems in process systems engineering: (i) pooling, (ii) process scheduling, and (iii) heat exchanger network synthesis. We survey relevant results and raise major open questions. Further, we present approximation algorithms applications which are relevant to process systems engineering: (i) better mathematical modeling, (ii) problem classification, (iii) designing solution methods, and (iv) dealing with uncertainty. This paper aims to motivate further research at the intersection of approximation algorithms and process systems engineering.


title: Scaling Local Control to Large-Scale Topological Navigation

abstract: Visual topological navigation has been revitalized recently thanks to the advancement of deep learning that substantially improves robot perception. However, the scalability and reliability issue remain challenging due to the complexity and ambiguity of real world images and mechanical constraints of real robots. We present an intuitive solution to show that by accurately measuring the capability of a local controller, large-scale visual topological navigation can be achieved while being scalable and robust. Our approach achieves state-of-the-art results in trajectory following and planning in large-scale environments. It also generalizes well to real robots and new environments without retraining or finetuning.


title: An Optimal Control Approach to Flocking

abstract: Flocking behavior has attracted considerable attention in multi-agent systems. The structure of flocking has been predominantly studied through the application of artificial potential fields coupled with velocity consensus. These approaches, however, do not consider the energy cost of the agents during flocking, which is especially important in large-scale robot swarms. This paper introduces an optimal control framework to induce flocking in a group of agents. Guarantees of energy minimization and safety are provided, along with a decentralized algorithm that satisfies the optimality conditions and can be realized in real time. The efficacy of the proposed control algorithm is evaluated through simulation in both MATLAB and Gazebo.


title: SIMEX Estimation in Parametric Modal Regression with Measurement Error

abstract: For a class of parametric modal regression models with measurement error, a simulation extrapolation estimation procedure is proposed in this paper for estimating the modal regression coefficients. Large sample properties of the proposed estimation procedure, including the consistency and asymptotic normality, are thoroughly investigated. Simulation studies are conducted to evaluate its robustness to potential outliers and the effectiveness in reducing the bias caused by the measurement error.


title: Investigating Particle Acceleration in the Wolf-Rayet Bubble G2.4+1.4

abstract: The supersonic winds produced by massive stars carry a large amount of kinetic power. In numerous scenarios such winds have been proven to produce shocks in which relativistic particles are accelerated emitting non-thermal radiation. Here, we report the first detection of non-thermal emission from a single stellar bubble, G2.4+1.4, associated with a WO star. We observed this source with the uGMRT in Band 4 ($550-850$ MHz) and Band 5 ($1050-1450$ MHz). We present intensity and spectral index maps for this source that are consistent with synchrotron emission (average spectral index, $\alpha = -0.83 \pm 0.10$). The fraction of the available kinetic wind power that is converted into cosmic ray acceleration is estimated to be of the order of a few per cent. This finding constitutes an observational breakthrough and gives new insight on the non-thermal physical processes taking place in the environments of isolated massive stars. In particular, our results show that non-runaway isolated massive stars are capable of accelerating relativistic particles and are therefore confirmed as sources of Galactic cosmic rays.


title: Cavity-enhanced Raman scattering for in situ alignment and characterization of solid-state microcavities

abstract: We report cavity-enhanced Raman scattering from a single-crystal diamond membrane embedded in a highly miniaturized fully-tunable Fabry-P\'{e}rot cavity. The Raman intensity is enhanced 58.8-fold compared to the corresponding confocal measurement. The strong signal amplification results from the Purcell effect. We show that the cavity-enhanced Raman scattering can be harnessed as a narrowband, high-intensity, internal light-source. The Raman process can be triggered in a simple way by using an optical excitation frequency outside the cavity stopband and is independent of the lateral positioning of the cavity mode with respect to the diamond membrane. The strong Raman signal emerging from the cavity output facilitates in situ mode-matching of the cavity mode to single-mode collection optics; it also represents a simple way of measuring the dispersion and spatial intensity-profile of the cavity modes. The optimization of the cavity performance via the strong Raman process is extremely helpful in achieving efficient cavity-outcoupling of the relatively weak emission of single color-centers such as nitrogen-vacancy centers in diamond or rare-earth ions in crystalline hosts with low emitter density.


title: Spectral decomposition of discrepancy kernels on the Euclidean ball, the special orthogonal group, and the Grassmannian manifold

abstract: To numerically approximate Borel probability measures by finite atomic measures, we study the spectral decomposition of discrepancy kernels when restricted to compact subsets of $\mathbb{R}^d$. For restrictions to the Euclidean ball in odd dimensions, to the rotation group $SO(3)$, and to the Grassmannian manifold $\mathcal{G}_{2,4}$, we compute the kernels' Fourier coefficients and determine their asymptotics. The $L_2$-discrepancy is then expressed in the Fourier domain that enables efficient numerical minimization based on the nonequispaced fast Fourier transform. For $SO(3)$, the nonequispaced fast Fourier transform is publicly available, and, for $\mathcal{G}_{2,4}$, the transform is derived here. We also provide numerical experiments for $SO(3)$ and $\mathcal{G}_{2,4}$.


title: Rethinking Text Attribute Transfer: A Lexical Analysis

abstract: Text attribute transfer is modifying certain linguistic attributes (e.g. sentiment, style, authorship, etc.) of a sentence and transforming them from one type to another. In this paper, we aim to analyze and interpret what is changed during the transfer process. We start from the observation that in many existing models and datasets, certain words within a sentence play important roles in determining the sentence attribute class. These words are referred to as \textit{the Pivot Words}. Based on these pivot words, we propose a lexical analysis framework, \textit{the Pivot Analysis}, to quantitatively analyze the effects of these words in text attribute classification and transfer. We apply this framework to existing datasets and models and show that: (1) the pivot words are strong features for the classification of sentence attributes; (2) to change the attribute of a sentence, many datasets only requires to change certain pivot words; (3) consequently, many transfer models only perform the lexical-level modification, while leaving higher-level sentence structures unchanged. Our work provides an in-depth understanding of linguistic attribute transfer and further identifies the future requirements and challenges of this task\footnote{Our code can be found at https://github.com/FranxYao/pivot_analysis}.


title: Pure point spectrum for the Maryland model: a constructive proof

abstract: We develop a constructive method to prove and study pure point spectrum for the Maryland model with Diophantine frequencies.


title: Viscosity solutions for controlled McKean--Vlasov jump-diffusions

abstract: We study a class of non linear integro-differential equations on the Wasserstein space related to the optimal control of McKean--Vlasov jump-diffusions. We develop an intrinsic notion of viscosity solutions that does not rely on the lifting to an Hilbert space and prove a comparison theorem for these solutions. We also show that the value function is the unique viscosity solution.


title: Hardware Design and Analysis of the ACE and WAGE Ciphers

abstract: This paper presents the hardware design and analysis of ACE and WAGE, two candidate ciphers for the NIST Lightweight Cryptography standardization. Both ciphers use sLiSCP's unified sponge duplex mode. ACE has an internal state of 320 bits, uses three 64 bit Simeck boxes, and implements both authenticated encryption and hashing. WAGE is based on the Welch-Gong stream cipher and provides authenticated encryption. WAGE has 259 bits of state, two 7 bit Welch-Gong permutations, and four lightweight 7 bit S-boxes. ACE and WAGE have the same external interface and follow the same I/O protocol to transition between phases. The paper illustrates how a hardware perspective influenced key aspects of the ACE and WAGE algorithms. The paper reports area, power, and energy results for both serial and parallel (unrolled) implementations using four different ASIC libraries: two 65 nm libraries, a 90 nm library, and a 130 nm library. ACE implementations range from a throughput of 0.5 bits-per-clock cycle (bpc) and an area of 4210 GE (averaged across the four ASIC libraries) up to 4 bpc and 7260 GE. WAGE results range from 0.57 bpc with 2920 GE to 4.57 bpc with 11080 GE.


title: Coin_flipper at eHealth-KD Challenge 2019: Voting LSTMs for Key Phrases and Semantic Relation Identification Applied to Spanish eHealth Texts

abstract: This paper describes our approach presented for the eHealth-KD 2019 challenge. Our participation was aimed at testing how far we could go using generic tools for Text-Processing but, at the same time, using common optimization techniques in the field of Data Mining. The architecture proposed for both tasks of the challenge is a standard stacked 2-layer bi-LSTM. The main particularities of our approach are: (a) The use of a surrogate function of F1 as loss function to close the gap between the minimization function and the evaluation metric, and (b) The generation of an ensemble of models for generating predictions by majority vote. Our system ranked second with an F1 score of 62.18% in the main task by a narrow margin with the winner that scored 63.94%.


title: At Stability's Edge: How to Adjust Hyperparameters to Preserve Minima Selection in Asynchronous Training of Neural Networks?

abstract: Background: Recent developments have made it possible to accelerate neural networks training significantly using large batch sizes and data parallelism. Training in an asynchronous fashion, where delay occurs, can make training even more scalable. However, asynchronous training has its pitfalls, mainly a degradation in generalization, even after convergence of the algorithm. This gap remains not well understood, as theoretical analysis so far mainly focused on the convergence rate of asynchronous methods. Contributions: We examine asynchronous training from the perspective of dynamical stability. We find that the degree of delay interacts with the learning rate, to change the set of minima accessible by an asynchronous stochastic gradient descent algorithm. We derive closed-form rules on how the learning rate could be changed, while keeping the accessible set the same. Specifically, for high delay values, we find that the learning rate should be kept inversely proportional to the delay. We then extend this analysis to include momentum. We find momentum should be either turned off, or modified to improve training stability. We provide empirical experiments to validate our theoretical findings.


title: Conserved Restricted Solid-on-Solid Model

abstract: We study a model for the movement of surfaces, namely the conserved, restricted solid-on-solid model. The surface configurations are restricted such that the difference between the heights at adjacent sites is no more than one. In addition the total number of particles is preserved by the dynamics of the model. Mean-field approximations are used to approximate the one-site probability function of the model.


title: Compressed Sensing Microscopy with Scanning Line Probes

abstract: In applications of scanning probe microscopy, images are acquired by raster scanning a point probe across a sample. Viewed from the perspective of compressed sensing (CS), this pointwise sampling scheme is inefficient, especially when the target image is structured. While replacing point measurements with delocalized, incoherent measurements has the potential to yield order-of-magnitude improvements in scan time, implementing the delocalized measurements of CS theory is challenging. In this paper we study a partially delocalized probe construction, in which the point probe is replaced with a continuous line, creating a sensor which essentially acquires line integrals of the target image. We show through simulations, rudimentary theoretical analysis, and experiments, that these line measurements can image sparse samples far more efficiently than traditional point measurements, provided the local features in the sample are enough separated. Despite this promise, practical reconstruction from line measurements poses additional difficulties: the measurements are partially coherent, and real measurements exhibit nonidealities. We show how to overcome these limitations using natural strategies (reweighting to cope with coherence, blind calibration for nonidealities), culminating in an end-to-end demonstration.


title: Line shapes in narrow-line Seyfert 1 galaxies: a tracer of physical properties?

abstract: Line profiles can provide fundamental information on the physics of active galactic nuclei (AGN). In the case of narrow-line Seyfert 1 galaxies (NLS1s) this is of particular importance since past studies revealed how their permitted line profiles are well reproduced by a Lorentzian function instead of a Gaussian. This has been explained with different properties of the broad-line region (BLR), which may present a more pronounced turbulent motions in NLS1s with respect to other AGN. We investigated the line profiles in a recent large NLS1 sample classified using SDSS, and we divided the sources into two subsamples according to their line shapes, Gaussian or Lorentzian. The line profiles clearly separate all the properties of NLS1s. Black hole mass, Eddington ratio, [O III], and Fe II strength are all very different in the Lorentzian and Gaussian samples. We interpret this in terms of evolution within the class of NLS1s. The Lorentzian sources may be the youngest objects, while Gaussian profiles may be typically associated to more evolved objects. Further detailed spectroscopic studies are needed to fully confirm our hypothesis.


title: Thermodynamics and reentrant phase transition for logarithmic nonlinear charged black holes in massive gravity

abstract: We investigate a new class of $(n+1)$-dimensional topological black hole solutions in the context of massive gravity and in the presence of logarithmic nonlinear electrodynamics. Exploring higher dimensional solutions in massive gravity coupled to nonlinear electrodynamics is motivated by holographic hypothesis as well as string theory. We first construct exact solutions of the field equations and then explore the behavior of the metric functions for different values of the model parameters. We observe that our black holes admit the multi-horizons caused by a quantum effect called anti-evaporation. Next, by calculating the conserved and thermodynamic quantities, we obtain a generalized Smarr formula. We find that the first law of black holes thermodynamics is satisfied on the black hole horizon. We study thermal stability of the obtained solutions in both canonical and grand canonical ensembles. We reveal that depending on the model parameters, our solutions exhibit a rich variety of phase structures. Finally, we explore, for the first time without extending thermodynamics phase space, the critical behavior and reentrant phase transition for black hole solutions in massive gravity theory. We realize that there is a zeroth order phase transition for a specified range of charge value and the system experiences a large/small/large reentrant phase transition due to the presence of nonlinear electrodynamics.


title: Optimizing simulation parameters for weak lensing analyses involving non-Gaussian observables

abstract: We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in raytracing DM-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light-cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counter-intuitively, we find that using thin lens planes ($< 60~h^{-1}$Mpc on a $240~h^{-1}$Mpc simulation box) suppresses the power spectrum over a broad range of scales beyond what would be acceptable for an LSST-type survey. A mass resolution of $7.2\times 10^{11}~h^{-1}\,M_{\odot}$ per DM particle (or 256$^3$ particles in a ($240~h^{-1}$Mpc)$^3$ box) is sufficient to extract information using the power spectrum and non-Gaussian statistics from weak lensing data at angular scales down to 1 arcmin with LSST-like levels of shape noise.


title: System-level, Input-output and New Parameterizations of Stabilizing Controllers, and Their Numerical Computation

abstract: It is known that the set of internally stabilizing controller $\mathcal{C}_{\text{stab}}$ is non-convex, but it admits convex characterizations using certain closed-loop maps: a classical result is the Youla parameterization, and two recent notions are the system-level parameterization (SLP) and the input-output parameterization (IOP). In this paper, we address the existence of new convex parameterizations and discuss potential tradeoffs of each parametrization in different scenarios. Our main contributions are: 1) We reveal that only four groups of stable closed-loop transfer matrices are equivalent to internal stability: one of them is used in the SLP, another one is used in the IOP, and the other two are new, leading to two new convex parameterizations of $\mathcal{C}_{\text{stab}}$. 2) We investigate the properties of these parameterizations after imposing the finite impulse response (FIR) approximation, revealing that the IOP has the best ability of approximating $\mathcal{C}_{\text{stab}}$ given FIR constraints. 3) These four parameterizations require no \emph{a priori} doubly-coprime factorization of the plant, but impose a set of equality constraints. However, these equality constraints will never be satisfied exactly in floating-point arithmetic computation and/or implementation. We prove that the IOP is numerically robust for open-loop stable plants, in the sense that small mismatches in the equality constraints do not compromise the closed-loop stability; but a direct IOP implementation will fail to stabilize open-loop unstable systems in practice. The SLP is known to enjoy numerical robustness in the state feedback case; here, we show that numerical robustness of the four-block SLP controller requires case-by-case analysis even the plant is open-loop stable.


title: The structure of the blue whirl revealed

abstract: The blue whirl is a small, stable, spinning blue flame that evolved spontaneously in laboratory experiments while studying, violent, turbulent fire whirls. The blue whirl cleanly burns heavy, liquid hydrocarbon fuels with no soot production, presenting a new potential way for low-emission combustion. It is reproducible, appears for a range of different fuels and initial conditions, is quiet, appears laminar, and has characteristics which led to the idea that it results from vortex breakdown in whirling, reacting fluid. Since its discovery, considerable effort has been put into measurements, which have shown its temperature structure and sensitivity to the boundary layer near the surface. This has led to considerable speculation about the type of flames that comprise it. Simultaneously, there was a numerical effort to study its structure by performing simulations of vortex breakdown in gaseous reactive flows. The simulations described in this paper show that the stable blue whirl is composed of three different flame structures -- a diffusion flame and a premixed rich and lean flame -- all of which meet in a fourth structure, a triple flame which appears as a whirling blue ring. In addition, the blue whirl structure emerges as the result of vortex breakdown in a swirling reactive flow, as evidenced in the simulation by a bubble mode that is usually invisible in the experiments but at the center of the whirl. This paper also presents the tool used for the study and discusses how this might be used for future investigations.


title: The Doubled Geometry of Nilmanifold Reductions

abstract: A class of special holonomy spaces arise as nilmanifolds fibred over a line interval and are dual to intersecting brane solutions of string theory. Further dualities relate these to T-folds, exotic branes, essentially doubled spaces and spaces with R-flux. We develop the doubled geometry of these spaces, with the various duals arising as different slices of the doubled space.


title: Model Predictive Control-Based Battery Scheduling and Incentives to Manipulate Demand Response Baselines

abstract: We study operations of a battery energy storage system under a baseline-based demand response (DR) program with an uncertain schedule of DR events. Baseline-based DR programs may provide undesired incentives to inflate baseline consumption in non-event days, in order to increase "apparent" DR reduction in event days and secure higher DR payments. Our goal is to identify and quantify such incentives. To understand customer decisions, we formulate the problem of determining hourly battery charging and discharge schedules to minimize expected net costs, defined as energy purchase costs minus energy export rebates and DR payments, over a sufficiently long time horizon (e.g., a year). The complexity of this stochastic optimization problem grows exponentially with the time horizon considered. To obtain computationally tractable solutions, we propose using multistage model predictive control with scenario sampling. Numerical results indicate that our solutions are near optimal (e.g., within 3% from the optimum in the test cases). Finally, we apply our solutions to study an example residential customer with solar photovoltaic and battery systems participating in a typical existing baseline-based DR program. Results reveal that over 66% of the average apparent load reduction during DR events could result from inflation of baseline consumption during non-event days.


title: Popular Differences for Corners in Abelian Groups

abstract: For a compact abelian group $G$, a corner in $G \times G$ is a triple of points $(x,y)$, $(x,y+d)$, $(x+d,y)$. The classical corners theorem of Ajtai and Szemer\'edi implies that for every $\alpha > 0$, there is some $\delta > 0$ such that every subset $A \subset G \times G$ of density $\alpha$ contains a $\delta$ fraction of all corners in $G \times G$, as $x,y,d$ range over $G$. Recently, Mandache proved a "popular differences" version of this result in the finite field case $G = \mathbb F_p^n$, showing that for any subset $A \subset G \times G$ of density $\alpha$, one can fix $d \neq 0$ such that $A$ contains a large fraction, now known to be approximately $\alpha^4$, of all corners with difference $d$, as $x,y$ vary over $G$. We generalize Mandache's result to all compact abelian groups $G$, as well as the case of corners in $\mathbb Z^2$.


title: All Trees are Seven-Cordial

abstract: For any integer $k>0$, a tree $T$ is $k$-cordial if there exists a labeling of the vertices of $T$ by $\mathbb{Z}_k$, inducing edge-weights as the sum modulo $k$ of the labels on incident vertices to a given edge, which furthermore satisfies the following conditions: (i) Each label appears on at most one more vertex than any other label. (ii) Each edge-weight appears on at most one more edge than any other edge-weight. Mark Hovey (1991) conjectured that all trees are $k$-cordial for any integer $k$. Cahit (1987) had shown earlier that all trees are $2$-cordial and Hovey proved that all trees are $3,4,$ and $5$-cordial. Driscoll, et. al. (2017), used an adjustment to Hovey's test to show that all trees are $6$-cordial. It is shown here that all trees are $7$-cordial by that same adjustment.


title: Witten-Veneziano mechanism and pseudoscalar glueball-meson mixing in holographic QCD

abstract: We revisit the $U(1)_A$ anomaly in the holographic model of low-energy QCD by Witten, Sakai, and Sugimoto, presenting a new and direct derivation of the Witten-Veneziano mechanism for generating the mass of the $\eta'$ through an anomalous mixing of the Ramond-Ramond $C_1$ field with the singlet component of the pseudoscalar mesons. The latter turns out to have a kinetic mixing with the normalizable modes of the $C_1$ field representing pseudoscalar glueballs, yielding additional vertices for their production and their decay that dominate over those of the unmixed case considered previously in the Witten-Sakai-Sugimoto model. The leading channel is predicted to be decay into two vector mesons, followed in importance by decay into three pseudoscalar mesons. The issue of production of pseudoscalar glueballs in radiative $J/\psi$ decays and in double diffractive processes is also discussed briefly.


title: Interactive Particle Systems on Hypergraphs, Drift Analysis and the WalkSAT algorithm

abstract: We analyze the expected running time of WalkSAT, a well-known local search procedure for satisfiability solving, on satisfiable instances of the k-XOR SAT problem. We obtain estimates of this expected running time by reducing the problem to a setting amenable to classical techniques from drift analysis. A crucial ingredient of this reduction is the definition of (new, explosive) hypergraph versions of interacting particle systems, notably of coalescing and annihilating random walks as well as the voter model. The use of these tools allows to show that the expected running time of WalkSAT depends on structural parameter (we call odd Cheeger drift) of the dual of the formula hypergraph.


title: Realtime Simulation of Thin-Shell Deformable Materials using CNN-Based Mesh Embedding

abstract: We address the problem of accelerating thin-shell deformable object simulations by dimension reduction. We present a new algorithm to embed a high-dimensional configuration space of deformable objects in a low-dimensional feature space, where the configurations of objects and feature points have approximate one-to-one mapping. Our key technique is a graph-based convolutional neural network (CNN) defined on meshes with arbitrary topologies and a new mesh embedding approach based on physics-inspired loss term. We have applied our approach to accelerate high-resolution thin shell simulations corresponding to cloth-like materials, where the configuration space has tens of thousands of degrees of freedom. We show that our physics-inspired embedding approach leads to higher accuracy compared with prior mesh embedding methods. Finally, we show that the temporal evolution of the mesh in the feature space can also be learned using a recurrent neural network (RNN) leading to fully learnable physics simulators. After training our learned simulator runs $500-10000\times$ faster and the accuracy is high enough for robot manipulation tasks.