title: A general theory for preferential sampling in environmental networks

abstract: This paper presents a general model framework for detecting the preferential sampling of environmental monitors recording an environmental process across space and/or time. This is achieved by considering the joint distribution of an environmental process with a site--selection process that considers where and when sites are placed to measure the process. The environmental process may be spatial, temporal or spatio--temporal in nature. By sharing random effects between the two processes, the joint model is able to establish whether site placement was stochastically dependent of the environmental process under study. The embedding into a spatio--temporal framework also allows for the modelling of the dynamic site---selection process itself. Real--world factors affecting both the size and location of the network can be easily modelled and quantified. Depending upon the choice of population of locations to consider for selection across space and time under the site--selection process, different insights about the precise nature of preferential sampling can be obtained. The general framework developed in the paper is designed to be easily and quickly fit using the R-INLA package. We apply this framework to a case study involving particulate air pollution over the UK where a major reduction in the size of a monitoring network through time occurred. It is demonstrated that a significant response--biased reduction in the air quality monitoring network occurred. We also show that the network was consistently unrepresentative of the levels of particulate matter seen across much of GB throughout the operating life of the network. Finally we show that this may have led to a severe over-reporting of the population--average exposure levels experienced across GB. This could have great impacts on estimates of the health effects of black smoke levels.


title: XNLI: Evaluating Cross-lingual Sentence Representations

abstract: State-of-the-art natural language processing systems rely on supervision in the form of annotated data to learn competent models. These models are generally trained on data in a single language (usually English), and cannot be directly used beyond that language. Since collecting data in every language is not realistic, there has been a growing interest in cross-lingual language understanding (XLU) and low-resource cross-language transfer. In this work, we construct an evaluation set for XLU by extending the development and test sets of the Multi-Genre Natural Language Inference Corpus (MultiNLI) to 15 languages, including low-resource languages such as Swahili and Urdu. We hope that our dataset, dubbed XNLI, will catalyze research in cross-lingual sentence understanding by providing an informative standard evaluation task. In addition, we provide several baselines for multilingual sentence understanding, including two based on machine translation systems, and two that use parallel data to train aligned multilingual bag-of-words and LSTM encoders. We find that XNLI represents a practical and challenging evaluation suite, and that directly translating the test data yields the best performance among available baselines.


title: IncSQL: Training Incremental Text-to-SQL Parsers with Non-Deterministic Oracles

abstract: We present a sequence-to-action parsing approach for the natural language to SQL task that incrementally fills the slots of a SQL query with feasible actions from a pre-defined inventory. To account for the fact that typically there are multiple correct SQL queries with the same or very similar semantics, we draw inspiration from syntactic parsing techniques and propose to train our sequence-to-action models with non-deterministic oracles. We evaluate our models on the WikiSQL dataset and achieve an execution accuracy of 83.7% on the test set, a 2.1% absolute improvement over the models trained with traditional static oracles assuming a single correct target SQL query. When further combined with the execution-guided decoding strategy, our model sets a new state-of-the-art performance at an execution accuracy of 87.1%.


title: Single-layer antiferromagnetic semiconductor CoS2 with the pentagonal structure

abstract: Structure-property relationships have always been guiding principles in discovering new materials. Here we explore the relationships to discover novel two-dimensional (2D) materials with the goal of identifying 2D magnetic semiconductors for spintronics applications. In particular, we report a density functional theory + $U$ study of single-layer antiferromagnetic (AFM) semiconductor CoS$_2$ with the pentagonal structure forming the so-called Cairo Tessellation. We find that this single-layer magnet exhibits an indirect bandgap of 1.06 eV with light electron and hole effective masses of 0.03 and 0.10 $m_0$, respectively, which may lead to high carrier mobilities. The hybrid density functional theory calculations correct the bandgap to 2.24 eV. We also compute the magnetocrystalline anisotropy energy (MAE), showing that the easy axis of the AFM ordering is out of plane with a sizable MAE of 153 $\mu$eV per Co ion. We further calculate the magnon frequencies at different spin-spiral vectors, based on which we estimate the N$\acute{e}$el temperatures to be 20.4 and 13.3 K using the mean field and random phase approximations, respectively. We then apply biaxial strains to tune the bandgap of single-layer pentagonal CoS$_2$. We find that the energy difference between the ferromagnetic and AFM structures strongly depends on the biaxial strain, but the ground state remains the AFM ordering. Although the low critical temperature prohibits the magnetic applications of single-layer pentagonal CoS$_2$ at room temperature, the excellent electrical properties may find this novel single-layer semiconductor applications in optoelectronic nanodevices.


title: Normal Modes in Magnetized Two-Fluid Spin Quantum Plasmas

abstract: We extend the classical two-fluid magnetohydrodynamic (MHD) formalism to include quantum effects such as electron Fermi pressure, Bohm pressure and spin couplings. At scales smaller than the electron skin-depth, the Hall effect and electron inertia must be taken into account, and can overlap with the quantum effects. We write down the full set of two-fluid quantum MHD (QMHD) and analyze the relative importance ofthese effects in the high density environments of neutron star atmospheres and white dwarf interiors, finding that for a broad range of parameters all these effects are operative. Of all spin interactions we analyze only the spin-magnetic coupling, as it is linear in $\hbar$ and consequently it is the strongest spin effect. We re-obtain the classical two-fluid MHD dispersion relations corresponding to the magnetosonic and Alfv\'en modes, modified by quantum effects. In the zero-spin case, for propagation parallel to the magnetic field, we find that the frequency of the fast mode is due to quantum effects modified by electron inertia, while the frequency of the Alfv\'en-slow sector has no quantum corrections. For perpendicular propagation, the fast-mode frequency is the same as for the parallel propagation plus a correction due only to classical two-fluid effects. When spin is considered, a whistler mode appears, which is due to two-fluid spin-magnetic interaction. There are no modifications due to spin for parallel propagation of magnetosonic and Alfv\'en waves, while for perpendicular prop agation a dispersive term due to spin arises in the two-fluid expression for the fast magnetosonic mode.


title: Gate-efficient simulation of molecular eigenstates on a quantum computer

abstract: A key requirement to perform simulations of large quantum systems on near-term quantum hardware is the design of quantum algorithms with short circuit depth that finish within the available coherence time. A way to stay within the limits of coherence is to reduce the number of gates by implementing a gate set that matches the requirements of the specific algorithm of interest directly in hardware. Here, we show that exchange-type gates are a promising choice for simulating molecular eigenstates on near-term quantum devices since these gates preserve the number of excitations in the system. Complementing the theoretical work by Barkoutsos et al. [PRA 98, 022322 (2018)], we report on the experimental implementation of a variational algorithm on a superconducting qubit platform to compute the eigenstate energies of molecular hydrogen. We utilize a parametrically driven tunable coupler to realize exchange-type gates that are configurable in amplitude and phase on two fixed-frequency superconducting qubits. With gate fidelities around 95% we are able to compute the eigenstates within an accuracy of 50 mHartree on average, a limit set by the coherence time of the tunable coupler.


title: Tire Noise Optimization Problem: a Mixed Integer Linear Program Approach

abstract: We present a Mixed Integer Linear Program (MILP) approach in order to model the nonlinear problem of minimizing the tire noise. We first take more industrial constraints into account than in a former work of the authors. Then, we associate a Branch-and-Cut algorithm to the MILP to obtain exact solutions. We compare our experimental results with those obtained by other methods.


title: Double Diffusion Encoding Prevents Degeneracy in Parameter Estimation of Biophysical Models in Diffusion MRI

abstract: Purpose: Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, the general Standard Model has recently shown that model parameter estimation from dMRI data is ill-posed unless very strong magnetic gradients are used. We analyse this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from Single Diffusion Encoding (SDE) to Double Diffusion Encoding (DDE) solves the ill-posedness and increases the accuracy of the parameter estimation. Methods: We analyse theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions. Results: We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space. Conclusions: DDE adds additional information for estimating the model parameters, unexplored by SDE, which is enough to solve the degeneracy in the NODDIDA model parameter estimation.


title: The NLS approximation for two dimensional deep gravity waves

abstract: This article is concerned with infinite depth gravity water waves in two space dimensions. We consider this system expressed in position-velocity potential holomorphic coordinates. Our goal is to study this problem with small wave packet data, and to show that this is well approximated by the cubic nonlinear Schr\"odinger equation (NLS) on the natural cubic time scale.


title: Kinetic Solvers with Adaptive Mesh in Phase Space for Low-Temperature Plasmas

abstract: We describe the implementation of 1d1v and 1d2v Vlasov and Fokker-Planck kinetic solvers with adaptive mesh refinement in phase space (AMPS) and coupling these kinetic solvers to Poisson equation solver for electric field. We demonstrate that coupling AMPS kinetic and electrostatic solvers can be done efficiently without splitting phase-space transport. We show that Eulerian fluid and kinetic solvers with dynamically adaptive Cartesian mesh can be used for simulations of collisionless plasma expansion into vacuum. The Vlasov-Fokker-Planck solver is demonstrated for the analysis of electron acceleration and scattering as well as the generation of runaway electrons in spatially inhomogeneous electric fields.


title: Non-Fickian current enabling the uphill diffusion of impurities diffusing by the mechanism of mobile impurity-defect pairs

abstract: Equations governing the diffusion mediated by the mobile pairs have been derived. In addition to known non-Fickian (n-F) behavior observable at small spatiotemporal scale under homogeneous defect distribution the equations predict another type of n-F diffusion that may take place only in the inhomogeneous case and can be operative at any scale. In the limit of slowly varying inhomogeneity the n-F diffusion current has been identified and shown to be able to induce the uphill diffusion. It has been argued that the latter should be observable in about one third of impurities in FCC metals.


title: Monitoring radiation damage in the LHCb Tracker Turicensis

abstract: This paper presents the techniques used to monitor radiation damage in the LHCb Tracker Turicensis during the LHC Runs 1 and 2. Bulk leakage currents in the silicon sensors were monitored continuously, while the full depletion voltages of the sensors were estimated at regular intervals by performing dedicated scans of the charge collection efficiency as a function of the applied bias voltage. Predictions of the expected leakage currents and full depletion voltages are extracted from the simulated radiation profile, the luminosity delivered by the LHC, and the thermal history of the silicon sensors. Good agreement between measurements and predictions is found.


title: PSR J2234+0611: A new laboratory for stellar evolution

abstract: We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.


title: Lyapunov analysis of multiscale dynamics: The slow bundle of the two-scale Lorenz 96 model

abstract: We investigate the geometrical structure of instabilities in the two-scales Lorenz 96 model through the prism of Lyapunov analysis. Our detailed study of the full spectrum of covariant Lyapunov vectors reveals the presence of a slow bundle in tangent space, composed by a set of vectors with a significant projection on the slow degrees of freedom; they correspond to the smallest (in absolute value) Lyapunov exponents and thereby to the longer time scales. We show that the dimension of the slow bundle is extensive in the number of both slow and fast degrees of freedom, and discuss its relationship with the results of a finite-size analysis of instabilities, supporting the conjecture that the slow-variable behavior is effectively determined by a non-trivial subset of degrees of freedom. More precisely, we show that the slow bundle corresponds to the Lyapunov spectrum region where fast and slow instability rates overlap, "mixing" their evolution into a set of vectors which simultaneously carry information on both scales. We suggest these results may pave the way for future applications to ensemble forecasting and data assimilations in weather and climate models.


title: Simulated Tempering Method in the Infinite Switch Limit with Adaptive Weight Learning

abstract: We investigate the theoretical foundations of the simulated tempering method and use our findings to design efficient algorithms. Employing a large deviation argument first used for replica exchange molecular dynamics [Plattner et al., J. Chem. Phys. 135:134111 (2011)], we demonstrate that the most efficient approach to simulated tempering is to vary the temperature infinitely rapidly. In this limit, we can replace the equations of motion for the temperature and physical variables by averaged equations for the latter alone, with the forces rescaled according to a position-dependent function defined in terms of temperature weights. The averaged equations are similar to those used in Gao's integrated-over-temperature method, except that we show that it is better to use a continuous rather than a discrete set of temperatures. We give a theoretical argument for the choice of the temperature weights as the reciprocal partition function, thereby relating simulated tempering to Wang-Landau sampling. Finally, we describe a self-consistent algorithm for simultaneously sampling the canonical ensemble and learning the weights during simulation. This algorithm is tested on a system of harmonic oscillators as well as a continuous variant of the Curie-Weiss model, where it is shown to perform well and to accurately capture the second-order phase transition observed in this model.


title: Seeing Tree Structure from Vibration

abstract: Humans recognize object structure from both their appearance and motion; often, motion helps to resolve ambiguities in object structure that arise when we observe object appearance only. There are particular scenarios, however, where neither appearance nor spatial-temporal motion signals are informative: occluding twigs may look connected and have almost identical movements, though they belong to different, possibly disconnected branches. We propose to tackle this problem through spectrum analysis of motion signals, because vibrations of disconnected branches, though visually similar, often have distinctive natural frequencies. We propose a novel formulation of tree structure based on a physics-based link model, and validate its effectiveness by theoretical analysis, numerical simulation, and empirical experiments. With this formulation, we use nonparametric Bayesian inference to reconstruct tree structure from both spectral vibration signals and appearance cues. Our model performs well in recognizing hierarchical tree structure from real-world videos of trees and vessels.


title: Learning Shape Priors for Single-View 3D Completion and Reconstruction

abstract: The problem of single-view 3D shape completion or reconstruction is challenging, because among the many possible shapes that explain an observation, most are implausible and do not correspond to natural objects. Recent research in the field has tackled this problem by exploiting the expressiveness of deep convolutional networks. In fact, there is another level of ambiguity that is often overlooked: among plausible shapes, there are still multiple shapes that fit the 2D image equally well; i.e., the ground truth shape is non-deterministic given a single-view input. Existing fully supervised approaches fail to address this issue, and often produce blurry mean shapes with smooth surfaces but no fine details. In this paper, we propose ShapeHD, pushing the limit of single-view shape completion and reconstruction by integrating deep generative models with adversarially learned shape priors. The learned priors serve as a regularizer, penalizing the model only if its output is unrealistic, not if it deviates from the ground truth. Our design thus overcomes both levels of ambiguity aforementioned. Experiments demonstrate that ShapeHD outperforms state of the art by a large margin in both shape completion and shape reconstruction on multiple real datasets.


title: Cwikel's bound reloaded

abstract: There are a couple of proofs by now for the famous Cwikel--Lieb--Rozenblum (CLR) bound, which is a semiclassical bound on the number of bound states for a Schr\"odinger operator, proven in the 1970s. Of the rather distinct proofs by Cwikel, Lieb, and Rozenblum, the one by Lieb gives the best constant, the one by Rozenblum does not seem to yield any reasonable estimate for the constants, and Cwikel's proof is said to give a constant which is at least about 2 orders of magnitude off the truth. This situation did not change much during the last 40+ years. It turns out that this common belief, i.e, Cwikel's approach yields bad constants, is not set in stone: We give a drastic simplification of Cwikel's original approach which leads to an astonishingly good bound for the constant in the CLR inequality. Our proof is also quite flexible and leads to rather precise bounds for a large class of Schr\"odinger-type operators with generalized kinetic energies. Moreover, it highlights a natural but overlooked connection of the CLR bound with bounds for maximal Fourier multipliers from harmonic analysis.


title: Physical Primitive Decomposition

abstract: Objects are made of parts, each with distinct geometry, physics, functionality, and affordances. Developing such a distributed, physical, interpretable representation of objects will facilitate intelligent agents to better explore and interact with the world. In this paper, we study physical primitive decomposition---understanding an object through its components, each with physical and geometric attributes. As annotated data for object parts and physics are rare, we propose a novel formulation that learns physical primitives by explaining both an object's appearance and its behaviors in physical events. Our model performs well on block towers and tools in both synthetic and real scenarios; we also demonstrate that visual and physical observations often provide complementary signals. We further present ablation and behavioral studies to better understand our model and contrast it with human performance.


title: Leading order CFT analysis of multi-scalar theories in d>2

abstract: We investigate multi-field multicritical scalar theories using CFT constraints on two- and three-point functions combined with the Schwinger-Dyson equation. This is done in general and without assuming any symmetry for the models, which we just define to admit a Landau-Ginzburg description that includes the most general critical interactions built from monomials of the form $\phi_{i_1} \cdots \phi_{i_m}$. For all such models we analyze to the leading order of the $\epsilon$-expansion the anomalous dimensions of the fields and those of the composite quadratic operators. For models with even $m$ we extend the analysis to an infinite tower of composite operators of arbitrary order. The results are supplemented by the computation of some families of structure constants. We also find the equations which constrain the nontrivial critical theories at leading order and show that they coincide with the ones obtained with functional perturbative RG methods. This is done for the case $m=3$ as well as for all the even models. We ultimately specialize to $S_q$ symmetric models, which are related to the $q$-state Potts universality class, and focus on three realizations appearing below the upper critical dimensions $6$, $4$ and $\frac{10}{3}$, which can thus be nontrivial CFTs in three dimensions.


title: A controlled-NOT gate for frequency-bin qubits

abstract: The realization of strong photon-photon interactions has presented an enduring challenge across photonics, particularly in quantum computing, where two-photon gates form essential components for scalable quantum information processing (QIP). While linear-optic schemes have enabled probabilistic entangling gates in spatio-polarization encoding, solutions for many other useful degrees of freedom remain missing. In particular, no two-photon gate for the important platform of frequency encoding has been experimentally demonstrated, due in large part to the additional challenges imparted by the mismatched wavelengths of the interacting photons. In this article, we design and implement the first entangling gate for frequency-bin qubits, a coincidence-basis controlled-NOT (CNOT), using line-by-line pulse shaping and electro-optic modulation. We extract a quantum gate fidelity of $0.91 \pm 0.01$ via a novel parameter inference approach based on Bayesian machine learning, which enables accurate gate reconstruction from measurements in the two-photon computational basis alone. Our CNOT imparts a single-photon frequency shift controlled by the frequency of another photon---an important capability in itself---and should enable new directions in fiber-compatible QIP.


title: Channel Polarization through the Lens of Blackwell Measures

abstract: Each memoryless binary-input channel (BIC) can be uniquely described by its Blackwell measure, which is a probability distribution on the unit interval $[0,1]$ with mean $1/2$. Conversely, any such probability distribution defines a BIC. The evolution of the Blackwell measure under Arikan's polar transform is derived for general BICs, and is analogous to density evolution as cited in the literature. The present analysis emphasizes functional equations. Consequently, the evolution of a variety of channel functionals is characterized, including the symmetric capacity, Bhattacharyya parameter, moments of information density, Hellinger affinity, Gallager's reliability function, the Hirschfeld-Gebelein-Renyi maximal correlation, and the Bayesian information gain. The evolution of measure is specialized for symmetric BICs according to their decomposition into binary symmetric (sub)-channels (BSCs), which simplifies iterative computations and the construction of polar codes. It is verified that, as a consequence of the Blackwell--Sherman--Stein theorem, all channel functionals $\mathrm{I}_f$ that can be expressed as an expectation of a convex function $f$ with respect to the Blackwell measure of a channel polarize in each iteration due to the polar transformation on the class of symmetric BICs. Moreover, for $f$ either convex or non-convex, a necessary and sufficient condition is established to determine whether the random process associated with each $\mathrm{I}_f$ is a martingale, submartingale, or supermartingale. Represented via functional inequalities in terms of $f$, this condition is numerically verifiable for all $\mathrm{I}_f$, and can generate analytical proofs. To exhibit one such proof, it is shown that the random process associated with the squared maximal correlation parameter is a supermartingale, and converges almost surely on the unit interval $[0,1]$.


title: Derivative-free online learning of inverse dynamics models

abstract: This paper discusses online algorithms for inverse dynamics modelling in robotics. Several model classes including rigid body dynamics (RBD) models, data-driven models and semiparametric models (which are a combination of the previous two classes) are placed in a common framework. While model classes used in the literature typically exploit joint velocities and accelerations, which need to be approximated resorting to numerical differentiation schemes, in this paper a new `derivative-free' framework is proposed that does not require this preprocessing step. An extensive experimental study with real data from the right arm of the iCub robot is presented, comparing different model classes and estimation procedures, showing that the proposed `derivative-free' methods outperform existing methodologies.


title: Where Does Haydn End and Mozart Begin? Composer Classification of String Quartets

abstract: For centuries, the history and music of Joseph Franz Haydn and Wolfgang Amadeus Mozart have been compared by scholars. Recently, the growing field of music information retrieval (MIR) has offered quantitative analyses to complement traditional qualitative analyses of these composers. In this MIR study, we classify the composer of Haydn and Mozart string quartets based on the content of their scores. Our contribution is an interpretable statistical and machine learning approach that provides high classification accuracies and musical relevance. We develop novel global features that are automatically computed from symbolic data and informed by musicological Haydn-Mozart comparative studies, particularly relating to the sonata form. Several of these proposed features are found to be important for distinguishing between Haydn and Mozart string quartets. Our Bayesian logistic regression model attains leave-one-out classification accuracies over 84%, higher than prior works and providing interpretations that could aid in assessing musicological claims. Overall, our work can help expand the longstanding dialogue surrounding Haydn and Mozart and exemplify the benefit of interpretable machine learning in MIR, with potential applications to music generation and classification of other classical composers.


title: Computer Vision-aided Atom Tracking in STEM Imaging

abstract: To address the SMC'17 data challenge -- "Data mining atomically resolved images for material properties", we first used the classic "blob detection" algorithms developed in computer vision to identify all atom centers in each STEM image frame. With the help of nearest neighbor analysis, we then found and labeled every atom center common to all the STEM frames and tracked their movements through the given time interval for both Molybdenum or Selenium atoms.


title: MSc Dissertation: Exclusive Row Biclustering for Gene Expression Using a Combinatorial Auction Approach

abstract: The availability of large microarray data has led to a growing interest in biclustering methods in the past decade. Several algorithms have been proposed to identify subsets of genes and conditions according to different similarity measures and under varying constraints. In this paper we focus on the exclusive row biclustering problem for gene expression data sets, in which each row can only be a member of a single bicluster while columns can participate in multiple ones. This type of biclustering may be adequate, for example, for clustering groups of cancer patients where each patient (row) is expected to be carrying only a single type of cancer, while each cancer type is associated with multiple (and possibly overlapping) genes (columns). We present a novel method to identify these exclusive row biclusters through a combination of existing biclustering algorithms and combinatorial auction techniques. We devise an approach for tuning the threshold for our algorithm based on comparison to a null model in the spirit of the Gap statistic approach. We demonstrate our approach on both synthetic and real-world gene expression data and show its power in identifying large span non-overlapping rows sub matrices, while considering their unique nature. The Gap statistic approach succeeds in identifying appropriate thresholds in all our examples.


title: The Free Fermion Anomaly and Representations of the Pin Groups

abstract: We consider a duality between a boson on a ring and a free fermion and show that they have an anomaly which corresponds to the states transforming under double covers of O(2). There are two (in general not isometric) double covers of O(2), known as Pin$_+(2)$ and Pin$_-(2)$. These can in general be distinguished at the group level by checking whether reflections square to $\pm 1$. We show that in irreducible representations in complex Hilbert spaces the commutation of time reversal and charge conjugation gives another method for distinguishing Pin$_+(2)$ and Pin$_-(2)$. While we only demonstrate the duality for a single fermion, the anomaly is also present in any number of free fermions. For an even number of fermions we show that the two double covers of O(2n) are isomorphic. For an odd number of fermions we show that the distinction between irreducible representations of Pin$_+$ and Pin$_-$ can still be detected by the commutation of time reversal and charge conjugation namely Pin$_\pm(2n)$ will have $TCT = \pm C$.


title: Competition and partnership between conformity and payoff-based imitations in social dilemmas

abstract: Learning from a partner who collects higher payoff is a frequently used working hypothesis in evolutionary game theory. One of the alternative dynamical rules is when the focal player prefers to follow the strategy choice of the majority in the local neighborhood, which is often called as conformity-driven strategy update. In this work we assume that both strategy learning methods are present and compete for space within the framework of a coevolutionary model. Our results reveal that the presence of payoff-driven strategy learning method becomes exclusive for high succer's payoff and/or high temptation values that represent a snowdrift game dilemma situation. In general, however, the competition of the mentioned strategy learning methods could be useful to enlarge the parameter space where only cooperators prevail. This success of cooperation is based on the enforced coordination of cooperator players which reveals the benefit of the latter strategy. Interestingly, the payoff-based and the conformity-based cooperator players can form an effective alliance against defectors that can also extend the parameter space of full cooperator solution in the stag-hunt game region. Our work highlights that the coevolution of strategies and individual features such as learning method can provide novel type of pattern formation mechanism that cannot be observed in a static model, hence remains hidden in traditional models.

categories: physics.soc-ph cond-mat.stat-mech cs.GT q-bio.PE


title: Large-scale motions from a direct numerical simulation of a turbulent boundary layer

abstract: A study of large-scale motions from a new direct numerical simulation database of the turbulent boundary layer up to Re{\theta} ~ 2500 is presented. The statistics of large-scale streamwise structures have been investigated using two-dimensional and three-dimensional extraction procedures. The large-scale structures are abstracted using a robust skeletonization method usually applied to other research domains to simplify complex 3D objects. Different structure parameters such as the length, shape or angle are investigated. The features of the detected structures are compared to their mean counterparts extracted from two-point correlations. Structures as large as 10 boundary layer thickness are observed. The streamwise length of these structures follows a -2 power law distribution, similar to the experimental findings at higher Reynolds numbers.


title: Demonstration of displacement sensing of a mg-scale pendulum for mm- and mg- scale gravity measurements

abstract: Gravity generated by large masses has been observed using a variety of probes from atomic interferometers to torsional balances. However, gravitational coupling between small masses has never been observed so far. Here, we demonstrate sensitive displacement sensing of the Brownian motion of an optically trapped 7-mg pendulum motion whose natural quality factor is increased to $10^8$ through dissipation dilution. The sensitivity for an integration time of one second corresponds to the displacement generated by the gravitational coupling between the probe and a mm separated 100 mg mass, whose position is modulated at the pendulum mechanical resonant frequency. Development of such a sensitive displacement sensor using a mg-scale device will pave the way for a new class of experiments where gravitational coupling between small masses in quantum regimes can be achieved.


title: Submodular Secretary Problem with Shortlists

abstract: In submodular $k$-secretary problem, the goal is to select $k$ items in a randomly ordered input so as to maximize the expected value of a given monotone submodular function on the set of selected items. In this paper, we introduce a relaxation of this problem, which we refer to as submodular $k$-secretary problem with shortlists. In the proposed problem setting, the algorithm is allowed to choose more than $k$ items as part of a shortlist. Then, after seeing the entire input, the algorithm can choose a subset of size $k$ from the bigger set of items in the shortlist. We are interested in understanding to what extent this relaxation can improve the achievable competitive ratio for the submodular $k$-secretary problem. In particular, using an $O(k)$ shortlist, can an online algorithm achieve a competitive ratio close to the best achievable online approximation factor for this problem? We answer this question affirmatively by giving a polynomial time algorithm that achieves a $1-1/e-\epsilon -O(k^{-1})$ competitive ratio for any constant $\epsilon > 0$, using a shortlist of size $\eta_\epsilon(k) = O(k)$. Also, for the special case of m-submodular functions, we demonstrate an algorithm that achieves a $1-\epsilon$ competitive ratio for any constant $\epsilon > 0$, using an $O(1)$ shortlist. Finally, we show that our algorithm can be implemented in the streaming setting using a memory buffer of size $\eta_\epsilon(k) = O(k)$ to achieve a $1 - 1/e - \epsilon-O(k^{-1})$ approximation for submodular function maximization in the random order streaming model. This substantially improves upon the previously best known approximation factor of $1/2 + 8 \times 10^{-14}$ [Norouzi-Fard et al. 2018] that used a memory buffer of size $O(k \log k)$.


title: Quotients of the magmatic operad: lattice structures and convergent rewrite systems

abstract: We study quotients of the magmatic operad, that is the free nonsymmetric operad over one binary generator. In the linear setting, we show that the set of these quotients admits a lattice structure and we show an analog of the Grassmann formula for the dimensions of these operads. In the nonlinear setting, we define comb associative operads, that are operads indexed by nonnegative integers generalizing the associative operad. We show that the set of comb associative operads admits a lattice structure, isomorphic to the lattice of nonnegative integers equipped with the division order. Driven by computer experimentations, we provide a finite convergent presentation for the comb associative operad in correspondence with~$3$. Finally, we study quotients of the magmatic operad by one cubic relation by expressing their Hilbert series and providing combinatorial realizations.


title: Convolution quadrature methods for time-domain scattering from unbounded penetrable interfaces

abstract: This paper presents a class of boundary integral equation methods for the numerical solution of acoustic and electromagnetic time-domain scattering problems in the presence of unbounded penetrable interfaces in two-spatial dimensions. The proposed methodology relies on Convolution Quadrature (CQ) methods in conjunction with the recently introduced Windowed Green Function (WGF) method. As in standard time-domain scattering from bounded obstacles, a CQ method of the user's choice is utilized to transform the problem into a finite number of (complex) frequency-domain problems posed on the domains involving penetrable unbounded interfaces. Each one of the frequency-domain transmission problems is then formulated as a second-kind integral equation that is effectively reduced to a bounded interface by means of the WGF method---which introduces errors that decrease super-algebraically fast as the window size increases. The resulting windowed integral equations can then be solved by means of any (accelerated or unaccelerated) off-the-shelf Helmholtz boundary integral equation solver capable of handling complex wavenumbers with a large imaginary part. A high-order Nystr\"om method based on Alpert quadrature rules is utilized here. A variety of numerical examples including wave propagation in open waveguides as well as scattering from multiply layered media demonstrate the capabilities of the proposed approach.


title: Work and Heat Value of Bound Entanglement

abstract: Entanglement has recently been recognized as an energy resource which can outperform classical resources if decoherence is relatively low. Multi-atom entangled states can mutate irreversibly to so called bound entangled (BE) states under noise. Resource value of BE states in information applications has been under critical study and a few cases where they can be useful have been identified. We explore the energetic value of typical BE states. Maximal work extraction is determined in terms of ergotropy. Since the BE states are non-thermal, extracting heat from them is less obvious. We compare single and repeated interaction schemes to operationally define and harvest heat from BE states. BE and free entangled (FE) states are compared in terms of their ergotropy and maximal heat values. Distinct roles of distillability in work and heat values of FE and BE states are pointed out. Decoherence effects in dynamics of ergotropy and mutation of FE states into BE states are examined to clarify significance of the work value of BE states. Thermometry of distillability of entanglement using micromaser cavity is proposed.


title: A mean-field limit of the Lohe matrix model and emergent dynamics

abstract: The Lohe matrix model is a continuous-time dynamical system describing the collective dynamics of group elements in the unitary group manifold, and it has been introduced as a toy model of a non abelian generalization of the Kuramoto phase model. In the absence of couplings, it reduces to the finite-dimensional decoupled free Schr\"{o}dinger equations with constant Hamiltonians. In this paper, we study a rigorous mean-field limit of the Lohe matrix model which results in a Vlasov type equation for the probability density function on the corresponding phase space. We also provide two different settings for the emergent synchronous dynamics of the kinetic Lohe equation in terms of the initial data and the coupling strength.


title: Whole-system ultrasound resonances as the basis for acoustophoresis in all-polymer microfluidic devices

abstract: Using a previously well-tested numerical model, we demonstrate theoretically that good acoustophoresis can be obtained in a microchannel embedded in an acoustically soft, all-polymer chip, by excitation of whole-system ultrasound resonances. In contrast to conventional techniques based on a standing bulk acoustic wave inside a liquid-filled microchannel embedded in an elastic, acoustically hard material, such as glass or silicon, the proposed whole-system resonance does not need a high acoustic contrast between the liquid and surrounding solid. Instead, it relies on the very high acoustic contrast between the solid and the surrounding air. In microchannels of usual dimensions, we demonstrate the existence of whole-system resonances in an all-polymer device, which support acoustophoresis of a quality fully comparable to that of a conventional hard-walled system. Our results open up for using cheap and easily processable polymers in a controlled manner to design and fabricate microfluidic devices for single-use acoustophoresis.


title: High Throughput Cryptocurrency Routing in Payment Channel Networks

abstract: Despite growing adoption of cryptocurrencies, making fast payments at scale remains a challenge. Payment channel networks (PCNs) such as the Lightning Network have emerged as a viable scaling solution. However, completing payments on PCNs is challenging: payments must be routed on paths with sufficient funds. As payments flow over a single channel (link) in the same direction, the channel eventually becomes depleted and cannot support further payments in that direction; hence, naive routing schemes like shortest-path routing can deplete key payment channels and paralyze the system. Today's PCNs also route payments atomically, worsening the problem. In this paper, we present Spider, a routing solution that "packetizes" transactions and uses a multi-path transport protocol to achieve high-throughput routing in PCNs. Packetization allows Spider to complete even large transactions on low-capacity payment channels over time, while the multi-path congestion control protocol ensures balanced utilization of channels and fairness across flows. Extensive simulations comparing Spider with state-of-the-art approaches shows that Spider requires less than 25% of the funds to successfully route over 95% of transactions on balanced traffic demands, and offloads 4x more transactions onto the PCN on imbalanced demands.


title: Bounds on extra dimensions from micro black holes in the context of the metastable Higgs vacuum

abstract: We estimate the rate at which collisions between ultra-high energy cosmic rays can form small black holes in models with extra dimensions. If recent conjectures about false vacuum decay catalyzed by black hole evaporation apply, the lack of vacuum decay events in our past light cone may place new bounds on the black hole formation rate and thus on the fundamental scale of gravity in these models. For theories with fundamental scale $E_{*}$ above the Higgs instability scale of the Standard Model, we find a lower bound on $E_{*}$ that is within about an order of magnitude of the energy where the cosmic ray spectrum begins to show suppression from the GZK effect. Otherwise, the abundant formation of semiclassical black holes with short lifetimes would likely initiate vacuum decay. Assuming a Higgs instability scale at the low end of the range compatible with experimental data, the excluded range is approximately $10^{17} \,\text{eV} \lesssim E_{*} \leq 10^{18.8}\,\text{eV}$ for theories with $n=1$ extra dimension, narrowing to $10^{17}\,\text{eV} \lesssim E_{*} \leq 10^{18.1}\,\text{eV}$ for $n=6$. These bounds rule out regions of parameter space that are inaccessible to collider experiments, small-scale gravity tests, or estimates of Kaluza-Klein processes in neutron stars and supernovae.


title: Stimulated Raman with Broadband LED Stokes Source for Analysis of Glucose

abstract: We demonstrate stimulated Raman gain using a broadband LED Stokes source to measure vibrational spectra of aqueous glucose solutions. This versatile and cost-effective method increases Raman signal for a variety of applications. We measured both stimulated Raman and spontaneous Raman spectra of glucose solutions with concentrations up to 10 mM with a photon counter and lock-in amplifier. We built partial least squares regression models based on both stimulated Raman and spontaneous Raman spectral data measured with each instrument for predicting concentrations of the glucose solutions. The stimulated Raman spectra measured with the lock-in amplifier based model had the strongest predictive power and predicted the concentrations of the test set of glucose solutions with a mean squared error value an order of magnitude lower than those of the spontaneous Raman based model.


title: Estimating Historical Functional Linear Models with a Nested Group Bridge Approach

abstract: We study a scalar-on-function historical linear regression model which assumes that the functional predictor does not influence the response when the time passes a certain cutoff point. We approach this problem from the perspective of locally sparse modeling, where a function is locally sparse if it is zero on a substantial portion of its defining domain. In the historical linear model, the slope function is exactly a locally sparse function that is zero beyond the cutoff time. A locally sparse estimate then gives rise to an estimate of the cutoff time. We propose a nested group bridge penalty that is able to specifically shrink the tail of a function. Combined with the B-spline basis expansion and penalized least squares, the nested group bridge approach can identify the cutoff time and produce a smooth estimate of the slope function simultaneously. The proposed locally sparse estimator is shown to be consistent, while its numerical performance is illustrated by simulation studies. The proposed method is demonstrated with an application of determining the effect of the past engine acceleration on the current particulate matter emission.


title: Polynomial mixing time of edge flips on quadrangulations

abstract: We establish the first polynomial upper bound for the mixing time of random edge flips on rooted quadrangulations: we show that the spectral gap of the edge flip Markov chain on quadrangulations with $n$ faces admits, up to constants, an upper bound of $n^{-5/4}$ and a lower bound of $n^{-11/2}$. In order to obtain the lower bound, we also consider a very natural Markov chain on plane trees (or, equivalently, on Dyck paths) and improve the previous lower bound for its spectral gap obtained by Shor and Movassagh.


title: Switching quantum reference frames in the N-body problem and the absence of global relational perspectives

abstract: Given the importance of quantum reference frames (QRFs) to both quantum and gravitational physics, it is pertinent to develop a systematic method for switching between the descriptions of physics relative to different choices of QRFs, which is valid in both fields. Here we continue with such a unifying approach, begun in arxiv:1809.00556, whose key ingredient is a symmetry principle, which enforces physics to be relational. Thanks to gauge related redundancies, this leads to a perspective-neutral structure which contains all frame choices at once and via which frame perspectives can be consistently switched. Formulated in the language of constrained systems, the perspective-neutral structure is the constraint surface classically and the gauge invariant Hilbert space in the Dirac quantized theory. By contrast, a perspective relative to a specific frame corresponds to a gauge choice and the associated reduced phase and Hilbert space. QRF changes thus amount to a gauge transformation. We show that they take the form of `quantum coordinate changes'. We illustrate this in a general mechanical model, namely the relational $N$-body problem in 3D space with rotational and translational symmetry. This model is especially interesting because it features the Gribov problem so that globally valid gauge fixing conditions, and hence relational frame perspectives, are absent. The constraint surface is topologically non-trivial and foliated by 3-, 5- and 6-dimensional gauge orbits, where the lower dimensional orbits are a set of measure zero. The $N$-body problem also does not admit globally valid canonically conjugate pairs of Dirac observables. These challenges notwithstanding, we exhibit how one can construct the QRF transformations for the 3-body problem. Our construction also sheds new light on the generic inequivalence of Dirac and reduced quantization through its interplay with QRF perspectives.


title: Cosmological Tests with the Joint Lightcurve Analysis

abstract: We examine whether a comparison between wCDM and R_h=ct using merged Type Ia SN catalogs produces results consistent with those based on a single homogeneous sample. Using the Betoule et al. (2014) joint lightcurve analysis (JLA) of a combined sample of 613 events from SNLS and SDSS-II, we estimate the parameters of the two models and compare them. We find that the improved statistics can alter the model selection in some cases, but not others. In addition, based on the model fits, we find that there appears to be a lingering systematic offset of ~0.04-0.08 mag between the SNLS and SDSS-II sources, in spite of the cross-calibration in the JLA. Treating wCDM, LCDM and R_h=ct as separate models, we find in an unbiased pairwise statistical comparison that the Bayes Information Criterion (BIC) favors the R_h=ct universe with a likelihood of 82.8% versus 17.2% for wCDM, but the ratio of likelihoods is reversed (16.2% versus 83.8%) when w_de=-1 (i.e., LCDM) and strongly reversed (1.0% versus 99.0%) if in addition k=0 (i.e., flat LCDM). We point out, however, that the value of k is a measure of the net energy (kinetic plus gravitational) in the Universe and is not constrained theoretically, though some models of inflation would drive k to 0 due to an expansion-enforced dilution. Since we here consider only the basic LCDM model, the value of k needs to be measured and, therefore, the pre-assumption of flatness introduces a significant bias into the BIC.


title: Effect of Blast Exposure on Gene-Gene Interactions

abstract: Repeated exposure to low-level blast may initiate a range of adverse health problem such as traumatic brain injury (TBI). Although many studies successfully identified genes associated with TBI, yet the cellular mechanisms underpinning TBI are not fully elucidated. In this study, we investigated underlying relationship among genes through constructing transcript Bayesian networks using RNA-seq data. The data for pre- and post-blast transcripts, which were collected on 33 individuals in Army training program, combined with our system approach provide unique opportunity to investigate the effect of blast-wave exposure on gene-gene interactions. Digging into the networks, we identified four subnetworks related to immune system and inflammatory process that are disrupted due to the exposure. Among genes with relatively high fold change in their transcript expression level, ATP6V1G1, B2M, BCL2A1, PELI, S100A8, TRIM58 and ZNF654 showed major impact on the dysregulation of the gene-gene interactions. This study reveals how repeated exposures to traumatic conditions increase the level of fold change of transcript expression and hypothesizes new targets for further experimental studies.


title: Negative Update Intervals in Deep Multi-Agent Reinforcement Learning

abstract: In Multi-Agent Reinforcement Learning (MA-RL), independent cooperative learners must overcome a number of pathologies to learn optimal joint policies. Addressing one pathology often leaves approaches vulnerable towards others. For instance, hysteretic Q-learning addresses miscoordination while leaving agents vulnerable towards misleading stochastic rewards. Other methods, such as leniency, have proven more robust when dealing with multiple pathologies simultaneously. However, leniency has predominately been studied within the context of strategic form games (bimatrix games) and fully observable Markov games consisting of a small number of probabilistic state transitions. This raises the question of whether these findings scale to more complex domains. For this purpose we implement a temporally extend version of the Climb Game, within which agents must overcome multiple pathologies simultaneously, including relative overgeneralisation, stochasticity, the alter-exploration and moving target problems, while learning from a large observation space. We find that existing lenient and hysteretic approaches fail to consistently learn near optimal joint-policies in this environment. To address these pathologies we introduce Negative Update Intervals-DDQN (NUI-DDQN), a Deep MA-RL algorithm which discards episodes yielding cumulative rewards outside the range of expanding intervals. NUI-DDQN consistently gravitates towards optimal joint-policies in our environment, overcoming the outlined pathologies.


title: The Supersonic Project: rotational effects of supersonic motions on the first structures in the Universe

abstract: We introduce the "Supersonic Project," aimed at investigating the effects of the supersonic relative velocity between dark matter (DM) and baryons at high redshift using a combination of analytical calculations and cosmological simulations. In this paper, we study the effect of this stream velocity on the angular momentum of the first structures in the early Universe using simulations. We focus on DM haloes and their gas component as well as the recently predicted supersonically-induced gas objects (SIGOs) that arise as a result of the stream velocity phase shift. We find that the spin parameter of the gas component in these first haloes is increased with the stream velocity. Moreover, we find that when the stream velocity is taken into account, the angular momentum vectors of the DM component and the gas component are typically misaligned and this misalignment angle has a nearly isotropic distribution. The spin parameter value of the gas component is higher than in the no stream velocity case, which even in the absence of cooling, may result in more prolate objects. We also generalize the spin parameter to the SIGOs and find that they typically have a larger spin parameter with respect to their dark matter counterparts and that there is no correlation of the spin parameter and the prolateness of such structures. We speculate that SIGOs may be observed as very low luminosity objects in the early Universe and may serve as potential progenitors of Little Blue Dot-like systems.


title: Colloquium: Ionic phenomena in nanoscale pores through 2D materials

abstract: Ion transport through nanopores permeates through many areas of science and technology, from cell behavior to sensing and separation to catalysis and batteries. Two-dimensional materials, such as graphene, molybdenum disulfide (MoS$_2$), and hexagonal boron nitride (hBN), are recent additions to these fields. Low-dimensional materials present new opportunities to develop filtration, sensing, and power technologies, encompassing ion exclusion membranes, DNA sequencing, single molecule detection, osmotic power generation, and beyond. Moreover, the physics of ionic transport through pores and constrictions within these materials is a distinct realm of competing many-particle interactions (e.g., solvation/dehydration, electrostatic blockade, hydrogen bond dynamics) and confinement. This opens up alternative routes to creating biomimetic pores and may even give analogues of quantum phenomena, such as quantized conductance, in the classical domain. These prospects make membranes of 2D materials -- i.e., 2D membranes -- fascinating. We will discuss the physics and applications of ionic transport through nanopores in 2D membranes.opores in 2D membranes.


title: EM counterparts of structured jets from 3D GRMHD simulations

abstract: GW170817/GRB170817A has offered unprecedented insight into binary neutron star post-merger systems. Its Prompt and afterglow emission imply the presence of a tightly collimated relativistic jet with a smooth transverse structure. However, it remains unclear whether and how the central engine can produce such structured jets. Here, we utilize 3D GRMHD simulations starting with a black hole surrounded by a magnetized torus with properties typically expected of a post-merger system. We follow the jet, as it is self-consistently launched, from the scale of the compact object out to more than 3 orders of magnitude in distance. We find that this naturally results in a structured jet, which is collimated by the disk wind into a half-opening angle of roughly 10 degrees, its emission can explain features of both the prompt and afterglow emission of GRB170817A for a 30 degree observing angle. Our work is the first to compute the afterglow, in the context of a binary merger, from a relativistic magnetized jet self-consistently generated by an accreting black hole, with the jet's transverse structure determined by the accretion physics and not prescribed at any point.


title: Absence of unidirectionally propagating surface plasmon-polaritons in nonreciprocal plasmonics

abstract: In the presence of an external magnetic field, the surface plasmon polariton that exists at the metal-dielectric interface is believed to support a unidirectional frequency range near the surface plasmon frequency, where the surface plasmon polariton propagates along one but not the opposite direction. Recent works have pointed to some of the paradoxical consequences of such a unidirectional range, including in particular the violation of the time-bandwidth product constraint that should otherwise apply in general in static systems. Here we show that such a unidirectional frequency range is nonphysical, using both a general thermodynamic argument, and a detailed calculation based on a nonlocal hydrodynamic Drude model for the metal permittivity. Our calculation reveals that the surface plasmon-polariton remains bidirectional for all frequencies. This work overturns a long-held belief in nonreciprocal photonics, and highlights the importance of quantum plasmonic concepts for the understanding of nonreciprocal plasmonic effects.


title: SymBuild: a package for the computation of integrable symbols in scattering amplitudes

abstract: The article presents and documents the Mathematica package SymBuild. This package implements the computation and manipulation of integrable symbols appearing in various calculations in high-energy scattering amplitudes. By using Gr\"obner bases, implementing various simplifications and by the potential utilization of the C++ program SpaSM, integrable symbols in a large class of alphabets (including roots) can be computed to high weights.


title: Resolving the Weinberg Paradox with Topology

abstract: Long ago Weinberg showed, from first principles, that the amplitude for a single photon exchange between an electric current and a magnetic current violates Lorentz invariance. The obvious conclusion at the time was that monopoles were not allowed in quantum field theory. Since the discovery of topological monopoles there has thus been a paradox. On the one hand, topological monopoles are constructed in Lorentz invariant quantum field theories, while on the other hand, the low-energy effective theory for such monopoles will reproduce Weinberg's result. We examine a toy model where both electric and magnetic charges are perturbatively coupled and show how soft-photon resummation for hard scattering exponentiates the Lorentz violating pieces to a phase that is the covariant form of the Aharonov-Bohm phase due to the Dirac string. The modulus of the scattering amplitudes (and hence observables) are Lorentz invariant, and when Dirac charge quantization is imposed the amplitude itself is also Lorentz invariant. For closed paths there is a topological component of the phase that relates to aspects of 4D topological quantum field theory.


title: The faint end of the Centaurus A satellite luminosity function

abstract: The Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS) is constructing a wide-field map of the resolved stellar populations in the extended halos of these two nearby, prominent galaxies. We present new Magellan/Megacam imaging of a $\sim3$ deg$^2$ area around Centaurus A (Cen A), which filled in much of our coverage to its south, leaving a nearly complete halo map out to a projected radius of $\sim$150 kpc and allowing us to identify two new resolved dwarf galaxies. We have additionally obtained deep Hubble Space Telescope (HST) optical imaging of eleven out of the thirteen candidate dwarf galaxies identified around Cen A and presented in Crnojevi\'c et al. (2016): seven are confirmed to be satellites of Cen A, while four are found to be background galaxies. We derive accurate distances, structural parameters, luminosities and photometric metallicities for the seven candidates confirmed by our HST/ACS imaging. We further study the stellar population along the $\sim$60 kpc long (in projection) stream associated with Dw3, which likely had an initial brightness of $M_{V}$$\sim$$-$15 and shows evidence for a metallicity gradient along its length. Using the total sample of eleven dwarf satellites discovered by the PISCeS survey, as well as thirteen brighter previously known satellites of Cen A, we present a revised galaxy luminosity function for the Cen A group down to a limiting magnitude of $M_V\sim-8$, which has a slope of $-1.14\pm0.17$, comparable to that seen in the Local Group and in other nearby groups of galaxies.


title: Stochastic cosmic ray sources and the TeV break in the all-electron spectrum

abstract: Despite significant progress over more than 100 years, no accelerator has been unambiguously identified as the source of the locally measured flux of cosmic rays. High-energy electrons and positrons are of particular importance in the search for nearby sources as radiative energy losses constrain their propagation to distances of about 1 kpc around 1 TeV. At the highest energies, the spectrum is therefore dominated and shaped by only a few sources whose properties can be inferred from the fine structure of the spectrum at energies currently accessed by experiments like AMS-02, CALET, DAMPE, Fermi-LAT, H.E.S.S. and ISS-CREAM. We present a stochastic model of the Galactic all-electron flux and evaluate its compatibility with the measurement recently presented by the H.E.S.S. collaboration. To this end, we have MC generated a large sample of the all-electron flux from an ensemble of random distributions of sources. We confirm the non-Gaussian nature of the probability density of fluxes at individual energies previously reported in analytical computations. For the first time, we also consider the correlations between the fluxes at different energies, treating the binned spectrum as a random vector and parametrising its joint distribution with the help of a pair-copula construction. We show that the spectral break observed in the all-electron spectrum by H.E.S.S. and DAMPE is statistically compatible with a distribution of astrophysical sources like supernova remnants or pulsars, but requires a rate smaller than the canonical supernova rate. This important result provides an astrophysical interpretation of the spectrum at TeV energies and allows differentiating astrophysical source models from exotic explanations, like dark matter annihilation. We also critically assess the reliability of using catalogues of known sources to model the electron-positron flux.


title: FliPer: A global measure of power density to estimate surface gravities of main-sequence Solar-like stars and red giants

abstract: Asteroseismology provides global stellar parameters such as masses, radii or surface gravities using the mean global seismic parameters as well as the effective temperature for thousands of low-mass stars $(0.8 M_\odot <M<3 M_\odot)$. This methodology has been successfully applied to stars in which acoustic modes excited by turbulent convection are measured. Other techniques such as the Flicker can also be used to determine stellar surface gravities, but only works for $\log{g}$ above $2.5$ dex. In this work, we present a new metric called FliPer (the acronym stands for Flicker in spectral power density, in opposition to the standard Flicker measurement which is computed in the time domain that is able to extend the range for which reliable surface gravities can be obtained ($0.1<\log{g}<4.6$ dex) without performing any seismic analysis for stars brighter than $\textit{Kp}$ $<$ 14. FliPer takes into account the average variability of a star measured in the power density spectrum in a given range of frequencies. However, FliPer values calculated on several ranges of frequency are required to better characterize a star. Using a large set of asteroseismic targets it is possible to calibrate the behavior of surface gravity with FliPer through machine learning. This calibration made with a random forest regressor covers a wide range of surface gravities from main-sequence stars to subgiants and red giants, with very small uncertainties from $0.04$ to $0.1$ dex. FliPer values can be inserted in automatic global seismic pipelines to either give an estimation of the stellar surface gravity or to assess the quality of the seismic results by detecting any outliers in the obtained $\nu_{max}$ values. FliPer also constrain the surface gravities of main-sequence dwarfs using only long cadence data for which the Nyquist frequency is too low to measure the acoustic-mode properties.


title: A Successful 3D Core-Collapse Supernova Explosion Model

abstract: In this paper, we present the results of our three-dimensional, multi-group, multi-neutrino-species radiation/hydrodynamic simulation using the state-of-the-art code F{\sc{ornax}} of the terminal dynamics of the core of a non-rotating 16-M$_{\odot}$ stellar progenitor. The calculation incorporates redistribution by inelastic scattering, a correction for the effect of many-body interactions on the neutrino-nucleon scattering rates, approximate general relativity (including the effects of gravitational redshifts), velocity-dependent frequency advection, and an implementation of initial perturbations in the progenitor core. The model explodes within $\sim$100 milliseconds of bounce (near when the silicon-oxygen interface is accreted through the temporarily-stalled shock) and by the end of the simulation (here, $\sim$677 milliseconds after bounce) is accumulating explosion energy at a rate of $\sim$2.5$\times$10$^{50}$ ergs s$^{-1}$. The supernova explosion resembles an asymmetrical multi-plume structure, with one hemisphere predominating. The gravitational mass of the residual proto-neutron star at $\sim$677 milliseconds is $\sim$1.42 M$_{\odot}$. Even at the end of the simulation, explosion in most of the solid angle is accompanied by some accretion in an annular fraction at the wasp-like waist of the debris field. The ejecta electron fraction (Y$_e$) is distributed from $\sim$0.48 to $\sim$0.56, with most of the ejecta mass proton-rich. This may have implications for supernova nucleosynthesis, and could have a bearing on the p- and $\nu$p-processes and on the site of the first peak of the r-process. The ejecta spatial distributions of both Y$_e$ and mass density are predominantly in wide-angle plumes and large-scale structures, but are nevertheless quite patchy.


title: NuSTAR reveals that the heavily obscured nucleus of NGC 2785 was the contaminant of IRAS 09104+4109 in the BeppoSAX/PDS hard X-rays

abstract: The search for heavily obscured active galactic nuclei (AGNs) has been revitalized in the last five years by NuSTAR, which has provided a good census and spectral characterization of a population of such objects, mostly at low redshift, thanks to its enhanced sensitivity above 10 keV compared to previous X-ray facilities, and its hard X-ray imaging capabilities. We aim at demonstrating how NGC2785, a local (z=0.009) star-forming galaxy, is responsible, in virtue of its heavily obscured active nucleus, for significant contamination in the non-imaging BeppoSAX/PDS data of the relatively nearby (~17 arcmin) quasar IRAS 09104+4109 (z=0.44), which was originally mis-classified as Compton thick. We analyzed ~71 ks NuSTAR data of NGC2785 using the MYTorus model and provided a physical description of the X-ray properties of the source for the first time. We found that NGC2785 hosts a heavily obscured (NH~3*10^{24} cm^{-2}) nucleus. The intrinsic X-ray luminosity of the source, once corrected for the measured obscuration (L(2-10 keV)~10^{42} erg/s), is consistent within a factor of a few with predictions based on the source mid-infrared flux using widely adopted correlations from the literature. Based on NuSTAR data and previous indications from the Neil Gehrels Swift Observatory (BAT instrument), we confirm that NGC2785, because of its hard X-ray emission and spectral shape, was responsible for at least one third of the 20-100 keV emission observed using the PDS instrument onboard BeppoSAX, originally completely associated with IRAS 09104+4109. Such emission led to the erroneous classification of this source as a Compton-thick quasar, while it is now recognized as Compton thin.


title: Scattering of scalar, electromagnetic and gravitational waves from binary systems

abstract: The direct detection of gravitational waves crowns decades of efforts in the modelling of sources and of increasing detectors' sensitivity. With future third-generation Earth-based detectors or space-based observatories, gravitational-wave astronomy will be at its full bloom. Previously brushed-aside questions on environmental or other systematic effects in the generation and propagation of gravitational waves are now begging for a systematic treatment. Here, we study how electromagnetic and gravitational radiation is scattered by a binary system. Scattering cross-sections, resonances and the effect of an impinging wave on a gravitational-bound binary are worked out for the first time. The ratio between the scattered-wave amplitude and the incident wave can be of order $10^{-5}$ for known pulsars, bringing this into the realm of future gravitational-wave observatories. For currently realistic distribution of compact-object binaries, the interaction cross-section is too small to be of relevance.


title: Strongly Coupled Quintessence

abstract: We present a family of consistent quantum field theories of monodromy quintessence in strong coupling, which can serve as benchmarks in modeling dark energy different from cosmological constant. These theories have discrete gauge symmetries which can protect them from quantum field theory and quantum gravity corrections, both perturbative and nonperturbative. The strong coupling effects, at scales $\ga {\rm mm}^{-1}$, flatten the potential and activate operators with higher powers of derivatives. The predicted equation of state is close to, but not exactly equal to $-1$, thus being within reach of the (near!) future programs to explore the nature of dark energy.


title: Bridging Hubbard Model Physics and Quantum Hall Physics in Trilayer Graphene/h-BN moir\'e superlattice

abstract: The moir\'e superlattice formed by ABC stacked trilayer graphene aligned with a hexagonal boron nitride substrate (TG/h-BN) provides an interesting system where both the bandwidth and the topology can be tuned by an applied perpendicular electric field $D$ . Thus the TG/h-BN system can simulate both Hubbard model physics and nearly flat Chern band physics within one sample. We derive lattice models for both signs of $D$ (which controls the band topology) separately through explicit Wannier orbital construction and mapping of Coulomb interaction. When the bands are topologically trivial, we discuss possible candidates for Mott insulators at integer number of holes per site (labeled as $\nu_T$). These include both broken symmetry states and quantum spin liquid insulators which may be particularly favorable in the vicinity of the Mott transition. We propose feasible experiments to study carefully the bandwidth tuned and the doping tuned Mott metal-insulator transition at both $\nu_T=1$ and $\nu_T=2$. We discuss the interesting possibility of probing experimentally a bandwidth (or doping) controlled continuous Mott transition between a Fermi liquid metal and a quantum spin liquid insulator. Finally we also show that the system has a large valley Zeeman coupling to a small out-of-plane magnetic field, which can be used to control the valley degree of freedom.


title: Harmonic Analysis and Mean Field Theory

abstract: We review some aspects of harmonic analysis for the Euclidean conformal group, including conformally-invariant pairings, the Plancherel measure, and the shadow transform. We introduce two efficient methods for computing these quantities: one based on weight-shifting operators, and another based on Fourier space. As an application, we give a general formula for OPE coefficients in Mean Field Theory (MFT) for arbitrary spinning operators. We apply this formula to several examples, including MFT for fermions and "seed" operators in 4d, and MFT for currents and stress-tensors in 3d.


title: Asteroseismic radii of dwarfs: New accuracy constraints from Gaia DR2 parallaxes

abstract: Precise stellar masses and radii can be determined using asteroseismology, but their accuracy must be tested against independent estimates. Using radii derived from Gaia DR2 parallaxes, we test the accuracy of asteroseismic radii for a sample of 93 dwarfs based on both individual frequency fitting and the seismic scaling relations. Radii from frequency fitting are about 1 per cent smaller than Gaia radii on average; however, this difference may be explained by a negative bias of 30 $\mu$as in the Gaia parallaxes. This indicates that the radii derived from frequency fitting are accurate to within 1 per cent. The scaling relations are found to overestimate radii by more than 5 per cent, compared to the Gaia radii, at the highest temperatures. We demonstrate that this offset is reduced to 3 per cent after applying corrections based on model frequencies to the scaling relation for $\Delta\nu$, but only when the model frequencies are corrected for the surface effect. With corrections to $\Delta\nu$, the scaling relation gives radii accurate to about 2--3 per cent for dwarfs in the temperature range $5400$--$6700$ K. The remaining offset at the highest temperatures may indicate the need for a correction to the scaling relation for $\nu_{\mathrm{max}}$.


title: The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad-Line Region

abstract: We present models of the H$\beta$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},\sigma})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio.


title: A Comparison of Stellar and Gas-Phase Chemical Abundances in Dusty Early-Type Galaxies

abstract: While we observe a large amount of cold interstellar gas and dust in a subset of the early-type galaxy (ETG) population, the source of this material remains unclear. The two main, competing scenarios are external accretion of lower mass, gas-rich dwarfs and internal production from stellar mass loss and/or cooling from the hot interstellar medium (ISM). We test these hypotheses with measurements of the stellar and nebular metallicities of three ETGs (NGC 2768, NGC 3245, and NGC 4694) from new long-slit, high signal-to-noise ratio spectroscopy from the Multi-Object Double Spectographs (MODs) on the Large Binocular Telescope (LBT). These ETGs have modest star formation rates and minimal evidence of nuclear activity. We model the stellar continuum to derive chemical abundances and measure gas-phase abundances with standard nebular diagnostics. We find that the stellar and gas-phase abundances are very similar, which supports internal production and is very inconsistent with the accretion of smaller, lower metallicity dwarfs. All three of these galaxies are also consistent with an extrapolation of the mass-metallicity relation to higher mass galaxies with lower specific star formation rates. The emission line flux ratios along the long-slit, as well as global line ratios clearly indicate that photoionization dominates and ionization by alternate sources including AGN activity, shocks, cosmic rays, dissipative magnetohydrodynamic waves, and single degenerate Type Ia supernovae progenitors do not significantly affect the line ratios.


title: Characterizing Circumgalactic Gas around Massive Ellipticals at z~0.4 - II. Physical Properties and Elemental Abundances

abstract: We present a systematic investigation of the circumgalactic medium (CGM) within projected distances d<160 kpc of luminous red galaxies (LRGs). The sample comprises 16 intermediate-redshift (z=0.21-0.55) LRGs of stellar mass M_star>1e11 M_sun. Combining far-ultraviolet Cosmic Origin Spectrograph spectra from the Hubble Space Telescope and optical echelle spectra from the ground enables a detailed ionization analysis based on resolved component structures of a suite of absorption transitions, including the full HI Lyman series and various ionic metal transitions. By comparing the relative abundances of different ions in individually-matched components, we show that cool gas (T~1e4 K) density and metallicity can vary by more than a factor of ten in in an LRG halo. Specifically, metal-poor absorbing components with <1/10 solar metallicity are seen in 50% of the LRG halos, while gas with solar and super-solar metallicity is also common. These results indicate a complex multiphase structure and poor chemical mixing in these quiescent halos. We calculate the total surface mass density of cool gas, \Sigma_cool, by applying the estimated ionization fraction corrections to the observed HI column densities. The radial profile of \Sigma_cool is best-described by a projected Einasto profile of slope \alpha=1 and scale radius r_s=48 kpc. We find that typical LRGs at z~0.4 contain cool gas mass of M_cool= (1-2) x1e10 M_sun at d<160 kpc (or as much as 4x1e10 M_sun at d<500 kpc), comparable to the cool CGM mass of star-forming galaxies. Furthermore, we show that high-ionization OVI and low-ionization absorption species exhibit distinct velocity profiles, highlighting their different physical origins. We discuss the implications of our findings for the origin and fate of cool gas in LRG halos.


title: Unistructurality of cluster algebras

abstract: We prove that any skew-symmetrizable cluster algebra is unistructural, which is a conjecture by Assem, Schiffler, and Shramchenko. As a corollary, we obtain that a cluster automorphism of a cluster algebra $\mathcal A(\mathcal S)$ is just an automorphism of the ambient field $\mathcal F$ which restricts to a permutation of cluster variables of $\mathcal A(S)$.


title: Sidon sets and 2-caps in $\mathbb{F}_3^n$

abstract: For each natural number $d$, we introduce the concept of a $d$-cap in $\mathbb{F}_3^n$. A subset of $\mathbb{F}_3^n$ is called a $d$-cap if, for each $k = 1, 2, \dots, d$, no $k+2$ of the points lie on a $k$-dimensional flat. This generalizes the notion of a cap in $\mathbb{F}_3^n$. We prove that the $2$-caps in $\mathbb{F}_3^n$ are exactly the Sidon sets in $\mathbb{F}_3^n$ and study the problem of determining the size of the largest $2$-cap in $\mathbb{F}_3^n$.


title: Rigidity of weighted composition operators on $H^p$

abstract: We show that every non-compact weighted composition operator $f \mapsto u\cdot (f\circ\phi)$ acting on a Hardy space $H^p$ for $1 \leq p < \infty$ fixes an isomorphic copy of the sequence space $\ell^p$ and therefore fails to be strictly singular. We also characterize those weighted composition operators on $H^p$ which fix a copy of the Hilbert space $\ell^2$. These results extend earlier ones obtained for unweighted composition operators.


title: Towards Secure Infrastructure-based Cooperative Adaptive Cruise Control

abstract: Cooperative Adaptive Cruise Control (CACC) is a pivotal vehicular application that would allow transportation field to achieve its goals of increased traffic throughput and roadway capacity. This application is of paramount interest to the vehicular technology community with a large body of literature dedicated to research within different aspects of CACC, including but not limited to security with CACC. Of all available literature, the overwhelming focus in on CACC utilizing vehicle-to-vehicle (V2V) communication. In this work, we assert that a qualitative increase in vehicle-to-infrastructure (V2I) and infrastructure-to-vehicle (I2V) involvement has the potential to add greater value to CACC. In this study, we developed a strategy for detection of a denial-of-service (DoS) attack on a CACC platoon where the system edge in the vehicular network plays a central role in attack detection. The proposed security strategy is substantiated with a simulation-based evaluation using the ns-3 discrete event network simulator. Empirical evidence obtained through simulation-based results illustrate successful detection of the DoS attack at four different levels of attack severity using this security strategy.


title: Time preference and information acquisition

abstract: I consider the sequential implementation of a target information structure. I characterize the set of decision time distributions induced by all signal processes that satisfy a per-period learning capacity constraint. I find that all decision time distributions have the same expectation, and the maximal and minimal elements by mean-preserving spread order are deterministic distribution and exponential distribution. The result implies that when time preference is risk loving (e.g. standard or hyperbolic discounting), Poisson signal is optimal since it induces the most risky exponential decision time distribution. When time preference is risk neutral (e.g. constant delay cost), all signal processes are equally optimal.


title: Singular Hochschild cohomology via the singularity category

abstract: We show that for a noetherian algebra $A$ whose bounded dg derived category is smooth, the singular Hochschild cohomology (=Tate--Hochschild cohomology) is isomorphic, as a graded algebra, to the Hochschild cohomology of the dg singularity category of $A$. The existence of such an isomorphism is suggested by recent work of Zhengfang Wang.


title: Holomorphic subgraph reduction of higher-point modular graph forms

abstract: Modular graph forms are a class of modular covariant functions which appear in the genus-one contribution to the low-energy expansion of closed string scattering amplitudes. Modular graph forms with holomorphic subgraphs enjoy the simplifying property that they may be reduced to sums of products of modular graph forms of strictly lower loop order. In the particular case of dihedral modular graph forms, a closed form expression for this holomorphic subgraph reduction was obtained previously by D'Hoker and Green. In the current work, we extend these results to trihedral modular graph forms. Doing so involves the identification of a modular covariant regularization scheme for certain conditionally convergent sums over discrete momenta, with some elements of the sum being excluded. The appropriate regularization scheme is identified for any number of exclusions, which in principle allows one to perform holomorphic subgraph reduction of higher-point modular graph forms with arbitrary holomorphic subgraphs.


title: A mechanism for holography for non-interacting fields on anti-de Sitter spacetimes

abstract: In the setting of non-interacting Klein-Gordon fields on asymptotically anti-de Sitter spacetimes, we show that algebras of observables localized in a neighborhood of the boundary are subalgebras of a boundary algebra. The underlying mechanism is directly related to holography for classical fields. In particular, the proof relies on unique continuation theorems at the conformal boundary.


title: Enhanced Network Embeddings via Exploiting Edge Labels

abstract: Network embedding methods aim at learning low-dimensional latent representation of nodes in a network. While achieving competitive performance on a variety of network inference tasks such as node classification and link prediction, these methods treat the relations between nodes as a binary variable and ignore the rich semantics of edges. In this work, we attempt to learn network embeddings which simultaneously preserve network structure and relations between nodes. Experiments on several real-world networks illustrate that by considering different relations between different node pairs, our method is capable of producing node embeddings of higher quality than a number of state-of-the-art network embedding methods, as evaluated on a challenging multi-label node classification task.


title: Dynamics of weakly interacting bosons in optical lattices with flux

abstract: Realization of strong synthetic magnetic fields in driven optical lattices has enabled implementation of topological bands in cold-atom setups. A milestone has been reached by a recent measurement of a finite Chern number based on the dynamics of incoherent bosonic atoms. The measurements of the quantum Hall effect in semiconductors are related to the Chern-number measurement in a cold-atom setup, however, the design and complexity of the two types of measurements are quite different. Motivated by these recent developments, we investigate the dynamics of weakly interacting incoherent bosons in a two-dimensional driven optical lattice exposed to an external force, which provides a direct probe of the Chern number. We consider a realistic driving protocol in the regime of high driving frequency and focus on the role of weak repulsive interactions. We find that interactions lead to the redistribution of atoms over topological bands both through the conversion of interaction energy into kinetic energy during the expansion of the atomic cloud and due to an additional heating. Remarkably, we observe that the moderate atomic repulsion facilitates the measurement by flattening the distribution of atoms in the quasi-momentum space. Our results also show that weak interactions can suppress the contribution of some higher-order non-topological terms in favor of the topological part of the effective model.


title: Control of long-range correlations in turbulence

abstract: The character of turbulence depends on where it develops. Turbulence near boundaries, for instance, is different than in a free stream. To elucidate the differences between flows, it is instructive to vary the structure of turbulence systematically, but there are few ways of stirring turbulence that make this possible. In other words, an experiment typically examines either a boundary layer or a free stream, say, and the structure of the turbulence is fixed by the geometry of the experiment. We introduce a new active grid with many more degrees of freedom than previous active grids. The additional degrees of freedom make it possible to control various properties of the turbulence. We show how long-range correlations in the turbulent velocity fluctuations can be shaped by changing the way the active grid moves. Specifically, we show how not only the correlation length but also the detailed shape of the correlation function depends on the correlations imposed in the motions of the grid. Until now, large-scale structure had not been adjustable in experiments. This new capability makes possible new systematic investigations into turbulence dissipation and dispersion, for example, and perhaps in flows that mimic features of boundary layers, free streams, and flows of intermediate character.


title: IL-Net: Using Expert Knowledge to Guide the Design of Furcated Neural Networks

abstract: Deep neural networks (DNN) excel at extracting patterns. Through representation learning and automated feature engineering on large datasets, such models have been highly successful in computer vision and natural language applications. Designing optimal network architectures from a principled or rational approach however has been less than successful, with the best successful approaches utilizing an additional machine learning algorithm to tune the network hyperparameters. However, in many technical fields, there exist established domain knowledge and understanding about the subject matter. In this work, we develop a novel furcated neural network architecture that utilizes domain knowledge as high-level design principles of the network. We demonstrate proof-of-concept by developing IL-Net, a furcated network for predicting the properties of ionic liquids, which is a class of complex multi-chemicals entities. Compared to existing state-of-the-art approaches, we show that furcated networks can improve model accuracy by approximately 20-35%, without using additional labeled data. Lastly, we distill two key design principles for furcated networks that can be adapted to other domains.


title: Nuclear Processes in Other Universes: Varying the Strength of the Weak Force

abstract: Motivated by the possibility that the laws of physics could be different in other regions of space-time, we consider nuclear processes in universes where the weak interaction is either stronger or weaker than observed. We focus on the physics of both Big Bang Nucleosynthesis (BBN) and stellar evolution. For sufficiently ineffective weak interactions, neutrons do not decay during BBN, and the baryon-to-photon ratio $\eta$ must be smaller in order for protons to survive without becoming incorporated into larger nuclei. For stronger weak interactions, neutrons decay before the onset of BBN, and the early universe is left with nearly a pure hydrogen composition. We then consider stellar structure and evolution for the different nuclear compositions resulting from BBN, a wide range of weak force strengths, and the full range of stellar masses for a given universe. We delineate the range of this parameter space that supports working stars, along with a determination of the dominant nuclear reactions over the different regimes. Deuterium burning dominates the energy generation in stars when the weak force is sufficiently weak, whereas proton-proton burning into helium-3 dominates for the regime where the weak force is much stronger than in our universe. Although stars in these universes are somewhat different, they have comparable surface temperatures, luminosities, radii, and lifetimes, so that a wide range of such universes remain potentially habitable.


title: The past, present and future of the heavier electroweakinos in the light of LHC and other data

abstract: The aim of this paper is to showcase the novel multilepton ($nl + E\!\!\!\!/_T$, n = 3 - 5) signals, hitherto unexplored at the LHC, arising from the heavier electroweakinos, in several generic pMSSMs at the upcoming LHC experiments. We first briefly review our old constraints on the full electroweakino sector of these models, containing both lighter and heavier sparticles, using the ATLAS trilepton data from the LHC Run I. Next we derive new stronger constraints on this sector for the first time using the ATLAS Run II data. We identify some benchmark points and explore the prospect of observing multilepton events in future LHC experiments. Our focus is on the channels with $n > 3$ which are the hallmarks of the heavier electroweakinos. If the spectrum of the lighter electroweakinos is compressed, these signals might very well be the discovery channels of the electroweakinos at the high luminosity LHC. We also discuss the implications of the new LHC constraints for the observed dark matter relic density of the universe, the measured value of the anomalous magnetic moment of the muon and the dark matter direct detection experiments.


title: Irrational Toric Varieties

abstract: Classical toric varieties are among the simplest objects in algebraic geometry. They arise in an elementary fashion as varieties parametrized by monomials whose exponents are a finite subset $\mathcal{A}$ of $\mathbb{Z}^n$. They may also be constructed from a rational fan $\Sigma$ in $\mathbb{R}^n$. The combinatorics of the set $\mathcal{A}$ or fan $\Sigma$ control the geometry of the associated toric variety. These toric varieties have an action of an algebraic torus with a dense orbit. Applications of algebraic geometry in geometric modeling and algebraic statistics have long studied the nonnegative real part of a toric variety as the main object, where the set $\mathcal{A}$ may be an arbitrary set in $\mathbb{R}^n$. These are called irrational affine toric varieties. This theory has been limited by the lack of a construction of an irrational toric variety from an arbitrary fan in $\mathbb{R}^n$. We construct a theory of irrational toric varieties associated to arbitrary fans. These are $(\mathbb{R}_>)^n$-equivariant cell complexes dual to the fan. Such an irrational toric variety is projective (may be embedded in a simplex) if and only if its fan is the normal fan of a polytope, and in that case, the toric variety is homeomorphic to that polytope. We use irrational toric varieties to show that the space of Hausdorff limits of translates an irrational toric variety associated to a finite subset $\mathcal{A}$ of $\mathbb{R}^n$ is homeomorphic to the secondary polytope of $\mathcal{A}$.