title: The Boundary between Gas-rich and Gas-poor Planets

abstract: Sub-Saturns straddle the boundary between gas-rich Jupiters and gas-poor super-Earths/sub-Neptunes. Their large radii (4--8$R_\oplus$) suggest that their gas-to-core mass ratios range $\sim$0.1--1.0. With their envelopes as massive as their cores, sub-Saturns are just on the verge of runaway gas accretion; they are expected to be significantly less populous than gas giants. Yet, the observed occurrence rates of sub-Saturns and Jupiters are comparable within $\sim$100 days. We show that in these inner regions of planetary systems, the growth of sub-Saturns/Jupiters is ultimately limited by local and global hydrodynamic flows---runaway accretion terminates and the formation of gas giants is suppressed. We derive a simple analytic formula for the local hydrodynamic accretion rate---an expression that has been previously reported only as an empirical fit to numerical simulations. Evolving simultaneously the background disk gas and the gas accretion onto planetary cores, we find that both the ubiquity of super-Earths/sub-Neptunes and the rarity of gas-rich planets are best explained if an underlying core mass distribution is peaked at $\sim$4.3$M_\oplus$. Within a finite disk lifetime $\sim$10 Myrs, massive cores ($\gtrsim 10M_\oplus$) can become either gas-poor or gas-rich depending on when they assemble but smaller cores ($\lesssim 10M_\oplus$) can only become gas-poor. This wider range of possible outcomes afforded by more massive cores may explain why metal-rich stars harbor a more diverse set of planets.


title: MIND THE GAP: Is The Too Big To Fail Problem Resolved?

abstract: The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $\Lambda$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $\Lambda$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected.


title: Topological phase diagram of a three-terminal Josephson junction: From the conventional to the Majorana regime

abstract: We study the evolution of averaged transconductances in three-terminal Josephson junctions when the superconducting leads are led throughout a topological phase transition from an $s$-wave to a $p$-wave (Majorana) phase by an in-plane magnetic field $B_x$. We provide a complete description of this transition as a function of $B_x$ and a magnetic flux $\Phi$ threading the junction. For that we use a spinful model within a formalism that allows to treat on equal footing the contribution to the transconductance from both the Andreev subgap levels and the continuum spectrum. We unveil a fractionalization in the quantization of the transconductance due to the presence of Majorana quasiparticles, reflecting the effective pumping of half a Cooper pair charge in the $p$-wave regime.


title: Empirical case for two pseudogaps in cuprate superconductors

abstract: Superconductivity in cuprates is achieved by doping holes into a correlated charge-transfer insulator. While the correlated character of the parent insulator is now understood, there is no accepted theory for the "normal" state of the doped insulator. I present a mostly empirical analysis of a large range of experimental characterizations, making the case for two pseudogaps: (1) a large pseudogap resulting from the competition between the energy of superexchange-coupled local Cu moments and the kinetic energy of doped holes; (2) a small pseudogap that results from dopant disorder and consequent variations in local charge density, leading to a distribution of local superconducting onset temperatures. The large pseudogap closes as hole kinetic energy dominates at higher doping and the dynamic antiferromagnetic correlations become overdamped. Establishing spatially-homogeneous $d$-wave superconductivity is limited by those regions with the weakest superconducting phase coherence, which tends to be limited by low-energy spin fluctuations. The magnitude of the small pseudogap is correlated with the doping-dependent energy $E_{\rm cross}$ associated with the neck of the hour-glass dispersion of spin excitations. The consequences of this picture are discussed.


title: Atmospheres and UV Environments of Earth-like Planets Throughout Post-Main Sequence Evolution

abstract: During the post-main sequence phase of stellar evolution the orbital distance of the habitable zone, which allows for liquid surface water on terrestrial planets, moves out past the system's original frost line, providing an opportunity for outer planetary system surface habitability. We use a 1D coupled climate/photochemistry code to study the impact of the stellar environment on the planetary atmospheres of Earth-like planets/moons throughout its time in the post-main sequence habitable zone. We also explore the ground UV environments of such planets/moons and compare them to Earth's. We model the evolution of star-planet systems with host stars ranging from 1.0 to 3.5 M$_\odot$ throughout the post-main sequence, calculating stellar mass loss and its effects on planetary orbital evolution and atmospheric erosion. The maximum amount of time a rocky planet can spend continuously in the evolving post-MS habitable zone ranges between 56 and 257 Myr for our grid stars. Thus, during the post-MS evolution of their host star, subsurface life on cold planets and moons could become remotely detectable once the initially frozen surface melts. Frozen planets or moons, like Europa in our Solar System, experience a relatively stable environment on the horizontal branch of their host stars' evolution for millions of years.


title: Translating solitons $C^1-$asymptotic to two half-hyperplanes

abstract: We prove that the hyperplanes parallel to ${\bf e}_{n+1}$ are the unique examples of translating solitons $C^1-$asymptotic to two half-hyperplanes outside a vertical cylinder in $\R^{n+1}$.


title: Tailored tunnel magnetoresistance response in three ultrathin chromium trihalides

abstract: Materials that demonstrate large magnetoresistance have attracted significant interest for many decades. Recently, extremely large tunnel magnetoresistance (TMR) has been reported by several groups across ultrathin CrI$_3$ by exploiting the weak antiferromagnetic coupling between adjacent layers. Here, we report a comparative study of TMR in all three chromium trihalides (CrX$_3$, X= Cl, Br, or I) in the two-dimensional limit. As the materials exhibit different transition temperatures and interlayer magnetic ordering in the ground state, tunneling measurements allow for an easy determination of the field-temperature phase diagram for the three systems. By changing sample thickness and biasing conditions, we then demonstrate how to maximize and further tailor the TMR response at different temperatures for each material. In particular, near the magnetic transition temperature, TMR is non-saturating up to the highest fields measured for all three compounds owing to the large, field-induced exchange coupling.


title: On the Ultrametricity Property in Random Field Ising Models

abstract: In this paper we show that the ultrametricity property remains valid in random field Ising models with independent disorder whenever the field strength is a small perturbation.

categories: math-ph cond-mat.stat-mech math.MP math.PR


title: New Duality Results for Evenly Convex Optimization Problems

abstract: We present new results on optimization problems where the involved functions are evenly convex. By means of a generalized conjugation scheme and the perturbation theory introduced by Rockafellar, we propose an alternative dual problem for a general optimization one defined on a separated locally convex topological space. Sufficient conditions for converse and total duality involving the even convexity of the perturbation function and $c$-subdifferentials are given. Formulae for the $c$-subdifferential and biconjugate of the objective function of a general optimization problem are provided, too. We also characterize the total duality also by means of the saddle-point theory for a notion of Lagrangian adapted to the considered framework.


title: Counting Induced Subgraphs: An Algebraic Approach to #W[1]-hardness

abstract: We study the problem #IndSub(P) of counting all induced subgraphs of size k in a graph G that satisfy the property P. This problem was introduced by Jerrum and Meeks and shown to be #W[1]-hard when parameterized by k for some families of properties P including, among others, connectivity [JCSS 15] and even- or oddness of the number of edges [Combinatorica 17]. Very recently [IPEC 18], two of the authors introduced a novel technique for the complexity analysis of #IndSub(P), inspired by the "topological approach to evasiveness" of Kahn, Saks and Sturtevant [FOCS 83] and the framework of graph motif parameters due to Curticapean, Dell and Marx [STOC 17], allowing them to prove hardness of a wide range of properties P. In this work, we refine this technique for graph properties that are non-trivial on edge-transitive graphs with a prime power number of edges. In particular, we fully classify the case of monotone bipartite graph properties: It is shown that, given any graph property P that is closed under the removal of vertices and edges, and that is non-trivial for bipartite graphs, the problem #IndSub(P) is #W[1]-hard and cannot be solved in time f(k)*n^{o(k)} for any computable function f, unless the Exponential Time Hypothesis fails. This holds true even if the input graph is restricted to be bipartite and counting is done modulo a fixed prime. A similar result is shown for properties that are closed under the removal of edges only.


title: Measurement of the Photon-Plasmon Coupling Phase

abstract: Scattering processes have played a crucial role in the development of quantum theory. In the field of optics, scattering phase shifts have been utilized to unveil interesting forms of light-matter interactions. Here, we investigate the mode-coupling phase of single photons to surface plasmon polaritons in a quantum plasmonic tritter. We observe that the coupling process induces a phase jump that occurs when photons scatter into surface plasmons and vice versa. This interesting coupling phase dynamics is of particular relevance for quantum plasmonic experiments. Furthermore, it is demonstrated that this photon-plasmon interaction can be modeled through a quantum-mechanical tritter. We show that the visibility of a double-slit and a triple-slit interference patterns are convenient observables to characterize the interaction at a slit and determine the coupling phase. Our accurate and simple model of the interaction, validated by simulations and experiments, has important implications not only for quantum plasmonic interference effects, but is also advantageous to classical applications.


title: ECG Reconstruction via PPG: A Pilot Study

abstract: In this paper, the relation between electrocardiogram (ECG) and photoplethysmogram (PPG) signals is studied, and the waveform of ECG is inferred via the PPG signals. In order to address this inverse problem, a transform is proposed to map the discrete cosine transform (DCT) coefficients of each PPG cycle to those of the corresponding ECG cycle. The resulting DCT coefficients of the ECG cycle are inversely transformed to obtain the reconstructed ECG waveform. The proposed method is evaluated on a benchmark dataset of subjects with a variety of combinations of age and weight. Experimental results show that the proposed method can achieve a high accuracy at 0.98 in averaged correlation.


title: On the quasi-isometric rigidity of graphs of surface groups

abstract: Let $\Gamma$ be a word hyperbolic group with a cyclic JSJ decomposition that has only rigid vertex groups, which are all fundamental groups of closed surface groups. We show that any group $H$ quasi-isometric to $\Gamma$ is abstractly commensurable with $\Gamma$.


title: Gravitational waves from fermion production during axion inflation

abstract: We present analytic results for the gravitational wave power spectrum induced in models where the inflaton is coupled to a fermionic pseudocurrent. We show that although such a coupling creates helically polarized fermions, the polarized component of the resulting gravitational waves is parametrically suppressed with respect to the non-polarized one. We also show that the amplitude of the gravitational wave signal associated to this production cannot exceed that generated by the standard mechanism of amplification of vacuum fluctuations. We previously found that this model allows for a regime in which the backreaction of the produced fermions allows for slow-roll inflation even for a steep inflaton potential, and still leads to Gaussian primordial scalar perturbations. The present analysis shows that this regime also results in a gravitational wave signal compatible with the current bounds.


title: TurbuStat: Turbulence Statistics in Python

abstract: We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions.


title: Frenesy: time-symmetric dynamical activity in nonequilibria

abstract: We review the concept of dynamical ensembles in nonequilibrium statistical mechanics as specified from an action functional or Lagrangian on spacetime. There, under local detailed balance, the breaking of time-reversal invariance is quantified via the entropy flux, and we revisit some of the consequences for fluctuation and response theory. Frenesy is the time-symmetric part of the path-space action with respect to a reference process. It collects the variable quiescence and dynamical activity as function of the system's trajectory, and as has been introduced under different forms in studies of nonequilibria. We discuss its various realizations for physically inspired Markov jump and diffusion processes and why it matters a good deal for nonequilibrium physics. This review then serves also as an introduction to the exploration of frenetic contributions in nonequilibrium phenomena.


title: Fisher-Information Efficient Metagratings for Transverse Displacement Metrology at the Measurement Independent Shot-Noise Limit

abstract: We derive the measurement independent precision bound for transverse displacement metrology using a classical source. Using tools from quantum estimation theory and metamaterial designs, we prove that the bound is achieved efficiently using optimized metagratings by detecting only a very small amount of photons used to probe the displacement, and we also reveal the direct link between the resonant property in the unit cell and the conditions to achieve the bound: one with the 0th dipole resonance and the other one with the anapole condition of the 1st dipole.


title: The Cauchy problem for a fractional conservation laws driven by L\'evy noise

abstract: In this article, we explore some of the main mathematical problems connected to multidimensional fractional conservation laws driven by L\'evy processes. Making use of an adapted entropy formulation, a result of existence and uniqueness of a solution is established. Moreover, using bounded variation (BV) estimates for vanishing viscosity approximations, we derive an explicit continuous dependence estimate on the nonlinearities of the entropy solutions under the assumption that L\'evy noise depends only on the solution. This result is used to show the error estimate for the stochastic vanishing viscosity method. Furthermore, we establish a result on vanishing non-local regularization of scalar stochastic conservation laws.


title: Gamma-Ray Bursts: Theoretical Issues and Developments

abstract: I discuss some aspects of the evolution of the standard GRB model, emphasizing various theoretical developments in the last decade, and review the impact of some of the most recent observational discoveries and the new challenges they pose in the expanding realm of multi-messenger astrophysics.


title: Chunkflow: Distributed Hybrid Cloud Processing of Large 3D Images by Convolutional Nets

abstract: It is now common to process volumetric biomedical images using 3D Convolutional Networks (ConvNets). This can be challenging for the teravoxel and even petavoxel images that are being acquired today by light or electron microscopy. Here we introduce chunkflow, a software framework for distributing ConvNet processing over local and cloud GPUs and CPUs. The image volume is divided into overlapping chunks, each chunk is processed by a ConvNet, and the results are blended together to yield the output image. The frontend submits ConvNet tasks to a cloud queue. The tasks are executed by local and cloud GPUs and CPUs. Thanks to the fault-tolerant architecture of Chunkflow, cost can be greatly reduced by utilizing cheap unstable cloud instances. Chunkflow currently supports PyTorch for GPUs and PZnet for CPUs. To illustrate its usage, a large 3D brain image from serial section electron microscopy was processed by a 3D ConvNet with a U-Net style architecture. Chunkflow provides some chunk operations for general use, and the operations can be composed flexibly in a command line interface.


title: The Masses of Supernova Remnant Progenitors in M83

abstract: We determine the ages of the young, resolved stellar populations at the locations of 237 optically-identified supernova remnants in M83. These age distributions put constraints on the progenitor masses of the supernovae that produced 199 of the remnants. The other 38 show no evidence for having a young progenitor and are therefore good Type Ia SNR candidates. Starting from Hubble Space Telescope broadband imaging, we measured resolved stellar photometry of seven archival WFC3/UVIS fields in F336W, F438W, and F814W. We generate color-magnitude diagrams of the stars within 50 pc of each SNR and fit them with stellar evolution models to obtain the population ages. From these ages we infer the progenitor mass that corresponds to the lifetime of the most prominent age that is $<$50 Myr. In this sample, there are 47 SNRs with best-fit progenitor masses $>$15 M$_{\odot}$, and 5 of these are $>$15 M$_{\odot}$ at 84% confidence. This is the largest collection of high-mass progenitors to date, including our highest-mass progenitor inference found so far, with a constraint of $<$8 Myr. Overall, the distribution of progenitor masses has a power-law index of $-3.0^{+0.2}_{-0.7}$, steeper than Salpeter initial mass function ($-2.35$). It remains unclear whether the reason for the low number of high-mass progenitors is due to the difficulty of finding and measuring such objects or because only a fraction of very massive stars produce supernovae.


title: Quantum geometry of moduli spaces of local systems and representation theory

abstract: Let G be a split semi-simple adjoint group, and S a colored decorated surface, given by an oriented surface with punctures, special boundary points, and a specified collection of boundary intervals. We introduce a moduli space P(G,S) parametrizing G-local system on S with some boundary data, and prove that it carries a cluster Poisson structure, equivariant under the action of the cluster modular group M(G,S), containing the mapping class group of S, the group of outer automorphisms of G, and the product of Weyl / braid groups over punctures / boundary components. We prove that the dual moduli space A(G,S) carries a M(G,S)-equivariant cluster structure, and the pair (A(G,S), P(G,S)) is a cluster ensemble. These results generalize the works of V. Fock & the first author, and of I. Le. We quantize cluster Poisson varieties X for any Planck constant h s.t. h>0 or |h|=1. First, we define a *-algebra structure on the Langlands modular double A(h; X) of the algebra of functions on X. We construct a principal series of representations of the *-algebra A(h; X), equivariant under a unitary projective representation of the cluster modular group M(X). This extends works of V. Fock and the first author when h>0. Combining this, we get a M(G,S)-equivariant quantization of the moduli space P(G,S), given by the *-algebra A(h; P(G,S)) and its principal series representations. We construct realizations of the principal series *-representations. In particular, when S is punctured disc with two special points, we get a principal series *-representations of the Langlands modular double of the quantum group Uq(g). We conjecture that there is a nondegenerate pairing between the local system of coinvariants of oscillatory representations of the W-algebra and the one provided by the projective representation of the mapping class group of S.


title: Photon echoes in optically dense media

abstract: Coherent nonlinear multi-pulse processes, nonlinear waves and echo effects in resonant media are the topical problems of modern optics and important tools of coherent spectroscopy and quantum information science. We generalize the McCall-Hahn area theorem to the formation of an arbitrary photon echo generated during the multi-pulse excitation of the optically dense resonant media. The derived theorem made it possible to reveal the nonlinear mechanism of generation and evolution of the photon echo signals inside the media after a two-pulse excitation. We find that a series of self-reviving echo signals with total area of $2\pi$ or $0\pi$ is excited and propagates in the media depth, with each pulse having an individual area less than $\pi$. The resulting echo pulse train is a new alternative to the well-known soliton or breather. The developed pulse-area approach paves the way for more precise coherent spectroscopy, studies of different photon echo signals and quantum control of light pulses in the optically dense media.


title: Counting perfect matchings and the eight-vertex model

abstract: We study the approximation complexity of the partition function of the eight-vertex model on general 4-regular graphs. For the first time, we relate the approximability of the eight-vertex model to the complexity of approximately counting perfect matchings, a central open problem in this field. Our results extend those in arXiv:1811.03126 [cs.CC]. In a region of the parameter space where no previous approximation complexity was known, we show that approximating the partition function is at least as hard as approximately counting perfect matchings via approximation-preserving reductions. In another region of the parameter space which is larger than the previously known FPRASable region, we show that computing the partition function can be reduced to (with or without approximation) counting perfect matchings. Moreover, we give a complete characterization of nonnegatively weighted (not necessarily planar) 4-ary matchgates, which has been open for several years. The key ingredient of our proof is a geometric lemma. We also identify a region of the parameter space where approximating the partition function on planar 4-regular graphs is feasible but on general 4-regular graphs is equivalent to approximately counting perfect matchings. To our best knowledge, these are the first problems of this kind.


title: Balanced Crossover Operators in Genetic Algorithms

abstract: In several combinatorial optimization problems arising in cryptography and design theory, the admissible solutions must often satisfy a balancedness constraint, such as being represented by bitstrings with a fixed number of ones. For this reason, several works in the literature tackling these optimization problems with Genetic Algorithms (GA) introduced new balanced crossover operators which ensure that the offspring has the same balancedness characteristics of the parents. However, the use of such operators has never been thoroughly motivated, except for some generic considerations about search space reduction. In this paper, we undertake a rigorous statistical investigation on the effect of balanced and unbalanced crossover operators against three optimization problems from the area of cryptography and coding theory: nonlinear balanced Boolean functions, binary Orthogonal Arrays (OA) and bent functions. In particular, we consider three different balanced crossover operators (each with two variants: "left-to-right" and "shuffled"), two of which have never been published before, and compare their performances with classic one-point crossover. We are able to confirm that the balanced crossover operators performs better than all three balanced crossover operators. Furthermore, in two out of three crossovers, the "left-to-right" version performs better than the "shuffled" version.


title: On laws exhibiting universal ordering under stochastic restart

abstract: For each of (i) arbitrary stochastic reset, (ii) deterministic reset with arbitrary period, (iii) reset at arbitrary constant rate, and then in the sense of either (a) first-order stochastic dominance or (b) expectation (i.e. for each of the six possible combinations of the preceding), those laws of random times are precisely characterized that are rendered no bigger [rendered no smaller; left invariant] by all possible restart laws (within the classes (i), (ii), (iii), as the case may be). Partial results in the same vein for reset with branching are obtained. In particular it is found that deterministic and arbitrary stochastic restart lead to the same characterizations, but this equivalence fails to persist for exponential (constant-rate) reset.


title: Modelling the radiation pattern of a dual circular polarization system

abstract: We present the electromagnetic model of a dual circular polarization antenna-feed system, consisting of a corrugated feedhorn, a polarizer and an orthomode transducer. This model was developed for the passive front-end implemented in the Q-band receivers of the STRIP instrument of the Large Scale Polarization Explorer experiment. Its applicability, however, is completely general. The model has been implemented by superposing the response of two linearly polarized feedhorns with a $\pi$/2 phase difference, thus taking into account for the effect of the polarizer, that behaves differently for the two polarizations of the incoming electric field. The model has been verified by means of radiation pattern measurements performed in the anechoic chamber at the Physics Department of the University of Milan. We measured both arms of the orthomode transducer, in order to check that the diagram at one port is the 90_ rotation of the diagram at the other port. Simulations and measurements show an agreement at the level of fraction of a dB up to the first sidelobes, thus confirming the model.


title: Springer theory for symplectic Galois groups

abstract: A classical and beautiful story in geometric representation theory is the construction by Springer of an action of the Weyl group on the cohomology of the fibres of the Springer resolution of the nilpotent cone. We establish a natural extension of Springer's theory to arbitrary symplectic resolutions of conical symplectic singularities. We analyse features of the action in the case of affine quiver varieties, constructing Weyl group actions on the cohomology of $ADE$ quiver varieties, and also consider "symplectically dual" examples arising from slices in the affine Grassmannian. Along the way, we document some basic features of the symplectic geometry of quiver varieties.


title: Information Based Centralization of Locomotion in Animals and Robots

abstract: Movement in biology is often achieved with distributed control of coupled subcomponents, e.g. muscles and limbs. Coupling could range from weak and local, i.e. decentralized, to strong and global, i.e. centralized. We developed a model-free measure of centralization that compares information shared between control signals and both global and local states. A second measure, called co-information, quantifies the net redundant information the control signal shares with both states. We first validate our measures through simulations of coupled oscillators and show that it successfully reconstructs the shift from low to high coupling strengths. We then measure centralization in freely running cockroaches. Surprisingly, extensor muscle activity in the middle leg is more informative of movements of all legs combined than the movements of that particular leg. Cockroach centralization successfully recapitulates a specific model of a strongly coupled oscillator network previously used to model cockroach leg kinematics. When segregated by stride frequency, slower cockroach strides exhibit more shared information per stride between control and output states than faster strides, indicative of an information bandwidth limitation. However, centralization remains consistent between the two groups. We then used a robotic model to show that centralization can be affected by mechanical coupling independent of neural coupling. The mechanically coupled bounding gait is decentralized and becomes more decentralized as mechanical coupling decreases while internal parameters of control remain constant. The results of these systems span a design space of centralization and co-information that can be used to test biological hypotheses and advise the design of robotic control.


title: Comparing Samples from the $\mathcal{G}^0$ Distribution using a Geodesic Distance

abstract: The $\mathcal{G}^0$ distribution is widely used for monopolarized SAR image modeling because it can characterize regions with different degree of texture accurately. It is indexed by three parameters: the number of looks (which can be estimated for the whole image), a scale parameter and a texture parameter. This paper presents a new proposal for comparing samples from the $\mathcal{G}^0$ distribution using a Geodesic Distance (GD) as a measure of dissimilarity between models. The objective is quantifying the difference between pairs of samples from SAR data using both local parameters (scale and texture) of the $\mathcal{G}^0$ distribution. We propose three tests based on the GD which combine the tests presented in~\cite{GeodesicDistanceGI0JSTARS}, and we estimate their probability distributions using permutation methods.


title: Natural Language Interactions in Autonomous Vehicles: Intent Detection and Slot Filling from Passenger Utterances

abstract: Understanding passenger intents and extracting relevant slots are important building blocks towards developing contextual dialogue systems for natural interactions in autonomous vehicles (AV). In this work, we explored AMIE (Automated-vehicle Multi-modal In-cabin Experience), the in-cabin agent responsible for handling certain passenger-vehicle interactions. When the passengers give instructions to AMIE, the agent should parse such commands properly and trigger the appropriate functionality of the AV system. In our current explorations, we focused on AMIE scenarios describing usages around setting or changing the destination and route, updating driving behavior or speed, finishing the trip and other use-cases to support various natural commands. We collected a multi-modal in-cabin dataset with multi-turn dialogues between the passengers and AMIE using a Wizard-of-Oz scheme via a realistic scavenger hunt game activity. After exploring various recent Recurrent Neural Networks (RNN) based techniques, we introduced our own hierarchical joint models to recognize passenger intents along with relevant slots associated with the action to be performed in AV scenarios. Our experimental results outperformed certain competitive baselines and achieved overall F1 scores of 0.91 for utterance-level intent detection and 0.96 for slot filling tasks. In addition, we conducted initial speech-to-text explorations by comparing intent/slot models trained and tested on human transcriptions versus noisy Automatic Speech Recognition (ASR) outputs. Finally, we compared the results with single passenger rides versus the rides with multiple passengers.


title: Weighted estimates for the Bergman projection on the Hartogs triangle

abstract: We apply modern techniques of dyadic harmonic analysis to obtain sharp estimates for the Bergman projection in weighted Bergman spaces. Our main theorem focuses on the Bergman projection on the Hartogs triangle. The estimates of the operator norm are in terms of a Bekoll\'e-Bonami type constant. As an application of the results obtained, we give, for example, an upper bound for the $L^p$ norm of the Bergman projection on the generalized Hartogs triangle $\mathbb H_{m/n}$ in $\mathbb C^2$.


title: Relative-error inertial-relaxed inexact versions of Douglas-Rachford and ADMM splitting algorithms

abstract: This paper derives new inexact variants of the Douglas-Rachford splitting method for maximal monotone operators and the alternating direction method of multipliers (ADMM) for convex optimization. The analysis is based on a new inexact version of the proximal point algorithm that includes both an inertial step and overrelaxation. We apply our new inexact ADMM method to LASSO and logistic regression problems and obtain somewhat better computational performance than earlier inexact ADMM methods.


title: Fine-Grained Named Entity Recognition using ELMo and Wikidata

abstract: Fine-grained Named Entity Recognition is a task whereby we detect and classify entity mentions to a large set of types. These types can span diverse domains such as finance, healthcare, and politics. We observe that when the type set spans several domains the accuracy of the entity detection becomes a limitation for supervised learning models. The primary reason being the lack of datasets where entity boundaries are properly annotated, whilst covering a large spectrum of entity types. Furthermore, many named entity systems suffer when considering the categorization of fine grained entity types. Our work attempts to address these issues, in part, by combining state-of-the-art deep learning models (ELMo) with an expansive knowledge base (Wikidata). Using our framework, we cross-validate our model on the 112 fine-grained entity types based on the hierarchy given from the Wiki(gold) dataset.


title: Understanding the efficacy, reliability and resiliency of computer vision techniques for malware detection and future research directions

abstract: My research lies in the intersection of security and machine learning. This overview summarizes one component of my research: combining computer vision with malware exploit detection for enhanced security solutions. I will present the perspectives of efficacy, reliability and resiliency to formulate threat detection as computer vision problems and develop state-of-the-art image-based malware classification. Representing malware binary as images provides a direct visualization of data samples, reduces the efforts for feature extraction, and consumes the whole binary for holistic structural analysis. Employing transfer learning of deep neural networks effective for large scale image classification to malware classification demonstrates superior classification efficacy compared with classical machine learning algorithms. To enhance reliability of these vision-based malware detectors, interpretation frameworks can be constructed on the malware visual representations and useful for extracting faithful explanation, so that security practitioners have confidence in the model before deployment. In cyber-security applications, we should always assume that a malware writer constantly modifies code to bypass detection. Addressing the resiliency of the malware detectors is equivalently important as efficacy and reliability. Via understanding the attack surfaces of machine learning models used for malware detection, we can greatly improve the robustness of the algorithms to combat malware adversaries in the wild. Finally I will discuss future research directions worth pursuing in this research community.


title: Neutrino Splitting for Lorentz-Violating Neutrinos: Detailed Analysis

abstract: Lorentz-violating neutrino parameters have been severely constrained on the basis of astrophysical considerations. In the high-energy limit, one generally assumes a superluminal dispersion relation of an incoming neutrino of the form E ~ |p|v, where E is the energy, p is the momentum and $v = sqrt(1 + delta) > 1. Lepton-pair creation due to a Cerenkov-radiation-like process (nu -> nu + e^- + e^+) becomes possible above a certain energy threshold, and bounds on the Lorentz-violating parameter delta can be derived. Here, we investigate a related process, nu_i -> nu_i + nu_f + bar_nu_f, where nu_i is an incoming neutrino mass eigenstate, while nu_f is the final neutrino mass eigenstate, with a superluminal velocity that is slightly slower than that of the initial state. This process is kinematically allowed if the Lorentz-violating parameters at high energy differ for the different neutrino mass eigenstates. Neutrino splitting is not subject to any significant energy threshold condition and could yield quite a substantial contribution to decay and energy loss processes at high energy, even if the differential Lorentz violation among neutrino flavors is severely constrained by other experiments. We also discuss the SU(2)-gauge invariance of the superluminal models and briefly discuss the use of a generalized vierbein formalism in the formulation of the Lorentz-violating Dirac equation.


title: Detailed Human Shape Estimation from a Single Image by Hierarchical Mesh Deformation

abstract: This paper presents a novel framework to recover detailed human body shapes from a single image. It is a challenging task due to factors such as variations in human shapes, body poses, and viewpoints. Prior methods typically attempt to recover the human body shape using a parametric based template that lacks the surface details. As such the resulting body shape appears to be without clothing. In this paper, we propose a novel learning-based framework that combines the robustness of parametric model with the flexibility of free-form 3D deformation. We use the deep neural networks to refine the 3D shape in a Hierarchical Mesh Deformation (HMD) framework, utilizing the constraints from body joints, silhouettes, and per-pixel shading information. We are able to restore detailed human body shapes beyond skinned models. Experiments demonstrate that our method has outperformed previous state-of-the-art approaches, achieving better accuracy in terms of both 2D IoU number and 3D metric distance. The code is available in https://github.com/zhuhao-nju/hmd.git


title: Fekete's lemma for componentwise subadditive functions of two or more real variables

abstract: We prove an analogue of Fekete's subadditivity lemma for functions of several real variables which are subadditive in each variable taken singularly. This extends both the classical case for subadditive functions of one real variable, and a result in a previous paper by the author. While doing so, we prove that the functions with the property mentioned above are bounded in every closed and bounded subset of their domain. The arguments follows those of Chapter 6 in E. Hille's 1948 textbook.


title: Quantum-Inspired Computing: Can it be a Microscopic Computing Model of the Brain?

abstract: Quantum computing and the workings of the brain have many aspects in common and have been attracting increasing attention in academia and industry. The computation in both is parallel and non-discrete. Though the underlying physical dynamics (e.g., equation of motion) may be deterministic, the observed or interpreted outcomes are often probabilistic. Consequently, various investigations have been undertaken to understand and reproduce the brain on the basis of quantum physics and computing. However, there have been arguments on whether the brain can and have to take advantage of quantum phenomena that need to survive in the macroscopic space-time region at room temperature. This paper presents a unique microscopic computational model for the brain based on an ansatz that the brain computes in a manner similar to quantum computing, but with classical waves. Log-scale encoding of information in the context of computing with waves is shown to play a critical role in bridging the computing models with classical and quantum waves. Our quantum-inspired computing model opens up a possibility of unifying the computing framework of artificial intelligence and quantum computing beyond quantum machine learning approaches.


title: Generating Long Sequences with Sparse Transformers

abstract: Transformers are powerful sequence models, but require time and memory that grows quadratically with the sequence length. In this paper we introduce sparse factorizations of the attention matrix which reduce this to $O(n \sqrt{n})$. We also introduce a) a variation on architecture and initialization to train deeper networks, b) the recomputation of attention matrices to save memory, and c) fast attention kernels for training. We call networks with these changes Sparse Transformers, and show they can model sequences tens of thousands of timesteps long using hundreds of layers. We use the same architecture to model images, audio, and text from raw bytes, setting a new state of the art for density modeling of Enwik8, CIFAR-10, and ImageNet-64. We generate unconditional samples that demonstrate global coherence and great diversity, and show it is possible in principle to use self-attention to model sequences of length one million or more.


title: Baryogenesis via Gravitational Spontaneous Symmetry Breaking

abstract: We study baryogenesis in effective field theories where a $\mathrm{U}(1)_{ B-L}$-charged scalar couples to gravity via curvature invariants. We analyze the general possibilities in such models, noting the relationships between some of them and existing models, such as Affleck-Dine baryogenesis. We then identify a novel mechanism in which $\mathrm{U}(1)_{ B-L}$ can be broken by couplings to the Gauss-Bonnet invariant during inflation and reheating. Using analytic methods, we demonstrate that this mechanism provides a new way to dynamically generate the net matter-anti-matter asymmetry observed today, and verify this numerically.


title: Large Scale Structure in CHILES

abstract: We demonstrate that the Discrete Persistent Source Extractor (DisPerSE) can be used with spectroscopic redshifts to define the cosmic web and its distance to galaxies in small area deepfields. Here we analyze the use of DisPerSE to identify structure in observational data. We apply DisPerSE to the distribution of galaxies in the COSMOS field and find the best parameters to identify filaments. We compile a catalog of 11500 spectroscopic redshifts from the Galaxy and Mass Assembly (GAMA) G10 data release. We analyze two-dimensional slices, extract filaments and calculate the distance for each galaxy to its nearest filament. We find that redder and more massive galaxies are closer to filaments. To study the growth of galaxies across cosmic time, and environment, we are carrying out an HI survey covering redshifts z = 0 - 0.45, the COSMOS HI Large Extragalactic Survey (CHILES). In addition we present the predicted HI mass fraction as a function of distance to filaments for the spectroscopically known galaxies in CHILES. Lastly, we discuss the cold gas morphology of a few individual galaxies and their positions with respect to the cosmic web. The identification of the cosmic web, and the ability of CHILES to study the resolved neutral hydrogen morphologies and kinematics of galaxies, will allow future studies of the properties of neutral hydrogen in different cosmic web environments across the redshift range z = 0.1 - 0.45.


title: Complex time, shredded propagator method for large-scale GW calculations

abstract: The GW method is a many-body electronic structure technique capable of generating accurate quasiparticle properties for realistic systems spanning physics, chemistry, and materials science. Despite its power, GW is not routinely applied to large complex assemblies due to its large computational overhead and quartic scaling with particle number. Here, the GW equations are recast, exactly, as Fourier-Laplace time integrals over complex time propagators. The propagators are then "shredded" via energy partitioning and the time integrals approximated in a controlled manner using generalized Gaussian quadrature(s) while discrete variable methods are employed to represent the required propagators in real-space. The resulting cubic scaling GW method has a sufficiently small prefactor to outperform standard quartic scaling methods on small systems ($\gtrapprox$ 10 atoms) and also represents a substantial improvement over other cubic methods tested for all system sizes studied. The approach can be applied to any theoretical framework containing large sums of terms with energy differences in the denominator.


title: Correlations between supermassive black holes, hot atmospheres, and the total masses of early type galaxies

abstract: We present a study of relations between the masses of the central supermassive black holes (SMBHs) and the atmospheric gas temperatures and luminosities measured within a range of radii between $R_{\rm e}$ and 5$R_{\rm e}$, for a sample of 47 early-type galaxies observed by the {\it Chandra X-ray Observatory}. We report the discovery of a tight correlation between the atmospheric temperatures of the brightest cluster/group galaxies (BCGs) and their central SMBH masses. Furthermore, our hydrostatic analysis reveals an approximately linear correlation between the total masses of BCGs ($M_{\rm tot}$) and their central SMBH masses ($M_{\rm BH}$). State-of-the-art cosmological simulations show that the SMBH mass could be determined by the binding energy of the halo through radiative feedback during the rapid black hole growth by accretion, while for the most massive galaxies mergers are the chief channel of growth. In the scenario of a simultaneous growth of central SMBHs and their host galaxies through mergers, the observed linear correlation could be a natural consequence of the central limit theorem.


title: A Fast and Accurate Algorithm for Spherical Harmonic Analysis on HEALPix Grids with Applications to the Cosmic Microwave Background Radiation

abstract: The Hierarchical Equal Area isoLatitude Pixelation (HEALPix) scheme is used extensively in astrophysics for data collection and analysis on the sphere. The scheme was originally designed for studying the Cosmic Microwave Background (CMB) radiation, which represents the first light to travel during the early stages of the universe's development and gives the strongest evidence for the Big Bang theory to date. Refined analysis of the CMB angular power spectrum can lead to revolutionary developments in understanding the nature of dark matter and dark energy. In this paper, we present a new method for performing spherical harmonic analysis for HEALPix data, which is a central component to computing and analyzing the angular power spectrum of the massive CMB data sets. The method uses a novel combination of a non-uniform fast Fourier transform, the double Fourier sphere method, and Slevinsky's fast spherical harmonic transform (Slevinsky, 2019). For a HEALPix grid with $N$ pixels (points), the computational complexity of the method is $\mathcal{O}(N\log^2 N)$, with an initial set-up cost of $\mathcal{O}(N^{3/2}\log N)$. This compares favorably with $\mathcal{O}(N^{3/2})$ runtime complexity of the current methods available in the HEALPix software when multiple maps need to be analyzed at the same time. Using numerical experiments, we demonstrate that the new method also appears to provide better accuracy over the entire angular power spectrum of synthetic data when compared to the current methods, with a convergence rate at least two times higher.

categories: math.NA astro-ph.CO cs.NA physics.comp-ph


title: Pointwise strong (H, Phi) approximation by Fourier series of L^{Psi} integrable functions

abstract: We essentially extend and improve the classical result of G. H. Hardy and J. E. Littlewood on strong summability of Fourier series. We will present an estimation of the generalized strong mean (H; Phi) as an approximation version of the Totik type generalization of the result of G. H. Hardy, J. E. Littlewood, in case of integrable functions from L^{Psi}. As a measure of such approximation we will use the function constructed by function Psi complementary to Phi on the base of defnition of the L^{Psi} points. Some corollary and remarks will also be given.


title: Abraham-Rubin-Shelah Open Colorings and a Large Continuum

abstract: We show that the Abraham-Rubin-Shelah Open Coloring Axiom is consistent with a large continuum, in particular, consistent with $2^{\aleph_0}=\aleph_3$. This answers one of the main open questions from the 1985 paper of Abraham-Rubin-Shelah. As in their paper, we need to construct names for so-called preassignments of colors in order to add the necessary homogeneous sets. However, these names are constructed over models satisfying the CH. In order to address this difficulty, we show how to construct such names with very strong symmetry conditions. This symmetry allows us to combine them in many different ways, using a new type of poset called a Partition Product, and thereby obtain a model of this axiom in which $2^{\aleph_0}=\aleph_3$.


title: Current-induced switching of YIG/Pt bilayers with in-plane magnetization due to Oersted fields

abstract: We report on the switching of the in-plane magnetization of thin yttrium iron garnet (YIG)/Pt bilayers induced by an electrical current. The switching is either field-induced and assisted by a dc current, or current-induced and assisted by a static magnetic field. The reversal of the magnetization occurs at a current density as low as $10^5$~A/cm$^{2}$ and magnetic fields of $\sim 40$~$\mu$T, two orders of magnitude smaller than in ferromagnetic metals, consistently with the weak uniaxial anisotropy of the YIG layers. We use the transverse component of the spin Hall magnetoresistance to sense the magnetic orientation of YIG while sweeping the current. Our measurements and simulations reveal that the current-induced effective field responsible for switching is due to the Oersted field generated by the current flowing in the Pt layer rather than by spin-orbit torques, and that the switching efficiency is influenced by pinning of the magnetic domains.


title: Flag-transitive block designs with prime replication number and almost simple groups

abstract: In this article, we study $2$-designs with prime replication number admitting a flag-transitive automorphism group. The automorphism groups of these designs are point-primitive of almost simple or affine type. We determine $2$-designs with prime replication number admitting an almost simple automorphism group.


title: Big images of two-dimensional pseudorepresentations

abstract: Bella\"iche has recently applied Pink-Lie theory to prove that, under mild conditions, the image of a continuous 2-dimensional pseudorepresentation $\rho$ of a profinite group on a local pro-$p$ domain $A$ contains a nontrivial congruence subgroup of ${\rm SL}_2(B)$ for a certain subring $B$ of $A$. We enlarge Bella\"iche's ring and give this new $B$ a conceptual interpretation in terms of conjugate self-twists of $\rho$, symmetries that naturally constrain its image. As a corollary, this new $B$ is optimal among congruence subgroups contained in the image. We also interpret the new $B$ vis-a-vis the adjoint trace ring of $\rho$, which we show is a more natural ring for these questions in general. Finally, we use our purely algebraic result to recover and extend a variety of arithmetic big-image results for ${\rm GL}_2$ Galois representations arising from elliptic, Hilbert, and Bianchi modular forms and $p$-adic Hida or Coleman families of elliptic and Hilbert modular forms.


title: Measurement of the mass difference and the binding energy of the hypertriton and antihypertriton

abstract: According to the CPT theorem, which states that the combined operation of charge conjugation, parity transformation and time reversal must be conserved, particles and their antiparticles should have the same mass and lifetime but opposite charge and magnetic moment. Here, we test CPT symmetry in a nucleus containing a strange quark, more specifically in the hypertriton. This hypernucleus is the lightest one yet discovered and consists of a proton, a neutron, and a $\Lambda$ hyperon. With data recorded by the STAR detector{\cite{TPC,HFT,TOF}} at the Relativistic Heavy Ion Collider, we measure the $\Lambda$ hyperon binding energy $B_{\Lambda}$ for the hypertriton, and find that it differs from the widely used value{\cite{B_1973}} and from predictions{\cite{2019_weak, 1995_weak, 2002_weak, 2014_weak}}, where the hypertriton is treated as a weakly bound system. Our results place stringent constraints on the hyperon-nucleon interaction{\cite{Hammer2002, STAR-antiH3L}}, and have implications for understanding neutron star interiors, where strange matter may be present{\cite{Chatterjee2016}}. A precise comparison of the masses of the hypertriton and the antihypertriton allows us to test CPT symmetry in a nucleus with strangeness for the first time, and we observe no deviation from the expected exact symmetry.


title: Accounting for Correlations When Fitting Extra Cosmological Parameters

abstract: Current cosmological tensions motivate investigating extensions to the standard $\Lambda$CDM model. Additional model parameters are typically varied one or two at a time, in a series of separate tests. The purpose of this paper is to highlight that information is lost by not also examining the correlations between these additional parameters, which arise when their effects on model predictions are similar, even if the parameters are not varied simultaneously. We show how these correlations can be quantified with simulations and Markov Chain Monte Carlo (MCMC) methods. As an example, we assume that $\Lambda$CDM is the true underlying model, and calculate the correlations expected between the phenomenological lensing amplitude parameter, $A_L$, the running of the spectral index, $n_{\rm run}$, and the primordial helium mass fraction, $Y_P$, when these parameters are varied one at a time along with the $\Lambda$CDM parameters in fits to the $\textit{Planck}$ 2015 temperature power spectrum. These correlations are not small, ranging from 0.31 ($A_L-n_{\rm run}$) to $-0.93$ ($n_{\rm run}-Y_P$). We find that the values of these three parameters from the $\textit{Planck}$ data are consistent with $\Lambda$CDM expectations within $0.9\sigma$ when the correlations are accounted for. This does not explain the 1.8-2.7$\sigma$ $\textit{Planck}$ preference for $A_L>1$, but provides an additional $\Lambda$CDM consistency test. For example, if $A_L>1$ was a symptom of an underlying systematic error or some real but unknown physical effect that also produced spurious correlations with $n_{\rm run}$ or $Y_P$ our test might have revealed this. We recommend that future cosmological analyses examine correlations between additional model parameters in addition to investigating them separately, one a time.


title: Block-distributed Gradient Boosted Trees

abstract: The Gradient Boosted Tree (GBT) algorithm is one of the most popular machine learning algorithms used in production, for tasks that include Click-Through Rate (CTR) prediction and learning-to-rank. To deal with the massive datasets available today, many distributed GBT methods have been proposed. However, they all assume a row-distributed dataset, addressing scalability only with respect to the number of data points and not the number of features, and increasing communication cost for high-dimensional data. In order to allow for scalability across both the data point and feature dimensions, and reduce communication cost, we propose block-distributed GBTs. We achieve communication efficiency by making full use of the data sparsity and adapting the Quickscorer algorithm to the block-distributed setting. We evaluate our approach using datasets with millions of features, and demonstrate that we are able to achieve multiple orders of magnitude reduction in communication cost for sparse data, with no loss in accuracy, while providing a more scalable design. As a result, we are able to reduce the training time for high-dimensional data, and allow more cost-effective scale-out without the need for expensive network communication.


title: A neural network-based framework for financial model calibration

abstract: A data-driven approach called CaNN (Calibration Neural Network) is proposed to calibrate financial asset price models using an Artificial Neural Network (ANN). Determining optimal values of the model parameters is formulated as training hidden neurons within a machine learning framework, based on available financial option prices. The framework consists of two parts: a forward pass in which we train the weights of the ANN off-line, valuing options under many different asset model parameter settings; and a backward pass, in which we evaluate the trained ANN-solver on-line, aiming to find the weights of the neurons in the input layer. The rapid on-line learning of implied volatility by ANNs, in combination with the use of an adapted parallel global optimization method, tackles the computation bottleneck and provides a fast and reliable technique for calibrating model parameters while avoiding, as much as possible, getting stuck in local minima. Numerical experiments confirm that this machine-learning framework can be employed to calibrate parameters of high-dimensional stochastic volatility models efficiently and accurately.


title: Defective hierarchical porous copper-based metal-organic frameworks synthesised via facile acid etching strategy

abstract: Introducing hierarchical pore structure to microporous materials such as metal-organic frameworks (MOFs) can be beneficial for reactions where the rate of reaction is limited by low rates of diffusion or high pressure drop. This advantageous pore structure can be obtained by defect formation, mostly via post-synthetic acid etching, which has been studied extensively on water-stable MOFs. Here we show that a water-unstable HKUST-1 MOF can also be modified in a corresponding manner by using phosphoric acid as a size-selective etching agent and a mixture of dimethyl sulfoxide and methanol as a dilute solvent. Interestingly, we demonstrate that the etching process which is time- and acidity- dependent, can result in formation of defective HKUST-1 with extra interconnected hexagonal macropores without compromising on the bulk crystallinity. These findings suggest an intelligent scalable synthetic method for formation of hierarchical porosity in MOFs that are prone to hydrolysis, for improved molecular accessibility and diffusion for catalysis.


title: A Rapid Fault Reconstruction Strategy Using a Bank of Sliding Mode Observers

abstract: This paper deals with the design of a model-based rapid fault detection and isolation strategy using sliding mode observers. To address this problem, a new scheme is proposed by adaptively combining the information provided by a bank of observers. In this regard, a new structure for sliding mode observers is considered. Then, the well-known recursive least square algorithm is utilized to merge individual state estimations suitably such that the system fault is detected faster. The required condition for enhancing perfect state estimation is derived, and the stability of the overall system is proven via Lyapunov's direct method. The supremacy of proposed scheme is fully discussed through mathematical analyses as well as simulations.


title: Upper Limits on Very-High-Energy Gamma-ray Emission from Core-Collapse Supernovae Observed with H.E.S.S

abstract: Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 hours to 53 hours. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 hours. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the $>1$ TeV gamma-ray flux of the order of $\sim$10$^{-13}$ cm$^{-2}$s$^{-1}$ are established, corresponding to upper limits on the luminosities in the range $\sim$2 $\times$ 10$^{39}$ erg s$^{-1}$ to $\sim$1 $\times$ 10$^{42}$ erg s$^{-1}$. These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between $\sim$2 $\times 10^{-5}$ and $\sim$2 $\times 10^{-3}$M$_{\odot}$yr$^{-1}$ under reasonable assumptions on the particle acceleration parameters.


title: Deconstructing the Filter Bubble: User Decision-Making and Recommender Systems

abstract: We study a model of user decision-making in the context of recommender systems via numerical simulation. Our model provides an explanation for the findings of Nguyen, et. al (2014), where, in environments where recommender systems are typically deployed, users consume increasingly similar items over time even without recommendation. We find that recommendation alleviates these natural filter-bubble effects, but that it also leads to an increase in homogeneity across users, resulting in a trade-off between homogenizing across-user consumption and diversifying within-user consumption. Finally, we discuss how our model highlights the importance of collecting data on user beliefs and their evolution over time both to design better recommendations and to further understand their impact.


title: Putting Mind Back into Nature: A Tribute to Henry P. Stapp

abstract: Henry Stapp has for 60 years been a leader, perhaps the leader, in exploring the role of mind (psyche, consciousness, experience) in the ontology of quantum mechanics. Henry's contention is that the very structure of quantum mechanics implies a central and irreducible role for mind: an experiential aspect of nature distinct from that of the physical matter and energy described by the dynamical equations of physics. The task then becomes generating interest across disciplines in the exploration of this thesis, and in particular to seek connections with neuroscience and with empirical psychology.


title: Approximate expression for the ground-state energy of the two- and three-dimensional Hubbard model at arbitrary filling obtained from dimensional scaling

abstract: We generalize the linear discrete dimensional scaling approach for the repulsive Hubbard model to obtain a nonlinear scaling relation that yields accurate approximations to the ground-state energy in both two and three dimensions, as judged by comparison to auxiliary-field quantum Monte Carlo (QMC) data. Predictions are made for the per-site ground-state energies in two and three dimensions for $n$ (filling factor) and $U$ (Coulomb interaction) values for which QMC data are currently unavailable.


title: Six Years of Sustained Activity from Active Asteroid (6478) Gault

abstract: We present archival observations demonstrating that main belt asteroid (6478) Gault has an extensive history of comet-like activity. Outbursts have taken place during multiple epochs since 2013 and at distances extending as far as 2.68 au, nearly aphelion. (6478) Gault is a member of the predominately S-type (i.e., volatile-poor) Phocaea family; no other main belt object of this type has ever shown more than a single activity outburst. Furthermore, our data suggest that this is the longest duration of activity caused by a body spinning near the rotational breakup barrier. If activity is indeed unrelated to volatiles, as appears to be the case, (6478) Gault represents a new class of object, perpetually active due to rotational spin-up.


title: Anisotropic Moser-Trudinger inequality involving $L^n$ norm

abstract: The paper is concerned about a sharp form of Anisotropic Moser-Trudinger inequality which involves $L^{n}$ norm. Let \begin{equation*} \lambda_{1}(\Omega) = \inf_{u\in W_0^{1,n}(\Omega),u\not\equiv 0} ||F(\nabla u)||_{L^n(\Omega)}^n / ||u||_{L^n(\Omega)}^n \end{equation*} be the first eigenvalue associated with $n$-Finsler-Laplacian. using blowing up analysis, we obtain that \begin{equation*} \sup_{u\in W_{0}^{1,n}(\Omega),||F(\nabla u)||_{L^n(\Omega)} = 1} \int_{\Omega}e^{\lambda_n (1+\alpha||u||_{L^n (\Omega)}^n)^{\frac{1}{n-1}} |u|^{\frac{n}{n-1}}}dx \end{equation*} is finite for any $0\leq \alpha<\lambda_{1}(\Omega)$,and the supremum is infinite for any $\alpha\geq \lambda_{1}(\Omega)$, where $\lambda_{n} = n^{\frac{n}{n-1}} \kappa_n^{\frac{1}{n-1}}$ ($\kappa_{n}$ is the volume of the unit wulff ball) and the function $F$ is positive,convex and homogeneous of degree $1$, and its polar $F^o$ represents a Finsler metric on $\mathbb{R}^n$. Furthermore, the supremum is attained for any $0\leq \alpha<\lambda_{1}(\Omega)$.


title: The Moore-Penrose inverses of split quaternions

abstract: In this paper, we find the roots of lightlike quaternions. By introducing the concept of the Moore-Penrose inverse in split quaternions, we solve the linear equations $axb=d$, $xa=bx$ and $xa=b\bar{x}$. Also we obtain necessary and sufficient conditions for two split quaternions to be similar or consimilar.


title: Estimating the size of the scatterer

abstract: Formula for the size of the scatterer is derived explicitly in terms of the scattering amplitude corresponding to this scatterer. By the scatterer either a bounded obstacle $D$ or the support of the compactly supported potential is meant


title: Global existence and uniqueness of the solution to a nonlinear parabolic equation

abstract: Consider the equation $$ u'(t)-\Delta u+|u|^\rho u=0, \quad u(0)=u_0(x), (1), $$ where $ u':=\frac {du}{dt}$, $ \rho=const >0, $ $x\in \mathbb{R}^3$, $t>0$. Assume that $u_0$ is a smooth and decaying function, $$\|u_0\|\:=\sup_{x\in \mathbb{R}^3, t\in \mathbb{R}_+} |u(x,t)|.$$ It is proved that problem (1) has a unique global solution and this solution satisfies the following estimate $$\|u(x,t)\|<c, $$ where $c>0$ does not depend on $x,t$.


title: Evaluation of MRI to ultrasound registration methods for brain shift correction: The CuRIOUS2018 Challenge

abstract: In brain tumor surgery, the quality and safety of the procedure can be impacted by intra-operative tissue deformation, called brain shift. Brain shift can move the surgical targets and other vital structures such as blood vessels, thus invalidating the pre-surgical plan. Intra-operative ultrasound (iUS) is a convenient and cost-effective imaging tool to track brain shift and tumor resection. Accurate image registration techniques that update pre-surgical MRI based on iUS are crucial but challenging. The MICCAI Challenge 2018 for Correction of Brain shift with Intra-Operative UltraSound (CuRIOUS2018) provided a public platform to benchmark MRI-iUS registration algorithms on newly released clinical datasets. In this work, we present the data, setup, evaluation, and results of CuRIOUS 2018, which received 6 fully automated algorithms from leading academic and industrial research groups. All algorithms were first trained with the public RESECT database, and then ranked based on test dataset of 10 additional cases with identical data curation and annotation protocols as the RESECT database. The article compares the results of all participating teams and discusses the insights gained from the challenge, as well as future work.


title: Frequency measurement of the $^{1}\mathrm{S}_{0}, F=5/2\leftrightarrow\,^{3}\mathrm{P}_{1}, F=7/2$ transition of $^{27}$Al$^{+}$ via quantum logic spectroscopy with $^{40}$Ca$^{+}$

abstract: We perform quantum logic spectroscopy with a $^{27}$Al$^{+}$/$^{40}$Ca$^{+}$ mixed ion crystal in a linear Paul trap for a measurement of the $(3s^{2})\,^{1}\mathrm{S}_{0} \leftrightarrow \, (3s3p)\,^{3}\mathrm{P}_{1}, F=7/2$ intercombination transition in $^{27}$Al$^{+}$. Towards this end, Ramsey spectroscopy is used for probing the transition in $^{27}$Al$^{+}$ and the $(4s^{2})\,\mathrm{S}_{1/2} \leftrightarrow \, (4s3d)\,\mathrm{D}_{5/2}$ clock transition in $^{40}$Ca$^{+}$ in interleaved measurements. By using the precisely measured frequency of the clock transition in $^{40}$Ca$^{+}$ as a frequency reference, we determine the frequency of the intercombination line to be $\nu_{^{1}\mathrm{S}_{0} \leftrightarrow \,^{3}\mathrm{P}_{1},F=7/2}$=1122 842 857 334 736(93) Hz and the Land\'e g-factor of the excited state to be $g_{^{3}\mathrm{P}_{1}, F=7/2}$=0.428132(2).


title: An Instrument-Free Demonstration of Quantum Key Distribution for High-School Students

abstract: It has become increasingly common for high-school students to see media reports on the importance of quantum mechanics in the development of next-generation industries such as drug development and secure communication, but few of them have been exposed to fundamental quantum mechanical concepts in a meaningful classroom activity. In order to bridge this gap, we design and test a low-cost 20-minute demonstration of the Bell test, which is used in several entanglement-based quantum key distribution protocols. The demonstration introduces ideas such as the quantum state, quantum measurement, spin quantization, cryptography, and entanglement; all without using concepts beyond the 9th grade of the Chilean high-school curriculum. The demonstration can serve to promote early exposure of the future adopters and developers of quantum technology with its conceptual building blocks, and also to educate the general public about the importance of quantum mechanics in modern industry


title: Extra invariance of principal shift invariant spaces and the Zak transform

abstract: We prove a necessary and sufficient condition for a principal shift invariant space of $L^2(\mathbb{R})$ to be invariant under translations by the subgroup $\frac{1}{N} \mathbb{Z}, N>1$. This condition is given in terms of the Zak transform of the group $\frac{1}{N} \mathbb{Z}.$ This result is extended to principal shift invariant spaces generated by a lattice in a general locally compact abelian (LCA) group.


title: Diversity in density profiles of self-interacting dark matter satellite halos

abstract: We present results from N-body simulations of self-interacting dark matter (SIDM) subhalos, which could host ultra-faint dwarf spheroidal galaxies, inside a Milky-Way-like main halo. We find that high-concentration subhalos are driven to gravothermal core collapse, while low-concentration subhalos develop large (kpc-sized) low-density cores, with both effects depending sensitively on the satellite's orbit and the self-interaction cross section over mass $\sigma/m$. The overall effect for $\sigma/m \gtrsim 3 \ \rm cm^2/g$ is to increase the range of inner densities, potentially explaining the observed diversity of Milky Way satellites, which include compact systems like Draco and Segue 1 that are dense in dark matter, and less dense, diffuse systems like Sextans and Crater II. We discuss possible ways of distinguishing SIDM models from collisionless dark matter models using the inferred dark matter densities and stellar sizes of the dwarf spheroidal galaxies.


title: Some Problems in Logic: Applications of Kripke's Notion of Fulfilment

abstract: This is a study of S. Kripke's notion of fulfilment. Motivated by Paris-Harrington statement, Kripke was looking for a proof of G\"odel's Incompleteness Theorem which was model-theoretic, natural (without self-reference), and easy. Fulfilment gives a versatile tool for both Proof and Model Theory. We begin with short proofs to a number of classical results. With two new results: there is an easily definable subring $R$ of the primitive recursive functions such that for any non-principal ultrafilter $D$ on $\omega$, $R/D$ is a recursively saturated model of Peano arithmetic; and for any r.e. theory $\textsf{T}$ and for any given r.e. set, we can feasibly find a $\Sigma_1^0$ formula which semi-represents it in $\textsf{T}$. We then give a version of Herbrand's Theorem, and of the Hilbert-Ackermann method of proving consistency, answering a problem of D. Guaspari:\[\{\ulcorner\phi\urcorner\in\Pi^0_k:\phi\text{ is }\Sigma^0_k\text{-conservative over $\textsf{PA}$}\}\]is a complete $\Pi^0_2$ set. We extend H. Friedman's method for results of such as $\Sigma^1_2\textsf{-AC}$ is $\Pi^1_3$-conservative over $(\Pi^1_1\textsf{-CA})_{<\varepsilon_0}\upharpoonright$, along with uniform versions\[\forall\alpha<\varepsilon_0\;(\Pi^1_1\textsf{-CA})_{\alpha}\upharpoonright\,\vdash\textsf{RFN}_{\Pi^1_3}({\Sigma^1_2}\textsf{-AC})\] Then there are some model-theoretic applications, starting with non-$\omega$-models. We extend the theorem of D. Scott involving Weak K\"onig's Lemma, and describe the order types of elementary initial segments of recursively saturated models. For $\omega$-models, we extend Friedman's theorem on minimal models of analysis, and develop indicators for countable fragments of $\mathcal L_{\infty\omega}$, with some representability results in $\omega$-logic. We close with an exposition of the Paris-Harrington statement.


title: Fixed-Depth Two-Qubit Circuits and the Monodromy Polytope

abstract: For a native gate set which includes all single-qubit gates, we apply results from symplectic geometry to analyze the spaces of two-qubit programs accessible within a fixed number of gates. These techniques yield an explicit description of this subspace as a convex polytope, presented by a family of linear inequalities themselves accessible via a finite calculation. We completely describe this family of inequalities in a variety of familiar example cases, and as a consequence we highlight a certain member of the "XY-family" for which this subspace is particularly large, i.e., for which many two-qubit programs admit expression as low-depth circuits.


title: Topological bound states for quantum charges

abstract: We discuss how, in appropriately designed configurations, solenoids carrying a semifluxon can be used as topological energy barriers for charged quantum systems. We interpret this phenomenon as a consequence of the fact that such solenoids induce nodal lines in the wave function describing the charge, which on itself is a consequence of the Aharonov-Bohm effect. Moreover, we present a thought experiment with a cavity where two solenoids are sufficient to create bound states.


title: Comment on "Search for an axion-induced oscillating electric dipole moment for electrons using atomic magnetometers"

abstract: In the recent work arXiv:1809.02446, the authors proposed a new method measuring the electron oscillating electric dipole moment (eOEDM) using atomic magnetomaters. This eOEDM is induced by the interaction between the electron magnetic dipole moment, electric field and axion field. The result is sensitive to the axion-photon coupling according to [Hill, PRD 91, 111702 (2015)]. Here we want to describe that the same experimental method can be also sensitive to the axion-electron coupling according to [Alexander and Sims, PRD 98, 015011 (2018)]. In this article, we will show the corresponding sensitivity plot and compare with other constraints.


title: Developments in Cosmic Growth and Gravitation

abstract: Cosmic surveys of large scale structure have imaged hundreds of millions of galaxies and mapped the 3D positions of over a million. Surveys starting over the next few years will increase these numbers more than tenfold. Simultaneously, developments in extracting information on dark energy, dark matter, neutrinos, and gravity on cosmic scales have advanced greatly, with many important works from Asia Pacific institutions.


title: Leptogenesis in $\Delta(27)$ with a Universal Texture Zero

abstract: We investigate the possibility of viable leptogenesis in an appealing $\Delta(27)$ model with a universal texture zero in the (1,1) entry. The model accommodates the mass spectrum, mixing and CP phases for both quarks and leptons and allows for grand unification. Flavoured Boltzmann equations for the lepton asymmetries are solved numerically, taking into account both $N_1$ and $N_2$ right-handed neutrino decays. The $N_1$-dominated scenario is successful and the most natural option for the model, with $M_1 \in [10^9, 10^{12}]$ GeV, and $M_1/M_2 \in [0.002, 0.1]$, which constrains the parameter space of the underlying model and yields lower bounds on the respective Yukawa couplings. Viable leptogenesis is also possible in the $N_2$-dominated scenario, with the asymmetry in the electron flavour protected from $N_1$ washout by the texture zero. However, this occurs in a region of parameter space which has a stronger mass hierarchy $M_1/M_2 < 0.002 $, and $M_2$ relatively close to $M_3$, which is not a natural expectation of the $\Delta(27)$ model.


title: Embedded nonlinear model predictive control for obstacle avoidance using PANOC

abstract: We employ the proximal averaged Newton-type method for optimal control (PANOC) to solve obstacle avoidance problems in real time. We introduce a novel modeling framework for obstacle avoidance which allows us to easily account for generic, possibly nonconvex, obstacles involving polytopes, ellipsoids, semialgebraic sets and generic sets described by a set of nonlinear inequalities. PANOC is particularly well-suited for embedded applications as it involves simple steps, its implementation comes with a low memory footprint and its fast convergence meets the tight runtime requirements of fast dynamical systems one encounters in modern mechatronics and robotics. The proposed obstacle avoidance scheme is tested on a lab-scale autonomous vehicle.


title: Vortex dynamics under pulsatile flow in axisymmetric constricted tubes

abstract: An improved understanding of how vortices develop and propagate under pulsatile flow can shed important light on the mixing and transport processes including the transition to turbulent regime occurring in such systems. For example, the characterization of pulsatile flows in obstructed artery models serves to encourage research into flow-induced phenomena associated with changes in morphology, blood viscosity, wall elasticity and flow rate. In this work, an axisymmetric rigid model was used to study the behaviour of the flow pattern with varying constriction degree ($d_0$), mean Reynolds number ($\bar{Re}$) and Womersley number ($\alpha$). Velocity fields were acquired experimentally using Digital Particle Image Velocimetry and generated numerically. For the acquisition of data, $\bar{Re}$ was varied from 385 to 2044, $d_0$ was 1.0 cm and 1.6 cm, and $\alpha$ was varied from 17 to 33 in the experiments and from 24 to 50 in the numerical simulations. Results for the considered Reynolds number, showed that the flow pattern consisted of two main structures: a central jet around the tube axis and a recirculation zone adjacent to the inner wall of the tube, where vortices shed. Using the vorticity fields, the trajectory of vortices was tracked and their displacement over their lifetime calculated. The analysis led to a scaling law equation for the maximum vortex displacement as a function of a dimensionless variable dependent on the system parameters Re and $\alpha$.


title: Uncertainty-aware demand management of water distribution networks in deregulated energy markets

abstract: We present an open-source solution for the operational control of drinking water distribution networks which accounts for the inherent uncertainty in water demand and electricity prices in the day-ahead market of a volatile deregulated economy. As increasingly more energy markets adopt this trading scheme, the operation of drinking water networks requires uncertainty-aware control approaches that mitigate the effect of volatility and result in an economic and safe operation of the network that meets the consumers' need for uninterrupted water supply. We propose the use of scenario-based stochastic model predictive control: an advanced control methodology which comes at a considerable computation cost which is overcome by harnessing the parallelization capabilities of graphics processing units (GPUs) and using a massively parallelizable algorithm based on the accelerated proximal gradient method.


title: Chandra and Hubble Space Telescope observations of dark gamma-ray bursts and their host galaxies

abstract: We present a study of 21 dark gamma-ray burst (GRB) host galaxies, predominantly using X-ray afterglows obtained with the Chandra X-Ray Observatory (CXO) to precisely locate the burst in deep Hubble Space Telescope (HST) imaging of the burst region. The host galaxies are well-detected in F160W in all but one case and in F606W imaging in approx 60 per cent of cases. We measure magnitudes and perform a morphological analysis of each galaxy. The asymmetry, concentration and ellipticity of the dark burst hosts are compared against the host galaxies of optically bright GRBs. In agreement with other studies, we find that dark GRB hosts are redder and more luminous than the bulk of the GRB host population. The distribution of projected spatial offsets for dark GRBs from their host galaxy centroids is comparable to that of optically-bright bursts. The dark GRB hosts are physically larger, more massive and redder, but are morphologically similar to the hosts of bright GRBs in terms of concentration and asymmetry. Our analysis constrains the fraction of high redshift (z greater than 5) GRBs in the sample to approx 14 per cent, implying an upper limit for the whole long-GRB population of less than 4.4 per cent. If dust is the primary cause of afterglow darkening amongst dark GRBs, the measured extinction may require a clumpy dust component in order to explain the observed offset and ellipticity distributions.


title: Sparse Nerves in Practice

abstract: Topological data analysis combines machine learning with methods from algebraic topology. Persistent homology, a method to characterize topological features occurring in data at multiple scales is of particular interest. A major obstacle to the wide-spread use of persistent homology is its computational complexity. In order to be able to calculate persistent homology of large datasets, a number of approximations can be applied in order to reduce its complexity. We propose algorithms for calculation of approximate sparse nerves for classes of Dowker dissimilarities including all finite Dowker dissimilarities and Dowker dissimilarities whose homology is Cech persistent homology. All other sparsification methods and software packages that we are aware of calculate persistent homology with either an additive or a multiplicative interleaving. In dowker_homology, we allow for any non-decreasing interleaving function $\alpha$. We analyze the computational complexity of the algorithms and present some benchmarks. For Euclidean data in dimensions larger than three, the sizes of simplicial complexes we create are in general smaller than the ones created by SimBa. Especially when calculating persistent homology in higher homology dimensions, the differences can become substantial.


title: CascadeML: An Automatic Neural Network Architecture Evolution and Training Algorithm for Multi-label Classification

abstract: Multi-label classification is an approach which allows a datapoint to be labelled with more than one class at the same time. A common but trivial approach is to train individual binary classifiers per label, but the performance can be improved by considering associations within the labels. Like with any machine learning algorithm, hyperparameter tuning is important to train a good multi-label classifier model. The task of selecting the best hyperparameter settings for an algorithm is an optimisation problem. Very limited work has been done on automatic hyperparameter tuning and AutoML in the multi-label domain. This paper attempts to fill this gap by proposing a neural network algorithm, CascadeML, to train multi-label neural network based on cascade neural networks. This method requires minimal or no hyperparameter tuning and also considers pairwise label associations. The cascade algorithm grows the network architecture incrementally in a two phase process as it learns the weights using adaptive first order gradient algorithm, therefore omitting the requirement of preselecting the number of hidden layers, nodes and the learning rate. The method was tested on 10 multi-label datasets and compared with other multi-label classification algorithms. Results show that CascadeML performs very well without hyperparameter tuning.


title: ML-KFHE: Multi-label ensemble classification algorithm exploiting sensor fusion properties of the Kalman filter

abstract: Despite the success of ensemble classification methods in multi-class classification problems, ensemble methods based on approaches other than bagging have not been widely explored for multi-label classification problems. The Kalman Filter-based Heuristic Ensemble (KFHE) is an ensemble method that exploits the sensor fusion properties of the Kalman filter to combine several classifier models, and that has been shown to be very effective. This work proposes a multi-label version of KFHE, ML-KFHE, demonstrating the effectiveness of the KFHE method on multi-label datasets. Two variants are introduced based on the underlying component classifier algorithm, ML-KFHE-HOMER, and ML-KFHE-CC which uses HOMER and Classifier Chain (CC) as the underlying multi-label algorithms respectively. ML-KFHE-HOMER and ML-KFHE-CC sequentially train multiple HOMER and CC multi-label classifiers and aggregate their outputs using the sensor fusion properties of the Kalman filter. Extensive experiments and detailed analysis were performed on thirteen multi-label datasets and eight other algorithms, which included state-of-the-art ensemble methods. The results show, for both versions, the ML-KFHE framework improves the predictive performance significantly with respect to bagged combinations of HOMER (named E-HOMER), also introduced in this paper, and bagged combination of CC, Ensemble Classifier Chains (ECC), thus demonstrating the effectiveness of ML-KFHE. Also, the ML-KFHE-HOMER variant was found to perform consistently and significantly better than the compared multi-label methods including existing approaches based on ensembles.


title: Quantum Adiabatic Doping with Incommensurate Optical Lattices

abstract: Quantum simulations of Fermi-Hubbard models have been attracting considerable efforts in the optical lattice research, with the ultracold anti-ferromagnetic atomic phase reached at half filling in recent years. An unresolved issue is to dope the system while maintaining the low thermal entropy. Here we propose to achieve the low temperature phase of the doped Fermi-Hubbard model using incommensurate optical lattices through adiabatic quantum evolution. In this theoretical proposal, we find that one major problem about the adiabatic doping that shows up is atomic localization in the incommensurate lattice, potentially causing exponential slowing down of the adiabatic procedure. We study both one- and two-dimensional incommensurate optical lattices, and find that the localization prevents efficient adiabatic doping in the strong lattice regime for both cases. With density matrix renormalization group calculation, we further show that the slowing down problem in one dimension can be circumvented by considering interaction induced many-body delocalization, which is experimentally feasible using Feshbach resonance techniques. This protocol is expected to be efficient as well in two dimensions where the localization phenomenon is less stable.


title: Deep Q-Learning for Nash Equilibria: Nash-DQN

abstract: Model-free learning for multi-agent stochastic games is an active area of research. Existing reinforcement learning algorithms, however, are often restricted to zero-sum games, and are applicable only in small state-action spaces or other simplified settings. Here, we develop a new data efficient Deep-Q-learning methodology for model-free learning of Nash equilibria for general-sum stochastic games. The algorithm uses a local linear-quadratic expansion of the stochastic game, which leads to analytically solvable optimal actions. The expansion is parametrized by deep neural networks to give it sufficient flexibility to learn the environment without the need to experience all state-action pairs. We study symmetry properties of the algorithm stemming from label-invariant stochastic games and as a proof of concept, apply our algorithm to learning optimal trading strategies in competitive electronic markets.


title: Sur l'indice des sous-alg\`ebres biparaboliques de $\mathfrak{gl}(n)$

abstract: We give some properties on the index of seaweed subalgebras of the complex Lie algebra $\mathfrak{gl}(n)$ which allow to obtain formulas for the index of some interesting classes of this family and to give new families of Frobenius Lie algebras.


title: Fractional Calculus in Pharmacokinetics

abstract: We are witnessing the birth of a new variety of pharmacokinetics where non-integer-order differential equations are employed to study the time course of drugs in the body: this is dubbed "fractional pharmacokinetics." The presence of fractional kinetics has important clinical implications such as the lack of a half-life, observed, for example with the drug amiodarone and the associated irregular accumulation patterns following constant and multiple-dose administration. Building models that accurately reflect this behaviour is essential for the design of less toxic and more effective drug administration protocols and devices. This article introduces the readers to the theory of fractional pharmacokinetics and the research challenges that arise. After a short introduction to the concepts of fractional calculus, and the main applications that have appeared in literature up to date, we address two important aspects. First, numerical methods that allow us to simulate fractional order systems accurately and second, optimal control methodologies that can be used to design dosing regimens to individuals and populations.


title: A short proof of the discontinuity of phase transition in the planar random-cluster model with $q>4$

abstract: The goal of this paper is to provide a short proof of the discontinuity of phase transition for the random-cluster model on the square lattice with parameter $q>4$. This result was recently shown via the so-called Bethe ansatz for the six-vertex model. Our proof also exploits the connection to the six-vertex model, but does not rely on the Bethe ansatz. Our argument is soft and only uses very basic properties of the random-cluster model (for example, we do not need the Russo--Seymour--Welsh theory).


title: 1H(n,el) as a Cross Section Reference in a White Source Neutron Beam With the fissionTPC

abstract: We provide a quantitative description of a method to measure neutron-induced fission cross sections in ratio to elastic hydrogen scattering in a white-source neutron beam with the fission Time Projection Chamber. This detector has measured precision fission cross section ratios using actinide references such as $^{235}$U(n,f) and $^{238}$U(n,f). However, by employing a more precise reference such as the H(n,el) cross section there is the potential to further reduce the evaluation uncertainties of the measured cross sections. In principle the fissionTPC could provide a unique measurement by simultaneously measuring both fission fragments and proton recoils over a large solid angle. We investigate one method with a hydrogenous gas target and with the neutron energy determined by the proton recoil kinematics. This method enables the measurement to be performed in a white-source neutron beam and with the current configuration of the fissionTPC. We show that while such a measurement is feasible in the energy range of 0.5 MeV to $\sim$10 MeV, uncertainties on the proton detection efficiency and the neutron energy resolution do not allow us to preform a fission ratio measurement to the desired precision. Utilizing either a direct measurement of the neutron time-of-flight for the recoil proton or a mono-energetic neutron source or some combination of both would provide a path to a sub-percent precision measurement.


title: Neutrino Oscillations in a Quantum Processor

abstract: Quantum computing technologies promise to revolutionize calculations in many areas of physics, chemistry, and data science. Their power is expected to be especially pronounced for problems where direct analogs of a quantum system under study can be encoded coherently within a quantum computer. A first step toward harnessing this power is to express the building blocks of known physical systems within the language of quantum gates and circuits. In this paper, we present a quantum calculation of an archetypal quantum system: neutrino oscillations. We define gate arrangements that implement the neutral lepton mixing operation and neutrino time evolution in two-, three-, and four-flavor systems. We then calculate oscillation probabilities by coherently preparing quantum states within the processor, time evolving them unitarily, and performing measurements in the flavor basis, with close analogy to the physical processes realized in neutrino oscillation experiments, finding excellent agreement with classical calculations. We provide recipes for modeling oscillation in the standard three-flavor paradigm as well as beyond-standard-model scenarios, including systems with sterile neutrinos, non-standard interactions, Lorentz symmetry violation, and anomalous decoherence.


title: Cold, old and metal-poor: New stellar substructures in the Milky Way's dwarf spheroidals

abstract: Dwarf spheroidal galaxies (dSph) orbiting the Milky Way are complex objects often with complicated star formation histories and internal dynamics. In this work, we search for stellar substructures in four of the classical dSph satellites of the Milky Way: Sextans, Carina, Leo I, and Leo II. We apply two methods to search for stellar substructure: the minimum spanning tree method, which helps us to find and quantify spatially connected structures, and the "brute-force" method, which is able to find elongated stellar substructures. We detected the previously known substructure in Sextans, and also found a new stellar substructure within Sextans. Furthermore, we identified a new stellar substructure close to the core radius of the Carina dwarf galaxy. We report a detection of one substructure in Leo I and two in Leo II, but we note that we are dealing with a low number of stars in the samples used. Such old stellar substructures in dSph galaxies could help us to shed light on the nature of the dark matter halos, within which such structures form, evolve, and survive.


title: A $p$-adic Stark conjecture in the rank one setting

abstract: We give a new definition of a $p$-adic $L$-function for a mixed signature character of a real quadratic field and for a nontrivial ray class character of an imaginary quadratic field. We then state a $p$-adic Stark conjecture for this $p$-adic $L$-function. We prove our conjecture in the case when $p$ is split in the imaginary quadratic field by relating our construction to Katz's $p$-adic $L$-function. We also provide numerical evidence for our conjecture in three examples.


title: A Non-Linear Roth Theorem for Fractals of Sufficiently Large Dimension

abstract: Suppose that $d \geq 2$, and that $A \subset [0,1]$ has sufficiently large dimension, $1 - \epsilon_d < \dim_H(A) < 1$. Then for any polynomial $P$ of degree $d$ with no constant term, there exists a point configuration $\{ x, x-t,x-P(t) \} \subset A$ with $t \approx_P 1$.


title: Thermal photon drag in many-body systems

abstract: We demonstrate the existence of a thermal analog of Coulomb drag in many-body systems which is driven by thermal photons. We show that this frictional effect can either be positive or negative depending on the separation distances within the system. Also we highlight that the persistent heat currents flowing in non-reciprocal systems at equilibrium are subject to this effect and the latter can even amplify these flows.