title: Coarse-graining Molecular Systems by Spectral Matching

abstract: Coarse-graining has become an area of tremendous importance within many different research fields. For molecular simulation, coarse-graining bears the promise of finding simplified models such that long-time simulations of large-scale systems become computationally tractable. While significant progress has been made in tuning thermodynamic properties of reduced models, it remains a key challenge to ensure that relevant kinetic properties are retained by coarse-grained dynamical systems. In this study, we focus on data-driven methods to preserve the rare-event kinetics of the original system, and make use of their close connection to the low-lying spectrum of the system's generator. Building on work by Crommelin and Vanden-Eijnden, SIAM Multiscale Model. Simul. (2011), we present a general framework, called spectral matching, which directly targets the generator's leading eigenvalue equations when learning parameters for coarse-grained models. We discuss different parametric models for effective dynamics and derive the resulting data-based regression problems. We show that spectral matching can be used to learn effective potentials which retain the slow dynamics, but also to correct the dynamics induced by existing techniques, such as force matching.

categories: physics.chem-ph math.DS math.PR physics.bio-ph


title: Curvature computations in Finsler Geometry using a distinguished class of anisotropic connections

abstract: We show how to compute tensor derivatives and curvature tensors using affine connections. This allows for all computations to be obtained without using coordinate systems, in a way that parallels the computations appearing in classical Riemannian Geometry. In particular, we obtain Bianchi identities for the curvature tensor of any anisotropic connection, we compare the curvature tensors of any two anisotropic connections, and we find a family of anisotropic connections which are well suited to study the geometry of Finsler metrics.


title: Mathematical analysis of weak and strong solutions to an evolutionary model for magnetoviscoelasticity

abstract: The paper is concerned with the analysis of an evolutionary model for magnetoviscoelastic materials in two dimensions. The model consists of a Navier-Stokes system featuring a dependence of the stress tensor on elastic and magnetic terms, a regularized system for the evolution of the deformation gradient and the Landau-Lifshitz-Gilbert system for the dynamics of the magnetization. First, we show that our model possesses global in time weak solutions, thus extending work by Bene\v{s}ov\'a et al. 2018. Compared to that work, we include the stray field energy and relax the assumptions on the elastic energy density. Second, we prove the local in time existence of strong solutions. Both existence results are based on the Galerkin method. Finally, we show a weak-strong uniqueness property.


title: Synthetic Neural Vision System Design for Motion Pattern Recognition in Dynamic Robot Scenes

abstract: Insects have tiny brains but complicated visual systems for motion perception. A handful of insect visual neurons have been computationally modeled and successfully applied for robotics. How different neurons collaborate on motion perception, is an open question to date. In this paper, we propose a novel embedded vision system in autonomous micro-robots, to recognize motion patterns in dynamic robot scenes. Here, the basic motion patterns are categorized into movements of looming (proximity), recession, translation, and other irrelevant ones. The presented system is a synthetic neural network, which comprises two complementary sub-systems with four spiking neurons -- the lobula giant movement detectors (LGMD1 and LGMD2) in locusts for sensing looming and recession, and the direction selective neurons (DSN-R and DSN-L) in flies for translational motion extraction. Images are transformed to spikes via spatiotemporal computations towards a switch function and decision making mechanisms, in order to invoke proper robot behaviors amongst collision avoidance, tracking and wandering, in dynamic robot scenes. Our robot experiments demonstrated two main contributions: (1) This neural vision system is effective to recognize the basic motion patterns corresponding to timely and proper robot behaviors in dynamic scenes. (2) The arena tests with multi-robots demonstrated the effectiveness in recognizing more abundant motion features for collision detection, which is a great improvement compared with former studies.


title: Spin-to-charge conversion in magnetic Weyl semimetals

abstract: Weyl semimetals (WSMs) are a newly discovered class of quantum materials which can host a number of exotic bulk transport properties, such as the chiral magnetic effect, negative magneto-resistance, and the anomalous Hall effect. In this work, we investigate theoretically the spin-to-charge conversion in a bilayer consisting of a magnetic WSM and a normal metal (NM), where a charge current can be induced in the WSM by an spin current injection at the interface. We show that the induced charge current exhibits a peculiar anisotropy: it vanishes along the magnetization orientation of the magnetic WSM, regardless of the direction of the injected spin. This anisotropy originates from the unique band structure of magnetic WSMs and distinguishes the spin-to-charge conversion effect in WSM/NM structures from that observed in other systems, such as heterostructures involving heavy metals or topological insulators. The induced charge current depends strongly on injected spin orientation, as well as on the position of the Fermi level relative to the Weyl nodes and the separation between them. These dependencies provide additional means to control and manipulate spin-charge conversion in these topological materials.


title: Momentum-space probability density of ${}^6$He in Halo Effective Field Theory

abstract: We compute the momentum-space probability density of ${}^6$He at leading order in Halo EFT. In this framework, the ${}^6$He nucleus is treated as a three-body problem with a ${}^4$He core ($c$) and two valence neutrons ($n$). This requires the $nn$ and $nc$ t-matrices as well as a $cnn$ force as input in the Faddeev equations. Since the $nc$ t-matrix corresponds to an energy-dependent potential, we consider the consequent modifications to the standard normalization and orthogonality conditions. We find that these are small for momenta within the domain of validity of Halo EFT. In this regime, the ${}^6$He probability density is regulator independent, provided the cutoff is significantly above the EFT breakdown scale.


title: Some cyclic properties of $L_1$-graphs

abstract: A graph $G$ is called an $L_1$-graph if $d(u)+d(v)\ge|N(u)\cup N(v)\cup N(w)|-1$ for every triple of vertices $u,v,w$ where $u$ and $v$ are at distance 2 and $w\in N(u)\cap N(v)$. Asratian et al. (1996) proved that all finite connected $L_1$-graphs on at least three vertices such that $|N(u)\cap N(v)|\ge2$ for each pair of vertices $u,v$ at distance 2 are Hamiltonian, except for a simple family $\mathcal{K}$ of exceptions. We show that not all such graphs are pancyclic, but that any non-Hamiltonian cycle in such a graph can be extended to a larger cycle containing all vertices of the original cycle and at most two other vertices. We also prove a similar result for paths whose endpoints do not have any common neighbors.


title: A monotone scheme for G-equations with application to the explicit convergence rate of robust central limit theorem

abstract: We propose a monotone approximation scheme for a class of fully nonlinear PDEs called G-equations. Such equations arise often in the characterization of G-distributed random variables in a sublinear expectation space. The proposed scheme is constructed recursively based on a piecewise constant approximation of the viscosity solution to the G-equation. We establish the convergence of the scheme and determine the convergence rate with an explicit error bound, using the comparison principles for both the scheme and the equation together with a mollification procedure. The first application is obtaining the convergence rate of Peng's robust central limit theorem with an explicit bound of Berry-Esseen type. The second application is an approximation scheme with its convergence rate for the Black-Scholes-Barenblatt equation.


title: Excited states in bilayer graphene quantum dots

abstract: We report on ground- and excited state transport through an electrostatically defined few-hole quantum dot in bilayer graphene in both parallel and perpendicular applied magnetic fields. A remarkably clear level scheme for the two-particle spectra is found by analyzing finite bias spectroscopy data within a two-particle model including spin and valley degrees of freedom. We identify the two-hole ground-state to be a spin-triplet and valley-singlet state. This spin alignment can be seen as Hund's rule for a valley-degenerate system, which is fundamentally different to quantum dots in carbon nano tubes and GaAs-based quantum dots. The spin-singlet excited states are found to be valley-triplet states by tilting the magnetic field with respect to the sample plane. We quantify the exchange energy to be 0.35meV and measure a valley and spin g-factor of 36 and 2, respectively.


title: Blow-up at space infinity for solutions of a system of non-autonomous semilinear heat equations

abstract: In this paper we will see that the global or local existence of solutions to \begin{eqnarray*} \dfrac{\partial u_{1}}{\partial t} & = & \mathit{k}_{1} (t) \Delta u_{1} + h_{1}(t) u_{1}^{p_{11}} u_{2}^{p_{12}},\\ \dfrac{\partial u_{2}}{\partial t} & = & \mathit{k}_{2} (t) \Delta u_{2} + h_{2}(t) u_{2}^{p_{22}} u_{1}^{p_{21}}, \end{eqnarray*} depends on the initial datums and the global or local existence of solutions to \begin{eqnarray*} \dfrac{dy_{1}}{dt} & = & h_{1}(t) y_{1}^{p_{11}}(t) y_{2}^{p_{12}}(t),\\ \dfrac{dy_{2}}{dt} & = & h_{2}(t) y_{2}^{p_{22}}(t) y_{1}^{p_{21}}(t). \end{eqnarray*} We also give some bounds for the maximal existence time of the partial differential system. Moreover, if such existence time is finite and $\min\{p_{11} + p_{12},p_{22}+p_{21}\} > 1$ then we will prove the partial differential system has solutions that blows-up at space infinite.


title: The Numerical Stability of Regularized Barycentric Interpolation Formulae for Interpolation and Extrapolation

abstract: The $\ell_2-$ and $\ell_1-$regularized modified Lagrange interpolation formulae over $[-1,1]$ are deduced in this paper. This paper mainly analyzes the numerical characteristics of regularized barycentric interpolation formulae, which are presented in [C. An and H.-N. Wu, 2019], and regularized modified Lagrange interpolation formulae for both interpolation and extrapolation. Regularized barycentric interpolation formulae can be carried out in $\mathcal{O}(N)$ operations based on existed algorithms [H. Wang, D. Huybrechs and S. Vandewalle, Math. Comp., 2014], and regularized modified Lagrange interpolation formulae can be realized in an algorithm of $\mathcal{O}(N\log N)$ operations. For interpolation, the regularized modified Lagrange interpolation formulae are blessed with backward stability and forward stability, whereas the regularized barycentric interpolation formulae are only provided with forward stability. For extrapolation, the regularized barycentric interpolation formulae meet loss of accuracy outside $[-1,1]$, but the regularized modified Lagrange interpolation formulae still work. Consistent results for extrapolation are also verified outside the Chebfun ellipse (a special Bernstein ellipse) in the complex plain. Finally, we illustrate that regularized interpolation formulae perform better than classical interpolation formulae without regularization in noise reduction.


title: Cyberbullying and Traditional Bullying in Greece: An Empirical Study

abstract: The rapid evolution of technology is a preferred method of interacting. A new world was created for young people who are sending emails, visiting websites, using webcams and chat rooms, and instant messaging through social media. In the past years, people used face-to-face communication though, in recent years, people are using internet technology in order to communicate with each other. As a consequence the communication change, created a new type of bullying, cyberbullying, in which bullying is taking place by using internet technology. Cyberbullying is a phenomenon which is increasing day by day over the world. This was one of the reasons that prompted us to do this survey. Furthermore, bullies aggressive reaction is affected by four factors which are conductive during childhood.


title: Reinforcement Learning with Probabilistic Guarantees for Autonomous Driving

abstract: Designing reliable decision strategies for autonomous urban driving is challenging. Reinforcement learning (RL) has been used to automatically derive suitable behavior in uncertain environments, but it does not provide any guarantee on the performance of the resulting policy. We propose a generic approach to enforce probabilistic guarantees on an RL agent. An exploration strategy is derived prior to training that constrains the agent to choose among actions that satisfy a desired probabilistic specification expressed with linear temporal logic (LTL). Reducing the search space to policies satisfying the LTL formula helps training and simplifies reward design. This paper outlines a case study of an intersection scenario involving multiple traffic participants. The resulting policy outperforms a rule-based heuristic approach in terms of efficiency while exhibiting strong guarantees on safety.


title: Explicit Spatial Encoding for Deep Local Descriptors

abstract: We propose a kernelized deep local-patch descriptor based on efficient match kernels of neural network activations. Response of each receptive field is encoded together with its spatial location using explicit feature maps. Two location parametrizations, Cartesian and polar, are used to provide robustness to a different types of canonical patch misalignment. Additionally, we analyze how the conventional architecture, i.e. a fully connected layer attached after the convolutional part, encodes responses in a spatially variant way. In contrary, explicit spatial encoding is used in our descriptor, whose potential applications are not limited to local-patches. We evaluate the descriptor on standard benchmarks. Both versions, encoding 32x32 or 64x64 patches, consistently outperform all other methods on all benchmarks. The number of parameters of the model is independent of the input patch resolution.


title: Radial oscillations of boson stars made of ultralight repulsive dark matter

abstract: We compute the lowest frequency radial oscillation modes of boson stars. It is assumed that the object is made of pseudo-Goldstone bosons subjected to a scalar potential that leads to a repulsive self-interaction force, and which is characterized by two unknown mass scales $m$ (mass of the particle) and $F$ (decay constant). First we integrate the Tolman-Oppenheimer-Volkoff equations for the hydrostatic equilibrium of the star, and then we solve the Sturm-Liouville boundary value problem for the perturbations using the shooting method. The effective potential that enters into the Schr{\"o}dinger-like equation as well as several associated eigenfunctions are shown as well. Moreover, we found that the large frequency separation, i.e. the difference between consecutive modes, is proportional to the square root of the mass of the star and the cube of the mass scale defined by $\Lambda \equiv \sqrt{m F}$.


title: Comparison of statistical post-processing methods for probabilistic NWP forecasts of solar radiation

abstract: The increased usage of solar energy places additional importance on forecasts of solar radiation. Solar panel power production is primarily driven by the amount of solar radiation and it is therefore important to have accurate forecasts of solar radiation. Accurate forecasts that also give information on the forecast uncertainties can help users of solar energy to make better solar radiation based decisions related to the stability of the electrical grid. To achieve this, we apply statistical post-processing techniques that determine relationships between observations of global radiation (made within the KNMI network of automatic weather stations in the Netherlands) and forecasts of various meteorological variables from the numerical weather prediction (NWP) model HARMONIE-AROME (HA) and the atmospheric composition model CAMS. Those relationships are used to produce probabilistic forecasts of global radiation. We compare 7 different statistical post-processing methods, consisting of two parametric and five non-parametric methods. We find that all methods are able to generate probabilistic forecasts that improve the raw global radiation forecast from HA according to the root mean squared error (on the median) and the potential economic value. Additionally, we show how important the predictors are in the different regression methods. We also compare the regression methods using various probabilistic scoring metrics, namely the continuous ranked probability skill score, the Brier skill score and reliability diagrams. We find that quantile regression and generalized random forests generally perform best. In (near) clear sky conditions the non-parametric methods have more skill than the parametric ones.


title: ACRONYM III: Radial Velocities for 336 Candidate Young Low-Mass Stars in the Solar Neighborhood, Including 78 Newly Confirmed Young Moving Group Members

abstract: Young, low-mass stars in the solar neighborhood are vital for completing the mass function for nearby, young coeval groups, establishing a more complete census for evolutionary studies, and providing targets for direct-imaging exoplanet and/or disk studies. We present properties derived from high-resolution optical spectra for 336 candidate young nearby, low-mass stars. These include measurements of radial velocities and age diagnostics such as H$\alpha$ and Li $\lambda$6707 equivalent widths. Combining our radial velocities with astrometry from Gaia DR2, we provide full 3D kinematics for the entire sample. We combine the measured spectroscopic youth information with additional age diagnostics (e.g., X-ray and UV fluxes, CMD positions) and kinematics to evaluate potential membership in nearby, young moving groups and associations. We identify 78 objects in our sample as bonafide members of 10 different moving groups, 15 of which are completely new members or have had their group membership reassigned. We also reject 44 previously proposed candidate moving group members. Furthermore, we have newly identified or confirmed the youth of numerous additional stars that do not belong to any currently known group, and find 69 co-moving systems using Gaia DR2 astrometry. We also find evidence that the Carina association is younger than previously thought, with an age similar to the $\beta$ Pictoris moving group ($\sim$22 Myr).


title: Strong Stability Preserving Integrating Factor Two-step Runge--Kutta Methods

abstract: Problems that feature significantly different time scales, where the stiff time-step restriction comes from a linear component, implicit-explicit (IMEX) methods alleviate this restriction if the concern is linear stability. However, where the SSP property is needed, IMEX SSP Runge-Kutta (SSP-IMEX) methods have very restrictive time-steps. An alternative to SSP-IMEX schemes is to adopt an integrating factor approach to handle the linear component exactly and step the transformed problem forward using some time-evolution method. The strong stability properties of integrating factor Runge--Kutta methods were previously established, where it was shown that it is possible to define explicit integrating factor Runge-Kutta methods that preserve strong stability properties satisfied by each of the two components when coupled with forward Euler time-stepping. It was proved that the solution will be SSP if the transformed problem is stepped forward with an explicit SSP Runge-Kutta method that has non-decreasing abscissas. However, explicit SSP Runge-Kutta methods have an order barrier of p=4, and sometimes higher order is desired. In this work we consider explicit SSP two-step Runge--Kutta integrating factor methods to raise the order. We show that strong stability is ensured if the two-step Runge-Kutta method used to evolve the transformed problem is SSP and has non-decreasing abscissas. We find such methods up to eighth order and present their SSP coefficients. Adding a step allows us to break the fourth order barrier on explicit SSP Runge-Kutta methods; furthermore, our explicit SSP two-step Runge--Kutta methods with non-decreasing abscissas typically have larger SSP coefficients than the corresponding one-step methods.


title: Acoustic modes of pulsating axion stars: Non-radial oscillations

abstract: We compute the lowest frequency non-radial oscillation modes of dilute axion stars. The effective potential that enters into the Schr{\"o}dinger-like equation, several associated eigenfunctions, and the large as well as the small frequency separations are shown as well.


title: On the equality problem of generalized Bajraktarevi\'c means

abstract: The purpose of this paper is to investigate the equality problem of generalized Bajraktarevi\'c means, i.e., to solve the functional equation \begin{equation}\label{E0}\tag{*} f^{(-1)}\bigg(\frac{p_1(x_1)f(x_1)+\dots+p_n(x_n)f(x_n)}{p_1(x_1)+\dots+p_n(x_n)}\bigg)=g^{(-1)}\bigg(\frac{q_1(x_1)g(x_1)+\dots+q_n(x_n)g(x_n)}{q_1(x_1)+\dots+q_n(x_n)}\bigg), \end{equation} which holds for all $x=(x_1,\dots,x_n)\in I^n$, where $n\geq 2$, $I$ is a nonempty open real interval, the unknown functions $f,g:I\to\mathbb{R}$ are strictly monotone, $f^{(-1)}$ and $g^{(-1)}$ denote their generalized left inverses, respectively, and $p=(p_1,\dots,p_n):I\to\mathbb{R}_{+}^n$ and $q=(q_1,\dots,q_n):I\to\mathbb{R}_{+}^n$ are also unknown functions. This equality problem in the symmetric two-variable (i.e., when $n=2$) case was already investigated and solved under sixth-order regularity assumptions by Losonczi in 1999. In the nonsymmetric two-variable case, assuming three times differentiability of $f$, $g$ and the existence of $i\in\{1,2\}$ such that either $p_i$ is twice continuously differentiable and $p_{3-i}$ is continuous on $I$, or $p_i$ is twice differentiable and $p_{3-i}$ is once differentiable on $I$, we prove that \eqref{E0} holds if and only if there exist four constants $a,b,c,d\in\mathbb{R}$ with $ad\neq bc$ such that \begin{equation*} cf+d>0,\qquad g=\frac{af+b}{cf+d},\qquad\mbox{and}\qquad q_\ell=(cf+d)p_\ell\qquad (\ell\in\{1,\dots,n\}). \end{equation*} In the case $n\geq 3$, we obtain the same conclusion with weaker regularity assumptions. Namely, we suppose that $f$ and $g$ are three times differentiable, $p$ is continuous and there exist $i,j,k\in\{1,\dots,n\}$ with $i\neq j\neq k\neq i$ such that $p_i,p_j,p_k$ are differentiable.


title: Identification of Parameters for Large-scale Kinetic Models

abstract: Inverse problem for the identification of the parameters for large-scale systems of nonlinear ordinary differential equations (ODEs) arising in systems biology is analyzed. In a recent paper in \textit{Mathematical Biosciences, 305(2018), 133-145}, the authors implemented the numerical method suggested by one of the authors in \textit{J. Optim. Theory Appl., 85, 3(1995), 509-526} for identification of parameters in moderate scale models of systems biology. This method combines Pontryagin optimization or Bellman's quasilinearization with sensitivity analysis and Tikhonov regularization. We suggest modification of the method by embedding a method of staggered corrector for sensitivity analysis and by enhancing multi-objective optimization which enables application of the method to large-scale models with practically non-identifiable parameters based on multiple data sets, possibly with partial and noisy measurements. We apply the modified method to a benchmark model of a three-step pathway modeled by 8 nonlinear ODEs with 36 unknown parameters and two control input parameters. The numerical results demonstrate geometric convergence with a minimum of five data sets and with minimum measurements per data set. Software package \textit{qlopt} is developed and posted in GitHub. MATLAB package AMIGO2 is used to demonstrate advantage of \textit{qlopt} over most popular methods/software such as \textit{lsqnonlin}, \textit{fmincon} and \textit{nl2sol}.


title: Total positivity in reductive groups, II

abstract: We continue the study of total positivity in reductive groups following my 1994 paper.


title: Exact Rate-Distortion in Autoencoders via Echo Noise

abstract: Compression is at the heart of effective representation learning. However, lossy compression is typically achieved through simple parametric models like Gaussian noise to preserve analytic tractability, and the limitations this imposes on learning are largely unexplored. Further, the Gaussian prior assumptions in models such as variational autoencoders (VAEs) provide only an upper bound on the compression rate in general. We introduce a new noise channel, \emph{Echo noise}, that admits a simple, exact expression for mutual information for arbitrary input distributions. The noise is constructed in a data-driven fashion that does not require restrictive distributional assumptions. With its complex encoding mechanism and exact rate regularization, Echo leads to improved bounds on log-likelihood and dominates $\beta$-VAEs across the achievable range of rate-distortion trade-offs. Further, we show that Echo noise can outperform flow-based methods without the need to train additional distributional transformations.


title: A Discussion on Solving Partial Differential Equations using Neural Networks

abstract: Can neural networks learn to solve partial differential equations (PDEs)? We investigate this question for two (systems of) PDEs, namely, the Poisson equation and the steady Navier--Stokes equations. The contributions of this paper are five-fold. (1) Numerical experiments show that small neural networks (< 500 learnable parameters) are able to accurately learn complex solutions for systems of partial differential equations. (2) It investigates the influence of random weight initialization on the quality of the neural network approximate solution and demonstrates how one can take advantage of this non-determinism using ensemble learning. (3) It investigates the suitability of the loss function used in this work. (4) It studies the benefits and drawbacks of solving (systems of) PDEs with neural networks compared to classical numerical methods. (5) It proposes an exhaustive list of possible directions of future work.


title: Mid-infrared (3-8 {\mu}m) $\text{Ge}_{1-y}\text{Sn}_y$ alloys (0.15 < $y$ < 0.30): synthesis, structural, and optical properties

abstract: $\text{Ge}_{1-y}\text{Sn}_y$ alloys with compositions in the 0.15 < $y$ < 0.30 range have been grown directly on Si substrates using a chemical vapor deposition approach that allows for growth temperatures as high as 290 $^{\circ}$C. The films show structural properties that are consistent with results from earlier materials with much lower Sn concentrations. These include the lattice parameter and the Ge-Ge Raman frequency, which are found to depend linearly on composition. The simplicity of the structures, directly grown on Si, makes it possible to carry out detailed optical studies. Sharp absorption edges are found, reaching 8 $\mu$m near $y$ =0.3. The compositional dependence of edge energies shows a cubic deviation from the standard quadratic alloy expression. The cubic term may dramatically impact the ability of the alloys to cover the long-wavelength (8-12 $\mu$m) mid-IR atmospheric window.


title: Neutrino oscillations and decoherence in short-GRB progenitors

abstract: Neutrinos are produced in cosmic accelerators, like active galactic nuclei (AGNs), blazars, supernova (SN) remnants and gamma-ray bursts (GRBs). On their way to the Earth they experience flavor-oscillations. The interactions of the neutrinos coming from the source with other particles, e.g. intergalactic primordial neutrinos or heavy-mass right-handed neutrinos, in their way to the detector may transform the original wave packet in pointer states. This phenomenon, known as decoherence, becomes important in the reconstruction of processes at the source. In this work we study neutrino emission in short GRBs by adopting the Fireshell Model. We consider $e^{-}e^{+}$-pair annihilation as the main channel for neutrino production. We compare the properties of the neutrino-flux with the characteristic photon-signal produced once the transparency condition is reached. We study the effects of flavor-oscillations and decoherence as neutrinos travel from the region near the black-hole (BH) event-horizon outwards. We consider the source to be in thermal equilibrium, and calculate energy distribution functions for electrons and neutrinos. To compute the effects of decoherence we use a Gaussian model. In this scenario the emitted electron-neutrinos transform into pointer states consisting of $67.8\%$ electron-neutrinos and $32.2\%$ as a combination of mu and tau neutrinos. We found that decoherence plays an important role in the evolution of the neutrino wave packet, leading to the detected pointer states on Earth.


title: Anisotropic upper critical field of pristine and proton-irradiated single crystals of the magnetically ordered superconductor RbEuFe$_4$As$_4$

abstract: We present a study of the upper critical field, H$_{c2}$, of pristine and proton-irradiated RbEuFe$_4$As$_4$ crystals in pulsed magnetic fields of up to 65 T. The data for H$_{c2}$ reveal pronounced downwards curvature, particularly for the in-plane field orientation, and a superconducting anisotropy that decreases with decreasing temperature. These features are indicative of Pauli paramagnetic limiting. For the interpretation of these data, we use a model of a clean single-band superconductor with an open Fermi surface in the shape of a warped cylinder, which includes strong paramagnetic limiting. Fits to the data reveal that the in-plane upper critical field is Pauli paramagnetic limited, while the out-of-plane upper critical field is orbitally limited and that the orbital and paramagnetic fields have opposite anisotropies. A consequence of this particular combination is the unusual inversion of the anisotropy, $H_{c2}^{ab} < H_{c2}^c$, of the irradiated sample at temperatures below 10 K. The fits also yield an in-plane Maki parameter, $\alpha_M^{110} \approx$ 2.6, exceeding the critical value for the formation of the Fulde-Ferrell-Larkin-Ovchinnikov state. Nevertheless, the current measurements did not reveal direct evidence for the occurrence of this state.


title: Tutorial: Safe and Reliable Machine Learning

abstract: This document serves as a brief overview of the "Safe and Reliable Machine Learning" tutorial given at the 2019 ACM Conference on Fairness, Accountability, and Transparency (FAT* 2019). The talk slides can be found here: https://bit.ly/2Gfsukp, while a video of the talk is available here: https://youtu.be/FGLOCkC4KmE, and a complete list of references for the tutorial here: https://bit.ly/2GdLPme.


title: Dynamics of racial segregation and gentrification in New York City

abstract: Racial residential segregation is interconnected with several other phenomena such as income inequalities, property values inequalities, and racial disparities in health and in education. Furthermore, recent literature suggests the phenomena of gentrification as a cause of perpetuation or increase of racial residential segregation in some American cities. In this paper, we analyze the dynamics of racial residential segregation for white, black, Asian, and Hispanic citizens in New York City in the years of 1990, 2000, and 2010. It was possible to observe that segregation between white and Hispanic citizens, and discrimination between white and Asian ones has grown, while segregation between white and black is quite stable. Furthermore, we analyzed the per capita income and the Gini coefficient in each segregated zone, showing that the highest inequalities occur in the zones where there is overlap of high-density zones of pair of races. Focusing on census tracts that have changed density of population during these twenty years, and, particularly, by analyzing white and black people's segregation, our analysis reveals that a positive flux of white (black) people is associated to a substantial increase (decrease) of the property values, as compared with the city mean. Furthermore, by clustering the region of high density of black citizens, we measured the variation of area and displacement of the four biggest clusters in the period from 1990 to 2010. The large displacements (~ 1.6 Km) observed for two of these clusters, namely, one in the neighborhood of Harlem and the other inside the borough of Brooklyn, led to the emergence of typically gentrified regions.


title: An LGMD Based Competitive Collision Avoidance Strategy for UAV

abstract: Building a reliable and efficient collision avoidance system for unmanned aerial vehicles (UAVs) is still a challenging problem. This research takes inspiration from locusts, which can fly in dense swarms for hundreds of miles without collision. In the locust's brain, a visual pathway of LGMD-DCMD (lobula giant movement detector and descending contra-lateral motion detector) has been identified as collision perception system guiding fast collision avoidance for locusts, which is ideal for designing artificial vision systems. However, there is very few works investigating its potential in real-world UAV applications. In this paper, we present an LGMD based competitive collision avoidance method for UAV indoor navigation. Compared to previous works, we divided the UAV's field of view into four subfields each handled by an LGMD neuron. Therefore, four individual competitive LGMDs (C-LGMD) compete for guiding the directional collision avoidance of UAV. With more degrees of freedom compared to ground robots and vehicles, the UAV can escape from collision along four cardinal directions (e.g. the object approaching from the left-side triggers a rightward shifting of the UAV). Our proposed method has been validated by both simulations and real-time quadcopter arena experiments.


title: Extreme wave statistics in co-propagating windsea and swell

abstract: We investigate how the extreme wave statistics of a windsea is modified by a following swell, by means of laboratory experiments and simulations using a High Order Spectral Method (HOSM) of long-crested sea. The windsea spectrum is kept equal for all cases, while the swell is altered. Analysis of the combined wave system gives the impression that the mixed sea has milder extreme wave statistics than the windsea alone, especially when the nonlinearities in the two systems clearly differ. Upon partitioning the mixed sea into windsea and swell, the windsea part is found to be nearly unaffected by the swell and to be governed by essentially the same statistics irrespective of a swell being present or not. This result is found for skewness, kurtosis and exceedance probability of envelope and crest height.


title: More testing or more disease? A counterfactual approach to explaining observed increases in positive tests over time

abstract: Observed gonorrhea case rates (number of positive tests per 100,000 individuals) increased by 75 percent in the United States between 2009 and 2017, predominantly among men. However, testing recommendations by the Centers for Disease Control and Prevention (CDC) have also changed over this period with more frequent screening for sexually transmitted infections (STIs) recommended among men who have sex with men (MSM) who are sexually active. In this and similar disease surveillance settings, a common question is whether observed increases in the overall proportion of positive tests over time is due only to increased testing of diseased individuals, increased underlying disease or both. By placing this problem within a counterfactual framework, we can carefully consider untestable assumptions under which this question may be answered and, in turn, a principled approach to statistical analysis. This report outlines this thought process.


title: Attention-Passing Models for Robust and Data-Efficient End-to-End Speech Translation

abstract: Speech translation has traditionally been approached through cascaded models consisting of a speech recognizer trained on a corpus of transcribed speech, and a machine translation system trained on parallel texts. Several recent works have shown the feasibility of collapsing the cascade into a single, direct model that can be trained in an end-to-end fashion on a corpus of translated speech. However, experiments are inconclusive on whether the cascade or the direct model is stronger, and have only been conducted under the unrealistic assumption that both are trained on equal amounts of data, ignoring other available speech recognition and machine translation corpora. In this paper, we demonstrate that direct speech translation models require more data to perform well than cascaded models, and while they allow including auxiliary data through multi-task training, they are poor at exploiting such data, putting them at a severe disadvantage. As a remedy, we propose the use of end-to-end trainable models with two attention mechanisms, the first establishing source speech to source text alignments, the second modeling source to target text alignment. We show that such models naturally decompose into multi-task-trainable recognition and translation tasks and propose an attention-passing technique that alleviates error propagation issues in a previous formulation of a model with two attention stages. Our proposed model outperforms all examined baselines and is able to exploit auxiliary training data much more effectively than direct attentional models.


title: Essential self-adjointness of perturbed quadharmonic operators on Riemannian manifolds with an application to the separation problem

abstract: We consider perturbed quadharmonic operators, $\Delta^4 + V$, acting on sections of a Hermitian vector bundle over a complete Riemannian manifold, with the potential $V$ satisfying a bound from below by a non-positive function depending on the distance from a point. Under a bounded geometry assumption on the Hermitian vector bundle and the underlying Riemannian manifold, we give a sufficient condition for the essential self-adjointness of such operators. We then apply this to prove the separation property in $L^2$ when the perturbed operator acts on functions.


title: On the Phases of a Complex Matrix

abstract: In this paper, we define the phases of a complex sectorial matrix to be its canonical angles, which are uniquely determined from a sectorial decomposition of the matrix. Various properties of matrix phases are studied, including those of compressions, Schur complements, matrix products, and sums. In particular, by exploiting a notion known as the compound numerical range, we establish a majorization relation between the phases of the eigenvalues of $AB$ and the phases of $A$ and $B$. This is then applied to investigate the rank of $I + AB$ with phase information of $A$ and $B$, which plays a crucial role in feedback stability analysis. A pair of problems: banded sectorial matrix completion and decomposition is studied. The phases of the Kronecker and Hadamard products are also discussed.


title: Time-delay between neutrinos and gamma-rays in short GRBs

abstract: In this work we make an estimate of the time-delay between signals, recorded at detectors on Earth, of neutrinos and photons originated in a short gamma-ray burst. We describe the geometry and dynamics of the system according to the Fireshell model. The delay in the photon's arrival time is produced because the system is originally opaque to radiation; thus, the photons remain trapped and thermalize until the transparency condition is reached. We calculate the time interval between neutrino- and photon- emission in the black hole frame and transform it to the observer-frame using Lorentz transformations. We obtain a difference in the arrival time at Earth of $\Delta t_{Earth} \approx 854.57$ s.


title: Comparison of two early warning systems for regional flash flood hazard forecasting

abstract: The anticipation of flash flood events is crucial to issue warnings to mitigate their impact. This work presents a comparison of two early warning systems for real-time flash flood hazard forecasting at regional scale. The two systems are based in a gridded drainage network and they use weather radar precipitation inputs to assess the hazard level in different points of the study area, considering the return period (in years) as the indicator of the flash flood hazard. The essential difference between the systems is that one is a rainfall-based system (ERICHA), using the upstream basin-aggregated rainfall as the variable to determine the hazard level, while the other (Flood-PROOFS) is a system based on a distributed rainfall-runoff model to compute the streamflows at pixel scale. The comparison has been done for three rainfall events in the autumn of 2014 that resulted in severe flooding in the Liguria region (Northwest of Italy). The results obtained by the two systems show many similarities, particularly for larger catchments and for large return periods (extreme floods).


title: New binary black hole mergers in the second observing run of Advanced LIGO and Advanced Virgo

abstract: We report the detection of new binary black hole merger events in the publicly available data from the second observing run of advanced LIGO and advanced Virgo (O2). The mergers were discovered using the new search pipeline described in Venumadhav et al. [Phys. Rev. D 100, 023011 (2019)], and are above the detection thresholds as defined in Abbott et al. (LIGO Scientific and Virgo Collaborations) [Phys. Rev. X 9, 031040 (2019).] Three of the mergers (GW170121, GW170304, GW170727) have inferred probabilities of being of astrophysical origin $p_{\rm astro} > 0.98$. The remaining three (GW170425, GW170202, GW170403) are less certain, with $p_{\rm astro}$ ranging from 0.5 to 0.8. The newly found mergers largely share the statistical properties of previously reported events, with the exception of GW170403, the least secure event, which has a highly negative effective spin parameter $\chi_{\rm eff}$ . The most secure new event, GW170121 ($p_{\rm astro} > 0.99$), is also notable due to its inferred negative value of $\chi_{\rm eff}$, which is inconsistent with being positive at the ~95.8% confidence level. The new mergers nearly double the sample of gravitational wave events reported from O2, and present a substantial opportunity to explore the statistics of the binary black hole population in the Universe. The number of detected events is not surprising since we estimate that the detection volume of our pipeline may be larger than that of other pipelines by as much as a factor of two (with significant uncertainties in the estimate). The increase in volume is larger when the constituent detectors of the network have very different sensitivities, as is likely to be the case in current and future runs.

categories: astro-ph.HE astro-ph.CO astro-ph.IM gr-qc


title: The surprisingly low carbon mass in the debris disk around HD 32297

abstract: Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks.


title: A Linear Upper Bound on the Weisfeiler-Leman Dimension of Graphs of Bounded Genus

abstract: The Weisfeiler-Leman (WL) dimension of a graph is a measure for the inherent descriptive complexity of the graph. While originally derived from a combinatorial graph isomorphism test called the Weisfeiler-Leman algorithm, the WL dimension can also be characterised in terms of the number of variables that is required to describe the graph up to isomorphism in first-order logic with counting quantifiers. It is known that the WL dimension is upper-bounded for all graphs that exclude some fixed graph as a minor (Grohe, JACM 2012). However, the bounds that can be derived from this general result are astronomic. Only recently, it was proved that the WL dimension of planar graphs is at most 3 (Kiefer, Ponomarenko, and Schweitzer, LICS 2017). In this paper, we prove that the WL dimension of graphs embeddable in a surface of Euler genus $g$ is at most $4g+3$. For the WL dimension of graphs embeddable in an orientable surface of Euler genus $g$, our approach yields an upper bound of $2g+3$.


title: Resurgence, a problem of missing exponential corrections in asymptotic expansions

abstract: It is well known that perturbative expansions of path integrals are divergent. These expansions are to be understood as asymptotic expansions, which encode the limiting behaviour of the path integral for positive small coupling. Conventionally, the method of Borel summation assigns a finite answer to the divergent expansion. Still, the Borel sum might not encode the full information of a function, because it misses exponentially small corrections. In the present work, we consider a slight variation of the conventional Borel summation, in which a generalised Borel transform (an inverse Laplace transform) is followed by a directional Laplace transform. These new tools will allow us to give perhaps better answers to typical problems in Borel summation: missing exponential corrections and ambiguities in the Borel summation. In addition, we will define resurgence as a connection between the discontinuity of a function and the coefficients of its asymptotic expansion. From this definition, we will be able to reduce resurgence to the problem of missing exponential corrections in asymptotic expansions and understand, within a unified framework, different approaches to resurgence found in the literature.


title: From Hall's Marriage Theorem to Boolean Satisfiability and Back

abstract: Motivated by the application of Hall's Marriage Theorem in various LP-rounding problems, we introduce a generalization of the classical marriage problem (CMP) that we call the Fractional Marriage Problem. We show that the Fractional Marriage Problem is NP-Complete by reduction from Boolean Satisfiability (SAT). We show that when we view the classical marriage problem (a.k.a. bipartite matching) as a sub-class of SAT we get a new class of polynomial-time satisfiable SAT instances that we call CMP-SAT, different from the classically known polynomial-time satisfiable SAT instances 2-SAT, Horn-SAT and XOR-SAT. We next turn to the problem of recognizing CMP-SAT instances, first using SAT embeddings, and then using their embeddings within the universe of Fractional Marriage Problems (FMPs). In the process we are led to another generalization of the CMP that we call the Symmetric Marriage Problem, which is polynomial time decidable and leads to a slight enlargement of the CMP-SAT class. We develop a framework for simplifying FMP problems to identify CMP instances that we call Fragment Logic. Finally we give a result that sheds light on how expressive the FMP need be to still be NP-Complete. The result gives a second NP-Complete reduction of the FMP, this time to Tripartite Matching. We conclude with a wide assortment of suggested additional problems.


title: The Tur\'an number of blow-ups of trees

abstract: A conjecture of Erd\H{o}s from 1967 asserts that any graph on $n$ vertices which does not contain a fixed $r$-degenerate bipartite graph $F$ has at most $Cn^{2-1/r}$ edges, where $C$ is a constant depending only on $F$. We show that this bound holds for a large family of $r$-degenerate bipartite graphs, including all $r$-degenerate blow-ups of trees. Our results generalise many previously proven cases of the Erd\H{o}s conjecture, including the related results of F\"uredi and Alon, Krivelevich and Sudakov. Our proof uses supersaturation and a random walk on an auxiliary graph.


title: Learning Discriminative Model Prediction for Tracking

abstract: The current strive towards end-to-end trainable computer vision systems imposes major challenges for the task of visual tracking. In contrast to most other vision problems, tracking requires the learning of a robust target-specific appearance model online, during the inference stage. To be end-to-end trainable, the online learning of the target model thus needs to be embedded in the tracking architecture itself. Due to the imposed challenges, the popular Siamese paradigm simply predicts a target feature template, while ignoring the background appearance information during inference. Consequently, the predicted model possesses limited target-background discriminability. We develop an end-to-end tracking architecture, capable of fully exploiting both target and background appearance information for target model prediction. Our architecture is derived from a discriminative learning loss by designing a dedicated optimization process that is capable of predicting a powerful model in only a few iterations. Furthermore, our approach is able to learn key aspects of the discriminative loss itself. The proposed tracker sets a new state-of-the-art on 6 tracking benchmarks, achieving an EAO score of 0.440 on VOT2018, while running at over 40 FPS. The code and models are available at https://github.com/visionml/pytracking.


title: The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

abstract: We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.


title: Asymmetric balance in symmetry breaking

abstract: Spontaneous symmetry breaking is central to our understanding of physics and explains many natural phenomena, from cosmic scales to subatomic particles. Its use for applications requires devices with a high level of symmetry, but engineered systems are always imperfect. Surprisingly, the impact of such imperfections has barely been studied, and restricted to a single asymmetry. Here, we experimentally study spontaneous symmetry breaking with two controllable asymmetries. We remarkably find that features typical of spontaneous symmetry breaking, while destroyed by one asymmetry, can be restored by introducing a second asymmetry. In essence, asymmetries are found to balance each other. Our study illustrates aspects of the universal unfolding of the pitchfork bifurcation, and provides new insights into a key fundamental process. It also has practical implications, showing that asymmetry can be exploited as an additional degree of freedom. In particular, it would enable sensors based on symmetry breaking or exceptional points to reach divergent sensitivity even in presence of imperfections. Our experimental implementation built around an optical fiber ring additionally constitutes the first observation of the polarization symmetry breaking of passive driven nonlinear resonators.


title: Joint Discriminative and Generative Learning for Person Re-identification

abstract: Person re-identification (re-id) remains challenging due to significant intra-class variations across different cameras. Recently, there has been a growing interest in using generative models to augment training data and enhance the invariance to input changes. The generative pipelines in existing methods, however, stay relatively separate from the discriminative re-id learning stages. Accordingly, re-id models are often trained in a straightforward manner on the generated data. In this paper, we seek to improve learned re-id embeddings by better leveraging the generated data. To this end, we propose a joint learning framework that couples re-id learning and data generation end-to-end. Our model involves a generative module that separately encodes each person into an appearance code and a structure code, and a discriminative module that shares the appearance encoder with the generative module. By switching the appearance or structure codes, the generative module is able to generate high-quality cross-id composed images, which are online fed back to the appearance encoder and used to improve the discriminative module. The proposed joint learning framework renders significant improvement over the baseline without using generated data, leading to the state-of-the-art performance on several benchmark datasets.


title: The 2019 Discovery of a Meteor of Interstellar Origin

abstract: The earliest confirmed interstellar object, `Oumuamua, was discovered in the Solar System by Pan-STARRS in 2017, allowing for a calibration of the abundance of interstellar objects of its size $\sim 100\;$ m. This was followed by the discovery of Borisov, which allowed for a similar calibration of its size $\sim 0.4 - 1 \mathrm{\; km}$. One would expect a much higher abundance of significantly smaller interstellar objects, with some of them colliding with Earth frequently enough to be noticeable. Based on the CNEOS catalog of bolide events, we identify the $\sim 0.45$m meteor detected at 2014-01-08 17:05:34 UTC as originating from an unbound hyperbolic orbit with 99.999\% confidence. The U.S. Department of Defense has since verified that "the velocity estimate reported to NASA is sufficiently accurate to indicate an interstellar trajectory." We infer that the meteor had an asymptotic speed of $v_{\infty} \sim 42.1 \pm 5.5\; \mathrm{km \; s^{-1}}$ outside of the solar system. Its origin is approximately towards R.A. $49.4 \pm 4.1^{\circ}$ and declination $11.2 \pm 1.8^{\circ}$, implying that its initial velocity vector was $58\pm6\; \mathrm{km\;s^{-1}}$ away from the velocity of the Local Standard of Rest (LSR). Its high LSR speed implies a possible origin from the deep interior of a planetary system or a star in the thick disk of the Milky Way galaxy. The local number density of its population is $10^{6{^{+0.75}_{-1.5}}} \; \mathrm{AU^{-3}}$ or $9 \times 10^{21{^{+0.75}_{-1.5}}} \; \mathrm{pc^{-3}}$ (necessitating 0.2 -- 20 Earth masses of material to be ejected per local star). We show that the detections of CNEOS 2014-01-08, `Oumuamua, and Borisov collectively imply that the differential size distribution in good agreement with a collisional distribution, with a power-law slope is $q \sim 3.6 \pm 0.5$, where the quoted uncertainty corresponds to $2 \sigma$.


title: Autonomous adaptive noise characterization in quantum computers

abstract: New quantum computing architectures consider integrating qubits as sensors to provide actionable information useful for decoherence mitigation on neighboring data qubits, but little work has addressed how such schemes may be efficiently implemented in order to maximize information utilization. Techniques from classical estimation and dynamic control, suitably adapted to the strictures of quantum measurement, provide an opportunity to extract augmented hardware performance through automation of low-level characterization and control. In this work, we present an autonomous learning framework, Noise Mapping for Quantum Architectures (NMQA), for adaptive scheduling of sensor-qubit measurements and efficient spatial noise mapping (prior to actuation) across device architectures. Via a two-layer particle filter, NMQA receives binary measurements and determines regions within the architecture that share common noise processes; an adaptive controller then schedules future measurements to reduce map uncertainty. Numerical analysis and experiments on an array of trapped ytterbium ions demonstrate that NMQA outperforms brute-force mapping by up-to $18$x ($3$x) in simulations (experiments), calculated as a reduction in the number of measurements required to map a spatially inhomogeneous magnetic field with a target error metric. As an early adaptation of robotic control to quantum devices, this work opens up exciting new avenues in quantum computer science.


title: From multi-dimensional black scholes to Hamilton jacobi

abstract: The first widely used financial model is linked to dynamical Hamilton jacobi model


title: First passage time for Slepian process with linear barrier

abstract: In this paper we extend results of L.A. Shepp by finding explicit formulas for the first passage probability $F_{a,b}(T\, |\, x)={\rm Pr}(S(t)<a+bt \text{ for all } t\in[0,T]\,\, | \,\,S(0)=x)$, for all $T>0$, where $S(t)$ is a Gaussian process with mean 0 and covariance $\mathbb{E} S(t)S(t')=\max\{0,1-|t-t'|\}\,.$ We then extend the results to the case of piecewise-linear barriers and outline applications to change-point detection problems. Previously, explicit formulas for $F_{a,b}(T\, |\, x)$ were known only for the cases $b=0$ (constant barrier) or $T\leq 1$ (short interval).


title: Multidimensional topological Galois theory

abstract: In this preprint we present an outline of the multidimensional version of topological Galois theory. The theory studies topological obstruction to solvability of equations "in finite terms" (i.e. to their solvability by radicals, by elementary functions, by quadratures and so on). This preprint is based on the author's book on topological Galois theory. It contains definitions, statements of results and comments to them. Basically no proofs are presented. This preprint was written as a part of the comments to a new edition (in preparation) of the classical book "Integration in finite terms" by J.F.~Ritt.


title: Quantum knots and knotted zeros

abstract: In this paper we show how to place Michael Berry's discovery of knotted zeros in the quantum states of hydrogen in the context of general knot theory and in the context of our formulations for quantum knots. Berry gave a time independent wave function for hydrogen, as a map from three space to the complex plane and such that the inverse image of zero in the complex plane contains a knotted curve in three space. We show that for knots in three space this is a generic situation in that every smooth knot K in three space has a smooth classifying map f from three space to the complex plane such that the inverse image of zero is the knot K. This leaves open the question of characterizing just when such f are wave-functions for quantum systems. One can compare this result with the work of Mark Dennis and his collaborators, with the work of Daniel Peralta-Salas and his collaborators, and with the work of Lee Rudolph. Our approach provides great generality to the structure of knotted zeros of a wavefunction and opens up many new avenues for research in the relationships of quantum theory and knot theory. We show how this classifying construction can be related our previous work on two dimensional and three dimensional mosaic and lattice quantum knots.


title: Diamond Twin

abstract: As noticed in 2006 by the author of the present article, the hypothetical crystal---described by crystallographer F. Laves (1932) for the first time and designated ``Laves' graph of girth ten" by geometer H. S. M. Coxeter (1955)---is a unique crystal net sharing a remarkable symmetric property with the diamond crystal, thus deserving to be called the diamond twin although their shapes look quite a bit different at first sight. In this short note, we shall provide an interesting mutual relationship between them, expressed in terms of ``building blocks" and ``period lattices." This may give further justification to employ the word ``twin." What is more, our discussion brings us to the notion of ``orthogonally symmetric lattice," a generalization of irreducible root lattices, which makes the diamond and its twin very distinct among all crystal structures.


title: Critical elements for connectivity analysis of brain networks

abstract: In recent years, new and important perspectives were introduced in the field of neuroimaging with the emergence of the connectionist approach. In this new context, it is important to know not only which brain areas are activated by a particular stimulus but, mainly, how these areas are structurally and functionally connected, distributed, and organized in relation to other areas. Additionally, the arrangement of the network elements, i.e., its topology, and the dynamics they give rise to are also important. This new approach is called connectomics. It brings together a series of techniques and methodologies capable of systematizing, from the different types of signals and images of the nervous system, how neuronal units to brain areas are connected. Through this approach, the different patterns of connectivity can be graphically and mathematically represented by the so-called connectomes. The connectome uses quantitative metrics to evaluate structural and functional information from images of neural tracts and pathways or signals from the metabolic and/or electrophysiologic activity of cell populations or brain areas. Besides, with adequate treatment of this information, it is also possible to infer causal relationships. In this way, structural and functional evaluations are complementary descriptions which, together, represent the anatomic and physiologic neural properties, establishing a new paradigm for understanding how the brain functions by looking at brain connections. Here, we highlight five critical elements of a network that allows an integrative analysis, focusing mainly on a functional description. These elements include; (i) the properties of its nodes; (ii) the metrics for connectivity and coupling between nodes; (iii) the network topologies; (iv) the network dynamics and (v) the interconnections between different domains and scales of network representations.


title: Ultra-compact spherically symmetric dark matter charged star objects

abstract: We study the properties of ultra-compact spherically symmetric dark matter sector star objects, being the solution of Einstein equations with two $U(1)$-gauge fields. One of them is the ordinary Maxwell field, while the auxiliary gauge field pertains to the hidden sector, and mimics the properties of dark matter. The visible and hidden sectors are coupled by a kinetic mixing term with a coupling constant $\alpha$. We also investigate the possibility of condensation of charged scalar field around the reflecting dark matter star object. It happens that dark matter sector, both dark matter charge and coupling constant, cause shrinking of the star radius for which the condensation may occur.


title: Estimation of Linear Motion in Dense Crowd Videos using Langevin Model

abstract: Crowd gatherings at social and cultural events are increasing in leaps and bounds with the increase in population. Surveillance through computer vision and expert decision making systems can help to understand the crowd phenomena at large gatherings. Understanding crowd phenomena can be helpful in early identification of unwanted incidents and their prevention. Motion flow is one of the important crowd phenomena that can be instrumental in describing the crowd behavior. Flows can be useful in understanding instabilities in the crowd. However, extracting motion flows is a challenging task due to randomness in crowd movement and limitations of the sensing device. Moreover, low-level features such as optical flow can be misleading if the randomness is high. In this paper, we propose a new model based on Langevin equation to analyze the linear dominant flows in videos of densely crowded scenarios. We assume a force model with three components, namely external force, confinement/drift force, and disturbance force. These forces are found to be sufficient to describe the linear or near-linear motion in dense crowd videos. The method is significantly faster as compared to existing popular crowd segmentation methods. The evaluation of the proposed model has been carried out on publicly available datasets as well as using our dataset. It has been observed that the proposed method is able to estimate and segment the linear flows in the dense crowd with better accuracy as compared to state-of-the-art techniques with substantial decrease in the computational overhead.


title: Bounded and Approximate Strong Satisfiability in Workflows

abstract: There has been a considerable amount of interest in recent years in the problem of workflow satisfiability, which asks whether the existence of constraints in a workflow specification makes it impossible to allocate authorized users to each step in the workflow. Recent developments have seen the workflow satisfiability problem (WSP) studied in the context of workflow specifications in which the set of steps may vary from one instance of the workflow to another. This, in turn, means that some constraints may only apply to certain workflow instances. Inevitably, WSP becomes more complex for such workflow specifications. Other approaches have considered the possibility of associating costs with the violation of `soft' constraints and authorizations. Workflow satisfiability in this context becomes a question of minimizing the cost of allocating users to steps in the workflow. In this paper, we introduce new problems, which we believe to be of practical relevance, that combine these approaches. In particular, we consider the question of whether, given a workflow specification with costs and a `budget', all possible workflow instances have an allocation of users to steps that does not exceed the budget. We design a fixed-parameter tractable algorithm to solve this problem parameterized by the total number of steps, release points and xor branchings.


title: Focus Is All You Need: Loss Functions For Event-based Vision

abstract: Event cameras are novel vision sensors that output pixel-level brightness changes ("events") instead of traditional video frames. These asynchronous sensors offer several advantages over traditional cameras, such as, high temporal resolution, very high dynamic range, and no motion blur. To unlock the potential of such sensors, motion compensation methods have been recently proposed. We present a collection and taxonomy of twenty two objective functions to analyze event alignment in motion compensation approaches (Fig. 1). We call them Focus Loss Functions since they have strong connections with functions used in traditional shape-from-focus applications. The proposed loss functions allow bringing mature computer vision tools to the realm of event cameras. We compare the accuracy and runtime performance of all loss functions on a publicly available dataset, and conclude that the variance, the gradient and the Laplacian magnitudes are among the best loss functions. The applicability of the loss functions is shown on multiple tasks: rotational motion, depth and optical flow estimation. The proposed focus loss functions allow to unlock the outstanding properties of event cameras.


title: On the Occurrence of Crossings Between the Angular Distributions of Electron Neutrinos and Antineutrinos in the Supernova Core

abstract: Neutrino fast pairwise conversions have been postulated to occur in the dense core of a core-collapse supernova (SN), possibly having dramatic consequences on the SN mechanism and the observable neutrino signal. One crucial condition favoring pairwise conversions is the presence of crossings between the electron neutrino and antineutrino angular distributions (i.e., electron neutrino lepton number crossings, ELN crossings). A stationary and spherically symmetric SN toy-model is constructed to reproduce the development of the neutrino angular distributions in the dense SN core in the absence of perturbations induced by hydrodynamical instabilities. By iteratively solving the neutrino Boltzmann equations including the collisional term, our model predicts that ELN crossings can develop only in the proximity of the decoupling region and for a sharp radial evolution of the baryon density, when the electron neutrino and antineutrino number densities are comparable. Such conditions are likely to occur only in the late SN stages. Interestingly, flavor instabilities induced by spatial or temporal perturbations are unlikely to generate ELN crossings dynamically within our simplified setup.


title: Time-delay cosmographic forecasts with strong lensing and JWST stellar kinematics

abstract: We present a joint strong lensing and stellar dynamical framework for future time-delay cosmography purposes. Based on a pixelated source reconstruction and the axisymmetric Jeans equations, we are capable of constraining cosmological distances and hence the current expansion rate of the Universe ($H_0$) to the few percent level per lens, when high signal-to-noise integral field unit (IFU) observations from the next generation of telescopes become available. For illustrating the power of this method, we mock up IFU stellar kinematic data of the prominent lens system RXJ1131-1231, given the specifications of the James Webb Space Telescope. Our analysis shows that the time-delay distance ($D_{\Delta t}$) can be constrained with 3.1% uncertainty at best, if future IFU stellar kinematics are included in the fit and if the set of candidate model parameterisations contains the true lens potential. These constraints would translate to a 3.2% precision measurement on $H_0$ in flat $\Lambda$CDM cosmology from the single lens RXJ1131-1231, and can be expected to yield an $H_0$ measure with < 2.0% uncertainty, if similar gains in precision can be reached for two additional lens systems. Moreover, the angular diameter distance ($D_\mathrm{d}$) to RXJ1131-1231 can be constrained with 2.4% precision, providing two distance measurements from a single lens system, which is extremely powerful to further constrain the matter density ($\Omega_\mathrm{m}$). The measurement accuracy of $D_\mathrm{d}$, however, is highly sensitive to any systematics in the measurement of the stellar kinematics. For both distance measurements, we strongly advise to probe a large set of physically motivated lens potentials in the future, to minimise the systematic errors associated with the lens mass parameterisation.


title: High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

abstract: The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $\Lambda$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$\alpha$ versus star formation rate relation, and the H$\alpha$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.


title: Generalised model-independent characterisation of strong gravitational lenses VI: the origin of the formalism intrinsic degeneracies and their influence on $H_0$

abstract: We give a physical interpretation of the formalism intrinsic degeneracies of the gravitational lensing formalism that we derived on a mathematical basis in part IV of this series. We find that all degeneracies occur due to the partition of the mass density along the line of sight. Usually, it is partitioned into a background (cosmic) density and a foreground deflecting object. The latter can be further partitioned into a main deflecting object and perturbers. Weak deflecting objects along the line of sight are also added, either to the deflecting object or as a correction of the angular diameter distances, perturbing the cosmological background density. A priori, this is an arbitrary choice of reference frame and partition. They can be redefined without changing the lensing observables which are sensitive to the integrated deflecting mass density along the entire line of sight. Reformulating the time delay equation such that this interpretation of the degeneracies becomes easily visible, we note that the source can be eliminated from this formulation, which simplifies reconstructions of the deflecting mass distribution or the inference of the Hubble constant, $H_0$. Subsequently, we list necessary conditions to break the formalism intrinsic degeneracies and discuss ways to break them by model choices or including non-lensing observables like velocity dispersions along the line of sight with their advantages and disadvantages. We conclude with a systematic summary of all formalism intrinsic degeneracies and possibilities to break them.


title: Solvable Strong-coupling Quantum Dot Model with a Non-Fermi-liquid Pairing Transition

abstract: We show that a random interacting model exhibits solvable non-Fermi liquid behavior and exotic pairing behavior. This model, dubbed as the Yukawa-SYK model, describes the random Yukawa coupling between $M$ quantum dots each hosting $N$ flavors of fermions and $N^2$ bosons that self-tunes to criticality at low energies. The diagrammatic expansion is controlled by $1/MN$, and the results become exact in a large-$M$, large-$N$ limit. We find that pairing only develops within a region of the $(M,N)$ plane --- even though the pairing interaction is strongly attractive, the incoherence of the fermions can spoil the forming of Cooper pairs, rendering the system a non-Fermi liquid down to zero temperature. By solving the Eliashberg equation and the renormalization group equation, we show that the transition into the pairing phase exhibits Kosterlitz-Thouless quantum-critical behavior.


title: Entanglement Renormalization for Interacting Field Theories

abstract: A general method to build the entanglement renormalization (cMERA) for interacting quantum field theories is presented. We improve upon the well-known Gaussian formalism used in free theories through a class of variational non-Gaussian wavefunctionals for which expectation values of local operators can be efficiently calculated analytically and in a closed form. The method consists of a series of scale-dependent nonlinear canonical transformations on the fields of the theory under consideration. Here, the $\lambda\, \phi^4$ and the sine-Gordon scalar theories are used to illustrate how non-perturbative effects far beyond the Gaussian approximation are obtained by considering the energy functional and the correlation functions of the theory.


title: Topological phases without crystalline counterparts

abstract: Recent years saw the complete classification of topological band structures, revealing an abundance of topological crystalline insulators. Here we theoretically demonstrate the existence of topological materials beyond this framework, protected by quasicrystalline symmetries. We construct a higher-order topological phase protected by a point group symmetry that is impossible in any crystalline system. Our tight-binding model describes a superconductor on a quasicrystalline Ammann-Beenker tiling which hosts localized Majorana zero modes at the corners of an octagonal sample. The Majorana modes are protected by particle-hole symmetry and by the combination of an 8-fold rotation and in-plane reflection symmetry. We find a bulk topological invariant associated with the presence of these zero modes, and show that they are robust against large symmetry preserving deformations, as long as the bulk remains gapped. The nontrivial bulk topology of this phase falls outside all currently known classification schemes.


title: A Mechanism for Triple-Ridge Emission Structure of AGN Jets

abstract: Recent radio VLBI observations of the relativistic jet in M87 radio galaxy have shown a triple-ridge structure that consists of the conventional limb-brightened feature and a central narrow ridge. Motivated by these observations, we examine a steady axisymmetric force-free model of a jet driven by the central black hole (BH) with its electromagnetic structure being consistent with general relativistic magnetohydrodynamic simulations, and find that it can produce triple-ridge images even if we assume a simple Gaussian distribution of emitting electrons at the base of the jet. We show that the fluid velocity field associated with the electromagnetic field produces the central ridge component due to the relativistic beaming effect, while the limb-brightened feature arises due to strong magnetic field around the jet edge which also induces the electrons to be dense there. We also show that the computed image strongly depends on the electromagnetic field structure, viewing angle, and parameters related to the electrons' spatial distribution at the jet base. This study will help constraining the non-thermal electron injection mechanism of BH jets and be complementary to theoretical analyses of the upcoming data of Event Horizon Telescope.


title: High efficiency photospheric emission entailed by formation of a collimation shock in gamma-ray bursts

abstract: The primary dissipation mechanism in jets of gamma-ray bursts (GRBs), and the high efficiency of the prompt emission are long standing issues. One possibility is strong collimation of a weakly magnetized relativistic jet by the surrounding medium, which can considerably enhance the efficiency of the photospheric emission. We derive a simple analytic criterion for the radiative efficiency of a collimated jet showing that it depends most strongly on the baryon loading. We confirm this analytic result by 3D numerical simulations, and further find that mixing of jet and cocoon material at the collimation throat leads to a substantial stratification of the outflow as well as sporadic loading, even if the injected jet is uniform and continuous. One consequence of this mixing is a strong angular dependence of the radiative efficiency. Another is large differences in the Lorentz factor of different fluid elements that lead to formation of internal shocks. Our analysis indicates that in both long and short GRBs a prominent photospheric component cannot be avoided when observed within an angle of a few degrees to the axis, unless the asymptotic Lorentz factor is limited by baryon loading at the jet base to $\Gamma_\infty <100$ (with a weak dependence on outflow power). Photon generation by newly created pairs behind the collimation shock regulates the observed temperature at $\sim 50~\theta_0^{-1}$ keV, where $\theta_0$ is the initial jet opening angle, in remarkable agreement with the observed peak energies of prompt emission spectra. Further consequences for the properties of the prompt emission are discussed at the end.


title: Axion production and detection with superconducting RF cavities

abstract: We propose a novel design of a laboratory search for axions based on photon regeneration with superconducting RF cavities. Our particular setup uses a toroid as a region of confined static magnetic field, while production and detection cavities are positioned in regions of vanishing external field. This permits cavity operation at quality factors of $Q \sim 10^{10} - 10^{12}$. The limitations due to fundamental issues such as signal screening and back-reaction are discussed, and the optimal sensitivity is calculated. This experimental design can potentially probe axion-photon couplings beyond astrophysical limits, comparable and complementary to next generation optical experiments.


title: Revealing Dust Obscured Star Formation in CLJ1449+0856, a Cluster at z=2

abstract: We present SCUBA-2 450$\mu$m and 850$\mu$m data of the mature redshift 2 cluster CLJ1449. We combine this with archival Herschel data to explore the star forming properties of CLJ1449. Using high resolution ALMA and JVLA data we identify potentially confused galaxies, and use the Bayesian inference tool XID+ to estimate fluxes for them. Using archival optical and near infrared data with the energy-balance code CIGALE we calculate star formation rates, and stellar masses for all our cluster members, and find the star formation rate varies between 20-1600M$_{\odot}$yr$^{-1}$ over the entire 3Mpc radial range. The central 0.5Mpc region itself has a total star formation rate of 800$\pm$200M$_{\odot}$yr$^{-1}$, which corresponds to a star formation rate density of (1.2$\pm$0.3)$\times$10$^{4}$M$_{\odot}$yr$^{-1}$Mpc$^{-3}$, which is approximately five orders of magnitude greater than expected field values. When comparing this cluster to those at lower redshifts we find that there is an increase in star formation rate per unit volume towards the centre of the cluster. This indicates that there is indeed a reversal in the star formation/density relation in CLJ1449. Based on the radial star-formation rate density profile, we see evidence for an elevation in the star formation rate density, even out to radii of 3Mpc. At these radii the elevation could be an order of magnitude greater than field values, but the exact number cannot be determined due to ambiguity in the redshift associations. If this is the case it would imply that this cluster is still accreting material which is possibly interacting and undergoing vigorous star-formation.


title: Out of equilibrium higher-order topological insulator: Floquet engineering and quench dynamics

abstract: Higher-order topological~(HOT) states,~hosting topologically protected modes on lower-dimensional boundaries,~such as hinges and corners, have recently extended the realm of the static topological phases.~Here we demonstrate the possibility of realizing a two-dimensional \emph{Floquet} second-order topological insulator, featuring corner-localized zero quasienergy modes and characterized by quantized Floquet qudrupolar moment $Q^{\rm Flq}_{xy}=0.5$, by periodically kicking a quantum spin Hall insulator (QSHI) with a discrete fourfold ($C_4$) and time-reversal (${\mathcal T}$) symmetry breaking Dirac mass perturbation.~Furthermore, we show that $Q^{\rm Flq}_{xy}$ becomes independent of the choice of origin as the system approaches the thermodynamic limit.~We also analyze the dynamics of a corner mode after a sudden quench, when the $C_4$ and ${\mathcal T}$ symmetry breaking perturbation is switched off, and find that the corresponding survival probability displays periodic appearances of complete, partial and no revival for long time, encoding the signature of corner modes in a QSHI.~Our protocol is sufficiently general to explore the territory of dynamical HOT phases in insulators and gapless systems.

categories: cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el quant-ph


title: Machine Learning in Astronomy: a practical overview

abstract: Astronomy is experiencing a rapid growth in data size and complexity. This change fosters the development of data-driven science as a useful companion to the common model-driven data analysis paradigm, where astronomers develop automatic tools to mine datasets and extract novel information from them. In recent years, machine learning algorithms have become increasingly popular among astronomers, and are now used for a wide variety of tasks. In light of these developments, and the promise and challenges associated with them, the IAC Winter School 2018 focused on big data in Astronomy, with a particular emphasis on machine learning and deep learning techniques. This document summarizes the topics of supervised and unsupervised learning algorithms presented during the school, and provides practical information on the application of such tools to astronomical datasets. In this document I cover basic topics in supervised machine learning, including selection and preprocessing of the input dataset, evaluation methods, and three popular supervised learning algorithms, Support Vector Machines, Random Forests, and shallow Artificial Neural Networks. My main focus is on unsupervised machine learning algorithms, that are used to perform cluster analysis, dimensionality reduction, visualization, and outlier detection. Unsupervised learning algorithms are of particular importance to scientific research, since they can be used to extract new knowledge from existing datasets, and can facilitate new discoveries.


title: Exploring the CMB Power Suppression in Canonical Inflation Models

abstract: There exists some evidence of a suppression in power in the CMB multipoles around $l \sim 20-30$. If taken seriously, this is in tension with the simplest inflationary models driven by a single scalar field with a standard type of slowly varying potential function $V(\phi)$. Such potential functions generate a nearly scale invariant spectrum and so they do not possess the requisite suppression in power. In this paper we explore if canonical two-derivative inflation models, with a step-like feature in the potential, can improve agreement with data. We find that improvement can be made when one utilizes the standard slow-roll approximation formula for the power spectrum. However, we find that in order to have a feature in the power spectrum that is sufficiently localized so as to not significantly disrupt higher $l$ or lower $l$, the potential's step-like feature must be so sharp that the standard slow-roll approximations break down. This leads us to perform an exact computation of the power spectrum by solving for the Bunch-Davies mode functions numerically. We find that the corresponding CMB multipoles do not provide a good agreement with the data. We conclude that, unless there is fine-tuning, canonical inflation models do not fit this suppression in the data.


title: Anomaly Inflow for M5-branes on Punctured Riemann Surfaces

abstract: We derive the anomaly polynomials of 4d $\mathcal{N}=2$ theories that are obtained by wrapping M5-branes on a Riemann surface with arbitrary regular punctures, using anomaly inflow in the corresponding M-theory setup. Our results match the known anomaly polynomials for the 4d $\mathcal{N}=2$ class $\mathcal{S}$ SCFTs. In our approach, the contributions to the 't Hooft anomalies due to boundary conditions at the punctures are determined entirely by $G_4$-flux in the 11d geometry. This computation provides a top-down derivation of these contributions that utilizes the geometric definition of the field theories, complementing the previous field-theoretic arguments.


title: Electric Field Response In Breathing Pyrochlores

abstract: We study the effects of a uniform electric field on the the ground state and excitations of the three-dimensional U(1) spin liquid phase of a breathing pyrochlore lattice, arising due to the coupling between the conventional (Maxwell) electric field and the emergent electrodynamics of the quantum spin ice material. This is an extension of the studies for isotropic pyrochlores in Phys. Rev. B 96, 125145 (2017) to the anisotropic case, as the lattice inversion symmetry is broken in breathing pyrochlores. The emergent photons are found to exhibit birefringence, analogous to the isotropic case. However, the system possesses a nonzero polarization even in the absence of an external electric field, unlike the isotropic pyrochlore. We also find that a sufficiently strong electric field triggers a quantum phase transition into new U(1) quantum spin liquid phases which trap $\pi$-fluxes of the emergent electric field. Such transitions are seen to occur even when the applied electric field is along a direction that does not show a phase transition in the isotropic limit.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech hep-th


title: Anderson Localization on the Bethe Lattice using Cages and the Wegner Flow

abstract: Anderson localization on tree-like graphs such as the Bethe lattice, Cayley tree, or random regular graphs has attracted attention due to its apparent mathematical tractability, hypothesized connections to many-body localization, and the possibility of non-ergodic extended regimes. This behavior has been conjectured to also appear in many-body localization as a "bad metal" phase, and constitutes an intermediate possibility between the extremes of ergodic quantum chaos and integrable localization. Despite decades of research, a complete consensus understanding of this model remains elusive. Here, we use cages, maximally tree-like structures from extremal graph theory; and numerical continuous unitary Wegner flows of the Anderson Hamiltonian to develop an intuitive picture which, after extrapolating to the infinite Bethe lattice, appears to capture ergodic, non-ergodic extended, and fully localized behavior.


title: Virgo: an unlikely cluster of galaxies because of its environment

abstract: Galaxy clusters constitute powerful cosmological probes thanks to comparisons between observed and simulated clusters. As such Virgo constitutes a formidable source for detailed observations facilitated by its proximity. However, the diversity of clusters complicates the comparisons on a one-to-one basis. Simulated clusters must be carefully selected, a daunting task since most properties are unknown. Alternatively, lookalikes produced in the proper large scale environment can be used. Additionally, their statistical study give access to the mean properties of the observed cluster including its most probable history as well as its deviation from an average cluster. This paper presents such a statistical study with 200 Virgo-like and 400+ cluster-size random dark matter halos. Only 18%(0.5%) of these random halos comply within 3(2)-sigma with the mean values (radius, velocity dispersion, number of substructures, spin, velocity, concentration, center of mass offset with respect to the spherical center) of Virgo halos at z=0 and abide by a similar merging history up to redshift 4. None are within 1-sigma because of environmentally induced properties (number of substructures and velocity). For further comparisons, random halos are selected to reproduce the mass distribution of the lookalikes to cancel mass bias effects. Redshift 1 appears then as a turning point: random to Virgo-like property ratios are alternatively smaller/larger than 1. This highlights the importance of studying clusters within their proper large scale environment: simulated galaxy population, grandly affected by the cluster history, can then be compared with the observed one in details. Direct lookalikes simplify grandly the challenge.


title: Linking gas and galaxies at high redshift: MUSE surveys the environments of six damped Lyman alpha galaxies at z~3

abstract: We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.


title: Discovery of a Third Transiting Planet in the Kepler-47 Circumbinary System

abstract: Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.


title: Infrared Galaxies in the Field of the Massive Cluster Abell S1063: Discovery of a Luminous Kiloparsec-Sized HII Region in a Gravitationally Lensed IR-Luminous Galaxy at $z=0.6$

abstract: Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $\mu$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly devoid of infrared sources, showing only a few detectable sources within the central r$\sim1^{\prime}$. There is, however, one particularly bright source (2.3 mJy at 24 $\mu$m; 106 mJy at 160 $\mu$m), which corresponds to a background galaxy at $z=0.61$. The modest magnification factor (4.0$\times$) implies that this galaxy is intrinsically IR-luminous (L$_{\rm IR}=3.1\times10^{11}\ \rm L_{\odot}$). What is particularly interesting about this galaxy is that HST optical/near-infrared images show a remarkably bright and large (1 kpc) clump at one edge of the disk. Our follow-up optical/near-infrared spectroscopy shows Balmer (H$\alpha$-H8) and forbidden emission from this clump ([OII] $\lambda$3727, [OIII] $\lambda\lambda$4959,5007, [NII] $\lambda\lambda$6548,6583), indicating that it is a HII region. The HII region appears to have formed in-situ, as kinematically it is part of a rotating disk, and there is no evidence of nearby interacting galaxies. With an extinction correction of A$_{\rm V}=1.5$ mag, the star formation rate of this giant HII region is $\sim$10 M$_{\odot}$ yr$^{-1}$, which is exceptionally large, even for high redshift HII regions. Such a large and luminous HII region is often seen at $z\sim2$ but quite rare in the nearby Universe.