title: Solar Flare Energy Partitioning and Transport -- the Gradual Phase (a Heliophysics 2050 White Paper)

abstract: Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare gradual phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, the flare gradual phase persists much longer than predicted so, by 2050, we must identify the characteristics of the significant energy deposition sustaining the gradual phase, and address the fundamental processes of turbulence and non-local heat flux.


title: Thermodynamics of Energy Magnetization

abstract: We construct the thermodynamics of energy magnetization in the presence of gravitomagnetic field. We show that the free energy must be modified to account for the modification of the energy current operator in the presence of a confining potential. The explicit expression of the energy magnetization is derived for a periodic system, and the Streda formula for the thermal Hall conductivity is rigorously established. We demonstrate our theory of the energy magnetization and the Streda formula in a Chern insulator.


title: Evaluating Low-Frequency Pulsar Observations to Monitor Dispersion with the Giant Metrewave Radio Telescope

abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project has the primary goal of detecting and characterizing low-frequency gravitational waves through high-precision pulsar timing. The mitigation of interstellar effects is crucial to achieve the necessary precision for gravitational wave detection. Effects like dispersion and scattering are more influential at lower observing frequencies, with the variation of these quantities over week-month timescales requiring high-cadence multi-frequency observations for pulsar timing projects. In this work, we utilize the dual-frequency observing capability of the Giant Metrewave Radio Telescope (GMRT) and evaluate the potential decrease in dispersion measure (DM) uncertainties when combined with existing pulsar timing array data. We present the timing analysis for four millisecond pulsars observed with the GMRT simultaneously at 322 and 607 MHz, and compare the DM measurements with those obtained through NANOGrav observations with the Green Bank Telescope (GBT) and Arecibo Observatory at 1400 to 2300 MHz frequencies. Measured DM values with the GMRT and NANOGrav program show significant offsets for some pulsars, which could be caused by pulse profile evolution in the two frequency bands. In comparison to the predicted DM uncertainties when incorporating these low-frequency data into the NANOGrav dataset, we find that higher-precision GMRT data is necessary to provide improved DM measurements. Through the detection and analysis of pulse profile baseline ripple in data on test pulsar B1929+10, we find that, while not important for this data, it may be relevant for other timing datasets. We discuss the possible advantages and challenges of incorporating GMRT data into NANOGrav and International Pulsar Timing Array datasets.


title: Population Mapping in Informal Settlements with High-Resolution Satellite Imagery and Equitable Ground-Truth

abstract: We propose a generalizable framework for the population estimation of dense, informal settlements in low-income urban areas--so called 'slums'--using high-resolution satellite imagery. Precise population estimates are a crucial factor for efficient resource allocations by government authorities and NGO's, for instance in medical emergencies. We utilize equitable ground-truth data, which is gathered in collaboration with local communities: Through training and community mapping, the local population contributes their unique domain knowledge, while also maintaining agency over their data. This practice allows us to avoid carrying forward potential biases into the modeling pipeline, which might arise from a less rigorous ground-truthing approach. We contextualize our approach in respect to the ongoing discussion within the machine learning community, aiming to make real-world machine learning applications more inclusive, fair and accountable. Because of the resource intensive ground-truth generation process, our training data is limited. We propose a gridded population estimation model, enabling flexible and customizable spatial resolutions. We test our pipeline on three experimental site in Nigeria, utilizing pre-trained and fine-tune vision networks to overcome data sparsity. Our findings highlight the difficulties of transferring common benchmark models to real-world tasks. We discuss this and propose steps forward.


title: Dynamic causal modelling of immune heterogeneity

abstract: An interesting inference drawn by some Covid-19 epidemiological models is that there exists a proportion of the population who are not susceptible to infection -- even at the start of the current pandemic. This paper introduces a model of the immune response to a virus. This is based upon the same sort of mean-field dynamics as used in epidemiology. However, in place of the location, clinical status, and other attributes of people in an epidemiological model, we consider the state of a virus, B and T-lymphocytes, and the antibodies they generate. Our aim is to formalise some key hypotheses as to the mechanism of resistance. We present a series of simple simulations illustrating changes to the dynamics of the immune response under these hypotheses. These include attenuated viral cell entry, pre-existing cross-reactive humoral (antibody-mediated) immunity, and enhanced T-cell dependent immunity. Finally, we illustrate the potential application of this sort of model by illustrating variational inversion (using simulated data) of this model to illustrate its use in testing hypotheses. In principle, this furnishes a fast and efficient immunological assay--based on sequential serology--that provides a (i) quantitative measure of latent immunological responses and (ii) a Bayes optimal classification of the different kinds of immunological response (c.f., glucose tolerance tests used to test for insulin resistance). This may be especially useful in assessing SARS-CoV-2 vaccines.


title: Solving the Optimal Trading Trajectory Problem Using Simulated Bifurcation

abstract: We use an optimization procedure based on simulated bifurcation (SB) to solve the integer portfolio and trading trajectory problem with an unprecedented computational speed. The underlying algorithm is based on a classical description of quantum adiabatic evolutions of a network of non-linearly interacting oscillators. This formulation has already proven to beat state of the art computation times for other NP-hard problems and is expected to show similar performance for certain portfolio optimization problems. Inspired by such we apply the SB approach to the portfolio integer optimization problem with quantity constraints and trading activities. We show first numerical results for portfolios of up to 1000 assets, which already confirm the power of the SB algorithm for its novel use-case as a portfolio and trading trajectory optimizer.


title: Not sure? Handling hesitancy of COVID-19 vaccines

abstract: From the moment the first COVID-19 vaccines are rolled out, there will need to be a large fraction of the global population ready in line. It is therefore crucial to start managing the growing global hesitancy to any such COVID-19 vaccine. The current approach of trying to convince the "no"s cannot work quickly enough, nor can the current policy of trying to find, remove and/or rebut all the individual pieces of COVID and vaccine misinformation. Instead, we show how this can be done in a simpler way by moving away from chasing misinformation content and focusing instead on managing the "yes--no--not-sure" hesitancy ecosystem.


title: Data-Driven Snapshot Calibration via Monotonic Feature Matching

abstract: Snapshot matrices of hyperbolic equations have a slow singular value decay, resulting in inefficient reduced-order models. We develop on the idea of inducing a faster singular value decay by computing snapshots on a transformed spatial domain, or the so-called snapshot calibration/transformation. We are particularly interested in problems involving shock collision, shock rarefaction-fan collision, shock formation, etc. For such problems, we propose a realizable algorithm to compute the spatial transform using monotonic feature matching. We consider discontinuities and kinks as features, and by carefully partitioning the parameter domain, we ensure that the spatial transform has properties that are desirable both from a theoretical and an implementation standpoint. We use these properties to prove that our method results in a fast m-width decay of a so-called calibrated manifold. A crucial observation we make is that due to calibration, the m-width does not only depend on m but also on the accuracy of the full order model, which is in contrast to elliptic and parabolic problems that do not need calibration. The method we propose only requires the solution snapshots and not the underlying partial differential equation (PDE) and is therefore, data-driven. We perform several numerical experiments to demonstrate the effectiveness of our method.


title: Steady state entropy production rate for scalar Langevin field theories

abstract: The entropy production rate (EPR) offers a quantitative measure of time reversal symmetry breaking in non-equilibrium systems. It can be defined either at particle level or at the level of coarse-grained fields such as density; the EPR for the latter quantifies the extent to which these coarse-grained fields behave irreversibly. In this work, we first develop a general method to compute the EPR of scalar Langevin field theories with additive noise. This large class of theories includes active versions of Model A (non-conserved density dynamics) and Model B (conserved) and also models where both types of dynamics are simultaneously present (such as Model AB). Treating the scalar field $\phi$ (and its time derivative $\dot\phi$) as the sole observable(s), we arrive at an expression for the EPR that is non-negative for every field configuration and is quadratic in the time-antisymmetric component of the dynamics. Our general expression is a function of the quasipotential, which determines the full probability distribution for configurations, and is not generally calculable. To alleviate this difficulty, we present a small-noise expansion of the EPR, which only requires knowledge of the deterministic (mean-field) solution for the scalar field in steady state, which generally is calculable, at least numerically. We demonstrate this calculation for the case of Model AB. We then present a similar EPR calculation for Model AB with the conservative and non-conservative contributions to $\dot\phi = \dot\phi_{\rm A} + \dot\phi_{\rm B}$ viewed as separately observable quantities. The results are qualitatively different, confirming that the field-level EPR depends on the choice of coarse-grained information retained within the dynamical description.


title: Near-Optimal Decremental Hopsets with Applications

abstract: Given a weighted undirected graph $G=(V,E,w)$, a hopset $H$ of hopbound $\beta$ and stretch $(1+\epsilon)$ is a set of edges such that for any pair of nodes $u, v \in V$, there is a path in $G \cup H$ of at most $\beta$ hops, whose length is within a $(1+\epsilon)$ factor from the distance between $u$ and $v$ in $G$. We show the first efficient decremental algorithm for maintaining hopsets with a polylogarithmic hopbound. The update time of our algorithm matches the best known static algorithm up to polylogarithmic factors. All the previous decremental hopset constructions had a superpolylogarithmic (but subpolynomial) hopbound of $2^{\log^{\Omega(1)} n}$ [Bernstein, FOCS'09; HKN, FOCS'14; Chechik, FOCS'18]. By applying our decremental hopset construction, we get improved or near optimal bounds for several distance problems. Most importantly, we show how to decrementally maintain $(2k-1)(1+\epsilon)$-approximate all-pairs shortest paths (for any constant $k \geq 2)$, in $\tilde{O}(n^{1/k})$ amortized update time and $O(k)$ query time. This improves (by a polynomial factor) over the update-time of the best previously known decremental algorithm in the constant query time regime. Moreover, it improves over the result of [Chechik, FOCS'18] that has a query time of $O(\log \log(nW))$, where $W$ is the aspect ratio, and the amortized update time is $n^{1/k}\cdot(\frac{1}{\epsilon})^{\tilde{O}(\sqrt{\log n})}$. For sparse graphs our construction nearly matches the best known static running time / query time tradeoff. We also obtain near-optimal bounds for maintaining approximate multi-source shortest paths and distance sketches, and get improved bounds for approximate single-source shortest paths. Our algorithms are randomized and our bounds hold with high probability against an oblivious adversary.


title: Gravitational radiation from birefringent matter dynamics

abstract: We present the second-order gravitational dynamics for a spacetime inhabited by matter fields which feature vacuum birefringence. The derivation follows a perturbative variant of the covariant constructive gravity program, ensuring diffeomorphism invariance of gravity and causal compatibility of matter theory and gravity. A subsequent spatio-temporal split of this theory reveals the presence of unphysical artifacts, which are cured by imposing constraints on the gravitational constants, reducing their number from ten to seven. Within this sector, we derive the gravitational radiation emitted by a binary system in circular motion. The system emits massless waves which correspond to the radiation predicted by Einstein gravity, but also massive waves, which are generated only above a certain angular frequency threshold and are unknown to Einstein gravity. A gravitational-wave detector modeled as sphere of freely falling test masses shows quantitatively and qualitatively new behavior under the influence of this radiation. The result is a prediction of gravitational self-coupling from first principles, demonstrating the predictive power of covariant constructive gravity for modified gravity research, especially in the era of gravitational-wave astronomy.


title: Regularisation by regular noise

abstract: We show that perturbing ill-posed differential equations with (potentially very) smooth random processes can restore well-posedness -- even if the perturbation is (potentially much) more regular than the drift component of the solution. The noise considered is of fractional Brownian type, and the familiar regularity condition $\alpha>1-1/(2H)$ is recovered for all non-integer $H>1$.


title: Moments and hybrid subconvexity for symmetric-square L-functions

abstract: We establish sharp bounds for the second moment of symmetric-square $L$-functions attached to Hecke Maass cusp forms $u_j$ with spectral parameter $t_j$, where the second moment is a sum over $t_j$ in a short interval. At the central point $s=1/2$ of the $L$-function, our interval is smaller than previous known results. More specifically, for $|t_j|$ of size $T$, our interval is of size $T^{1/5}$, while the previous best was $T^{1/3}$ from work of Lam. A little higher up on the critical line, our second moment yields a subconvexity bound for the symmetric-square $L$-function. More specifically, we get subconvexity at $s=1/2+it$ provided $|t_j|^{6/7+\delta}\le |t| \le (2-\delta)|t_j|$ for any fixed $\delta>0$. Since $|t|$ can be taken significantly smaller than $|t_j|$, this may be viewed as an approximation to the notorious subconvexity problem for the symmetric-square $L$-function in the spectral aspect at $s=1/2$.


title: Utilizing remote sensing data in forest inventory sampling via Bayesian optimization

abstract: In large-area forest inventories a trade-off between the amount of data to be sampled and the costs of collecting the data is necessary. It is not always possible to have a very large data sample when dealing with sampling-based inventories. It is therefore necessary to optimize the sampling design in order to achieve optimal population parameter estimation. On the contrary, the availability of remote sensing (RS) data correlated with the forest inventory variables is usually much higher. The combination of RS and the sampled field measurement data is often used for improving the forest inventory parameter estimation. In addition, it is also reasonable to study the utilization of RS data in inventory sampling, which can further improve the estimation of forest variables. In this study, we propose a data sampling method based on Bayesian optimization which uses RS data in forest inventory sample selection. The presented method applies the learned functional relationship between the RS and inventory data in new sampling decisions. We evaluate our method by conducting simulated sampling experiments with both synthetic data and measured data from the Aland region in Finland. The proposed method is benchmarked against two baseline methods: simple random sampling and the local pivotal method. The results of the simulated experiments show the best results in terms of MSE values for the proposed method when the functional relationship between RS and inventory data is correctly learned from the available training data.


title: Self-Adaptive Amplified Spontaneous Emission Suppression with a Photorefractive Two-Beam Coupling Filter

abstract: Amplified spontaneous emission is a source of broadband noise that parasitically limits the achievable gain in laser amplifiers. While optical bandpass filtering elements can suppress these broadband noise contributions, such filters are typically designed around specific frequencies or require manual tuning, rendering them less compatible with tunable laser systems. Here, we introduce a nonlinear self-adaptive filter and demonstrate the suppression of amplified spontaneous emission surrounding the lasing mode of a tunable 780 nm external cavity diode laser, using the two-beam coupling interaction in photorefractive BaTiO$_3$. A peak suppression of $-$10 dB is observed $\pm$2.5 nm from the lasing mode, with an overall 50% filter power throughput. The dynamic photorefractive filter is automatically centered on the peak frequency due to the continuous writing and readout of the volume holographic grating and can thereby also automatically adapt to frequency tuning, drift, or mode hopping with an estimated auto-tuning rate of 100 GHz/s under typical conditions. Additionally, we present opportunities for enhancing filter suppression characteristics via the input intensity ratio and tuning the bandwidth via the coupling angle, toward versatile, self-adaptive optical filtering.


title: Elastica: A compliant mechanics environment for soft robotic control

abstract: Soft robots are notoriously hard to control. This is partly due to the scarcity of models able to capture their complex continuum mechanics, resulting in a lack of control methodologies that take full advantage of body compliance. Currently available simulation methods are either too computational demanding or overly simplistic in their physical assumptions, leading to a paucity of available simulation resources for developing such control schemes. To address this, we introduce Elastica, a free, open-source simulation environment for soft, slender rods that can bend, twist, shear and stretch. We demonstrate how Elastica can be coupled with five state-of-the-art reinforcement learning algorithms to successfully control a soft, compliant robotic arm and complete increasingly challenging tasks.


title: Constraining the orientation of the spin axes of extrasolar minor bodies 1I/2017 U1 (`Oumuamua) and 2I/Borisov

abstract: Context. The orientation of the spin axis of a comet is defined by the values of its equatorial obliquity and its cometocentric longitude of the Sun. These parameters can be computed from the components of the nongravitational force caused by outgassing, if the cometary activity is well characterized. The trajectories of known interstellar bodies passing through the Solar System show nongravitational accelerations. Aims. The spin-axis orientation of 1I/2017 U1 (`Oumuamua) remains to be determined; for 2I/Borisov, the already released results are mutually exclusive. Here, we investigate the orientation of the spin axes of `Oumuamua and 2I/Borisov using public orbit determinations that consider outgassing. Methods. We applied a Monte Carlo simulation using the covariance matrix method together with Monte Carlo random search techniques to compute the distributions of equatorial obliquities and cometocentric longitudes of the Sun at perihelion of `Oumuamua and 2I/Borisov from the values of the nongravitational parameters. Results. We find that the equatorial obliquity of `Oumuamua could be about 93 deg, if it has a very prolate (fusiform) shape, or close to 16 deg, if it is very oblate (disk-like). Different orbit determinations of 2I/Borisov gave obliquity values of 59 deg and 90 deg. The distributions of cometocentric longitudes were in general multimodal. Conclusions. The most probable spin-axis direction of `Oumuamua in equatorial coordinates is (280 degrees, +46 degrees) if very prolate or (312 deg, -50 deg) if very oblate. For the orbit determinations of 2I/Borisov used here, we find most probable poles pointing near (275 deg, +65 deg) and (231 deg, +30 deg), respectively. Although our analysis favors an oblate shape for 2I/Borisov, a prolate one cannot be ruled out.


title: Modeling Task Effects on Meaning Representation in the Brain via Zero-Shot MEG Prediction

abstract: How meaning is represented in the brain is still one of the big open questions in neuroscience. Does a word (e.g., bird) always have the same representation, or does the task under which the word is processed alter its representation (answering "can you eat it?" versus "can it fly?")? The brain activity of subjects who read the same word while performing different semantic tasks has been shown to differ across tasks. However, it is still not understood how the task itself contributes to this difference. In the current work, we study Magnetoencephalography (MEG) brain recordings of participants tasked with answering questions about concrete nouns. We investigate the effect of the task (i.e. the question being asked) on the processing of the concrete noun by predicting the millisecond-resolution MEG recordings as a function of both the semantics of the noun and the task. Using this approach, we test several hypotheses about the task-stimulus interactions by comparing the zero-shot predictions made by these hypotheses for novel tasks and nouns not seen during training. We find that incorporating the task semantics significantly improves the prediction of MEG recordings, across participants. The improvement occurs 475-550ms after the participants first see the word, which corresponds to what is considered to be the ending time of semantic processing for a word. These results suggest that only the end of semantic processing of a word is task-dependent, and pose a challenge for future research to formulate new hypotheses for earlier task effects as a function of the task and stimuli.


title: Superuniversality of superdiffusion

abstract: Anomalous finite-temperature transport has recently been observed in numerical studies of various integrable models in one dimension; these models share the feature of being invariant under a continuous non-abelian global symmetry. This work offers a comprehensive group-theoretic account of this elusive phenomenon. For an integrable quantum model invariant under a global non-abelian simple Lie group $G$, we find that finite-temperature transport of Noether charges associated with symmetry $G$ in thermal states that are invariant under $G$ is universally superdiffusive and characterized by dynamical exponent $z = 3/2$. This conclusion holds regardless of the Lie algebra symmetry, local degrees of freedom (on-site representations), Lorentz invariance, or particular realization of microscopic interactions: we accordingly dub it as superuniversal. The anomalous transport behavior is attributed to long-lived giant quasiparticles dressed by thermal fluctuations. We provide an algebraic viewpoint on the corresponding dressing transformation and elucidate formal connections to fusion identities amongst the quantum-group characters. We identify giant quasiparticles with nonlinear soliton modes of classical field theories that describe low-energy excitations above ferromagnetic vacua. Our analysis of these field theories also provides a complete classification of the low-energy (i.e., Goldstone-mode) spectra of quantum isotropic ferromagnetic chains.

categories: cond-mat.stat-mech cond-mat.str-el math-ph math.MP


title: On the cohomology of Lie algebras with an invariant inner product

abstract: In this work we consider all metric Lie algebras, having a nondegenerate symmetric invariant bilinear form, over \C and \R up to dimension 5 and all metric Lie algebras over \C in dimension 6. We introduce cyclic and reduced cyclic cohomology to identify their metric deformations.


title: Discovering Dynamic Salient Regions for Spatio-Temporal Graph Neural Networks

abstract: Graph Neural Networks are perfectly suited to capture latent interactions between various entities in the spatio-temporal domain (e.g. videos). However, when an explicit structure is not available, it is not obvious what atomic elements should be represented as nodes. Current works generally use pre-trained object detectors or fixed, predefined regions to extract graph nodes. Improving upon this, our proposed model learns nodes that dynamically attach to well-delimited salient regions, which are relevant for a higher-level task, without using any object-level supervision. Constructing these localized, adaptive nodes gives our model inductive bias towards object-centric representations and we show that it discovers regions that are well correlated with objects in the video. In extensive ablation studies and experiments on two challenging datasets, we show superior performance to previous graph neural networks models for video classification.


title: Radar-Camera Sensor Fusion for Joint Object Detection and Distance Estimation in Autonomous Vehicles

abstract: In this paper we present a novel radar-camera sensor fusion framework for accurate object detection and distance estimation in autonomous driving scenarios. The proposed architecture uses a middle-fusion approach to fuse the radar point clouds and RGB images. Our radar object proposal network uses radar point clouds to generate 3D proposals from a set of 3D prior boxes. These proposals are mapped to the image and fed into a Radar Proposal Refinement (RPR) network for objectness score prediction and box refinement. The RPR network utilizes both radar information and image feature maps to generate accurate object proposals and distance estimations. The radar-based proposals are combined with image-based proposals generated by a modified Region Proposal Network (RPN). The RPN has a distance regression layer for estimating distance for every generated proposal. The radar-based and image-based proposals are merged and used in the next stage for object classification. Experiments on the challenging nuScenes dataset show our method outperforms other existing radar-camera fusion methods in the 2D object detection task while at the same time accurately estimates objects' distances.


title: Sensitivity of steady states in a degenerately-damped stochastic Lorenz system

abstract: We study stability of solutions for a randomly driven and degenerately damped version of the Lorenz '63 model. Specifically, we prove that when damping is absent in one of the temperature components, the system possesses a unique invariant probability measure if and only if noise acts on the convection variable. On the other hand, if there is a positive growth term on the vertical temperature profile, we prove that there is no normalizable invariant state. Our approach relies on the derivation and analysis of non-trivial Lyapunov functions which ensure positive recurrence or null-recurrence/transience of the dynamics.


title: Analytical Approximation of the Second-Harmonic Conversion Efficiency

abstract: The second-harmonic generation process of a focused laser beam inside a nonlinear crystal is described by the Boyd-Kleinman theory. Calculating the actual conversion efficiency and upconverted power requires the solution of a double integral that is analytically intractable. We provide an expression that predicts the exact gain coefficient within an error margin of less than 2% over several orders of magnitude of the confocal parameter and as a function of the walk-off parameter. Our result allows for readily tuning the beam parameters to optimize the performance of the upconversion process and improve optical system designs.


title: Integrals of functions containing parameters

abstract: This paper offers what seems at first to be a minor technical correction to the current practice of computing indefinite integrals, and introduces the idea of a "Kahanian constant of integration". However, the total impact of this minor correction is potentially large because the current practice is taught early at the university level and to very many students---most of whom do not go on to become mathematics majors. Moreover, computer algebra systems have become widespread, including good free ones, some of which are available for smartphones. Most current computer algebra systems apply current textbook rules and amplify the effects of fundamental "minor" errors such as the error in continuity that we address in this article. So in practice, the correction we present is important.


title: A Time To Event Framework For Multi-touch Attribution

abstract: Multi-touch attribution (MTA) estimates the relative contributions of the multiple ads a user may see prior to any observed conversions. Increasingly, advertisers also want to base budget and bidding decisions on these attributions, spending more on ads that drive more conversions. We describe two requirements for an MTA system to be suitable for this application: First, it must be able to handle continuously updated and incomplete data. Second, it must be sufficiently flexible to capture that an ad's effect will change over time. We describe an MTA system, consisting of a model for user conversion behavior and a credit assignment algorithm, that satisfies these requirements. Our model for user conversion behavior treats conversions as occurrences in an inhomogeneous Poisson process, while our attribution algorithm is based on iteratively removing the last ad in the path.


title: On the global controllability of scalar conservation laws with boundary and source controls

abstract: We provide global and semi-global controllability results for hyperbolic conservation laws on a bounded domain, with a general (not necessarily convex)flux and a time-dependent source term acting as a control. The results are achieved for, possibly critical, both continuously differentiable states and BV states. The proofs are based on a combination of the return method and on the analysis of the Riccati equaiton for the space derivative of the solution.


title: Resource distillation in convex Gaussian resource theories

abstract: It is known that distillation in continuous variable resource theories is impossible when restricted to Gaussian states and operations. To overcome this limitation, we enlarge the theories to include convex mixtures of Gaussian states and operations. This extension is operationally well-motivated since classical randomness is easily accessible. We find that resource distillation becomes possible for convex Gaussian resource theories-albeit in a limited fashion. We derive this limitation by studying the convex roof extension of a Gaussian resource measure and then go on to show that our bound is tight by means of example protocols for the distillation of squeezing and entanglement.


title: Large Norms of CNN Layers Do Not Hurt Adversarial Robustness

abstract: Since the Lipschitz properties of convolutional neural networks (CNNs) are widely considered to be related to adversarial robustness, we theoretically characterize the $\ell_1$ norm and $\ell_\infty$ norm of 2D multi-channel convolutional layers and provide efficient methods to compute the exact $\ell_1$ norm and $\ell_\infty$ norm. Based on our theorem, we propose a novel regularization method termed norm decay, which can effectively reduce the norms of convolutional layers and fully-connected layers. Experiments show that norm-regularization methods, including norm decay, weight decay, and singular value clipping, can improve generalization of CNNs. However, they can slightly hurt adversarial robustness. Observing this unexpected phenomenon, we compute the norms of layers in the CNNs trained with three different adversarial training frameworks and surprisingly find that adversarially robust CNNs have comparable or even larger layer norms than their non-adversarially robust counterparts. Furthermore, we prove that under a mild assumption, adversarially robust classifiers can be achieved using neural networks, and an adversarially robust neural network can have an arbitrarily large Lipschitz constant. For this reason, enforcing small norms on CNN layers may be neither necessary nor effective in achieving adversarial robustness. The code is available at https://github.com/youweiliang/norm_robustness.


title: A global microscopic description of nucleon-nucleus scattering with quantified uncertainties

abstract: We develop for the first time a microscopic global nucleon-nucleus optical potential with quantified uncertainties suitable for analyzing nuclear reaction experiments at next-generation rare-isotope beam facilities. Within the improved local density approximation and without any adjustable parameters, we begin by computing proton-nucleus and neutron-nucleus optical potentials from a set of five nuclear forces from chiral effective field theory for 1800 target nuclei in the mass range 12 $\leq$ A $\leq$ 242 for energies between 0 MeV $<$ E $\lesssim$ 150 MeV. We then parameterize a global optical potential for each chiral force that depends smoothly on the projectile energy as well as the target nucleus mass number and isospin asymmetry. Uncertainty bands for elastic scattering observables are generated from a full covariance analysis of the parameters entering in the description of our global optical potential and benchmarked against existing experimental data for stable target nuclei. Since our approach is purely microscopic, we anticipate a similar quality of the model for nucleon scattering on unstable isotopes.


title: FIGARO: Improving System Performance via Fine-Grained In-DRAM Data Relocation and Caching

abstract: DRAM Main memory is a performance bottleneck for many applications due to the high access latency. In-DRAM caches work to mitigate this latency by augmenting regular-latency DRAM with small-but-fast regions of DRAM that serve as a cache for the data held in the regular-latency region of DRAM. While an effective in-DRAM cache can allow a large fraction of memory requests to be served from a fast DRAM region, the latency savings are often hindered by inefficient mechanisms for relocating copies of data into and out of the fast regions. Existing in-DRAM caches have two sources of inefficiency: (1) the data relocation granularity is an entire multi-kilobyte row of DRAM; and (2) because the relocation latency increases with the physical distance between the slow and fast regions, multiple fast regions are physically interleaved among slow regions to reduce the relocation latency, resulting in increased hardware area and manufacturing complexity. We propose a new substrate, FIGARO, that uses existing shared global buffers among subarrays within a DRAM bank to provide support for in-DRAM data relocation across subarrays at the granularity of a single cache block. FIGARO has a distance-independent latency within a DRAM bank, and avoids complex modifications to DRAM. Using FIGARO, we design a fine-grained in-DRAM cache called FIGCache. The key idea of FIGCache is to cache only small, frequently-accessed portions of different DRAM rows in a designated region of DRAM. By caching only the parts of each row that are expected to be accessed in the near future, we can pack more of the frequently-accessed data into FIGCache, and can benefit from additional row hits in DRAM. Our evaluations show that FIGCache improves the average performance of a system using DDR4 DRAM by 16.3% and reduces average DRAM energy consumption by 7.8% for 8-core workloads, over a conventional system without in-DRAM caching.


title: Competitiveness of MAP-Elites against Proximal Policy Optimization on locomotion tasks in deterministic simulations

abstract: The increasing importance of robots and automation creates a demand for learnable controllers which can be obtained through various approaches such as Evolutionary Algorithms (EAs) or Reinforcement Learning (RL). Unfortunately, these two families of algorithms have mainly developed independently and there are only a few works comparing modern EAs with deep RL algorithms. We show that Multidimensional Archive of Phenotypic Elites (MAP-Elites), which is a modern EA, can deliver better-performing solutions than one of the state-of-the-art RL methods, Proximal Policy Optimization (PPO) in the generation of locomotion controllers for a simulated hexapod robot. Additionally, extensive hyper-parameter tuning shows that MAP-Elites displays greater robustness across seeds and hyper-parameter sets. Generally, this paper demonstrates that EAs combined with modern computational resources display promising characteristics and have the potential to contribute to the state-of-the-art in controller learning.


title: A Cantor-Bernstein theorem for infinite matroids

abstract: We give a common matroidal generalisation of `A Cantor-Bernstein theorem for paths in graphs' by Diestel and Thomassen and `A Cantor-Bernstein-type theorem for spanning trees in infinite graphs' by ourselves.


title: Metrological complementarity reveals the Einstein-Podolsky-Rosen paradox

abstract: The Einstein-Podolsky-Rosen (EPR) paradox plays a fundamental role in our understanding of quantum mechanics, and is associated with the possibility of predicting the results of non-commuting measurements with a precision that seems to violate the uncertainty principle. This apparent contradiction to complementarity is made possible by nonclassical correlations stronger than entanglement, called steering. Quantum information recognises steering as an essential resource for a number of tasks but, contrary to entanglement, its role for metrology has so far remained unclear. Here, we formulate the EPR paradox in the framework of quantum metrology, showing that it enables the precise estimation of a local phase shift and of its generating observable. Employing a stricter formulation of quantum complementarity, we derive a criterion based on the quantum Fisher information that detects steering in a larger class of states than well-known uncertainty-based criteria. Our result identifies useful steering for quantum-enhanced precision measurements and allows one to uncover steering of non-Gaussian states in state-of-the-art experiments.


title: A Computational Approach to Understanding Empathy Expressed in Text-Based Mental Health Support

abstract: Empathy is critical to successful mental health support. Empathy measurement has predominantly occurred in synchronous, face-to-face settings, and may not translate to asynchronous, text-based contexts. Because millions of people use text-based platforms for mental health support, understanding empathy in these contexts is crucial. In this work, we present a computational approach to understanding how empathy is expressed in online mental health platforms. We develop a novel unifying theoretically-grounded framework for characterizing the communication of empathy in text-based conversations. We collect and share a corpus of 10k (post, response) pairs annotated using this empathy framework with supporting evidence for annotations (rationales). We develop a multi-task RoBERTa-based bi-encoder model for identifying empathy in conversations and extracting rationales underlying its predictions. Experiments demonstrate that our approach can effectively identify empathic conversations. We further apply this model to analyze 235k mental health interactions and show that users do not self-learn empathy over time, revealing opportunities for empathy training and feedback.


title: On the Cauchy problem for the Muskat equation. II: Critical initial data

abstract: We prove that the Cauchy problem for the Muskat equation is well-posed locally in time for any initial data in the critical space of Lipschitz functions with three-half derivative in $L^2$. Moreover, we prove that the solution exists globally in time under a smallness assumption.


title: Tropical time series, iterated-sums signatures and quasisymmetric functions

abstract: Aiming for a systematic feature-extraction from time series, we introduce the iterated-sums signature over arbitrary commutative semirings. The case of the tropical semiring is a central, and our motivating example. It leads to features of (real-valued) time series that are not easily available using existing signature-type objects. We demonstrate how the signature extracts chronological aspects of a time series, and that its calculation is possible in linear time. We identify quasisymmetric expressions over semirings as the appropriate framework for iterated-sums signatures over semiring-valued time series.


title: A five distance theorem for Kronecker sequences

abstract: The three distance theorem (also known as the three gap theorem or Steinhaus problem) states that, for any given real number $\alpha$ and integer $N$, there are at most three values for the distances between consecutive elements of the Kronecker sequence $\alpha, 2\alpha,\ldots, N\alpha$ mod 1. In this paper we consider a natural generalisation of the three distance theorem to the higher dimensional Kronecker sequence $\vec\alpha, 2\vec\alpha,\ldots, N\vec\alpha$ modulo an integer lattice. We prove that in two dimensions there are at most five values that can arise as a distance between nearest neighbors, for all choices of $\vec\alpha$ and $N$. Furthermore, for almost every $\vec\alpha$, five distinct distances indeed appear for infinitely many $N$ and hence five is the best possible general upper bound. In higher dimensions we have similar explicit, but less precise, upper bounds. For instance in three dimensions our bound is 13, though we conjecture the truth to be 9. We furthermore study the number of possible distances from a point to its nearest neighbor in a restricted cone of directions. This may be viewed as a generalisation of the gap length in one dimension. For large cone angles we use geometric arguments to produce explicit bounds directly analogous to the three distance theorem. For small cone angles we use ergodic theory of homogeneous flows in the space of unimodular lattices to show that the number of distinct lengths is (a) unbounded for almost all $\vec\alpha$ and (b) bounded for $\vec\alpha$ that satisfy certain Diophantine conditions.


title: Self-Supervised Meta-Learning for Few-Shot Natural Language Classification Tasks

abstract: Self-supervised pre-training of transformer models has revolutionized NLP applications. Such pre-training with language modeling objectives provides a useful initial point for parameters that generalize well to new tasks with fine-tuning. However, fine-tuning is still data inefficient -- when there are few labeled examples, accuracy can be low. Data efficiency can be improved by optimizing pre-training directly for future fine-tuning with few examples; this can be treated as a meta-learning problem. However, standard meta-learning techniques require many training tasks in order to generalize; unfortunately, finding a diverse set of such supervised tasks is usually difficult. This paper proposes a self-supervised approach to generate a large, rich, meta-learning task distribution from unlabeled text. This is achieved using a cloze-style objective, but creating separate multi-class classification tasks by gathering tokens-to-be blanked from among only a handful of vocabulary terms. This yields as many unique meta-training tasks as the number of subsets of vocabulary terms. We meta-train a transformer model on this distribution of tasks using a recent meta-learning framework. On 17 NLP tasks, we show that this meta-training leads to better few-shot generalization than language-model pre-training followed by finetuning. Furthermore, we show how the self-supervised tasks can be combined with supervised tasks for meta-learning, providing substantial accuracy gains over previous supervised meta-learning.


title: VALES VII: Molecular and ionized gas properties in pressure balanced interstellar medium of starburst galaxies at z ~ 0.15

abstract: Context. Spatially resolved observations of the ionized and molecular gas are critical for understanding the physical processes that govern the interstellar medium (ISM) in galaxies. Aims. To study the morpho-kinematic properties of the ionized and molecular gas in three dusty starburst galaxies at $z = 0.12-0.17$ to explore the relation between molecular ISM gas phase dynamics and the star-formation activity. Methods. We analyse $\sim$kpc-scale ALMA CO(1--0) and seeing limited SINFONI Paschen-$\alpha$ observations. We use a dynamical mass model, which accounts for beam-smearing effects, to constrain the CO-to-H$_2$ conversion factor. Results. One starburst galaxy shows irregular morphology which may indicate a major merger, while the other two systems show disc-like morpho-kinematics. The two disc-like starbursts show molecular gas velocity dispersion values comparable with that seen in local LIRG/ULIRGs, but in an ISM with molecular gas fraction and surface density values consistent to that reported for local star-forming galaxies. These molecular gas velocity dispersion values can be explained by assuming vertical pressure equilibrium. The star-formation activity is correlated with the molecular gas content suggesting depletion times of the order of $\sim 0.1-1$ Gyr. The star formation rate surface density ($\Sigma_{\rm SFR}$) correlates with the ISM pressure set by self-gravity ($P_{\rm grav}$) following a power law with an exponent close to 0.8. Conclusions. In dusty disc-like starburst galaxies, our data support the scenario in which the molecular gas velocity dispersion values are driven by the ISM pressure set by self-gravity, responsible to maintain the vertical pressure balance. The correlation between $\Sigma_{\rm SFR}$ and $P_{\rm grav}$ suggests that, in these dusty starbursts galaxies, the star formation activity arises as a consequence of the ISM pressure balance.


title: Coordinate Methods for Matrix Games

abstract: We develop primal-dual coordinate methods for solving bilinear saddle-point problems of the form $\min_{x \in \mathcal{X}} \max_{y\in\mathcal{Y}} y^\top A x$ which contain linear programming, classification, and regression as special cases. Our methods push existing fully stochastic sublinear methods and variance-reduced methods towards their limits in terms of per-iteration complexity and sample complexity. We obtain nearly-constant per-iteration complexity by designing efficient data structures leveraging Taylor approximations to the exponential and a binomial heap. We improve sample complexity via low-variance gradient estimators using dynamic sampling distributions that depend on both the iterates and the magnitude of the matrix entries. Our runtime bounds improve upon those of existing primal-dual methods by a factor depending on sparsity measures of the $m$ by $n$ matrix $A$. For example, when rows and columns have constant $\ell_1/\ell_2$ norm ratios, we offer improvements by a factor of $m+n$ in the fully stochastic setting and $\sqrt{m+n}$ in the variance-reduced setting. We apply our methods to computational geometry problems, i.e. minimum enclosing ball, maximum inscribed ball, and linear regression, and obtain improved complexity bounds. For linear regression with an elementwise nonnegative matrix, our guarantees improve on exact gradient methods by a factor of $\sqrt{\mathrm{nnz}(A)/(m+n)}$.


title: Motivic Atiyah-Segal completion theorem

abstract: Let T be a torus, X a smooth quasi-compact separated scheme equipped with a T-action, and [X/T] the associated quotient stack. Given any localizing A1-homotopy invariant of dg categories E, we prove that the derived completion of E([X/T]) at the augmentation ideal I of the representation ring R(T) of T agrees with the Borel construction associated to the T-action on X. Moreover, for certain localizing A1-homotopy invariants, we extend this result to the case of a linearly reductive group scheme G. As a first application, we obtain an alternative proof of Krishna's completion theorem in algebraic K-theory, of Thomason's completion theorem in \'etale K-theory with coefficients, and also of Atiyah-Segal's completion theorem in topological K-theory. These alternative proofs lead to a spectral enrichment of the corresponding completion theorems and also to the following improvements: in the case of Thomason's completion theorem the base field no longer needs to be separably closed, and in the case of Atiyah-Segal's completion theorem the topological spaces no longer needs to be compact and the equivariant topological K-theory groups no longer need to be finitely generated over the representation ring. As a second application, we obtain new completion theorems in l-adic \'etale K-theory, in (real) semi-topological K-theory and also in periodic cyclic homology. As a third application, we obtain a purely algebraic description of the different equivariant cohomology groups in the literature (motivic, l-adic, (real) morphic, Betti, de Rham, etc). Finally, in two appendixes of independent interest, we extend a result of Weibel on homotopy K-theory from the realm of schemes to the broad setting of quotient stacks and establish some useful properties of (real) semi-topological K-theory.


title: 'Less Than One'-Shot Learning: Learning N Classes From M<N Samples

abstract: Deep neural networks require large training sets but suffer from high computational cost and long training times. Training on much smaller training sets while maintaining nearly the same accuracy would be very beneficial. In the few-shot learning setting, a model must learn a new class given only a small number of samples from that class. One-shot learning is an extreme form of few-shot learning where the model must learn a new class from a single example. We propose the `less than one'-shot learning task where models must learn $N$ new classes given only $M<N$ examples and we show that this is achievable with the help of soft labels. We use a soft-label generalization of the k-Nearest Neighbors classifier to explore the intricate decision landscapes that can be created in the `less than one'-shot learning setting. We analyze these decision landscapes to derive theoretical lower bounds for separating $N$ classes using $M<N$ soft-label samples and investigate the robustness of the resulting systems.


title: Fragmentation production of fully-charmed tetraquarks at LHC

abstract: The $X(6900)$ resonance, very recently discovered in the double-$J/\psi$ channel at LHCb experiment, has spurred intensive interest in unravelling the nature of the fully charmed tetraquark state. The aim of this paper is to present a model-independent theoretical framework to study the inclusive production of this novel species of exotic hadrons, the resonances composed of four heavy quark (commonly referred to as $T_{4c}$), at large $p_T$ in hadron collision experiments. Appealing to asymptotic freedom and the fact $m_c\gg \Lambda_{\rm QCD}$, we propose that the nonpertubative yet universal gluon-to-$T_{4c}$ fragmentation function, can be decomposed into the product of the perturbatively calculable short-distance coefficient and the long-distance NRQCD matrix elements. We compute the short-distance coefficient at lowest-order in $\alpha_s$ and velocity expansion. Adopting the diquark ansatz to roughly estimate those not-yet-known NRQCD matrix elements, together with the standard QCD factorization theorem, we predict the differential production rates for the $T_{4c/4b}(0^{++})$ and $T_{4c/4b}(2^{++})$ at large $p_T$ in $pp$ collision, which eagerly awaits the confrontation with the future LHC experiments.


title: MSTREAM: Fast Anomaly Detection in Multi-Aspect Streams

abstract: Given a stream of entries in a multi-aspect data setting i.e., entries having multiple dimensions, how can we detect anomalous activities in an unsupervised manner? For example, in the intrusion detection setting, existing work seeks to detect anomalous events or edges in dynamic graph streams, but this does not allow us to take into account additional attributes of each entry. Our work aims to define a streaming multi-aspect data anomaly detection framework, termed MSTREAM which can detect unusual group anomalies as they occur, in a dynamic manner. MSTREAM has the following properties: (a) it detects anomalies in multi-aspect data including both categorical and numeric attributes; (b) it is online, thus processing each record in constant time and constant memory; (c) it can capture the correlation between multiple aspects of the data. MSTREAM is evaluated over the KDDCUP99, CICIDS-DoS, UNSW-NB 15 and CICIDS-DDoS datasets, and outperforms state-of-the-art baselines.


title: Real-Time Anomaly Detection in Edge Streams

abstract: Given a stream of graph edges from a dynamic graph, how can we assign anomaly scores to edges in an online manner, for the purpose of detecting unusual behavior, using constant time and memory? Existing approaches aim to detect individually surprising edges. In this work, we propose MIDAS, which focuses on detecting microcluster anomalies, or suddenly arriving groups of suspiciously similar edges, such as lockstep behavior, including denial of service attacks in network traffic data. We further propose MIDAS-F, to solve the problem by which anomalies are incorporated into the algorithm's internal states, creating a `poisoning' effect that can allow future anomalies to slip through undetected. MIDAS-F introduces two modifications: 1) We modify the anomaly scoring function, aiming to reduce the `poisoning' effect of newly arriving edges; 2) We introduce a conditional merge step, which updates the algorithm's data structures after each time tick, but only if the anomaly score is below a threshold value, also to reduce the `poisoning' effect. Experiments show that MIDAS-F has significantly higher accuracy than MIDAS. MIDAS has the following properties: (a) it detects microcluster anomalies while providing theoretical guarantees about its false positive probability; (b) it is online, thus processing each edge in constant time and constant memory, and also processes the data orders-of-magnitude faster than state-of-the-art approaches; (c) it provides up to 62% higher ROC-AUC than state-of-the-art approaches.


title: MEAL V2: Boosting Vanilla ResNet-50 to 80%+ Top-1 Accuracy on ImageNet without Tricks

abstract: We introduce a simple yet effective distillation framework that is able to boost the vanilla ResNet-50 to 80%+ Top-1 accuracy on ImageNet without tricks. We construct such a framework through analyzing the problems in the existing classification system and simplify the base method ensemble knowledge distillation via discriminators by: (1) adopting the similarity loss and discriminator only on the final outputs and (2) using the average of softmax probabilities from all teacher ensembles as the stronger supervision. Intriguingly, three novel perspectives are presented for distillation: (1) weight decay can be weakened or even completely removed since the soft label also has a regularization effect; (2) using a good initialization for students is critical; and (3) one-hot/hard label is not necessary in the distillation process if the weights are well initialized. We show that such a straight-forward framework can achieve state-of-the-art results without involving any commonly-used techniques, such as architecture modification; outside training data beyond ImageNet; autoaug/randaug; cosine learning rate; mixup/cutmix training; label smoothing; etc. Our method obtains 80.67% top-1 accuracy on ImageNet using a single crop-size of 224x224 with vanilla ResNet-50, outperforming the previous state-of-the-arts by a significant margin under the same network structure. Our result can be regarded as a strong baseline using knowledge distillation, and to our best knowledge, this is also the first method that is able to boost vanilla ResNet-50 to surpass 80% on ImageNet without architecture modification or additional training data. On smaller ResNet-18, our distillation framework consistently improves from 69.76% to 73.19%, which shows tremendous practical values in real-world applications. Our code and models are available at: https://github.com/szq0214/MEAL-V2.


title: ExGAN: Adversarial Generation of Extreme Samples

abstract: Mitigating the risk arising from extreme events is a fundamental goal with many applications, such as the modelling of natural disasters, financial crashes, epidemics, and many others. To manage this risk, a vital step is to be able to understand or generate a wide range of extreme scenarios. Existing approaches based on Generative Adversarial Networks (GANs) excel at generating realistic samples, but seek to generate typical samples, rather than extreme samples. Hence, in this work, we propose ExGAN, a GAN-based approach to generate realistic and extreme samples. To model the extremes of the training distribution in a principled way, our work draws from Extreme Value Theory (EVT), a probabilistic approach for modelling the extreme tails of distributions. For practical utility, our framework allows the user to specify both the desired extremeness measure, as well as the desired extremeness probability they wish to sample at. Experiments on real US Precipitation data show that our method generates realistic samples, based on visual inspection and quantitative measures, in an efficient manner. Moreover, generating increasingly extreme examples using ExGAN can be done in constant time (with respect to the extremeness probability $\tau$), as opposed to the $\mathcal{O}(\frac{1}{\tau})$ time required by the baseline approach.


title: Recombination-independent determination of the sound horizon and the Hubble constant from BAO

abstract: The Hubble tension and attempts to resolve it by modifying the physics of (or at) recombination motivate finding ways to determine $H_0$ and the sound horizon at the epoch of baryon decoupling $r_{\rm d}$ in ways that neither rely on a recombination model nor on late-time Hubble data. In this work, we investigate what one can learn from the current and future BAO data when treating $r_{\rm d}$ and $H_0$ as independent free parameters. It is well known that BAO gives exquisite constraints on the product $r_{\rm d}H_0$. We show here that imposing a moderate prior on $\Omega_{\rm m} h^2$ breaks the degeneracy between $r_{\rm d}$ and $H_0$. Using the latest BAO data, including the recently released eBOSS DR16, along with a $\Omega_{\rm m} h^2$ prior based on the Planck best fit $\Lambda$CDM model, we find $r_{\rm d} =143.7 \pm 2.7$ Mpc and $H_0 = 69.6 \pm 1.8$ km/s/Mpc. BAO data therefore prefers somewhat lower $r_{\rm d}$ and higher $H_0$ than those inferred from Planck data in a $\Lambda$CDM model. We find similar values when combing BAO with the Pantheon supernovae, DES galaxy weak lensing, Planck or SPTPol CMB lensing and the cosmic chronometers data. We perform a forecast for DESI and find that, when aided with a moderate prior on $\Omega_{\rm m} h^2$, DESI will measure $r_{\rm d}$ and $H_0$ without assuming a recombination model with an accuracy surpassing the current best estimates from Planck.


title: Capturing Richer Information -- On Establishing the Validity of an Interval-Valued Survey Response Mode

abstract: Obtaining quantitative survey responses that are both accurate and informative is crucial to a wide range of fields. Traditional and ubiquitous response formats such as Likert and Visual Analogue Scales require condensation of responses into discrete point values - but sometimes a range of options may better represent the correct answer. In this paper, we propose an efficient interval-valued response mode, whereby responses are made by marking an ellipse along a continuous scale. We discuss its potential to capture and quantify valuable information that would be lost using conventional approaches, while preserving a high degree of response-efficiency. The information captured by the response interval may represent a possible response range - i.e., a conjunctive set, such as the real numbers between three and six. Alternatively, it may reflect uncertainty in respect to a distinct response - i.e., a disjunctive set, such as a confidence interval. We then report a validation study, utilizing our recently introduced open-source software (DECSYS) to explore how interval-valued survey responses reflect experimental manipulations of several factors hypothesised to influence interval width, across multiple contexts. Results consistently indicate that respondents used interval widths effectively, and subjective participant feedback was also positive. We present this as initial empirical evidence for the efficacy and value of interval-valued response capture. Interestingly, our results also provide insight into respondents' reasoning about the different aforementioned types of intervals - we replicate a tendency towards overconfidence for those representing epistemic uncertainty (i.e., disjunctive sets), but find intervals representing inherent range (i.e., conjunctive sets) to be well-calibrated.


title: Online Semi-Supervised Learning in Contextual Bandits with Episodic Reward

abstract: We considered a novel practical problem of online learning with episodically revealed rewards, motivated by several real-world applications, where the contexts are nonstationary over different episodes and the reward feedbacks are not always available to the decision making agents. For this online semi-supervised learning setting, we introduced Background Episodic Reward LinUCB (BerlinUCB), a solution that easily incorporates clustering as a self-supervision module to provide useful side information when rewards are not observed. Our experiments on a variety of datasets, both in stationary and nonstationary environments of six different scenarios, demonstrated clear advantages of the proposed approach over the standard contextual bandit. Lastly, we introduced a relevant real-life example where this problem setting is especially useful.


title: The MUSE Deep Lensed Field on the Hubble Frontier Field MACS~J0416

abstract: Context: A census of faint and tiny star forming complexes at high redshift is key to improving our understanding of reionizing sources, galaxy growth and the formation of globular clusters. Aims: We present the MUSE Deep Lensed Field (MDLF) program. Methods: We describe Deep MUSE observations of 17.1 hours integration on a single pointing over the Hubble Frontier Field galaxy cluster MACS~J0416. Results: We confirm spectroscopic redshifts for all 136 multiple images of 48 source galaxies at 0.9<z<6.2. Within those galaxies, we securely identify 182 multiple images of 66 galaxy components that we use to constrain our lens model. We identify 116 clumps belonging to background high-z galaxies; the majority of them are multiple images and span magnitude, size and redshift intervals of [-18,-10], [~400-3] parsec and 1<z<6.6, respectively, with the most magnified ones probing possible single gravitationally bound star clusters. The depth of the MDLF combined with lensing magnification lead us to reach a detection limit for unresolved emission lines of a few 10$^{-20}$ erg/s/cm2, after correction for lensing magnification. Ultraviolet high-ionization metal lines (and HeII1640) are detected with S/N>10 for individual objects down to de-lensed magnitude 28-30 suggesting that they are common in such faint sources. Conclusions:Deep MUSE observations, in combination with existing HST imaging, allowed us to:(1) confirm redshifts for extremely faint high-z sources;(2) peer into their internal clumps (down to 100-200 pc scale);(3) in some cases break down such clumps into bound star clusters (<20 pc scale);(4) double the number of constraints for the lens model,reaching an unprecedented set of 182 bona-fide multiple images and confirming up to 213 galaxy cluster members. These results demonstrate the power that JWST and future ELTs will have when combined to study gravitational telescopes.[abridged]


title: Likelihood-free inference with neural compression of DES SV weak lensing map statistics

abstract: In many cosmological inference problems, the likelihood (the probability of the observed data as a function of the unknown parameters) is unknown or intractable. This necessitates approximations and assumptions, which can lead to incorrect inference of cosmological parameters, including the nature of dark matter and dark energy, or create artificial model tensions. Likelihood-free inference covers a novel family of methods to rigorously estimate posterior distributions of parameters using forward modelling of mock data. We present likelihood-free cosmological parameter inference using weak lensing maps from the Dark Energy Survey (DES) SV data, using neural data compression of weak lensing map summary statistics. We explore combinations of the power spectra, peak counts, and neural compressed summaries of the lensing mass map using deep convolution neural networks. We demonstrate methods to validate the inference process, for both the data modelling and the probability density estimation steps. Likelihood-free inference provides a robust and scalable alternative for rigorous large-scale cosmological inference with galaxy survey data (for DES, Euclid and LSST). We have made our simulated lensing maps publicly available.


title: Causality Constraints in Large $N$ QCD Coupled to Gravity

abstract: Confining gauge theories contain glueballs and mesons with arbitrary spin, and these particles become metastable at large $N$. However, metastable higher spin particles, when coupled to gravity, are in conflict with causality. This tension can be avoided only if the gravitational interaction is accompanied by interactions involving other higher spin states well below the Planck scale $M_{\rm pl}$. These higher spin states can come from either the QCD sector or the gravity sector, but both these resolutions have some surprising implications. For example, QCD states can resolve the problem since there is a non-trivial mixing between the QCD sector and the gravity sector, requiring all particles to interact with glueballs at tree-level. If gravity sector states restore causality, any weakly coupled UV completion of the gravity-sector must have many stringy features, with an upper bound on the string scale. Under the assumption that gravity is weakly coupled, both scenarios imply that the theory has a stringy description above $N\gtrsim \frac{M_{\rm pl}}{\Lambda_{\rm QCD}}$, where $\Lambda_{\rm QCD}$ is the confinement scale.


title: MOA-2009-BLG-319Lb: A Sub-Saturn Planet Inside the Predicted Mass Desert

abstract: We present an adaptive optics (AO) analysis of images from the Keck-II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The $\sim$10 year baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of $66.5\pm 1.7\,$mas from the source star, consistent with the light curve model prediction. The combination of the host star brightness and light curve parameters yield host star and planet masses of M_host = 0.514 $\pm$ 0.063M_Sun and m_p = 66.0 $\pm$ 8.1M_Earth at a distance of $D_L = 7.0 \pm 0.7\,$kpc. The star-planet projected separation is $2.03 \pm 0.21\,$AU. The planet-star mass ratio of this system, $q = (3.857 \pm 0.029)\times 10^{-4}$, places it in the predicted "planet desert" at $10^{-4} < q < 4\times 10^{-4}$ according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. (2016) sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 $\leq$ M_host/M_Sun $\leq$ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M_Sun. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts.


title: Understanding Q-Balls Beyond the Thin-Wall Limit

abstract: Complex scalar fields charged under a global U(1) symmetry can admit non-topological soliton configurations called Q-balls which are stable against decay into individual particles or smaller Q-balls. These Q-balls are interesting objects within quantum field theory, but are also of phenomenological interest in several cosmological and astrophysical contexts. The Q-ball profiles are determined by a nonlinear differential equation, and so generally require solution by numerical methods. In this work, we derive analytical approximations for the Q-ball profile in a polynomial potential and obtain simple expressions for the important Q-ball properties of charge, energy, and radius. These results improve significantly on the often-used thin-wall approximation and make it possible to describe Q-balls to excellent precision without having to solve the underlying differential equation.


title: The spatially resolved broad line region of IRAS 09149-6206

abstract: We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$\gamma$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ~0.5 degree between the line and continuum. This represents an offset of ~120 $\mu$as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3 $\mu$m continuum. The offset is well within the dust sublimation region, which matches the measured ~0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br$\gamma$ line. We infer the radius of the BLR to be ~65 $\mu$as (0.075 pc), which is consistent with the radius-luminosity relation of nearby active galactic nuclei derived based on the time lag of the H$\beta$ line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is $\sim 1\times10^8\,M_\odot$, which is a little below, but consistent with, the standard $M_{\rm BH}$-$\sigma_*$ relation.


title: MusE GAs FLOw and Wind (MEGAFLOW) V. The dust/metallicity-anisotropy of the Circum-Galactic Medium

abstract: We investigate whether the dust content of the circum-galactic medium (CGM) depends on the location of the quasar sightline with respect to the galaxy major-axis using 13 galaxy-MgII absorber pairs (9 - 81 kpc distance) from the MusE GAs FLOw and Wind (MEGAFLOW) survey at 0.4 < z < 1.4. The dust content of the CGM is obtained from [Zn/Fe] using Ultraviolet and Visual Echelle Spectrograph (UVES) data. When a direct measurement of [Zn/Fe] is unavailable, we estimate the dust depletion from a method which consists in solving for the depletion from multiple singly ionized ions (e.g. MnII, CrII, ZnII) since each ion depletes on dust grains at different rates. We find a positive correlation between the azimuthal angle and [Zn/Fe] with a Pearson's r = 0.70 +/- 0.14. The sightlines along the major axis show [Zn/Fe] < 0.5, whereas the [Zn/Fe] is > 0.8 along the minor axis. These results suggest that the CGM along the minor axis is on average more metal enriched (by ~ 1 dex) than the gas located along the major axis of galaxies provided that dust depletion is a proxy for metallicity. This anisotropic distribution is consistent with recent results on outflow and accretion in hydro-dynamical simulations.


title: Dying of the Light: An X-ray Fading Cold Quasar at z ~ 0.405

abstract: Cold quasars are a rare subpopulation observed to host unobscured, X-ray luminous active galactic nuclei (AGN) while also retaining a cold gas supply fueling high star formation rates. These objects are interpreted as AGN early in their evolution. We present new SOFIA HAWC+ far-infrared observations, FUV-FIR photometry, and optical spectroscopy to characterize the accretion and star formation behavior in a cold quasar at z ~ 0.405 (CQ 4479). CQ 4479 is a starburst galaxy with a predominantly young stellar population and a high gas mass fraction of ~50-70%. The AGN component has yet to become the dominant component of the FIR emission. We also find AGN bolometric luminosity that varies as a function of observation method and AGN region probed. Finally, we identify a candidate outflow feature corroborating the hypothesis that cold quasars have energetic feedback. This object presents an intriguing look into the early stages of AGN feedback and probes the rare phase where an AGN and cold gaseous component co-exist.


title: Quantized surface magnetism and higher-order topology: Application to the Hopf insulator

abstract: We identify topological aspects of the subextensive magnetic moment contributed by the surfaces of a three-dimensional crystallite -- assumed to be insulating in the bulk as well as on all surface facets, with trivial Chern invariants in the bulk. The geometric component of this subextensive moment is given by its derivative with respect to the chemical potential, at zero temperature and zero field, per unit surface area, and hence corresponds to the surface magnetic compressibility. The sum of the surface compressibilities contributed by two opposite facets of a cube-shaped crystallite is quantized to an integer multiple of the fundamental constant $e/h c$; this integer is in one-to-one correspondence with the net chirality of hinge modes on the surface of the crystallite, manifesting a link with higher-order topology. The contribution by a single facet to the magnetic compressibility need not be quantized to integers; however, symmetry and/or Hilbert-space constraints can fix the single-facet compressibility to half-integer multiples of $e/hc$, as will be exemplified by the Hopf insulator.


title: Waveform systematics in the gravitational-wave inference of tidal parameters and equation of state from binary neutron star signals

abstract: Gravitational-wave signals from binary neutron star coalescences carry information about the star's equation of state in their tidal signatures. A major issue in the inference of the tidal parameters (or directly of the equation of state) is the systematic error introduced by the waveform approximants. We use a bottom-up approach based on gauge-invariant phase analysis and the Fisher information matrix to investigate waveform systematics and help identifying biases in parameter estimation. A mock analysis of 15 different binaries indicates that systematics in current waveform models dominate over statistical errors at signal-to-noise ratio (SNR) ${\gtrsim} 80$. This implies biases in the inference of the reduced tidal parameter that are are larger than the statistical $90\%$ credible-intervals. For example, while the neutron-star radius could be constrained at ${\sim} 5\%$ level at SNR 80, systematics can be at the ${\sim} 10\%$ level. We apply our approach to GW170817 (SNR ${\sim}30$) and confirm that no significant systematic effects are present. Using an optimal frequency range for the analysis, we estimate a neutron-star radius of $12.5^{+1.1}_{-1.8}\,$km. The latter is consistent with an electromagnetic-informed prior and the recent NICER measurement. Exploring SNR ${\gtrsim}100$ in view of third-generation detectors, we find that all the current waveform models lead to differences of at least 1-sigma in the inference of the reduced tidal parameter (for any value of the latter). We conclude that current waveform models, including those from numerical relativity, are insufficient to infer the equation of state in the loudest (and potentially most informative) events that will be observed by advanced and third generation detectors.


title: Black Hole Mergers from Hierarchical Triples in Dense Star Clusters

abstract: Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical triples, formed via binary--binary encounters in the CMC Cluster Catalog, a suite of cluster simulations with present-day properties representative of the Milky Way's globular clusters. We compare the properties of the mergers from triples to the other merger channels in dense star clusters, and show that triple systems do not produce significant differences in terms of mass and effective spin distribution. However, they represent an important pathway for forming eccentric mergers, which could be detected by LIGO--Virgo/KAGRA (LVK), and future missions such as LISA and DECIGO. We derive a conservative lower limit for the merger rate from this channel of $0.35$ Gpc$^{-3}$yr$^{-1}$ in the local Universe and up to $\sim9\%$ of these events may have a detectable eccentricity at LVK design sensitivity. Additionally, we find that triple systems could play an important role in retaining second-generation BHs, which can later merge again in the core of the host cluster.


title: Protecting Quantum Information in Quantum Dot Spin Chains by Driving Exchange Interactions Periodically

abstract: Recent work has demonstrated a new route to discrete time crystal physics in quantum spin chains by periodically driving nearest-neighbor exchange interactions in gate-defined quantum dot arrays [arXiv:2006.10913]. Here, we present a detailed analysis of exchange-driven Floquet physics in small arrays of GaAs quantum dots, including phase diagrams and additional diagnostics. We also show that emergent time-crystalline behavior can benefit the protection and manipulation of multi-spin states. For typical levels of nuclear spin noise in GaAs, the combination of driving and interactions protects spin-singlet states beyond what is possible in the absence of exchange interactions. We further show how to construct a time-crystal-inspired CZ gate between singlet-triplet qubits with high fidelity. These results show that periodically driving exchange couplings can enhance the performance of quantum dot spin systems for quantum information applications.


title: AstroVaDEr: Astronomical Variational Deep Embedder for Unsupervised Morphological Classification of Galaxies and Synthetic Image Generation

abstract: We present AstroVaDEr, a variational autoencoder designed to perform unsupervised clustering and synthetic image generation using astronomical imaging catalogues. The model is a convolutional neural network that learns to embed images into a low dimensional latent space, and simultaneously optimises a Gaussian Mixture Model (GMM) on the embedded vectors to cluster the training data. By utilising variational inference, we are able to use the learned GMM as a statistical prior on the latent space to facilitate random sampling and generation of synthetic images. We demonstrate AstroVaDEr's capabilities by training it on gray-scaled \textit{gri} images from the Sloan Digital Sky Survey, using a sample of galaxies that are classified by Galaxy Zoo 2. An unsupervised clustering model is found which separates galaxies based on learned morphological features such as axis ratio, surface brightness profile, orientation and the presence of companions. We use the learned mixture model to generate synthetic images of galaxies based on the morphological profiles of the Gaussian components. AstroVaDEr succeeds in producing a morphological classification scheme from unlabelled data, but unexpectedly places high importance on the presence of companion objects---demonstrating the importance of human interpretation. The network is scalable and flexible, allowing for larger datasets to be classified, or different kinds of imaging data. We also demonstrate the generative properties of the model, which allow for realistic synthetic images of galaxies to be sampled from the learned classification scheme. These can be used to create synthetic image catalogs or to perform image processing tasks such as deblending.


title: A Galactic Dust Devil: far-infrared observations of the Tornado Supernova Remnant candidate

abstract: We present complicated dust structures within multiple regions of the candidate supernova remnant (SNR) the `Tornado' (G357.7-0.1) using observations with Spitzer and Herschel. We use Point Process Mapping, PPMAP, to investigate the distribution of dust in the Tornado at a resolution of 8", compared to the native telescope beams of 5-36". We find complex dust structures at multiple temperatures within both the head and the tail of the Tornado, ranging from 15 to 60K. Cool dust in the head forms a shell, with some overlap with the radio emission, which envelopes warm dust at the X-ray peak. Akin to the terrestrial sandy whirlwinds known as `Dust Devils', we find a large mass of dust contained within the Tornado. We derive a total dust mass for the Tornado head of 16.7 solar masses, assuming a dust absorption coefficient of kappa_300 =0.56m^2 kg^1, which can be explained by interstellar material swept up by a SNR expanding in a dense region. The X-ray, infra-red, and radio emission from the Tornado head indicate that this is a SNR. The origin of the tail is more unclear, although we propose that there is an X-ray binary embedded in the SNR, the outflow from which drives into the SNR shell. This interaction forms the helical tail structure in a similar manner to that of the SNR W50 and microquasar SS433.


title: ACCESS: Confirmation of no potassium in the atmosphere of WASP-31b

abstract: We present a new optical (400-950nm) transmission spectrum of the hot Jupiter WASP-31b (M=0.48 MJ; R= 1.54 RJ; P=3.41 days), obtained by combining four transits observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na I and K I, with particular focus on K I, for which there have been two previously published disagreeing results. Observations with HST/STIS detected K I, whereas ground-based low- and high-resolution observations did not. We use equilibrium and non-equilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (alpha=5.3+2.9-3.1), high-altitude clouds at a log cloud top pressure of -3.6+2.7-2.1 bars, and possible muted H2O features. We find that our observations support other ground-based claims of no K I. Clouds are likely why signals like H2O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (alpha = 8.3+/-1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.


title: Towards simulating a realistic data analysis with an optimised angular power spectrum of spectroscopic galaxy surveys

abstract: The angular power spectrum is a natural tool to analyse the observed galaxy number count fluctuations. In a standard analysis, the angular galaxy distribution is sliced into concentric redshift bins and all correlations of its harmonic coefficients between bin pairs are considered---a procedure referred to as `tomography'. However, the unparalleled quality of data from oncoming spectroscopic galaxy surveys for cosmology will render this method computationally unfeasible, given the increasing number of bins. Here, we put to test against synthetic data a novel method proposed in a previous study to save computational time. According to this method, the whole galaxy redshift distribution is subdivided into thick bins, neglecting the cross-bin correlations among them; each of the thick bin is, however, further subdivided into thinner bins, considering in this case all the cross-bin correlations. We create a simulated data set that we then analyse in a Bayesian framework. We confirm that the newly proposed method saves computational time and gives results that surpass those of the standard approach.


title: Hierarchical Multi-Grained Generative Model for Expressive Speech Synthesis

abstract: This paper proposes a hierarchical generative model with a multi-grained latent variable to synthesize expressive speech. In recent years, fine-grained latent variables are introduced into the text-to-speech synthesis that enable the fine control of the prosody and speaking styles of synthesized speech. However, the naturalness of speech degrades when these latent variables are obtained by sampling from the standard Gaussian prior. To solve this problem, we propose a novel framework for modeling the fine-grained latent variables, considering the dependence on an input text, a hierarchical linguistic structure, and a temporal structure of latent variables. This framework consists of a multi-grained variational autoencoder, a conditional prior, and a multi-level auto-regressive latent converter to obtain the different time-resolution latent variables and sample the finer-level latent variables from the coarser-level ones by taking into account the input text. Experimental results indicate an appropriate method of sampling fine-grained latent variables without the reference signal at the synthesis stage. Our proposed framework also provides the controllability of speaking style in an entire utterance.


title: Hall coefficient of semimetals

abstract: A recently developed formula for the Hall coefficient [A. Auerbach, Phys. Rev. Lett. 121, 66601 (2018)] is applied to nodal line and Weyl semimetals (including graphene), and to spin-orbit split semiconductor bands in two and three dimensions. The calculation reduces to a ratio of two equilibrium susceptibilities, where corrections are negligible at weak disorder. Deviations from Drude's inverse carrier density are associated with band degeneracies, Fermi surface topology, and interband scattering. Experiments which can measure these deviations are proposed.


title: Congruences of algebraic automorphic forms and supercuspidal representations

abstract: Let $G$ be a connected reductive group over a totally real field $F$ which is compact modulo center at archimedean places. We find congruences modulo an arbitrary power of p between the space of arbitrary automorphic forms on $G(\mathbb A_F)$ and that of automorphic forms with supercuspidal components at p, provided that p is larger than the Coxeter number of the absolute Weyl group of $G$. We illustrate how such congruences can be applied in the construction of Galois representations. Our proof is based on type theory for representations of p-adic groups, generalizing the prototypical case of GL(2) in [arXiv:1506.04022, Section 7] to general reductive groups. We exhibit a plethora of new supercuspidal types consisting of arbitrarily small compact open subgroups and characters thereof. We expect these results of independent interest to have further applications. For example, we extend the result by Emerton--Pa\v{s}k\=unas on density of supercuspidal points from definite unitary groups to general $G$ as above.


title: Illuminating the electronic properties of WS$_2$ polytypism with electron microscopy

abstract: Tailoring the specific stacking sequences (polytypes) of layered materials represents a powerful strategy to identify and design novel physical properties. While nanostructures built upon transition-metal dichalcogenides (TMDs) with either the 2H or 3R crystalline phases have been routinely studied, our knowledge of those based on mixed 2H/3R polytypes is far more limited. Here we report on the characterization of mixed 2H/3R free-standing WS$_2$ nanostructures displaying a flower-like configuration by means of advanced transmission electron microscopy. We correlate their rich variety of shape-morphology combinations with relevant local electronic properties such as their edge, surface, and bulk plasmons. Electron energy-loss spectroscopy combined with machine learning reveals that the 2H/3R polytype displays an indirect band gap with $E_{\rm BG}=1.6^{+0.3}_{-0.2}~{\rm eV}$. Further, we identify the presence of energy-gain peaks in the EEL spectra characterized by a gain-to-loss ratio $I_G/I_L>1$. Such property could be exploited to develop novel cooling strategies for atomically thin TMD nanostructures and devices built upon them. Our results represent a stepping stone towards an improved understanding of TMDs nanostructures based on mixed crystalline phases.


title: PhenoTagger: A Hybrid Method for Phenotype Concept Recognition using Human Phenotype Ontology

abstract: Automatic phenotype concept recognition from unstructured text remains a challenging task in biomedical text mining research. Previous works that address the task typically use dictionary-based matching methods, which can achieve high precision but suffer from lower recall. Recently, machine learning-based methods have been proposed to identify biomedical concepts, which can recognize more unseen concept synonyms by automatic feature learning. However, most methods require large corpora of manually annotated data for model training, which is difficult to obtain due to the high cost of human annotation. In this paper, we propose PhenoTagger, a hybrid method that combines both dictionary and machine learning-based methods to recognize Human Phenotype Ontology (HPO) concepts in unstructured biomedical text. We first use all concepts and synonyms in HPO to construct a dictionary, which is then used to automatically build a distantly supervised training dataset for machine learning. Next, a cutting-edge deep learning model is trained to classify each candidate phrase (n-gram from input sentence) into a corresponding concept label. Finally, the dictionary and machine learning-based prediction results are combined for improved performance. Our method is validated with two HPO corpora, and the results show that PhenoTagger compares favorably to previous methods. In addition, to demonstrate the generalizability of our method, we retrained PhenoTagger using the disease ontology MEDIC for disease concept recognition to investigate the effect of training on different ontologies. Experimental results on the NCBI disease corpus show that PhenoTagger without requiring manually annotated training data achieves competitive performance as compared with state-of-the-art supervised methods.


title: Deterministic Algorithms for Decremental Shortest Paths via Layered Core Decomposition

abstract: In the decremental single-source shortest paths (SSSP) problem, the input is an undirected graph $G=(V,E)$ with $n$ vertices and $m$ edges undergoing edge deletions, together with a fixed source vertex $s\in V$. The goal is to maintain a data structure that supports shortest-path queries: given a vertex $v\in V$, quickly return an (approximate) shortest path from $s$ to $v$. The decremental all-pairs shortest paths (APSP) problem is defined similarly, but now the shortest-path queries are allowed between any pair of vertices of $V$. Both problems have been studied extensively since the 80's, and algorithms with near-optimal total update time and query time have been discovered for them. Unfortunately, all these algorithms are randomized and, more importantly, they need to assume an oblivious adversary. Our first result is a deterministic algorithm for the decremental SSSP problem on weighted graphs with $O(n^{2+o(1)})$ total update time, that supports $(1+\epsilon)$-approximate shortest-path queries, with query time $O(|P|\cdot n^{o(1)})$, where $P$ is the returned path. This is the first $(1+\epsilon)$-approximation algorithm against an adaptive adversary that supports shortest-path queries in time below $O(n)$, that breaks the $O(mn)$ total update time bound of the classical algorithm of Even and Shiloah from 1981. Our second result is a deterministic algorithm for the decremental APSP problem on unweighted graphs that achieves total update time $O(n^{2.5+\delta})$, for any constant $\delta>0$, supports approximate distance queries in $O(\log\log n)$ time; the algorithm achieves an $O(1)$-multiplicative and $n^{o(1)}$-additive approximation on the path length. All previous algorithms for APSP either assume an oblivious adversary or have an $\Omega(n^{3})$ total update time when $m=\Omega(n^{2})$.


title: A Clockwork Solution to the Doublet-Triplet Splitting Problem

abstract: Maybe the biggest puzzle in grand unified theories (GUTs) is the apparent large splitting of the doublet and triplet Higgs masses. We suggest a novel mechanism to solve this puzzle, which relies on the clockwork mechanism to generate large hierarchies from order-one numbers. The tension between gauge coupling unification and proton lifetime from minimal SU(5) GUTs is also removed in this scenario, and the theory remains perturbative until the Planck scale.


title: Excitation and acceleration of molecular outflows in LIRGs: The extended ESO 320-G030 outflow on 200-pc scales

abstract: We used high-spatial resolution (70 pc; 0.3") CO multi-transition (1-0, 2-1, 4-3, and 6-5) ALMA data to study the physical conditions and kinematics of the cold molecular outflow in the local LIRG ESO320-G030 (d=48 Mpc, log LIR/Lsun=11.3). ESO320-G030 is a double-barred isolated spiral, but its compact and obscured nuclear starburst (SFR~15 Msun/yr; Av~40 mag) resembles those of more luminous ULIRGs. In the outflow, the 1-0/2-1 ratio is enhanced with respect to the rest of the galaxy and the CO(4-3) transition is undetected. This indicates that the outflowing molecular gas is less excited than the gas in the nuclear starburst (launching site) and the galaxy disk. Non-LTE radiative transfer modeling reveals that the properties of the outflow molecular clouds differ from those of the nuclear and disk clouds: The kinetic temperature is lower (~9 K) in the outflow, and the outflowing clouds have lower column densities. Assuming a 10^-4 CO abundance, the large internal velocity gradients, 60^+250_-45 km/s/pc, imply that the outflowing molecular clouds are not bound by self-gravity. All this suggests that the life-cycle (formation, collapse, dissipation) of the disk clouds might differ from that of the outflowing clouds which might not be able to form stars. The low Tkin of the molecular outflow remains constant up to 1.7 kpc. This indicates that the heating by the hotter ionized outflow phase is not efficient and may favor the survival of the outflow molecular phase. The velocity structure of the outflow shows a 0.8 km/s/pc velocity gradient between 190-560 pc and then a constant maximum velocity (~750 km/s) up to 1.7 kpc. This is compatible with a pure gravitational evolution of the outflow under certain mass outflow rate and launching velocity variations. Alternatively, ram pressure acceleration and cloud evaporation could explain the observed kinematics and size of the molecular phase.