title: Exact Computation of Maximum Rank Correlation Estimator

abstract: In this paper we provide a computation algorithm to get a global solution for the maximum rank correlation estimator using the mixed integer programming (MIP) approach. We construct a new constrained optimization problem by transforming all indicator functions into binary parameters to be estimated and show that it is equivalent to the original problem. We also consider an application of the best subset rank prediction and show that the original optimization problem can be reformulated as MIP. We derive the non-asymptotic bound for the tail probability of the predictive performance measure. We investigate the performance of the MIP algorithm by an empirical example and Monte Carlo simulations.


title: Existence, nonexistence, and asymptotic behavior of solutions for $N$-Laplacian equations involving critical exponential growth in the whole $\mathbb{R}^N$

abstract: In this paper, we are interested in studying the existence or non-existence of solutions for a class of elliptic problems involving the $N$-Laplacian operator in the whole space. The nonlinearity considered involves critical Trudinger-Moser growth. Our approach is non-variational, and in this way, we can address a wide range of problems not yet contained in the literature. Even $W^{1,N}(\mathbb{R}^N)\hookrightarrow L^\infty(\mathbb{R}^N)$ failing, we establish $\|u_{\lambda}\|_{L^\infty(\mathbb{R}^N)} \leq C \|u\|_{W^{1,N}(\mathbb{R}^N)}^{\Theta}$ (for some $\Theta>0$), when $u$ is a solution. To conclude, we explore some asymptotic properties.


title: On Architecture to Architecture Mapping for Concurrency

abstract: Mapping programs from one architecture to another plays a key role in technologies such as binary translation, decompilation, emulation, virtualization, and application migration. Although multicore architectures are ubiquitous, the state-of-the-art translation tools do not handle concurrency primitives correctly. Doing so is rather challenging because of the subtle differences in the concurrency models between architectures. In response, we address various aspects of the challenge. First, we develop correct and efficient translations between the concurrency models of two mainstream architecture families: x86 and ARM (versions 7 and 8). We develop direct mappings between x86 and ARMv8 and ARMv7, and fence elimination algorithms to eliminate redundant fences after direct mapping. Although our mapping utilizes ARMv8 as an intermediate model for mapping between x86 and ARMv7, we argue that it should not be used as an intermediate model in a decompiler because it disallows common compiler transformations. Second, we propose and implement a technique for inserting memory fences for safely migrating programs between different architectures. Our technique checks robustness against x86 and ARM, and inserts fences upon robustness violations. Our experiments demonstrate that in most of the programs both our techniques introduce significantly fewer fences compared to naive schemes for porting applications across these architectures.


title: Ultrahigh flexoelectric effect of 3D interconnected porous polymers: modelling and verification

abstract: Non-conductive materials like rubbers, plastics, ceramics, and even semiconductors have the property of flexoelectricity, which means that they can generate electricity when bent and twisted. However, an irregular shape or a peculiar load has been the necessary condition to realize flexoelectricity, and the weight and deformability specific ratios of flexoelectricity of solids are limited. In this work, we develop a theoretical model of flexoelectricity of three-dimensional interconnected porous materials. Compared to the solid materials, porous materials can exhibit flexoelectricity under arbitrary loading forms due to their complex microstructures, and the weight and deformability specific flexoelectric output is much higher than that of the solids. Then, we verify the model by measuring the flexoelectric response of polydimethylsiloxane (PDMS) and porous polyvinylidene fluoride (PVDF). The porous PDMS with 3D micron-scale interconnected structures exhibits two orders of magnitude higher weight and deformability specific flexoelectric output than that of the solid truncated pyramid PDMS. The flexoelectric signal is found to be linearly proportional to the applied strain, the microstructural size and the frequency. Finally, we apply the theory to a more practical bending sensor, and demonstrate its stable functioning and accurate response. Our model can be applied to other porous materials, and the results highlight the new potential of porous micro-structured materials with a significant flexoelectric effect in the fields of mechanical sensing, actuating, energy harvesting, and biomimetics as light-weight materials.


title: Topological metals and finite-momentum superconductors

abstract: We show that Zeeman field can induce a topological transition in two-dimensional spin-orbit coupled metals, and concomitantly, a first-order phase transition in the superconducting state involving a discontinuous change of Cooper pair momentum. Depending on the spin-orbit coupling strength, we find different phase diagrams of 2D superconductors under in-plane magnetic field.


title: LynyrdSkynyrd at WNUT-2020 Task 2: Semi-Supervised Learning for Identification of Informative COVID-19 English Tweets

abstract: We describe our system for WNUT-2020 shared task on the identification of informative COVID-19 English tweets. Our system is an ensemble of various machine learning methods, leveraging both traditional feature-based classifiers as well as recent advances in pre-trained language models that help in capturing the syntactic, semantic, and contextual features from the tweets. We further employ pseudo-labelling to incorporate the unlabelled Twitter data released on the pandemic. Our best performing model achieves an F1-score of 0.9179 on the provided validation set and 0.8805 on the blind test-set.


title: Maximizing Privacy in MIMO Cyber-Physical Systems Using the Chapman-Robbins Bound

abstract: Privacy breaches of cyber-physical systems could expose vulnerabilities to an adversary. Here, privacy leaks of step inputs to linear-time-invariant systems are mitigated through additive Gaussian noise. Fundamental lower bounds on the privacy are derived, which are based on the variance of any estimator that seeks to recreate the input. Fully private inputs are investigated and related to transmission zeros. Thereafter, a method to increase the privacy of optimal step inputs is presented and a privacy-utility trade-off bound is derived. Finally, these results are verified on data from the KTH Live-In Lab Testbed, showing good correspondence with theoretical results.


title: Referenced Thermodynamic Integration for Bayesian Model Selection: Application to COVID-19 Model Selection

abstract: Model selection is a fundamental part of the applied Bayesian statistical methodology. Metrics such as the Akaike Information Criterion are commonly used in practice to select models but do not incorporate the uncertainty of the models' parameters and can give misleading choices. One approach that uses the full posterior distribution is to compute the ratio of two models' normalising constants, known as the Bayes factor. Often in realistic problems, this involves the integration of analytically intractable, high-dimensional distributions, and therefore requires the use of stochastic methods such as thermodynamic integration (TI). In this paper we apply a variation of the TI method, referred to as referenced TI, which computes a single model's normalising constant in an efficient way by using a judiciously chosen reference density. The advantages of the approach and theoretical considerations are set out, along with explicit pedagogical 1 and 2D examples. Benchmarking is presented with comparable methods and we find favourable convergence performance. The approach is shown to be useful in practice when applied to a real problem - to perform model selection for a semi-mechanistic hierarchical Bayesian model of COVID-19 transmission in South Korea involving the integration of a 200D density.


title: Application of The Wind-Driven Model to A Sample of Tidal Disruption Events

abstract: An origin of the Optical/UV radiation from tidal disruption events (TDEs) has recently been discussed for different scenarios, but observational support is generally missing. In this Letter, we test applicability of the `Wind-Driven model' (Uno & Maeda 2020) to a sample of UV/Optical TDEs. With the model, we aim to derive the physical properties of the Optical/UV TDEs, such as mass-loss rates and characteristic radii. The model assumes optically thick continuous outflows like stellar winds, and one key question is how the wind-launched radius is connected to physical processes in TDEs. As one possibility, through a comparison between the escape velocities estimated from their black-hole masses and the wind velocities estimated from observed line widths, we propose that the outflow is launched from the self-interaction radius ($R_{\rm SI}$) where the stellar debris stretched by the tidal force intersects; we show that the escape velocities at $R_{\rm SI}$ are roughly consistent with the wind velocities. By applying the model to a sample of Optical/UV TDE candidates, we find that explosive mass ejections ($\gtrsim 10 ~M_{\odot}{\rm yr^{-1}}$) from $R_{\rm SI}$ ($\sim 10^{14}{\rm ~cm}$) can explain the observed properties of TDEs around peak luminosity. We also apply the same framework to a peculiar transient, AT2018cow. The model suggests that AT2018cow is likely a TDE induced by an intermediate-mass black hole ($M_{\rm BH} \sim 10^{4}~M_{\odot}$).


title: Novel Invariants for Almost Geodesic Mappings of the Third Type

abstract: Two kinds of invariance for geometrical objects under transformations are involved in this paper. With respect to these kinds, we obtained novel invariants for almost geodesic mappings of the third type of a non-symmetric affine connection space in this paper. Our results are presented in two sections. In the Section 3, we obtained the invariants for the equitorsion almost geodesic mappings which do not have the property of reciprocity.


title: GW190521 and the GWTC-1 Events: Implication on the Black Hole Mass Function of Coalescing Binary Black Hole Systems

abstract: With the black hole mass function (BHMF; assuming an exponential cutoff at a mass of $\sim 40\,M_\odot$) of coalescing binary black hole systems constructed with the events detected in the O1 run of the advanced LIGO/Virgo network, Liang et al.(2017) predicted that the birth of the lightest intermediate mass black holes (LIMBHs; with a final mass of $\gtrsim 100\,M_\odot$) is very likely to be caught by the advanced LIGO/Virgo detectors in their O3 run. The O1 and O2 observation run data, however, strongly favor a cutoff of the BHMF much sharper than the exponential one. In this work we show that a power-law function followed by a sudden drop at $\sim 40\,M_\odot$ by a factor of $\sim $a few tens and then a new power-law component extending to $\geq 100M_\odot$ are consistent with the O1 and O2 observation run data. With this new BHMF, quite a few LIMBH events can be detected in the O3 observation run of advanced LIGO/Virgo. The first LIMBH born in GW190521, an event detected in the early stage of the O3 run of advanced LIGO/Virgo network, provides additional motivation for our hypothesis.


title: Induction and Exploitation of Subgoal Automata for Reinforcement Learning

abstract: In this paper we present ISA, an approach for learning and exploiting subgoals in episodic reinforcement learning (RL) tasks. ISA interleaves reinforcement learning with the induction of a subgoal automaton, an automaton whose edges are labeled by the task's subgoals expressed as propositional logic formulas over a set of high-level events. A subgoal automaton also consists of two special states: a state indicating the successful completion of the task, and a state indicating that the task has finished without succeeding. A state-of-the-art inductive logic programming system is used to learn a subgoal automaton that covers the traces of high-level events observed by the RL agent. When the currently exploited automaton does not correctly recognize a trace, the automaton learner induces a new automaton that covers that trace. The interleaving process guarantees the induction of automata with the minimum number of states, and applies a symmetry breaking mechanism to shrink the search space whilst remaining complete. We evaluate ISA in several gridworld and continuous state space problems using different RL algorithms that leverage the automaton structures. We provide an in-depth empirical analysis of the automaton learning performance in terms of the traces, the symmetry breaking and specific restrictions imposed on the final learnable automaton. For each class of RL problem, we show that the learned automata can be successfully exploited to learn policies that reach the goal, achieving an average reward comparable to the case where automata are not learned but handcrafted and given beforehand.


title: Leggett-Garg tests for macrorealism: interference experiments and the simple harmonic oscillator

abstract: Leggett-Garg (LG) tests for macrorealism were originally designed to explore quantum coherence on the macroscopic scale. Interference experiments and systems modelled by harmonic oscillators provide useful examples of situations in which macroscopicity has been approached experimentally and are readily turned into LG tests for a single dichotomic variable Q. Applying this approach to the double-slit experiment in which a non-invasive measurement at the slits is included, we exhibit LG violations. We find that these violations are always accompanied by destructive interference. The converse is not true in general and we find that there are non-trivial regimes in which there is destructive interference but the two-time LG inequalities are satisfied which implies that it is in fact often possible to assign (indirectly determined) probabilities for the interferometer paths. Similar features have been observed in recent work involving a LG analysis of a Mach-Zehnder interferometer and we compare with those results. We extend the analysis to the triple-slit experiment again finding LG violations, and we also exhibit examples of some surprising relationships between LG inequalities and NSIT conditions that do not exist for dichotomic variables. For the simple harmonic oscillator, we find an analytically tractable example showing a two-time LG violation with a gaussian initial state, echoing in simpler form recent results of Bose et al (Phys. Rev. Lett. 120, 210402 (2018)).


title: Nonlinear computations in spiking neural networks through multiplicative synapses

abstract: The brain efficiently performs nonlinear computations through its intricate networks of spiking neurons, but how this is done remains elusive. While nonlinear computations can be implemented successfully in spiking neural networks, this requires supervised training and the resulting connectivity can be hard to interpret. In contrast, the required connectivity for any computation in the form of a linear dynamical system can be directly derived and understood with the spike coding network (SCN) framework. These networks also have biologically realistic activity patterns and are highly robust to cell death. Here we extend the SCN framework to directly implement any polynomial dynamical system, without the need for training. This results in networks requiring a mix of synapse types (fast, slow, and multiplicative), which we term multiplicative spike coding networks (mSCNs). Using mSCNs, we demonstrate how to directly derive the required connectivity for several nonlinear dynamical systems. We also show how to carry out higher-order polynomials with coupled networks that use only pair-wise multiplicative synapses, and provide expected numbers of connections for each synapse type. Overall, our work demonstrates a novel method for implementing nonlinear computations in spiking neural networks, while keeping the attractive features of standard SCNs (robustness, realistic activity patterns, and interpretable connectivity). Finally, we discuss the biological plausibility of our approach, and how the high accuracy and robustness of the approach may be of interest for neuromorphic computing.


title: Dynamic Max-Consensus and Size Estimation of Anonymous Multi-Agent Networks

abstract: In this paper we propose a novel consensus protocol for discrete-time multi-agent systems (MAS), which solves the dynamic consensus problem on the max value, i.e., the dynamic max-consensus problem. In the dynamic max-consensus problem to each agent is fed a an exogenous reference signal, the objective of each agent is to estimate the instantaneous and time-varying value of the maximum among the signals fed to the network, by exploiting only local and anonymous interactions among the agents. The absolute and relative tracking error of the proposed distributed control protocol is theoretically characterized and is shown to be bounded and by tuning its parameters it is possible to trade-off convergence time for steady-state error. The dynamic Max-consensus algorithm is then applied to solve the distributed size estimation problem in a dynamic setting where the size of the network is time-varying during the execution of the estimation algorithm. Numerical simulations are provided to corroborate the theoretical analysis.


title: Trajectory Based Podcast Recommendation

abstract: Podcast recommendation is a growing area of research that presents new challenges and opportunities. Individuals interact with podcasts in a way that is distinct from most other media; and primary to our concerns is distinct from music consumption. We show that successful and consistent recommendations can be made by viewing users as moving through the podcast library sequentially. Recommendations for future podcasts are then made using the trajectory taken from their sequential behavior. Our experiments provide evidence that user behavior is confined to local trends, and that listening patterns tend to be found over short sequences of similar types of shows. Ultimately, our approach gives a450%increase in effectiveness over a collaborative filtering baseline.


title: Designing Transportable Experiments

abstract: We consider the problem of designing a randomized experiment on a source population to estimate the Average Treatment Effect (ATE) on a target population. We propose a novel approach which explicitly considers the target when designing the experiment on the source. Under the covariate shift assumption, we design an unbiased importance-weighted estimator for the target population's ATE. To reduce the variance of our estimator, we design a covariate balance condition (Target Balance) between the treatment and control groups based on the target population. We show that Target Balance achieves a higher variance reduction asymptotically than methods that do not consider the target population during the design phase. Our experiments illustrate that Target Balance reduces the variance even for small sample sizes.


title: On the Nature of the T Tauri Triple System

abstract: We present a multi-wavelength analysis to reveal the nature of the enigmatic T Tauri triple star system. New optical and infrared measurements are coupled with archival X-ray, UV and mm datasets to show morphologies of disk material and outflow kinematics. A dark lane of obscuring material is seen in silhouette in several emission lines and in model-subtracted ALMA mm continuum dust residuals near the position of T Tau Sa+Sb, revealing the attenuating circumbinary ring around T Tau S. The flux variability of T Tau S is linked in part to the binary orbit; T Tau Sb brightens near orbital apastron as it emerges from behind circumbinary material. Outflow diagnostics confirm that T Tau N powers the blue-shifted western outflow, and the T Tau S binary drives the northwest-southeastern flow. Analysis of the southern outflow shows periodic arcs ejected from the T Tau system. Correlation of these arc locations and tangential kinematics with the orbit timing suggests that launch of the last four southern outflow ejections is contemporaneous with, and perhaps triggered by, the T Tau Sa+Sb binary periastron passage. We present a geometry of the T Tau triple that has the southern components foreground to T Tau N, obscured by a circumbinary ring, with mis-aligned disks and interacting outflows. Particularly, a wind from T Tauri Sa that is perpendicular to its circumstellar disk might interact with the circumbinary material, which may explain conflicting high contrast measurements of the system outflows in the literature. T Tauri is an important laboratory to understand early dynamical processes in young multiple systems. We discuss the historical and future characteristics of the system in this context.


title: A brief review on cosmological analysis of galaxy surveys with multiple tracers

abstract: Galaxy redshift surveys are one of the key probes in modern cosmology. In the data analysis of galaxy surveys, the precision of the statistical measurement is primarily limited by the cosmic variance on large scales. Fortunately, this limitation can in principle be evaded by observing multiple types of biased tracers. In this brief review, we present the idea of the multi-tracer method, outline key steps in the data analysis, and show several worked examples based on the GAMA, BOSS and eBOSS galaxy surveys.


title: TanhSoft -- a family of activation functions combining Tanh and Softplus

abstract: Deep learning at its core, contains functions that are composition of a linear transformation with a non-linear function known as activation function. In past few years, there is an increasing interest in construction of novel activation functions resulting in better learning. In this work, we propose a family of novel activation functions, namely TanhSoft, with four undetermined hyper-parameters of the form tanh({\alpha}x+{\beta}e^{{\gamma}x})ln({\delta}+e^x) and tune these hyper-parameters to obtain activation functions which are shown to outperform several well known activation functions. For instance, replacing ReLU with xtanh(0.6e^x)improves top-1 classification accuracy on CIFAR-10 by 0.46% for DenseNet-169 and 0.7% for Inception-v3 while with tanh(0.87x)ln(1 +e^x) top-1 classification accuracy on CIFAR-100 improves by 1.24% for DenseNet-169 and 2.57% for SimpleNet model.


title: Contraction $\mathcal{L}_1$-Adaptive Control using Gaussian Processes

abstract: We present $\mathcal{CL}_1$-$\mathcal{GP}$, a control framework that enables safe simultaneous learning and control for systems subject to uncertainties. The two main constituents are contraction theory-based $\mathcal{L}_1$ ($\mathcal{CL}_1$) control and Bayesian learning in the form of Gaussian process (GP) regression. The $\mathcal{CL}_1$ controller ensures that control objectives are met while providing safety certificates. Furthermore, $\mathcal{CL}_1$-$\mathcal{GP}$ incorporates any available data into a GP model of uncertainties, which improves performance and enables the motion planner to achieve optimality safely. This way, the safe operation of the system is always guaranteed, even during the learning transients. We provide a few illustrative examples for the safe learning and control of planar quadrotor systems in a variety of environments.


title: Quasi-isometry invariants of weakly special square complexes

abstract: We define the intersection complex for the universal cover of a compact weakly special square complex and show that it is a quasi-isometry invariant. By using this quasi-isometry invariant, we study the quasi-isometric classification of 2-dimensional right-angled Artin groups and planar graph 2-braid groups. Our results cover two well-known cases of 2-dimensional right-angled Artin groups: (1) those whose defining graphs are trees and (2) those whose outer automorphism groups are finite. Finally, we show that there are infinitely many graph 2-braid groups which are quasi-isometric to right-angled Artin groups and infinitely many which are not.


title: Revisiting the Cosmological Constant Problem within Quantum Cosmology

abstract: A new perspective on the Cosmological Constant Problem (CCP) is proposed and discussed within the multiverse approach of Quantum Cosmology. It is assumed that each member of the ensemble of universes has a characteristic scale $a$ that can be used as integration variable in the partition function. An averaged characteristic scale of the ensemble is estimated by using only members that satisfy the Einstein field equations. The averaged characteristic scale is compatible with the Planck length when considering an ensemble of solutions to the Einstein field equations with an effective cosmological constant. The multiverse ensemble is split in Planck-seed universes with vacuum energy density of order one; thus, $ \tilde{\Lambda}\approx 8\pi$ in Planck units and $a$-derivable universes. For~$a$-derivable universe with a characteristic scale of the order of the observed Universe $a\approx 8\times10^{60}$, the cosmological constant $\Lambda=\tilde{\Lambda}/a^{2}$ is in the range $10^{-121}$--$10^{-122}$, which is close in magnitude to the observed value $10^{-123}$. We point out that the smallness of $\Lambda$ can be viewed to be natural if its value is associated with the entropy of the Universe. This approach to the CCP reconciles the Planck-scale huge vacuum energy--density predicted by QFT considerations, as valid for Planck-seed universes, with the observed small value of the cosmological constant as relevant to an $a$-derivable universe as~observed.


title: Resolvent estimates and resonance free domains for Schr\"odinger operators with matrix-valued potentials

abstract: We establish semiclassical resolvent estimates for Schr\"odinger operators with long-range matrix-valued potentials. As an application we prove resonance free domains both in trapping and non-trapping situations. Our results generalize the well-known results of \cite{Burq1, Martinez1} in the case of scalar Schr\"odinger operators.


title: Procedural Generation of STEM Quizzes

abstract: Electronic quizzes are used extensively for summative and formative assessment. Current Learning Management Systems (LMS) allow instructors to create quizzes through a Graphical User Interface. Despite having a smooth learning curve, question generation/editing process with such interfaces is often slow and the creation of question variants is mostly limited to random parameters. In this paper we argue that procedural question generation greatly facilitates the task of creating varied, formative, up-to-date, adaptive question banks for STEM quizzes. We present and evaluate a proof-of-concept Python API for script-based question generation, and propose different question design patterns that greatly facilitate question authoring. The API supports questions including mathematical formulas, dynamically generated images and videos, as well as interactive content such as 3D model viewers. Output questions can be imported in major LMS. For basic usage, the required programming skills are minimal. More advanced uses do require some programming knowledge, but at a level that is common in STEM instructors. A side advantage of our system is that the question bank is actually embedded in Python code, making collaboration, version control, and maintenance tasks very easy. We demonstrate the benefits of script-based generation over traditional GUI-based approaches, in terms of question richness, authoring speed and content re-usability.


title: Free Convolution Powers via Roots of Polynomials

abstract: Let $\mu$ be a compactly supported probability measure on the real line. Bercovici-Voiculescu and Nica-Speicher proved the existence of a free convolution power $\mu^{\boxplus k}$ for any real $k \geq 1$. The purpose of this short note is to give an elementary description of $\mu^{\boxplus k}$ in terms of of polynomials and roots of their derivatives. This bridge allows us to switch back and forth between free probability and the asymptotic behavior of polynomials.


title: Production of massive bosons from the decay of a massless particle beam

abstract: Taking a two interacting scalar toy model with interaction term $g\phi\chi^2$, we study the production of $\chi$ particles coming from the decay of an asymptotic and highly occupied beam of $\phi$ particles. We perform a nonperturbative analysis coming from parametric resonant instabilities and investigate the possibility that massive $\chi$ particles are produced from decays of massless $\phi$ particles from the beam. Although this process is not present in a perturbative analysis, our nonperturbative approach allows it to happen under certain conditions. For a momentum $p$ of the beam particles and a mass $m_\chi$ of the produced ones, we find that the decay is allowed if the energy density of the beam exceeds the instability threshold $p^2 m_\chi^4/(2g^2)$. We also provide an analytical expression for the spontaneous decay rate at the earliest time.


title: Intraoperative Liver Surface Completion with Graph Convolutional VAE

abstract: In this work we propose a method based on geometric deep learning to predict the complete surface of the liver, given a partial point cloud of the organ obtained during the surgical laparoscopic procedure. We introduce a new data augmentation technique that randomly perturbs shapes in their frequency domain to compensate the limited size of our dataset. The core of our method is a variational autoencoder (VAE) that is trained to learn a latent space for complete shapes of the liver. At inference time, the generative part of the model is embedded in an optimisation procedure where the latent representation is iteratively updated to generate a model that matches the intraoperative partial point cloud. The effect of this optimisation is a progressive non-rigid deformation of the initially generated shape. Our method is qualitatively evaluated on real data and quantitatively evaluated on synthetic data. We compared with a state-of-the-art rigid registration algorithm, that our method outperformed in visible areas.


title: A KELT-TESS Eclipsing Binary in a Young Triple System Associated with a "Stellar String" Theia 301

abstract: HD 54236 is a nearby, wide common-proper-motion visual pair that has been previously identified as likely being very young by virtue of strong X-ray emission and lithium absorption. Here we report the discovery that the brighter member of the wide pair, HD~54236A, is itself an eclipsing binary (EB), comprising two near-equal solar-mass stars on a 2.4 d orbit. It represents a potentially valuable opportunity to expand the number of benchmark-grade EBs at young stellar ages. Using new observations of Ca2H&K emission and lithium absorption in the wide K-dwarf companion, HD 54236B, we obtain a robust age estimate of 225 +/- 50 Myr for the system. This age estimate and Gaia proper motions show HD 54236 is associated with Theia~301, a newly discovered local "stellar string", which itself may be related to the AB Dor moving group through shared stellar members. Applying this age estimate to AB~Dor itself alleviates reported tension between observation and theory that arises for the luminosity of the 90M_Jup star/brown dwarf AB Dor C when younger age estimates are used.


title: Machine Intelligence for Outcome Predictions of Trauma Patients During Emergency Department Care

abstract: Trauma mortality results from a multitude of non-linear dependent risk factors including patient demographics, injury characteristics, medical care provided, and characteristics of medical facilities; yet traditional approach attempted to capture these relationships using rigid regression models. We hypothesized that a transfer learning based machine learning algorithm could deeply understand a trauma patient's condition and accurately identify individuals at high risk for mortality without relying on restrictive regression model criteria. Anonymous patient visit data were obtained from years 2007-2014 of the National Trauma Data Bank. Patients with incomplete vitals, unknown outcome, or missing demographics data were excluded. All patient visits occurred in U.S. hospitals, and of the 2,007,485 encounters that were retrospectively examined, 8,198 resulted in mortality (0.4%). The machine intelligence model was evaluated on its sensitivity, specificity, positive and negative predictive value, and Matthews Correlation Coefficient. Our model achieved similar performance in age-specific comparison models and generalized well when applied to all ages simultaneously. While testing for confounding factors, we discovered that excluding fall-related injuries boosted performance for adult trauma patients; however, it reduced performance for children. The machine intelligence model described here demonstrates similar performance to contemporary machine intelligence models without requiring restrictive regression model criteria or extensive medical expertise.


title: High-Bandwidth Spatial Equalization for mmWave Massive MU-MIMO with Processing-In-Memory

abstract: All-digital basestation (BS) architectures enable superior spectral efficiency compared to hybrid solutions in massive multi-user MIMO systems. However, supporting large bandwidths with all-digital architectures at mmWave frequencies is challenging as traditional baseband processing would result in excessively high power consumption and large silicon area. The recently-proposed concept of finite-alphabet equalization is able to address both of these issues by using equalization matrices that contain low-resolution entries to lower the power and complexity of high-throughput matrix-vector products in hardware. In this paper, we explore two different finite-alphabet equalization hardware implementations that tightly integrate the memory and processing elements: (i) a parallel array of multiply-accumulate (MAC) units and (ii) a bit-serial processing-in-memory (PIM) architecture. Our all-digital VLSI implementation results in 28nm CMOS show that the bit-serial PIM architecture reduces the area and power consumption up to a factor of 2x and 3x, respectively, when compared to a parallel MAC array that operates at the same throughput.


title: A Structure Theorem on Intersections of General Doubling Measures and Its Applications

abstract: We unite two themes in dyadic analysis and number theory by studying an analogue of the failure of the Hasse principle in harmonic analysis. Explicitly, we construct an explicit family of measures on the real line that are $p$-adic and $q$-adic doubling for any distinct primes $p$ and $q$, yet not doubling, and we apply these results to show analogous statements about the reverse H\"older and Muckenhoupt $A_p$ classes of weights. The proofs involve a delicate interplay among several geometric and number theoretic properties.


title: Medium induced QCD cascades: broadening and rescattering during branching

abstract: We study evolution equations describing jet propagation through quark--gluon plasma (QGP). In particular we investigate the contribution of momentum transfer during branching and find that such a contribution is sizeable. Furthermore, we study various approximations, such as the Gaussian approximation and the diffusive approximation to the jet-broadening term. We notice that in order to reproduce the BDIM equation (without the momentum transfer in the branching) the diffusive approximation requires a very large value of the jet-quenching parameter $\hat q$.


title: Modified lepton couplings and the Cabibbo-angle anomaly

abstract: Significant discrepancies between the different determinations of the Cabibbo angle have been observed. Here, we point out that this "Cabibbo-angle anomaly" can be explained by lepton flavour universality (LFU) violating New Physics (NP) in the neutrino sector. However, modified neutrino couplings to Standard Model gauge bosons also affect many other observables sensitive to LFU violation, which have to be taken into account in order to assess the viability of this explanation. Therefore, we perform a model-independent Bayesian global analysis and find that non-zero modifications of electron and muon neutrino couplings are preferred at more than $99.99\%\,{\rm C.L.}$ (corresponding to more than $4\,\sigma$). Our results show that constructive effects in the muon sector are necessary, meaning simple models with right-handed neutrinos are discarded and more sophisticated NP models required.


title: Convolution Neural Networks for diagnosing colon and lung cancer histopathological images

abstract: Lung and Colon cancer are one of the leading causes of mortality and morbidity in adults. Histopathological diagnosis is one of the key components to discern cancer type. The aim of the present research is to propose a computer aided diagnosis system for diagnosing squamous cell carcinomas and adenocarcinomas of lung as well as adenocarcinomas of colon using convolutional neural networks by evaluating the digital pathology images for these cancers. Hereby, rendering artificial intelligence as useful technology in the near future. A total of 2500 digital images were acquired from LC25000 dataset containing 5000 images for each class. A shallow neural network architecture was used classify the histopathological slides into squamous cell carcinomas, adenocarcinomas and benign for the lung. Similar model was used to classify adenocarcinomas and benign for colon. The diagnostic accuracy of more than 97% and 96% was recorded for lung and colon respectively.


title: Robust Bayesianism and Likelihoodism

abstract: We defend a new theory of statistical evidence, which we call Robust Bayesianism (RB). We prove that, under widely accepted assumptions, RB entails the law of likelihood [Royall, 1997], the likelihood principle [Berger and Wolpert, 1988], and a variety of other widely-accepted "statistical principles", e.g., the sufficiency principle [Birnbaum, 1962, 1972] and stopping-rule principle [Berger and Wolpert, 1988]. The main technical contribution of this paper is to extend some of those results to a qualitative framework in which experimenters are justified only in making comparative, non-numerical judgments of the form "A given B is more likely than C given D."


title: A fast algorithm for quadratic resource allocation problems with nested constraints

abstract: We study the quadratic resource allocation problem and its variant with lower and upper constraints on nested sums of variables. This problem occurs in many applications, in particular battery scheduling within decentralized energy management (DEM) for smart grids. We present an algorithm for this problem that runs in $O(n \log n)$ time and, in contrast to existing algorithms for this problem, achieves this time complexity using relatively simple and easy-to-implement subroutines and data structures. This makes our algorithm very attractive for real-life adaptation and implementation. Numerical comparisons of our algorithm with a subroutine for battery scheduling within an existing tool for DEM research indicates that our algorithm significantly reduces the overall execution time of the DEM system, especially when the battery is expected to be completely full or empty multiple times in the optimal schedule. Moreover, computational experiments with synthetic data show that our algorithm outperforms the currently most efficient algorithm by more than one order of magnitude. In particular, our algorithm is able to solves all considered instances with up to one million variables in less than 17 seconds on a personal computer.


title: Identification and characterization of current sheets in collisionless plasma turbulence

abstract: The properties of current sheets forming in a ion-kinetically turbulent collisionless plasma are investigated by utilizing the results of two-dimensional hybrid-kinetic numerical simulations. For this sake the algorithm proposed by Zhdankin et al. (2013) for the analysis of current sheets forming in MHD-turbulent plasmas, was extended to analyse the role and propertes of current sheets formating in a much noisier kinetically turbulent plasma. The applicability of this approach to the analysis of kinetically-turbulent plasmas is verified. Invesigated are, e.g., the effects of the choice of parameters on the current sheet recognition, viz. the threshold current density, the minimum current density and of the local regions around current density peaks. The main current sheet properties are derived, their peak current density, the peak current carrier velocity (mainly electrons), the thickness and length of the current sheets, i.e. also their aspect ratio (length/thickness). By varying the grid resolution of the simulations it is shown that, as long as the electron inertia is not taken into account, the current sheets thin down well below ion inertial length scale until numerical (grid-resolution based) dissipation stops any the further thinning.


title: Recovering Gardner Restacking with Purely Diffusive Operations

abstract: The maximum particle kinetic energy that can be extracted from an initial six-dimensional phase space distribution motivates the concept of free or available energy. The free energy depends on the allowed operations that can be performed. A key concept underlying the theoretical treatment of plasmas is the Gardner free energy, where the exchange of the contents of equal phase volumes is allowed. A second free energy concept is the diffusive free energy, in which the contents of volumes are instead averaged. For any finite discretization of phase space, the diffusive free energy is known to be less than the Gardner free energy. However, despite the apparent fundamental differences between these free energies, it is demonstrated here that the Gardner free energy may be recovered from the continuous limit of the diffusive free energy, leading to the surprise that macroscopic phase-space conservation can be achieved by ostensibly entropy-producing microscopic operations.


title: Massive Celestial Fermions

abstract: In an effort to further the study of amplitudes in the celestial CFT (CCFT), we construct conformal primary wavefunctions for massive fermions. Upon explicitly calculating the wavefunctions for Dirac fermions, we deduce the corresponding transformation of momentum space amplitudes to celestial amplitudes. The shadow wavefunctions are shown to have opposite spin and conformal dimension $2-\Delta$. The Dirac conformal primary wavefunctions are delta function normalizable with respect to the Dirac inner product provided they lie on the principal series with conformal dimension $\Delta = 1+i\lambda$ for $\lambda\in\mathbb{R}$. It is shown that there are two choices of a complete basis: single spin $J=\frac{1}{2}$ or $J=-\frac{1}{2}$ and $\lambda\in\mathbb{R}$ or multiple spin $J=\pm\frac{1}{2}$ and $\lambda\in\mathbb{R}_{+\cup 0}$. The massless limit of the Dirac conformal primary wavefunctions is shown to agree with previous literature. The momentum generators on the celestial sphere are derived and, along with the Lorentz generators, form a representation of the Poincare algebra. Finally, we show that the massive spin-$1$ conformal primary wavefunctions can be constructed from the Dirac conformal primary wavefunctions using the standard Clebsch-Gordan coefficients. We use this procedure to write the massive spin-$\frac{3}{2}$, Rarita-Schwinger, conformal primary wavefunctions. This provides a prescription for constructing all massive fermionic and bosonic conformal primary wavefunctions starting from spin-$\frac{1}{2}$.


title: Edge Selection in Bilinear Dynamical Networks

abstract: We develop some basic principles for the design and robustness analysis of a continuous-time bilinear dynamical network, where an attacker can manipulate the strength of the interconnections/edges between some of the agents/nodes. We formulate the edge protection optimization problem of picking a limited number of attack-free edges and minimizing the impact of the attack over the bilinear dynamical network.In particular, the $\mathcal{H}_2$-norm of bilinear systems is known to capture robustness and performance properties analogous to its linear counterpart and provides valuable insights for identifying which edges are most sensitive to attacks. The exact optimization problem is combinatorial in the number of edges, and brute-force approaches show poor scalability. However, we show that the $\mathcal{H}_2$-norm as a cost function is supermodular and, therefore, allows for efficient greedy approximations of the optimal solution. We illustrate and compare the effectiveness of our theoretical findings via numerical simulations.


title: Physical Conditions of Five O VI Absorption Systems Towards PG $1522+101$

abstract: We present the analysis of five O VI absorbers identified across a redshift path of z $\sim (0.6 - 1.3)$ towards the background quasar PG $1522+101$ with information on five consecutive ionization stages of oxygen from O II to O VI. The combined $HST$ and $Keck$ spectra cover UV, redshifted EUV, and optical transitions from a multitude of ions spanning ionization energies in the range of $\sim (13 - 300)$ eV. Low ionization (C II, O II, Si II, Mg II) and very high ionization species (Ne VIII, Mg X) are non-detections in all the absorbers. Three of the absorbers have coverage of He I, in one of which it is a $> 3 \sigma$ detection. The kinematic structures of these absorbers are extracted from C IV detected in $HIRES$ spectra. The farthest absorber in our sample also contains the detections of Ne V and Ne VI. Assuming co-spatial absorbing components, the ionization models show the medium to be multiphased with small-scale density-temperature inhomogeneities that are sometimes kinematically unresolved. In two of the absorbers, there is an explicit indication of the presence of a warm gas phase ($T \gtrsim 10^5$ K) traced by O VI. In the remaining absorbers, the column densities of the ions are consistent with a non-uniform photoionized medium. The sub-solar [C/O] relative abundances inferred for the absorbers point at enrichment from massive Type II supernovae. Despite metal enrichment, the inferred wide range for [O/H] $\sim$ [$-2.1, +0.2$] amongst the absorbers along with their anti-correlation with the observed H I suggest poor small-scale mixing of metals with hydrogen in the regions surrounding galaxies and the IGM.


title: The COMBS Survey -- II. Distinguishing the Metal-Poor Bulge from the Halo Interlopers

abstract: The metal-poor stars in the bulge are important relics of the Milky Way's formation history, as simulations predict that they are some of the oldest stars in the Galaxy. In order to determine if they are truly ancient stars, we must understand the origins of this population. Currently, it is unclear if the metal-poor stars in the bulge ([Fe/H] < -1 dex) are merely halo interlopers, a unique accreted population, part of the boxy/peanut-shaped (B/P) bulge or a classical bulge population. In this work, we use spectra from the VLT/FLAMES spectrograph to obtain metallicity estimates using the Ca-II triplet (CaT) of 473 bulge stars (187 of which have [Fe/H]<-1 dex), targeted using SkyMapper photometry. We also use Gaia DR2 parallaxes and proper motions to infer the Galactic positions and velocities along with orbital properties for 523 bulge stars. We employ a probabilistic orbit analysis and find that about half of our sample has a > 50\% probability of being bound to the bulge, and half are halo interlopers. We also see that the occurrence rate of halo interlopers increases steadily with decreasing metallicity across the full range of our sample (-3 < [Fe/H] < 0.5). Our examination of the kinematics of the confined compared to the unbound stars indicates the metal-poor bulge ([Fe/H] < -1 dex) comprises at least two populations; those confined to the boxy/peanut bulge and halo stars passing through the inner galaxy. We conclude that an orbital analysis approach, as we have employed, is important to uncover and understand the composite nature of the metal-poor stars in the inner region.


title: Low-Rank Training of Deep Neural Networks for Emerging Memory Technology

abstract: The recent success of neural networks for solving difficult decision tasks has incentivized incorporating smart decision making "at the edge." However, this work has traditionally focused on neural network inference, rather than training, due to memory and compute limitations, especially in emerging non-volatile memory systems, where writes are energetically costly and reduce lifespan. Yet, the ability to train at the edge is becoming increasingly important as it enables real-time adaptability to device drift and environmental variation, user customization, and federated learning across devices. In this work, we address two key challenges for training on edge devices with non-volatile memory: low write density and low auxiliary memory. We present a low-rank training scheme that addresses these challenges while maintaining computational efficiency. We then demonstrate the technique on a representative convolutional neural network across several adaptation problems, where it out-performs standard SGD both in accuracy and in number of weight writes.


title: Simulating gamma-ray production from cosmic rays interacting with the solar atmosphere in the presence of coronal magnetic fields

abstract: Cosmic rays can interact with the solar atmosphere and produce a slew of secondary messengers, making the Sun a bright gamma-ray source in the sky. Detailed observations with Fermi-LAT have shown that these interactions must be strongly affected by solar magnetic fields in order to produce the wide range of observational features, such as high flux and hard spectrum. However, the detailed mechanisms behind these features are still a mystery. In this work, we tackle this problem by performing particle interaction simulations in the solar atmosphere in the presence of coronal magnetic fields modelled using the potential field source surface (PFSS) model. We find that the low-energy (~GeV) gamma-ray production is significantly enhanced by the coronal magnetic fields, but the enhancement decreases rapidly with energy. The enhancement is directly correlated with the production of gamma rays with large deviation angles relative to the input cosmic-ray direction. We conclude that coronal magnetic fields are essential for correctly modeling solar disk gamma rays below 10GeV, but above that the effect of coronal magnetic fields diminishes. Other magnetic field structures are needed to explain the high-energy disk emission.


title: $q$-Supercongruences from the $q$-Saalsch\"{u}tz identity

abstract: In terms of the $q$-Saalsch\"{u}tz identity and the Chinese remainder theorem for coprime polynomials, we establish some $q$-supercongruences modulo the third power of a cyclotomic polynomial. In particular, we give a $q$-analogue of a formula due to Long and Ramakrishna [Adv. Math. 290 (2016), 773--808].


title: Dark Matter = Modified Gravity? Scrutinising the spacetime-matter distinction through the modified gravity/ dark matter lens

abstract: This paper scrutinises the tenability of a strict conceptual distinction between space(time) and matter via the lens of the debate between modified gravity and dark matter. In particular, we consider Berezhiani and Khoury's novel 'superfluid dark matter theory' (SFDM) as a case study. Two families of criteria for being matter and being spacetime, respectively, are extracted from the literature. Evaluation of the new scalar field postulated by SFDM according to these criteria reveals that it is as much (dark) matter as anything could possibly be, but also$-$below the critical temperature for superfluidity$-$as much (of a modification of) spacetime as anything could possibly be. A sequel paper examines possible interpretations of SFDM in light of this result, as well as the consequences for our understanding of (the importance of) the modified gravity/ dark matter distinction and the broader spacetime-matter distinction.

categories: physics.hist-ph astro-ph.CO gr-qc hep-th


title: Cartography of the space of theories: an interpretational chart for fields that are both (dark) matter and spacetime

abstract: This paper pushes back against the Democritean-Newtonian tradition of assuming a strict conceptual dichotomy between spacetime and matter. Our approach proceeds via the more narrow distinction between modified gravity/spacetime (MG) and dark matter (DM). A prequel paper argued that the novel field $\Phi$ postulated by Berezhiani and Khoury's 'superfluid dark matter theory' is as much (dark) matter as anything could possibly be, but also$-$below the critical temperature for superfluidity$-$as much (of a modification of) spacetime as anything could possibly be. Here we introduce and critically evaluate three groups of interpretations that one should consider for such Janus-faced theories. The consubstantiality interpretation holds that $\Phi$ is both (dark) matter and a modification of spacetime, analogously to the sense in which Jesus (according to catholicism) is both human and god. The fundamendalist interpretations consider for each of these roles whether they are instantiated fundamentally or emergently. The breakdown interpretations focus on the question of whether $\Phi$ signals the breakdown, in some sense to be specified, of the MG-DM dichotomy and perhaps even the broader spacetime-matter distinction. More generally, it is argued that hybrid theories urge a move towards a single space of theories, rather than two separate spaces of spacetime theories and matter theories, respectively.


title: Neural-PDE: A RNN based neural network for solving time dependent PDEs

abstract: Partial differential equations (PDEs) play a crucial role in studying a vast number of problems in science and engineering. Numerically solving nonlinear and/or high-dimensional PDEs is often a challenging task. Inspired by the traditional finite difference and finite elements methods and emerging advancements in machine learning, we propose a sequence deep learning framework called Neural-PDE, which allows to automatically learn governing rules of any time-dependent PDE system from existing data by using a bidirectional LSTM encoder, and predict the next n time steps data. One critical feature of our proposed framework is that the Neural-PDE is able to simultaneously learn and simulate the multiscale variables.We test the Neural-PDE by a range of examples from one-dimensional PDEs to a high-dimensional and nonlinear complex fluids model. The results show that the Neural-PDE is capable of learning the initial conditions, boundary conditions and differential operators without the knowledge of the specific form of a PDE system.In our experiments the Neural-PDE can efficiently extract the dynamics within 20 epochs training, and produces accurate predictions. Furthermore, unlike the traditional machine learning approaches in learning PDE such as CNN and MLP which require vast parameters for model precision, Neural-PDE shares parameters across all time steps, thus considerably reduces the computational complexity and leads to a fast learning algorithm.


title: Quantum aspects of chaos and complexity from bouncing cosmology: A study with two-mode single field squeezed state formalism

abstract: $Circuit~ Complexity$, a well known computational technique has recently become the backbone of the physics community to probe the chaotic behaviour and random quantum fluctuations of quantum fields. This paper is devoted to the study of out-of-equilibrium aspects and quantum chaos appearing in the universe from the paradigm of two well known bouncing cosmological solutions viz. $Cosine~ hyperbolic$ and $Exponential$ models of scale factors. Besides $circuit~ complexity$, we use the $Out-of-Time~ Ordered~ correlation~ (OTOC)$ functions for probing the random behaviour of the universe both at early and the late times. In particular, we use the techniques of well known two-mode squeezed state formalism in cosmological perturbation theory as a key ingredient for the purpose of our computation. To give an appropriate theoretical interpretation that is consistent with the observational perspective we use the scale factor and the number of e-foldings as a dynamical variable instead of conformal time for this computation. From this study, we found that the period of post bounce is the most interesting one. Though it may not be immediately visible, but an exponential rise can be seen in the $complexity$ once the post bounce feature is extrapolated to the present time scales. We also find within the very small acceptable error range a universal connecting relation between Complexity computed from two different kinds of cost functionals-$linearly~ weighted$ and $geodesic~ weighted$ with the OTOC. Furthermore, from the $complexity$ computation obtained from both the cosmological models and also using the well known MSS bound on quantum Lyapunov exponent, $\lambda\leq 2\pi/\beta$ for the saturation of chaos, we estimate the lower bound on the equilibrium temperature of our universe at late time scale. Finally, we provide a rough estimation of the scrambling time in terms of the conformal time.

categories: hep-th cond-mat.dis-nn gr-qc nlin.CD quant-ph


title: Solving two-dimensional non-relativistic electronic and muonic atoms governed by Chern-Simon potential

abstract: Bidimensional muonic and electronic atoms, with nuclei composed of a proton, deuteron, and triton, and governed by Chern-Simons potential, are numerically solved. Their eigenvalues and eigenfunctions are determined with a slightly modified Numerov method. Results are compared with those assuming that the same atoms are governed by the usual $1/r$ potential even in a two-dimensional space, as well as with its three-dimensional analogs.


title: Viscous evolution of a massive disk surrounding stellar-mass black holes in full general relativity

abstract: Long-term viscous neutrino-radiation hydrodynamics simulations in full general relativity are performed for a massive disk surrounding spinning stellar-mass black holes with mass $M_{\rm BH}=4$, $6$, and $10M_\odot$ and initial dimensionless spin $\chi \approx 0.8$. The initial disk is chosen to have mass $M_{\rm disk}\approx 0.1$ or $3M_\odot$ as plausible models of the remnants for the merger of black hole-neutron star binaries or the stellar core collapse from a rapidly rotating progenitor, respectively. For $M_{\rm disk} \approx 0.1M_\odot$ with the outer disk edge initially located at $r_{\rm out} \sim 200$ km, we find that $15$%-$20$% of $M_{\rm disk}$ is ejected and the average electron fraction of the ejecta is $\langle Y_e \rangle = 0.30$-$0.35$ as found in the previous study. For $M_{\rm disk} \approx 3M_\odot$, we find that $\approx 10$%-$20$% of $M_{\rm disk}$ is ejected for $r_{\rm out}\approx 200$-$1000$ km. In addition, $\langle Y_e \rangle$ of the ejecta can be enhanced to be $\gtrsim 0.4$ because the electron fraction is increased significantly during the long-term viscous expansion of the disk with high neutrino luminosity until the mass ejection sets in. Our results suggest that not heavy $r$-process elements but light trans-iron elements would be synthesized in the matter ejected from a massive torus surrounding stellar-mass black holes. We also find that the outcomes of the viscous evolution for the high-mass disk case is composed of a rapidly spinning black hole surrounded by a torus with a narrow funnel, which appears to be suitable for generating gamma-ray bursts.


title: Reducing orbital eccentricity in initial data of black hole--neutron star binaries in the puncture framework

abstract: We develop a method to compute low-eccentricity initial data of black hole--neutron star binaries in the puncture framework extending previous work on other types of compact binaries. In addition to adjusting the orbital angular velocity of the binary, the approaching velocity of a neutron star is incorporated by modifying the helical Killing vector used to derive equations of the hydrostationary equilibrium. The approaching velocity of the black hole is then induced by requiring the vanishing of the total linear momentum of the system, differently from the case of binary black holes in the puncture framework where the linear momentum of each black hole is specified explicitly. We successfully reduce the orbital eccentricity to <~0.001 by modifying the parameters iteratively using simulations of ~3 orbits both for nonprecessing and precessing configurations. We find that empirical formulas for binary black holes derived in the excision framework do not reduce the orbital eccentricity to ~0.001 for black hole--neutron star binaries in the puncture framework, although they work for binary neutron stars.


title: Quantifying the Causal Effects of Conversational Tendencies

abstract: Understanding what leads to effective conversations can aid the design of better computer-mediated communication platforms. In particular, prior observational work has sought to identify behaviors of individuals that correlate to their conversational efficiency. However, translating such correlations to causal interpretations is a necessary step in using them in a prescriptive fashion to guide better designs and policies. In this work, we formally describe the problem of drawing causal links between conversational behaviors and outcomes. We focus on the task of determining a particular type of policy for a text-based crisis counseling platform: how best to allocate counselors based on their behavioral tendencies exhibited in their past conversations. We apply arguments derived from causal inference to underline key challenges that arise in conversational settings where randomized trials are hard to implement. Finally, we show how to circumvent these inference challenges in our particular domain, and illustrate the potential benefits of an allocation policy informed by the resulting prescriptive information.


title: Structured, relativistic jets driven by radiation

abstract: Relativistic jets, or highly collimated and fast-moving outflows, are endemic to many astrophysical phenomena. The jets produced by gamma-ray bursts and tidal disruption events are accompanied by the accretion of material onto a black hole or neutron star, with the accretion rate exceeding the Eddington limit of the compact object by orders of magnitude. In such systems, radiation dominates the energy-momentum budget of the outflow, and the dynamical evolution of the jet is governed by the equations of radiation hydrodynamics. Here we show that there are analytic solutions to the equations of radiation hydrodynamics in the viscous (i.e., diffusive) regime that describe structured, relativistic jets, which consist of a fast-moving, highly relativistic core surrounded by a slower-moving, less relativistic sheath. In these solutions, the slower-moving, outer sheath contains most of the mass, and the jet structure is mediated by local anisotropies in the radiation field. We show that, depending on the pressure and density profile of the ambient medium, the angular profile of the jet Lorentz factor is Gaussian or falls off even more steeply with angle. These solutions have implications for the nature of jet production and evolution in hyperaccreting systems, and demonstrate that such jets -- and the corresponding jet structure -- can be sustained entirely by radiative processes. We discuss the implications of these findings in the context of jetted tidal disruption events and short and long gamma-ray bursts.


title: Edges and Endpoints in 21-cm Observations from Resonant Photon Production

abstract: We introduce a novel class of signatures---spectral edges and endpoints---in 21-cm measurements resulting from interactions between the standard and dark sectors. Within the context of a kinetically mixed dark photon, we demonstrate how resonant dark photon-to-photon conversions can imprint distinctive spectral features in the observed 21-cm brightness temperature, with implications for current, upcoming, and proposed experiments targeting the cosmic dawn and the dark ages. These signatures open up a qualitatively new way to look for physics beyond the Standard Model using 21-cm observations.


title: A novel framework to study the impact of binding energy distributions on the chemistry of dust grains

abstract: The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determined by evaporation, with important consequences for the formation of planets. The amount of molecules that are desorbed depends on the interaction between the species and the grain surface, which is controlled by the binding energy. Recent theoretical and experimental works point towards a distribution of values for this parameter instead of the single value often employed in astrochemical models.We present here a new "multi-binding energy" framework, to assess the effects that a distribution of binding energies has on the amount of species bound to the grains. We find that the efficiency of the surface chemistry is significantly influenced by this process with crucial consequences on the theoretical estimates of the desorbed species.


title: Which feedback mechanisms dominate in the high-pressure environment of the Central Molecular Zone?

abstract: Supernovae (SNe) dominate the energy and momentum budget of stellar feedback, but the efficiency with which they couple to the interstellar medium (ISM) depends strongly on how effectively early, pre-SN feedback clears dense gas from star-forming regions. There are observational constraints on the magnitudes and timescales of early stellar feedback in low ISM pressure environments, yet no such constraints exist for more cosmologically typical high ISM pressure environments. In this paper, we determine the mechanisms dominating the expansion of HII regions as a function of size-scale and evolutionary time within the high-pressure ($P/k_\rm{B}$~$10^{7-8}$K cm$^{-3}$) environment in the inner 100pc of the Milky Way. We calculate the thermal pressure from the warm ionised ($P_\rm{HII}$; 10$^{4}$K) gas, direct radiation pressure ($P_\rm{dir}$), and dust processed radiation pressure ($P_\rm{IR}$). We find that (1) $P_\rm{dir}$ dominates the expansion on small scales and at early times (0.01-0.1pc; $<$0.1Myr); (2) the expansion is driven by $P_\rm{HII}$ on large scales at later evolutionary stages ($>0.1$pc; $>1$Myr); (3) during the first ~1Myr of growth, but not thereafter, either $P_{\rm IR}$ or stellar wind pressure likely make a comparable contribution. Despite the high confining pressure of the environment, natal star-forming gas is efficiently cleared to radii of several pc within ~2Myr, i.e. before the first SNe explode. This `pre-processing' means that subsequent SNe will explode into low density gas, so their energy and momentum will efficiently couple to the ISM. We find the HII regions expand to a radius of 3pc, at which point they have internal pressures equal with the surrounding external pressure. A comparison with HII regions in lower pressure environments shows that the maximum size of all HII regions is set by pressure equilibrium with the ambient ISM.


title: Unboxing Quantum Black Box Models: Learning Non-Markovian Dynamics

abstract: Characterizing the memory properties of the environment has become critical for the high-fidelity control of qubits and other advanced quantum systems. However, current non-Markovian tomography techniques are either limited to discrete superoperators, or they employ machine learning methods, neither of which provide physical insight into the dynamics of the quantum system. To circumvent this limitation, we design learning architectures that explicitly encode physical constraints like the properties of completely-positive trace-preserving maps in a differential form. This method preserves the versatility of the machine learning approach without sacrificing the efficiency and fidelity of traditional parameter estimation methods. Our approach provides the physical interpretability that machine learning and opaque superoperators lack. Moreover, it is aware of the underlying continuous dynamics typically disregarded by superoperator-based tomography. This paradigm paves the way to noise-aware optimal quantum control and opens a path to exploiting the bath as a control and error mitigation resource.


title: Effective models of a semi-quark gluon plasma

abstract: In the deconfined regime of a non-Abelian gauge theory at nonzero temperature, previously it was argued that if a (gauge invariant) source is added to generate nonzero holonomy, that this source must be linear for small holonomy. The simplest example of this is the second Bernoulli polynomial. However, then there is a conundrum in computing the free energy to $\sim g^3$ in the coupling constant $g$, as part of the free energy is discontinuous as the holonomy vanishes. In this paper we investigate two ways of generating the second Bernoulli polynomial dynamically: as a mass derivative of an auxiliary field, and from two dimensional ghosts embedded isotropically in four dimensions. Computing the holonomous hard thermal loop (HHTL) in the gluon self-energy, we find that the limit of small holonomy is only well behaved for two dimensional ghosts, with a free energy which to $\sim g^3$ is continuous as the holonomy vanishes.


title: Spin-induced scalarized black holes

abstract: It was recently shown that a scalar field suitably coupled to the Gauss-Bonnet invariant $\mathcal{G}$ can undergo a spin-induced linear tachyonic instability near a Kerr black hole. This instability appears only once the dimensionless spin $j$ is sufficiently large, that is, $j \gtrsim 0.5$. A tachyonic instability is the hallmark of spontaneous scalarization. Focusing, for illustrative purposes, on a class of theories that do exhibit this instability, we show that stationary, rotating black hole solutions do indeed have scalar hair once the spin-induced instability threshold is exceeded, while black holes that lie below the threshold are described by the Kerr solution. Our results provide strong support for spin-induced black hole scalarization.


title: Spin-induced black hole scalarization in Einstein-scalar-Gauss-Bonnet theory

abstract: We construct black hole solutions with spin-induced scalarization in a class of models where a scalar field is quadratically coupled to the topological Gauss-Bonnet term. Starting from the tachyonically unstable Kerr solutions, we obtain families of scalarized black holes such that the scalar field has either even or odd parity, and we investigate their domain of existence. The scalarized black holes can violate the Kerr rotation bound. We identify "critical" families of scalarized black hole solutions such that the expansion of the metric functions and of the scalar field at the horizon no longer allows for real coefficients. For the quadratic coupling considered here, solutions with spin-induced scalarization are entropically favored over Kerr solutions with the same mass and angular momentum.


title: Deviations from the Infrared-Radio Correlation in Massive, Ultra-compact Starburst Galaxies

abstract: Feedback through energetic outflows has emerged as a key physical process responsible for transforming star-forming galaxies into the quiescent systems observed in the local universe. To explore this process, this paper focuses on a sample of massive and compact merger remnant galaxies hosting high-velocity gaseous outflows ($|v| \gtrsim 10^{3}$ km s$^{-1}$), found at intermediate redshift ($z \sim 0.6$). From their mid-infrared emission and compact morphologies, these galaxies are estimated to have exceptionally large star formation rate (SFR) surface densities ($\Sigma_{SFR} \sim 10^{3}$ $\mathrm{M_{\odot}}$ yr$^{-1}$ kpc$^{-2}$), approaching the Eddington limit for radiation pressure on dust grains. This suggests that star formation feedback may be driving the observed outflows. However, these SFR estimates suffer from significant uncertainties. We therefore sought an independent tracer of star formation to probe the compact starburst activity in these systems. In this paper, we present SFR estimates calculated using 1.5 GHz continuum Jansky Very Large Array observations for 19 of these galaxies. We also present updated infrared (IR) SFRs calculated from WISE survey data. We estimate SFRs from the IR to be larger than those from the radio for 16 out of 19 galaxies by a median factor of 2.5. We find that this deviation is maximized for the most compact galaxies hosting the youngest stellar populations, suggesting that compact starbursts deviate from the IR-radio correlation. We suggest that this deviation stems either from free-free absorption of synchrotron emission, a difference in the timescale over which each indicator traces star formation, or exceptionally hot IR-emitting dust in these ultra-dense galaxies.


title: Elliptic genera of pure gauge theories in two dimensions with semisimple non-simply-connected gauge groups

abstract: In this paper we describe a systematic method to compute elliptic genera of (2,2) supersymmetric gauge theories in two dimensions with gauge group G/Gamma (for G semisimple and simply-connected, Gamma a subgroup of the center of G) with various discrete theta angles. We apply the technique to examples of pure gauge theories with low-rank gauge groups. Our results are consistent with expectations from decomposition of two-dimensional theories with finite global one-form symmetries and with computations of supersymmetry breaking for some discrete theta angles in pure gauge theories. Finally, we make predictions for the elliptic genera of all the other remaining pure gauge theories by applying decomposition and matching to known supersymmetry breaking patterns.


title: Multi-timescale reverberation mapping of Mrk 335

abstract: Time lags due to X-ray reverberation have been detected in several Seyfert galaxies. The different travel time between reflected and directly observed rays naturally causes this type of lag, which depends directly on the light-crossing timescale of the system and hence scales with the mass of the central black hole. Featureless `hard lags' not associated with reverberation, and often interpreted as propagating mass accretion rate fluctuations, dominate the longer timescale variability. Here we fit our reltrans model simultaneously to the time-averaged energy spectrum and the lag-energy spectra of the Seyfert galaxy Mrk 335 over two timescales (Fourier frequency ranges). We model the hard lags as fluctuations in the slope and strength of the illuminating spectrum, and self-consistently account for the effects that these fluctuations have on the reverberation lags. The resulting mass estimate is $1.1^{+2.0}_{-0.7} \times 10^6~M_\odot$, which is significantly lower than the mass measured with the optical reverberation mapping technique (14 - 26 million $M_\odot$). When we add the correlated variability amplitudes to the time lags by fitting the full complex cross-spectra, the model is unable to describe the characteristic reverberation Fe K$\alpha$ line and cannot constrain the black hole mass. This may be due to the assumption that the direct radiation is emitted by a point-like source.


title: Flavorgenesis

abstract: We present a model where all fermions are contained in a single irreducible representation of an SU(19) gauge symmetry group. If there is only one scalar field, Yukawa interactions are controlled by a single number rather than by one or more $3\times3$ matrices of couplings. The low-energy concept of flavor emerges entirely from the scalar-sector parameters; more specifically, entries of the Standard Model Yukawa matrices are controlled by several vacuum expectation values.


title: Evolution of Retrograde Orbiters in an AGN Disk

abstract: AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifying the retrograde contribution is important, since roughly half of all orbiters should initially be on retrograde orbits when the disk forms. We perform an analytic calculation of the evolution of ROs in an AGN disk. Because this evolution could cause the orbits of ROs to cross those of prograde BBHs, we derive the collision rate between a given RO and a given BBH orbiting in the prograde direction. In the examples given here, ROs in the inner region of the disk experience a rapid decrease in the semimajor axis of their orbits while also becoming highly eccentric in less than a million years. This rapid orbital evolution could lead to extreme mass ratio inspirals detectable by the Laser Interferometer Space Antenna. The collision rates of our example ROs with prograde BBHs in the migration trap depend strongly on the volume of the inner radiation-pressure-dominated region which depends on the mass of the supermassive black hole (SMBH). Rates are lowest for larger mass SMBHs, which dominate the AGN merger channel, suggesting that merger rates for this channel may not be significantly altered by ROs.


title: The cosmic merger rate density of compact objects: impact of star formation, metallicity, initial mass function and binary evolution

abstract: We evaluate the redshift distribution of binary black hole (BBH), black hole - neutron star binary (BHNS) and binary neutron star (BNS) mergers, exploring the main sources of uncertainty: star formation rate (SFR) density, metallicity evolution, common envelope, mass transfer via Roche lobe overflow, natal kicks, core-collapse supernova model and initial mass function. Among binary evolution processes, uncertainties on common envelope ejection have a major impact: the local merger rate density of BNSs varies from $\sim{}10^3$ to $\sim{}20$ Gpc$^{-3}$ yr$^{-1}$ if we change the common envelope efficiency parameter from $\alpha_{\rm CE}=7$ to 0.5, while the local merger rates of BBHs and BHNSs vary by a factor of $\sim{}2-3$. The BBH merger rate changes by one order of magnitude, when $1 \sigma$ uncertainties on metallicity evolution are taken into account. In contrast, the BNS merger rate is almost insensitive to metallicity. Hence, BNSs are the ideal test bed to put constraints on uncertain binary evolution processes, such as common envelope and natal kicks. Only models assuming values of $\alpha_{\rm CE}\gtrsim{}2$ and moderately low natal kicks (depending on the ejected mass and the SN mechanism), result in a local BNS merger rate density within the 90% credible interval inferred from the second gravitational-wave transient catalogue.

categories: astro-ph.HE astro-ph.CO astro-ph.GA astro-ph.SR


title: The stellar populations of high-redshift dwarf galaxies

abstract: We use high-resolution ($\approx 10$ pc), zoom-in simulations of a typical (stellar mass $M_\star\simeq10^{10}M_\odot$) Lyman Break Galaxy (LBG) at $z\simeq 6$ to investigate the stellar populations of its six dwarf galaxy satellites, whose stellar [gas] masses are in the range $\log (M_\star/M_\odot) \simeq 6-9$ [$\log (M_{gas}/M_\odot) \simeq4.3-7.75$]. The properties and evolution of satellites show no dependence on the distance from the central massive LBG ($< 11.5$ kpc). Instead, their star formation and chemical enrichment histories are tightly connected their stellar (and sub-halo) mass. High-mass dwarf galaxies ($\rm M_\star \gtrsim 5\times 10^8 M_\odot$) experience a long history of star formation, characterised by many merger events. Lower-mass systems go through a series of short star formation episodes, with no signs of mergers; their star formation activity starts relatively late ($z\approx 7$), and it is rapidly quenched by internal stellar feedback. In spite of the different evolutionary patterns, all satellites show a spherical morphology, with ancient and more metal-poor stars located towards the inner regions. All six dwarf satellites experienced high star formation rate ($\rm >5\,M_\odot yr ^{-1}$) bursts, which can be detected by JWST while targeting high-$z$ LBGs.


title: Spectral index-flux relation for investigating the origins of steep decay in $\gamma$-ray bursts

abstract: $\gamma$-ray bursts (GRBs) are short-lived transients releasing a large amount of energy ($10^{51}-10^{53} $ erg) in the keV-MeV energy range. GRBs are thought to originate from internal dissipation of the energy carried by ultra-relativistic jets launched by the remnant of a massive star's death or a compact binary coalescence. While thousands of GRBs have been observed over the last thirty years, we still have an incomplete understanding of where and how the radiation is generated in the jet. Here we show a relation between the spectral index and the flux found by investigating the X-ray tails of bright GRB pulses via time-resolved spectral analysis. This relation is incompatible with the long standing scenario which invokes the delayed arrival of photons from high-latitude parts of the jet. While the alternative scenarios cannot be firmly excluded, the adiabatic cooling of the emitting particles is the most plausible explanation for the discovered relation, suggesting a proton-synchrotron origin of the GRB emission.


title: Control issues of KKLT

abstract: We analyze to which extent the KKLT proposal for the construction of de Sitter vacua in string theory is quantitatively controlled. Our focus is on the quality of the 10d supergravity approximation. As our main finding, we uncover and quantify an issue which one may want to call the "singular-bulk problem". In particular, we show that, requiring the curvature to be small in the conifold region, one is generically forced into a regime where the warp factor becomes negative in a significant part of the Calabi-Yau orientifold. This implies true singularities, independent of the familiar, string-theoretically controlled singularities of this type in the vicinity of O-planes. We also discuss possible escape routes as well as other control issues, related to the need for a large tadpole and hence for a complicated topology.


title: UV upturn versus UV weak galaxies: differences and similarities of their stellar populations unveiled by a de-biased sample

abstract: The ultraviolet (UV) upturn is characterised by an unexpected up-rise of the UV flux in quiescent galaxies between the Lyman limit and 2500\AA. By making use of colour-colour diagrams, one can subdivide UV bright red-sequence galaxies in two groups: UV weak and upturn. With these two groups, we propose a comparison between their stellar population properties with the goal of establishing differences and similarities between them. We make use of propensity score matching (PSM) to mitigate potential biases between the two samples, by selecting similar objects in terms of redshift and stellar mass. Also, we take advantage of spectral energy distribution (SED) fitting results from magphys made available by the GAMA collaboration. The analyses are made by comparing the distributions from the SED fitting directly, as well as investigating the differences in correlations between their parameters, and finally by using principal component analysis (PCA). We explore important differences and similarities between UV weak and upturn galaxies in terms of several parameters, such as: metallicity, age, specific star formation rate, time of last burst of star-formation, to mention a few. Notable differences are those concerning (g-r) colour, metallicity, and time since last burst of star-formation: UV upturn are redder in the optical, more metallic, and their last burst of star-formation happened earlier in time. These differences suggest that UV upturn systems have shorter star-formation histories (i.e. have been evolving more passively) when compared to UV weak galaxies. Consequently, these last seem to have a higher diversity of stellar populations.


title: A unified scenario for the origin of spiral and elliptical galaxy structural scaling laws

abstract: Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws that link their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP). Over the years, a number of studies have attempted to understand these empirical relations, usually in terms of variations of the virial theorem for E galaxies and in terms of the scaling relations of dark matter halos for spirals. We use Lambda cold dark matter (LCDM) cosmological hydrodynamical simulations to show that the scaling relations of both ellipticals and spirals arise as the result of (i) a tight galaxy mass-dark halo mass relation and (ii) the self-similar mass profile of cold dark matter halos. In this interpretation, E and S galaxies of a given stellar mass inhabit halos of similar masses, and their different scaling laws result from the varying amounts of dark matter enclosed within their luminous radii. This scenario suggests a new galaxy distance indicator applicable to galaxies of all morphologies and provides simple and intuitive explanations for long-standing puzzles, such as why the TFR is independent of surface brightness, or what causes the "tilt" in the FP. Our results provide strong support for the predictions of LCDM in the strongly non-linear regime, as well as guidance for further improvements to cosmological simulations of galaxy formation.


title: The Axion Quality Problem: Global Symmetry Breaking and Wormholes

abstract: Continuous global symmetries are expected to be broken by gravity, which can lead to important phenomenological consequences. A prime example is the threat that this poses to the viability of the Peccei-Quinn solution to the strong CP problem. In this paper, we explore the impact of wormholes as a source of global symmetry breaking by gravity. We review the current status of wormholes and global symmetries and note that, surprisingly, the axion has a quality problem within non-perturbative Einstein gravity. Although these wormholes lead to a large breaking of global symmetries, we show that their effect is nonetheless relevant for the model building of gauge protected axions. We also find wormhole solutions within two scenarios: (i) an extended global symmetry group within Einstein gravity, and (ii) U(1) wormholes within the low-energy limit of an open String Theory. The former allows us to show that the concept of a global symmetry in General Relativity is somewhat ill-defined. The latter illustrates that for motivated values of the string coupling constant, axions appear to have a quality problem within the open String Theory we consider.


title: Quantum steering with vector vortex photon states with the detection loophole closed

abstract: Violating a nonlocality inequality enables the most powerful remote quantum information tasks and fundamental tests of quantum physics. Loophole-free photonic verification of nonlocality has been achieved with polarization-entangled photon pairs, but not with states entangled in other degrees of freedom. Here we demonstrate completion of the quantum steering nonlocality task, with the detection loophole closed, when entanglement is distributed by transmitting a photon in an optical vector vortex state, formed by optical orbital angular momentum (OAM) and polarization. As well as opening up a high-efficiency encoding beyond polarization, the critically-important demonstration of vector vortex steering opens the door to new free-space and satellite-based secure quantum communication devices and device-independent protocols.


title: Constraining the entropy of formation from young transiting planets

abstract: Recently K2 and TESS have discovered transiting planets with radii between $\sim$ 5-10 R$_\oplus$ around stars with ages $<100$ Myr. These young planets are likely to be the progenitors of the ubiquitous super-earths/sub-neptunes, that are well studied around stars with ages $\gtrsim 1$ Gyr. The formation and early evolution of super-earths/sub-neptunes are poorly understood. Various planetary origin scenarios predict a wide range of possible formation entropies. We show how the formation entropies of young ($\sim$20-60 Myr), highly irradiated planets can be constrained if their mass, radius and age are measured. This method works by determining how low-mass a H/He envelope a planet can retain against mass-loss. This lower bound on the H/He envelope mass can then be converted into an upper bound on the entropy. If planet mass measurements with errors $\lesssim 20$\% can be achieved for the discovered young planets around DS Tuc A and V1298 Tau, then insights into their origins can be obtained. For these planets, higher measured planet masses would be consistent with standard core-accretion theory. In contrast, lower planet masses ($\lesssim 6-7$ M$_\oplus$) would require a "boil-off" phase during protoplanetary disc dispersal to explain.


title: An unsupervised method for identifying $X$-enriched stars directly from spectra: Li in LAMOST

abstract: Stars with peculiar element abundances are important markers of chemical enrichment mechanisms. We present a simple method, tangent space projection (TSP), for the detection of $X$-enriched stars, for arbitrary elements $X$, even from blended lines. Our method does not require stellar labels, but instead directly estimates the counterfactual unrenriched spectrum from other unlabelled spectra. As a case study, we apply this method to the $6708~$\AA{} Li doublet in LAMOST DR5, identifying 8,428 Li-enriched stars seamlessly across evolutionary state. We comment on the explanation for Li-enrichement for different subpopulations, including planet accretion, nonstandard mixing, and youth.


title: Fiber Bundle Codes: Breaking the $N^{1/2} \operatorname{polylog}(N)$ Barrier for Quantum LDPC Codes

abstract: We present a quantum LDPC code family that has distance $\Omega(N^{3/5}/\operatorname{polylog}(N))$ and $\tilde\Theta(N^{3/5})$ logical qubits. This is the first quantum LDPC code construction which achieves distance greater than $N^{1/2} \operatorname{polylog}(N)$. The construction is based on generalizing the homological product of codes to a fiber bundle.


title: Optimal portfolio choice with path dependent benchmarked labor income: a mean field model

abstract: We consider the life-cycle optimal portfolio choice problem faced by an agent receiving labor income and allocating her wealth to risky assets and a riskless bond subject to a borrowing constraint. In this paper, to reflect a realistic economic setting, we propose a model where the dynamics of the labor income has two main features. First, labor income adjust slowly to financial market shocks, a feature already considered in Biffis, E., Gozzi, F. and Prosdocimi, C. (2020) - "Optimal portfolio choice with path dependent labor income: the infinite horizon case". Second, the labor income $y_i$ of an agent $i$ is benchmarked against the labor incomes of a population $y^n:=(y_1,y_2,\ldots,y_n)$ of $n$ agents with comparable tasks and/or ranks. This last feature has not been considered yet in the literature and is faced taking the limit when $n\to +\infty$ so that the problem falls into the family of optimal control of infinite dimensional McKean-Vlasov Dynamics, which is a completely new and challenging research field. We study the problem in a simplified case where, adding a suitable new variable, we are able to find explicitly the solution of the associated HJB equation and find the optimal feedback controls.


title: Quantum many-body scars with chiral topological order in 2D and critical properties in 1D

abstract: We construct few-body, interacting, nonlocal Hamiltonians with a quantum scar state in an otherwise thermalizing many-body spectrum. In one dimension, the embedded state is a critical state, and in two dimensions, the embedded state is a chiral topologically ordered state. The models are defined on slightly disordered lattices, and the scar state appears independent of the precise realization of the disorder. A parameter allows the scar state to be placed at any position in the spectrum. We show that the level spacing distributions are Wigner-Dyson and that the entanglement entropies of the states in the middle of the spectrum are close to the Page value. Finally, we confirm the topological order in the scar state by showing that one can insert anyons into the state.


title: Variational wavefunctions for Sachdev-Ye-Kitaev models

abstract: Given a class of $q$-local Hamiltonians, is it possible to find a simple variational state whose energy is a finite fraction of the ground state energy in the thermodynamic limit? Whereas product states often provide an affirmative answer in the case of bosonic (or qubit) models, we show that Gaussian states fail dramatically in the fermionic case, like for the Sachdev-Ye-Kitaev (SYK) models. This prompts us to propose a new class of wavefunctions for SYK models inspired by the variational coupled cluster algorithm. We introduce a static ("0+0D") large-$N$ field theory to study the energy, two-point correlators, and entanglement properties of these states. Most importantly, we demonstrate a finite disorder-averaged approximation ratio of $r \approx 0.62$ between the variational and ground state energy of SYK for $q=4$. Moreover, the variational states provide an exact description of spontaneous symmetry breaking in a related two-flavor SYK model.


title: The interacting nature of dwarf galaxies hosting superluminous supernovae

abstract: (Abridged) Type I superluminous supernovae (SLSNe I) are rare, powerful explosions whose mechanism and progenitors remain elusive. SLSNe I show a preference for low-metallicity, actively star-forming dwarf galaxies. We investigate whether the hosts of SLSNe I show increased evidence for interaction. We use a sample of 42 SLSN I images obtained with $\textit{HST}$ and measure the number of companion galaxies by counting the objects detected within a given radius from the host. As a comparison, we used two Monte Carlo-based methods to estimate the expected average number of companion objects in the same images, as well as a sample of 32 galaxies that have hosted long gamma-ray bursts (GRBs). About 50% of SLSN I hosts have at least one major companion (within a flux ratio of 1:4) within 5 kpc. The average number of major companions per SLSN I host galaxy is $0.70^{+0.19}_{-0.14}$. Our Monte Carlo comparison methods yield a lower number of companions for random objects of similar brightness in the same image or for the SLSN host after randomly redistributing the sources in the same image. The Anderson-Darling test shows that this difference is statistically significant independent of the redshift range. The same is true for the projected distance distribution of the companions. The SLSN I hosts are, thus, found in areas of their images, where the object number density is greater than average. SLSN I hosts have more companions than GRB hosts ($0.44^{+0.25}_{-0.13}$ companions per host distributed over 25% of the hosts) but the difference is not statistically significant. The difference between their separations is, however, marginally significant. The dwarf galaxies hosting SLSNe I are often part of interacting systems. This suggests that SLSNe I progenitors are formed after a recent burst of star formation. Low metallicity alone cannot explain this tendency.


title: The MOSDEF Survey: An Improved Voronoi Binning Technique on Spatially Resolved Stellar Populations at z~2

abstract: We use a sample of 350 star-forming galaxies at $1.25<z<2.66$ from the MOSFIRE Deep Evolution Field survey to demonstrate an improved Voronoi binning technique that we use to study the properties of resolved stellar populations in $z\sim2$ galaxies. Stellar population and dust maps are constructed from the high-resolution CANDELS/3D-HST multi-band imaging. Rather than constructing the layout of resolved elements (i.e., Voronoi bins) from the S/N distribution of the $H_{160}$-band alone, we introduce a modified Voronoi binning method that additionally incorporates the S/N distribution of several resolved filters. The SED-derived resolved E(B-V)$_{\text{stars}}$, stellar population ages, SFRs, and stellar masses that are inferred from the Voronoi bins constructed from multiple filters are generally consistent with the properties inferred from the integrated photometry within the uncertainties, with the exception of the inferred E(B-V)$_{\text{stars}}$ from our $z\sim1.5$ sample due to their UV slopes being unconstrained by the resolved photometry. The results from our multi-filter Voronoi binning technique are compared to those derived from a "traditional" single-filter Voronoi binning approach. We find that single-filter binning produces inferred E(B-V)$_{\text{stars}}$ that are systematically redder by 0.02 mag on average, but could differ by up to 0.20 mag, and could be attributed to poorly constrained resolved photometry covering the UV slope. Overall, we advocate that our methodology produces more reliable SED-derived parameters due to the best-fit resolved SEDs being better constrained at all resolved wavelengths--particularly those covering the UV slope.