title: Single Image Dehazing Algorithm Based on Sky Region Segmentation

abstract: In this paper a hybrid image defogging approach based on region segmentation is proposed to address the dark channel priori algorithm's shortcomings in de-fogging the sky regions. The preliminary stage of the proposed approach focuses on the segmentation of sky and non-sky regions in a foggy image taking the advantageous of Meanshift and edge detection with embedded confidence. In the second stage, an improved dark channel priori algorithm is employed to defog the non-sky region. Ultimately, the sky area is processed by DehazeNet algorithm, which relies on deep learning Convolutional Neural Networks. The simulation results show that the proposed hybrid approach in this research addresses the problem of color distortion associated with sky regions in foggy images. The approach greatly improves the image quality indices including entropy information, visibility ratio of the edges, average gradient, and the saturation percentage with a very fast computation time, which is a good indication of the excellent performance of this model.


title: HSD Shared Task in VLSP Campaign 2019:Hate Speech Detection for Social Good

abstract: The paper describes the organisation of the "HateSpeech Detection" (HSD) task at the VLSP workshop 2019 on detecting the fine-grained presence of hate speech in Vietnamese textual items (i.e., messages) extracted from Facebook, which is the most popular social network site (SNS) in Vietnam. The task is organised as a multi-class classification task and based on a large-scale dataset containing 25,431 Vietnamese textual items from Facebook. The task participants were challenged to build a classification model that is capable of classifying an item to one of 3 classes, i.e., "HATE", "OFFENSIVE" and "CLEAN". HSD attracted a large number of participants and was a popular task at VLSP 2019. In particular, there were 71 teams signed up for the task, 14 of them submitted results with 380 valid submissions from 20th September 2019 to 4th October 2019.


title: Relaxation in a Phase-separating Two-dimensional Active Matter System with Alignment Interaction

abstract: Via computer simulations we study kinetics of pattern formation in a two-dimensional active matter system. Self-propulsion in our model is incorporated via the Vicsek-like activity, i.e., particles have the tendency of aligning their velocities with the average directions of motion of their neighbors. In addition to this dynamic or active interaction, there exists passive inter-particle interaction in the model for which we have chosen the standard Lennard-Jones form. Following quenches of homogeneous configurations to a point deep inside the region of coexistence between high and low density phases, as the systems exhibit formation and evolution of particle-rich clusters, we investigate properties related to the morphology, growth and aging. A focus of our study is on the understanding of the effects of structure on growth and aging. To quantify the latter we use the two-time order-parameter autocorrelation function. This correlation, as well as the growth, is observed to follow power-law time dependence, qualitatively similar to the scaling behavior reported for passive systems. The values of the exponents have been estimated and discussed by comparing with the previously obtained numbers for other dimensions as well as with the new results for the passive limit of the considered model. We have also presented results on the effects of temperature on the activity mediated phase separation.


title: The $\sigma$-meson: four-quark vs. two-quark components and decay width in a Bethe-Salpeter approach

abstract: We study the dynamical generation of resonances in isospin singlet channels with mixing between two- and four-quark states. To this end we generalise a Bethe-Salpeter approach to four-quark states employed previously \cite{Heupel:2012ua} to accommodate for mixing diagrams. The $q\bar{q}q\bar{q}$ and $q\bar{q}$ components of the Bethe-Salpeter wave function (with light quarks $q\in\{u,d\}$) are determined consistently in a symmetry-preserving truncation of the underlying Dyson-Schwinger equations. As a prominent example we deal with the isospin-singlet $0^{++}$ meson with light quark content. We find that the $\pi\pi$ contribution of the four-quark component is mainly responsible for the low (real part of the) mass of the resulting state. We also study the analytic structure in the complex momentum plane and find a branch cut at the two-pion threshold and a singularity in the second Riemann sheet indicating a considerable decay width. Our findings are in excellent qualitative agreement with the general picture for the $\sigma/f_0(500)$ that emerged in the past two decades from dispersive approaches \cite{Pelaez:2015qba}.


title: Flexible Coaxial Ribbon Cable for High-Density Superconducting Microwave Device Arrays

abstract: Superconducting electronics often require high-density microwave interconnects capable of transporting signals between temperature stages with minimal loss, cross talk, and heat conduction. We report the design and fabrication of superconducting 53 wt% Nb-47 wt% Ti (Nb47Ti) FLexible coAXial ribbon cables (FLAX). The ten traces each consist of a 0.076 mm O.D. NbTi inner conductor insulated with PFA (0.28 mm O.D.) and sheathed in a shared 0.025 mm thick Nb47Ti outer conductor. The cable is terminated with G3PO coaxial push-on connectors via stainless steel capillary tubing (1.6 mm O.D., 0.13 mm thick) soldered to a coplanar wave guide transition board. The 30 cm long cable has 1 dB of loss at 8 GHz with -60 dB nearest-neighbor forward cross talk. The loss is 0.5 dB more than commercially available superconducting coax likely due to impedance mismatches caused by manufacturing imperfections in the cable. The reported cross talk is 30 dB lower than previously developed laminated NbTi-onKapton microstrip cables. We estimate the heat load from 1 K to 90 mK to be 20 nW per trace, approximately half the computed load from the smallest commercially available superconducting coax from CryoCoax


title: Intermittency acceleration of water droplet population dynamics inside the interfacial layer between cloudy and clear air environments

abstract: We use direct numerical simulation to study the temporal evolution of a perturbation localized on the turbulent layer that typically separates a cloud from the surrounding clear air. Across this shearless layer, a turbulent kinetic energy gradient naturally forms. Here, a finite perturbation in the form of a local initial temperature fluctuation is applied to simulate a hydrodynamic instability inside the background turbulent air flow. A numerical initial value problem for two diametrically opposite types of drop population distributions is then solved. Specifically, we consider a mono-disperse population of droplets of 15 $\mu$m of radius and a poly-disperse distribution with radii in the range 0.6 - 30 $\mu$m. For both distributions, it is observed that the evaporation and condensation have a dramatically different weight inside the homogeneous cloudy region and the interfacial anisotropic mixing region. It is observed that the dynamics of drop collisions is highly effected by the turbulence structure of the host region. The two populations show a common aspect during their energy decay transient. That is the increased probability of collisions in the interfacial layer hat houses intense anisotropic velocity fluctuations. This layer, in fact, induces an enhanced differentiation on droplets kinetic energy and sizes. Both polydisperse and monodisperse initial particle distributions contain $10^7$ droplets, matching an initial liquid water content of $0.8 g/m^3$. An estimate of the turbulent collision kernel for geometric collisions used in the population balance equations is given. A preliminary discussion is presented on the structure of the two unsteady non ergodic collision kernels obtained inside the cloud interface region.


title: Wafer-scale low-loss lithium niobate photonic integrated circuits

abstract: Thin-film lithium niobate (LN) photonic integrated circuits (PICs) could enable ultrahigh performance in electro-optic and nonlinear optical devices. To date, realizations have been limited to chip-scale proof-of-concepts. Here we demonstrate monolithic LN PICs fabricated on 4- and 6-inch wafers with deep ultraviolet lithography and show smooth and uniform etching, achieving 0.27 dB/cm optical propagation loss on wafer-scale. Our results show that LN PICs are fundamentally scalable and can be highly cost-effective.


title: Effect of interactions and non-uniform magnetic states on the magnetization reversal of iron nanowire arrays

abstract: Ordered ferromagnetic nanowire arrays are widely studied due to the diversity of possible applications. However, there is still no complete understanding of the relation between the array's parameters and its magnetic behavior. The effect of vortex states on the magnetization reversal of large-diameter nanowires is of particular interest. Here, we compare analytical and micromagnetic models with experimental results for three arrays of iron nanowires with diameters of 33, 52 and 70 nm in order to find the balance between the number of approximations and resources used for the calculations. The influence of the vortex states and the effect of interwire interactions on the remagnetization curves is discussed. It has been found that 7 nanowires treated by a mean field model are able to reproduce well the reversal behaviour of the whole array in the case of large diameter nanowires. Vortex states tend to decrease the influence of the structural inhomogeneities on reversal process and thus lead to the increased predictability of the system.


title: Phonon Feedback Effects on Dynamics of Phase Slip Centers in Finite Gap Superconductors

abstract: The results on the behavior of phase-slip centers in thin superconducting wires based on finite-element modeling and time-dependent Ginzburg-Landau (TDGL) equations are discussed. For closer relationship with experiments, we used finite-gap formulation of the TDGL system. Both the dynamic equation for the Cooper-pair condensate wave-function and the expression for the electric current are more complex than in the gapless case. On this basis, the influence of nonequilibrium phonons is explored. These phonons can essentially change the location of geometrical points in which the phase slippage takes place. They also affect the frequency of phase-slip oscillations. The reported effects are experimentally detectable and can be used in practical devices.


title: Impact of COVID-19 on Air Quality in Israel

abstract: The COVID-19 pandemic has caused, in general, a sharp reduction in traffic and industrial activities. This in turn leaded to a reduction in air pollution around the world. It is important to quantity the amount of that reduction in order to estimate the influence weight of traffic and industrial activities over the total variation of air quality. The aim of this paper is to evaluate the impact of the COVID-19 outbreak on air pollution in Israel, which is considered one of the countries with a higher air pollution than other Western countries. The results reveal two main findings: 1. During the COVID-19 outbreak, relative to its earlier closest period, the pollution from transport, based on Nitrogen oxides, had reduced by 40$\%$ on average, whereas the pollution from industrial, based on Grand-level ozone, had increased by 34$\%$ on average. Relative to 2019, the COVID-19 outbreak caused a reduction in air pollution from transport and industrial as well. 2. The explanation percent of the time period of COVID-19 is at most 22$\%$ over the total variation of each pollutant amount.


title: Five-Dimensional Yang-Mills Black Holes in Massive Gravity's Rainbow

abstract: In this paper, we will analyze a five-dimensional Yang-Mills black hole solution in massive gravity's rainbow. We will also investigate the flow of such a solution with scale. Then, we will discuss the scale dependence of the thermodynamics for this black hole. Also, we study the criticality in the extended phase space by treating the cosmological constant as the thermodynamics pressure of this black hole solution. Moreover, we will use the partition function for this solution to obtaining corrections to the thermodynamics of this system and examine their key role on the behavior of corrected solutions.


title: PRI-VAE: Principle-of-Relevant-Information Variational Autoencoders

abstract: Although substantial efforts have been made to learn disentangled representations under the variational autoencoder (VAE) framework, the fundamental properties to the dynamics of learning of most VAE models still remain unknown and under-investigated. In this work, we first propose a novel learning objective, termed the principle-of-relevant-information variational autoencoder (PRI-VAE), to learn disentangled representations. We then present an information-theoretic perspective to analyze existing VAE models by inspecting the evolution of some critical information-theoretic quantities across training epochs. Our observations unveil some fundamental properties associated with VAEs. Empirical results also demonstrate the effectiveness of PRI-VAE on four benchmark data sets.


title: Towards Practical Lipreading with Distilled and Efficient Models

abstract: Lipreading has witnessed a lot of progress due to the resurgence of neural networks. Recent works have placed emphasis on aspects such as improving performance by finding the optimal architecture or improving generalization. However, there is still a significant gap between the current methodologies and the requirements for an effective deployment of lipreading in practical scenarios. In this work, we propose a series of innovations that significantly bridge that gap: first, we raise the state-of-the-art performance by a wide margin on LRW and LRW-1000 to 88.5% and 46.6%, respectively using self-distillation. Secondly, we propose a series of architectural changes, including a novel Depthwise Separable Temporal Convolutional Network (DS-TCN) head, that slashes the computational cost to a fraction of the (already quite efficient) original model. Thirdly, we show that knowledge distillation is a very effective tool for recovering performance of the lightweight models. This results in a range of models with different accuracy-efficiency trade-offs. However, our most promising lightweight models are on par with the current state-of-the-art while showing a reduction of 8.2x and 3.9x in terms of computational cost and number of parameters, respectively, which we hope will enable the deployment of lipreading models in practical applications.


title: Lyapunov exponents of the SHE for general initial data

abstract: We consider the $(1+1)$-dimensional stochastic heat equation (SHE) with multiplicative white noise and the Cole-Hopf solution of the Kardar-Parisi-Zhang (KPZ) equation. We show an exact way of computing the Lyapunov exponents of the SHE for a large class of initial data which includes any bounded deterministic positive initial data and the stationary initial data. As a consequence, we derive exact formulas for the upper tail large deviation rate functions of the KPZ equation for general initial data.


title: Mobility patterns of the Portuguese population during the COVID-19 pandemic

abstract: SARS-CoV-2 emerged in late 2019. Since then, it has spread to several countries, becoming classified as a pandemic. So far, there is no definitive treatment or vaccine, so the best solution is to prevent transmission between individuals through social distancing. However, it is difficult to measure the effectiveness of these distance measures. Therefore, this study uses data from Google COVID-19 Community Mobility Reports to try to understand the mobility patterns of the Portuguese population during the COVID-19 pandemic. In this study, the Rt value was modeled for Portugal. Also, the changepoint was calculated for the population mobility patterns. Thus, the change in the mobility pattern was used to understand the impact of social distance measures on the dissemination of COVID-19. As a result, it can be stated that the initial Rt value in Portugal was very close to 3, falling to values close to 1 after 25 days. Social isolation measures were adopted quickly. Furthermore, it was observed that public transport was avoided during the pandemic. Finally, until the emergence of a vaccine or an effective treatment, this is the new normal, and it must be understood that new patterns of mobility, social interaction, and hygiene must be adapted to this reality.


title: Industry Adoption Scenarios for Authoritative Data Stores using the ISDA Common Domain Model

abstract: In this paper we explore opportunities for the post-trade industry to standardize and simplify in order to significantly increase efficiency and reduce costs. We start by summarizing relevant industry problems (inconsistent processes, inconsistent data and duplicated data) and then present the corresponding potential industry solutions (process standardization, data standardization and authoritative data stores). This includes transitioning to the International Swaps and Derivatives Association Common Domain Model (CDM) as a standard set of digital representations for the business events and processes throughout the life cycle of a trade. We then explore how financial market infrastructures could operate authoritative data stores that make CDM business events available to broker-dealers, considering both traditional centralized models and potential decentralized models. For both types of model, there are many possible adoption scenarios (depending on each broker-dealer's degree of integration with the authoritative data store and usage of the CDM), and we identify some of the key scenarios.


title: A Lagrangian Perturbation Theory in the presence of massive neutrinos

abstract: We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using Convolution-LPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote simulation suite, we find good accuracy down to $r=20 \,\text{Mpc} \, h^{-1}$ at redshift $z=0.5$, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to $r=20 \,\text{Mpc} \, h^{-1}$ when comparing to simulated halos.


title: Terahertz spectroscopy evidence of possible 40 K superconductivity in rhenium-doped strontium ruthenates

abstract: Strontium ruthenates have many similarities with copper oxide superconductors and are of particular interest for the investigation of the mechanisms and conditions which lead to high-temperature superconductivity. We report here on multiple experimental indications of superconductivity with onset at 40 K in strontium ruthenate doped by rhenium and selenium with chlorine used as the flux. The main experimental evidence arises from terahertz spectroscopy of this material followed by AC and DC magnetization, as well as measurements of its heat capacity and magnetoresistance. Structural and morphological studies revealed the heterophase nature of this polycrystalline material as well as the changes of lattice parameters relative to the original phases. Experimental data show a higher critical temperature on the surface compared to that of the bulk of the sample.


title: Mean-variance-utility portfolio selection with time and state dependent risk aversion

abstract: Under mean-variance-utility framework, we propose a new portfolio selection model, which allows wealth and time both have influences on risk aversion in the process of investment. We solved the model under a game theoretic framework and analytically derived the equilibrium investment (consumption) policy. The results conform with the facts that optimal investment strategy heavily depends on the investor's wealth and future income-consumption balance as well as the continuous optimally consumption process is highly dependent on the consumption preference of the investor.


title: An Enhanced Text Classification to Explore Health based Indian Government Policy Tweets

abstract: Government-sponsored policy-making and scheme generations is one of the means of protecting and promoting the social, economic, and personal development of the citizens. The evaluation of effectiveness of these schemes done by government only provide the statistical information in terms of facts and figures which do not include the in-depth knowledge of public perceptions, experiences and views on the topic. In this research work, we propose an improved text classification framework that classifies the Twitter data of different health-based government schemes. The proposed framework leverages the language representation models (LR models) BERT, ELMO, and USE. However, these LR models have less real-time applicability due to the scarcity of the ample annotated data. To handle this, we propose a novel GloVe word embeddings and class-specific sentiments based text augmentation approach (named Mod-EDA) which boosts the performance of text classification task by increasing the size of labeled data. Furthermore, the trained model is leveraged to identify the level of engagement of citizens towards these policies in different communities such as middle-income and low-income groups.


title: Deep Learning for Distributed Channel Feedback and Multiuser Precoding in FDD Massive MIMO

abstract: This paper shows that deep neural network (DNN) can be used for efficient and distributed channel estimation, quantization, feedback, and downlink multiuser precoding for a frequency-division duplex massive multiple-input multiple-output system in which a base station (BS) serves multiple mobile users, but with rate-limited feedback from the users to the BS. A key observation is that the multiuser channel estimation and feedback problem can be thought of as a distributed source coding problem. In contrast to the traditional approach where the channel state information (CSI) is estimated and quantized at each user independently, this paper shows that a joint design of pilots and a new DNN architecture, which maps the received pilots directly into feedback bits at the user side then maps the feedback bits from all the users directly into the precoding matrix at the BS, can significantly improve the overall performance. This paper further proposes robust design strategies with respect to channel parameters and also a generalizable DNN architecture for varying number of users and number of feedback bits. Numerical results show that the DNN-based approach with short pilot sequences and very limited feedback overhead can already approach the performance of conventional linear precoding schemes with full CSI.


title: Kolmogorov flow: Linear Stability and Energy Transfers in a minimal low-dimensional model

abstract: In this paper, we derive a four-mode model for the Kolmogorov flow by employing Galerkin truncation and Craya-Herring basis for the decomposition of velocity field. After this, we perform a bifurcation analysis of the model. Though our low-dimensional model has fewer modes than the past models, it captures the essential features of the primary bifurcation of the Kolmogorov flow. For example, it reproduces the critical Reynolds number for the supercritical pitchfork bifurcation and the flow structures of the past works. We also demonstrate energy transfers from intermediate scales to large scales. We perform direct numerical simulations of the Kolmogorov flow and show that our model predictions match with the numerical simulations very well.


title: Discovery of a hot ultramassive rapidly rotating DBA White Dwarf

abstract: We report the discovery of a nearby massive white dwarf with He-H atmosphere. The white dwarf is located at a distance of $74.5\pm0.9$ pc. Its radius, mass, effective temperature, H/He ratio and age are $R=2500\pm100$ km, $M=1.33\pm0.01$ $\rm M_{\odot}$, $T_{\mathrm{eff}}=31200\pm 1200$ K, ${\mathrm{H/He}}\sim0.1$ and $330\pm40$ Myr, respectively. The observed spectrum is redshifted by $V_{\mathrm{r}}=+240\pm15$ km s$^{-1},$ which is mostly attributed to the gravitational redshift. The white dwarf shows a regular stable photometric variability with amplitude $\Delta g\approx0.06^{\mathrm{m}}$ and period $P=353.456$ s suggesting rapid rotation. This massive, hot and rapidly rotating white dwarf is likely to originate from the merging of close binary white dwarf system that avoided explosion in a thermonuclear type Ia supernova at the Carboniferous Period of the Earth history.


title: Theory of the thickness dependence of the charge density wave transition in 1T-TiTe$_2$

abstract: Most metallic transition metal dichalcogenides undergo charge density wave (CDW) instabilities with similar or identical ordering vectors in bulk and in single layer, albeit with different critical temperatures. Metallic 1T-TiTe$_2$ is a remarkable exception as it shows no evidence of charge density wave formation in bulk, but it displays a stable $2\times2$ reconstruction in single-layer form. The mechanism for this 3D-2D crossover of the transition is still unclear, although strain from the substrate and the exchange interaction have been pointed out as possible formation mechanisms. Here, by performing non-perturbative anharmonic calculations with gradient corrected and hybrid functionals, we explain the thickness behaviour of the transition in 1T-TiTe$_2$. We demonstrate that the occurrence of the CDW in single-layer TiTe$_2$ occurs from the interplay of non-perturbative anharmonicity and an exchange enhancement of the electron-phonon interaction, larger in the single layer than in the bulk. Finally, we study the electronic and structural properties of the single-layer CDW phase and provide a complete description of its electronic structure, phonon dispersion as well as infrared and Raman active phonon modes.


title: Uncertain-DeepSSM: From Images to Probabilistic Shape Models

abstract: Statistical shape modeling (SSM) has recently taken advantage of advances in deep learning to alleviate the need for a time-consuming and expert-driven workflow of anatomy segmentation, shape registration, and the optimization of population-level shape representations. DeepSSM is an end-to-end deep learning approach that extracts statistical shape representation directly from unsegmented images with little manual overhead. It performs comparably with state-of-the-art shape modeling methods for estimating morphologies that are viable for subsequent downstream tasks. Nonetheless, DeepSSM produces an overconfident estimate of shape that cannot be blindly assumed to be accurate. Hence, conveying what DeepSSM does not know, via quantifying granular estimates of uncertainty, is critical for its direct clinical application as an on-demand diagnostic tool to determine how trustworthy the model output is. Here, we propose Uncertain-DeepSSM as a unified model that quantifies both, data-dependent aleatoric uncertainty by adapting the network to predict intrinsic input variance, and model-dependent epistemic uncertainty via a Monte Carlo dropout sampling to approximate a variational distribution over the network parameters. Experiments show an accuracy improvement over DeepSSM while maintaining the same benefits of being end-to-end with little pre-processing.


title: TITUS: Visualization of Neutrino Events in Liquid Argon Time Projection Chambers

abstract: The amount and complexity of data recorded by high energy physics experiments are rapidly growing, and with these grow the difficulties in visualizing such data. To study the physics of neutrinos, a type of elementary particle, scientists use liquid argon time projection chamber (LArTPC) detectors, among other technologies. LArTPCs have a very high spatial resolution and resolve many of the elementary particles that come out of a neutrino interacting within the argon in the detector. Visualizing these neutrino interactions is of fundamental importance to understanding the properties of neutrinos, but also monitoring and checking on the detector conditions and operations. From these ideas, we have developed TITUS, an event display that shows images recorded by these neutrino detectors. TITUS is a piece of software that reads data coming from LArTPC detectors (as well as the corresponding simulation) and allows users to explore such data in multiple ways. TITUS is flexible to enable fast prototyping and customization.


title: NNLO QCD$\times$EW corrections to on-shell $Z$ production

abstract: We present the first analytical results for the ${\cal O}(\alpha\alpha_s)$ corrections to the total cross section for the inclusive production of an on-shell $Z$ boson at hadron colliders. We include the complete set of contributions, with photonic and massive weak gauge boson effects, which have been computed in analytical form and expressed in terms of polylogarithmic and elliptic integrals. We present numerical results, relevant for the precision studies at the LHC. These corrections increase the accuracy of the predictions and contribute to the reduction of the QCD component of the theoretical uncertainty.


title: An integral version of Zariski decompositions on normal surfaces

abstract: We show that any pseudo-effective divisor on a normal surface decomposes uniquely into its "integral positive" part and "integral negative" part, which is an integral analog of Zariski decompositions. By using this decomposition, we give three applications: a vanishing theorem of divisors on surfaces (a generalization of Kawamata-Viehweg and Miyaoka vanishing theorems), Reider-type theorems of adjoint linear systems on surfaces (including a log version and a relative version of the original one) and extension theorems of morphisms defined on curves on surfaces (generalizations of Serrano and Paoletti's results).


title: Representation Formula for Viscosity Solution to a PDE Problem involving Pucci's Extremal Operator

abstract: We provide a representation formula for viscosity solutions to an elliptic Dirichlet problem involving Pucci's extremal operators. This is done through a dynamic programming principle derived from Denis, Hu and Peng (2010). The formula can be seen as a nonlinear extension of the Feynman-Kac formula.


title: Predicting the long-term stability of compact multiplanet systems

abstract: We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.


title: Undular diffusion in nonlinear sigma models

abstract: We discuss general features of charge transport in non-relativistic classical field theories invariant under non-abelian unitary Lie groups by examining the full structure of two-point dynamical correlation functions in grand-canonical ensembles at finite charge densities (polarized ensembles). Upon explicit breaking of non-abelian symmetry, two distinct transport laws characterized by dynamical exponent $z=2$ arise. While in the unbroken symmetry sector the Cartan fields exhibit normal diffusion, the transversal sectors governed by the nonlinear analogues of Goldstone modes disclose an unconventional law of diffusion characterized by a complex diffusion constant and undulating patterns in the spatiotemporal correlation profiles. In the limit of strong polarization, one retrieves the imaginary-time diffusion for uncoupled linear Goldstone modes, whereas for weak polarizations the imaginary component of the diffusion constant becomes small. In models of higher rank symmetry, we prove absence of dynamical correlations among distinct transversal sectors.


title: The Calder\'on problem in the $L^p$ framework on Riemann surfaces

abstract: The purpose of this article is to extend the uniqueness results for the two dimensional Calder\'on problem to unbounded potentials on general geometric settings. We prove that the Cauchy data sets for Schr\"odinger equations uniquely determines potentials in $L^{p}$ for $p> 4/3$. In doing so, we first recover singularities of the potential, from which point a $L^2$-based method of stationary phase can be applied. Both steps are done via constructions of complex geometric optic solutions and Carleman estimates.


title: Tensor-based techniques for fast discretization and solution of 3D elliptic equations with random coefficients

abstract: In this paper, we propose and analyze the numerical algorithms for fast solution of periodic elliptic problems in random media in $\mathbb{R}^d$, $d=2,3$. We consider the stochastic realizations using checkerboard configuration of the equation coefficients built on a large $L \times L \times L$ lattice, where $L$ is the size of representative volume elements. The Kronecker tensor product scheme is introduced for fast generation of the stiffness matrix for FDM discretization on a tensor grid. We describe tensor techniques for the construction of the low Kronecker rank spectrally equivalent preconditioner in periodic setting to be used in the framework of PCG iteration. In our construction the diagonal matrix of the discrete Laplacian inverse represented in the Fourier basis is reshaped into a 3D tensor, which is then approximated by a low-rank canonical tensor, calculated by the multigrid Tucker-to-canonical tensor transform. The FDM discretization scheme on a tensor grid is described in detail, and the computational characteristics in terms of $L$ for the 3D Matlab implementation of the PCG iteration are illustrated. The present work continues the developments in [22], where the numerical primer to study the asymptotic convergence rate vs. $L$ for the homogenized matrix for 2D elliptic PDEs with random coefficients was investigated numerically. The presented elliptic problem solver can be applied for calculation of long sequences of stochastic realizations in numerical analysis of 3D stochastic homogenization problems, for solving 3D quasi-periodic geometric homogenization problems, as well as in the numerical simulation of dynamical many body interaction processes and multi-particle electrostatics.


title: Split Higgs triplet

abstract: We find if a Higgs triplet with hypercharge has a special dimension-6 operator with the standard model Higgs doublet, i.e. a square of the trilinear triplet-doublet coupling, its scalar(pseudo-scalar) component can obtain a small quadratic term while its pseudo-scalar(scalar) and charged-scalar components can hold their masses heavy enough. Such split Higgs triplet can spontaneously develop a small vacuum expectation value to realize a Majorana neutrino mass generation without causing any high-dimensional lepton number violations including the well-known Weinberg dimension-5 operator. Alternatively it can mediate a non-standard neutrino self-interaction motivated by resolving the tension in Hubble constant measurements. This effective theory with rich observable phenomena can be induced by the Georgi-Machacek Higgs triplets at tree level or some dark matter fields at one-loop order.


title: Radius and equation of state constraints from massive neutron stars and GW190814

abstract: Motivated by the unknown nature of the $2.50-2.67\,M_\odot$ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. We demonstrate that a modification of the highly uncertain supra-saturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of $2.6\,M_\odot$ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of $\sim 2.0\,M_\odot$ neutron stars (but not necessarily the radii of $1.4\,M_\odot$ neutron stars), which further motivates future NICER radius measurements of PSR J1614-2230 and PSR J0740+6620.


title: Unique Continuation Property for Biharmonic Hypersurfaces in Spheres

abstract: We study properties of non-minimal biharmonic hypersurfaces of spheres. The main result is a CMC Unique Continuation Theorem for biharmonic hypersurfaces of spheres. We then deduce new rigidity theorems to support the Conjecture that biharmonic submanifolds of Euclidean spheres must be of constant mean curvature.


title: Random extrapolation for primal-dual coordinate descent

abstract: We introduce a randomly extrapolated primal-dual coordinate descent method that adapts to sparsity of the data matrix and the favorable structures of the objective function. Our method updates only a subset of primal and dual variables with sparse data, and it uses large step sizes with dense data, retaining the benefits of the specific methods designed for each case. In addition to adapting to sparsity, our method attains fast convergence guarantees in favorable cases \textit{without any modifications}. In particular, we prove linear convergence under metric subregularity, which applies to strongly convex-strongly concave problems and piecewise linear quadratic functions. We show almost sure convergence of the sequence and optimal sublinear convergence rates for the primal-dual gap and objective values, in the general convex-concave case. Numerical evidence demonstrates the state-of-the-art empirical performance of our method in sparse and dense settings, matching and improving the existing methods.


title: Interpreting deep learning models for weak lensing

abstract: Deep Neural Networks (DNNs) are powerful algorithms that have been proven capable of extracting non-Gaussian information from weak lensing (WL) data sets. Understanding which features in the data determine the output of these nested, non-linear algorithms is an important but challenging task. We analyze a DNN that has been found in previous work to accurately recover cosmological parameters in simulated maps of the WL convergence ($\kappa$). We derive constraints on the cosmological parameter pair $(\Omega_m,\sigma_8)$ from a combination of three commonly used WL statistics (power spectrum, lensing peaks, and Minkowski functionals), using ray-traced simulated $\kappa$ maps. We show that the network can improve the inferred parameter constraints relative to this combination by $20\%$ even in the presence of realistic levels of shape noise. We apply a series of well established saliency methods to interpret the DNN and find that the most relevant pixels are those with extreme $\kappa$ values. For noiseless maps, regions with negative $\kappa$ account for $86-69\%$ of the attribution of the DNN output, defined as the square of the saliency in input space. In the presence of shape nose, the attribution concentrates in high convergence regions, with $36-68\%$ of the attribution in regions with $\kappa > 3 \sigma_{\kappa}$.


title: Race and gender income inequality in the USA: black women vs. white men

abstract: Income inequality between different races in the U.S. is especially large. This difference is even larger when gender is involved. In a complementary study, we have developed a dynamic microeconomic model accurately describing the evolution of male and female incomes since 1930. Here, we extend our analysis and model the disparity between black and white population in the U.S., separately for males and females. Unfortunately, income microdata provided by the U.S. Census Bureau for other races and ethnic groups are not time compatible or too short for modelling purposes. We are forced to constrain our analysis to black and white population, but all principal results can be extrapolated to other races and ethnicities. Our analysis shows that black females and white males are two poles of the overall income inequality. The prediction of income distribution for two extreme cases with one model is the main challenge of this study.


title: A Robotic Framework for Making Eye Contact with Humans

abstract: Meeting eye contact is the essential prerequisite skill of a human to initiate any conversation with others. However, it is not an easy task for a robot to meet eye contact with a human if they are not facing each other initially or the human is intensely engaged his or her task. If the robot would like to start communication with a particular person, it should turn its gaze to that person first. However, only such a turning action alone is not always enough to set up eye contact. Sometimes, the robot should perform some strong actions so that it can capture the human's attention toward it. In this paper, we proposed a computational model for robots that can proactively capture human attention and makes eye contact with him or her. Evaluation experiments by using a robotic head reveal the effectiveness of the proposed model in different viewing situations.


title: Graphene Domain Signature of Raman Spectra of sp2 Amorphous Carbons

abstract: The paper presents a joint consideration of Raman spectra of sp2 amorphous carbons alongside with the nature and type of their amorphicity. The latter was attributed to the enforced fragmentation. The fragments, presented with size-restricted graphene domains with heteroatom necklaces in the circumference, are the basic structural units (BSUs) of the solids, determining them as amorphics with molecular structure. The standard G-D-2D pattern of Raman spectra of polycyclic aromatic hydrocarbons, sp2 amorphous carbons, graphene and/or graphite crystal is attributed to BSUs graphene domains. The molecular approximation allows connecting the G-D spectra image of one-phonon spectra with a considerable dispersion of the C=C bond lengths within graphene domains, governed by size, heteroatom necklace of BSUs as well as BSUs packing. The interpretation of 2D two-phonon spectra reveals a particular role of electrical anharmonicity in the spectra formation and attributes this effect to a high degree of the electron density delocalization in graphene domains.


title: S2RMs: Spatially Structured Recurrent Modules

abstract: Capturing the structure of a data-generating process by means of appropriate inductive biases can help in learning models that generalize well and are robust to changes in the input distribution. While methods that harness spatial and temporal structures find broad application, recent work has demonstrated the potential of models that leverage sparse and modular structure using an ensemble of sparingly interacting modules. In this work, we take a step towards dynamic models that are capable of simultaneously exploiting both modular and spatiotemporal structures. We accomplish this by abstracting the modeled dynamical system as a collection of autonomous but sparsely interacting sub-systems. The sub-systems interact according to a topology that is learned, but also informed by the spatial structure of the underlying real-world system. This results in a class of models that are well suited for modeling the dynamics of systems that only offer local views into their state, along with corresponding spatial locations of those views. On the tasks of video prediction from cropped frames and multi-agent world modeling from partial observations in the challenging Starcraft2 domain, we find our models to be more robust to the number of available views and better capable of generalization to novel tasks without additional training, even when compared against strong baselines that perform equally well or better on the training distribution.


title: Energy-Momentum portal to dark matter and emergent gravity

abstract: We propose a new scenario where dark matter belongs to a secluded sector coupled to the Standard Model through energy--momentum tensors. Our model is motivated by constructions where gravity {\it emerges} from a hidden sector, the graviton being identified by the kinetic term of the fields in the secluded sector. Supposing that the lighter particle of the secluded sector is the dark component of the Universe, we show that we can produce it in a sufficiently large amount despite the suppressed couplings of the theory, thanks to large temperatures of the thermal bath in the early stage of the Universe.


title: The New Digital Platforms: Merger Control in Pakistan

abstract: The Pakistan competition policy, as in many other countries, was originally designed to regulate business conduct in traditional markets and for tangible goods and services. However, the development and proliferation of the internet has led to the emergence of digital companies which have disrupted many sectors of the economy. These platforms provide digital infrastructure for a range of services including search engines, marketplaces, and social networking sites. The digital economy poses a myriad of challenges for competition authorities worldwide, especially with regard to digital mergers and acquisitions (M&As). While some jurisdictions such as the European Union and the United States have taken significant strides in regulating technological M&As, there is an increasing need for developing countries such as Pakistan to rethink their competition policy tools. This paper investigates whether merger reviews in the Pakistan digital market are informed by the same explanatory variables as in the traditional market, by performing an empirical comparative analysis of the Competition Commission of Pakistan's (CCP's) M&A decisions between 2014 and 2019. The findings indicate the CCP applies the same decision factors in reviewing both traditional and digital M&As. As such, this paper establishes a basis for igniting the policy and economic debate of regulating the digital platform industry in Pakistan.


title: Co-t-structures, cotilting and cotorsion pairs

abstract: Let $\mathsf{T}$ be a triangulated category with shift functor $\Sigma \colon \mathsf{T} \to \mathsf{T}$. Suppose $(\mathsf{A},\mathsf{B})$ is a co-t-structure with coheart $\mathsf{S} = \Sigma \mathsf{A} \cap \mathsf{B}$ and extended coheart $\mathsf{C} = \Sigma^2 \mathsf{A} \cap \mathsf{B} = \mathsf{S} * \Sigma \mathsf{S}$, which is an extriangulated category. We show that there is a bijection between co-t-structures $(\mathsf{A}',\mathsf{B}')$ in $\mathsf{T}$ such that $\mathsf{A} \subseteq \mathsf{A}' \subseteq \Sigma \mathsf{A}$ and complete cotorsion pairs in the extended coheart $\mathsf{C}$. In the case that $\mathsf{T}$ is Hom-finite, $\mathbf{k}$-linear and Krull-Schmidt, we show further that there is a bijection between complete cotorsion pairs in $\mathsf{C}$ and functorially finite torsion pairs in $\mathsf{mod}\, \mathsf{S}$.


title: Blockchain-Federated-Learning and Deep Learning Models for COVID-19 detection using CT Imaging

abstract: With the increase of COVID-19 cases worldwide, an effective way is required to diagnose COVID-19 patients. The primary problem in diagnosing COVID-19 patients is the shortage and reliability of testing kits, due to the quick spread of the virus, medical practitioners are facing difficulty identifying the positive cases. The second real-world problem is to share the data among the hospitals globally while keeping in view the privacy concerns of the organizations. Building a collaborative model and preserving privacy are major concerns for training a global deep learning model. This paper proposes a framework that collects a small amount of data from different sources (various hospitals) and trains a global deep learning model using blockchain based federated learning. Blockchain technology authenticates the data and federated learning trains the model globally while preserving the privacy of the organization. First, we propose a data normalization technique that deals with the heterogeneity of data as the data is gathered from different hospitals having different kinds of CT scanners. Secondly, we use Capsule Network-based segmentation and classification to detect COVID-19 patients. Thirdly, we design a method that can collaboratively train a global model using blockchain technology with federated learning while preserving privacy. Additionally, we collected real-life COVID-19 patients data, which is, open to the research community. The proposed framework can utilize up-to-date data which improves the recognition of computed tomography (CT) images. Finally, our results demonstrate a better performance to detect COVID-19 patients.


title: Nichols algebras over classical Weyl groups

abstract: We show that except in several cases conjugacy classes of classical Weyl groups $W(B_n)$ and $W(D_n)$ are of type {\rm D}. We prove that except in three cases Nichols algebras of irreducible Yetter-Drinfeld ({\rm YD} in short )modules over the classical Weyl groups are infinite dimensional.


title: Benchmarking 16-element quantum search algorithms on superconducting quantum processors

abstract: We present experimental results on running 4-qubit unstructured search on IBM quantum processors. Our best attempt attained probability of success around 24.5%. We try several algorithms and use the most recent developments in quantum search to reduce the number of entangling gates that are currently considered the main source of errors in quantum computations. Comparing theoretical expectations of an algorithm performance with the actual data, we explore the hardware limits, showing sharp, phase-transition-like degradation of performance on quantum processors. We conclude that it is extremely important to design hardware-aware algorithms and to include any other low level optimizations on NISQ devices.


title: Open Data Quality

abstract: The research discusses how (open) data quality could be described, what should be considered developing a data quality management solution and how it could be applied to open data to check its quality. The proposed approach focuses on development of data quality specification which can be executed to get data quality evaluation results, find errors in data and possible problems which must be solved. The proposed approach is applied to several open data sets to evaluate their quality. Open data is very popular, free available for every stakeholder - it is often used to make business decisions. It is important to be sure that this data is trustable and error-free as its quality problems can lead to huge losses.


title: Bayesian Modeling of COVID-19 Positivity Rate -- the Indiana experience

abstract: In this short technical report we model, within the Bayesian framework, the rate of positive tests reported by the the State of Indiana, accounting also for the substantial variability (and overdispeartion) in the daily count of the tests performed. The approach we take, results with a simple procedure for prediction, a posteriori, of this rate of 'positivity' and allows for an easy and a straightforward adaptation by any agency tracking daily results of COVID-19 tests. The numerical results provided herein were obtained via an updatable R Markdown document.


title: Loss Function Search for Face Recognition

abstract: In face recognition, designing margin-based (e.g., angular, additive, additive angular margins) softmax loss functions plays an important role in learning discriminative features. However, these hand-crafted heuristic methods are sub-optimal because they require much effort to explore the large design space. Recently, an AutoML for loss function search method AM-LFS has been derived, which leverages reinforcement learning to search loss functions during the training process. But its search space is complex and unstable that hindering its superiority. In this paper, we first analyze that the key to enhance the feature discrimination is actually \textbf{how to reduce the softmax probability}. We then design a unified formulation for the current margin-based softmax losses. Accordingly, we define a novel search space and develop a reward-guided search method to automatically obtain the best candidate. Experimental results on a variety of face recognition benchmarks have demonstrated the effectiveness of our method over the state-of-the-art alternatives.


title: Dynamics of ternary statistical experiments with equilibrium state

abstract: We study the scenarios of the dynamics of ternary statistical experiments, modeled employing difference equations. The important features are a balance condition and the existence of a steady-state (equilibrium). We give a classification of scenarios of the model evolution which are significantly different between them, depending on the domain of the values of the model basic parameters.


title: Free-running SIMilarity-Based Angiography (SIMBA) for simplified anatomical MR imaging of the heart

abstract: Purpose: Whole-heart MRA techniques typically target pre-determined motion states and address cardiac and respiratory dynamics independently. We propose a novel fast reconstruction algorithm, applicable to ungated free-running sequences, that leverages inherent similarities in the acquired data to avoid such physiological constraints. Theory and Methods: The proposed SIMilarity-Based Angiography (SIMBA) method clusters the continuously acquired k-space data in order to find a motion-consistent subset that can be reconstructed into a motion-suppressed whole-heart MRA. Free-running 3D radial datasets from six ferumoxytol-enhanced scans of pediatric cardiac patients and twelve non-contrast scans of healthy volunteers were reconstructed with a non-motion-suppressed regridding of all the acquired data (All Data), our proposed SIMBA method, and a previously published free-running framework (FRF) that uses cardiac and respiratory self-gating and compressed sensing. Images were compared for blood-myocardium interface sharpness, contrast ratio, and visibility of coronary artery ostia. Results: Both the fast SIMBA reconstruction (~20s) and the FRF provided significantly higher blood-myocardium sharpness than All Data (P<0.001). No significant difference was observed among the former two. Significantly higher blood-myocardium contrast ratio was obtained with SIMBA compared to All Data and FRF (P<0.01). More coronary ostia could be visualized with both SIMBA and FRF than with All Data (All Data: 4/36, SIMBA: 30/36, FRF: 33/36, both P<0.001) but no significant difference was found between the first two. Conclusion: The combination of free-running sequences and the fast SIMBA reconstruction, which operates without a priori assumptions related to physiological motion, forms a simple workflow for obtaining whole-heart MRA with sharp anatomical structures.


title: Stability of Parker's Solar Wind Solution Near the Solar Surface

abstract: Stability of Parker's steady solar wind solution near the solar surface is systematically investigated by posing a Sturm-Liouville problem for this solution. Parker's solar wind solution, whether it turns into a breeze or a supersonic wind depending on the pressure in the interstellar medium, is shown to describe the solar wind to start in a stable way as it leaves the solar surface from a state of: rest, co-rotation with the Sun, slow motion. The isothermal gas flow assumption in Parker's solar wind model is then relaxed, and more realistic barotropic fluid and diabatic models are used for the gas flow. The stability of the solar wind flow, as it starts from a state of rest at the solar surface, is shown to continue to hold. Parker's solar wind solution therefore appears to be a stable attractor of this dynamical system.


title: On the bordism group for group actions on the torus

abstract: In this short note, we study the bordism problem for group actions on the torus and give examples of groups acting on the torus by diffeomorphisms isotopic to the identity that cannot be extended to an action on a bounding 3-manifold. This solves a question raised in the previous work of the authors.


title: Constraining annihilating dark matter mass by the radio continuum spectral data of a high-redshift galaxy cluster

abstract: In the past decade, the properties of annihilating dark matter models were examined by various kinds of data, including the data of gamma rays, radio waves, X-ray, positrons, electrons, antiprotons and neutrinos. In particular, most of the studies focus on the data of our Galaxy, nearby galaxies (e.g. M31 galaxy) or nearby galaxy clusters (e.g. Fornax cluster). In this article, we examine the archival radio continuum spectral data of a relatively high-redshift galaxy cluster (A697 cluster) to constrain the properties of annihilating dark matter. We find that leptophilic annihilation channels ($e^+e^-$, $\mu^+\mu^-$ and $\tau^+\tau^-$) can give very good fits to the radio continuum spectrum of the A697 cluster.


title: Relations between scaling exponents in unimodular random graphs

abstract: We investigate the validity of the "Einstein relations" in the general setting of unimodular random networks. These are equalities relating scaling exponents: $d_w = d_f + \tilde{\zeta}$ and $d_s = 2 d_f/d_w$, where $d_w$ is the walk dimension, $d_f$ is the fractal dimension, $d_s$ is the spectral dimension, and $\tilde{\zeta}$ is the resistance exponent. Roughly speaking, this relates the mean displacement and return probability of a random walker to the density and conductivity of the underlying medium. We show that if $d_f$ and $\tilde{\zeta} \geq 0$ exist, then $d_w$ and $d_s$ exist, and the aforementioned equalities hold. Moreover, our primary new estimate is the relation $d_w \geq d_f + \tilde{\zeta}$, which is established for all $\tilde{\zeta} \in \mathbb{R}$. For the uniform infinite planar triangulation (UIPT), this yields the consequence $d_w=4$ using $d_f=4$ (Angel 2003) and $\tilde{\zeta}=0$ (established here as a consequence of the Liouville Quantum Gravity theory, following Gwynne-Miller 2017 and Ding-Gwynne 2020). The conclusion $d_w=4$ had been previously established by Gwynne and Hutchcroft (2018) using more elaborate methods. A new consequence is that $d_w = d_f$ for the uniform infinite Schnyder-wood decorated triangulation, implying that the simple random walk is subdiffusive, since $d_f > 2$ (Ding and Gwynne 2020). For the random walk on $\mathbb{Z}^2$ driven by conductances from an exponentiated Gaussian free field with exponent $\gamma > 0$, one has $d_f = d_f(\gamma)$ and $\tilde{\zeta}=0$ (Biskup, Ding, and Goswami 2020). This yields $d_s=2$ and $d_w = d_f$, confirming two predictions of those authors.


title: Preferences Single-Peaked on a Tree: Multiwinner Elections and Structural Results

abstract: A preference profile is single-peaked on a tree if the candidate set can be equipped with a tree structure so that the preferences of each voter are decreasing from their top candidate along all paths in the tree. This notion was introduced by Demange (1982), and subsequently Trick (1989) described an efficient algorithm for deciding if a given profile is single-peaked on a tree. We study the complexity of multiwinner elections under several variants of the Chamberlin-Courant rule for preferences single-peaked on trees. We show that the egalitarian version of this problem admits a polynomial-time algorithm. For the utilitarian version, we prove that winner determination remains NP-hard, even for the Borda scoring function; however, a winning committee can be found in polynomial time if either the number of leaves or the number of internal vertices of the underlying tree is bounded by a constant. To benefit from these positive results, we need a procedure that can determine whether a given profile is single-peaked on a tree that has additional desirable properties (such as, e.g., a small number of leaves). To address this challenge, we develop a structural approach that enables us to compactly represent all trees with respect to which a given profile is single-peaked. We show how to use this representation to efficiently find the best tree for a given profile for use with our winner determination algorithms: Given a profile, we can efficiently find a tree with the minimum number of leaves, or a tree with the minimum number of internal vertices among trees on which the profile is single-peaked. We also consider several other optimization criteria for trees: for some we obtain polynomial-time algorithms, while for others we show NP-hardness results.


title: Reconstruction in one dimension from unlabeled Euclidean lengths

abstract: Let $G$ be a $3$-connected ordered graph with $n$ vertices and $m$ edges. Let $\mathbf{p}$ be a randomly chosen mapping of these $n$ vertices to the integer range $\{1, 2,3, \ldots, 2^b\}$ for $b\ge m^2$. Let $\ell$ be the vector of $m$ Euclidean lengths of $G$'s edges under $\mathbf{p}$. In this paper, we show that, with high probability over $\mathbf{p}$, we can efficiently reconstruct both $G$ and $\mathbf{p}$ from $\ell$. This reconstruction problem is NP-HARD in the worst case, even if both $G$ and $\ell$ are given. We also show that our results stand in the presence of small amounts of error in $\ell$, and in the real setting, with sufficiently accurate length measurements. Our method combines lattice reduction, which has previously been used to solve random subset sum problems, with an algorithm of Seymour that can efficiently reconstruct an ordered graph given an independence oracle for its matroid.


title: Critical Islands

abstract: We discuss a doubly-holographic prescription for black holes in braneworlds with a vanishing cosmological constant. It involves calculating Ryu-Takayanagi surfaces in AdS black funnel spacetimes attached to braneworld black holes in the $ critical$ Randall-Sundrum II model. Critical braneworlds have the virtue of having massless gravitons. Our approach should be useful when the braneworld is a cosmological black hole interacting with deconfined, large-$N$ matter. In higher dimensions, explicit funnel metrics will have to be constructed numerically -- but based on the general structure of the geometry, we present a natural guess for where one might find the semi-classical island. In a 3-dimensional example where a toy analytic black funnel is known, we can check our guess by direct calculation. We argue that this resolves a version of the information paradox in these braneworld systems, by finding strong evidence for "cosmological islands". Comoving Ryu-Takayanagi surfaces and associated UV cut-offs on the brane, play natural roles.


title: Relaxing the I.I.D. Assumption: Adaptively Minimax Optimal Regret via Root-Entropic Regularization

abstract: We consider prediction with expert advice when data are generated from distributions varying arbitrarily within an unknown constraint set. This semi-adversarial setting includes (at the extremes) the classical i.i.d. setting, when the unknown constraint set is restricted to be a singleton, and the unconstrained adversarial setting, when the constraint set is the set of all distributions. The Hedge algorithm -- long known to be minimax (rate) optimal in the adversarial regime -- was recently shown to be simultaneously minimax optimal for i.i.d. data. In this work, we propose to relax the i.i.d. assumption by seeking adaptivity at all levels of a natural ordering on constraint sets. We provide matching upper and lower bounds on the minimax regret at all levels, show that Hedge with deterministic learning rates is suboptimal outside of the extremes, and prove that one can adaptively obtain minimax regret at all levels. We achieve this optimal adaptivity using the follow-the-regularized-leader (FTRL) framework, with a novel adaptive regularization scheme that implicitly scales as the square root of the entropy of the current predictive distribution, rather than the entropy of the initial predictive distribution. Finally, we provide novel technical tools to study the statistical performance of FTRL along the semi-adversarial spectrum.


title: Radial oscillations and gravitational wave echoes of strange stars for various equations of state

abstract: We study the radial oscillations of non-rotating strange stars and their characteristic echo frequencies for three Equations of state (EoSs), viz., MIT Bag model EoS, linear EoS and polytropic EoS. The frequencies of radial oscillations of these compact stars are computed for these EoSs. 22 lowest radial frequencies for each of these three EoSs have been computed. First, for each EoS, we have integrated Tolman-Oppenheimer-Volkoff (TOV) equations numerically to calculate the radial and pressure perturbations of strange stars. Next, the mass-radius relationships for these stars are obtained using these three EoSs. Then the radial frequencies of oscillations for these EoSs are calculated. Further, the characteristic gravitational wave echo frequencies and the repetition of echo frequencies of strange stars are computed for these EoSs. Our numerical results show that the radial frequencies and also echo frequencies vastly depend on the model and on the value of the model parameter. Our results also show that, the radial frequencies of strange stars are maximum for polytropic EoS in comparison to MIT Bag model EoS and linear EoS. Moreover, strange stars with MIT Bag model EoS and linear EoS are found to emit gravitational wave echoes. Whereas, strange stars with polytropic EoS are not emitting gravitational wave echoes.


title: Two-Dimensional Quantum Walk of Correlated Photons

abstract: Quantum walks in an elaborately designed graph, is a powerful tool simulating physical and topological phenomena, constructing analog quantum algorithms and realizing universal quantum computing. Integrated photonics technology has emerged as a versatile platform to implement various quantum information tasks and a promising candidate to perform large-scale quantum walks. Both extending physical dimensions and involving more particles will increase the complexity of the evolving systems and the desired quantum resources. Pioneer works have demonstrated single particle walking on two-dimensional (2D) lattices and multiple walkers interfering on a one-dimensional structure. However, 2D multi-particle quantum walk, genuinely being not classically simulatable, has been a vacancy for nearly ten years. Here, we present a genuine 2D quantum walk with correlated photons on a triangular photonic lattice, which can be mapped to a state space up to 37X37 dimensions. This breaks through the physically restriction of single-particle evolution, which can encode information in a large space and constitute high-dimensional graphs indeed beneficial to quantum information processing. A site-by-site addressing between the chip facet and the 2D fanout interface enables an observation of over 600 non-classical interferences simultaneously, violating a classical limit up to 57 standard deviations. Our platform offers a promising prospect for multi-photon quantum walks in a large-scale 2D arrangement, paving the way for practical quantum simulation and quantum computation beyond classical regime.


title: Adversarial robustness via robust low rank representations

abstract: Adversarial robustness measures the susceptibility of a classifier to imperceptible perturbations made to the inputs at test time. In this work we highlight the benefits of natural low rank representations that often exist for real data such as images, for training neural networks with certified robustness guarantees. Our first contribution is for certified robustness to perturbations measured in $\ell_2$ norm. We exploit low rank data representations to provide improved guarantees over state-of-the-art randomized smoothing-based approaches on standard benchmark datasets such as CIFAR-10 and CIFAR-100. Our second contribution is for the more challenging setting of certified robustness to perturbations measured in $\ell_\infty$ norm. We demonstrate empirically that natural low rank representations have inherent robustness properties, that can be leveraged to provide significantly better guarantees for certified robustness to $\ell_\infty$ perturbations in those representations. Our certificate of $\ell_\infty$ robustness relies on a natural quantity involving the $\infty \to 2$ matrix operator norm associated with the representation, to translate robustness guarantees from $\ell_2$ to $\ell_\infty$ perturbations. A key technical ingredient for our certification guarantees is a fast algorithm with provable guarantees based on the multiplicative weights update method to provide upper bounds on the above matrix norm. Our algorithmic guarantees improve upon the state of the art for this problem, and may be of independent interest.


title: Colossal magnetoresistance in a nonsymmorphic antiferromagnetic insulator

abstract: Here we investigate antiferromagnetic Eu$_{5}$In$_{2}$Sb$_{6}$, a nonsymmorphic Zintl phase. Our electrical transport data show that Eu$_{5}$In$_{2}$Sb$_{6}$ is remarkably insulating and exhibits an exceptionally large negative magnetoresistance, which is consistent with the presence of magnetic polarons. From {\it ab initio} calculations, the paramagnetic state of Eu$_{5}$In$_{2}$Sb$_{6}$ is a topologically nontrivial semimetal within the generalized gradient approximation (GGA), whereas an insulating state with trivial topological indices is obtained using a modified Becke-Johnson potential. Notably, GGA+U calculations suggest that the antiferromagnetic phase of Eu$_{5}$In$_{2}$Sb$_{6}$ may host an axion insulating state. Our results provide important feedback for theories of topological classification and highlight the potential of realizing clean magnetic narrow-gap semiconductors in Zintl materials.


title: Prediction-Centric Learning of Independent Cascade Dynamics from Partial Observations

abstract: Spreading processes play an increasingly important role in modeling for diffusion networks, information propagation, marketing and opinion setting. We address the problem of learning of a spreading model such that the predictions generated from this model are accurate and could be subsequently used for the optimization, and control of diffusion dynamics. We focus on a challenging setting where full observations of the dynamics are not available, and standard approaches such as maximum likelihood quickly become intractable for large network instances. We introduce a computationally efficient algorithm, based on a scalable dynamic message-passing approach, which is able to learn parameters of the effective spreading model given only limited information on the activation times of nodes in the network. The popular Independent Cascade model is used to illustrate our approach. We show that tractable inference from the learned model generates a better prediction of marginal probabilities compared to the original model. We develop a systematic procedure for learning a mixture of models which further improves the prediction quality.


title: Fast Global Convergence of Natural Policy Gradient Methods with Entropy Regularization

abstract: Natural policy gradient (NPG) methods are among the most widely used policy optimization algorithms in contemporary reinforcement learning. This class of methods is often applied in conjunction with entropy regularization -- an algorithmic scheme that encourages exploration -- and is closely related to soft policy iteration and trust region policy optimization. Despite the empirical success, the theoretical underpinnings for NPG methods remain limited even for the tabular setting. This paper develops $\textit{non-asymptotic}$ convergence guarantees for entropy-regularized NPG methods under softmax parameterization, focusing on discounted Markov decision processes (MDPs). Assuming access to exact policy evaluation, we demonstrate that the algorithm converges linearly -- or even quadratically once it enters a local region around the optimal policy -- when computing optimal value functions of the regularized MDP. Moreover, the algorithm is provably stable vis-\`a-vis inexactness of policy evaluation. Our convergence results accommodate a wide range of learning rates, and shed light upon the role of entropy regularization in enabling fast convergence.


title: Graph Structure of Neural Networks

abstract: Neural networks are often represented as graphs of connections between neurons. However, despite their wide use, there is currently little understanding of the relationship between the graph structure of the neural network and its predictive performance. Here we systematically investigate how does the graph structure of neural networks affect their predictive performance. To this end, we develop a novel graph-based representation of neural networks called relational graph, where layers of neural network computation correspond to rounds of message exchange along the graph structure. Using this representation we show that: (1) a "sweet spot" of relational graphs leads to neural networks with significantly improved predictive performance; (2) neural network's performance is approximately a smooth function of the clustering coefficient and average path length of its relational graph; (3) our findings are consistent across many different tasks and datasets; (4) the sweet spot can be identified efficiently; (5) top-performing neural networks have graph structure surprisingly similar to those of real biological neural networks. Our work opens new directions for the design of neural architectures and the understanding on neural networks in general.


title: Mass Ratios of Merging Double Neutron Stars as Implied by the Milky Way Population

abstract: Of the seven known double neutron stars (DNS) with precisely measure masses in the Milky Way that will merge within a Hubble time, all but one has a mass ratio, $q$, close to unity. Recently, precise measurements of three post-Keplerian parameters in the DNS J1913$+$1102 constrain this system to have a significantly non-unity mass ratio of 0.78$\pm$0.03. One may be tempted to conclude that approximately one out of seven (14\%) DNS mergers detected by gravitational wave observatories will have mass ratios significantly different from unity. However J1913$+$1102 has a relatively long merger time of 470 Myr. We show that when merger times and observational biases are taken into account, the population of Galactic DNSs imply that $\simeq98\%$ of all merging DNSs will have $q>$0.9. We then apply two separate fitting formulas informed by 3D hydrodynamic simulations of DNS mergers to our results on Galactic DNS masses, finding that either $\simeq$0.004 ${M_{\odot}}$ or $\simeq$0.010 ${M_{\odot}}$ of material will be ejected at merger, depending on which formula is used. These ejecta masses have implications for both the peak bolometric luminosities of electromagnetic counterparts (which we find to be $\sim$10$^{41}$ erg s$^{-1}$) as well as the $r$-process enrichment of the Milky Way.


title: Noncollinear antiferromagnetic order in the buckled honeycomb lattice of magnetoelectric Co4Ta2O9 determined by single-crystal neutron diffraction

abstract: Co4Ta2O9 exhibits a three-dimensional magnetic lattice based on the buckled honeycomb motif. It shows unusual magnetoelectric effects, including the sign change and non-linearity. These effects cannot be understood without the detailed knowledge of the magnetic structure. Herein, we report neutron diffraction and direction-dependent magnetic susceptibility measurements on Co4Ta2O9 single crystals. Below 20.3 K, we find a long-range antiferromagnetic order in the alternating buckled and flat honeycomb layers of Co2+ ions stacked along the c axis. Within experimental accuracy, the magnetic moments lie in the ab plane. They form a canted antiferromagnetic structure with a tilt angle of ~ 14 degrees at 15 K in the buckled layers, while the magnetic moments in each flat layer are collinear. This is directly evidenced by a finite (0, 0, 3) magnetic Bragg peak intensity, which would be absent in the collinear magnetic order. The magnetic space group is C2'/c. It is different from the previously reported C2/c' group, also found in the isostructural Co4Nb2O9. The revised magnetic structure successfully explains the major features of the magnetoelectric tensor of Co4Ta2O9 within the framework of the spin-flop model.


title: The role of clouds on the depletion of methane and water dominance in the transmission spectra of irradiated exoplanets

abstract: Observations suggest an abundance of water and paucity of methane in the majority of observed exoplanetary atmospheres. We isolate the effect of atmospheric processes to investigate possible causes. Previously, we studied the effect of effective temperature, surface gravity, metallicity, carbon-to-oxygen ratio, and stellar type assuming cloud-free thermochemical equilibrium and disequilibrium chemistry. However, under these assumptions, methane remains a persisting spectral feature in the transmission spectra of exoplanets over a certain parameter space, the Methane Valley. In this work we investigate the role of clouds on this domain and we find that clouds change the spectral appearance of methane in two direct ways: 1) by heating-up the photosphere of colder planets, and 2) by obscuring molecular features. The presence of clouds also affects methane features indirectly: 1) cloud heating results in more evaporation of condensates and hence releases additional oxygen, causing water dominated spectra of colder carbon-poor exoplanets, and 2) HCN/CO production results in a suppression of depleted methane features by these molecules. The presence of HCN/CO and a lack of methane could be an indication of cloud formation on hot exoplanets. Cloud heating can also deplete ammonia. Therefore, a simultaneous depletion of methane and ammonia is not unique to photochemical processes. We propose that the best targets for methane detection are likely to be massive but smaller planets with a temperature around 1450 K orbiting colder stars. We also construct Spitzer synthetic color-maps and find that clouds can explain some of the high contrast observations by IRAC's channel 1 and 2.


title: HOBFLOPS CNNs: Hardware Optimized Bitslice-Parallel Floating-Point Operations for Convolutional Neural Networks

abstract: Convolutional neural networks (CNNs) are typically trained using 16- or 32-bit floating-point (FP) and researchers show that low-precision floating-point (FP) can be highly effective for inference. Low-precision FP can be implemented in field programmable gate array (FPGA) and application-specific integrated circuit (ASIC) accelerators, but existing processors do not generally support custom precision FP. We propose hardware optimized bitslice-parallel floating-point operators (HOBFLOPS), a method of generating efficient custom-precision emulated bitslice-parallel software FP arithmetic. We generate custom-precision FP routines optimized using a hardware synthesis design flow to create circuits. We provide standard cell libraries matching the bitwise operations on the target microprocessor architecture, and a code-generator to translate the hardware circuits to bitslice software equivalents. We exploit bitslice parallelism to create a very wide (32-512 element) vectorized convolutional neural network (CNN) convolution. Hardware optimized bitslice-parallel floating-point operators (HOBFLOPS) multiply-accumulate (MAC) performance in CNN convolution on Arm and Intel processors are compared to Berkeley's SoftFP16 equivalent MAC. HOBFLOPS16 outperforms SoftFP16 by 8x on Intel AVX512. HOBFLOPS offers arbitrary-precision FP with custom range and precision e.g., HOBFLOPS9 performs at 6x the performance of HOBFLOPS16 on Arm Neon. HOBFLOPS allows researchers to prototype different levels of custom FP precision in the arithmetic of software CNN accelerators. Furthermore, HOBFLOPS fast custom-precision FP CNNs may be valuable in cases where memory bandwidth is limited.


title: Uncertainty relations in terms of the Gini index for finite quantum systems

abstract: Lorenz values and the Gini index are popular quantities in Mathematical Economics, and are used here in the context of quantum systems with finite-dimensional Hilbert space. They quantify the uncertainty in the probability distribution related to an orthonormal basis. It is shown that Lorenz values are superadditive functions and the Gini indices are subadditive functions. The supremum over all density matrices of the sum of the two Gini indices with respect to position and momentum states, is used to define an uncertainty coefficient which quantifies the uncertainty in the quantum system. It is shown that the uncertainty coefficient is positive, and an upper bound for it is given. Various examples demonstrate these ideas.


title: FocusLiteNN: High Efficiency Focus Quality Assessment for Digital Pathology

abstract: Out-of-focus microscopy lens in digital pathology is a critical bottleneck in high-throughput Whole Slide Image (WSI) scanning platforms, for which pixel-level automated Focus Quality Assessment (FQA) methods are highly desirable to help significantly accelerate the clinical workflows. Existing FQA methods include both knowledge-driven and data-driven approaches. While data-driven approaches such as Convolutional Neural Network (CNN) based methods have shown great promises, they are difficult to use in practice due to their high computational complexity and lack of transferability. Here, we propose a highly efficient CNN-based model that maintains fast computations similar to the knowledge-driven methods without excessive hardware requirements such as GPUs. We create a training dataset using FocusPath which encompasses diverse tissue slides across nine different stain colors, where the stain diversity greatly helps the model to learn diverse color spectrum and tissue structures. In our attempt to reduce the CNN complexity, we find with surprise that even trimming down the CNN to the minimal level, it still achieves a highly competitive performance. We introduce a novel comprehensive evaluation dataset, the largest of its kind, annotated and compiled from TCGA repository for model assessment and comparison, for which the proposed method exhibits superior precision-speed trade-off when compared with existing knowledge-driven and data-driven FQA approaches.


title: A unified machine learning approach to time series forecasting applied to demand at emergency departments

abstract: There were 25.6 million attendances at Emergency Departments (EDs) in England in 2019 corresponding to an increase of 12 million attendances over the past ten years. The steadily rising demand at EDs creates a constant challenge to provide adequate quality of care while maintaining standards and productivity. Managing hospital demand effectively requires an adequate knowledge of the future rate of admission. Using 8 years of electronic admissions data from two major acute care hospitals in London, we develop a novel ensemble methodology that combines the outcomes of the best performing time series and machine learning approaches in order to make highly accurate forecasts of demand, 1, 3 and 7 days in the future. Both hospitals face an average daily demand of 208 and 106 attendances respectively and experience considerable volatility around this mean. However, our approach is able to predict attendances at these emergency departments one day in advance up to a mean absolute error of +/- 14 and +/- 10 patients corresponding to a mean absolute percentage error of 6.8% and 8.6% respectively. Our analysis compares machine learning algorithms to more traditional linear models. We find that linear models often outperform machine learning methods and that the quality of our predictions for any of the forecasting horizons of 1, 3 or 7 days are comparable as measured in MAE. In addition to comparing and combining state-of-the-art forecasting methods to predict hospital demand, we consider two different hyperparameter tuning methods, enabling a faster deployment of our models without compromising performance. We believe our framework can readily be used to forecast a wide range of policy relevant indicators.


title: Lossless Compression of Structured Convolutional Models via Lifting

abstract: Lifting is an efficient technique to scale up graphical models generalized to relational domains by exploiting the underlying symmetries. Concurrently, neural models are continuously expanding from grid-like tensor data into structured representations, such as various attributed graphs and relational databases. To address the irregular structure of the data, the models typically extrapolate on the idea of convolution, effectively introducing parameter sharing in their, dynamically unfolded, computation graphs. The computation graphs themselves then reflect the symmetries of the underlying data, similarly to the lifted graphical models. Inspired by lifting, we introduce a simple and efficient technique to detect the symmetries and compress the neural models without loss of any information. We demonstrate through experiments that such compression can lead to significant speedups of structured convolutional models, such as various Graph Neural Networks, across various tasks, such as molecule classification and knowledge-base completion.


title: Impact of PAH photodissociation on the formation of small hydrocarbons in the Orion Bar and the Horsehead PDRs

abstract: We study whether polycyclic aromatic hydrocarbons (PAHs) can be a weighty source of small hydrocarbons in photo-dissociation regions (PDRs). We modeled the evolution of 20 specific PAH molecules in terms of dehydrogenation and destruction of the carbon skeleton under the physical conditions of two well-studied PDRs, the Orion Bar and the Horsehead nebula which represent prototypical examples of PDRs irradiated by "high" and "low" ultraviolet radiation field. PAHs are described as microcanonical systems. The acetylene molecule is considered as the main carbonaceous fragment of the PAH dissociation as it follows from laboratory experiments and theory. We estimated the rates of acetylene production in gas phase chemical reactions and compared them with the rates of the acetylene production through the PAH dissociation. It is found that the latter rates can be higher than the former rates in the Orion Bar at $A_{\rm V}<1$ and also at $A_{\rm V}>3.5$. In the Horsehead nebula, the chemical reactions provide more acetylene than the PAH dissociation. The produced acetylene participate in the reactions of the formation of small hydrocarbons (C$_2$H, C$_3$H, C$_3$H$^{+}$, C$_3$H$_2$, C$_4$H). Acetylene production via the PAH destruction may increase the abundances of small hydrocarbons produced in gas phase chemical reactions in the Orion Bar only at $A_{\rm V}>3.5$. In the Horsehead nebula, the contribution of PAHs to the abundances of the small hydrocarbons is negligible. We conclude that the PAHs are not a major source of small hydrocarbons in both PDRs except some locations in the Orion Bar.


title: Conformal mapping in linear time

abstract: Given any $\epsilon >0$ and any planar region $\Omega$ bounded by a simple n-gon $P$ we construct a ($1 + \epsilon)$-quasiconformal map between $\Omega$ and the unit disk in time $C(\epsilon)n$. One can take $ C(\epsilon) = C + C \log (1/\epsilon) \log \log (1/\epsilon)$.


title: Towards causal benchmarking of bias in face analysis algorithms

abstract: Measuring algorithmic bias is crucial both to assess algorithmic fairness, and to guide the improvement of algorithms. Current methods to measure algorithmic bias in computer vision, which are based on observational datasets, are inadequate for this task because they conflate algorithmic bias with dataset bias. To address this problem we develop an experimental method for measuring algorithmic bias of face analysis algorithms, which manipulates directly the attributes of interest, e.g., gender and skin tone, in order to reveal causal links between attribute variation and performance change. Our proposed method is based on generating synthetic ``transects'' of matched sample images that are designed to differ along specific attributes while leaving other attributes constant. A crucial aspect of our approach is relying on the perception of human observers, both to guide manipulations, and to measure algorithmic bias. Besides allowing the measurement of algorithmic bias, synthetic transects have other advantages with respect to observational datasets: they sample attributes more evenly allowing for more straightforward bias analysis on minority and intersectional groups, they enable prediction of bias in new scenarios, they greatly reduce ethical and legal challenges, and they are economical and fast to obtain, helping make bias testing affordable and widely available. We validate our method by comparing it to a study that employs the traditional observational method for analyzing bias in gender classification algorithms. The two methods reach different conclusions. While the observational method reports gender and skin color biases, the experimental method reveals biases due to gender, hair length, age, and facial hair.


title: Inverse Cubic Iteration

abstract: There are thousands of papers on rootfinding for nonlinear scalar equations. Here is one more, to talk about an apparently new method, which I call ``Inverse Cubic Iteration'' (ICI) in analogy to the Inverse Quadratic Iteration in Richard Brent's zeroin method. The possibly new method is based on a cubic blend of tangent-line approximations for the inverse function. We rewrite this iteration for numerical stability as an average of two Newton steps and a secant step: only one new function evaluation and derivative evaluation is needed for each step. The total cost of the method is therefore only trivially more than Newton's method, and we will see that it has order $1+\sqrt{3} = 2.732...$, thus ensuring that to achieve a given accuracy it usually takes fewer steps than Newton's method while using essentially the same effort per step.


title: A physical model for radiative, convective dusty disk in AGN

abstract: An accretion disk in an Active Galactic Nucleus (AGN) harbors and shields dust from external illumination: at the mid-plane of the disk around a $M_{{\rm BH}}=10^{7}M_{\odot}$ black hole, dust can exist at $0.1$pc from the black hole, compared to 0.5pc outside of the disk. We construct a physical model of a disk region approximately located between the radius of dust sublimation at the disk mid-plane and the radius at which dust sublimes at the disk surface. Our main conclusion is that for a wide range of model parameters such as local accretion rate and/or opacity, the accretion disk's own radiation pressure on dust significantly influences its vertical structure. In addition to being highly convective, such a disk can transform from geometrically thin to slim. Our model fits into the narrative of a "failed wind" scenario of Czerny & Hryniewicz (2011) and the "compact torus" model of Baskin & Laor (2018), incorporating them as variations of the radiative dusty disk model.


title: The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300 au with SPHERE

abstract: The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies.


title: Scattering Equations in AdS: Scalar Correlators in Arbitrary Dimensions

abstract: We introduce a bosonic ambitwistor string theory in AdS space. Even though the theory is anomalous at the quantum level, one can nevertheless use it in the classical limit to derive a novel formula for correlation functions of boundary CFT operators in arbitrary space-time dimensions. The resulting construction can be treated as a natural extension of the CHY formalism for the flat-space S-matrix, as it similarly expresses tree-level amplitudes in AdS as integrals over the moduli space of Riemann spheres with punctures. These integrals localize on an operator-valued version of scattering equations, which we derive directly from the ambitwistor string action on a coset manifold. As a testing ground for this formalism we focus on the simplest case of ambitwistor string coupled to two current algebras, which gives bi-adjoint scalar correlators in AdS. In order to evaluate them directly, we make use of a series of contour deformations on the moduli space of punctured Riemann spheres and check that the result agrees with tree level Witten diagram computations to all multiplicity. We also initiate the study of eigenfunctions of scattering equations in AdS, which interpolate between conformal partial waves in different OPE channels, and point out a connection to an elliptic deformation of the Calogero-Sutherland model.