title: Enhancing the discovery prospects for SUSY-like decays with a forgotten kinematic variable

abstract: The lack of a new physics signal thus far at the Large Hadron Collider motivates us to consider how to look for challenging final states, with large Standard Model backgrounds and subtle kinematic features, such as cascade decays with compressed spectra. Adopting a benchmark SUSY-like decay topology with a four-body final state proceeding through a sequence of two-body decays via intermediate resonances, we focus our attention on the kinematic variable $\Delta_{4}$ which previously has been used to parameterize the boundary of the allowed four-body phase space. We highlight the advantages of using $\Delta_{4}$ as a discovery variable, and present an analysis suggesting that the pairing of $\Delta_{4}$ with another invariant mass variable leads to a significant improvement over more conventional variable choices and techniques.


title: Kempf-Ness type theorems and Nahm equations

abstract: We prove a version of the affine Kempf-Ness theorem for non-algebraic symplectic structures and shifted moment maps, and use it to describe hyperkahler quotients of T*G, where G is a complex reductive group.


title: A Simple Elementary Proof of P=NP based on the Relational Model of E. F. Codd

abstract: The P versus NP problem is studied under the relational model of E. F. Codd. I found that the term "complete configuration" is unnecessary and harmful in computational complexity theory because of excessive symbol redundancy. For an input, its valid sequences of complete configurations are normalized into a relational model of shared trichoices with no redundancy. To simplify the problem, a polynomial time nondeterministic Turing machine is polynomially reduced to a periodic machine, which only reverses its tape head displacement at the tape ends. By enumerating all the O(p(n)) shared trichoices, a polynomial time p(n) periodic machine is simulated in time O((p(n))^4) under logarithmic cost. A simple elementary proof of P=NP is obtained.


title: Genetic algorithms with DNN-based trainable crossover as an example of partial specialization of general search

abstract: Universal induction relies on some general search procedure that is doomed to be inefficient. One possibility to achieve both generality and efficiency is to specialize this procedure w.r.t. any given narrow task. However, complete specialization that implies direct mapping from the task parameters to solutions (discriminative models) without search is not always possible. In this paper, partial specialization of general search is considered in the form of genetic algorithms (GAs) with a specialized crossover operator. We perform a feasibility study of this idea implementing such an operator in the form of a deep feedforward neural network. GAs with trainable crossover operators are compared with the result of complete specialization, which is also represented as a deep neural network. Experimental results show that specialized GAs can be more efficient than both general GAs and discriminative models.


title: Quantum walks on hypergraphs

abstract: In this work we introduce the concept of a quantum walk on a hypergraph. We show that the staggered quantum walk model is a special case of a quantum walk on a hypergraph.


title: Are Starspots and Plages Co-Located on Active G and K Stars?

abstract: We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation consistent with a bright hemisphere followed by a relatively dark hemisphere, generating photometric variability with a semi-amplitude of 4%. We find that KIC 9652680 is darkest when its $S$-index of Ca II H & K emission is at its maximum. We interpret this anti-correlation between flux and $S$ to indicate that dark starspots in the photosphere are co-located with the bright plages in the chromosphere, as they are on the Sun. Moving to lower masses and slower rotators, we present K2 observations with simultaneous spectroscopy of EPIC 211928486 (K5V) and EPIC 211966629 (K4V), two active stars in the 650 Myr-old open cluster Praesepe. The K2 photometry reveals that both stars have rotation periods of 11.7 days; while their flux varies by 1 and 2% respectively, their Ca II H & K $S$-indices seem to hold relatively constant as a function of rotational phase. This suggests that extended chromospheric networks of plages are not concentrated into regions of emission centered on the starspots that drive rotational modulation, unlike KIC 9652680. We also note that the Ca II emission of EPIC 211928486 dipped and recovered suddenly over the duration of one rotation, suggesting that the evolution timescale of plages may be of order the rotation period.


title: Highly-efficient second and third harmonic generation in a monocrystalline lithium niobate microresonator

abstract: Nonlinear optics in whispering-gallery-mode (WGM) microresonators have attracted much attention. Owing to strong confinement of the light in a small volume, a WGM microresonator can dramatically boost the strength of light field, giving rise to enhancement of the nonlinear interactions of light with the resonator material. Here, we demonstrate highly efficient second harmonic generation (SHG) and third harmonic generation (THG) in an on-chip monocrystalline lithium niobate (LN) microresonator. Benefitting from a cyclic phase matching scheme for the transverse-electric WGMs, nonlinear wavelength conversion utilizing the largest second-order nonlinear coefficient d33, which has not been demonstrated until now despite of its obvious advantage, is successfully achieved in an X-cut LN microresonator. We obtain high conversion efficiencies in not only the SHG (3.8% mW^-1) but also the THG (0.3% mW^-2) realized through a cascaded sum-frequency process. Our results represent a major step toward the classical and quantum photonic integrated circuits.


title: Stochastic and deterministic switches in a bistable polariton micropillar under short optical pulses

abstract: Optical bistability of exciton polaritons in semiconductor microcavities is a promising platform for digital optical devices. Steady states of coherently driven polaritons can be toggled in tens of picoseconds by a short external pulse of appropriate amplitude and phase. We have analyzed the switching behavior of polaritons depending on the pulse amplitude, phase, and duration. The switches are found to change dramatically when the inverse pulse duration becomes comparable to the frequency detuning between the driving field and polariton resonance. If the detuning is large compared to the polariton linewidth, the system becomes extremely sensitive to initial conditions and thus can serve as a fast random-number generator.


title: Label-less Learning for Traffic Control in an Edge Network

abstract: With the development of intelligent applications (e.g., self-driving, real-time emotion recognition, etc), there are higher requirements for the cloud intelligence. However, cloud intelligence depends on the multi-modal data collected by user equipments (UEs). Due to the limited capacity of network bandwidth, offloading all data generated from the UEs to the remote cloud is impractical. Thus, in this article, we consider the challenging issue of achieving a certain level of cloud intelligence while reducing network traffic. In order to solve this problem, we design a traffic control algorithm based on label-less learning on the edge cloud, which is dubbed as LLTC. By the use of the limited computing and storage resources at edge cloud, LLTC evaluates the value of data, which will be offloaded. Specifically, we first give a statement of the problem and the system architecture. Then, we design the LLTC algorithm in detail. Finally, we set up the system testbed. Experimental results show that the proposed LLTC can guarantee the required cloud intelligence while minimizing the amount of data transmission.


title: Measuring monopole and dipole polarizability of acoustic meta-atoms

abstract: We present a method to extract monopole and dipole polarizability from experimental measurements of two-dimensional acoustic meta-atoms. In contrast to extraction from numerical results, this enables all second-order effects and uncertainties in material properties to be accounted for. We apply the technique to 3D-printed labyrinthine meta-atoms of a variety of geometries. We show that the polarizability of structures with shorter acoustic path length agrees well with numerical results. However, those with longer path lengths suffer strong additional damping, which we attribute to the strong viscous and thermal losses in narrow channels.


title: On three-loop RGE for the Higgs sector of 2HDM

abstract: We discuss renormalization group equations (RGE) for the parameters of the Higgs sector in general Two-Higgs-Doublet Model (2HDM). We present the three-loop results but consider only contributions due to self-couplings of the Higgs doublets. We study the structure of RGE and express beta-functions in terms of reparametrization invariants with respect to higgs-basis rotations. The Cayley-Hamilton theorem is utilized to reduce both the number of independent tensor structures in matrix RGE and the number of invariants to a minimal set. As a by-product of our calculation we discovered that two-loop RGE of the scalar sector in general QFT with multiple higgses were not properly implemented in a number of public packages. The latter give rise to a wrong result when mixing in the scalar sector is allowed.


title: On Joint Distributions, Counterfactual Values, and Hidden Variables in Relation to Contextuality

abstract: This paper deals with three traditional ways of defining contextuality: (C1) in terms of (non)existence of certain joint distributions involving measurements made in several mutually exclusive contexts; (C2) in terms of relationship between factual measurements in a given context and counterfactual measurements that could be made if one used other contexts; and (C3) in terms of (non)existence of ``hidden variables'' that determine the outcomes of all factually performed measurements. It is generally believed that the three meanings are equivalent, but the issues involved are not entirely transparent. Thus, arguments have been offered that C2 may have nothing to do with C1, and the traditional formulation of C1 itself encounters difficulties when measurement outcomes in a contextual system are treated as random variables. I show that if C1 is formulated within the framework of the Contextuality-by-Default (CbD) theory, the notion of a probabilistic coupling, the core mathematical tool of CbD, subsumes both counterfactual values and ``hidden variables''. In the latter case, a coupling itself can be viewed as a maximally parsimonious choice of a hidden variable.


title: Image contrast enhancement using fuzzy logic

abstract: Image enhancement is a method of improving the quality of an image and contrast is a major aspect. Traditional methods of contrast enhancement like histogram equalization results in over/under enhancement of the image especially a lower resolution one. This paper aims at developing a new Fuzzy Inference System to enhance the contrast of the low resolution images overcoming the shortcomings of the traditional methods. Results obtained using both the approaches are compared.


title: Steklov Regularization and Trajectory Methods for Univariate Global Optimization

abstract: We introduce a new regularization technique, using what we refer to as the Steklov regularization function, and apply this technique to devise an algorithm that computes a global minimizer of univariate coercive functions. First, we show that the Steklov regularization convexifies a given univariate coercive function. Then, by using the regularization parameter as the independent variable, a trajectory is constructed on the surface generated by the Steklov function. For monic quartic polynomials, we prove that this trajectory does generate a global minimizer. In the process, we derive some properties of quartic polynomials. Comparisons are made with a previous approach which uses a quadratic regularization function. We carry out numerical experiments to illustrate the working of the new method on polynomials of various degree as well as a non-polynomial function.


title: Line shapes of the magnesium resonance lines in cool DZ white dwarf atmospheres

abstract: Line shapes of the magnesium resonance lines in white dwarf spectra are determined by the properties of magnesium atoms and the structure of the white dwarf atmosphere. Through their blanketing effect, these lines have a dominant influence on the model structure and thus on the determination from the spectra of other physical parameters that describe the stellar atmosphere and elemental abundances. In continuation of previous work on Mg$^+$He lines in the UV, we present theoretical profiles of the resonance line of neutral Mg perturbed by He at the extreme density conditions found in the cool largely transparent atmosphere of DZ white dwarfs. We accurately determined the broadening of Mg by He in a unified theory of collisional line profiles using ab initio calculations of MgHe potential energies and transition matrix elements among the singlet electronic states that are involved for the observable spectral lines. We computed the shapes and line parameters of the Mg lines and studied their dependence on helium densities and temperatures. We present results over the full range of temperatures from 4000 to 12000 K needed for input to stellar spectra models. Atmosphere models were constructed for a range of effective temperatures and surface gravities typical for cool DZ white dwarfs. We present synthetic spectra tracing the behavior of the Mg resonance line profiles under the low temperatures and high gas pressures prevalent in these atmospheres. The determination of accurate opacity data of magnesium resonance lines together with an improved atmosphere model code lead to a good fit of cool DZ white dwarf stars. The broadening of spectral lines by helium needs to be understood to accurately determine the H/He and Mg/He abundance ratio in DZ white dwarf atmospheres. We emphasize that no free potential parameters or ad hoc adjustments were used to calculate the line profiles.


title: Characterizing the learning dynamics in extremum seeking

abstract: We consider perturbation-based extremum seeking, which recovers an approximate gradient of an analytically unknown objective function through measurements. Using classical needle variation analysis, we are able to explicitly quantify the recovered gradient in the scalar case. We reveal that it corresponds to an averaged gradient of the objective function, even for very general extremum seeking systems. From this, we create a recursion which represents the learning dynamics along the recovered gradient. These results give rise to the interpretation that extremum seeking actually optimizes a function other than the original one. From this insight, a new perspective on global optimization of functions with local extrema emerges: because the gradient is averaged over a certain time period, local extrema might be evened out in the learning dynamics. Moreover, a multidimensional extension of the scalar results is given.


title: SETI Detection Strategies for Single Dish Radio Telescopes

abstract: Radio Searches for Extra Terrestrial Intelligence aim at detecting artificial transmissions from extra terrestrial communicative civilizations. The lack of prior knowledge concerning these potential transmissions increase the search parameter space. Ground-based single dish radio telescopes offer high sensitivity, but standard data products are limited to power spectral density estimates. To overcome important classical energy detector limitations, two detection strategies based on asynchronous ON and OFF astronomical target observations are proposed. Statistical models are described to enable threshold selection and detection performance assessment.


title: Relaxion Dark Matter

abstract: We highlight a new connection between the Standard Model hierarchy problem and the dark matter sector. The key piece is the relaxion field, which besides scanning the Higgs mass and setting the electroweak scale, also constitutes the observed dark matter abundance of the universe. The relaxation mechanism is realized during inflation, and the necessary friction is provided by particle production. Using this framework we show that the relaxion is a phenomenologically viable dark matter candidate in the keV mass range.


title: A coupled bulk-surface model for cell polarisation

abstract: Several cellular activities, such as directed cell migration, are coordinated by an intricate network of biochemical reactions which lead to a polarised state of the cell, in which cellular symmetry is broken, causing the cell to have a well defined front and back. Recent work on balancing biological complexity with mathematical tractability resulted in the proposal and formulation of a famous minimal model for cell polarisation, known as the wave pinning model. In this study, we present a three-dimensional generalisation of this mathematical framework through the maturing theory of coupled bulk-surface semilinear partial differential equations in which protein compartmentalisation becomes natural. We show how a local perturbation over the surface can trigger propagating reactions, eventually stopped in a stable profile by the interplay with the bulk component. We describe the behavior of the model through asymptotic and local perturbation analysis, in which the role of the geometry is investigated. The bulk-surface finite element method is used to generate numerical simulations over simple and complex geometries, which confirm our analysis, showing pattern formation due to propagation and pinning dynamics. The generality of our mathematical and computational framework allows to study more complex biochemical reactions and biomechanical properties associated with cell polarisation in multi-dimensions.


title: Unpaired Brain MR-to-CT Synthesis using a Structure-Constrained CycleGAN

abstract: The cycleGAN is becoming an influential method in medical image synthesis. However, due to a lack of direct constraints between input and synthetic images, the cycleGAN cannot guarantee structural consistency between these two images, and such consistency is of extreme importance in medical imaging. To overcome this, we propose a structure-constrained cycleGAN for brain MR-to-CT synthesis using unpaired data that defines an extra structure-consistency loss based on the modality independent neighborhood descriptor to constrain structural consistency. Additionally, we use a position-based selection strategy for selecting training images instead of a completely random selection scheme. Experimental results on synthesizing CT images from brain MR images demonstrate that our method is better than the conventional cycleGAN and approximates the cycleGAN trained with paired data.


title: Topological phases in pseudospin-1 Fermi gases with two-dimensional spin-orbit coupling

abstract: The recent experimental realization of spin-orbit (SO) coupling for ultracold bosons and fermions opens an exciting avenue for engineering quantum matter that may be challenging to realize in solid state materials such as SO coupled spin-1 fermions. While one-dimensional SO coupling for spin-1 bosons has been experimentally realized, the generation of two-dimensional (2D) SO coupling and its topological properties are largely unexplored. Here we propose an experimental scheme for realizing a 2D Dresselhaus-type SO coupling in a square lattice for spin-1 Fermi gases. Because of the extended spin degree of freedom, many interesting topological phases could exist without relying on lattice point group symmetries that are crucial in solid state materials. These exotic phases include triply-degenerate points, quadratic band touching, a large Chern number ($C=5$) superfluid with 5 Majorana modes, triple-Weyl fermions, etc. Our scheme can be generalized to larger spins and provides a new route for engineering topological quantum matter by utilizing large spin degrees of freedom, instead of specific lattice symmetries.


title: Poisson brackets after Jacobi and Plucker

abstract: We construct a symplectic realization and a bi-hamiltonian formulation of a 3-dimensional system whose solution are the Jacobi elliptic functions. We generalize this system and the related Poisson brackets to higher dimensions. These more general systems are parametrized by lines in projective space. For these rank 2 Poisson brackets the Jacobi identity is satisfied only when the Pl\" ucker relations hold. Two of these Poisson brackets are compatible if and only if the corresponding lines in projective space intersect. We present several examples of such systems.


title: Frequency-Aware Model Predictive Control

abstract: Transferring solutions found by trajectory optimization to robotic hardware remains a challenging task. When the optimization fully exploits the provided model to perform dynamic tasks, the presence of unmodeled dynamics renders the motion infeasible on the real system. Model errors can be a result of model simplifications, but also naturally arise when deploying the robot in unstructured and nondeterministic environments. Predominantly, compliant contacts and actuator dynamics lead to bandwidth limitations. While classical control methods provide tools to synthesize controllers that are robust to a class of model errors, such a notion is missing in modern trajectory optimization, which is solved in the time domain. We propose frequency-shaped cost functions to achieve robust solutions in the context of optimal control for legged robots. Through simulation and hardware experiments we show that motion plans can be made compatible with bandwidth limits set by actuators and contact dynamics. The smoothness of the model predictive solutions can be continuously tuned without compromising the feasibility of the problem. Experiments with the quadrupedal robot ANYmal, which is driven by highly-compliant series elastic actuators, showed significantly improved tracking performance of the planned motion, torque, and force trajectories and enabled the machine to walk robustly on terrain with unmodeled compliance.


title: Ligand dissociation mechanisms from all-atom simulations: Are we there yet?

abstract: Large parallel gains in the development of both computational resources as well as sampling methods have now made it possible to simulate dissociation events in ligand-protein complexes with all--atom resolution. Such encouraging progress, together with the inherent spatiotemporal resolution associated with molecular simulations, has left their use for investigating dissociation processes brimming with potential, both in rational drug design, where it can be an invaluable tool for determining the mechanistic driving forces behind dissociation rate constants, as well as in force-field development, where it can provide a catalog of transient molecular structures on which to refine force-fields. Although much progress has been made in making force-fields more accurate, reducing their error for transient structures along a transition path could yet prove to be a critical development helping to make kinetic predictions much more accurate. In what follows we will provide a state-of-the-art compilation of the molecular dynamics (MD) methods used to investigate the kinetics and mechanisms of ligand-protein dissociation processes. Due to the timescales of such processes being slower than what is accessible using straightforward MD simulations, several ingenious schemes are being devised at a rapid rate to overcome this obstacle. Here we provide an up-to-date compendium of such methods and their achievements/shortcomings in extracting mechanistic insight into ligand-protein dissociation. We conclude with a critical and provocative appraisal attempting to answer the title of this review.


title: Lugsail lag windows for estimating time-average covariance matrices

abstract: Lag windows are commonly used in time series, econometrics, steady-state simulation, and Markov chain Monte Carlo to estimate time-average covariance matrices. In the presence of positive correlation of the underlying process, estimators of this matrix almost always exhibit significant negative bias, leading to undesirable finite-sample properties. We propose a new family of lag windows specifically designed to improve finite-sample performance by offsetting this negative bias. Any existing lag window can be adapted into a lugsail equivalent with no additional assumptions. We use these lag windows within spectral variance estimators and demonstrate its advantages in a linear regression model with autocorrelated and heteroskedastic residuals. We further employ the lugsail lag windows in weighted batch means estimators due to their computational efficiency on large simulation output. We obtain bias and variance results for these multivariate estimators and significantly weaken the mixing condition on the process. Superior finite-sample properties are illustrated in a vector autoregressive process and a Bayesian logistic regression model.


title: The Inductive Bias of Restricted f-GANs

abstract: Generative adversarial networks are a novel method for statistical inference that have achieved much empirical success; however, the factors contributing to this success remain ill-understood. In this work, we attempt to analyze generative adversarial learning -- that is, statistical inference as the result of a game between a generator and a discriminator -- with the view of understanding how it differs from classical statistical inference solutions such as maximum likelihood inference and the method of moments. Specifically, we provide a theoretical characterization of the distribution inferred by a simple form of generative adversarial learning called restricted f-GANs -- where the discriminator is a function in a given function class, the distribution induced by the generator is restricted to lie in a pre-specified distribution class and the objective is similar to a variational form of the f-divergence. A consequence of our result is that for linear KL-GANs -- that is, when the discriminator is a linear function over some feature space and f corresponds to the KL-divergence -- the distribution induced by the optimal generator is neither the maximum likelihood nor the method of moments solution, but an interesting combination of both.


title: Quantum control of photoelectron circular dichroism

abstract: We demonstrate coherent control over the photoelectron circular dichroism in randomly oriented chiral molecules, based on quantum interference between multiple photoionization pathways. To significantly enhance the chiral signature, we use a finite manifold of indistinguishable $(1+1^\prime)$ REMPI pathways interfering at a common photoelectron energy but probing different intermediate states. We show that this coherent control mechanism maximizes the number of molecular states that constructively contribute to the dichroism at an optimal photoelectron energy and thus outperforms other schemes, including interference between opposite-parity pathways driven by bichromatic $({\omega}, 2{\omega})$ fields as well as sequential pump-probe ionization.


title: Propagation of Moments and Semiclassical Limit from Hartree to Vlasov Equation

abstract: In this paper, we prove a quantitative version of the semiclassical limit from the Hartree to the Vlasov equation with singular interaction, including the Coulomb potential. To reach this objective, we also prove the propagation of velocity moments and weighted Schatten norms which implies the boundedness of the space density of particles uniformly in the Planck constant.


title: Photochemical Oxygen in Non-1 Bar CO2 Atmospheres of Terrestrial Exoplanets

abstract: Atmospheric chemistry models have shown molecular oxygen can build up in CO2-dominated atmospheres on potentially habitable exoplanets without an input of life. Existing models typically assume a surface pressure of 1 bar. Here we present model scenarios of CO2-dominated atmospheres with the surface pressure ranging from 0.1 to 10 bars, while keeping the surface temperature at 288 K. We use a one-dimensional photochemistry model to calculate the abundance of O2 and other key species, for outgassing rates ranging from a Venus-like volcanic activity up to 20x Earth-like activity. The model maintains the redox balance of the atmosphere and the ocean, and includes the pressure dependency of outgassing on the surface pressure. Our calculations show that the surface pressure is a controlling parameter in the photochemical stability and oxygen buildup of CO2-dominated atmospheres. The mixing ratio of O2 monotonically decreases as the surface pressure increases at the very high outgassing rates, whereas it increases as the surface pressure increases at the lower-than-Earth outgassing rates. Abiotic O2 can only build up to the detectable level, defined as 1e-3 in volume mixing ratio, in 10-bar atmospheres with the Venus-like volcanic activity rate and the reduced outgassing rate of H2 due to the high surface pressure. Our results support the search for biological activities and habitability via atmospheric O2 on terrestrial planets in the habitable zone of Sun-like stars.


title: Dedicated symplectic integrators for rotation motions

abstract: We propose to use the properties of the Lie algebra of the angular momentum to build symplectic integrators dedicated to the Hamiltonian of the free rigid body. By introducing a dependence of the coefficients of integrators on the moments of inertia of the integrated body, we can construct symplectic dedicated integrators with fewer stages than in the general case for a splitting in three parts of the Hamiltonian. We perform numerical tests to compare the developed dedicated 4th-order integrators to the existing reference integrators for the water molecule. We also estimate analytically the accuracy of these new integrators for the set of the rigid bodies and conclude that they are more accurate than the existing ones only for very asymmetric bodies.


title: Using the Tsetlin Machine to Learn Human-Interpretable Rules for High-Accuracy Text Categorization with Medical Applications

abstract: Medical applications challenge today's text categorization techniques by demanding both high accuracy and ease-of-interpretation. Although deep learning has provided a leap ahead in accuracy, this leap comes at the sacrifice of interpretability. To address this accuracy-interpretability challenge, we here introduce, for the first time, a text categorization approach that leverages the recently introduced Tsetlin Machine. In all brevity, we represent the terms of a text as propositional variables. From these, we capture categories using simple propositional formulae, such as: if "rash" and "reaction" and "penicillin" then Allergy. The Tsetlin Machine learns these formulae from a labelled text, utilizing conjunctive clauses to represent the particular facets of each category. Indeed, even the absence of terms (negated features) can be used for categorization purposes. Our empirical comparison with Na\"ive Bayes, decision trees, linear support vector machines (SVMs), random forest, long short-term memory (LSTM) neural networks, and other techniques, is quite conclusive. The Tsetlin Machine either performs on par with or outperforms all of the evaluated methods on both the 20 Newsgroups and IMDb datasets, as well as on a non-public clinical dataset. On average, the Tsetlin Machine delivers the best recall and precision scores across the datasets. Finally, our GPU implementation of the Tsetlin Machine executes 5 to 15 times faster than the CPU implementation, depending on the dataset. We thus believe that our novel approach can have a significant impact on a wide range of text analysis applications, forming a promising starting point for deeper natural language understanding with the Tsetlin Machine.


title: Classification of Simple Cuspidal Modules over a Lattice Lie Algebra of Witt type

abstract: Let $W_\pi$ be the lattice Lie algebra of Witt type associated with an additive inclusion $\pi: \mathbb{Z}^N \hookrightarrow \mathbb{C}^2$ with $N>1$. In this article, the classification of simple $\mathbb{Z}^N$-graded $W_\pi$-modules, whose multiplicities are uniformly bounded, is given.


title: Driving Skill Modeling Using Neural Networks for Performance-based Haptic Assistance

abstract: This paper addresses a data-driven framework, modeling expert driving skills for performance-based haptic assistance using neural networks (NNs). We have built a haptic driving training simulator to collect expert driving data and to provide proper haptic feedback. We establish an expert driving skill model by training NNs with the collected data. Then, the skill model is applied to performance-based haptic assistance to provide optimized references of the steering/pedaling movements. We evaluate the skill model and its application to performance-based haptic assistance in two user experiments. The results of the first experiment demonstrate that our skill model has appropriately captured experts' steering/pedaling skills. The results of the second experiment show that our performance-based haptic assistance can help novice drivers perform steering as expert drivers, but cannot assist their pedaling performance.


title: Cross-correlating Carbon Monoxide Line-intensity Maps with Spectroscopic and Photometric Galaxy Surveys

abstract: Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$\alpha$ emitters). We find that obtaining $\sigma_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.


title: Constraining cosmology with the cosmic microwave and infrared backgrounds correlation

abstract: We explore the use of the cosmic infrared background as a tracer of the LSS for cross-correlating with the CMB and exploit the ISW. We use the improved linear CIB model of Maniyar et al (2018) and derive the theoretical CIBxISW cross-correlation for Planck HFI frequencies 217, 353, 545 and 857 GHz and IRAS 3000 GHz. We predict a positive cross-correlation between the CIB and CMB whose amplitude decreases rapidly at small scales. We perform a signal-to-noise ratio (SNR) analysis on this cross-correlation. In the ideal case of the cross-correlation obtained over 70% (40%) of the sky with no residual contaminants (e.g. galactic dust) in maps, the SNR ranges from 4.2-5.6 (3.2-4.3) with the highest for 857 GHz. A Fisher matrix analysis shows that an ISW signal detected with such high SNR on the 40% sky can improve the constraints on the cosmological parameters considerably; constraints on the equation of state of the dark energy are improved by 80%. We then perform a more realistic analysis with the effect of residual galactic dust contamination in CIB maps. We calculate the dust power spectra for several frequencies and sky fractions which dominate over CIB at lower multipoles we are interested in. Considering conservative 10% residual level of galactic dust in the CIB power spectra, we find that the SNR drops drastically making ISW detection difficult. To check the capability of current maps to detect ISW via this method, we measure the cross-correlation of the CIB and the CMB Planck maps on so called GASS field covering an area of 11% in the southern hemisphere. We find that with such a small sky fraction and dust residuals present in CIB maps, we do not detect ISW signal and the measured signal is consistent with zero. In order not to degrade the SNR for the ISW measurement by more than 10% on the 40% sky, we find that the dust needs to be cleaned up to 0.01% level on the power spectrum.


title: Optimal probabilistic storage and retrieval of unitary channels

abstract: We address the question of a quantum memory storage of quantum dynamics. In particular, we design an optimal protocol for $N\to 1$ probabilistic storage-and-retrieval of unitary channels on $d$-dimensional quantum systems. If we may access the unknown unitary gate only $N$-times, the optimal success probability of perfect retrieval of its single use is $N/(N-1+d^2)$. The derived size of the memory system exponentially improves the known upper bound on the size of the program register needed for probabilistic programmable quantum processors. Our results are closely related to probabilistic perfect alignment of reference frames and probabilistic port-based teleportation.


title: End-to-end Audiovisual Speech Activity Detection with Bimodal Recurrent Neural Models

abstract: Speech activity detection (SAD) plays an important role in current speech processing systems, including automatic speech recognition (ASR). SAD is particularly difficult in environments with acoustic noise. A practical solution is to incorporate visual information, increasing the robustness of the SAD approach. An audiovisual system has the advantage of being robust to different speech modes (e.g., whisper speech) or background noise. Recent advances in audiovisual speech processing using deep learning have opened opportunities to capture in a principled way the temporal relationships between acoustic and visual features. This study explores this idea proposing a \emph{bimodal recurrent neural network} (BRNN) framework for SAD. The approach models the temporal dynamic of the sequential audiovisual data, improving the accuracy and robustness of the proposed SAD system. Instead of estimating hand-crafted features, the study investigates an end-to-end training approach, where acoustic and visual features are directly learned from the raw data during training. The experimental evaluation considers a large audiovisual corpus with over 60.8 hours of recordings, collected from 105 speakers. The results demonstrate that the proposed framework leads to absolute improvements up to 1.2% under practical scenarios over a VAD baseline using only audio implemented with deep neural network (DNN). The proposed approach achieves 92.7% F1-score when it is evaluated using the sensors from a portable tablet under noisy acoustic environment, which is only 1.0% lower than the performance obtained under ideal conditions (e.g., clean speech obtained with a high definition camera and a close-talking microphone).


title: Proceedings 5th Workshop on Horn Clauses for Verification and Synthesis

abstract: Many Program Verification and Synthesis problems of interest can be modeled directly using Horn clauses and many recent advances in the CLP and CAV communities have centered around efficiently solving problems presented as Horn clauses. The HCVS series of workshops aims to bring together researchers working in the two communities of Constraint/Logic Programming (e.g., ICLP and CP), Program Verification (e.g., CAV, TACAS, and VMCAI), and Automated Deduction (e.g., CADE, IJCAR), on the topic of Horn clause based analysis, verification, and synthesis. Horn clauses for verification and synthesis have been advocated by these communities in different times and from different perspectives and HCVS is organized to stimulate interaction and a fruitful exchange and integration of experiences.


title: A rapid and well-conditioned algorithm for the Helmholtz--Hodge decomposition of vector fields on the sphere

abstract: A rapid algorithm is derived for the Helmholtz--Hodge decomposition on the surface of the sphere in spherical coordinates. The algorithm uncouples modes of spherical harmonics with different absolute order, writes the conversion as barely-overdetermined banded linear systems, and solves them with banded $QR$ decompositions that factor and execute in optimal complexity. Rigorous upper bounds on the $2$-norm relative condition number of the banded linear systems support the observable low error growth with respect to truncation degree.


title: Unsupervised Controllable Text Formalization

abstract: We propose a novel framework for controllable natural language transformation. Realizing that the requirement of parallel corpus is practically unsustainable for controllable generation tasks, an unsupervised training scheme is introduced. The crux of the framework is a deep neural encoder-decoder that is reinforced with text-transformation knowledge through auxiliary modules (called scorers). The scorers, based on off-the-shelf language processing tools, decide the learning scheme of the encoder-decoder based on its actions. We apply this framework for the text-transformation task of formalizing an input text by improving its readability grade; the degree of required formalization can be controlled by the user at run-time. Experiments on public datasets demonstrate the efficacy of our model towards: (a) transforming a given text to a more formal style, and (b) introducing appropriate amount of formalness in the output text pertaining to the input control. Our code and datasets are released for academic use.


title: Solving Sinhala Language Arithmetic Problems using Neural Networks

abstract: A methodology is presented to solve Arithmetic problems in Sinhala Language using a Neural Network. The system comprises of (a) keyword identification, (b) question identification, (c) mathematical operation identification and is combined using a neural network. Naive Bayes Classification is used in order to identify keywords and Conditional Random Field to identify the question and the operation which should be performed on the identified keywords to achieve the expected result. "One vs. all Classification" is done using a neural network for sentences. All functions are combined through the neural network which builds an equation to solve the problem. The paper compares each methodology in ARIS and Mahoshadha to the method presented in the paper. Mahoshadha2 learns to solve arithmetic problems with the accuracy of 76%.


title: Coordinated Heterogeneous Distributed Perception based on Latent Space Representation

abstract: We investigate a reinforcement approach for distributed sensing based on the latent space derived from multi-modal deep generative models. Our contribution provides insights to the following benefits: Detections can be exchanged effectively between robots equipped with uni-modal sensors due to a shared latent representation of information that is trained by a Variational Auto Encoder (VAE). Sensor-fusion can be applied asynchronously due to the generative feature of the VAE. Deep Q-Networks (DQNs) are trained to minimize uncertainty in latent space by coordinating robots to a Point-of-Interest (PoI) where their sensor modality can provide beneficial information about the PoI. Additionally, we show that the decrease in uncertainty can be defined as the direct reward signal for training the DQN.


title: Benchmarking and Optimization of Gradient Boosting Decision Tree Algorithms

abstract: Gradient boosting decision trees (GBDTs) have seen widespread adoption in academia, industry and competitive data science due to their state-of-the-art performance in many machine learning tasks. One relative downside to these models is the large number of hyper-parameters that they expose to the end-user. To maximize the predictive power of GBDT models, one must either manually tune the hyper-parameters, or utilize automated techniques such as those based on Bayesian optimization. Both of these approaches are time-consuming since they involve repeatably training the model for different sets of hyper-parameters. A number of software GBDT packages have started to offer GPU acceleration which can help to alleviate this problem. In this paper, we consider three such packages: XGBoost, LightGBM and Catboost. Firstly, we evaluate the performance of the GPU acceleration provided by these packages using large-scale datasets with varying shapes, sparsities and learning tasks. Then, we compare the packages in the context of hyper-parameter optimization, both in terms of how quickly each package converges to a good validation score, and in terms of generalization performance.


title: Game-Based Video-Context Dialogue

abstract: Current dialogue systems focus more on textual and speech context knowledge and are usually based on two speakers. Some recent work has investigated static image-based dialogue. However, several real-world human interactions also involve dynamic visual context (similar to videos) as well as dialogue exchanges among multiple speakers. To move closer towards such multimodal conversational skills and visually-situated applications, we introduce a new video-context, many-speaker dialogue dataset based on live-broadcast soccer game videos and chats from Twitch.tv. This challenging testbed allows us to develop visually-grounded dialogue models that should generate relevant temporal and spatial event language from the live video, while also being relevant to the chat history. For strong baselines, we also present several discriminative and generative models, e.g., based on tridirectional attention flow (TriDAF). We evaluate these models via retrieval ranking-recall, automatic phrase-matching metrics, as well as human evaluation studies. We also present dataset analyses, model ablations, and visualizations to understand the contribution of different modalities and model components.


title: Constant Amortized RMR Complexity Deterministic Abortable Mutual Exclusion Algorithm for CC and DSM Models

abstract: The abortable mutual exclusion problem was introduced by Scott and Scherer to meet a need that arises in database and real time systems, where processes sometimes have to abandon their attempt to acquire a mutual exclusion lock to initiate recovery from a potential deadlock or to avoid overshooting a deadline. Algorithms of O(1) RMR complexity have been known for the standard mutual exclusion problem for both the Cache-Coherent (CC) and Distributed Shared Memory (DSM) models of multiprocessors, but whether O(1) RMR complexity is also achievable for abortable mutual exclusion has remained open for the 18 years that this problem has been investigated. Jayanti gives a Theta(log n) worst case RMR complexity solution for both models, where n is the maximum number of processes that execute the algorithm concurrently. Giakouppis and Woelfel's algorithm, presented at PODC last year, is an O(1) amortized complexity algorithm, but it works only for the CC model, uses randomization, does not satisfy Starvation Freedom, and the O(1) amortized bound holds only in expectation and is proven for the a weak (oblivious) adversary model. We design an algorithm that is free of these limitations: our algorithm is deterministic, supports fast aborts (a process completes an abort in O(1) steps), has a small space complexity of O(n), requires hardware support for only the Fetch&Store instruction, satisfies a novely defined First Come First Served for abortable locks, and most importantly, has O(1) amortized RMR complexity for both the CC and DSM models. Our algorithm is short and practical with fewer than a dozen lines of code, and is accompanied by a rigorous proof of mutual exclusion through invariants and of starvation-freedom and complexity analysis through distance and potential functions. Thus, modulo amortization, our result answers affirmatively the long standing open question described above.


title: Signatures of moir\'e-trapped valley excitons in MoSe$_2$/WSe$_2$ heterobilayers

abstract: The creation of moir\'e patterns in crystalline solids is a powerful approach to manipulate their electronic properties, which are fundamentally influenced by periodic potential landscapes. In 2D materials, a moir\'e pattern with a superlattice potential can form by vertically stacking two layered materials with a twist and/or finite lattice constant difference. This unique approach has led to emergent electronic phenomena, including the fractal quantum Hall effect, tunable Mott insulators, and unconventional superconductivity. Furthermore, theory predicts intriguing effects on optical excitations by a moir\'e potential in 2D valley semiconductors, but these signatures have yet to be experimentally detected. Here, we report experimental evidence of interlayer valley excitons trapped in a moir\'e potential in MoSe$_2$/WSe$_2$ heterobilayers. At low temperatures, we observe photoluminescence near the free interlayer exciton energy but with over 100 times narrower linewidths. The emitter g-factors are homogeneous across the same sample and only take two values, -15.9 and 6.7, in samples with twisting angles near 60{\deg} and 0\deg, respectively. The g-factors match those of the free interlayer exciton, which is determined by one of two possible valley pairing configurations. At a twist angle near 20\deg, the emitters become two orders of magnitude dimmer, but remarkably, they possess the same g-factor as the heterobilayer near 60\deg. This is consistent with the Umklapp recombination of interlayer excitons near the commensurate 21.8{\deg} twist angle. The emitters exhibit strong circular polarization, which implies the preservation of three-fold rotation symmetry by the trapping potential. Together with the power and excitation energy dependence, all evidence points to their origin as interlayer excitons trapped in a smooth moir\'e potential with inherited valley-contrasting physics.


title: Finding Higher Order Mutants Using Variational Execution

abstract: Mutation testing is an effective but time consuming method for gauging the quality of a test suite. It functions by repeatedly making changes, called mutants, to the source code and checking whether the test suite fails (i.e., whether the mutant is killed). Recent work has shown cases in which applying multiple changes, called a higher order mutation, is more difficult to kill than a single change, called a first order mutation. Specifically, a special kind of higher order mutation, called a strongly subsuming higher order mutation (SSHOM), can enable equivalent accuracy in assessing the quality of the test suite with fewer executions of tests. Little is known about these SSHOMs, as they are difficult to find. Our goal in this research is to identify a faster, more reliable method for finding SSHOMs in order to characterize them in the future. We propose an approach based on variational execution to find SSHOMs. Preliminary results indicate that variational execution performs better than the existing genetic algorithm in terms of speed and completeness of results. Out of a set of 33 first order mutations, our variational execution approach finds all 38 SSHOMs in 4.5 seconds, whereas the genetic algorithm only finds 36 of the 38 SSHOMs in 50 seconds.


title: On the Generalization of Stochastic Gradient Descent with Momentum

abstract: While momentum-based accelerated variants of stochastic gradient descent (SGD) are widely used when training machine learning models, there is little theoretical understanding on the generalization error of such methods. In this work, we first show that there exists a convex loss function for which the stability gap for multiple epochs of SGD with standard heavy-ball momentum (SGDM) becomes unbounded. Then, for smooth Lipschitz loss functions, we analyze a modified momentum-based update rule, i.e., SGD with early momentum (SGDEM) under a broad range of step-sizes, and show that it can train machine learning models for multiple epochs with a guarantee for generalization. Finally, for the special case of strongly convex loss functions, we find a range of momentum such that multiple epochs of standard SGDM, as a special form of SGDEM, also generalizes. Extending our results on generalization, we also develop an upper bound on the expected true risk, in terms of the number of training steps, sample size, and momentum. Our experimental evaluations verify the consistency between the numerical results and our theoretical bounds. SGDEM improves the generalization error of SGDM when training ResNet-18 on ImageNet in practical distributed settings.


title: Optimization-Based Bound Tightening using a Strengthened QC-Relaxation of the Optimal Power Flow Problem

abstract: This article develops a strengthened convex quadratic convex (QC) relaxation of the AC Optimal Power Flow (AC-OPF) problem and presents an optimization-based bound-tightening (OBBT) algorithm to compute tight, feasible bounds on the voltage magnitude variables for each bus and the phase angle difference variables for each branch in the network. Theoretical properties of the strengthened QC relaxation that show its dominance over the other variants of the QC relaxation studied in the literature are also derived. The effectiveness of the strengthened QC relaxation is corroborated via extensive numerical results on benchmark AC-OPF test networks. In particular, the results demonstrate that the proposed relaxation consistently provides the tightest variable bounds and optimality gaps with negligible impacts on runtime performance.


title: Streaming enhanced flow-mediated transport

abstract: We investigate the capability of an active body (master) to manipulate a passive object (slave) purely via contactless flow-mediated mechanisms, motivated by potential applications in microfluidic devices and medicine (drug delivery purposes). We extend prior works on active-passive cylinder pairs by superimposing periodic oscillations to the master's linear motion. In a viscous fluid, such oscillations produce an additional viscous streaming field, which is leveraged for enhancing slave transport. We see that superimposing oscillations robustly improves transport across a range of Reynolds numbers. Comparison with results without oscillations highlights the flow mechanisms at work, which we capitalize on to design (master) geometries for augmented transport. These principles are found to extend to three-dimensional active-passive shapes as well.


title: Stable and Unstable Vortex Knots in Excitable Media

abstract: We study the dynamics of knotted vortices in a bulk excitable medium using the FitzHugh-Nagumo model. From a systematic survey of all knots of at most eight crossings we establish that the generic behaviour is of unsteady, irregular dynamics, with prolonged periods of expansion of parts of the vortex. The mechanism for the length expansion is a long-range `wave slapping' interaction, analogous to that responsible for the annihilation of small vortex rings by larger ones. We also show that there are stable vortex geometries for certain knots; in addition to the unknot, trefoil and figure eight knots reported previously, we have found stable examples of the Whitehead link and $6_2$ knot. We give a thorough characterisation of their geometry and steady state motion. For the unknot, trefoil and figure eight knots we greatly expand previous evidence that FitzHugh-Nagumo dynamics untangles initially complex geometries while preserving topology.


title: A computational model for exploring particle acceleration during reconnection in macro-scale systems

abstract: A new computational model is presented suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfven waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.


title: Multicellular actomyosin cables in epithelia under external anisotropic stress

abstract: The alignment of cell-cell junctions and associated cortical actomyosin across multiple cells to form supracellular cables in an epithelium is an example of the long range tissue organization that drives morphogenesis. Here we demonstrate that the ability of tissues to assemble these parallel cables depends on the initial packing topology of the cells in the epithelium. Using a computational vertex model we develop two methods of measuring a disordered tissue's favorability to forming cables under an external stress. These measures quantify the deformation of cells and the distribution of tension in the tissue under stress. Using these measures we show that passive stress-induced cell flow reduces a tissue's ability to form cables, whereas oriented divisions create a packing which can sustain multiple parallel cables. These measures are applied to a region of the the Drosophila demonstrating a shift to a more cable-friendly packing after a wave of oriented divisions in the region.


title: FINN-R: An End-to-End Deep-Learning Framework for Fast Exploration of Quantized Neural Networks

abstract: Convolutional Neural Networks have rapidly become the most successful machine learning algorithm, enabling ubiquitous machine vision and intelligent decisions on even embedded computing-systems. While the underlying arithmetic is structurally simple, compute and memory requirements are challenging. One of the promising opportunities is leveraging reduced-precision representations for inputs, activations and model parameters. The resulting scalability in performance, power efficiency and storage footprint provides interesting design compromises in exchange for a small reduction in accuracy. FPGAs are ideal for exploiting low-precision inference engines leveraging custom precisions to achieve the required numerical accuracy for a given application. In this article, we describe the second generation of the FINN framework, an end-to-end tool which enables design space exploration and automates the creation of fully customized inference engines on FPGAs. Given a neural network description, the tool optimizes for given platforms, design targets and a specific precision. We introduce formalizations of resource cost functions and performance predictions, and elaborate on the optimization algorithms. Finally, we evaluate a selection of reduced precision neural networks ranging from CIFAR-10 classifiers to YOLO-based object detection on a range of platforms including PYNQ and AWS\,F1, demonstrating new unprecedented measured throughput at 50TOp/s on AWS-F1 and 5TOp/s on embedded devices.


title: Efficient generation of random derangements with the expected distribution of cycle lengths

abstract: We show how to generate random derangements efficiently by two different techniques: random restricted transpositions and sequential importance sampling. The algorithm employing restricted transpositions can also be used to generate random fixed-point-free involutions only, a.k.a. random perfect matchings on the complete graph. Our data indicate that the algorithms generate random samples with the expected distribution of cycle lengths, which we derive, and for relatively small samples, which can actually be very large in absolute numbers, we argue that they generate samples indistinguishable from the uniform distribution. Both algorithms are simple to understand and implement and possess a performance comparable to or better than those of currently known methods. Simulations suggest that the mixing time of the algorithm based on random restricted transpositions (in the total variance distance with respect to the distribution of cycle lengths) is $O(n^{a}\log{n}^{2})$ with $a \simeq \frac{1}{2}$ and $n$ the length of the derangement. We prove that the sequential importance sampling algorithm generates random derangements in $O(n)$ time with probability $O(1/n)$ of failing.


title: Edge universality of separable covariance matrices

abstract: In this paper, we prove the edge universality of largest eigenvalues for separable covariance matrices of the form $\mathcal Q :=A^{1/2}XBX^*A^{1/2}$. Here $X=(x_{ij})$ is an $n\times N$ random matrix with $x_{ij}=N^{-1/2}q_{ij}$, where $q_{ij}$ are $i.i.d.$ random variables with zero mean and unit variance, and $A$ and $B$ are respectively $n \times n$ and $N\times N$ deterministic non-negative definite symmetric (or Hermitian) matrices. We consider the high-dimensional case, i.e. ${n}/{N}\to d \in (0, \infty)$ as $N\to \infty$. Assuming $\mathbb E q_{ij}^3=0$ and some mild conditions on $A$ and $B$, we prove that the limiting distribution of the largest eigenvalue of $\mathcal Q$ coincide with that of the corresponding Gaussian ensemble (i.e. the $\mathcal Q$ with $X$ being an $i.i.d.$ Gaussian matrix) as long as we have $\lim_{s \rightarrow \infty}s^4 \mathbb{P}(\vert q_{ij} \vert \geq s)=0$, which is a sharp moment condition for edge universality. If we take $B=I$, then $\mathcal Q$ becomes the normal sample covariance matrix and the edge universality holds true without the vanishing third moment condition. So far, this is the strongest edge universality result for sample covariance matrices with correlated data (i.e. non-diagonal $A$) and heavy tails, which improves the previous results in \cite{BPZ1,LS} (assuming high moments and diagonal $A$), \cite{Anisotropic} (assuming high moments) and \cite{DY} (assuming diagonal $A$).


title: Design and Performance of an Interferometric Trigger Array for Radio Detection of High-Energy Neutrinos

abstract: Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.


title: From Topological to Quantum Entanglement

abstract: Entanglement is a special feature of the quantum world that reflects the existence of subtle, often non-local, correlations between local degrees of freedom. In topological theories such non-local correlations can be given a very intuitive interpretation: quantum entanglement of subsystems means that there are "strings" connecting them. More generally, an entangled state, or similarly, the density matrix of a mixed state, can be represented by cobordisms of topological spaces. Using a formal mathematical definition of TQFT we construct basic examples of entangled states and compute their von Neumann entropy.


title: Ferroelectric switching of a two-dimensional metal

abstract: A ferroelectric is a material with a polar structure whose polarity can be reversed by applying an electric field. In metals, the itinerant electrons tend to screen electrostatic forces between ions, helping to explain why polar metals are very rare. Screening also excludes external electric fields, apparently ruling out the possibility of polarity reversal and thus ferroelectric switching. In principle, however, a thin enough polar metal could be penetrated by an electric field sufficiently to be switched. Here we show that the layered topological semimetal WTe2 provides the first embodiment of this principle. Although monolayer WTe2 is centrosymmetric and thus nonpolar, the stacked bulk structure is polar. We find that two- or three-layer WTe2 exhibits a spontaneous out-of-plane electric polarization which can be switched using gate electrodes. We directly detect and quantify the polarization using graphene as an electric field sensor. Moreover, the polarization states can be differentiated by conductivity, and the carrier density can be varied to modify the properties. The critical temperature is above 350 K, and even when WTe2 is sandwiched in graphene it retains its switching capability at room temperature, demonstrating a robustness suitable for applications in combination with other two-dimensional materials.


title: Rainbow numbers for $x_1+x_2=kx_3$ in $\mathbb{Z}_n$

abstract: In this work, we investigate the fewest number of colors needed to guarantee a rainbow solution to the equation $x_1 + x_2 = k x_3$ in $\mathbb{Z}_n$. This value is called the Rainbow number and is denoted by $rb(\mathbb{Z}_n, k)$ for positive integer values of $n$ and $k$. We find that $rb(\mathbb{Z}_p, 1) = 4$ for all primes greater than $3$ and that $rb(\mathbb{Z}_n, 1)$ can be deterimined from the prime factorization of $n$. Furthermore, when $k$ is prime, $rb(\mathbb{Z}_n, k)$ can be determined from the prime factorization of $n$.


title: A Short Note on a Weighted Friedrichs Inequality

abstract: In this short note we derive, for bounded domains, an upper bound for a Friedrichs type constant in a weighted Friedrichs type inequality. This upper bound generalizes a well known upper bound of the Friedrichs constant. This upper bound is also used to improve an upper bound of a Maxwell type constant for convex domains in $\mathbb{R}^3$. A simple numerical application is also given: we apply the main result in a posteriori error estimation for an elliptic problem.


title: Simplicity Creates Inequity: Implications for Fairness, Stereotypes, and Interpretability

abstract: Algorithms are increasingly used to aid, or in some cases supplant, human decision-making, particularly for decisions that hinge on predictions. As a result, two additional features in addition to prediction quality have generated interest: (i) to facilitate human interaction and understanding with these algorithms, we desire prediction functions that are in some fashion simple or interpretable; and (ii) because they influence consequential decisions, we also want them to produce equitable allocations. We develop a formal model to explore the relationship between the demands of simplicity and equity. Although the two concepts appear to be motivated by qualitatively distinct goals, we show a fundamental inconsistency between them. Specifically, we formalize a general framework for producing simple prediction functions, and in this framework we establish two basic results. First, every simple prediction function is strictly improvable: there exists a more complex prediction function that is both strictly more efficient and also strictly more equitable. Put another way, using a simple prediction function both reduces utility for disadvantaged groups and reduces overall welfare relative to other options. Second, we show that simple prediction functions necessarily create incentives to use information about individuals' membership in a disadvantaged group --- incentives that weren't present before simplification, and that work against these individuals. Thus, simplicity transforms disadvantage into bias against the disadvantaged group. Our results are not only about algorithms but about any process that produces simple models, and as such they connect to the psychology of stereotypes and to an earlier economics literature on statistical discrimination.


title: Privacy-preserving mHealth Data Release with Pattern Consistency

abstract: Mobile healthcare system integrating wearable sensing and wireless communication technologies continuously monitors the users' health status. However, the mHealth system raises a severe privacy concern as the data it collects are private information, such as heart rate and blood pressure. In this paper, we propose an efficient and privacy-preserving mHealth data release approach for the statistic data with the objectives to preserve the unique patterns in the original data bins. The proposed approach adopts the bucket partition algorithm and the differential privacy algorithm for privacy preservation. A customized bucket partition algorithm is proposed to combine the database value bins into buckets according to certain conditions and parameters such that the patterns are preserved. The differential privacy algorithm is then applied to the buckets to prevent an attacker from being able to identify the small changes at the original data. We prove that the proposed approach achieves differential privacy. We also show the accuracy of the proposed approach through extensive simulations on real data. Real experiments show that our partitioning algorithm outperforms the state-of-the-art in preserving the patterns of the original data by a factor of 1.75.


title: Distorting an Adversary's View in Cyber-Physical Systems

abstract: In Cyber-Physical Systems (CPSs), inference based on communicated data is of critical significance as it can be used to manipulate or damage the control operations by adversaries. This calls for efficient mechanisms for secure transmission of data since control systems are becoming increasingly distributed over larger geographical areas. Distortion based security, recently proposed as one candidate for CPSs security, is not only more appropriate for these applications but also quite frugal in terms of prior requirements on shared keys. In this paper, we propose distortion-based metrics to protect CPSs communication and show that it is possible to confuse adversaries with just a few bits of pre-shared keys.


title: Going Viral: Stability of Consensus-Driven Adoptive Spread

abstract: The spread of new products in a networked population is often modeled as an epidemic. However, in the case of `complex' contagion, these models {do not capture nuanced, dynamic social reinforcement effects in adoption behavior}. In this paper, we investigate a model of complex contagion which allows a coevolutionary interplay between adoption, modeled as an SIS epidemic spreading process, and social reinforcement effects, modeled as consensus opinion dynamics. Asymptotic stability analysis of the all-adopt as well as the none-adopt equilibria of the combined opinion-adoption model is provided through the use of Lyapunov arguments. In doing so, sufficient conditions are provided which determine the stability of the `flop' state, where no one adopts the product and everyone's opinion of the product is least favorable, and the `hit' state, where everyone adopts and their opinions are most favorable. These conditions are shown to extend to the bounded confidence opinion dynamic under a stronger assumption on the model parameters. To conclude, numerical simulations demonstrate behavior of the model which reflect findings from the sociology literature on adoption behavior.


title: Elastoresistive and Elastocaloric Anomalies at Magnetic and Electronic-Nematic Critical Points

abstract: Using Ba(Fe$_{0.975}$Co$_{0.025}$)$_2$As$_{2}$ as an exemplar material exhibiting second order electronic-nematic and antiferromagnetic transitions, we present measurements that reveal anomalies in the elastoresistance $\left(\frac{\partial \rho_{ij}}{\partial\varepsilon_{kl}}\right)$ and elastocaloric effect $\left(\frac{\partial T}{\partial\varepsilon_{kl}}\right)$ at both phase transitions. Both effects are understood to arise from the effect of strain on the transition temperatures; in the region close to the phase transitions this leads to (1) similarity between the strain and temperature derivatives of the resistivity and (2) similarity between the elastocaloric effect and the singular part of the specific heat. These mechanisms for elastoresistance and elastocaloric effect should be anticipated for any material in which mechanical deformation changes the transition temperature. Furthermore, these measurements provide evidence that the Fisher-Langer relation $\rho^{(c)} \propto U^{(c)}$ between the scattering from critical degrees of freedom and their energy-density, respectively, holds near each of the Ne\'{e}l and electronic nematic transitions in the material studied.


title: Efficient and Privacy-preserving Voice-based Search over mHealth Data

abstract: In-home IoT devices play a major role in healthcare systems as smart personal assistants. They usually come with a voice-enabled feature to add an extra level of usability and convenience to elderly, disabled people, and patients. In this paper, we propose an efficient and privacy-preserving voice-based search scheme to enhance the efficiency and the privacy of in-home healthcare applications. We consider an application scenario where patients use the devices to record and upload their voice to servers and the caregivers search the interested voices of their patient's based on the voice content, mood, tone and background sound. Our scheme preserves the richness and privacy of voice data and enables accurate and efficient voice-based search, while in current systems that use speech recognition the richness and privacy of voice data are compromised. Specifically, our scheme achieves the privacy by employing a homomorphic encryption; only encrypted voice data is uploaded to the server who is unable to access the original voice data. In addition, our scheme enables the server to selectively and accurately respond to caregiver's queries on the voice data based on voice's feature similarity. We evaluate our scheme through real experiments and show that our scheme even with privacy preservation can successfully match similar voice data at an average accuracy of 80.8%.


title: Gate-Based High Fidelity Spin Read-out in a CMOS Device

abstract: The engineering of electron spin qubits in a compact unit cell embedding all quantum functionalities is mandatory for large scale integration. In particular, the development of a high-fidelity and scalable spin readout method remains an open challenge. Here we demonstrate high-fidelity and robust spin readout based on gate reflectometry in a CMOS device comprising one qubit dot and one ancillary dot coupled to an electron reservoir to perform readout. This scalable method allows us to read out a spin with a fidelity above 99% for 1 ms integration time. To achieve such fidelity, we exploit a latched spin blockade mechanism that requires electron exchange between the ancillary dot and the reservoir. We show that the demonstrated high read-out fidelity is fully preserved up to 0.5 K. This results holds particular relevance for the future co-integration of spin qubits and classical control electronics.


title: Closed-Book Training to Improve Summarization Encoder Memory

abstract: A good neural sequence-to-sequence summarization model should have a strong encoder that can distill and memorize the important information from long input texts so that the decoder can generate salient summaries based on the encoder's memory. In this paper, we aim to improve the memorization capabilities of the encoder of a pointer-generator model by adding an additional 'closed-book' decoder without attention and pointer mechanisms. Such a decoder forces the encoder to be more selective in the information encoded in its memory state because the decoder can't rely on the extra information provided by the attention and possibly copy modules, and hence improves the entire model. On the CNN/Daily Mail dataset, our 2-decoder model outperforms the baseline significantly in terms of ROUGE and METEOR metrics, for both cross-entropy and reinforced setups (and on human evaluation). Moreover, our model also achieves higher scores in a test-only DUC-2002 generalizability setup. We further present a memory ability test, two saliency metrics, as well as several sanity-check ablations (based on fixed-encoder, gradient-flow cut, and model capacity) to prove that the encoder of our 2-decoder model does in fact learn stronger memory representations than the baseline encoder.


title: The Bernstein problem for Lipschitz intrinsic graphs in the Heisenberg group

abstract: We prove that, in the first Heisenberg group $\mathbb{H}$, an entire locally Lipschitz intrinsic graph admitting vanishing first variation of its sub-Riemannian area and non-negative second variation must be an intrinsic plane, i.e., a coset of a two dimensional subgroup of $\mathbb{H}$. Moreover two examples are given for stressing result's sharpness.


title: Distributed Chernoff Test: Optimal decision systems over networks

abstract: We study "active" decision making over sensor networks where the sensors' sequential probing actions are actively chosen by continuously learning from past observations. We consider two network settings: with and without central coordination. In the first case, the network nodes interact with each other through a central entity, which plays the role of a fusion center. In the second case, the network nodes interact in a fully distributed fashion. In both of these scenarios, we propose sequential and adaptive hypothesis tests extending the classic Chernoff test. We compare the performance of the proposed tests to the optimal sequential test. In the presence of a fusion center, our test achieves the same asymptotic optimality of the Chernoff test, minimizing the risk, expressed by the expected cost required to reach a decision plus the expected cost of making a wrong decision, when the observation cost per unit time tends to zero. The test is also asymptotically optimal in the higher moments of the time required to reach a decision. Additionally, the test is parsimonious in terms of communications, and the expected number of channel uses per network node tends to a small constant. In the distributed setup, our test achieves the same asymptotic optimality of Chernoff's test, up to a multiplicative constant in terms of both risk and the higher moments of the decision time. Additionally, the test is parsimonious in terms of communications in comparison to state-of-the-art schemes proposed in the literature. The analysis of these tests is also extended to account for message quantization and communication over channels with random erasures.

categories: stat.ME cs.IT cs.MA math.IT stat.AP stat.ML


title: Closed geodesics on connected sums and 3-manifolds

abstract: We study the asymptotics of the number N(t) of geometrically distinct closed geodesics of a Riemannian or Finsler metric on a connected sum of two compact manifolds of dimension at least three with non-trivial fundamental groups and apply this result to the prime decomposition of a three-manifold. In particular we show that the function N(t) grows at least like the prime numbers on a compact 3-manifold with infinite fundamental group. It follows that a generic Riemannian metric on a compact 3-manifold has infinitely many geometrically distinct closed geodesics. We also consider the case of a connected sum of a compact manifold with positive first Betti number and a simply-connected manifold which is not homeomorphic to a sphere.


title: Transition between two metallic ferroelectric orders in multiferroic Ca$_3$Ru$_2$O$_7$, induced by magnetism-mediated orbital re-polarization

abstract: For the past decades, the low-temperature phase of Ca$_3$Ru$_2$O$_7$ below the 48K first-order phase transition remains a puzzle with controversial suggestions involving metallic ferroelectric, orbital or magnetic ordering. Through analysis of experimental lattice structure, density functional theory calculation, and effective model analysis, we propose that the 48K phase transition is a bond formation transition promoted by the magnetism mediated orbital re-polarization. Most interestingly, this transition is accompanied by a switch of \textit{two} metallic ferroelectric orders from a $xy+y$ symmetry to $xz+z$. Our study not only resolves a long-standing puzzle of this phase transition in this material, but also demonstrates perhaps the first example of transition between multiple emergent ferroelectric orders in bad metals, resulting from interplay between multiferroic orders.


title: Topological Transformation Monoids

abstract: We investigate semigroup topologies on the full transformation monoid T(X) of an infinite set X. We show that the standard pointwise topology is the weakest Hausdorff semigroup topology on T(X), show that the pointwise topology is the unique Hausdorff semigroup topology on T(X) that induces the pointwise topology on the group of all permutations of X, and construct |X| distinct Hausdorff semigroup topologies on T(X). In the case where X is countable, we prove that the pointwise topology is the only Polish semigroup topology on T(X). We also show that every separable semigroup topology on T(X) is perfect, describe the compact sets in an arbitrary Hausdorff semigroup topology on T(X), and show that there are no locally compact perfect Hausdorff semigroup topologies on T(X) when |X| has uncountable cofinality.


title: On a Diophantine equation with five prime variables

abstract: Let $[x]$ denote the integral part of the real number $x$, and $N$ be a sufficiently large integer. In this paper, it is proved that, for $1<c<\frac{4109054}{1999527}, c\not=2$, the Diophantine equation $N=[p_1^c]+[p_2^c]+[p_3^c]+[p_4^c]+[p_5^c]$ is solvable in prime variables $p_1,p_2,p_3,p_4,p_5$.


title: Theory of bifurcation amplifiers utilizing the nonlinear dynamical response of an optically damped mechanical oscillator

abstract: We consider a standard optomechanical system where a mechanical oscillator is coupled to a cavity mode through the radiation pressure interaction. The oscillator is coherently driven at its resonance frequency, whereas the cavity mode is driven below its resonance, providing optical damping of the mechanical oscillations. We study the nonlinear coherent response of the mechanical oscillator in this setup. For large mechanical amplitudes, we find that the system can display dynamical multistability if the optomechanical cooperativity exceeds a critical value. This analysis relates standard optomechanical damping to the dynamical attractors known from the theory of optomechanical self-sustained oscillations. We also investigate the effect of thermal and quantum noise and estimate the noise-induced switching rate between the stable states of the system. We then consider applications of this system and primarily focus on how it can be used as bifurcation amplifiers for the detection of small mechanical or optical signals. Finally, we show that in a related but more complicated setup featuring resonant optomechanical interactions, the same effects can be realized with a relaxed requirement on the size of the mechanical oscillations.


title: Broadening of Ly alpha by neutral helium in DBA white dwarfs

abstract: Traces of photospheric hydrogen are detected in at least half of all white dwarfs with helium-dominated atmospheres through the presence of H alpha in high-quality optical spectroscopy. Previous studies have noted significant discrepancies between the hydrogen abundances derived from H alpha and Ly alpha for a number of stars where ultraviolet spectroscopy is also available. We demonstrate that this discrepancy is caused by inadequate treatment of the broadening of Ly alpha by neutral helium. When fitting Hubble Space Telescope COS spectroscopy of 17 DB white dwarfs using our new line profile calculations, we find good agreement between log(H/He) measured from Ly alpha and H alpha. Larger values of log(H/He) based on Ly alpha are still found for three stars, which are among the most distant in our sample, and we show that a small amount of interstellar absorption from neutral hydrogen can account for this discrepancy.


title: A Comparison of the Effects of Neutron and Gamma Radiation in Silicon Photomultipliers

abstract: The effects of radiation damage in silicon photomultipliers (SiPMs) from gamma rays have been measured and compared with the damage produced by neutrons. Several types of MPPCs from Hamamatsu were exposed to gamma rays and neutrons at the Solid State Gamma Ray Irradiation Facility (SSGRIF) at Brookhaven National Lab and the Institute for Nuclear Research (Atomki) in Debrecen, Hungary. The gamma ray exposures ranged from 1 krad to 1 Mrad and the neutron exposures ranged from 10$^8$ n/cm$^2$ to 10$^{12}$ n/cm$^2$. The main effect of gamma ray damage is an increase in the noise and leakage current in the irradiated devices, similar to what is seen from neutron damage, but the level of damage is considerably less at comparable high levels of exposure. In addition, the damage from gamma rays saturates after a few hundred krad, while the damage from neutrons shows no sign of saturation, suggestive of different damage mechanisms in the two cases. The change in optical absorption in the window material of the SiPMs due to radiation was also measured. This study was carried out in order to evaluate the use of SiPMs for particle physics applications with moderate levels of radiation exposures.


title: Quantum criticality in Ising chains with random hyperuniform couplings

abstract: We study quantum phase transitions in transverse-field Ising spin chains in which the couplings are random but hyperuniform, in the sense that their large-scale fluctuations are suppressed. We construct a one-parameter family of disorder models in which long-wavelength fluctuations are increasingly suppressed as a parameter $\alpha$ is tuned. For $\alpha = 0$, one recovers the familiar infinite-randomness critical point. For $0 < \alpha < 1$, we find a line of infinite-randomness critical points with continuously varying critical exponents; however, the Griffiths phases that flank the critical point at $\alpha = 0$ are absent at any $\alpha > 0$. When $\alpha > 1$, randomness is a dangerously irrelevant perturbation at the clean Ising critical point, leading to a state we call the critical Ising insulator. In this state, thermodynamics and equilibrium correlation functions behave as in the clean system. However, all finite-energy excitations are localized, thermal transport vanishes, and autocorrelation functions remain finite in the long-time limit. We characterize this line of hyperuniform critical points using a combination of perturbation theory, renormalization-group methods, and exact diagonalization.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.quant-gas cond-mat.str-el quant-ph


title: Very massive stars in not so massive clusters

abstract: Very young star clusters in the Milky Way exhibit a well-defined relation between their maximum stellar mass, m_max, and their mass in stars, M_ecl. A recent study shows that the young intermediate-mass star cluster VVV CL041 possibly hosts a >80 Msun star, WR62-2, which appears to violate the existence of the m_max-M_ecl relation since the mass of the star is almost two times higher than that expected from the relation. By performing direct N-body calculations with the same mass as the cluster VVV CL041 (3000 Msun), we study whether such a very massive star can be formed via dynamically induced stellar collisions in a binary-rich star cluster that initially follows the m_max-M_ecl relation. Eight out of 100 clusters form a star more massive than 80 Msun through multiple stellar collisions. This suggests that the VVV CL041 cluster may have become an outlier of the relation because of its early-dynamical evolution, even if the cluster followed the relation at birth. We find that more than half of our model clusters host a merger product as its most massive member within the first 5 Myr of cluster evolution. Thus, the existence of stars more massive than the m_max-M_ecl relation in some young clusters is expected due to dynamical processes despite the validity of the m_max-M_ecl relation. We briefly discuss evolution of binary populations in our model.


title: Neutron Star Mergers Are the Dominant Source of the r-process in the Early Evolution of Dwarf Galaxies

abstract: There are many candidate sites of the r-process: core-collapse supernovae (including rare magnetorotational core-collapse supernovae), neutron star mergers, and neutron star/black hole mergers. The chemical enrichment of galaxies---specifically dwarf galaxies---helps distinguish between these sources based on the continual build-up of r-process elements. This technique can distinguish between the r-process candidate sites by the clearest observational difference---how quickly these events occur after the stars are created. The existence of several nearby dwarf galaxies allows us to measure robust chemical abundances for galaxies with different star formation histories. Dwarf galaxies are especially useful because simple chemical evolution models can be used to determine the sources of r-process material. We have measured the r-process element barium with Keck/DEIMOS medium-resolution spectroscopy. We present the largest sample of barium abundances (almost 250 stars) in dwarf galaxies ever assembled. We measure [Ba/Fe] as a function of [Fe/H] in this sample and compare with existing [alpha/Fe] measurements. We have found that a large contribution of barium needs to occur at more delayed timescales than core-collapse supernovae in order to explain our observed abundances, namely the significantly more positive trend of the r-process component of [Ba/Fe] vs. [Fe/H] seen for [Fe/H] <~ -1.6 when compared to the [Mg/Fe] vs. [Fe/H] trend. We conclude that neutron star mergers are the most likely source of r-process enrichment in dwarf galaxies at early times.


title: Bayesian sparse reconstruction: a brute-force approach to astronomical imaging and machine learning

abstract: We present a principled Bayesian framework for signal reconstruction, in which the signal is modelled by basis functions whose number (and form, if required) is determined by the data themselves. This approach is based on a Bayesian interpretation of conventional sparse reconstruction and regularisation techniques, in which sparsity is imposed through priors via Bayesian model selection. We demonstrate our method for noisy 1- and 2-dimensional signals, including astronomical images. Furthermore, by using a product-space approach, the number and type of basis functions can be treated as integer parameters and their posterior distributions sampled directly. We show that order-of-magnitude increases in computational efficiency are possible from this technique compared to calculating the Bayesian evidences separately, and that further computational gains are possible using it in combination with dynamic nested sampling. Our approach can also be readily applied to neural networks, where it allows the network architecture to be determined by the data in a principled Bayesian manner by treating the number of nodes and hidden layers as parameters.


title: Dark Energy Survey Year 1 Results: The effect of intracluster light on photometric redshifts for weak gravitational lensing

abstract: We study the effect of diffuse intracluster light on the critical surface mass density estimated from photometric redshifts of lensing source galaxies, and the resulting bias in a weak lensing measurement of galaxy cluster mass. Under conservative assumptions, we find the bias to be negligible for imaging surveys like the Dark Energy Survey (DES) with a recommended scale cut of >=200 kpc distance from cluster centers. For significantly deeper source catalogs from present and future surveys like the Large Synoptic Survey Telescope (LSST) program, more conservative scale and source magnitude cuts or a correction of the effect may be necessary to achieve per-cent level lensing measurement accuracy, especially at the massive end of the cluster population.


title: Fast Semistochastic Heat-Bath Configuration Interaction

abstract: This paper presents in detail our fast semistochastic heat-bath configuration interaction (SHCI) method for solving the many-body Schrodinger equation. We identify and eliminate computational bottlenecks in both the variational and perturbative steps of the SHCI algorithm. We also describe the parallelization and the key data structures in our implementation, such as the distributed hash table. The improved SHCI algorithm enables us to include in our variational wavefunction two orders of magnitude more determinants than has been reported previously with other selected configuration interaction methods. We use our algorithm to calculate an accurate benchmark energy for the chromium dimer with the X2C relativistic Hamiltonian in the cc-pVDZ-DK basis, correlating 28 electrons in 76 spatial orbitals. Our largest calculation uses two billion Slater determinants in the variational space, and semistochastically includes perturbative contributions from at least trillions of additional determinants with better than 10 microhartree statistical uncertainty.


title: Quantum effective action for degenerate vector field theories

abstract: We calculate the divergent part of the one-loop effective action in curved spacetime for a particular class of second-order vector field operators with a degenerate principal part. The principal symbol of these operators has the structure of a longitudinal projector. In this case, standard heat-kernel techniques are not directly applicable. We present a method which reduces the problem to a nondegenerate scalar operator for which standard heat-kernel techniques are available. Interestingly, this method leads to the identification of an effective metric structure in the longitudinal sector. The one-loop divergences are compactly expressed in terms of invariants constructed from this metric.


title: Computing $N$-subjettiness for boosted jets

abstract: Jet substructure tools have proven useful in a number of high-energy particle-physics studies. A particular case is the discrimination, or tagging, between a boosted jet originated from an electroweak boson (signal), and a standard QCD parton (background). A common way to achieve this is to cut on a measure of the radiation inside the jet, i.e. a jet shape. Over the last few years, analytic calculations of jet substructure have allowed for a deeper understanding of these tools and for the development of more efficient ones. However, analytic calculations are often limited to the region where the jet shape is small. In this paper we introduce a new approach in perturbative QCD to compute jet shapes for a generic boosted jets, waiving the above limitation. We focus on an example common in the substructure literature: the jet mass distribution after a cut on the $N$-subjettiness $\tau_{21}$ ratio, extending previous works to the region relevant for phenomenology. We compare our analytic predictions to Monte Carlo simulations for both plain and SoftDrop-groomed jets. We use our results to construct analytically a decorrelated tagger.