title: Mutated hilltop inflation revisited

abstract: In this work we re-investigate pros and cons of mutated hilltop inflation. Applying Hamilton-Jacobi formalism we solve inflationary dynamics and find that inflation goes on along the ${\cal W}_{-1}$ branch of the Lambert function. Depending on the model parameter mutated hilltop model renders two types of inflationary solutions: one corresponds to small inflaton excursion during observable inflation and the other describes large field inflation. The inflationary observables from curvature perturbation are in tune with the current data for a wide range of the model parameter. The small field branch predicts negligible amount of tensor to scalar ratio $r\sim \mathcal{O}(10^{-4})$, while the large field sector is capable of generating high amplitude for tensor perturbations, $r\sim \mathcal{O}(10^{-1})$. Also, the spectral index is almost independent of the model parameter along with a very small negative amount of scalar running. Finally we find that the mutated hilltop inflation closely resembles the $\alpha$-attractor class of inflationary models in the limit of $\alpha\phi\gg 1$.


title: Invariant Solutions for the Einstein Field Equation

abstract: In this paper we provide a method capable of producing an infinite number of solutions for Einstein's equation on static spacetimes with perfect fluid as a matter field. All spacetimes of this type which are symmetric with respect to a given group of translations and whose spatial factor is conformally flat, are characterized. We use this method to give some exact solutions of the referred equation.


title: Few-mode fibers and AO-assisted high resolution spectroscopy: coupling efficiency and modal noise mitigation

abstract: NIRPS (Near Infra-Red Planet Searcher) is an AO-assisted and fiber-fed spectrograph for high precision radial velocity measurements that will operate in the YJH-bands. While using an AO system in such instrument is generally considered to feed a single-mode fiber, NIRPS is following a different path by using a small multi-mode fiber (more specifically called "few-mode fiber"). This choice offers an excellent trade-off by allowing to design a compact cryogenic spectrograph, while maintaining a high coupling efficiency under bad seeing conditions and for faint stars. The main drawback resides in a much more important modal-noise, a problem that has to be tackled for allowing 1m/s precision radial velocity measurements. We study the impact of using an AO system to couple light into few-mode fibers. We focus on two aspects: the coupling efficiency into few-mode fibers and the question of modal noise and scrambling. We show first that NIRPS can reach coupling >= 50% up to magnitude I=12, and offer a gain of 1-2 magnitudes over a single-mode solution. We finally show that the best strategy to mitigate modal noise with the AO system is among the simplest: a continuous tip-tilt scanning of the fiber core.


title: Random walk on random planar maps: spectral dimension, resistance, and displacement

abstract: We study simple random walk on the class of random planar maps which can be encoded by a two-dimensional random walk with i.i.d. increments or a two-dimensional Brownian motion via a "mating-of-trees" type bijection. This class includes the uniform infinite planar triangulation (UIPT), the infinite-volume limits of random planar maps weighted by the number of spanning trees, bipolar orientations, or Schnyder woods they admit, and the $\gamma$-mated-CRT map for $\gamma \in (0,2)$. For each of these maps, we obtain an upper bound for the Green's function on the diagonal, an upper bound for the effective resistance to the boundary of a metric ball, an upper bound for the return probability of the random walk to its starting point after $n$ steps, and a lower bound for the graph-distance displacement of the random walk, all of which are sharp up to polylogarithmic factors. When combined with work of Lee (2017), our bound for the return probability shows that the spectral dimension of each of these random planar maps is a.s. equal to 2, i.e., the (quenched) probability that the simple random walk returns to its starting point after $2n$ steps is $n^{-1+o_n(1)}$. Our results also show that the amount of time that it takes a random walk to exit a metric ball is at least its volume (up to a polylogarithmic factor). In the special case of the UIPT, this implies that random walk typically travels at least $n^{1/4 - o_n(1)}$ units of graph distance in $n$ units of time. The matching upper bound for the displacement is proven by Gwynne and Hutchcroft (2018). These two works together resolve a conjecture of Benjamini and Curien (2013) in the UIPT case. Our proofs are based on estimates for the mated-CRT map (which come from its relationship to SLE-decorated Liouville quantum gravity) and a strong coupling of the mated-CRT map with the other random planar map models.


title: Oversampling for Imbalanced Learning Based on K-Means and SMOTE

abstract: Learning from class-imbalanced data continues to be a common and challenging problem in supervised learning as standard classification algorithms are designed to handle balanced class distributions. While different strategies exist to tackle this problem, methods which generate artificial data to achieve a balanced class distribution are more versatile than modifications to the classification algorithm. Such techniques, called oversamplers, modify the training data, allowing any classifier to be used with class-imbalanced datasets. Many algorithms have been proposed for this task, but most are complex and tend to generate unnecessary noise. This work presents a simple and effective oversampling method based on k-means clustering and SMOTE oversampling, which avoids the generation of noise and effectively overcomes imbalances between and within classes. Empirical results of extensive experiments with 71 datasets show that training data oversampled with the proposed method improves classification results. Moreover, k-means SMOTE consistently outperforms other popular oversampling methods. An implementation is made available in the python programming language.


title: A Systems Approach for Eliciting Mission-Centric Security Requirements

abstract: The security of cyber-physical systems is first and foremost a safety problem, yet it is typically handled as a traditional security problem, which means that solutions are based on defending against threats and are often implemented too late. This approach neglects to take into consideration the context in which the system is intended to operate, thus system safety may be compromised. This paper presents a systems-theoretic analysis approach that combines stakeholder perspectives with a modified version of Systems-Theoretic Accident Model and Process (STAMP) that allows decision-makers to strategically enhance the safety, resilience, and security of a cyber-physical system against potential threats. This methodology allows the capture of vital mission-specific information in a model, which then allows analysts to identify and mitigate vulnerabilities in the locations most critical to mission success. We present an overview of the general approach followed by a real example using an unmanned aerial vehicle conducting a reconnaissance mission.


title: Valence Quark and Gluon Distributions of Kaon from J/Psi Production

abstract: The only experimental information on the parton distribution of kaons was obtained from the scarce kaon-induced Drell-Yan cross section data. From these data, evidence for a valence $u$ quark distribution of kaon softer than that of pion was found. We study the feasibility to extract kaon's valence quark distribution via the existing kaon-induced J/psi production data. We compare the NA3 data on the K- / pi- and K+ / pi+ ratios for J/psi production on a platinum target with theoretical calculations based on the color-evaporation model. We show that the NA3 data on J/psi production provide an independent evidence for a valence $u$ quark distribution softer in kaon than in pion. The data also suggest different gluon distributions for kaon and pion.


title: Self-assembly of a model supramolecular polymer studied by replica exchange with solute tempering

abstract: Conventional molecular-dynamics (cMD) simulation has a well-known limitation in accessible time and length scales, and thus various enhanced sampling techniques have been proposed to alleviate the problem. In this paper we explore the utility of replica exchange with solute tempering (REST) (i.e., a variant of Hamiltonian replica exchange methods) to simulate the self-assembly of a supramolecular polymer in explicit solvent, and compare the performance with temperature-based replica exchange MD (T-REMD) as well as cMD. As a test system, we consider a relatively simple all-atom model of supramolecular polymerization (namely, benzene-1,3,5-tricarboxamides in methylcyclohexane solvent). Our results show that both REST and T-REMD are able to predict highly ordered polymer structures with helical H-bonding patterns, in contrast to cMD which completely fails to obtain such a structure for the present model. At the same time, we have also experienced some technical challenge (i.e., aggregation-dispersion transition and the resulting bottleneck for replica traversal), which is illustrated numerically. Since the computational cost of REST scales more moderately than T-REMD, we expect that REST will be useful for studying the self-assembly of larger systems in solution with enhanced rearrangement of monomers.


title: Lower Bounds for Finding Stationary Points II: First-Order Methods

abstract: We establish lower bounds on the complexity of finding $\epsilon$-stationary points of smooth, non-convex high-dimensional functions using first-order methods. We prove that deterministic first-order methods, even applied to arbitrarily smooth functions, cannot achieve convergence rates in $\epsilon$ better than $\epsilon^{-8/5}$, which is within $\epsilon^{-1/15}\log\frac{1}{\epsilon}$ of the best known rate for such methods. Moreover, for functions with Lipschitz first and second derivatives, we prove no deterministic first-order method can achieve convergence rates better than $\epsilon^{-12/7}$, while $\epsilon^{-2}$ is a lower bound for functions with only Lipschitz gradient. For convex functions with Lipschitz gradient, accelerated gradient descent achieves the rate $\epsilon^{-1}\log\frac{1}{\epsilon}$, showing that finding stationary points is easier given convexity.


title: Correlation between heavy flavour production and multiplicity in pp and p-Pb collisions at high energy in the multi-pomeron exchange model

abstract: The multiplicity dependence of heavy flavour production in pp-collisions at LHC energies is studied in the framework of the multi-pomeron exchange model. The model is introducing the string-string interaction collectivity effects in pp collisions, which modifies multiplicity and transverse momenta, leading to the non-trivial mean pt vs. multiplicity (<pt>_Nch - Nch) correlation. The string collectivity strength parameter is fixed by experimental data on multiplicity and transverse momentum correlation in a wide energy range (from ISR to LHC). The particles discrimination is implemented according to Schwinger mechanism taking into account the strong decays of hadron resonances. We demonstrate, that the faster-than-linear growth of the open charm production with the event charged particle multiplicity, observed in experimental pp high energy collisions, can be explained by the modification of the string tension due to the increasing overlap and interaction of quark-gluon strings. The model is extended for p-A interactions and the calculations for p-Pb collisions are performed.


title: Generalized Probabilistic Bisection for Stochastic Root-Finding

abstract: We consider numerical schemes for root finding of noisy responses through generalizing the Probabilistic Bisection Algorithm (PBA) to the more practical context where the sampling distribution is unknown and location-dependent. As in standard PBA, we rely on a knowledge state for the approximate posterior of the root location. To implement the corresponding Bayesian updating, we also carry out inference of oracle accuracy, namely learning the probability of correct response. To this end we utilize batched querying in combination with a variety of frequentist and Bayesian estimators based on majority vote, as well as the underlying functional responses, if available. For guiding sampling selection we investigate both Information Directed sampling, as well as Quantile sampling. Our numerical experiments show that these strategies perform quite differently; in particular we demonstrate the efficiency of randomized quantile sampling which is reminiscent of Thompson sampling. Our work is motivated by the root-finding sub-routine in pricing of Bermudan financial derivatives, illustrated in the last section of the paper.


title: Chromatic homotopy theory is asymptotically algebraic

abstract: Inspired by the Ax--Kochen isomorphism theorem, we develop a notion of categorical ultraproducts to capture the generic behavior of an infinite collection of mathematical objects. We employ this theory to give an asymptotic solution to the approximation problem in chromatic homotopy theory. More precisely, we show that the ultraproduct of the $E(n,p)$-local categories over any non-prinicipal ultrafilter on the set of prime numbers is equivalent to the ultraproduct of certain algebraic categories introduced by Franke. This shows that chromatic homotopy theory at a fixed height is asymptotically algebraic.


title: Anyonic quasiparticles of hardcore anyons

abstract: Strongly interacting topologically ordered many-body systems consisting of fermions or bosons can host exotic quasiparticles with anyonic statistics. This raises the question whether many-body systems of anyons can also form anyonic quasiparticles. Here, we show that one can, indeed, construct many-anyon wavefunctions with anyonic quasiparticles. The braiding statistics of the emergent anyons are different from those of the original anyons. We investigate hole type and particle type anyonic quasiparticles in Abelian systems on a two-dimensional lattice and compute the density profiles and braiding properties of the emergent anyons by employing Monte Carlo simulations.


title: Autoequivalences of twisted K3 surfaces

abstract: Derived equivalences of twisted K3 surfaces induce twisted Hodge isometries between them; that is, isomorphisms of their cohomologies which respect certain natural lattice structures and Hodge structures. We prove a criterion for when a given Hodge isometry arises in this way. In particular, we describe the image of the representation which associates to any autoequivalence of a twisted K3 surface its realization in cohomology: this image is a subgroup of index one or two in the group of all Hodge isometries of the twisted K3 surface. We show that both indices can occur.


title: A single-mode phononic wire

abstract: Photons and electrons transmit information to form complex systems and networks. Phonons on the other hand, the quanta of mechanical motion, are often considered only as carriers of thermal energy. Nonetheless, their flow can also be molded in fabricated nanoscale circuits. We design and experimentally demonstrate wires for phonons that transmit information with little loss or scattering across a chip. By patterning the surface of a silicon chip, we completely eliminate all but one channel of phonon conduction. We observe the emergence of low-loss standing waves in millimeter long phononic wires that we address and cool optically. Coherent transport and strong optical coupling to a phononic wire enables new phononic technologies to manipulate information and energy on a chip.


title: Variational Inference of Disentangled Latent Concepts from Unlabeled Observations

abstract: Disentangled representations, where the higher level data generative factors are reflected in disjoint latent dimensions, offer several benefits such as ease of deriving invariant representations, transferability to other tasks, interpretability, etc. We consider the problem of unsupervised learning of disentangled representations from large pool of unlabeled observations, and propose a variational inference based approach to infer disentangled latent factors. We introduce a regularizer on the expectation of the approximate posterior over observed data that encourages the disentanglement. We also propose a new disentanglement metric which is better aligned with the qualitative disentanglement observed in the decoder's output. We empirically observe significant improvement over existing methods in terms of both disentanglement and data likelihood (reconstruction quality).


title: Dissipation Factors of Spherical Current Modes on Multiple Spherical Layers

abstract: Radiation efficiencies of modal current densities distributed on a spherical shell are evaluated in terms of dissipation factor. The presented approach is rigorous, yet simple and straightforward, leading to closed-form expressions. The same approach is utilized for a two-layered shell and the results are compared with other models existing in the literature. Discrepancies in this comparison are reported and reasons are analyzed. Finally, it is demonstrated that radiation efficiency potentially benefits from the use of internal volume which contrasts with the case of the radiation Q-factor.


title: Combinatorial identities associated with new families of the numbers and polynomials and their approximation values

abstract: Recently, the numbers $Y_{n}(\lambda )$ and the polynomials $Y_{n}(x,\lambda)$ have been introduced by the second author [22]. The purpose of this paper is to construct higher-order of these numbers and polynomials with their generating functions. By using these generating functions with their functional equations and derivative equations, we derive various identities and relations including two recurrence relations, Vandermonde type convolution formula, combinatorial sums, the Bernstein basis functions, and also some well known families of special numbers and their interpolation functions such as the Apostol--Bernoulli numbers, the Apostol--Euler numbers, the Stirling numbers of the first kind, and the zeta type function. Finally, by using Stirling's approximation for factorials, we investigate some approximation values of the special case of the numbers $Y_{n}\left( \lambda \right) $.


title: Provable defenses against adversarial examples via the convex outer adversarial polytope

abstract: We propose a method to learn deep ReLU-based classifiers that are provably robust against norm-bounded adversarial perturbations on the training data. For previously unseen examples, the approach is guaranteed to detect all adversarial examples, though it may flag some non-adversarial examples as well. The basic idea is to consider a convex outer approximation of the set of activations reachable through a norm-bounded perturbation, and we develop a robust optimization procedure that minimizes the worst case loss over this outer region (via a linear program). Crucially, we show that the dual problem to this linear program can be represented itself as a deep network similar to the backpropagation network, leading to very efficient optimization approaches that produce guaranteed bounds on the robust loss. The end result is that by executing a few more forward and backward passes through a slightly modified version of the original network (though possibly with much larger batch sizes), we can learn a classifier that is provably robust to any norm-bounded adversarial attack. We illustrate the approach on a number of tasks to train classifiers with robust adversarial guarantees (e.g. for MNIST, we produce a convolutional classifier that provably has less than 5.8% test error for any adversarial attack with bounded $\ell_\infty$ norm less than $\epsilon = 0.1$), and code for all experiments in the paper is available at https://github.com/locuslab/convex_adversarial.


title: Quenched invariance principles for the maximal particle in branching random walk in random environment and the parabolic Anderson model

abstract: We consider branching random walk in spatial random branching environment (BRWRE) in dimension one, as well as related differential equations: the Fisher-KPP equation with random branching and its linearized version, the parabolic Anderson model (PAM). When the random environment is bounded, we show that after recentering and scaling, the position of the maximal particle of the BRWRE, the front of the solution of the PAM, as well as the front of the solution of the randomized Fisher-KPP equation fulfill quenched invariance principles. In addition, we prove that at time t the distance between the median of the maximal particle of the BRWRE and the front of the solution of the PAM is in O(ln t). This partially transfers results from Bramson [Comm. Pure Appl. Math. 31 (1978), no. 5, 531--581] to the setting of BRWRE.


title: Using Bernstein-Vazirani Algorithm to Attack Block Ciphers

abstract: In this paper, we study applications of Bernstein-Vazirani algorithm and present several new methods to attack block ciphers. Specifically, we first present a quantum algorithm for finding the linear structures of a function. Based on it, we propose new quantum distinguishers for the 3-round Feistel scheme and a new quantum algorithm to recover partial key of the Even-Mansour construction. Afterwards, by observing that the linear structures of a encryption function are actually high probability differentials of it, we apply our algorithm to differential analysis and impossible differential cryptanalysis respectively. We also propose a new kind of differential cryptanalysis, called quantum small probability differential cryptanalysis, based on the fact that the linear structures found by our algorithm are also the linear structure of each component function. To our knowledge, no similar method was proposed before. The efficiency and success probability of all attacks are analyzed rigorously. Since our algorithm treats the encryption function as a whole, it avoid the disadvantage of traditional differential cryptanalysis that it is difficult to extending the differential path.


title: Accretion of a Symmetry Breaking Scalar Field by a Schwarzschild Black Hole

abstract: We simulate the behaviour of a Higgs-like field in the vicinity of a Schwarzschild black hole using a highly accurate numerical framework. We consider both the limit of the zero-temperature Higgs potential, and a toy model for the time-dependent evolution of the potential when immersed in a slowly cooling radiation bath. Through these numerical investigations, we aim to improve our understanding of the non-equilibrium dynamics of a symmetry breaking field (such as the Higgs) in the vicinity of a compact object such as a black hole. Understanding this dynamics may suggest new approaches for studying properties of scalar fields using black holes as a laboratory.


title: Critical Exponents of Extremal Kerr Perturbations

abstract: We show that scalar, electromagnetic, and gravitational perturbations of extremal Kerr black holes are asymptotically self-similar under the near-horizon, late-time scaling symmetry of the background metric. This accounts for the Aretakis instability (growth of transverse derivatives) as a critical phenomenon associated with the emergent symmetry. We compute the critical exponent of each mode, which is equivalent to its decay rate. It follows from symmetry arguments that, despite the growth of transverse derivatives, all generally covariant scalar quantities decay to zero.


title: Does the existence of a plane of satellites constrain properties of the Milky Way?

abstract: According to the hierarchical model of galaxy formation underlying our current understanding of cosmology, the Milky Way (MW) has continued to accrete smaller-sized dwarf galaxies since its formation. Remnants of this process surround the MW as debris streams and satellite galaxies, and provide information that is complementary to studies of the Galaxy itself. The satellite system thus has the potential to teach us about the formation and evolution of the MW. Can the existence of a narrow, co-rotating plane of satellite galaxies (the Vast Polar Structure, VPOS) put constraints on our Galaxy's properties? Are such satellite galaxy planes more narrow around less massive hosts, more abundant around more concentrated hosts, more kinematically coherent around more early-forming halos? To address such questions, we have looked for correlations between properties of satellite galaxy planes fitted to cosmological simulations in the ELVIS suite and properties of their host dark matter halos, while accounting for realistic observational biases such as the obscuration by the disk of the MW. We find no evidence for strong correlations that would allow conclusions on the host halo properties from the mere existence of the VPOS around our Galaxy.


title: Is Self-Interacting Dark Matter Undergoing Dark Fusion?

abstract: We suggest that two-to-two dark matter fusion may be the relaxation process that resolves the small-scale structure problems of the cold collisionless dark matter paradigm. In order for the fusion cross section to scale correctly across many decades of astrophysical masses from dwarf galaxies to galaxy clusters, we require the fractional binding energy released to be greater than v^n ~ [10^{-(2-3)}]^n, where n=1,2 depends on local dark sector chemistry. The size of the dark-sector interaction cross sections must be sigma ~ 0.1-1 barn, moderately larger than for Standard Model deuteron fusion, indicating a dark nuclear scale Lambda ~ O(100 MeV). Dark fusion firmly predicts constant sigma v below the characteristic velocities of galaxy clusters. Observations of the inner structure of galaxy groups with velocity dispersion of several hundred kilometer per second, of which a handful have been identified, could differentiate dark fusion from a dark photon model.

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


title: The Nelson-Barr Relaxion

abstract: Cosmological relaxation models in which the relaxion is identified with the QCD axion, generically fail to account for the smallness of the strong CP phase. We present a simple alternative solution to this "relaxion CP problem" based on the Nelson-Barr mechanism. We take CP to be a symmetry of the UV theory, and the relaxion to have no anomalous coupling with QCD. The non-zero vacuum expectation value of the relaxion breaks CP spontaneously, and the resulting phase is mapped to the Cabibbo-Kobayashi-Maskawa phase of the Standard Model. The extended Nelson-Barr quark sector generates the relaxion "rolling" potential radiatively, relating the new physics scale with the relaxion decay constant. With no new states within the reach of the LHC, our relaxion can still be probed in a variety of astrophysical and cosmological processes, as well as in flavor experiments.


title: Hubble PanCET: An isothermal day-side atmosphere for the bloated gas-giant HAT-P-32Ab

abstract: We present a thermal emission spectrum of the bloated hot Jupiter HAT-P-32Ab from a single eclipse observation made in spatial scan mode with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST). The spectrum covers the wavelength regime from 1.123 to 1.644 microns which is binned into 14 eclipse depths measured to an averaged precision of 104 parts-per million. The spectrum is unaffected by a dilution from the close M-dwarf companion HAT-P-32B, which was fully resolved. We complemented our spectrum with literature results and performed a comparative forward and retrieval analysis with the 1D radiative-convective ATMO model. Assuming solar abundance of the planet atmosphere, we find that the measured spectrum can best be explained by the spectrum of a blackbody isothermal atmosphere with Tp = 1995 +/- 17K, but can equally-well be described by a spectrum with modest thermal inversion. The retrieved spectrum suggests emission from VO at the WFC3 wavelengths and no evidence of the 1.4 micron water feature. The emission models with temperature profiles decreasing with height are rejected at a high confidence. An isothermal or inverted spectrum can imply a clear atmosphere with an absorber, a dusty cloud deck or a combination of both. We find that the planet can have continuum of values for the albedo and recirculation, ranging from high albedo and poor recirculation to low albedo and efficient recirculation. Optical spectroscopy of the planet's day-side or thermal emission phase curves can potentially resolve the current albedo with recirculation degeneracy.


title: Hunting Faint Dwarf Galaxies in the Field Using Integrated Light Surveys

abstract: We discuss the approach of searching low mass dwarf galaxies, $\lesssim10^6\textrm{ M}_{\odot}$, in the general field, using integrated light surveys. By exploring the limiting surface brightness-spatial resolution ($\mu_{\textrm{eff,lim}}-\theta$) parameter space, we suggest that faint field dwarfs in the Local Volume, between $3$ and $10 \textrm{ Mpc}$, are expected to be detected effectively and in large numbers using integrated light photometric surveys, complementary to the classical star counts method. We use a sample of Local Group dwarf galaxies to construct relations between their photometric and structural parameters, $\textrm{M}_{*}$-$\mu_{\textrm{eff,V}}$ and $\textrm{M}_{*}$-$\textrm{R}_{\textrm{eff}}$. We use these relations, along with assumed functional forms for the halo mass function and the stellar mass-halo mass relation, to calculate the lowest detectable stellar masses in the Local Volume and the expected number of galaxies as a function of the limiting surface brightness and spatial resolution. The number of detected galaxies depends mostly on the limiting surface brightness for distances $>3\textrm{ Mpc}$ while spatial resolution starts to play a role at distances $>8\textrm{ Mpc}$. Surveys with $\mu_{\textrm{eff,lim}}\sim30\textrm{ mag arcsec}^{-2}$ should be able to detect galaxies with stellar masses down to $\sim10^4 \textrm{ M}_{\odot}$ in the Local Volume. Depending on the assumed stellar mass-halo mass relation, the expected number of galaxies between $3$ and $10\textrm{ Mpc}$ is $0.04-0.35\textrm{ deg}^{-2}$, assuming a limiting surface brightness of $\sim29-30\textrm{ mag arcsec}^{-2}$ and a spatial resolution $<4''$. We currently look for field dwarf galaxies by performing a blank wide-field survey with the Dragonfly Telephoto Array, optimized for the detection of ultra-low surface brightness structures.


title: Theory of Large Intrinsic Spin Hall Effect in Iridate Semimetals

abstract: We theoretically investigate the mechanism to generate large intrinsic spin Hall effect in iridates or more broadly in 5d transition metal oxides with strong spin-orbit coupling. We demonstrate such a possibility by taking the example of orthorhombic perovskite iridate with nonsymmorphic lattice symmetry, SrIrO$_3$, which is a three-dimensional semimetal with nodal line spectrum. It is shown that large intrinsic spin Hall effect arises in this system via the spin-Berry curvature originating from the nearly degenerate electronic spectra surrounding the nodal line. This effect exists even when the nodal line is gently gapped out, due to the persistent nearly degenerate electronic structure, suggesting a distinct robustness. The magnitude of the spin Hall conductivity is shown to be comparable to the best known example such as doped topological insulators and the biggest in any transition metal oxides. To gain further insight, we compute the intrinsic spin Hall conductivity in both of the bulk and thin film systems. We find that the geometric confinement in thin films leads to significant modifications of the electronic states, leading to even bigger spin Hall conductivity in certain cases. We compare our findings with the recent experimental report on the discovery of large spin Hall effect in SrIrO$_3$ thin films.


title: The Highly Luminous Type Ibn Supernova ASASSN-14ms

abstract: We present photometric and spectroscopic follow-up observations of the highly luminous Type Ibn supernova ASASSN-14ms, which was discovered on UT 2014-12-26.61 at $m_V \sim 16.5$. With a peak absolute $V$-band magnitude brighter than $-20.5$, a peak bolometric luminosity of $1.7 \times 10^{44}$ ergs s$^{-1}$, and a total radiated energy of $2.1 \times 10^{50}$ ergs, ASASSN-14ms is one of the most luminous Type Ibn supernovae yet discovered. In simple models, the most likely power source for this event is a combination of the radioactive decay of $^{56}$Ni and $^{56}$Co at late times and the interaction of supernova ejecta with the progenitor's circumstellar medium at early times, although we cannot rule out the possibility of a magnetar-powered light curve. The presence of a dense circumstellar medium is indicated by the intermediate-width He I features in the spectra. The faint ($m_g \sim 21.6$) host galaxy SDSS J130408.52+521846.4 has an oxygen abundance below $12+\log(O/H) \lesssim 8.3$, a stellar mass of $M_* \sim 2.6 \times 10^8 M_{\odot}$, and a star formation rate of $\textrm{SFR} \sim 0.02$ $M_{\odot}$ yr$^{-1}$.


title: From Bosonic Topological Transition to Symmetric Fermion Mass Generation

abstract: The bosonic topological transition (BTT) is a quantum critical point between the bosonic symmetry protected topological phase and the trivial phase. In this work, we derive a description of this transition in terms of compact quantum electrodynamics (QED) with four fermion flavors ($N_f=4$). This allows us to describe the transition in a lattice model with the maximal microscopic symmetry: an internal SO(4) symmetry. Within a systematic renormalization group analysis, we identify the critical point with the desired O(4) emergent symmetry and all expected deformations. By lowering the microscopic symmetry we recover the previous $N_f=2$ non-compact QED description of the BTT. Finally, by merging two BTTs we recover a previously discussed theory of symmetric mass generation, as an SU(2) quantum chromodynamics-Higgs theory with $N_f=4$ flavors of SU(2) fundamental fermions and one SU(2) fundamental Higgs boson. This provides a consistency check on both theories.


title: The String Landscape, the Swampland, and the Missing Corner

abstract: We give a brief overview of the string landscape and techniques used to construct string compactifications. We then explain how this motivates the notion of the swampland and review a number of conjectures that attempt to characterize theories in the swampland. We also compare holography in the context of superstrings with the similar, but much simpler case of topological string theory. For topological strings, there is a direct definition of topological gravity based on a sum over a "quantum gravitational foam." In this context, holography is the statement of an identification between a gravity and gauge theory, both of which are defined independently of one another. This points to a missing corner in string dualities which suggests the search for a direct definition of quantum theory of gravity rather than relying on its strongly coupled holographic dual as an adequate substitute (Based on TASI 2017 lectures given by C. Vafa).


title: Firefly (Fitting IteRativEly For Likelihood analYsis): a full spectral fitting code

abstract: We present a new spectral fitting code, Firefly, for deriving the stellar population properties of stellar systems. Firefly is a chi-squared minimisation fitting code that fits combinations of single-burst stellar population models to spectroscopic data, following an iterative best-fitting process controlled by the Bayesian Information Criterion. No priors are applied, rather all solutions within a statistical cut are retained with their weight. Moreover, no additive or multiplicative polynomia are employed to adjust the spectral shape. This fitting freedom is envisaged in order to map out the effect of intrinsic spectral energy distribution (SED) degeneracies, such as age, metallicity, dust reddening on galaxy properties, and to quantify the effect of varying input model components on such properties. Dust attenuation is included using a new procedure, which was tested on Integral Field Spectroscopic (IFS) data in a previous paper. The fitting method is extensively tested with a comprehensive suite of mock galaxies, real galaxies from the Sloan Digital Sky Survey and Milky Way globular clusters. We also assess the robustness of the derived properties as a function of signal-to-noise ratio and adopted wavelength range. We show that \FF\ is able to recover age, metallicity, stellar mass and even the star formation history remarkably well down to a $S/N\sim5$, for moderately dusty systems. Code and results are publicly available at www.icg.port.ac.uk/firefly.


title: Foreground and sensitivity analysis for broad band (2D) 21cm--Ly-alpha and 21cm--H-alpha correlation experiments probing the Epoch of Reionization

abstract: A detection of the predicted anticorrelation between 21cm and either Ly-alpha or H-alpha from the Epoch of Reionization (EOR) would be a powerful probe of the first galaxies. While 3D intensity maps isolate foregrounds in low k_\parallel modes, infrared surveys cannot yet match the field of view and redshift resolution of radio intensity mapping experiments. In contrast, 2D (i.e., broad band) infrared intensity maps can be measured with current experiments and are limited by foregrounds instead of photon or thermal noise. We show 2D experiments can measure most of the 3D fluctuation power at k<0.2 Mpc^-1 while preserving its correlation properties. However, we show foregrounds pose two challenges: (1) simple geometric effects produce percent-level correlations between radio and infrared fluxes, even if their luminosities are uncorrelated; and (2) radio and infrared foreground residuals contribute sample variance noise to the cross spectrum. The first challenge demands better foreground masking and subtraction, while the second demands large fields of view to average away uncorrelated radio and infrared power. Using radio observations from the Murchison Widefield Array and near-infrared observations from the Asteroid Terrestrial-impact Last Alert System, we set an upper limit on residual foregrounds of the 21cm--Ly-alpha cross power spectrum at z\sim7 of \Delta^2<181 kJy/sr * mK (95\%) at \ell\sim800. We predict levels of foreground correlation and sample variance noise in future experiments, showing that higher resolution surveys such as LOFAR, SKA-LOW, and the Dark Energy Survey can start to probe models of the 21cm--Ly\alpha EOR cross spectrum.


title: The (Un)reliability of saliency methods

abstract: Saliency methods aim to explain the predictions of deep neural networks. These methods lack reliability when the explanation is sensitive to factors that do not contribute to the model prediction. We use a simple and common pre-processing step ---adding a constant shift to the input data--- to show that a transformation with no effect on the model can cause numerous methods to incorrectly attribute. In order to guarantee reliability, we posit that methods should fulfill input invariance, the requirement that a saliency method mirror the sensitivity of the model with respect to transformations of the input. We show, through several examples, that saliency methods that do not satisfy input invariance result in misleading attribution.


title: Three-dimensional GRMHD simulations of neutrino-cooled accretion disks from neutron star mergers

abstract: Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from these neutrino-cooled accretion disks represent an important site for $r$-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales ($380\,\mathrm{ms}$). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an $\sim\!20\,\mathrm{ms}$ cycle, and the generation of a `hot' disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into $\alpha$-particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low ($Y_e\sim0.1$) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of $Y_e\approx 0.2$ to robustly synthesize third-peak $r$-process elements. The quasi-spherical outflows are projected to unbind $40\%$ of the initial disk mass with typical asymptotic escape velocities of $0.1c$, and may thus represent the dominant mass ejection mechanism in NS-NS mergers. Including neutrino absorption, our findings agree with previous hydrodynamical $\alpha-$disk simulations that the entire range of $r$-process nuclei from the first to the third $r$-process peak can be synthesized in the outflows, in good agreement with observed solar system abundances. The asymptotic escape velocities and the quantity of ejecta, when extrapolated to moderately higher disk masses, are consistent with those needed to explain the red kilonova emission following the NS merger GW170817.


title: Introducing a model of pairing based on Base pair specific interactions between identical DNA sequences

abstract: At present, there have been suggested two types of physical mechanism that may facilitate preferential pairing between DNA molecules, with identical or similar base pair texts, without separation of base pairs. One solely relies on base pair specific patterns of helix distortion being the same on the two molecules, discussed extensively in the past. The other mechanism proposes that there are preferential interactions between base pairs of the same composition. We introduce a model, built on this second mechanism, where both thermal stretching and twisting fluctuations are included, as well as the base pair specific helix distortions. Firstly, we consider an approximation for weak pairing interactions, or short molecules. This yields a dependence of the energy on the square root of the molecular length, which could explain recent experimental data. However, analysis suggests that this approximation is no-longer valid at large DNA lengths. In a second approximation, for long molecules, we define two adaptation lengths for twisting and stretching, over which the pairing interaction can limit the accumulation of helix disorder. When the pairing interaction is sufficiently strong, both adaptation lengths are finite; however, as we reduce pairing strength, the stretching adaptation length remains finite but the torsional one becomes infinite. This second state persists to arbitrarily weak values of the pairing strength; suggesting that, if the molecules are long enough, the pairing energy scales as length. To probe differences between the two pairing mechanisms, we also construct a model of similar form. However, now, pairing between identical sequences solely relies on the intrinsic helix distortion patterns. Between the two models, we see interesting qualitative differences. We discuss our findings, and suggest new work to distinguish between the two mechanisms.

categories: physics.bio-ph cond-mat.soft cond-mat.stat-mech physics.chem-ph


title: On the UV/optical variation in NGC 5548: new evidence against the reprocessing diagram

abstract: The reprocessing scenario is widely adopted in literature to explain the observed tight inter-band correlation and short lags in the UV/optical variations of active galactic nuclei (AGNs). In this work we look into the color variability of the famous Seyfert galaxy NGC 5548 with high quality Swift multi-band UV/optical light curves. We find the color variation of NGC 5548 is clearly timescale-dependent, in a way that it is more prominent at shorter timescales. This is similar to that previously detected in quasar samples, but for the first time in an individual AGN. We show that while a reprocessing model with strict assumptions on the driving source and the disk size can apparently match the observed light curves and inter-band lags, it fails to reproduce the observed timescale dependency in the color variation. Such discrepancy raises a severe challenge to, and can hardly be reconciled under the widely accepted reprocessing diagram. It also demonstrates that the timescale dependency of the color variation is uniquely powerful in probing the physics behind AGN UV/optical variations.


title: Revealing missing charges with generalised quantum fluctuation relations

abstract: The non-equilibrium dynamics of quantum many-body systems is one of the most fascinating problems in physics. Open questions range from how they relax to equilibrium to how to extract useful work from them. A critical point lies in assessing whether a system has conserved quantities (or 'charges'), as these can drastically influence its dynamics. Here we propose a general protocol to reveal the existence of charges based on a set of exact relations between out-of-equilibrium fluctuations and equilibrium properties of a quantum system. We apply these generalised quantum fluctuation relations to a driven quantum simulator, demonstrating their relevance to obtain unbiased temperature estimates from non-equilibrium measurements. Our findings will help guide research on the interplay of quantum and thermal fluctuations in quantum simulation, in studying the transition from integrability to chaos and in the design of new quantum devices.


title: Necessary and sufficient state condition for two-qubit steering using two measurement settings per party and monogamy of steering

abstract: We consider the Cavalcanti-Foster-Fuwa-Wiseman inequality~\cite{achsh} which is a necessary and sufficient steerability condition for two-qubit states with two measurement settings on each side. We derive the criterion which an arbitrary two-qubit state must satisfy in order to violate this inequality, and obtain its maximum attainable violation in quantum mechanics. The derived condition on the state parameters enables us to establish a tight monogamy relation for two-qubit steering.


title: Hydrodynamics of the interacting Bose gas in the Quantum Newton Cradle setup

abstract: Describing and understanding the motion of quantum gases out of equilibrium is one of the most important modern challenges for theorists. In the groundbreaking Quantum Newton Cradle experiment [Kinoshita, Wenger and Weiss, Nature 440, 900, 2006], quasi-one-dimensional cold atom gases were observed with unprecedented accuracy, providing impetus for many developments on the effects of low dimensionality in out-of-equilibrium physics. But it is only recently that the theory of generalized hydrodynamics has provided the adequate tools for a numerically efficient description. Using it, we give a complete numerical study of the time evolution of an ultracold atomic gas in this setup, in an interacting parameter regime close to that of the original experiment. We evaluate the full evolving phase-space distribution of particles. We simulate oscillations due to the harmonic trap, the collision of clouds without thermalization, and observe a small elongation of the actual oscillation period and cloud deformations due to many-body dephasing. We also analyze the effects of weak anharmonicity. In the experiment, measurements are made after release from the one-dimensional trap. We evaluate the gas density curves after such a release, characterizing the actual time necessary for reaching the asymptotic state where the integrable quasi-particle momentum distribution function emerges.


title: Optical spectroscopy of the blue supergiant Sk-69 279 and its circumstellar shell with SALT

abstract: We report the results of optical spectroscopy of the blue supergiant Sk-69 279 and its circular shell in the Large Magellanic Cloud (LMC) with the Southern African Large Telescope (SALT). We classify Sk-69 279 as an O9.2 Iaf star and analyse its spectrum by using the stellar atmosphere code CMFGEN, obtaining a stellar temperature of \approx 30 kK, a luminosity of log(L_*/Lsun)=5.54, a mass-loss rate of log(\dot{M}/(Msun/yr))=-5.26, and a wind velocity of 800 km/s. We found also that Sk-69 279 possesses an extended atmosphere with an effective temperature of \approx 24 kK and that its surface helium and nitrogen abundances are enhanced, respectively, by factors of \approx 2 and 20--30. This suggests that either Sk-69 279 was initially a (single) fast-rotating (\ga 400 km/s) star, which only recently evolved off the main sequence, or that it is a product of close binary evolution. The long-slit spectroscopy of the shell around Sk-69 279 revealed that its nitrogen abundance is enhanced by the same factor as the stellar atmosphere, which implies that the shell is composed mostly of the CNO processed material lost by the star. Our findings support previous propositions that some massive stars can produce compact circumstellar shells and, presumably, appear as luminous blue variables while they are still on the main sequence or have only recently left it.


title: Direction Dependent Corrections in Polarimetric Radio Imaging II: A-Solver Methodology A low-order solver for the A-Term of the A-Projection algorithm

abstract: The effects of the antenna far-field power pattern limits the imaging performance of modern wide-bandwidth, high-sensitivity interferometric radio telescopes. Given a model for the aperture illumination pattern (AIP) of the antenna, referred to as the A-term, the wide-band (WB) A-Projection algorithm corrects for the effects of its time, frequency, and polarization structure. The level to which this correction is possible depends on how accurately the A-term, represents the true AIP. In this paper, we describe the A-Solver methodology that combines physical modeling with optimization to holographic measurements to build an accurate model for the AIP. Using a parametrized ray-tracing code as the predictor, we solve for the frequency dependence of the antenna optics and show that the resulting low-order model for the Karl G. Jansky Very Large Array (VLA) antenna captures the dominant frequency-dependent terms. The A-Solver methodology described here is generic and can be adapted for other types of antennas as well. The parameterization is based on the physical characteristics of the antenna structure and optics and is therefore arguably a compact representation (minimized degrees of freedom) of the frequency-dependent structure of the antenna A-term. In this paper, we also show that the parameters derived from A-Solver methodology are expected to improve sensitivity and imaging performance out to the first side-lobe of the antenna.


title: No surviving companion in Kepler's supernova

abstract: We have surveyed Kepler's supernova remnant in search of the companion star of the explosion. We have gone as deep as 2.6 solar luminosities in the stars within 20% of the radius of the remnant. We use FLAMES at the VLT-UT2 telescope to obtain high resolution spectra of the stellar candidates selected from HST images. The resulting set of stellar parameters suggests that these stars come from a rather ordinary mixture of field stars (mostly giants). A few of the stars seem to have low [Fe/H] (< -1) and they are consistent with being metal-poor giants. The radial velocities and rotational velocities vrot sin i are very well determined. There are no fast rotating stars as vrot sin i < 20 km/s. The radial velocities from the spectra and the proper motions determined from HST images are compatible with those expected from the Besan\c{c}on model of the Galaxy. The strong limits placed on luminosity suggest that this supernova could have arisen either from the core-degenerate scenario or from the double-degenerate scenario.


title: Using Bayesian latent Gaussian graphical models to infer symptom associations in verbal autopsies

abstract: Learning dependence relationships among variables of mixed types provides insights in a variety of scientific settings and is a well-studied problem in statistics. Existing methods, however, typically rely on copious, high quality data to accurately learn associations. In this paper, we develop a method for scientific settings where learning dependence structure is essential, but data are sparse and have a high fraction of missing values. Specifically, our work is motivated by survey-based cause of death assessments known as verbal autopsies (VAs). We propose a Bayesian approach to characterize dependence relationships using a latent Gaussian graphical model that incorporates informative priors on the marginal distributions of the variables. We demonstrate such information can improve estimation of the dependence structure, especially in settings with little training data. We show that our method can be integrated into existing probabilistic cause-of-death assignment algorithms and improves model performance while recovering dependence patterns between symptoms that can inform efficient questionnaire design in future data collection.


title: On Polymorphic Sessions and Functions: A Tale of Two (Fully Abstract) Encodings

abstract: This work exploits the logical foundation of session types to determine what kind of type discipline for the pi-calculus can exactly capture, and is captured by, lambda-calculus behaviours. Leveraging the proof theoretic content of the soundness and completeness of sequent calculus and natural deduction presentations of linear logic, we develop the first mutually inverse and fully abstract processes-as-functions and functions-as-processes encodings between a polymorphic session pi-calculus and a linear formulation of System F. We are then able to derive results of the session calculus from the theory of the lambda-calculus: (1) we obtain a characterisation of inductive and coinductive session types via their algebraic representations in System F; and (2) we extend our results to account for value and process passing, entailing strong normalisation.


title: Weighted a priori estimates for elliptic equations

abstract: We give a simpler proof of the a priori estimates obtained in the paper by Duran, Sanmartino and Toschi for solutions of elliptic problems in weighted Sobolev norms when the weights belong to the Muckenhoupt class $A_p$. The argument is a generalization to bounded domains of the one used in $\mathbb{R}^n$ to prove the continuity of singular integral operators in weighted norms. In the case of singular integral operators it is known that the $A_p$ condition is also necessary for the continuity. We do not know whether this is also true for the a priori estimates in bounded domains but we are able to prove a weaker result when the operator is the Laplacian or a power of it. We prove that a necessary condition is that the weight belongs to the local $A_p$ class.


title: Entanglement Distribution Statistic in Andreev Billiards

abstract: We investigate statistical aspects of the entanglement production for open chaotic mesoscopic billiards in contact with superconducting parts, known as Andreev billiards. The complete distributions of concurrence and entanglement of formation are obtained by using the Altland-Zirnbauer symmetry classes of circular ensembles of scattering matrices, which complements previous studies in chaotic universal billiards belonging to other classes of random matrix theory. Our results show a unique and very peculiar behavior: the realization of entanglement in a Andreev billiard always results in non-separable state, regardless of the time reversal symmetry. The analytical calculations are supported by a numerical Monte Carlo simulation.


title: On the Steady State of Continuous Time Stochastic Opinion Dynamics with Power Law Confidence

abstract: This paper introduces a class of non-linear and continuous-time opinion dynamics model with additive noise and state dependent interaction rates between agents. The model features interaction rates which are proportional to a negative power of opinion distances. We establish a non-local partial differential equation for the distribution of opinion distances and use Mellin transforms to provide an explicit formula for the stationary solution of the latter, when it exists. Our approach leads to new qualitative and quantitative results on this type of dynamics. To the best of our knowledge these Mellin transform results are the first quantitative results on the equilibria of opinion dynamics with distance-dependent interaction rates. The closed form expressions for this class of dynamics are obtained for the two agent case. However the results can be used in mean-field models featuring several agents whose interaction rates depend on the empirical average of their opinions. The technique also applies to linear dynamics, namely with a constant interaction rate, on an interaction graph.


title: Correcting Nuisance Variation using Wasserstein Distance

abstract: Profiling cellular phenotypes from microscopic imaging can provide meaningful biological information resulting from various factors affecting the cells. One motivating application is drug development: morphological cell features can be captured from images, from which similarities between different drug compounds applied at different doses can be quantified. The general approach is to find a function mapping the images to an embedding space of manageable dimensionality whose geometry captures relevant features of the input images. An important known issue for such methods is separating relevant biological signal from nuisance variation. For example, the embedding vectors tend to be more correlated for cells that were cultured and imaged during the same week than for those from different weeks, despite having identical drug compounds applied in both cases. In this case, the particular batch in which a set of experiments were conducted constitutes the domain of the data; an ideal set of image embeddings should contain only the relevant biological information (e.g. drug effects). We develop a general framework for adjusting the image embeddings in order to `forget' domain-specific information while preserving relevant biological information. To achieve this, we minimize a loss function based on distances between marginal distributions (such as the Wasserstein distance) of embeddings across domains for each replicated treatment. For the dataset we present results with, the only replicated treatment happens to be the negative control treatment, for which we do not expect any treatment-induced cell morphology changes. We find that for our transformed embeddings (i) the underlying geometric structure is not only preserved but the embeddings also carry improved biological signal; and (ii) less domain-specific information is present.


title: Comparison of the DeWitt metric in general relativity with the fourth-rank constitutive tensors in electrodynamics and in elasticity theory

abstract: We perform a short comparison between the local and linear constitutive tensor $\chi^{\lambda\nu\sigma\kappa}$ in four-dimensional electrodynamics (Sec.2), the elasticity tensor $c^{ijkl}$ in three-dimensional elasticity theory (Sec.3), and the DeWitt metric $G^{abcd}$ in general relativity, with ${a,b,\dots=1,2,3}$ (Sec.4). We find that the DeWitt metric has only six independent components.


title: An algebraic formulation of the locality principle in renormalisation

abstract: We study the mathematical structure underlying the concept of locality which lies at the heart of classical and quantum field theory, and develop a machinery used to preserve locality during the renormalisation procedure. Viewing renormalisation in the framework of Connes and Kreimer as the algebraic Birkhoff factorisation of characters on a Hopf algebra with values in a Rota-Baxter algebra, we build locality variants of these algebraic structures, leading to a locality variant of the algebraic Birkhoff factorisation. This provides an algebraic formulation of the conservation of locality while renormalising. As an application in the context of the Euler-Maclaurin formula on cones, we renormalise the exponential generating function which sums over the lattice points in convex cones. For a suitable multivariate regularisation, renormalisation from the algebraic Birkhoff factorisation amounts to composition by a projection onto holomorphic multivariate functions.


title: Using Deep Learning to Examine the Association between the Built Environment and Neighborhood Adult Obesity Prevalence

abstract: More than one-third of the adult population in the United States is obese. Obesity has been linked to factors such as, genetics, diet, physical activity and the environment. However, evidence indicating associations between the built environment and obesity has varied across studies and geographical contexts. Here, we used deep learning and approximately 150,000 high resolution satellite images to extract features of the built environment. We then developed linear regression models to consistently quantify the association between the extracted features and obesity prevalence at the census tract level for six cities in the United States. The extracted features of the built environment explained 72% to 90% of the variation in obesity prevalence across cities. Outof-sample predictions were considerably high with correlations greater than 80% between predicted and true obesity prevalence across all census tracts. This study supports a strong association between the built environment and obesity prevalence. Additionally, it also illustrates that features of the built environment extracted from satellite images can be useful for studying health indicators, such as obesity. Understanding the association between specific features of the built environment and obesity prevalence can lead to structural changes that could encourage physical activity and decreases in obesity prevalence.


title: Invariant ergodic measures and the classification of crossed product $C^\ast$-algebras

abstract: Let $\alpha: G\curvearrowright X$ be a minimal free continuous action of an infinite countable amenable group on an infinite compact metrizable space. In this paper, under the hypothesis that the invariant ergodic probability Borel measure space $E_G(X)$ is compact and zero-dimensional, we show that the action $\alpha$ has the small boundary property. This partially answers an open problem in dynamical systems that asks whether a minimal free action of an amenable group has the small boundary property if its space $M_G(X)$ of invariant Borel probability measures forms a Bauer simplex. In addition, under the same hypothesis, we show that dynamical comparison implies almost finiteness, which was shown by Kerr to imply that the crossed product is $\mathcal{Z}$-stable. Finally, we discuss some rank properties and provide two classifiability results for crossed products, one of which is based on the work of Elliott and Niu.


title: Probing quench dynamics across a quantum phase transition into a 2D Ising antiferromagnet

abstract: Simulating the real-time evolution of quantum spin systems far out of equilibrium poses a major theoretical challenge, especially in more than one dimension. We experimentally explore the dynamics of a two-dimensional Ising spin system with transverse and longitudinal fields as we quench it across a quantum phase transition from a paramagnet to an antiferromagnet. We realize the system with a near unit-occupancy atomic array of over 200 atoms obtained by loading a spin-polarized band insulator of fermionic lithium into an optical lattice and induce short-range interactions by direct excitation to a low-lying Rydberg state. Using site-resolved microscopy, we probe the correlations in the system after a sudden quench from the paramagnetic state and compare our measurements to exact calculations in the regime where it is possible. We achieve many-body states with longer-range antiferromagnetic correlations by implementing a near-adiabatic quench and study the buildup of correlations as we cross the quantum phase transition at different rates.


title: Set-to-Set Hashing with Applications in Visual Recognition

abstract: Visual data, such as an image or a sequence of video frames, is often naturally represented as a point set. In this paper, we consider the fundamental problem of finding a nearest set from a collection of sets, to a query set. This problem has obvious applications in large-scale visual retrieval and recognition, and also in applied fields beyond computer vision. One challenge stands out in solving the problem---set representation and measure of similarity. Particularly, the query set and the sets in dataset collection can have varying cardinalities. The training collection is large enough such that linear scan is impractical. We propose a simple representation scheme that encodes both statistical and structural information of the sets. The derived representations are integrated in a kernel framework for flexible similarity measurement. For the query set process, we adopt a learning-to-hash pipeline that turns the kernel representations into hash bits based on simple learners, using multiple kernel learning. Experiments on two visual retrieval datasets show unambiguously that our set-to-set hashing framework outperforms prior methods that do not take the set-to-set search setting.


title: Structured Generative Adversarial Networks

abstract: We study the problem of conditional generative modeling based on designated semantics or structures. Existing models that build conditional generators either require massive labeled instances as supervision or are unable to accurately control the semantics of generated samples. We propose structured generative adversarial networks (SGANs) for semi-supervised conditional generative modeling. SGAN assumes the data x is generated conditioned on two independent latent variables: y that encodes the designated semantics, and z that contains other factors of variation. To ensure disentangled semantics in y and z, SGAN builds two collaborative games in the hidden space to minimize the reconstruction error of y and z, respectively. Training SGAN also involves solving two adversarial games that have their equilibrium concentrating at the true joint data distributions p(x, z) and p(x, y), avoiding distributing the probability mass diffusely over data space that MLE-based methods may suffer. We assess SGAN by evaluating its trained networks, and its performance on downstream tasks. We show that SGAN delivers a highly controllable generator, and disentangled representations; it also establishes start-of-the-art results across multiple datasets when applied for semi-supervised image classification (1.27%, 5.73%, 17.26% error rates on MNIST, SVHN and CIFAR-10 using 50, 1000 and 4000 labels, respectively). Benefiting from the separate modeling of y and z, SGAN can generate images with high visual quality and strictly following the designated semantic, and can be extended to a wide spectrum of applications, such as style transfer.


title: Multivariate stochastic integrals with respect to independently scattered random measures on {\delta}-rings

abstract: In this paper we construct general vector-valued infinite-divisible independently scattered random measures with values in $\mathbb{R}^m$ and their corresponding stochastic integrals. Moreover, given such a random measure, the class of all integrable matrix-valued deterministic functions is characterized in terms of certain characteristics of the random measure. In addition a general construction principle is presented.


title: Balas formulation for the union of polytopes is optimal

abstract: A celebrated theorem of Balas gives a linear mixed-integer formulation for the union of two nonempty polytopes whose relaxation gives the convex hull of this union. The number of inequalities in Balas formulation is linear in the number of inequalities that describe the two polytopes and the number of variables is doubled. In this paper we show that this is best possible: in every dimension there exist two nonempty polytopes such that if a formulation for the convex hull of their union has a number of inequalities that is polynomial in the number of inequalities that describe the two polytopes, then the number of additional variables is at least linear in the dimension of the polytopes. We then show that this result essentially carries over if one wants to approximate the convex hull of the union of two polytopes and also in the more restrictive setting of lift-and-project.


title: Improved Adams-type inequalities and their extremals in dimension 2m

abstract: In this paper we prove the existence of extremal functions for the Adams-Moser-Trudinger inequality on the Sobolev space $H^{m}(\Omega)$, where $\Omega$ is any bounded, smooth, open subset of $\mathbb{R}^{2m}$, $m\ge 1$. Moreover, we extend this result to improved versions of Adams' inequality of Adimurthi-Druet type. Our strategy is based on blow-up analysis for sequences of subcritical extremals and introduces several new techniques and constructions. The most important one is a new procedure for obtaining capacity-type estimates on annular regions.


title: Reply to Comment [arXiv:1708.09341] by Birse and McGovern on "Nucleon spin-averaged forward virtual Compton tensor at large $Q^2$"

abstract: We reply to misleading claims made in a comment to our work by Birse and McGovern in arXiv:1708.09341.


title: Multi-Mention Learning for Reading Comprehension with Neural Cascades

abstract: Reading comprehension is a challenging task, especially when executed across longer or across multiple evidence documents, where the answer is likely to reoccur. Existing neural architectures typically do not scale to the entire evidence, and hence, resort to selecting a single passage in the document (either via truncation or other means), and carefully searching for the answer within that passage. However, in some cases, this strategy can be suboptimal, since by focusing on a specific passage, it becomes difficult to leverage multiple mentions of the same answer throughout the document. In this work, we take a different approach by constructing lightweight models that are combined in a cascade to find the answer. Each submodel consists only of feed-forward networks equipped with an attention mechanism, making it trivially parallelizable. We show that our approach can scale to approximately an order of magnitude larger evidence documents and can aggregate information at the representation level from multiple mentions of each answer candidate across the document. Empirically, our approach achieves state-of-the-art performance on both the Wikipedia and web domains of the TriviaQA dataset, outperforming more complex, recurrent architectures.


title: Intersections of multicurves from Dynnikov coordinates

abstract: We present an algorithm for calculating the geometric intersection number of two multicurves on the $n$-punctured disk, taking as input their Dynnikov coordinates. The algorithm has complexity $O(m^2n^4)$, where $m$ is the sum of the absolute values of the Dynnikov coordinates of the two multicurves. The main ingredient is an algorithm due to Cumplido for relaxing a multicurve.


title: Notes on the starlike log--harmonic mappings of order alpha

abstract: Let $h$ and $g$ be two analytic functions in the unit disc $\Delta$ that $g(0)=1$. Also let $\beta$ be a complex number with ${\rm Re}\{\beta\}>-1/2$. A function $f$ is said to be log--harmonic mapping if it has the following representation \begin{equation*} f(z)=z |z|^{2\beta} h(z)\overline{g(z)}\quad (z\in \Delta). \end{equation*} A log--harmonic mapping $f$ is said to be starlike log--harmonic mapping of order $\alpha$, where $0\leq \alpha<1$, if \begin{equation*} {\rm Re}\left\{\frac{zf_z -\overline{z}f_{\overline{z}}}{f}\right\}>\alpha\quad(z\in \Delta). \end{equation*} In this paper, by use of the subordination principle, we study some geometric properties of the starlike log--harmonic mappings of order $\alpha$. Also, we estimate the Jacobian of log--harmonic mappings.


title: Hard Substructure of Quenched Jets: a Monte Carlo Study

abstract: Modification of the hard jet substructure in terms of the Soft Drop jet grooming algorithm observables is studied for three different scenarios of jet quenching in a quark-gluon plasma: i) an explicit enhancement of the parton splitting functions, ii) increased soft gluon emissions induced by an in-medium virtuality gain, and iii) energy loss due to a drag force. Despite the fact that first two scenarios both correspond to a radiative energy loss mechanism and lead to similar modifications of parton showers, they are shown to have very different impacts on the momentum balance of hard subjets. Simulations for heavy-ion collisions based on the second scenario are presented and found to be in a good agreement with the experimental data.


title: The role of temperature and Coulomb correlation in stabilization of CsCl-phase in FeS under pressure

abstract: The iron-sulfur system is important for planetary interiors and is intensely studied, particularly for better understanding of the cores of Mars and Earth. Yet, there is a paradox about high-pressure stability of FeS: ab initio global optimization (at DFT level) predicts a Pmmn phase (with a distorted rocksalt structure) to be stable at pressures above ~120 GPa, which has not yet been observed in the experiments that instead revealed a CsCl-type phase which, according to density functional calculations, should not be stable. Using quasiharmonic free energy calculations and the dynamical mean field theory, we show that this apparent discrepancy is removed by proper account of electron correlations and entropic effects.


title: A Model of Polarisation Rotations in Blazars from Kink Instabilities in Relativistic Jets

abstract: This paper presents a simple model of polarisation rotation in optically thin relativistic jets of blazars. The model is based on the development of helical (kink) mode of current-driven instability. A possible explanation is suggested for the observational connection between polarisation rotations and optical/gamma-ray flares in blazars, if the current-driven modes are triggered by secular increases of the total jet power. The importance of intrinsic depolarisation in limiting the amplitude of coherent polarisation rotations is demonstrated. The polarisation rotation amplitude is thus very sensitive to the viewing angle, which appears to be inconsistent with the observational estimates of viewing angles in blazars showing polarisation rotations. Overall, there are serious obstacles to explaining large-amplitude polarisation rotations in blazars in terms of current-driven kink modes.


title: Partial constraint singularities in elastic rods

abstract: We present a unified classical treatment of partially constrained elastic rods. Partial constraints often entail singularities in both shapes and reactions. Our approach encompasses both sleeve and adhesion problems, and provides simple and unambiguous derivations of counterintuitive results in the literature. Relationships between reaction forces and moments, geometry, and adhesion energies follow from the balance of energy during quasistatic motion. We also relate our approach to the balance of material momentum and the concept of a driving traction. The theory is generalizable and can be applied to a wide array of contact, adhesion, gripping, and locomotion problems.


title: Imaging a Nitrogen-Vacancy Center with a Diamond Immersion Metalens

abstract: Solid-state quantum emitters have emerged as robust single-photon sources and addressable spins: key components in rapidly developing quantum technologies for broadband magnetometry, biological sensing, and quantum information science. Performance in these applications, be it magnetometer sensitivity or quantum key generation rate, is limited by the number of photons detected. However, efficient collection of a quantum emitter's photoluminescence (PL) is challenging as its atomic scale necessitates diffraction-limited imaging with nanometer-precision alignment, oftentimes at cryogenic temperatures. In this letter, we image an individual quantum emitter, an isolated nitrogen-vacancy (NV) center in diamond, using a dielectric metalens composed of subwavelength pillars etched into the diamond's surface. The metalens eliminates the need for an objective by operating as a high-transmission-efficiency immersion lens with a numerical aperture (NA) greater than 1.0. This design provides a scalable approach for fiber coupling solid-state quantum emitters that will enable the development of deployable quantum devices.


title: Moderate maximal inequalities for the Ornstein-Uhlenbeck process

abstract: The maximal inequalities for diffusion processes have drawn increasing attention in recent years. However, the existing proof of the $L^p$ maximum inequalities for the Ornstein-Uhlenbeck process was dubious. Here we give a rigorous proof of the moderate maximum inequalities for the Ornstein-Uhlenbeck process, which include the $L^p$ maximum inequalities as special cases and generalize the remarkable $L^1$ maximum inequalities obtained by Graversen and Peskir [P. Am. Math. Soc., 128(10):3035-3041, 2000]. As a corollary, we also obtain a new moderate maximal inequality for continuous local martingales, which can be viewed as a supplement of the classical Burkholder-Davis-Gundy inequality.


title: Acceleration of tensor-product operations for high-order finite element methods

abstract: This paper is devoted to GPU kernel optimization and performance analysis of three tensor-product operators arising in finite element methods. We provide a mathematical background to these operations and implementation details. Achieving close-to-the-peak performance for these operators requires extensive optimization because of the operators' properties: low arithmetic intensity, tiered structure, and the need to store intermediate results inside the kernel. We give a guided overview of optimization strategies and we present a performance model that allows us to compare the efficacy of these optimizations against an empirically calibrated roofline.


title: Quadratic homogeneous polynomial maps $H$ and Keller maps $x+H$ with $3 \le {\rm rk} J H \le 4$

abstract: We compute by hand all quadratic homogeneous polynomial maps $H$ and all Keller maps of the form $x + H$, for which ${\rm rk} J H = 3$, over a field of arbitrary characteristic. Furthermore, we use computer support to compute Keller maps of the form $x + H$ with ${\rm rk} J H = 4$, namely: $\bullet$ all such maps in dimension $5$ over fields with $\frac12$; $\bullet$ all such maps in dimension $6$ over fields without $\frac12$. We use these results to prove the following over fields of arbitrary characteristic: for Keller maps $x + H$ for which ${\rm rk} J H \le 4$, the rows of $J H$ are dependent over the base field.


title: Sparse-View X-Ray CT Reconstruction Using $\ell_1$ Prior with Learned Transform

abstract: A major challenge in X-ray computed tomography (CT) is reducing radiation dose while maintaining high quality of reconstructed images. To reduce the radiation dose, one can reduce the number of projection views (sparse-view CT); however, it becomes difficult to achieve high-quality image reconstruction as the number of projection views decreases. Researchers have applied the concept of learning sparse representations from (high-quality) CT image dataset to the sparse-view CT reconstruction. We propose a new statistical CT reconstruction model that combines penalized weighted-least squares (PWLS) and $\ell_1$ prior with learned sparsifying transform (PWLS-ST-$\ell_1$), and a corresponding efficient algorithm based on Alternating Direction Method of Multipliers (ADMM). To moderate the difficulty of tuning ADMM parameters, we propose a new ADMM parameter selection scheme based on approximated condition numbers. We interpret the proposed model by analyzing the minimum mean square error of its ($\ell_2$-norm relaxed) image update estimator. Our results with the extended cardiac-torso (XCAT) phantom data and clinical chest data show that, for sparse-view 2D fan-beam CT and 3D axial cone-beam CT, PWLS-ST-$\ell_1$ improves the quality of reconstructed images compared to the CT reconstruction methods using edge-preserving regularizer and $\ell_2$ prior with learned ST. These results also show that, for sparse-view 2D fan-beam CT, PWLS-ST-$\ell_1$ achieves comparable or better image quality and requires much shorter runtime than PWLS-DL using a learned overcomplete dictionary. Our results with clinical chest data show that, methods using the unsupervised learned prior generalize better than a state-of-the-art deep "denoising" neural network that does not use a physical imaging model.


title: Variance-Aware Optimal Power Flow

abstract: The incorporation of stochastic loads and generation into the operation of power grids gives rise to an exposure to stochastic risk. This risk has been addressed in prior work through a variety of mechanisms, such as scenario generation or chance constraints, that can be incorporated into OPF computations. Nevertheless, numerical experiments reveal that the resulting operational decisions can produce power flows with very high variance. In this paper we introduce a variety of convex variants of OPF that explicitly address the interplay of (power flow) variance with cost minimization, and present numerical experiments that highlight our contributions.


title: Analysis of a gauged model with a spin-$\frac{1}{2}$ field directly coupled to a Rarita-Schwinger spin-$\frac{3}{2}$ field

abstract: We give a detailed analysis of an abelianized gauge field model in which a Rarita-Schwinger spin-$\frac{3}{2}$ field is directly coupled to a spin-$\frac{1}{2}$ field. The model permits a perturbative expansion in powers of the gauge field coupling, and from the Feynman rules for the model we calculate the chiral anomaly.


title: Feynman-Kac formula for the stochastic Bessel operator

abstract: We introduce a stochastic process and functional that should describe the semigroup generated by the stochastic Bessel operator. Recently Gorin and Shkolnikov showed that the largest eigenvalues for certain random matrix ensembles with soft edge behavior can be understood by analyzing large powers of tridiagonal matrices, which converge to operators in the stochastic Airy semigroup. In this article we make some progress towards realizing Gorin and Shkolnikov's program at the random matrix hard edge. We analyze large powers of a suitable tridiagonal matrix model (a slight modification of the $\beta$-Laguerre ensemble). For finite $n$ we represent the matrix powers using Feynman-Kac type formulas, which identifies a sequence of stochastic processes $X^n$ and functionals $\Phi_n$. We show that $\Phi_n(X^n)$ converges in probability to the limiting functional $\Phi(X)$ for our proposed stochastic Bessel semigroup. We also discuss how the semigroup method may be used to understand transitions from a hard edge to a soft edge in the $\beta$-Laguerre models.


title: Weight-Based Variable Ordering in the Context of High-Level Consistencies

abstract: Dom/wdeg is one of the best performing heuristics for dynamic variable ordering in backtrack search [Boussemart et al., 2004]. As originally defined, this heuristic increments the weight of the constraint that causes a domain wipeout (i.e., a dead-end) when enforcing arc consistency during search. "The process of weighting constraints with dom/wdeg is not defined when more than one constraint lead to a domain wipeout [Vion et al., 2011]." In this paper, we investigate how weights should be updated in the context of two high-level consistencies, namely, singleton (POAC) and relational consistencies (RNIC). We propose, analyze, and empirically evaluate several strategies for updating the weights. We statistically compare the proposed strategies and conclude with our recommendations.


title: Nonrelativistic trace and diffeomorphism anomalies in particle number background

abstract: Using the heat kernel method, we compute nonrelativistic trace anomalies for Schr\"odinger theories in flat spacetime, with a generic background gauge field for the particle number symmetry, both for a free scalar and a free fermion. The result is genuinely nonrelativistic, and it has no counterpart in the relativistic case. Contrary to the naive expectations, the anomaly is not gauge-invariant; this is similar to the non-gauge covariance of the non-abelian relativistic anomaly. We also show that, in the same background, the gravitational anomaly for a nonrelativistic scalar vanishes.


title: Non-standard interactions and neutrinos from dark matter annihilation in the Sun

abstract: We perform an analysis of the influence of non-standard neutrino interactions (NSI) on neutrino signal from dark matter annihilations in the Sun. Taking experimentally allowed benchmark values for the matter NSI parameters we show that the evolution of such neutrinos with energies at GeV scale can be considerably modified. We simulate propagation of neutrinos from the Sun to the Earth for realistic dark matter annihilation channels and find that the matter NSI can result in at most 30% correction to the signal rate of muon track events at neutrino telescopes. Still present experimental bounds on dark matter from these searches are robust in the presence of NSI within considerable part of their allowed parameter space. At the same time electron neutrino flux from dark matter annihilation in the Sun can be changed by a factor of few.


title: Conditional fiducial models

abstract: The fiducial is not unique in general, but we prove that in a restricted class of models it is uniquely determined by the sampling distribution of the data. It depends in particular not on the choice of a data generating model. The arguments lead to a generalization of the classical formula found by Fisher (1930). The restricted class includes cases with discrete distributions, the case of the shape parameter in the Gamma distribution, and also the case of the correlation coefficient in a bivariate Gaussian model. One of the examples can also be used in a pedagogical context to demonstrate possible difficulties with likelihood-, Bayesian-, and bootstrap-inference. Examples that demonstrate non-uniqueness are also presented. It is explained that they can be seen as cases with restrictions on the parameter space. Motivated by this the concept of a conditional fiducial model is introduced. This class of models includes the common case of iid samples from a one-parameter model investigated by Hannig (2013), the structural group models investigated by Fraser (1968), and also certain models discussed by Fisher (1973) in his final writing on the subject.


title: Does Phase Matter For Monaural Source Separation?

abstract: The "cocktail party" problem of fully separating multiple sources from a single channel audio waveform remains unsolved. Current biological understanding of neural encoding suggests that phase information is preserved and utilized at every stage of the auditory pathway. However, current computational approaches primarily discard phase information in order to mask amplitude spectrograms of sound. In this paper, we seek to address whether preserving phase information in spectral representations of sound provides better results in monaural separation of vocals from a musical track by using a neurally plausible sparse generative model. Our results demonstrate that preserving phase information reduces artifacts in the separated tracks, as quantified by the signal to artifact ratio (GSAR). Furthermore, our proposed method achieves state-of-the-art performance for source separation, as quantified by a mean signal to interference ratio (GSIR) of 19.46.


title: Spinor frame formalism for amplitudes and constrained superamplitudes of 10D SYM and 11D supergravity

abstract: We show that 10D spinor helicity formalism can be understood as spinor moving frame approach to supersymmetric particles extended to the description of amplitudes. This allows us to develop the spinor helicity formalism for 11D supergravity and a new constrained superfield formalism for 10D SYM and 11D SUGRA amplitudes. We show how the constrained on-shell superfields, one-particle counterparts of the superamplitudes, can be obtained by quantization of massless superparticle mechanics. We make some stages towards the calculation of amplitudes of 10D SYM and 11D SUGRA in this framework. In particular we have found supersymmetric Ward identities for constrained amplitudes and an especially convenient gauge, fixed on the spinor frame variables corresponding to scattered particles, which promises to be an extremely useful tool for further development of our approach. We also discuss a candidate for generalization of the BCFW recurrent relations for the constrained tree superamplitudes, indicate and discuss a problem of dependence of the expressions obtained with it on a deformation vector, which is not fixed uniquely in higher dimensional $D>4$ cases.


title: Star formation is boosted (and quenched) from the inside out: radial star formation profiles from MaNGA

abstract: The tight correlation between total galaxy stellar mass and star formation rate (SFR) has become known as the star forming main sequence. Using ~487,000 spaxels from galaxies observed as part of the Sloan Digital Sky Survey Mapping Galaxies at Apache Point Observatory (MaNGA) survey, we confirm previous results that a correlation also exists between the surface densities of star formation (Sigma_SFR) and stellar mass (Sigma_mass) on kpc scales, representing a `resolved' main sequence. Using a new metric (Delta Sigma_SFR), which measures the relative enhancement or deficit of star formation on a spaxel-by-spaxel basis relative to the resolved main sequence, we investigate the SFR profiles of 864 galaxies as a function of their position relative to the global star forming main sequence (Delta SFR). For galaxies above the global main sequence (positive Delta SFR) Delta Sigma_SFR is elevated throughout the galaxy, but the greatest enhancement in star formation occurs at small radii (< 3 kpc, or 0.5 R_e). Moreover, galaxies that are at least a factor of three above the main sequence show diluted gas phase metallicities out to 2 R_e, indicative of metal-poor gas inflows accompanying the starbursts. For quiescent/passive galaxies that lie at least a factor of 10 below the star forming main sequence there is an analogous deficit of star formation throughout the galaxy with the lowest values of Delta Sigma_SFR in the central 3 kpc. Our results are in qualitative agreement with the `compaction' scenario in which a central starburst leads to mass growth in the bulge and may ultimately precede galactic quenching from the inside-out.


title: Converged GW quasiparticle energies for transition metal oxide perovskites

abstract: The ab initio calculation of quasiparticle (QP) energies is a technically and computationally challenging problem. In condensed matter physics the most widely used approach to determine QP energies is the GW approximation. Although the GW method has been widely applied to many typical semiconductors and insulators, its application to more complex compounds such as transition metal oxide perovskites has been comparatively rare, and its proper use is not well established from a technical point of view. In this work, we have applied the single-shot G0W0 method to a representative set of transition metal oxide perovskites including 3d (SrTiO3, LaScO3, SrMnO3, LaTiO3, LaVO3, LaCrO3, LaMnO3, and LaFeO3), 4d (SrZrO3, SrTcO3, and Ca2RuO4) and 5d (SrHfO3, KTaO3 and NaOsO3) compounds with different electronic configurations, magnetic orderings, structural characteristics and bandgaps ranging from 0.1 to 6.1 eV. We discuss the proper procedure to obtain well converged QP energies and accurate bandgaps within single-shot G0W0 by comparing the conventional approach based on an incremental variation of a specific set of parameters (number of bands, energy cutoff for the plane-wave expansion and number of k-points and the basis-set extrapolation scheme [Phys. Rev. B 90, 075125 (2014)]. In addition, we have inspected the difference between the adoption of norm-conserving and ultrasoft potentials in GW calculations. A minimal statistical analysis indicates that the correlation of the GW data with the DFT gap is more robust than the correlation with the experimental gaps; moreover we identify the static dielectric constant as alternative useful parameter for the approximation of GW gap in high-throughput automatic procedures. Finally, we compute the QP band structure and spectra within the random phase approximation and compare the results with available experimental data.


title: Baryon Content in a Sample of 91 Galaxy Clusters Selected by the South Pole Telescope at 0.2 < z < 1.25

abstract: We estimate total mass ($M_{500}$), intracluster medium (ICM) mass ($M_{\mathrm{ICM}}$) and stellar mass ($M_{\star}$) in a Sunyaev-Zel'dovich effect (SZE) selected sample of 91 galaxy clusters with masses $M_{500}\gtrsim2.5\times10^{14}M_{\odot}$ and redshift $0.2 < z < 1.25$ from the 2500 deg$^2$ South Pole Telescope SPT-SZ survey. The total masses $M_{500}$ are estimated from the SZE observable, the ICM masses $M_{\mathrm{ICM}}$ are obtained from the analysis of $Chandra$ X-ray observations, and the stellar masses $M_{\star}$ are derived by fitting spectral energy distribution templates to Dark Energy Survey (DES) $griz$ optical photometry and $WISE$ or $Spitzer$ near-infrared photometry. We study trends in the stellar mass, the ICM mass, the total baryonic mass and the cold baryonic fraction with cluster mass and redshift. We find significant departures from self-similarity in the mass scaling for all quantities, while the redshift trends are all statistically consistent with zero, indicating that the baryon content of clusters at fixed mass has changed remarkably little over the past $\approx9$ Gyr. We compare our results to the mean baryon fraction (and the stellar mass fraction) in the field, finding that these values lie above (below) those in cluster virial regions in all but the most massive clusters at low redshift. Using a simple model of the matter assembly of clusters from infalling groups with lower masses and from infalling material from the low density environment or field surrounding the parent halos, we show that the measured mass trends without strong redshift trends in the stellar mass scaling relation could be explained by a mass and redshift dependent fractional contribution from field material. Similar analyses of the ICM and baryon mass scaling relations provide evidence for the so-called "missing baryons" outside cluster virial regions.


title: Quotients of graph operators by symmetry representations

abstract: A finite dimensional operator that commutes with some symmetry group admits quotient operators, which are determined by the choice of associated representation. Taking the quotient isolates the part of the spectrum supporting the chosen representation and reduces the complexity of the spectral problem. Yet, such a quotient operator is not uniquely defined. Here we present a computationally simple way of choosing a special basis for the space of intertwiners, allowing us to construct a quotient that reflects the structure of the original operator. This quotient construction generalizes previous definitions for discrete graphs, which either dealt with restricted group actions or only with the trivial representation. We also extend the method to quantum graphs, which simplifies previous constructions within this context, answers an open question regarding self-adjointness and offers alternative viewpoints in terms of a scattering approach. Applications to isospectrality are discussed, together with numerous examples and comparisons with previous results.


title: Local Noncollinear Spin Analysis

abstract: In this work, we generalize the local spin analysis of Clark and Davidson [J. Chem. Phys. 115(16), 7382 (2001)] for the partitioning of the expectation value of the molecular spin square operator, $\langle S^2 \rangle$, into atomic contributions, $\langle S_A \cdot S_B \rangle$, to the noncollinear spin case in the framework of density functional theory (DFT). We derive the working equations and we show applications to the analysis of the noncollinear spin solutions of typical spin-frustrated systems and to the calculation of magnetic exchange couplings. In the former case, we employ the triangular H$_3$He$_3$ test molecule and a Mn$_3$ complex to show that the local spin analysis provides additional information that complements the standard one-particle spin population analysis. For the calculation of magnetic exchange couplings, $J_{AB}$, we employ the local spin partitioning to extract $\langle S_A \cdot S_B \rangle$ as a function of the interatomic spin orientation given by the angle $\theta$. This, combined with the dependence of the electronic energy with $\theta$, provides a methodology to extract $J_{AB}$ from DFT calculations that, in contrast to conventional energy differences based methods, does not require the use of $ad-hoc$ $S_A$ and $S_B$ values.


title: A triple boundary lemma for surface homeomorphisms

abstract: Given an orientation-preserving and area-preserving homeomorphism $f$ of the sphere, we prove that every point which is in the common boundary of three pairwise disjoint invariant open topological disks must be a fixed point. As an application, if $K$ is an invariant Wada type continuum, then $f^n|_K$ is the identity for some $n>0$. Another application is an elementary proof of the fact that invariant disks for a nonwandering homeomorphisms homotopic to the identity in an arbitrary surface are homotopically bounded if the fixed point set is inessential. The main results in this article are self-contained.