title: Long-period variables in the Gaia era

abstract: The second Gaia data release (DR2, spring 2018) included a unique all-sky catalogue of large-amplitude long-period variables (LPVs) containing Miras and semi-regular variables. These stars are on the Asymptotic Giant Branch (AGB), and are characterized by high luminosity, changing surface composition, and intense mass loss, that make them of paramount importance for stellar, galactic, and extra-galactic studies. An initial investigation of LPVs in the Large Magellanic Cloud (LMC) from the DR2 catalog of LPVs has revealed the possibility to disentangle O-rich and C-rich stars using a combination of optical Gaia and infrared 2MASS photometry. The so-called Gaia-2MASS diagram constructed to achieve this has further been shown to enable the identification of sub-groups of AGB stars among the O-rich and C-rich LPVs. Here, we extend this initial study of the Gaia-2MASS diagram to the Small Magellanic Cloud and the Galaxy, and use a variability amplitude proxy to identify LPVs from the full Gaia DR2 archive. We show that the remarkable properties found in the LMC also apply to these other stellar systems. Interesting features, moreover, emerge as a result of the different metallicities between the three stellar environments, which we highlight in this exploratory presentation of Gaia's potential to study stellar populations harboring LPVs. Finally, we look ahead to the future, and highlight the power of the exploitation of Gaia RP spectra for the identification of carbon stars using solely Gaia data in forthcoming data releases, as revealed in an Image of the Week published by the Gaia consortium on the European Space Agency's web site. These proceedings include three animated images that can be used as outreach material.


title: Hodge-theoretic analysis on manifolds with boundary, heatable currents, and Onsager's conjecture in fluid dynamics

abstract: We use Hodge theory and functional analysis to develop a clean approach to heat flows and Onsager's conjecture on Riemannian manifolds with boundary, where the weak solution lies in the trace-critical Besov space $B_{3,1}^{\frac{1}{3}}$. We also introduce heatable currents as the natural analogue to tempered distributions and justify their importance in Hodge theory.


title: Gravitational Waves from Black Holes in Merging Ultra-Dwarf Galaxies

abstract: The origin of the black-hole:black-hole mergers discovered through gravitational waves with for example the LIGO/Virgo collaboration are a mystery. We investigate the idea that some of these black holes originate from the centers of extremely low-mass ultra-dwarf galaxies that have merged together in the distant past at $z>1$. Extrapolating the central black hole to stellar mass ratio suggests that the black holes in these mergers could arise from galaxies of masses $\sim 10^{5} - 10^{6}$ M$_{\odot}\,$. We investigate whether these galaxies merge enough, or too much, to be consistent with the observed GW rate of $\sim 9.7-101$ Gpc$^{-3}$ yr$^{-1}$ using the latest LIGO/Virgo results. We show that in the nearby universe the merger rate and number densities of ultra-dwarf galaxies are too low, by an order or magnitude, to produce these black hole mergers. However, by considering that the merger fraction, merger-time scales, and the number densities of low-mass galaxies all conspire at $z>1-1.5$ to increase the merger rate for these galaxies at higher redshifts we argue that it is possible that some of the observed GW events arise from BHs in the centers of low-mass galaxies. The major uncertainty in this calculation is the dynamical time-scales for black holes in low-mass galaxies. Our results however suggest a very long BH merger time-scale of 4-7 Gyr, consistent with an extended black hole merger history. Further simulations are needed to verify this possibility, however our theory can be tested by searching for host galaxies of gravitational wave events. Results from these searches would will put limits on dwarf galaxy mergers and/or the presence and formation mechanisms of black holes through PopIII stars in the lowest mass galaxies.


title: Continuous changes of variables and the Magnus expansion

abstract: In this paper, we are concerned with a formulation of Magnus and Floquet-Magnus expansions for general nonlinear differential equations. To this aim, we introduce suitable continuous variable transformations generated by operators. As an application of the simple formulas so-obtained, we explicitly compute the first terms of the Floquet-Magnus expansion for the Van der Pol oscillator and the nonlinear Schr\"odinger equation on the torus.

categories: math-ph math.CA math.MP physics.class-ph


title: Warfarin dose estimation on multiple datasets with automated hyperparameter optimisation and a novel software framework

abstract: Warfarin is an effective preventative treatment for arterial and venous thromboembolism, but requires individualised dosing due to its narrow therapeutic range and high individual variation. Many machine learning techniques have been demonstrated in this domain. This study evaluated the accuracy of the most promising algorithms on the International Warfarin Pharmacogenetics Consortium dataset and a novel clinical dataset of South African patients. Support vectors and linear regression were amongst the top performers in both datasets and performed comparably to recent stacked ensemble approaches, whilst neural networks were one of the worst performers in both datasets. We also introduced genetic programming to automatically optimise model architectures and hyperparameters without human guidance. Remarkably, the generated models were found to match the performance of the best models hand-crafted by human experts. Finally, we present a novel software framework (Warfit-learn) for warfarin dosing research. It leverages the most successful techniques in preprocessing, imputation, and parallel evaluation, with the goal of accelerating research and making results in this domain more reproducible.


title: Performance Boundary Identification for the Evaluation of Automated Vehicles using Gaussian Process Classification

abstract: Safety is an essential aspect in the facilitation of automated vehicle deployment. Current testing practices are not enough, and going beyond them leads to infeasible testing requirements, such as needing to drive billions of kilometres on public roads. Automated vehicles are exposed to an indefinite number of scenarios. Handling of the most challenging scenarios should be tested, which leads to the question of how such corner cases can be determined. We propose an approach to identify the performance boundary, where these corner cases are located, using Gaussian Process Classification. We also demonstrate the classification on an exemplary traffic jam approach scenario, showing that it is feasible and would lead to more efficient testing practices.


title: Nature of the Doping-Driven Metal-Insulator Transition in Correlated Electron Systems with strong Hund's Exchange Coupling

abstract: We study the doping-driven Mott metal-insulator transition for multi-orbital Hubbard models with Hund's exchange coupling at finite temperatures. As in the single-orbital Hubbard model, the transition is of first-order within dynamical mean field theory, with a coexistence region where two solutions can be stabilized. We find, that in the presence of finite Hund's coupling, the insulating phase is connected to a badly metallic phase, which extends to surprisingly large dopings. While fractional power-law behavior of the self-energies on the Matsubara axis is found on both sides of the transition, a regime with frozen local moments develops only on the branch connected to the insulating phase.


title: Upper bounds for the regularity of symbolic powers of certain classes of edge ideals

abstract: Let $G$ be a finite simple graph and $I(G)$ denote the corresponding edge ideal in a polynomial ring over a field $\mathbb{K}$. In this paper, we obtain upper bounds for the Castelnuovo-Mumford regularity of symbolic powers of certain classes of edge ideals. We also prove that for several classes of graphs, the regularity of symbolic powers of their edge ideals coincides with that of their ordinary powers.


title: The sounds of science: a symphony for many instruments and voices

abstract: This paper is a celebration of the frontiers of science. Goodenough, the maestro who transformed energy usage and technology through the invention of the lithium ion battery, opens the programme, reflecting on the ultimate limits of battery technology. This applied theme continues through the subsequent pieces on energy related topics (the sodium ion battery and artificial fuels, by Mansson) and the ultimate challenge for 3 dimensional printing the eventual production of life, by Atala. A passage by Alexander follows, reflecting on a related issue: How might an artificially produced human being behave? Next comes a consideration of consiousness and free will by Allen and Lidstrom. Further voices and new instruments enter as Bowen, Mauranyapin and Madsen discuss whether dynamical processes of single molecules might be observed in their native state. The exploitation of chaos in science and technology, applications of Bose Einstein condensates and a consideration of the significance of entropy follow in pieces by Reichl, Rasel and Allen, respectively. Katsnelson and Koonin then discuss the potential generalisation of thermodynamic concepts in the context of biological evolution. Entering with the music of the cosmos, Yasskin discusses whether we might be able to observe torsion in the geometry of the universe. The crescendo comes with the crisis of singularities, their nature and whether they can be resolved through quantum effects, in the composition of Coley. The climax is Krenn, Melvin and Zeilinger consideration of how computer code can be autonomously surprising and creative. In a harmonious counterpoint, Yampolskiy concludes that such code is not yet able to take responsibility for coauthoring a paper.


title: On The Structure of Dyck Languages

abstract: We prove that the closure of the one-sided Dyck language in a free monoid is a two-sided Dyck language.


title: Abelian-square factors and binary words

abstract: In this work, we affirm the conjecture proposed by Gabriele Fici and Filippo Mignosi at the 10th Conference on Combinatorics on Words.


title: The Selfish Higgs

abstract: We propose a mechanism to solve the Higgs naturalness problem through a cosmological selection process. The discharging of excited field configurations through membrane nucleation leads to discrete jumps of the cosmological constant and the Higgs mass, which vary in a correlated way. The resulting multitude of universes are all empty, except for those in which the cosmological constant and the Higgs mass are both nearly vanishing. Only under these critical conditions can inflation be activated and create a non-empty universe.


title: Analysis and data-driven reconstruction of bivariate jump-diffusion processes

abstract: We introduce the bivariate jump-diffusion process, comprising two-dimensional diffusion and two-dimensional jumps, that can be coupled to one another. We present a data-driven, non-parametric estimation procedure of higher-order (up to 8) Kramers-Moyal coefficients that allows one to reconstruct relevant aspects of the underlying jump-diffusion processes and to recover the underlying parameters. The procedure is validated with numerically integrated data using synthetic bivariate time series from continuous and discontinuous processes. We further evaluate the possibility of estimating the parameters of the jump-diffusion model via data-driven analyses of the higher-order Kramers-Moyal coefficients, and the limitations arising from the scarcity of points in the data or disproportionate parameters in the system.


title: Repeated Measurements with Minimally Destructive Partial-Transfer Absorption Imaging

abstract: We demonstrate partial-transfer absorption imaging as a technique for repeatedly imaging an ultracold atomic ensemble with minimal perturbation. We prepare an atomic cloud in a state that is dark to the imaging light. We then use a microwave pulse to coherently transfer a small fraction of the ensemble to a bright state, which we image using in situ absorption imaging. The amplitude or duration of the microwave pulse controls the fractional transfer from the dark to the bright state. For small transfer fractions, we can image the atomic cloud up to 50 times before it is depleted. As a sample application, we repeatedly image an atomic cloud oscillating in a dipole trap to measure the trap frequency.


title: Torus Spectroscopy of the Gross-Neveu-Yukawa Quantum Field Theory: Free Dirac versus Chiral Ising Fixed Point

abstract: We establish the universal torus low-energy spectra at the free Dirac fixed point and at the strongly coupled chiral Ising fixed point and their subtle crossover behaviour in the Gross-Neuveu-Yukawa field theory with ${n_\text{D}=4}$ component Dirac spinors in $D=(2+1)$ dimensions. These fixed points and the field theories are directly relevant for the long-wavelength physics of certain interacting Dirac systems, such as repulsive spinless fermions on the honeycomb lattice or $\pi$-flux square lattice. The torus energy spectrum has been shown previously to serve as a characteristic fingerprint of relativistic fixed points and is a powerful tool to discriminate quantum critical behaviour in numerical simulations. Here, we use a combination of exact diagonalization and quantum Monte Carlo simulations of strongly interacting fermionic lattice models, to compute the critical torus energy spectrum on finite-size clusters with periodic boundaries and extrapolate them to the thermodynamic limit. Additionally, we compute the torus energy spectrum analytically using the perturbative expansion in ${\epsilon = 4 - D}$, which is in good agreement with the numerical results, thereby validating the presence of the chiral Ising fixed point in the lattice models at hand. We show that the strong interaction between the spinor field and the scalar order-parameter field strongly influences the critical torus energy spectrum and we observe prominent multiplicity features related to an emergent symmetry predicted from the quantum field theory. Building on these results we are able to address the subtle crossover physics of the low-energy spectrum flowing from the chiral Ising fixed point to the Dirac fixed point, and analyze earlier flawed attempts to extract Fermi velocity renormalizations from the low-energy spectrum.


title: Universal behavior in non stationary Mean Field Games

abstract: Mean Field Games provide a powerful framework to analyze the dynamics of a large number of controlled objects in interaction. Though these models are much simpler than the underlying differential games they describe in some limit, their behavior is still far from being fully understood. When the system is confined, a notion of "ergodic state" has been introduced that characterizes most of the dynamics for long optimization times. Here we consider a class of models without such an ergodic state, and show the existence of a scaling solution that plays similar role. Its universality and scaling behavior can be inferred from a mapping to an electrostatic problem.


title: Online Inference and Detection of Curbs in Partially Occluded Scenes with Sparse LIDAR

abstract: Road boundaries, or curbs, provide autonomous vehicles with essential information when interpreting road scenes and generating behaviour plans. Although curbs convey important information, they are difficult to detect in complex urban environments (in particular in comparison to other elements of the road such as traffic signs and road markings). These difficulties arise from occlusions by other traffic participants as well as changing lighting and/or weather conditions. Moreover, road boundaries have various shapes, colours and structures while motion planning algorithms require accurate and precise metric information in real-time to generate their plans. In this paper, we present a real-time LIDAR-based approach for accurate curb detection around the vehicle (360 degree). Our approach deals with both occlusions from traffic and changing environmental conditions. To this end, we project 3D LIDAR pointcloud data into 2D bird's-eye view images (akin to Inverse Perspective Mapping). These images are then processed by trained deep networks to infer both visible and occluded road boundaries. Finally, a post-processing step filters detected curb segments and tracks them over time. Experimental results demonstrate the effectiveness of the proposed approach on real-world driving data. Hence, we believe that our LIDAR-based approach provides an efficient and effective way to detect visible and occluded curbs around the vehicles in challenging driving scenarios.


title: Monocular 3D Sway Tracking for Assessing Postural Instability in Cerebral Hypoperfusion During Quiet Standing

abstract: Postural instability is prevalent in aging and neurodegenerative disease, decreasing quality of life and independence. Quantitatively monitoring balance control is important for assessing treatment efficacy and rehabilitation progress. However, existing technologies for assessing postural sway are complex and expensive, limiting their widespread utility. Here, we propose a monocular imaging system capable of assessing sub-millimeter 3D sway dynamics during quiet standing. Two anatomical targets with known feature geometries were placed on the lumbar and shoulder. Upper and lower trunk 3D kinematic motion was automatically assessed from a set of 2D frames through geometric feature tracking and an inverse motion model. Sway was tracked in 3D and compared between control and hypoperfusion conditions in 14 healthy young adults. The proposed system demonstrated high agreement with a commercial motion capture system (error $1.5 \times 10^{-4}~\text{mm}$, [$-0.52$, $0.52$]). Between-condition differences in sway dynamics were observed in anterior-posterior sway during early and mid stance, and medial-lateral sway during mid stance commensurate with decreased cerebral perfusion, followed by recovered sway dynamics during late stance with cerebral perfusion recovery. This inexpensive single-camera system enables quantitative 3D sway monitoring for assessing neuromuscular balance control in weakly constrained environments.


title: Method of moments for 3-D single particle ab initio modeling with non-uniform distribution of viewing angles

abstract: Single-particle reconstruction in cryo-electron microscopy (cryo-EM) is an increasingly popular technique for determining the 3-D structure of a molecule from several noisy 2-D projections images taken at unknown viewing angles. Most reconstruction algorithms require a low-resolution initialization for the 3-D structure, which is the goal of ab initio modeling. Suggested by Zvi Kam in 1980, the method of moments (MoM) offers one approach, wherein low-order statistics of the 2-D images are computed and a 3-D structure is estimated by solving a system of polynomial equations. Unfortunately, Kam's method suffers from restrictive assumptions, most notably that viewing angles should be distributed uniformly. Often unrealistic, uniformity entails the computation of higher-order correlations, as in this case first and second moments fail to determine the 3-D structure. In the present paper, we remove this hypothesis, by permitting an unknown, non-uniform distribution of viewing angles in MoM. Perhaps surprisingly, we show that this case is statistically easier than the uniform case, as now first and second moments generically suffice to determine low-resolution expansions of the molecule. In the idealized setting of a known, non-uniform distribution, we find an efficient provable algorithm inverting first and second moments. For unknown, non-uniform distributions, we use non-convex optimization methods to solve for both the molecule and distribution.


title: Quantum and Classical Algorithms for Approximate Submodular Function Minimization

abstract: Submodular functions are set functions mapping every subset of some ground set of size $n$ into the real numbers and satisfying the diminishing returns property. Submodular minimization is an important field in discrete optimization theory due to its relevance for various branches of mathematics, computer science and economics. The currently fastest strongly polynomial algorithm for exact minimization [LSW15] runs in time $\widetilde{O}(n^3 \cdot \mathrm{EO} + n^4)$ where $\mathrm{EO}$ denotes the cost to evaluate the function on any set. For functions with range $[-1,1]$, the best $\epsilon$-additive approximation algorithm [CLSW17] runs in time $\widetilde{O}(n^{5/3}/\epsilon^{2} \cdot \mathrm{EO})$. In this paper we present a classical and a quantum algorithm for approximate submodular minimization. Our classical result improves on the algorithm of [CLSW17] and runs in time $\widetilde{O}(n^{3/2}/\epsilon^2 \cdot \mathrm{EO})$. Our quantum algorithm is, up to our knowledge, the first attempt to use quantum computing for submodular optimization. The algorithm runs in time $\widetilde{O}(n^{5/4}/\epsilon^{5/2} \cdot \log(1/\epsilon) \cdot \mathrm{EO})$. The main ingredient of the quantum result is a new method for sampling with high probability $T$ independent elements from any discrete probability distribution of support size $n$ in time $O(\sqrt{Tn})$. Previous quantum algorithms for this problem were of complexity $O(T\sqrt{n})$.


title: Analytic and Probabilistic Problems in Discrete Geometry

abstract: The thesis concentrates on two problems in discrete geometry, whose solutions are obtained by analytic, probabilistic and combinatoric tools. The first chapter deals with the strong polarization problem. This states that for any sequence $u_1,\dots, u_n$ of norm 1 vectors in a real Hilbert space $\mathscr H$, there exists a unit vector $v \in \mathscr H$, such that $$ \sum \frac{1}{\langle u_i, v \rangle^2} \leq n^2. $$ The 2-dimensional case is proved by complex analytic methods. For the higher dimensional extremal cases, we prove a tensorisation result that is similar to F. John's theorem about characterisation of ellipsoids of maximal volume. From this, we deduce that the only full dimensional locally extremal system is the orthonormal system. We also obtain the same result for the weaker, original polarization problem. The second chapter investigates a problem in probabilistic geometry. Take $n$ independent, uniform random points in a triangle $T$. Convex chains between two fixed vertices of $T$ are defined naturally. Let $L_n$ denote the maximal size of a convex chain. We prove that the expectation of $L_n$ is asymptotically $\alpha \, n^{1/3}$, where $\alpha$ is a constant between 1.5 and 3.5 -- we conjecture that the correct value is 3. We also prove strong concentration results for $L_n$, which, in turn, imply a limit shape result for the longest convex chains.


title: Single Image Super-Resolution via CNN Architectures and TV-TV Minimization

abstract: Super-resolution (SR) is a technique that allows increasing the resolution of a given image. Having applications in many areas, from medical imaging to consumer electronics, several SR methods have been proposed. Currently, the best performing methods are based on convolutional neural networks (CNNs) and require extensive datasets for training. However, at test time, they fail to impose consistency between the super-resolved image and the given low-resolution image, a property that classic reconstruction-based algorithms naturally enforce in spite of having poorer performance. Motivated by this observation, we propose a new framework that joins both approaches and produces images with superior quality than any of the prior methods. Although our framework requires additional computation, our experiments on Set5, Set14, and BSD100 show that it systematically produces images with better peak signal to noise ratio (PSNR) and structural similarity (SSIM) than the current state-of-the-art CNN architectures for SR.


title: Adaptive Pricing in Insurance: Generalized Linear Models and Gaussian Process Regression Approaches

abstract: We study the application of dynamic pricing to insurance. We view this as an online revenue management problem where the insurance company looks to set prices to optimize the long-run revenue from selling a new insurance product. We develop two pricing models: an adaptive Generalized Linear Model (GLM) and an adaptive Gaussian Process (GP) regression model. Both balance between exploration, where we choose prices in order to learn the distribution of demands & claims for the insurance product, and exploitation, where we myopically choose the best price from the information gathered so far. The performance of the pricing policies is measured in terms of regret: the expected revenue loss caused by not using the optimal price. As is commonplace in insurance, we model demand and claims by GLMs. In our adaptive GLM design, we use the maximum quasi-likelihood estimation (MQLE) to estimate the unknown parameters. We show that, if prices are chosen with suitably decreasing variability, the MQLE parameters eventually exist and converge to the correct values, which in turn implies that the sequence of chosen prices will also converge to the optimal price. In the adaptive GP regression model, we sample demand and claims from Gaussian Processes and then choose selling prices by the upper confidence bound rule. We also analyze these GLM and GP pricing algorithms with delayed claims. Although similar results exist in other domains, this is among the first works to consider dynamic pricing problems in the field of insurance. We also believe this is the first work to consider Gaussian Process regression in the context of insurance pricing. These initial findings suggest that online machine learning algorithms could be a fruitful area of future investigation and application in insurance.

categories: econ.EM math.OC q-fin.MF q-fin.RM stat.ML


title: Magnetization Dynamics in 1D Chains of Ferromagnetic Nanoparticles Coupled with Dipolar Interactions: Blocking Temperature

abstract: There is so far no clear-cut experimental analysis that can determine whether dipole-dipole interactions enhance or reduce the blocking temperature $T_{B}$ of nanoparticle assemblies. It seems that the samples play a central role in the problem and therefore, their geometry should most likely be the key factor in this issue. Yet, in a previous work, J\"onsson and Garcia-Palacios did investigate theoretically this problem in a weak-interaction limit and without the presence of an external DC field. Based on symmetry arguments they reached the conclusion that the variation of the relaxation rate is monotonous. In the presence of an external magnetic field we show that these arguments may no longer hold depending on the experimental geometry. Therefore, the aim of this paper is to evaluate the variation of $T_{B}$ for a model system consisting of a chain of ferromagnetic nanoparticles coupled with long-range dipolar interaction with two different geometries. Rather than addressing a quantitative analysis, we focus on the qualitative variation of $T_{B}$ as a function of the interparticle distance a and of the external field $h$. The two following situations are investigated: a linear chain with a longitudinal axial anisotropy in a longitudinal DC field and a linear chain with a longitudinal axial anisotropy in a transverse field.


title: IDEAS: Immersive Dome Experiences for Accelerating Science

abstract: Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.


title: Visualiza\c{c}\~ao e anima\c{c}\~ao de aut\'omatos em Ocsigen Framework

abstract: Formal Languages and Automata Theory are important foundational topics in Computer Science. Their rigorous and formal characteristics make their learning them demanding. An important support for the assimilation of concepts is the possibility of interactively visualizing concrete examples of these computational models, facilitating understanding them. The tools available are neither complete nor fully support the interactive aspect. This project aims at the development of an interactive web tool in Portuguese to help in an assisted and intuitive way to understand the concepts and algorithms in question, seeing them work step-by-step, through typical examples preloaded or built by the user (an original aspect of our platform). The tool should therefore enable the creation and edition of an automata, as well as execute the relevant classical algorithms such as word acceptance, model conversions, etc. It is also intended to visualize not only the process of construction of the automaton, but also all the steps of applying the given algorithm. This tool uses the Ocsigen Framework because it provides the development of complete and interactive web tools written in OCaml, a functional language with a strong type checking system and therefore perfect for a web page without errors. Ocsigen was also chosen because it allows the creation of dynamic pages with a singular client-server system. This article presents the first phase of the development of the project. It is already possible to create automata, check the nature of its states and verify step-by-step (with undo) the acceptance of a word.


title: Brand vs. Generic: Addressing Non-Adherence, Secular Trends, and Non-Overlap

abstract: While generic drugs offer a cost-effective alternative to brand name drugs, regulators need a method to assess therapeutic equivalence in a post market setting. We develop such a method in the context of assessing the therapeutic equivalence of immediate release (IM) venlafaxine, based on a large insurance claims dataset provided by OptumLabs\textsuperscript{\textregistered}. To properly address this question, our methodology must deal with issues of non-adherence, secular trends in health outcomes, and lack of treatment overlap due to sharp uptake of the generic once it becomes available. We define, identify (under assumptions) and estimate (using G-computation) a causal effect for a time-to-event outcome by extending regression discontinuity to survival curves. We do not find evidence for a lack of therapeutic equivalence of brand and generic IM venlafaxine.


title: Testing independence between two random sets for the analysis of colocalization in bio-imaging

abstract: Colocalization aims at characterizing spatial associations between two fluorescently-tagged biomolecules by quantifying the co-occurrence and correlation between the two channels acquired in fluorescence microscopy. Colocalization is presented either as the degree of overlap between the two channels or the overlays of the red and green images, with areas of yellow indicating colocalization of the molecules. This problem remains an open issue in diffraction-limited microscopy and raises new challenges with the emergence of super-resolution imaging, a microscopic technique awarded by the 2014 Nobel prize in chemistry. We propose GcoPS, for Geo-coPositioning System, an original method that exploits the random sets structure of the tagged molecules to provide an explicit testing procedure. Our simulation study shows that GcoPS unequivocally outperforms the best competitive methods in adverse situations (noise, irregularly shaped fluorescent patterns, different optical resolutions). GcoPS is also much faster, a decisive advantage to face the huge amount of data in super-resolution imaging. We demonstrate the performances of GcoPS on two biological real datasets, obtained by conventional diffraction-limited microscopy technique and by super-resolution technique, respectively.


title: An alternative for the average income estimation using small area methods

abstract: The average household income is one of the most important indexes for decision making and the modelling of economic inequity and poverty. In this work we propose a practical procedure to estimate the average income using small area methods. We illustrate our proposal using information from a multipurpose survey and suitable economic and demographic variables such as the multidimensional poverty and the valorization indexes and the official population projections. We find that the standard relative errors for the income average estimates improve substantially when the proposed methodology is implemented.


title: Provably Efficient Reinforcement Learning with Linear Function Approximation

abstract: Modern Reinforcement Learning (RL) is commonly applied to practical problems with an enormous number of states, where function approximation must be deployed to approximate either the value function or the policy. The introduction of function approximation raises a fundamental set of challenges involving computational and statistical efficiency, especially given the need to manage the exploration/exploitation tradeoff. As a result, a core RL question remains open: how can we design provably efficient RL algorithms that incorporate function approximation? This question persists even in a basic setting with linear dynamics and linear rewards, for which only linear function approximation is needed. This paper presents the first provable RL algorithm with both polynomial runtime and polynomial sample complexity in this linear setting, without requiring a "simulator" or additional assumptions. Concretely, we prove that an optimistic modification of Least-Squares Value Iteration (LSVI)---a classical algorithm frequently studied in the linear setting---achieves $\tilde{\mathcal{O}}(\sqrt{d^3H^3T})$ regret, where $d$ is the ambient dimension of feature space, $H$ is the length of each episode, and $T$ is the total number of steps. Importantly, such regret is independent of the number of states and actions.


title: Non-conventional topological band properties and gapless helical edge states in elastic phononic waveguides with Kekul\'e distortion

abstract: This study investigates the topological behavior of a continuous elastic phononic structure characterized by a 3-term Kekul\'e distortion. The elastic waveguide consists of a hexagonal unit cell whose geometric dimensions are intentionally perturbed according to a generalized Kekul\'e scheme. The resulting structure exhibits an effective Hamiltonian that resembles a quantum spin Hall system, hence suggesting that the waveguide can support helical topological edge states. Using first-principles calculations, we show the existence of nondegenerate pseudospin states and a very peculiar 6-lobe pseudospin texture and Berry curvature pattern. Important insights are also provided concerning the topological states in the Kekul\'e lattice, so far considered indistinguishable, and their critical role in enabling unique gapless edge states, typically not achievable in phononic systems.


title: Reward Advancement: Transforming Policy under Maximum Causal Entropy Principle

abstract: Many real-world human behaviors can be characterized as a sequential decision making processes, such as urban travelers choices of transport modes and routes (Wu et al. 2017). Differing from choices controlled by machines, which in general follows perfect rationality to adopt the policy with the highest reward, studies have revealed that human agents make sub-optimal decisions under bounded rationality (Tao, Rohde, and Corcoran 2014). Such behaviors can be modeled using maximum causal entropy (MCE) principle (Ziebart 2010). In this paper, we define and investigate a general reward trans-formation problem (namely, reward advancement): Recovering the range of additional reward functions that transform the agent's policy from original policy to a predefined target policy under MCE principle. We show that given an MDP and a target policy, there are infinite many additional reward functions that can achieve the desired policy transformation. Moreover, we propose an algorithm to further extract the additional rewards with minimum "cost" to implement the policy transformation.


title: Walking Randomly, Massively, and Efficiently

abstract: We introduce a set of techniques that allow for efficiently generating many independent random walks in the Massive Parallel Computation (MPC) model with space per machine strongly sublinear in the number of vertices. In this space-per-machine regime, many natural approaches to graph problems struggle to overcome the $\Theta(\log n)$ MPC round complexity barrier. Our techniques enable breaking this barrier for PageRank---one of the most important applications of random walks---even in more challenging directed graphs, and for approximate bipartiteness and expansion testing. In the undirected case, we start our random walks from the stationary distribution, which implies that we approximately know the empirical distribution of their next steps. This allows for preparing continuations of random walks in advance and applying a doubling approach. As a result we can generate multiple random walks of length $l$ in $\Theta(\log l)$ rounds on MPC. Moreover, we show that under the popular 1-vs.-2-Cycles conjecture, this round complexity is asymptotically tight. For directed graphs, our approach stems from our treatment of the PageRank Markov chain. We first compute the PageRank for the undirected version of the input graph and then slowly transition towards the directed case, considering convex combinations of the transition matrices in the process. For PageRank, we achieve the following round complexities for damping factor equal to $1 - \epsilon$: * in $O(\log \log n + \log 1 / \epsilon)$ rounds for undirected graphs (with $\tilde O(m / \epsilon^2)$ total space), * in $\tilde O(\log^2 \log n + \log^2 1/\epsilon)$ rounds for directed graphs (with $\tilde O((m+n^{1+o(1)}) / poly\, \epsilon)$ total space).


title: Strain controlled superconductivity in few-layer NbSe2

abstract: The controlled tunability of superconductivity in low-dimensional materials may enable new quantum devices. Particularly in triplet or topological superconductors, tunneling devices such as Josephson junctions etc. can demonstrate exotic functionalities. The tunnel barrier, an insulating or normal material layer separating two superconductors, is a key component for the junctions. Thin layers of NbSe2 have been shown as a superconductor with strong spin orbit coupling, which can give rise to topological superconductivity if driven by a large magnetic exchange field. Here we demonstrate the superconductor-insulator transitions in epitaxially grown few-layer NbSe2 with wafer-scale uniformity on insulating substrates. We provide the electrical transport, Raman spectroscopy, cross-sectional transmission electron microscopy, and X-ray diffraction characterizations of the insulating phase. We show that the superconductor-insulator transition is driven by strain, which also causes characteristic energy shifts of the Raman modes. Our observation paves the way for high quality hetero-junction tunnel barriers to be seamlessly built into epitaxial NbSe2 itself, thereby enabling highly scalable tunneling devices for superconductor-based quantum electronics.


title: On The Evolution Of Operator Complexity Beyond Scrambling

abstract: We study operator complexity on various time scales with emphasis on those much larger than the scrambling period. We use, for systems with a large but finite number of degrees of freedom, the notion of K-complexity employed in arXiv:1812.08657 for infinite systems. We present evidence that K-complexity of ETH operators has indeed the character associated with the bulk time evolution of extremal volumes and actions. Namely, after a period of exponential growth during the scrambling period the K-complexity increases only linearly with time for exponentially long times in terms of the entropy, and it eventually saturates at a constant value also exponential in terms of the entropy. This constant value depends on the Hamiltonian and the operator but not on any extrinsic tolerance parameter. Thus K-complexity deserves to be an entry in the AdS/CFT dictionary. Invoking a concept of K-entropy and some numerical examples we also discuss the extent to which the long period of linear complexity growth entails an efficient randomization of operators.


title: The constructive Kan-Quillen model structure: two new proofs

abstract: We present two new proofs of Simon Henry's result that the category of simplicial sets admits a constructive counterpart of the classical Kan-Quillen model structure. Our proofs are entirely self-contained and avoid complex combinatorial arguments on anodyne extensions. We also give new constructive proofs of the left and right properness of the model structure.


title: Cortical Surface Parcellation using Spherical Convolutional Neural Networks

abstract: We present cortical surface parcellation using spherical deep convolutional neural networks. Traditional multi-atlas cortical surface parcellation requires inter-subject surface registration using geometric features with high processing time on a single subject (2-3 hours). Moreover, even optimal surface registration does not necessarily produce optimal cortical parcellation as parcel boundaries are not fully matched to the geometric features. In this context, a choice of training features is important for accurate cortical parcellation. To utilize the networks efficiently, we propose cortical parcellation-specific input data from an irregular and complicated structure of cortical surfaces. To this end, we align ground-truth cortical parcel boundaries and use their resulting deformation fields to generate new pairs of deformed geometric features and parcellation maps. To extend the capability of the networks, we then smoothly morph cortical geometric features and parcellation maps using the intermediate deformation fields. We validate our method on 427 adult brains for 49 labels. The experimental results show that our method out-performs traditional multi-atlas and naive spherical U-Net approaches, while achieving full cortical parcellation in less than a minute.


title: All-optical quantum sensing of rotational Brownian motion of magnetic molecules

abstract: Sensing local environment through the motional response of small molecules lays the foundation of many fundamental technologies. The information of local viscosity, for example, is contained in the random rotational Brownian motions of molecules. However, detection of the motions is challenging for molecules with sub-nanometer scale or high motional rates. Here we propose and experimentally demonstrate a novel method of detecting fast rotational Brownian motions of small magnetic molecules. With electronic spins as sensors, we are able to detect changes in motional rates, which yield different noise spectra and therefore different relaxation signals of the sensors. As a proof-of-principle demonstration, we experimentally implemented this method to detect the motions of gadolinium (Gd) complex molecules with nitrogen-vacancy (NV) centers in nanodiamonds. With all-optical measurements of the NV centers' longitudinal relaxation, we distinguished binary solutions with varying viscosities. Our method paves a new way for detecting fast motions of sub-nanometer sized magnetic molecules with better spatial resolution than conventional optical methods. It also provides a new tool in designing better contrast agents in magnetic resonance imaging.


title: Origin of gamma-ray families accompanied by halos and detected in experiments with x-ray emulsion chambers

abstract: The phenomenon of gamma-ray families featuring halos that is observed in an experiment with x-ray emulsion chambers (XREC) in the Pamir experiment and in other XREC experiments is explained. The experimental properties of halos are analyzed via a comparison with the results of their simulation. It is shown that gamma-ray families featuring halos are predominantly produced (more than 96 % of them) by protons and heliumnuclei. This makes it possible to employ the experimental properties of halos to estimate the fraction of protons and helium nuclei in the mass composition of primary cosmic radiation.


title: Analytic results for decays of color singlets to $gg$ and $q \bar q$ final states at NNLO QCD with the nested soft-collinear subtraction scheme

abstract: We present compact analytic formulas that describe the decay of colorless particles to both $q \bar q$ and $gg$ final states through next-to-next-to-leading order in perturbative QCD in the context of the nested soft-collinear subtraction scheme. In addition to their relevance for the description of decays like $V \to q \bar q'$, $V= Z,W$, $H \to b \bar b$ and $H \to gg$, these results provide an important building block for calculating NNLO QCD corrections to arbitrary processes at colliders within the nested soft-collinear subtraction scheme.


title: Optimal destabilization of K-unstable Fano varieties via stability thresholds

abstract: We show that for a K-unstable Fano variety, any divisorial valuation computing its stability threshold induces a non-trivial special test configuration preserving the stability threshold. When such a divisorial valuation exists, we show that the Fano variety degenerates to a uniquely determined twisted K-polystable Fano variety. We also show that the stability threshold can be approximated by divisorial valuations induced by special test configurations. As an application of the above results and the analytic work of Datar, Sz\'ekelyhidi, and Ross, we deduce that greatest Ricci lower bounds of Fano manifolds of fixed dimension form a finite set of rational numbers. As a key step in the proofs, we adapt the process of Li and Xu producing special test configurations to twisted K-stability in the sense of Dervan.


title: A modified micromorphic model based on micromechanics for granular materials

abstract: The purpose of this study is to propose a modified micromorphic continuum model for granular materials based on a micromechanics approach. In this model, Cauchy stress and the couple stress are symmetric conjugated with the symmetric strain and the symmetric curvature respectively, and the relative stress measures are asymmetric conjugated with the asymmetric relative strain measures. This modified micromorphic model considers a continuum material point as a granular volume element whose deformation behavior is influenced by the translation and the rotation of particles. And this study proposes that the microscopic actual motion is decomposed into a macroscopic motion and a fluctuation between the macro-micro motion. Based on this decomposition, the micromorphic constitutive relationships are derived for granular materials. In the constitutive relationships, the macroscopic constitutive relationships are first-order because of the introduce of the independent rotation of particle instead of the second-order micro-deformation gradient. Furthermore, the complex constitutive moduli in the micromorphic model are obtained in the expressions of the microstructural information such as the contact stiffness and the internal length.


title: Computational Concentration of Measure: Optimal Bounds, Reductions, and More

abstract: Product measures of dimension $n$ are known to be concentrated in Hamming distance: for any set $S$ in the product space of probability $\epsilon$, a random point in the space, with probability $1-\delta$, has a neighbor in $S$ that is different from the original point in only $O(\sqrt{n\ln(1/(\epsilon\delta))})$ coordinates. We obtain the tight computational version of this result, showing how given a random point and access to an $S$-membership oracle, we can find such a close point in polynomial time. This resolves an open question of [Mahloujifar and Mahmoody, ALT 2019]. As corollaries, we obtain polynomial-time poisoning and (in certain settings) evasion attacks against learning algorithms when the original vulnerabilities have any cryptographically non-negligible probability. We call our algorithm MUCIO ("MUltiplicative Conditional Influence Optimizer") since proceeding through the coordinates, it decides to change each coordinate of the given point based on a multiplicative version of the influence of that coordinate, where influence is computed conditioned on previously updated coordinates. We also define a new notion of algorithmic reduction between computational concentration of measure in different metric probability spaces. As an application, we get computational concentration of measure for high-dimensional Gaussian distributions under the $\ell_1$ metric. We prove several extensions to the results above: (1) Our computational concentration result is also true when the Hamming distance is weighted. (2) We obtain an algorithmic version of concentration around mean, more specifically, McDiarmid's inequality. (3) Our result generalizes to discrete random processes, and this leads to new tampering algorithms for collective coin tossing protocols. (4) We prove exponential lower bounds on the average running time of non-adaptive query algorithms.


title: Scattering amplitudes and contour deformations

abstract: We employ a scalar model to exemplify the use of contour deformations when solving Lorentz-invariant integral equations for scattering amplitudes. In particular, we calculate the onshell 2 -> 2 scattering amplitude for the scalar system. The integrals produce branch cuts in the complex plane of the integrand which prohibit a naive Euclidean integration path. By employing contour deformations, we can also access the kinematical regions associated with the scattering amplitude in Minkowski space. We show that in principle a homogeneous Bethe-Salpeter equation, together with analytic continuation methods such as the Resonances-via-Pad\'e method, is sufficient to determine the resonance pole locations on the second Riemann sheet. However, the scalar model investigated here does not produce resonance poles above threshold but instead virtual states on the real axis of the second sheet, which pose difficulties for analytic continuation methods. To address this, we calculate the scattering amplitude on the second sheet directly using the two-body unitarity relation which follows from the scattering equation.


title: Incrementalizing RASA's Open-Source Natural Language Understanding Pipeline

abstract: As spoken dialogue systems and chatbots are gaining more widespread adoption, commercial and open-sourced services for natural language understanding are emerging. In this paper, we explain how we altered the open-source RASA natural language understanding pipeline to process incrementally (i.e., word-by-word), following the incremental unit framework proposed by Schlangen and Skantze. To do so, we altered existing RASA components to process incrementally, and added an update-incremental intent recognition model as a component to RASA. Our evaluations on the Snips dataset show that our changes allow RASA to function as an effective incremental natural language understanding service.


title: Fibers add Flavor, Part I: Classification of 5d SCFTs, Flavor Symmetries and BPS States

abstract: We propose a graph-based approach to 5d superconformal field theories (SCFTs) based on their realization as M-theory compactifications on singular elliptic Calabi--Yau threefolds. Field-theoretically, these 5d SCFTs descend from 6d $\mathcal{N}=(1,0)$ SCFTs by circle compactification and mass deformations. We derive a description of these theories in terms of graphs, so-called Combined Fiber Diagrams, which encode salient features of the partially resolved Calabi--Yau geometry, and provides a combinatorial way of characterizing all 5d SCFTs that descend from a given 6d theory. Remarkably, these graphs manifestly capture strongly coupled data of the 5d SCFTs, such as the superconformal flavor symmetry, BPS states, and mass deformations. The capabilities of this approach are demonstrated by deriving all rank one and rank two 5d SCFTs. The full potential, however, becomes apparent when applied to theories with higher rank. Starting with the higher rank conformal matter theories in 6d, we are led to the discovery of previously unknown flavor symmetry enhancements and new 5d SCFTs.


title: Simulation of 3D elasto-acoustic wave propagation based on a Discontinuous Galerkin Spectral Element method

abstract: In this paper we present a numerical discretization of the coupled elasto-acoustic wave propagation problem based on a Discontinuous Galerkin Spectral Element (DGSE) approach in a three-dimensional setting. The unknowns of the coupled problem are the displacement field and the velocity potential, in the elastic and the acoustic domains, respectively, thereby resulting in a symmetric formulation. After stating the main theoretical results, we assess the performance of the method by convergence tests carried out on both matching and non-matching grids, and we simulate realistic scenarios where elasto-acoustic coupling occurs. In particular, we consider the case of Scholte waves and the scattering of elastic waves by an underground acoustic cavity. Numerical simulations are carried out by means of the code SPEED, available at http://speed.mox.polimi.it.


title: Using Chinese Characters To Generate Text-Based Passwords For Information Security

abstract: Graphical passwords (GPWs) are in many areas of the current world. Topological graphic passwords (Topsnut-gpws) are a new type of cryptography, and they differ from the existing GPWs. A Topsnut-gpw consists of two parts: one is a topological structure (graph), and one is a set of discrete elements (a graph labelling, or coloring), the topological structure connects these discrete elements together to form an interesting "story". Our idea is to transform Chinese characters into computer and electronic equipments with touch screen by speaking, writing and keyboard for forming Hanzi-graphs and Hanzi-gpws. We will use Hanzigpws to produce text-based passwords (TB-paws). We will introduce flawed graph labellings on disconnected Hanzi-graphs.


title: FARSIDE: A Low Radio Frequency Interferometric Array on the Lunar Farside

abstract: FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood.


title: Timely Cloud Computing: Preemption and Waiting

abstract: The notion of timely status updating is investigated in the context of cloud computing. Measurements of a time-varying process of interest are acquired by a sensor node, and uploaded to a cloud server to undergo some required computations. These computations consume random amounts of service time that are independent and identically distributed across different uploads. After the computations are done, the results are delivered to a monitor, constituting an update. The goal is to keep the monitor continuously fed with fresh updates over time, which is assessed by an age-of-information (AoI) metric. A scheduler is employed to optimize the measurement acquisition times. Following an update, an idle waiting period may be imposed by the scheduler before acquiring a new measurement. The scheduler also has the capability to preempt a measurement in progress if its service time grows above a certain cutoff time, and upload a fresher measurement in its place. Focusing on stationary deterministic policies, in which waiting times are deterministic functions of the instantaneous AoI and the cutoff time is fixed for all uploads, it is shown that the optimal waiting policy that minimizes the long term average AoI has a threshold structure, in which a new measurement is uploaded following an update only if the AoI grows above a certain threshold that is a function of the service time distribution and the cutoff time. The optimal cutoff is then found for standard and shifted exponential service times. While it has been previously reported that waiting before updating can be beneficial for AoI, it is shown in this work that preemption of late updates can be even more beneficial.


title: Change point detection for graphical models in the presence of missing values

abstract: We propose estimation methods for change points in high-dimensional covariance structures with an emphasis on challenging scenarios with missing values. We advocate three imputation like methods and investigate their implications on common losses used for change point detection. We also discuss how model selection methods have to be adapted to the setting of incomplete data. The methods are compared in a simulation study and applied to a time series from an environmental monitoring system. An implementation of our proposals within the R-package hdcd is available via the Supplementary materials.


title: On generalized Macdonald polynomials

abstract: Generalized Macdonald polynomials (GMP) are eigenfunctions of specifically-deformed Ruijsenaars Hamiltonians and are built as triangular polylinear combinations of Macdonald polynomials. They are orthogonal with respect to a modified scalar product, which could be constructed with the help of an increasingly important triangular perturbation theory, showing up in a variety of applications. A peculiar feature of GMP is that denominators in this expansion are fully factorized, which is a consequence of a hidden symmetry resulting from the special choice of the Hamiltonian deformation. We introduce also a simplified but deformed version of GMP, which we call generalized Schur functions. Our basic examples are bilinear in Macdonald polynomials.


title: Proof Theory for Positive Logic with Weak Negation

abstract: Proof-theoretic methods are developed for subsystems of Johansson's logic obtained by extending the positive fragment of intuitionistic logic with weak negations. These methods are exploited to establish properties of the logical systems. In particular, cut-free complete sequent calculi are introduced and used to provide a proof of the fact that the systems satisfy the Craig interpolation property. Alternative versions of the calculi are later obtained by means of an appropriate loop-checking history mechanism. Termination of the new calculi is proved, and used to conclude that the considered logical systems are PSPACE-complete.


title: A mathematical objection to the existence of relativistic mechanical systems of several particles

abstract: We will prove that, in general, a system formed by several particles moving along relativistic trajectories can not be described by a mechanical system. The contradiction that leads to the previous assertion is due to the fact that a mechanical system defines a second order differential equation and this, in turn, induces an absolute time that will generally be incompatible with the proper times of the different particles.


title: Two Neutron Correlation Study in Photofission of Actinides

abstract: It is well known that two fission fragments (FF's) are emitted essentially back to back in the laboratory frame. That can be used widely in many applications as a unique signature of fissionable materials. However, such fission fragments are difficult to detect. The energy and angular distributions of neutrons, on the other hand, are easy to measure, and that distribution will carry information about the fission fragment's energy and angular spectra, as well as the neutron spectra in the fission fragment rest frame. We propose to investigate the two neutron correlation yield resulting from two FF's as a function of different targets, the angle between the two neutrons and the neutron energies. The preliminary calculation of the two neutron correlation shows a huge asymmetry effect: many more neutrons are emitted anti-parallel to each other than parallel to each other. That asymmetry becomes even more if the energy cut on each neutron is done. This study will potentially permit a new technique for actinide detection for homeland security and safeguards applications as well as improve our knowledge of correlated neutron emission.


title: Variational principle for weakly dependent random fields

abstract: Using an alternative notion of entropy introduced by Datta, the max-entropy, we present a new simplified framework to study the minimizers of the specific free energy for random fields which are weakly dependent in the sense of Lewis, Pfister, and Sullivan. The framework is then applied to derive the variational principle for the loop $O(n)$ model and the Ising model in a random percolation environment in the nonmagnetic phase, and we explain how to extend the variational principle to similar models. To demonstrate the generality of the framework, we indicate how to naturally fit into it the variational principle for models with an absolutely summable interaction potential, and for the random-cluster model.


title: Learning to learn with quantum neural networks via classical neural networks

abstract: Quantum Neural Networks (QNNs) are a promising variational learning paradigm with applications to near-term quantum processors, however they still face some significant challenges. One such challenge is finding good parameter initialization heuristics that ensure rapid and consistent convergence to local minima of the parameterized quantum circuit landscape. In this work, we train classical neural networks to assist in the quantum learning process, also know as meta-learning, to rapidly find approximate optima in the parameter landscape for several classes of quantum variational algorithms. Specifically, we train classical recurrent neural networks to find approximately optimal parameters within a small number of queries of the cost function for the Quantum Approximate Optimization Algorithm (QAOA) for MaxCut, QAOA for Sherrington-Kirkpatrick Ising model, and for a Variational Quantum Eigensolver for the Hubbard model. By initializing other optimizers at parameter values suggested by the classical neural network, we demonstrate a significant improvement in the total number of optimization iterations required to reach a given accuracy. We further demonstrate that the optimization strategies learned by the neural network generalize well across a range of problem instance sizes. This opens up the possibility of training on small, classically simulatable problem instances, in order to initialize larger, classically intractably simulatable problem instances on quantum devices, thereby significantly reducing the number of required quantum-classical optimization iterations.


title: Topologically nontrivial Andreev bound states

abstract: Andreev bound states are low energy excitations appearing below the particle-hole gap of superconductors, and are expected to be topologically trivial. Here, we report the theoretical prediction of topologically $nontrivial$ Andreev bound states in one-dimensional superconductors. These states correspond to another topological invariant defined in a synthetic two-dimensional space, the particle-hole Chern number, which we construct in analogy to the spin Chern number in quantum spin Hall systems. Nontrivial Andreev bound states have distinct features and are topologically nonequivalent to Majorana bound states. Yet, they can coexist in the same system, have similar spectral signatures, and materialize with the concomitant opening of the particle-hole gap. The coexistence of Majorana and nontrivial Andreev bound state is the direct consequence of "double dimensionality", i.e., the dimensional embedding of the one-dimensional system in a synthetic two-dimensional space, which allows the definition of two distinct topological invariants ($\mathbb{Z}_2$ and $\mathbb{Z}$) in different dimensionalities.


title: Detecting composite orders in layered models via machine learning

abstract: Determining the phase diagram of systems consisting of smaller subsystems 'connected' via a tunable coupling is a challenging task relevant for a variety of physical settings. A general question is whether new phases, not present in the uncoupled limit, may arise. We use machine learning to study layered spin models, in which the spin variables constituting each of the uncoupled systems (to which we refer as layers) are coupled to each other via an interlayer coupling. In such systems, in general, composite order parameters involving spins of different layers may emerge as a consequence of the interlayer coupling. We focus on the layered Ising and Ashkin-Teller models as a paradigmatic case study, determining their phase diagram via the application of a machine learning algorithm to the Monte Carlo data. Remarkably our technique is able to correctly characterize all the system phases also in the case of hidden order parameters, i.e., order parameters whose expression in terms of the microscopic configurations would require additional preprocessing of the data fed to the algorithm. We correctly retrieve the three known phases of the Ashkin-Teller model with ferromagnetic couplings, including the phase described by a composite order parameter. For the bilayer and trilayer Ising models the phases we find are only the ferromagnetic and the paramagnetic ones. Within the approach we introduce, owing to the construction of convolutional neural networks, naturally suitable for layered image-like data with arbitrary number of layers, no preprocessing of the Monte Carlo data is needed, also with regard to its spatial structure. The physical meaning of our results is discussed and compared with analytical data, where available. Yet, the method can be used without any a priori knowledge of the phases one seeks to find and can be applied to other models and structures.

categories: cond-mat.dis-nn cond-mat.stat-mech physics.comp-ph physics.data-an quant-ph


title: Adversarial Objects Against LiDAR-Based Autonomous Driving Systems

abstract: Deep neural networks (DNNs) are found to be vulnerable against adversarial examples, which are carefully crafted inputs with a small magnitude of perturbation aiming to induce arbitrarily incorrect predictions. Recent studies show that adversarial examples can pose a threat to real-world security-critical applications: a "physical adversarial Stop Sign" can be synthesized such that the autonomous driving cars will misrecognize it as others (e.g., a speed limit sign). However, these image-space adversarial examples cannot easily alter 3D scans of widely equipped LiDAR or radar on autonomous vehicles. In this paper, we reveal the potential vulnerabilities of LiDAR-based autonomous driving detection systems, by proposing an optimization based approach LiDAR-Adv to generate adversarial objects that can evade the LiDAR-based detection system under various conditions. We first show the vulnerabilities using a blackbox evolution-based algorithm, and then explore how much a strong adversary can do, using our gradient-based approach LiDAR-Adv. We test the generated adversarial objects on the Baidu Apollo autonomous driving platform and show that such physical systems are indeed vulnerable to the proposed attacks. We also 3D-print our adversarial objects and perform physical experiments to illustrate that such vulnerability exists in the real world. Please find more visualizations and results on the anonymous website: https://sites.google.com/view/lidar-adv.


title: Classical Limit of Large N Gauge Theories with Conformal Symmetry

abstract: In this paper we study classical limit of conformal field theories realized by large N gauge theories using the generalized coherent states. For generic large N gauge theories with conformal symmetry, we show that the classical limit of them is described by the classical Einstein gravity. This may be regarded as a kind of derivation of the AdS/CFT correspondence.


title: Structure Exploiting Interior Point Methods

abstract: Interior point methods are among the most popular techniques for large scale nonlinear optimization, owing to their intrinsic ability of scaling to arbitrary large problem sizes. Their efficiency has attracted in recent years a lot of attention due to increasing demand for large scale optimization in industry and engineering. A parallel interior point method is discussed that exploits the intrinsic structure of large-scale nonlinear optimization problems so that the solution process can employ massively parallel high-performance computing infastructures. Since the overall performance of interior point methods relies heavily on scalable sparse linear algebra solvers, particular emphasis is given to the underlying algorithms for the distributed solution of the associated sparse linear systems obtained at each iteration from the linearization of the optimality conditions. The interior point algorithm is implemented in a object-oriented parallel IPM solver and applied for the solution of large scale optimal control problems solved in a daily basis for the secure transmission and distribution of electricity in modern power grids.


title: A DNS Study of entrainment in an axisymmetric turbulent jet as an episodic process

abstract: This investigation is based on a DNS of a steady self-preserving incompressible axisymmetric turbulent jet at a Reynolds number of 2400. The DNS data enable accurate maps of the outer irrotational flow field, and also the vorticity field in the turbulent core of the jet. It is found necessary to define two separate boundaries of the jet. The first is an inner boundary (turbulent/nonturbulent, T/NT), from where vorticity rises steeply towards to the core. The second is an outer rotational/irrotational boundary, beyond which the flow may be considered irrotational. The velocity field beyond the outer boundary often has ordered, nearly irrotational circulatory motions. These can be shown, in simpler cases, to be the velocity field induced by one or more vorticity elements in a coherent structure in the turbulent core. A detailed examination of axial and diametral sections indicates that there are periods when there is a large inrush of ambient fluid into parts of the T/NT interface, which gets distorted into a gulf or well that can be both twisted and deep. Sections of these wells often appear as what may be called as lakes of irrotational fluid in diametral sections of the jet flow. Part of the inrushing fluid crosses the T/NT interface within the well and is entrained into the turbulent core, by a process that can legitimately be called nibbling. The duration of such an inrush process can be of the order up to 20 flow units and suggests that entrainment can be an episodic process in which an inrush event accelerates ambient fluid even as it is pushed into a narrowing gulf, where it penetrates the T/NT interface of the gulf by nibbling. In the turbulent round jet, the entrainment burstiness is found to be of order 0.75, comparable to the momentum flux burstiness found in a turbulent boundary layer.


title: In Pursuit of Galactic Archaeology: Astro2020 Science White Paper

abstract: The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.


title: Heisenberg-Kitaev Models on Hyperhoneycomb and Stripyhoneycomb Lattices: 3D-2D Equivalence of Ordered States and Phase Diagrams

abstract: We discuss magnetically ordered states, arising in Heisenberg-Kitaev and related spin models, on three-dimensional (3D) harmonic honeycomb lattices. For large classes of ordered states, we show that they can be mapped onto two-dimensional (2D) counterparts on the honeycomb lattice, with the classical energetics being identical in the 2D and 3D cases. As an example, we determine the phase diagram of the classical nearest-neighbor Heisenberg-Kitaev model on the hyperhoneycomb lattice in a magnetic field: This displays rich and complex behavior akin to its 2D counterpart, with most phases and phase boundaries coinciding exactly. To make contact with the physics of the hyperhoneycomb iridate $\beta$-Li$_2$IrO$_3$, we also include a symmetric off-diagonal $\Gamma$ interaction, discuss its 3D-2D mapping, and determine the relevant phase diagrams. In particular, we demonstrate explicitly the adiabatic equivalence of the spiral magnetic orders in $\alpha$- and $\beta$-Li$_2$IrO$_3$. Our results pave the way to a systematic common understanding of 2D and 3D Kitaev materials.


title: On the prospect of using the maximum circular velocity of halos to encapsulate assembly bias in the galaxy-halo connection

abstract: We investigate a conceptual modification of the halo occupation distribution approach, using the halos' present-day maximal circular velocity, $\vmax$, as an alternative to halo mass. In particular, using a semi-analytic galaxy formation model applied to the Millennium WMAP7 simulation, we explore the extent that switching to $\vmax$ as the primary halo property incorporates the effects of assembly bias into the formalism. We consider fixed number density galaxy samples ranked by stellar mass and examine the variations in the halo occupation functions with either halo concentration or formation time. We find that using $\vmax$ results in a significant reduction in the occupancy variation of the central galaxies, particularly for concentration. The satellites occupancy variation on the other hand increases in all cases. We find effectively no change in the halo clustering dependence on concentration, for fixed bins of $\vmax$ compared to fixed halo mass. Most crucially, we calculate the impact of assembly bias on galaxy clustering by comparing the amplitude of clustering to that of a shuffled galaxy sample, finding that the level of galaxy assembly bias remains largely unchanged. Our results suggest that while using $\vmax$ as a proxy for halo mass diminishes some of occupancy variations exhibited in the galaxy-halo relation, it is not able to encapsulate the effects of assembly bias potentially present in galaxy clustering. The use of other more complex halo properties, such as $\vpeak$, the peak value of $\vmax$ over the assembly history, provides some improvement and warrants further investigation.


title: Cosmological Constraints on Invisible Neutrino Decays Revisited

abstract: Invisible neutrino decay modes are difficult to target at laboratory experiments, and current bounds on such decays from solar neutrino and neutrino oscillation experiments are somewhat weak. It has been known for some time that Cosmology can serve as a powerful probe of invisible neutrino decays. In this work, we show that in order for Big Bang Nucleosynthesis to be successful, the invisible neutrino decay lifetime is bounded to be $\tau_\nu > 10^{-3}\,\text{s}$ at 95\% CL. We revisit Cosmic Microwave Background constraints on invisible neutrino decays, and by using Planck2018 observations we find the following bound on the neutrino lifetime: $\tau_\nu > (1.3-0.3)\times 10^{9}\,\text{s} \, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$ at $95\%$ CL. We show that this bound is robust to modifications of the cosmological model, in particular that it is independent of the presence of dark radiation. We find that lifetimes relevant for Supernova observations ($\tau_\nu \sim 10^{5}\,\text{s}\, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$) are disfavoured at more than $5\,\sigma$ with respect to $\Lambda$CDM given the latest Planck CMB observations. Finally, we show that when including high-$\ell$ Planck polarization data, neutrino lifetimes $\tau_\nu = (2-16)\times 10^{9}\,\text{s} \, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$ are mildly preferred -- with a 1-2 $\sigma$ significance -- over neutrinos being stable.


title: Measuring the Obscuring Column of a Disk Megamaser AGN in a Nearby Merger

abstract: Active Galactic Nuclei (AGNs) hosting disk water megamasers are well known to be obscured by large amounts of gas, likely due to the presence along the line of sight of an almost edge-on disky structure orbiting the supermassive black hole. Correcting for the high obscuration is crucial to infer parameters intrinsic to the source, like its luminosity. We present a broadband X-ray spectral analysis of a water megamaser AGN in an early merger (NGC 5765B), combining Chandra and NuSTAR data. NGC 5765B is highly Compton-thick and reflection-dominated, following the general trend among disk megamasers. Combining the exquisite black hole mass from masers with our X-ray spectroscopy, the Eddington ratio of the megamaser is estimated to be in the $2-14\%$ range, and its robustness is confirmed through SED fitting.


title: Searching the Moon for Extrasolar Material and the Building Blocks of Extraterrestrial Life

abstract: Due to its absence of an atmosphere and relative geological inertness, the Moon's surface records past impacts of objects from the Solar system and beyond. We examine the prospects for discovering extrasolar material near the lunar surface and predict that its abundance is $\mathcal{O}(10)$ parts-per-million (ppm). The abundances of extrasolar organic carbon and biomolecular building blocks (e.g., amino acids) are estimated to be on the order of $0.1$ ppm and $< 0.1$ parts-per-billion (ppb), respectively. We describe strategies for identifying extrasolar material and potentially detecting extrasolar biomolecular building blocks as well as molecular biosignatures of extinct extraterrestrial life. Thus, viewed collectively, we argue that \emph{in situ} lunar exploration can provide vital new clues for astrobiology.


title: $\pi$-Corrected Heisenberg Limit

abstract: We consider the precision $\Delta \varphi$ with which the parameter $\varphi$, appearing in the unitary map $U_\varphi = e^{ i \varphi \Lambda}$ acting on some type of probe system, can be estimated when there is a finite amount of prior information about $\varphi$. We show that, if $U_\varphi$ acts $n$ times in total, then, asymptotically in $n$, there is a tight lower bound $\Delta \varphi \geq \frac{\pi}{n (\lambda_+ - \lambda_-)}$, where $\lambda_+$, $\lambda_-$ are the extreme eigenvalues of the generator $\Lambda$. This is greater by a factor of $\pi$ than the conventional Heisenberg limit, derived from the properties of the quantum Fisher information. That is, the conventional bound is never saturable. Our result makes no assumptions on the measurement protocol, and is relevant not only in the noiseless case but also if noise can be eliminated using quantum error correction techniques.


title: Spherical Contours, IR Divergences and the geometry of Feynman parameter integrands at one loop

abstract: Spherical contours introduced in \cite{SphericalContours} translate the concept of "discontinuity across a branch cut" to Feynman parameter space. In this paper, we further explore spherical contours and connect them to the computation of leading IR divergences of 1 loop graphs directly in Feynman parameter space. These spherical contours can be used to develop a Feynman parameter space analog of "Leading Singularities" of loop integrands which allows us to develop a method of determining Feynman parameter integrands with no reference to the momentum space loop integrand. Finally, we explore some interesting features of Feynman parameter integrands in $\mathcal{N}=4$ SYM.


title: Two-parameter integrable deformations of the $AdS_3 \times S^3 \times T^4$ superstring

abstract: For supercosets with isometry group of the form $\hat{G} \times \hat{G}$, the eta-deformation can be generalised to a two-parameter integrable deformation with independent $q$-deformations of the two copies. We study its kappa-symmetry and write down a formula for the Ramond-Ramond fluxes. We then focus on $\hat{G}= PSU(1,1|2)$ and construct two supergravity backgrounds for the two-parameter integrable deformation of the $AdS_3 \times S^3 \times T^4$ superstring, as well as explore their limits. We also construct backgrounds that are solutions of the weaker generalised supergravity equations of motion and compare them to the literature.


title: Imitation-Projected Programmatic Reinforcement Learning

abstract: We study the problem of programmatic reinforcement learning, in which policies are represented as short programs in a symbolic language. Programmatic policies can be more interpretable, generalizable, and amenable to formal verification than neural policies; however, designing rigorous learning approaches for such policies remains a challenge. Our approach to this challenge -- a meta-algorithm called PROPEL -- is based on three insights. First, we view our learning task as optimization in policy space, modulo the constraint that the desired policy has a programmatic representation, and solve this optimization problem using a form of mirror descent that takes a gradient step into the unconstrained policy space and then projects back onto the constrained space. Second, we view the unconstrained policy space as mixing neural and programmatic representations, which enables employing state-of-the-art deep policy gradient approaches. Third, we cast the projection step as program synthesis via imitation learning, and exploit contemporary combinatorial methods for this task. We present theoretical convergence results for PROPEL and empirically evaluate the approach in three continuous control domains. The experiments show that PROPEL can significantly outperform state-of-the-art approaches for learning programmatic policies.


title: The Kennicutt-Schmidt Law and Gas Scale Height in Luminous and Ultra-Luminous Infrared Galaxies

abstract: A new analysis of high-resolution data from the Atacama Large Millimeter/submillimeter Array (ALMA) for 5 luminous or ultra-luminous infrared galaxies gives a slope for the Kennicutt-Schmidt (KS) relation equal to $1.74^{+0.09}_{\rm -0.07}$ for gas surface densities $\Sigma_{\rm mol}>10^3\;M_\odot$ pc$^{-2}$ and an assumed constant CO-to-H$_2$ conversion factor. The velocity dispersion of the CO line, $\sigma_v$, scales approximately as the inverse square root of $\Sigma_{\rm mol}$, making the empirical gas scale height determined from $H\sim0.5\sigma^2/(\pi G\Sigma_{\rm mol})$ nearly constant, 150-190 pc, over 1.5 orders of magnitude in $\Sigma_{\rm mol}$. This constancy of $H$ implies that the average midplane density, which is presumably dominated by CO-emitting gas for these extreme star-forming galaxies, scales linearly with the gas surface density, which, in turn, implies that the gas dynamical rate (the inverse of the free-fall time) varies with $\Sigma_{\rm mol}^{1/2}$, thereby explaining most of the super-linear slope in the KS relation. Consistent with these relations, we also find that the mean efficiency of star formation per free-fall time is roughly constant, 5%-7%, and the gas depletion time decreases at high $\Sigma_{\rm mol}$, reaching only $\sim 16$ Myr at $\Sigma_{\rm mol}\sim10^4\;M_\odot$ pc$^{-2}$. The variation of $\sigma_v$ with $\Sigma_{\rm mol}$ and the constancy of $H$ are in tension with some feedback-driven models, which predict $\sigma_v$ to be more constant and $H$ to be more variable. However, these results are consistent with simulations in which large-scale gravity drives turbulence through a feedback process that maintains an approximately constant Toomre $Q$ instability parameter.


title: Optical conductivity in graphene: hydrodynamic regime

abstract: Experimental investigation of hydrodynamics in electron fluids is a highly topical research area that emerged during the last few years. A recent measurement of the optical conductivity in graphene [P. Gallagher et.al, Science 364, 158 (2019)] offers a possibility of experimental determination of microscopic time scales describing scattering processes in the electronic fluid. In this paper, I report a theoretical calculation of the optical conductivity in graphene at arbitrary doping levels, within the whole "hydrodynamic" temperature range, and for arbitrary, non-quantizing magnetic fields. The obtained results are in good agreement with the available experimental data.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el physics.flu-dyn physics.plasm-ph


title: Polarizability of the nucleon

abstract: The status of the experimental and theoretical investigations on the polarizabilities of the nucleon is presented. This includes a confirmation of the validity of the previously introduced recommended values of the polarizabilities [1,2]. It is shown that the only meaningful approach to a prediction of the polarizabilities is obtained from the nonsubtracted dispersion theory, where the appropriate degrees of freedom taken from other precise experimental data are taken in account. The present values of the recommended polarizabilities are $\alpha_p= 12.0 \pm 0.5$, $\beta_p= 1.9 \mp 0.5$, $\alpha_n= 12.6 \pm 1.2$, $\beta_n= 2.6 \mp 1.2$ in units of $10^{-4}$fm$^3$ and $\gamma^{(p)}_\pi= -36.4 \pm 1.5$, $\gamma^{(n)}_\pi =+58.6 \pm 4.0$, $\gamma^{(p)}_0=-0.58 \pm 0.20$, $\gamma^{(n)}_0= +0.38\pm 0.22$ in units of $10^{-4}$fm$^4$.


title: Nonlinear Perturbations of a periodic magnetic Choquard equation with Hardy-Littlewood-Sobolev critical exponent

abstract: In this paper, we consider the following magnetic nonlinear Choquard equation \[-(\nabla+iA(x))^2u+ V(x)u = \left(\frac{1}{|x|^{\alpha}}*|u|^{2_{\alpha}^*}\right) |u|^{2_{\alpha}^*-2} u + \lambda f(u)\ \textrm{ in }\ \R^N,\] where $2_{\alpha}^{*}=\frac{2N-\alpha}{N-2}$ is the critical exponent in the sense of the Hardy-Littlewood-Sobolev inequality, $\lambda>0$, $N\geq 3$, $0<\alpha< N$, $A: \mathbb{R}^{N}\rightarrow \mathbb{R}^{N}$ is an $C^1$, $\mathbb{Z}^N$-periodic vector potential and $V$ is a continuous scalar potential given as a perturbation of a periodic potential. Under suitable assumptions on different types of nonlinearities $f$, namely, $f(x,u)=\left(\frac{1}{|x|^{\alpha}}*|u|^{p}\right)|u|^{p-2} u$ for $(2N-\alpha)/N<p<2^{*}_{\alpha}$, then $f(u)=|u|^{p-1} u$ for $1<p<2^*-1$ and $f(u)=|u|^{2^* - 2}u$ (where $2^*=2N/(N-2)$), we prove the existence of at least one ground state solution for this equation by variational methods if $p$ belongs to some intervals depending on $N$ and $\lambda$.


title: Two-dimensional Dirac operators with singular interactions supported on closed curves

abstract: We study the two-dimensional Dirac operator with an arbitrary combination of electrostatic and Lorentz scalar $\delta$-interactions of constant strengths supported on a smooth closed curve. For any combination of the coupling constants a rigorous description of the self-adjoint realizations of the operators is given and the qualitative spectral properties are described. The analysis covers also all so-called critical combinations of coupling constants, for which there is a loss of regularity in the operator domain. In this case, if the mass is non-zero, the resulting operator has an additional point in the essential spectrum, and the position of this point inside the central gap can be made arbitrary by a suitable choice of the coupling constants. The analysis is based on a combination of the extension theory of symmetric operators with a detailed study of boundary integral operators viewed as periodic pseudodifferential operators.


title: Circumbinary disk inner radius as a diagnostic for disk--binary misalignment

abstract: We investigate the misalignment of the circumbinary disk around the binary HD 98800 BaBb with eccentricity $e\simeq 0.8$. \cite{Kennedy2019} observed the disk to be either at an inclination of $48^{\circ}$ or polar aligned to the binary orbital plane. Their simulations showed that alignment from $48^\circ$ to a polar configuration can take place on a shorter timescale than the age of this system. We perform hydrodynamical numerical simulations in order to estimate the cavity size carved by the eccentric binary for different disk inclinations as an independent check of polar alignment. Resonance theory suggests that torques on the inner parts of a polar disk around such a highly eccentric binary are much weaker than in the coplanar case, indicating a significantly smaller central cavity than in the coplanar case. We show that the inferred inner radius (from carbon monoxide measurements) of the accretion disk around BaBb can exclude the possibility of it being mildly inclined with respect to the binary orbital plane and therefore confirm the polar configuration. This study constitutes an important diagnostic for misaligned circumbinary disks, since it potentially allows us to infer the disk inclination from observed gas disk inner radii.


title: Subleading soft dressings of asymptotic states in QED and perturbative quantum gravity

abstract: We construct Faddeev-Kulish states in QED and perturbative quantum gravity to subleading order in the soft momentum expansion and to first order in the coupling constant, using the charge conservation formula of asymptotic symmetries associated with the tree-level subleading soft theorems. We demonstrate that the emission and absorption of soft photons/gravitons in dressed amplitudes vanish. The fact that no additional soft radiation may be added to a dressed amplitude supports the claim that, in the dressed state formalism, the soft and hard sectors of scattering processes are correlated. We also show that the dressed virtual amplitudes are equivalent to the infrared-finite part of the traditional amplitudes constructed using Fock states. Since there is no real soft radiation in the asymptotic Hilbert space, the dressed state formalism gives the same cross sections as the Bloch-Nordsieck method.


title: Commutant lifting and Nevanlinna-Pick interpolation in several variables

abstract: This paper concerns a commutant lifting theorem and a Nevanlinna-Pick type interpolation result in the setting of multipliers from vector-valued Drury-Arveson space to a large class of vector-valued reproducing kernel Hilbert spaces over the unit ball in $\mathbb{C}^n$. The special case of reproducing kernel Hilbert spaces includes all natural examples of Hilbert spaces like Hardy space, Bergman space and weighted Bergman spaces over the unit ball.


title: Quantum discord-breaking and discord-annihilating channels

abstract: Quantum discord-breaking channels were previously defined as the local channels that act on subsystem $A$ to produce classical-quantum states across system $AB$. However, unlike entanglement, discord is asymmetric. Here, we characterise the discord-breaking channels that act on subsystem $B$, thus completing the overall description of discord-breaking channels. We then introduce the notion of discord-annihilating channels, which act globally on system $AB$ to destroy quantum discord, and find their closed form. Discord-annihilating channels have clear operational description, involving subspace projections on subsystem $A$ and conditional preparation of fixed states on subsystem $B$.


title: Understanding Leakage Currents through $Al_2O_3$ on $SrTiO_3$

abstract: Leakage currents through insulators received continuous attention for decades, owing to their importance for a wide range of technologies, and interest in their fundamental mechanisms. This work investigates the leakage currents through atomic layer deposited (ALD) $Al_2O_3$, grown on $SrTiO_3$. This combination is not only a key building block of oxide electronics, but also a clean system for studying the leakage mechanisms without interfacial layers that form on most of the conventional bottom electrodes. We show how tiny differences in the deposition process can have a dramatic effect on the leakage behavior. Detailed analysis of the leakage behavior rules out Fowler-Nordheim tunneling (FNT) and thermionic emission, and leaves the trap-related mechanisms of trap-assisted tunneling (TAT) and Poole-Frenkel as the likely mechanisms. After annealing the sample in air, the currents are reduced, which is ascribed to transition from trap-based mechanism to FNT, due to the elimination of the traps. The dramatic role of the assumptions regarding the flat-band voltage used for analysis is critically discussed, and the sensitivity of the extracted parameters on this magnitude is quantitatively described. We show that field effect devices based on structures similar to those described here, should be able to modulate $>10^{13} cm^{-2}$ electrons. These results provide general guidelines for reducing and analyzing leakage currents in insulators, and highlight some of the possible approaches and pitfalls in their analysis.


title: Predicting engagement in online social networks: Challenges and opportunities

abstract: Since the introduction of social media, user participation or engagement has received little research attention. In this survey article, we establish the notion of participation in social media and main challenges that researchers may face while exploring this phenomenon. We surveyed a handful of research articles that had been done in this area, and tried to extract, analyze and summarize the techniques performed by the researchers. We classified these works based on our task definitions, and explored the machine learning models that have been used for any kind of participation prediction. We also explored the vast amount of features that have been proven useful, and classified them into categories for better understanding and ease of re-implementation. We have found that the success of a technique mostly depends on the type of the network that has been researched on, and there is no universal machine learning algorithm or feature sets that works reasonably well in all types of social media. There is a lack of attempts in implementing state-of-the-art machine learning techniques like neural networks, and the possibility of transfer learning and domain adaptation has not been explored.