title: Mixed Quantum States with Variable Planck Constant

abstract: Recent cosmological measurements tend to confirm that the fine structure constant {\alpha} is not immutable and has undergone a tiny variation since the Big Bang. Choosing adequate units, this could also reflect a variation of Planck's constant h. The aim of this Letter is to explore some consequences of such a possible change of h for the pure and mixed states of quantum mechanics. Surprisingly enough it is found that not only is the purity of a state extremely sensitive to such changes, but that quantum states can evolve into classical states, and vice versa. A complete classification of such transitions is however not possible for the moment being because of yet unsolved mathematical difficulties related to the study of positivity properties of trace class operators.


title: Free products in AQFT

abstract: We apply the free product construction to various local algebras in algebraic quantum field theory. If we take the free product of infinitely many identical half-sided modular inclusions with ergodic canonical endomorphism, we obtain a half-sided modular inclusion with ergodic canonical endomorphism and trivial relative commutant. On the other hand, if we take M\"obius covariant nets with trace class property, we are able to construct an inclusion of free product von Neumann algebras with large relative commutant, by considering either a finite family of identical inclusions or an infinite family of inequivalent inclusions. In two dimensional spacetime, we construct Borchers triples with trivial relative commutant by taking free products of infinitely many, identical Borchers triples. Free products of finitely many Borchers triples are possibly associated with Haag-Kastler net having S-matrix which is nontrivial and non asymptotically complete, yet the nontriviality of double cone algebras remains open.


title: Orbiton-magnon interplay in the spin-orbital polarons of KCuF3 and LaMnO3

abstract: We present a quasi-analytical solution of a spin-orbital model of KCuF$_{3}$, using the variational method for Green's functions. By analyzing the spectra for different partial bosonic compositions as well as the full solution, we show that hole propagation needs both orbiton and magnon excitations to develop, but the orbitons dominate the picture. We further elucidate the role of the different bosons by analyzing the self-energies for simplified models, establishing that because of the nature of the spin-orbital ground state, magnons alone do not produce a full quasiparticle band, in contrast to orbitons. Finally, using the electron-hole transformation between the $e_g$ states of KCuF$_3$ and LaMnO$_3$ we suggest the qualitative scenario for photoemission experiments in LaMnO$_3$.


title: Greenlees-May Duality in a Nutshell

abstract: This expository article delves into the Greenlees-May Duality Theorem which is widely thought of as a far-reaching generalization of the Grothendieck's Local Duality Theorem. This theorem is not addressed in the literature as it merits and its proof is indeed a tangled web in a series of scattered papers. By carefully scrutinizing the requisite tools, we present a clear-cut well-documented proof of this theorem for the sake of bookkeeping.


title: Influence of electrical and magnetic fields on the excess current in HTSC - normal metal point contacts

abstract: A polarity-dependent reversible change in the current-voltage characteristics between states corresponding to different values of the excess current Iexc is observed for bias voltages of several hundred millivolts in YBaCuO-Ag point contacts in the current-carrying state.


title: Evidence for Majorana bound state in an iron-based superconductor

abstract: The search for Majorana bound state (MBS) has recently emerged as one of the most active research areas in condensed matter physics, fueled by the prospect of using its non-Abelian statistics for robust quantum computation. A highly sought-after platform for MBS is two-dimensional topological superconductors, where MBS is predicted to exist as a zero-energy mode in the core of a vortex. A clear observation of MBS, however, is often hindered by the presence of additional low-lying bound states inside the vortex core. By using scanning tunneling microscope on the newly discovered superconducting Dirac surface state of iron-based superconductor FeTe1-xSex (x = 0.45, superconducting transition temperature Tc = 14.5 K), we clearly observe a sharp and non-split zero-bias peak inside a vortex core. Systematic studies of its evolution under different magnetic fields, temperatures, and tunneling barriers strongly suggest that this is the case of tunneling to a nearly pure MBS, separated from non-topological bound states which is moved away from the zero energy due to the high ratio between the superconducting gap and the Fermi energy in this material. This observation offers a new, robust platform for realizing and manipulating MBSs at a relatively high temperature.


title: Quantum transduction with adaptive control

abstract: Quantum transducers play a crucial role in hybrid quantum networks. A good quantum transducer can faithfully convert quantum signals from one mode to another with minimum decoherence. Most investigations of quantum transduction are based on the protocol of direct mode conversion. However, the direct protocol requires the matching condition, which in practice is not always feasible. Here we propose an adaptive protocol for quantum transducers, which can convert quantum signals without requiring the matching condition. The adaptive protocol only consists of Gaussian operations, feasible in various physical platforms. Moreover, we show that the adaptive protocol can be robust against imperfections associated with finite squeezing, thermal noise, and homodyne detection. It can be implemented to realize quantum state transfer between microwave and optical modes.


title: MHD simulations of the eruption of coronal flux ropes under coronal streamers

abstract: Using three-dimensional magnetohydrodynamic (MHD) simulations, we investigate the eruption of coronal flux ropes underlying coronal streamers and the development of a prominence eruption. We initialize a quasi-steady solution of a coronal helmet streamer, into which we impose at the lower boundary the slow emergence of a part of a twisted magnetic torus. As a result a quasi-equilibrium flux rope is built up under the streamer. With varying sizes of the streamer and the different length and total twist of the emerged flux rope, we found different scenarios for the evolution from quasi-equilibrium to eruption. In the cases with a broad streamer, the flux rope remains well confined until there is sufficient twist such that it first develops the kink instability and evolves through a sequence of kinked, confined states with increasing height until it eventually develops a "hernia-like" ejective eruption. For the significantly twisted flux ropes, prominence condensations form in the dips of the twisted field lines due to run-away radiative cooling. Once formed, the prominence carrying field becomes significantly non-force-free due to the prominence weight despite being low plasma $\beta$. As the flux rope erupts, we obtain the eruption of the prominence, which shows substantial draining along the legs of the erupting flux rope. The prominence may not show a kinked morphology even though the flux rope becomes kinked. On the other hand, in the case with a narrow streamer, the flux rope with less than 1 wind of twist can erupt via the onset of the torus instability.


title: Explosive Percolation on Directed Networks Due to Monotonic Flow of Activity

abstract: An important class of real-world networks have directed edges, and in addition, some rank ordering on the nodes, for instance the "popularity" of users in online social networks. Yet, nearly all research related to explosive percolation has been restricted to undirected networks. Furthermore, information on such rank ordered networks typically flows from higher ranked to lower ranked individuals, such as follower relations, replies and retweets on Twitter. Here we introduce a simple percolation process on an ordered, directed network where edges are added monotonically with respect to the rank ordering. We show with a numerical approach that the emergence of a dominant strongly connected component appears to be discontinuous. Large scale connectivity occurs at very high density compared with most percolation processes, and this holds not just for the strongly connected component structure but for the weakly connected component structure as well. We present analysis with branching processes which explains this unusual behavior and gives basic intuition for the underlying mechanisms. We also show that before the emergence of a dominant strongly connected component, multiple giant strongly connected components may exist simultaneously. By adding a competitive percolation rule with a small bias to link uses of similar rank, we show this leads to formation of two distinct components, one of high ranked users, and one of low ranked users, with little flow between the two components.


title: Correlators in the $\mathcal{N}=2$ Supersymmetric SYK Model

abstract: We study correlation functions in the one-dimensional $\mathcal{N}=2$ supersymmetric SYK model. The leading order 4-point correlation functions are computed by summing over ladder diagrams expanded in a suitable basis of conformal eigenfunctions. A novelty of the $\mathcal{N}=2$ model is that both symmetric and antisymmetric eigenfunctions are required. Although we use a component formalism, we verify that the operator spectrum and 4-point functions are consistent with $\mathcal{N}=2$ supersymmetry. We also confirm the maximally chaotic behavior of this model and comment briefly on its 6-point functions.


title: On Horizontal Recurrent Finsler Connections

abstract: In this paper we adopt the pullback approach to global Finsler geometry. We investigate horizontally recurrent Finsler connections. We prove that for each scalar ($\pi$)1-form $A$, there exists a unique horizontally recurrent Finsler connection whose $h$-recurrence form is $A$. This result generalizes the existence and uniqueness theorem of Cartan connection. We then study some properties of a special kind of horizontally recurrent Finsler connection, which we call special HRF-connection.


title: On the Raman Instability in Degenerate Relativistic Plasmas

abstract: The stimulated Raman scattering instability in a fully degenerate electron plasma is studied applying relativistic hydrodynamic and Maxwell equations. We demonstrated that the instability develops for weakly as well as strongly relativistic degenerate plasma. It is shown that in the field of strong radiation relativistically degenerate plasma effectively responses as in the case of weak degeneracy.


title: Endoscopic Depth Measurement and Super-Spectral-Resolution Imaging

abstract: Intra-operative measurements of tissue shape and multi/ hyperspectral information have the potential to provide surgical guidance and decision making support. We report an optical probe based system to combine sparse hyperspectral measurements and spectrally-encoded structured lighting (SL) for surface measurements. The system provides informative signals for navigation with a surgical interface. By rapidly switching between SL and white light (WL) modes, SL information is combined with structure-from-motion (SfM) from white light images, based on SURF feature detection and Lucas-Kanade (LK) optical flow to provide quasi-dense surface shape reconstruction with known scale in real-time. Furthermore, "super-spectral-resolution" was realized, whereby the RGB images and sparse hyperspectral data were integrated to recover dense pixel-level hyperspectral stacks, by using convolutional neural networks to upscale the wavelength dimension. Validation and demonstration of this system is reported on ex vivo/in vivo animal/ human experiments.


title: Preuve simplifi\'ee du th\'eor\`eme de Serret sur les nombres \'equivalents

abstract: If the continued fractions of two irrational numbers have a common complete quotient, then these two numbers are in the same orbit under the action of $\mathrm{PGL}(2,\mathbb{Z})$. The converse is Serret's well-known theorem, but we give a simpler and purely algebraic proof of it.


title: Towards Deep Learning Models Resistant to Adversarial Attacks

abstract: Recent work has demonstrated that deep neural networks are vulnerable to adversarial examples---inputs that are almost indistinguishable from natural data and yet classified incorrectly by the network. In fact, some of the latest findings suggest that the existence of adversarial attacks may be an inherent weakness of deep learning models. To address this problem, we study the adversarial robustness of neural networks through the lens of robust optimization. This approach provides us with a broad and unifying view on much of the prior work on this topic. Its principled nature also enables us to identify methods for both training and attacking neural networks that are reliable and, in a certain sense, universal. In particular, they specify a concrete security guarantee that would protect against any adversary. These methods let us train networks with significantly improved resistance to a wide range of adversarial attacks. They also suggest the notion of security against a first-order adversary as a natural and broad security guarantee. We believe that robustness against such well-defined classes of adversaries is an important stepping stone towards fully resistant deep learning models. Code and pre-trained models are available at https://github.com/MadryLab/mnist_challenge and https://github.com/MadryLab/cifar10_challenge.


title: Solving Integer Linear Programs with a Small Number of Global Variables and Constraints

abstract: Integer Linear Programming (ILP) has a broad range of applications in various areas of artificial intelligence. Yet in spite of recent advances, we still lack a thorough understanding of which structural restrictions make ILP tractable. Here we study ILP instances consisting of a small number of "global" variables and/or constraints such that the remaining part of the instance consists of small and otherwise independent components; this is captured in terms of a structural measure we call fracture backdoors which generalizes, for instance, the well-studied class of N -fold ILP instances. Our main contributions can be divided into three parts. First, we formally develop fracture backdoors and obtain exact and approximation algorithms for computing these. Second, we exploit these backdoors to develop several new parameterized algorithms for ILP; the performance of these algorithms will naturally scale based on the number of global variables or constraints in the instance. Finally, we complement the developed algorithms with matching lower bounds. Altogether, our results paint a near-complete complexity landscape of ILP with respect to fracture backdoors.


title: On the accuracy of ancestral sequence reconstruction for ultrametric trees with parsimony

abstract: We examine a mathematical question concerning the reconstruction accuracy of the Fitch algorithm for reconstructing the ancestral sequence of the most recent common ancestor given a phylogenetic tree and sequence data for all taxa under consideration. In particular, for the symmetric 4-state substitution model which is also known as Jukes-Cantor model, we answer affirmatively a conjecture of Li, Steel and Zhang which states that for any ultrametric phylogenetic tree and a symmetric model, the Fitch parsimony method using all terminal taxa is more accurate, or at least as accurate, for ancestral state reconstruction than using any particular terminal taxon or any particular pair of taxa. This conjecture had so far only been answered for two-state data by Fischer and Thatte. Here, we focus on answering the biologically more relevant case with four states, which corresponds to ancestral sequence reconstruction from DNA or RNA data.


title: An exponential lower bound for cut sparsifiers in planar graphs

abstract: Given an edge-weighted graph $G$ with a set $Q$ of $k$ terminals, a mimicking network is a graph with the same set of terminals that exactly preserves the sizes of minimum cuts between any partition of the terminals. A natural question in the area of graph compression is to provide as small mimicking networks as possible for input graph $G$ being either an arbitrary graph or coming from a specific graph class. In this note we show an exponential lower bound for cut mimicking networks in planar graphs: there are edge-weighted planar graphs with $k$ terminals that require $2^{k-2}$ edges in any mimicking network. This nearly matches an upper bound of $O(k 2^{2k})$ of Krauthgamer and Rika [SODA 2013, arXiv:1702.05951] and is in sharp contrast with the $O(k^2)$ upper bound under the assumption that all terminals lie on a single face [Goranci, Henzinger, Peng, arXiv:1702.01136]. As a side result we show a hard instance for the double-exponential upper bounds given by Hagerup, Katajainen, Nishimura, and Ragde~[JCSS 1998], Khan and Raghavendra~[IPL 2014], and Chambers and Eppstein~[JGAA 2013].


title: Aztec: A Platform to Render Biomedical Software Findable, Accessible, Interoperable, and Reusable

abstract: Precision medicine and health requires the characterization and phenotyping of biological systems and patient datasets using a variety of data formats. This scenario mandates the centralization of various tools and resources in a unified platform to render them Findable, Accessible, Interoperable, and Reusable (FAIR Principles). Leveraging these principles, Aztec provides the scientific community with a new platform that promotes a long-term, sustainable ecosystem of biomedical research software. Aztec is available at https://aztec.bio and its source code is hosted at https://github.com/BD2K-Aztec.


title: Inter-regional ECoG correlations predicted by communication dynamics, geometry, and correlated gene expression

abstract: Electrocorticography (ECoG) provides direct measurements of synchronized postsynaptic potentials at the exposed cortical surface. Patterns of signal covariance across ECoG sensors have been associated with diverse cognitive functions and remain a critical marker of seizure onset, progression, and termination. Yet, a systems level understanding of these patterns (or networks) has remained elusive, in part due to variable electrode placement and sparse cortical coverage. Here, we address these challenges by constructing inter-regional ECoG networks from multi-subject recordings, demonstrate similarities between these networks and those constructed from blood-oxygen-level-dependent signal in functional magnetic resonance imaging, and predict network topology from anatomical connectivity, interregional distance, and correlated gene expression patterns. Our models accurately predict out-of-sample ECoG networks and perform well even when fit to data from individual subjects, suggesting shared organizing principles across persons. In addition, we identify a set of genes whose brain-wide co-expression is highly correlated with ECoG network organization. Using gene ontology analysis, we show that these same genes are enriched for membrane and ion channel maintenance and function, suggesting a molecular underpinning of ECoG connectivity. Our findings provide fundamental understanding of the factors that influence interregional ECoG networks, and open the possibility for predictive modeling of surgical outcomes in disease.


title: Bayesian multi--dipole localization and uncertainty quantification from simultaneous EEG and MEG recordings

abstract: We deal with estimation of multiple dipoles from combined MEG and EEG time--series. We use a sequential Monte Carlo algorithm to characterize the posterior distribution of the number of dipoles and their locations. By considering three test cases, we show that using the combined data the method can localize sources that are not easily (or not at all) visible with either of the two individual data alone. In addition, the posterior distribution from combined data exhibits a lower variance, i.e. lower uncertainty, than the posterior from single device.


title: On the Bhattacharya-Mesner rank of third order hypermatrices

abstract: We introduce the Bhattacharya-Mesner rank of third order hypermatrices as a relaxation to the tensor rank and devise from it some bounds for the tensor rank. We use the Bhattacharya-Mesner rank to extend to third order hypermatrices the connection relating the rank to a notion of linear dependence. We also derive explicit necessary and sufficient conditions for the existence of third order hypermatrix inverse pair. Finally we use inverse pair to extend to third order hypermatrices the formulation and proof of the matrix rank-nullity theorem.


title: Counting Markov Equivalence Classes for DAG models on Trees

abstract: DAG models are statistical models satisfying a collection of conditional independence relations encoded by the nonedges of a directed acyclic graph (DAG) $\mathcal{G}$. Such models are used to model complex cause-effect systems across a variety of research fields. From observational data alone, a DAG model $\mathcal{G}$ is only recoverable up to Markov equivalence. Combinatorially, two DAGs are Markov equivalent if and only if they have the same underlying undirected graph (i.e. skeleton) and the same set of the induced subDAGs $i\to j \leftarrow k$, known as immoralities. Hence it is of interest to study the number and size of Markov equivalence classes (MECs). In a recent paper, the authors introduced a pair of generating functions that enumerate the number of MECs on a fixed skeleton by number of immoralities and by class size, and they studied the complexity of computing these functions. In this paper, we lay the foundation for studying these generating functions by analyzing their structure for trees and other closely related graphs. We describe these polynomials for some important families of graphs including paths, stars, cycles, spider graphs, caterpillars, and complete binary trees. In doing so, we recover important connections to independence polynomials, and extend some classical identities that hold for Fibonacci numbers. We also provide tight lower and upper bounds for the number and size of MECs on any tree. Finally, we use computational methods to show that the number and distribution of high degree nodes in a triangle-free graph dictates the number and size of MECs.


title: Evidence of an Internal Dissipation Origin for the High-energy Prompt Emission of GRB 170214A

abstract: The origin of the prompt high-energy ($>100$MeV) emission of Gamma-ray Bursts (GRBs), detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, is still under debate, for which both the external shock origin and internal dissipation origin have been suggested. In the internal dissipation scenario, the high energy emission is expected to exhibit significant temporal variability, tracking the keV/MeV fast variable behavior. Here, we report a detailed analysis on the Fermi data of GRB~170214A, which is sufficiently bright in the high energy to enable a quantitative analysis of the correlation between high-energy emission and keV/MeV emission with high statistics. Our result shows a clear temporal correlation between high-energy and keV/MeV emission in the whole prompt emission phase as well as in two decomposed short time intervals. Such correlation behavior is also found in some other bright LAT GRBs, i.e., GRB 080916C, 090902B and 090926A. For these GRBs as well as GRB 090510, we also find the rapid temporal variability in the high-energy emission. We thus conclude that the prompt high-energy emission in these bright LAT GRBs should be due to internal origin.


title: Extraterrestrial sedimentary rocks on Earth

abstract: This concept article discusses the possibilities for identifying sedimentary-origin meteorites. The paper concerns (i) the macroscopic candidate for sedimentary meteorite in the epicenter of the 1908 Tunguska catastrophe; (ii) potential parent bodies for sedimentary meteorites; (iii) isotopic heterogeneity of unmixed silicate reservoirs on Mars; (iv) possible terrestrial loss or contamination in the noble gas signatures in new type meteorites that spent time in extreme weather conditions; (v) cosmogenic isotopes and shielding; and (vi) pseudo meteorites. We conclude that the list of candidate parent bodies for sedimentary meteorites includes, but is not limited by the Earth, Mars, Enceladus, Ganymede, Europa, and hypothetical planets that could exist between orbits of Mars and Jupiter in the past. A parent body for extraterrestrial sedimentary rocks on the Earth should be identified based on the entire body of evidence which is not limited solely by tests of oxygen and noble gas isotopes whose signatures may undergo terrestrial contamination and may exhibit significant heterogeneity within the parent bodies. Observed fall of cosmic body, evidence of hypervelocity impact complying with the criteria of impact structures, and the presence of fusion crust on the fragments should be considered as priority signs of meteoritic origin.


title: Real-time digital signal recovery for a multi-pole low-pass transfer function system

abstract: In order to solve the problems of waveform distortion and signal delay by many physical and electrical systems with multi-pole linear low-pass transfer characteristics, a simple digital-signal-processing (DSP)-based method of real-time recovery of the original source waveform from the distorted output waveform is proposed. A mathematical analysis on the convolution kernel representation of the single-pole low-pass transfer function shows that the original source waveform can be accurately recovered in real time using a particular moving average algorithm applied on the input stream of the distorted waveform, which can also significantly reduce the overall delay time constant. This method is generalized for multi-pole low-pass systems and has noise characteristics of the inverse of the low-pass filter characteristics. This method can be applied to most sensors and amplifiers operating close to their frequency response limits to improve the overall performances of data acquisition systems and digital feedback control systems.


title: Block partitions: an extended view

abstract: Given a sequence $S=(s_1,\dots,s_m) \in [0, 1]^m$, a block $B$ of $S$ is a subsequence $B=(s_i,s_{i+1},\dots,s_j)$. The size $b$ of a block $B$ is the sum of its elements. It is proved in [1] that for each positive integer $n$, there is a partition of $S$ into $n$ blocks $B_1, \dots , B_n$ with $|b_i - b_j| \le 1$ for every $i, j$. In this paper, we consider a generalization of the problem in higher dimensions.


title: Inferring directed climatic interactions with renormalized partial directed coherence and directed partial correlation

abstract: Inferring interactions between processes promises deeper insight into mechanisms underlying network phenomena. Renormalised partial directed coherence (rPDC) is a frequency-domain representation of the concept of Granger causality while directed partial correlation (DPC) is an alternative approach for quantifying Granger causality in the time domain. Both methodologies have been successfully applied to neurophysiological signals for detecting directed relationships. This paper introduces their application to climatological time series. We first discuss the application to ENSO -- Monsoon interaction, and then apply the methodologies to the more challenging air-sea interaction in the South Atlantic Convergence Zone (SACZ). While in the first case the results obtained are fully consistent with present knowledge in climate modeling, in the second case the results are, as expected, less clear, and to fully elucidate the SACZ air-sea interaction, further investigations on the specificity and sensitivity of these methodologies are needed.


title: Control of accuracy in the Wang-Landau algorithm

abstract: The Wang-Landau (WL) algorithm has been widely used for simulations in many areas of physics. Our analysis of the WL algorithm explains its properties and shows that the difference of the largest eigenvalue of the transition matrix in the energy space from unity can be used to control the accuracy of estimating the density of states. Analytic expressions for the matrix elements are given in the case of the one-dimensional Ising model. The proposed method is further confirmed by numerical results for the one-dimensional and two-dimensional Ising models and also the two-dimensional Potts model.


title: How SU(2)$_4$ Anyons are Z$_3$ Parafermions

abstract: We consider the braid group representation which describes the non-abelian braiding statistics of the spin 1/2 particle world lines of an SU(2)$_4$ Chern-Simons theory. Up to an abelian phase, this is the same as the non-Abelian statistics of the elementary quasiparticles of the $k=4$ Read-Rezayi quantum Hall state. We show that these braiding statistics are identical to those of Z$_3$ Parafermions.


title: Multiwavelength characterisation of an ACT-selected, lensed dusty star-forming galaxy at z=2.64

abstract: We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow.


title: Flavourful $Z'$ models for $R_{K^{(*)}}$

abstract: We show how any flavour conserving $Z'$ model can be made flavour violating and non-universal by introducing mass mixing of quarks and leptons with a fourth family of vector-like fermions with non-universal $Z'$ couplings. After developing a general formalism, we focus on two concrete examples, namely a fermiophobic model, and an $SO(10)$ GUT model, and show how they can account for the anomalous $B$ decay ratios $R_K$ and $R_{K^*}$. A similar analysis could be performed for $B-L$ models, $E_6$ models, composite models, and so on.


title: Collapse and Nonlinear Instability of AdS with Angular Momenta

abstract: We present a numerical study of rotational dynamics in AdS$_5$ with equal angular momenta in the presence of a complex doublet scalar field. We determine that the endpoint of gravitational collapse is a Myers-Perry black hole for high energies and a hairy black hole for low energies. We investigate the timescale for collapse at low energies $E$, keeping the angular momenta $J\propto E$ in AdS length units. We find that the inclusion of angular momenta delays the collapse time, but retains a $t\sim1/E$ scaling. We perturb and evolve rotating boson stars, and find that boson stars near AdS appear stable, but those sufficiently far from AdS are unstable. We find that the dynamics of the boson star instability depend on the perturbation, resulting either in collapse to a Myers-Perry black hole, or development towards a stable oscillating solution.


title: Formation and Incidence of Shell Galaxies in the Illustris Simulation

abstract: Shells are low surface brightness tidal debris that appear as interleaved caustics with large opening angles, often situated on both sides of the galaxy center. In this paper, we study the incidence and formation processes of shell galaxies in the cosmological gravity+hydrodynamics Illustris simulation. We identify shells at redshift z=0 using stellar surface density maps, and we use stellar history catalogs to trace the birth, trajectory and progenitors of each individual star particle contributing to the tidal feature. Out of a sample of the 220 most massive galaxies in Illustris ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_{\odot}$), $18\%\pm3\%$ of the galaxies exhibit shells. This fraction increases with increasing mass cut: higher mass galaxies are more likely to have stellar shells. Furthermore, the fraction of massive galaxies that exhibit shells decreases with increasing redshift. We find that shell galaxies observed at redshift $z=0$ form preferentially through relatively major mergers ($\gtrsim$1:10 in stellar mass ratio). Progenitors are accreted on low angular momentum orbits, in a preferred time-window between $\sim$4 and 8 Gyrs ago. Our study indicates that, due to dynamical friction, more massive satellites are allowed to probe a wider range of impact parameters at accretion time, while small companions need almost purely radial infall trajectories in order to produce shells. We also find a number of special cases, as a consequence of the additional complexity introduced by the cosmological setting. These include galaxies with multiple shell-forming progenitors, satellite-of-satellites also forming shells, or satellites that fail to produce shells due to multiple major mergers happening in quick succession.


title: Probing motion of fast radio burst sources by timing strongly lensed repeaters

abstract: Given the possible repetitive nature of fast radio bursts (FRBs), their cosmological origin, and their high occurrence, detection of strongly lensed sources due to intervening galaxy lenses is possible with forthcoming radio surveys. We show that if multiple images of a repeating source are resolved with VLBI, using a method independent of lens modeling, accurate timing could reveal non-uniform motion, either physical or apparent, of the emission spot. This can probe the physical nature of FRBs and their surrounding environments, constraining scenarios including orbital motion around a stellar companion if FRBs require a compact star in a special system, and jet-medium interactions for which the location of the emission spot may randomly vary. The high timing precision possible for FRBs ($\sim {\rm ms}$) compared to the typical time delays between images in galaxy lensing ($\gtrsim 10\, {\rm days}$) enables the measurement of tiny fractional changes in the delays ($\sim 10^{-9}$), and hence the detection of time-delay variations induced by relative motions between the source, the lens, and the Earth. We show that uniform cosmic peculiar velocities only cause the delay time to drift linearly, and that the effect from the Earth's orbital motion can be accurately subtracted, thus enabling a search for non-trivial source motion. For a timing accuracy of $\sim 1\,$ms and a repetition rate (of detected bursts) $\sim 0.05$ per day of a single FRB source, non-uniform displacement $\gtrsim 0.1 - 1\,$AU of the emission spot perpendicular to the line of sight is detectable if repetitions are seen over a period of hundreds of days.


title: Anomaly constraints on deconfinement and chiral phase transition

abstract: We study constraints on thermal phase transitions of ${\rm SU}(N_c)$ gauge theories by using the 't Hooft anomaly involving the center symmetry and chiral symmetry. We consider two cases of massless fermions: (i) adjoint fermions, and (ii) $N_f$ flavors of fundamental fermions with a nontrivial greatest common divisor ${\rm gcd}(N_c,N_f) \neq 1$. For the first case (i), we show that the chiral symmetry restoration in terms of the standard Landau-Ginzburg effective action is impossible at a temperature lower than that of deconfinement. For the second case (ii), we introduce a modified version of the center symmetry which we call center-flavor symmetry, and draw similar conclusions under a certain definition of confinement. Moreover, at zero temperature, our results give a partial explanation of the appearance of dual magnetic gauge group in (supersymmetric) QCD when ${\rm gcd}(N_c,N_f) \neq 1$.


title: Angular Momentum Evolution of Young Stars in the Nearby Scorpius-Centaurus OB Association

abstract: We report the results of a study of archival SuperWASP light curves for stars in Scorpius-Centaurus (Sco-Cen), the nearest OB association. We use SuperWASP time-series photometry to extract rotation periods for 189 candidate members of the Sco-Cen complex, and verify that 162 of those are members of the classic Sco-Cen subgroups of Upper Scorpius (US), Upper Centaurus-Lupus (UCL), and Lower Centaurus-Crux (LCC). This study provides the first measurements of rotation periods explicitly for large samples of pre-main sequence (pre-MS) stars spanning the UCL and LCC subgroups. Our final sample of 157 well-characterized pre-MS stars spans ages of $\sim$ 10 -- 20 Myr, spectral types of $\sim$ F3 -- M0, and masses of $M \simeq$ 0.3 -- 1.5 $\mathcal{M}^{\rm N}_{\odot}$. For this sample, we find a distribution of stellar rotation periods with a median of $P_{\textrm{rot}} \simeq$ 2.4 days, overall range of 0.2 $<$ $P_{\textrm{rot}}$ $<$ 8 days, and a fairly well-defined mass-dependent upper envelope of rotation periods. This distribution of periods is consistent with recently developed stellar angular momentum evolution models. These data are significant because they represent an undersampled age range and the number of measurable rotation periods is large compared to recent studies of other regions. We also search for new examples of eclipsing disk or ring systems analogous to 1SWASP J140747.93-394542.6 ("J1407", V1400 Cen), but find none. Our survey yielded five eclipsing binaries, but only one appears to be physically associated with the Sco-Cen complex. V2394 Oph is a heavily reddened (A$_V$ $\simeq$ 5 mag) massive contact binary in the LDN 1689 cloud whose Gaia astrometry is clearly consistent with kinematic membership with the Ophiuchus star-forming region.


title: Capability of Detecting Ultra-Violet Counterparts of Gravitational Waves with GLUV

abstract: With the discovery of gravitational waves (GW), attention has turned towards detecting counterparts to these sources. In discussions on counterpart signatures and multi-messenger follow-up strategies to GW detections, ultra-violet (UV) signatures have largely been neglected, due to UV facilities being limited to SWIFT, which lacks high-cadence UV survey capabilities. In this paper, we examine the UV signatures from merger models for the major GW sources, highlighting the need for further modelling, while presenting requirements and a design for an effective UV survey telescope. Using $u'$-band models as an analogue, we find that a UV survey telescope requires a limiting magnitude of m$_{u'}\rm (AB)\approx 24$ to fully complement the aLIGO range and sky localisation. We show that a network of small, balloon-based UV telescopes with a primary mirror diameter of 30~cm could be capable of covering the aLIGO detection distance from $\sim$60--100\% for BNS events and $\sim$40\% for BHNS events. The sensitivity of UV emission to initial conditions suggests that a UV survey telescope would provide a unique dataset, that can act as an effective diagnostic to discriminate between models.


title: Pixel Color-Magnitude Diagram Analysis of the Brightest Cluster Galaxies in Dynamically Young and Old Clusters, Abell 1139 and Abell 2589

abstract: As a case study to understand the coevolution of Brightest Cluster Galaxies (BCGs) and their host clusters, we investigate the BCGs in dynamically young and old clusters, Abell 1139 (A1139) and Abell 2589 (A2589). We analyze the pixel color-magnitude diagrams (pCMDs) using deep g- and r-band images, obtained from the Canada-France-Hawaii Telescope observations. After masking foreground/background objects and smoothing pixels in consideration of the observational seeing size, detailed pCMD features are compared between the two BCGs. (1) While the overall shapes of the pCMDs are similar to those of typical early-type galaxies, the A2589-BCG tends to have redder mean pixel color and smaller pixel color deviation at given surface brightness than the A1139-BCG. (2) The mean pixel color distribution as a function of pixel surface brightness (pCMD backbone) indicates that the A2589-BCG formed a larger central body (~ 2.0 kpc in radius) by major dry mergers at an early epoch than the A1139-BCG (a central body ~ 1.3 kpc in radius), while they have grown commonly by subsequent minor mergers. (3) The spatial distributions of the pCMD outliers reveal that the A1139-BCG experienced considerable tidal events more recently than the A2589-BCG, whereas the A2589-BCG has an asymmetric compact core possibly resulting from major dry merger at an early epoch. (4) The A2589-BCG shows a very large faint-to-bright pixel number ratio compared to early-type non-BCGs, whereas the ratio for the A1139-BCG is not distinctively large. These results are consistent with the idea that the BCG in the dynamically older cluster (A2589) formed earlier and is relaxed better.


title: Microflare Heating of a Solar Active Region Observed with NuSTAR, Hinode/XRT, and SDO/AIA

abstract: NuSTAR is a highly sensitive focusing hard X-ray (HXR) telescope and has observed several small microflares in its initial solar pointings. In this paper, we present the first joint observation of a microflare with NuSTAR and Hinode/XRT on 2015 April 29 at ~11:29 UT. This microflare shows heating of material to several million Kelvin, observed in Soft X-rays (SXRs) with Hinode/XRT, and was faintly visible in Extreme Ultraviolet (EUV) with SDO/AIA. For three of the four NuSTAR observations of this region (pre-, decay, and post phases) the spectrum is well fitted by a single thermal model of 3.2-3.5 MK, but the spectrum during the impulsive phase shows additional emission up to 10 MK, emission equivalent to A0.1 GOES class. We recover the differential emission measure (DEM) using SDO/AIA, Hinode/XRT, and NuSTAR, giving unprecedented coverage in temperature. We find the pre-flare DEM peaks at ~3 MK and falls off sharply by 5 MK; but during the microflare's impulsive phase the emission above 3 MK is brighter and extends to 10 MK, giving a heating rate of about $2.5 \times 10^{25}$ erg s$^{-1}$. As the NuSTAR spectrum is purely thermal we determined upper-limits on the possible non-thermal bremsstrahlung emission. We find that for the accelerated electrons to be the source of the heating requires a power-law spectrum of $\delta \ge 7$ with a low energy cut-off $E_{c} \lesssim 7$ keV. In summary, this first NuSTAR microflare strongly resembles much more powerful flares.


title: Very massive stars, pair-instability supernovae and intermediate-mass black holes with the SEVN code

abstract: Understanding the link between massive ($\gtrsim 30$ M$_{\odot{}}$) stellar black holes (BHs) and their progenitor stars is a crucial step to interpret observations of gravitational-wave events. In this paper, we discuss the final fate of very massive stars (VMSs), with zero-age main sequence (ZAMS) mass $>150$ M$_{\odot}$, accounting for pulsational pair-instability supernovae (PPISNe) and for pair-instability supernovae (PISNe). We describe an updated version of our population synthesis code SEVN, in which we added stellar evolution tracks for VMSs with ZAMS mass up to $350 $M$_{\odot{}}$ and we included analytical prescriptions for PPISNe and PISNe. We use the new version of SEVN to study the BH mass spectrum at different metallicity $Z$, ranging from $Z=2.0\times 10^{-4}$ to $Z=2.0\times 10^{-2}$. The main effect of PPISNe and PISNe is to favour the formation of BHs in the mass range of the first gravitational-wave event (GW150914), while they prevent the formation of remnants with mass 60 - 120 M$_{\odot{}}$. In particular, we find that PPISNe significantly enhance mass loss of metal-poor ($Z\leq 2.0\times 10^{-3}$) stars with ZAMS mass $60\leq $M$_{\mathrm{ZAMS}}/$M$_{\odot{}}\leq 125$. In contrast, PISNe become effective only for moderately metal-poor ($Z<8.0\times 10^{-3}$) VMSs. VMSs with M$_{\rm ZAMS}\gtrsim{}220$ M$_\odot$ and $Z<10^{-3}$ do not undergo PISNe and form intermediate-mass BHs (IMBHs, with mass $\gtrsim 200 $M$_{\odot{}}$) via direct collapse.


title: Detection of faint broad emission lines in type 2 AGN: III. On the $M_{BH} - \sigma_\star$ relation of type 2 AGN

abstract: Type 2 active galactic nuclei (AGN) represent the majority of the AGN population. However, due to the difficulties in measuring their black hole (BH) masses, it is still unknown whether they follow the same BH mass-host galaxy scaling relations valid for quiescent galaxies and type 1 AGN. Here we present the locus of type 2 AGN having virial BH mass estimates in the $M_{BH}-\sigma_\star$ plane. Our analysis shows that the BH masses of type 2 AGN are $\sim0.9$ dex smaller than type 1 AGN at $\sigma_\star\sim 185$ km s$^{-1}$, regardless of the (early/late) AGN host galaxy morphology. Equivalently, type 2 AGN host galaxies have stellar velocity dispersions $\sim 0.2$ dex higher than type 1 AGN hosts at $M_{BH}\sim10^7$ M$_\odot$.


title: Probing the small-scale structure in strongly lensed systems via transdimensional inference

abstract: Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.


title: Double arc instability in the solar corona

abstract: The stability of the magnetic field in the solar corona is important for understanding the causes of solar eruptions. Although various scenarios have been suggested to date, the tether-cutting reconnection scenario proposed by Moore et al.(2001) is one of the widely accepted models to explain the onset process of solar eruptions. Although the tether-cutting reconnection scenario proposed that sigmoidal field formed by the internal reconnection is the magnetic field in pre-eruptive state, the stability of the sigmoidal field has not yet been investigated quantitatively. In this paper, in order to elucidate the stability problem of pre-eruptive state, we developed a simple numerical analysis, in which the sigmoidal field is modeled by a double arc electric current loop and its stability is analyzed. As a result, we found that the double arc loop is more easily destabilized than the axisymmetric torus, and it becomes unstable even if the external field does not decay with altitude, which is in contrast to the axisymmetric torus instability. This suggests that the tether-cutting reconnection may well work as the onset mechanism of solar eruptions, and if so the critical condition for eruption under certain geometry may be determined by a new type of instability rather than the torus instability. Based on them, we propose a new type of instability called double arc instability (DAI). We discuss the critical conditions for DAI and derive a new parameter $\kappa$ defined as the product of the magnetic twist and the normalized flux of tether-cutting reconnection.


title: Models and Materials for Generalized Kitaev Magnetism

abstract: The exactly solvable Kitaev model on the honeycomb lattice has recently received enormous attention linked to the hope of achieving novel spin-liquid states with fractionalized Majorana-like excitations. In this review, we analyze the mechanism proposed by G. Jackeli and G. Khaliullin to identify Kitaev materials based on spin-orbital dependent bond interactions and provide a comprehensive overview of its implications in real materials. We set the focus on experimental results and current theoretical understanding of planar honeycomb systems (Na$_2$IrO$_3$, $\alpha$-Li$_2$IrO$_3$, and $\alpha$-RuCl$_3$), three-dimensional Kitaev materials ($\beta$- and $\gamma$-Li$_2$IrO$_3$), and other potential candidates, completing the review with the list of open questions awaiting new insights.


title: Laser-irradiated Kondo insulators : Controlling Kondo effect and topological phases

abstract: We theoretically investigate the nature of laser-irradiated Kondo insulators. Using Floquet theory and slave boson approach, we study a periodic Anderson model and derive an effective model which describes the laser-irradiated Kondo insulators. In this model, we find two generic effects induced by laser light. One is the dynamical localization, which suppresses hopping and hybridization. The other is the laser-induced hopping and hybridization, which can be interpreted as a synthetic spin-orbit coupling or magnetic field. The first effect drastically changes the behavior of the Kondo effect. Especially, the Kondo effect under laser light qualitatively changes its character depending on whether the hybridization is on-site or off-site. The second effect triggers topological phase transitions. In topological Kondo insulators, linearly polarized laser light realizes phase transitions between trivial, weak topological, and strong topological Kondo insulators. Moreover, circularly polarized laser light breaks time-reversal symmetry and induces Weyl semimetallic phases. Our results pave the new way to dynamically control the Kondo effect and topological phases in heavy fermion systems. We also discuss experimental setups to detect the signatures.


title: Correlation Energy of Proton-Neutron Subsystem in Valence Orbit

abstract: Deuteron correlation energy (DCE) of the valence proton-neutron subsystem is evaluated by utilizing a simple three-body model. We focus on the $^6$Li and $^{18}$F nuclei assuming the doubly-closed core and the valence proton and neutron. Two interaction models, schematic density-dependent contact (SDDC) and Minnesota potentials, are utilized to describe the proton-neutron interaction. Evaluating DCE, we conclude that the proton-neutron binding in $^6$Li can be stronger than its counterpart of a deuteron in vacuum. On the other hand, in $^{18}$F, the energetic correlation is remarkably weak, and does not favor the bound, deuteron-like configuration. This significant difference between two systems can be understood from a competition between the proton-neutron kinetic and pairing energies, which are sensitive to the spatial extension of the wave function. This result indicates a remarkable dependence of the deuteron correlation to its environment and the valence orbits.


title: Observation of fractional Chern insulators in a van der Waals heterostructure

abstract: Topologically ordered phases are characterized by long-range quantum entanglement and fractional statistics rather than by symmetry breaking. First observed in a fractionally filled continuum Landau level, topological order has since been proposed to arise more generally at fractional filling of topologically non-trivial "Chern" bands. Here, we report the observation of gapped states at fractional filling of Harper-Hofstadter bands arising from the interplay of a magnetic field and a superlattice potential in a bilayer graphene/hexagonal boron nitride heterostructure. We observe new phases at fractional filling of bands with Chern indices $\mathcal{C} = -1, \pm 2,$ and $\pm 3$. Some of these, in $\mathcal{C}=-1$ and $\mathcal{C}=2$ bands, are characterized by fractional Hall conductance---they are `fractional Chern insulators' and constitute a new example of topological order beyond Landau levels.


title: Large, nonsaturating thermopower in a quantizing magnetic field

abstract: The thermoelectric effect is the generation of an electrical voltage from a temperature gradient in a solid material due to the diffusion of free charge carriers from hot to cold. Identifying materials with large thermoelectric response is crucial for the development of novel electric generators and coolers. In this paper we consider theoretically the thermopower of Dirac/Weyl semimetals subjected to a quantizing magnetic field. We contrast their thermoelectric properties with those of traditional heavily-doped semiconductors and we show that, under a sufficiently large magnetic field, the thermopower of Dirac/Weyl semimetals grows linearly with the field without saturation and can reach extremely high values. Our results suggest an immediate pathway for achieving record-high thermopower and thermoelectric figure of merit, and they compare well with a recent experiment on Pb$_{1-x}$Sn$_x$Se.


title: Three-body Unitarity with Isobars Revisited

abstract: The particle exchange model of hadron interactions can be used to describe three-body scattering under the isobar assumption. In this study we start from the 3->3 scattering amplitude for spinless particles, which contains an isobar-spectator scattering amplitude. Using a Bethe-Salpeter Ansatz for the latter, we derive a relativistic three-dimensional scattering equation that manifestly fulfills three-body unitarity and two-body unitarity for the sub-amplitudes. This property holds for energies above breakup and also in the presence of resonances in the sub-amplitudes.


title: The spatially-resolved star formation histories of CALIFA galaxies: Implications for galaxy formation

abstract: This paper presents the spatially resolved star formation history (SFH) of nearby galaxies with the aim of furthering our understanding of the different processes involved in the formation and evolution of galaxies. To this end, we apply the fossil record method of stellar population synthesis to a rich and diverse data set of 436 galaxies observed with integral field spectroscopy in the CALIFA survey. The sample covers a wide range of Hubble types, with stellar masses ranging from $M_\star \sim 10^9$ to $7 \times 10^{11} M_\odot$. Spectral synthesis techniques are applied to the datacubes to retrieve the spatially resolved time evolution of the star formation rate (SFR), its intensity ($\Sigma_{\rm SFR}$), and other descriptors of the 2D-SFH in seven bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc, and Sd), and five bins of stellar mass. Our main results are: a) Galaxies form very fast independently of their current stellar mass, with the peak of star formation at high redshift ($z > 2$). Subsequent star formation is driven by $M_\star$ and morphology, with less massive and later type spirals showing more prolonged periods of star formation. b) At any epoch in the past the SFR is proportional to $M_\star$, with most massive galaxies having the highest absolute (but lowest specific) SFRs. c) While nowadays $\Sigma_{\rm SFR}$ is similar for all spirals, and significantly lower in early type galaxies (ETG), in the past $\Sigma_{\rm SFR}$ scales well with morphology. The central regions of today's ETGs are where $\Sigma_{\rm SFR}$ reached the highest values ($> 10^3 \,M_\odot\,$Gyr$^{-1}\,$pc$^{-2}$), similar to those measured in high redshift star forming galaxies. d) The evolution of $\Sigma_{\rm SFR}$ in Sbc systems matches that of models for Milky-Way-like galaxies, suggesting that the formation of a thick disk may be a common phase in spirals at early epochs.


title: Multi-Label Annotation Aggregation in Crowdsourcing

abstract: As a means of human-based computation, crowdsourcing has been widely used to annotate large-scale unlabeled datasets. One of the obvious challenges is how to aggregate these possibly noisy labels provided by a set of heterogeneous annotators. Another challenge stems from the difficulty in evaluating the annotator reliability without even knowing the ground truth, which can be used to build incentive mechanisms in crowdsourcing platforms. When each instance is associated with many possible labels simultaneously, the problem becomes even harder because of its combinatorial nature. In this paper, we present new flexible Bayesian models and efficient inference algorithms for multi-label annotation aggregation by taking both annotator reliability and label dependency into account. Extensive experiments on real-world datasets confirm that the proposed methods outperform other competitive alternatives, and the model can recover the type of the annotators with high accuracy.


title: High-Precision Calculations in Strongly Coupled Quantum Field Theory with Next-to-Leading-Order Renormalized Hamiltonian Truncation

abstract: Hamiltonian Truncation (a.k.a. Truncated Spectrum Approach) is an efficient numerical technique to solve strongly coupled QFTs in d=2 spacetime dimensions. Further theoretical developments are needed to increase its accuracy and the range of applicability. With this goal in mind, here we present a new variant of Hamiltonian Truncation which exhibits smaller dependence on the UV cutoff than other existing implementations, and yields more accurate spectra. The key idea for achieving this consists in integrating out exactly a certain class of high energy states, which corresponds to performing renormalization at the cubic order in the interaction strength. We test the new method on the strongly coupled two-dimensional quartic scalar theory. Our work will also be useful for the future goal of extending Hamiltonian Truncation to higher dimensions d >= 3.


title: VAIN: Attentional Multi-agent Predictive Modeling

abstract: Multi-agent predictive modeling is an essential step for understanding physical, social and team-play systems. Recently, Interaction Networks (INs) were proposed for the task of modeling multi-agent physical systems, INs scale with the number of interactions in the system (typically quadratic or higher order in the number of agents). In this paper we introduce VAIN, a novel attentional architecture for multi-agent predictive modeling that scales linearly with the number of agents. We show that VAIN is effective for multi-agent predictive modeling. Our method is evaluated on tasks from challenging multi-agent prediction domains: chess and soccer, and outperforms competing multi-agent approaches.


title: Absorption systems at $z\sim 2$ as a probe of the circum galactic medium: a probabilistic approach

abstract: We characterize the properties of the intergalactic medium (IGM) around a sample of galaxies extracted from state-of-the-art hydrodynamical simulations of structure formation in a cosmological volume of 25 Mpc comoving at $z\sim 2$. The simulations are based on two different subresolution schemes for star formation and supernova feedback: the MUlti-Phase Particle Integrator (MUPPI) scheme and the Effective Model. We develop a quantitative and probabilistic analysis based on the apparent optical depth method of the properties of the absorbers as a function of impact parameter from their nearby galaxies: in such a way we probe different environments from circumgalactic medium to low-density filaments. Absorbers' properties are then compared with a spectroscopic observational data set obtained from high-resolution quasar spectra. Our main focus is on the N$_{\rm{ CIV}}$-N$_{\rm { HI}}$ relation around simulated galaxies: the results obtained with MUPPI and the Effective model are remarkably similar, with small differences only confined to regions at impact parameters $b = [1-3] \times r_{\rm {vir}}$. Using $\mbox{C IV}$ as a tracer of the metallicity, we obtain evidence that the observed metal absorption systems have the highest probability to be confined in a region of 150-400 kpc around galaxies. Near-filament environments have instead metallicities too low to be probed by present-day telescopes, but could be probed by future spectroscopical studies. Finally we compute $\mbox{C IV}$ covering fractions which are in agreement with observational data.


title: Ground-based photometry of the 21-day Neptune HD106315c

abstract: Space-based transit surveys such as K2 and TESS allow the detection of small transiting planets with orbital periods beyond 10 days. Few of these warm Neptunes are currently known around stars bright enough to allow for detailed follow-up observations dedicated to their atmospheric characterization. The 21-day period and 3.95 $R_\oplus$ planet HD106315c has been discovered based on the observation of two of its transits by K2. We have observed HD106315 using the 1.2m Euler telescope equipped with the EulerCam camera on two instances to confirm the transit using broad band photometry and refine the planetary period. Based on two observed transits of HD106315c, we detect its $\sim$1 mmag transit and obtain a precise measurement of the planetary ephemerids, which are critical for planning further follow-up observations. We have used the attained precision together with the predicted yield from the TESS mission to evaluate the potential for ground-based confirmation of Neptune-sized planets found by TESS. We find that 1-meter-class telescopes on the ground equipped with precise photometers could substantially contribute to the follow-up of 162 TESS candidates orbiting stars with magnitudes of $V \leq 14$. Out of these, 74 planets orbit stars with $V \leq 12$ and 12 planets orbit $V \leq 10$, which makes these candidates high-priority objects for atmospheric characterization with high-end instrumentation.


title: The galaxy-subhalo connection in low-redshift galaxy clusters from weak gravitational lensing

abstract: We measure the gravitational lensing signal around satellite galaxies in a sample of galaxy clusters at $z<0.15$ by combining high-quality imaging data from the Canada-France-Hawaii Telescope with a large sample of spectroscopically-confirmed cluster members. We use extensive image simulations to assess the accuracy of shape measurements of faint, background sources in the vicinity of bright satellite galaxies. We find a small but significant bias, as light from the lenses makes the shapes of background galaxies appear radially aligned with the lens. We account for this bias by applying a correction that depends on both lens size and magnitude. We also correct for contamination of the source sample by cluster members. We use a physically-motivated definition of subhalo mass, namely the mass bound to the subhalo, $m_\mathrm{bg}$, similar to definitions used by common subhalo finders in numerical simulations. Binning the satellites by stellar mass we provide a direct measurement of the subhalo-to-stellar-mass relation, $\log m_\mathrm{bg}/\mathrm{M}_\odot = (11.54\pm0.05) + (0.95\pm0.10)\log[m_\star/(2\times10^{10}\mathrm{M}_\odot)]$. This best-fitting relation implies that, at a stellar mass $m_\star\sim3\times10^{10}\,\mathrm{M}_\odot$, subhalo masses are roughly 50 per cent of those of central galaxies, and this fraction decreases at higher stellar masses. We find some evidence for a sharp change in the total-to-stellar mass ratio around the clusters' scale radius, which could be interpreted as galaxies within the scale radius having suffered more strongly from tidal stripping, but remain cautious regarding this interpretation.


title: Analytic approximation of solutions of parabolic partial differential equations with variable coefficients

abstract: A complete family of solutions for the one-dimensional reaction-diffusion equation \[ u_{xx}(x,t)-q(x)u(x,t) = u_t(x,t) \] with a coefficient $q$ depending on $x$ is constructed. The solutions represent the images of the heat polynomials under the action of a transmutation operator. Their use allows one to obtain an explicit solution of the noncharacteristic Cauchy problem for the considered equation with sufficiently regular Cauchy data as well as to solve numerically initial boundary value problems. In the paper the Dirichlet boundary conditions are considered however the proposed method can be easily extended onto other standard boundary conditions. The proposed numerical method is shown to reveal good accuracy.


title: Multiple Chirality in Nuclear Rotation: A Microscopic View

abstract: Covariant density functional theory and three-dimensional tilted axis cranking are used to investigate multiple chirality in nuclear rotation for the first time in a fully self-consistent and microscopic way. Two distinct sets of chiral solutions with negative and positive parities, respectively, are found in the nucleus 106Rh. The negative-parity solutions reproduce well the corresponding experimental spectrum as well as the B(M1)/B(E2) ratios of the transition strengths. This indicates that a predicted positive-parity chiral band should also exist. Therefore, it provides a further strong hint that multiple chirality is realized in nuclei.


title: Global Orientifolded Quivers with Inflation

abstract: We describe global embeddings of fractional D3 branes at orientifolded singularities in type IIB flux compactifications. We present an explicit Calabi-Yau example where the chiral visible sector lives on a local orientifolded quiver while non-perturbative effects, $\alpha'$ corrections and a T-brane hidden sector lead to full closed string moduli stabilisation in a de Sitter vacuum. The same model can also successfully give rise to inflation driven by a del Pezzo divisor. Our model represents the first explicit Calabi-Yau example featuring both an inflationary and a chiral visible sector.


title: Multiplexing 200 modes on a single digital hologram

abstract: The on-demand tailoring of light's spatial shape is of great relevance in a wide variety of research areas. Computer-controlled devices, such as Spatial Light Modulators (SLMs) or Digital Micromirror Devices (DMDs), offer a very accurate, flexible and fast holographic means to this end. Remarkably, digital holography affords the simultaneous generation of multiple beams (multiplexing), a tool with numerous applications in many fields. Here, we provide a self-contained tutorial on light beam multiplexing. Through the use of several examples, the readers will be guided step by step in the process of light beam shaping and multiplexing. Additionally, on the multiplexing capabilities of SLMs to provide a quantitative analysis on the maximum number of beams that can be multiplexed on a single SLM, showing approximately 200 modes on a single hologram.


title: Translation-invariant Gibbs measures for the Blum-Kapel model on a Cayley tree

abstract: In this paper we consider translation-invariant Gibbs measures for the Blum-Kapel model on a Cayley tree of order k. An approximate critical temperature T_{cr} is found such that for T\geq T_{cr} there exists a unique translation-invariant Gibbs measure and for 0<T<T_{cr} there are exactly three translation-invariant Gibbs measures. In addition, we studied the problem of (not) extremality for the unique Gibbs measure.


title: A bug's eye view: the Riemannian exponential map on polyhedral surfaces

abstract: We explore the perspective of a bug living on the two-dimensional surface of a polyhedron. Images of various kinds of effects like lensing and cloaking are shown via color pictures of three viewpoints: the first person perspective of the bug, a map of the bug's viewpoint, and a look at the bug on the embedded polyhedron from a three-dimensional exterior viewer. The pictures were constructed by computing the exponential map of a polyhedron by cutting and rotating faces into the tangent plane of the bug.


title: A Simple Depth of Search Metric for Exoplanet Imaging Surveys

abstract: We present a procedure for calculating expected exoplanet imaging yields, which explicitly separates the effects of instrument performance from assumptions of planet distributions. This `depth of search' approach allows for fast recalculation of yield values for variations in instrument parameters. We also describe a new target star selection metric with no dependence on an assumed planet population that can be used as a proxy for single-visit completeness. This approach allows for the recovery of the total mission completeness via convolution of the depth of search grid with an equivalent grid of assumed occurrence rates and integration over the part of the grid representing the population of interest (e.g., Earth-like planets on habitable zone orbits, etc.). In this work, we discuss the practical details of calculating the depth of search and present results of such calculations for one design iteration of the WFIRST coronagraphs.


title: Scalable Co-Optimization of Morphology and Control in Embodied Machines

abstract: Evolution sculpts both the body plans and nervous systems of agents together over time. In contrast, in AI and robotics, a robot's body plan is usually designed by hand, and control policies are then optimized for that fixed design. The task of simultaneously co-optimizing the morphology and controller of an embodied robot has remained a challenge. In psychology, the theory of embodied cognition posits that behavior arises from a close coupling between body plan and sensorimotor control, which suggests why co-optimizing these two subsystems is so difficult: most evolutionary changes to morphology tend to adversely impact sensorimotor control, leading to an overall decrease in behavioral performance. Here, we further examine this hypothesis and demonstrate a technique for "morphological innovation protection", which temporarily reduces selection pressure on recently morphologically-changed individuals, thus enabling evolution some time to "readapt" to the new morphology with subsequent control policy mutations. We show the potential for this method to avoid local optima and converge to similar highly fit morphologies across widely varying initial conditions, while sustaining fitness improvements further into optimization. While this technique is admittedly only the first of many steps that must be taken to achieve scalable optimization of embodied machines, we hope that theoretical insight into the cause of evolutionary stagnation in current methods will help to enable the automation of robot design and behavioral training -- while simultaneously providing a testbed to investigate the theory of embodied cognition.


title: Using Social Media to Predict the Future: A Systematic Literature Review

abstract: Social media (SM) data provides a vast record of humanity's everyday thoughts, feelings, and actions at a resolution previously unimaginable. Because user behavior on SM is a reflection of events in the real world, researchers have realized they can use SM in order to forecast, making predictions about the future. The advantage of SM data is its relative ease of acquisition, large quantity, and ability to capture socially relevant information, which may be difficult to gather from other data sources. Promising results exist across a wide variety of domains, but one will find little consensus regarding best practices in either methodology or evaluation. In this systematic review, we examine relevant literature over the past decade, tabulate mixed results across a number of scientific disciplines, and identify common pitfalls and best practices. We find that SM forecasting is limited by data biases, noisy data, lack of generalizable results, a lack of domain-specific theory, and underlying complexity in many prediction tasks. But despite these shortcomings, recurring findings and promising results continue to galvanize researchers and demand continued investigation. Based on the existing literature, we identify research practices which lead to success, citing specific examples in each case and making recommendations for best practices. These recommendations will help researchers take advantage of the exciting possibilities offered by SM platforms.


title: Localization for the one-dimensional Anderson model via positivity and large deviations for the Lyapunov exponent

abstract: We provide a complete and self-contained proof of spectral and dynamical localization for the one-dimensional Anderson model, starting from the positivity of the Lyapunov exponent provided by F\"urstenberg's theorem. That is, a Schr\"odinger operator in $\ell^2(\mathbb{Z})$ whose potential is given by independent identically distributed (i.i.d.) random variables almost surely has pure point spectrum with exponentially decaying eigenfunctions and its unitary group exhibits exponential off-diagonal decay, uniformly in time. This is achieved by way of a new result: for the Anderson model, one typically has Lyapunov behavior for all generalized eigenfunctions. We also explain how to obtain analogous statements for extended CMV matrices whose Verblunsky coefficients are i.i.d., as well as for half-line analogs of these models.


title: Element-centric clustering comparison unifies overlaps and hierarchy

abstract: Clustering is one of the most universal approaches for understanding complex data. A pivotal aspect of clustering analysis is quantitatively comparing clusterings; clustering comparison is the basis for many tasks such as clustering evaluation, consensus clustering, and tracking the temporal evolution of clusters. In particular, the extrinsic evaluation of clustering methods requires comparing the uncovered clusterings to planted clusterings or known metadata. Yet, as we demonstrate, existing clustering comparison measures have critical biases which undermine their usefulness, and no measure accommodates both overlapping and hierarchical clusterings. Here we unify the comparison of disjoint, overlapping, and hierarchically structured clusterings by proposing a new element-centric framework: elements are compared based on the relationships induced by the cluster structure, as opposed to the traditional cluster-centric philosophy. We demonstrate that, in contrast to standard clustering similarity measures, our framework does not suffer from critical biases and naturally provides unique insights into how the clusterings differ. We illustrate the strengths of our framework by revealing new insights into the organization of clusters in two applications: the improved classification of schizophrenia based on the overlapping and hierarchical community structure of fMRI brain networks, and the disentanglement of various social homophily factors in Facebook social networks. The universality of clustering suggests far-reaching impact of our framework throughout all areas of science.


title: A Low-power Reversible Alkali Atom Source

abstract: An electrically-controllable, solid-state, reversible device for sourcing and sinking alkali vapor is presented. When placed inside an alkali vapor cell, both an increase and decrease of the rubidium vapor density by a factor of two are demonstrated through laser absorption spectroscopy on 10 to 15 s time scales. The device requires low voltage (5 V), low power (<3.4 mW peak power), and low energy (<10.7 mJ per 10 s pulse). The absence of oxygen emission during operation is shown through residual gas analysis, indicating Rb is not lost through chemical reaction but rather by ion transport through the designed channel. This device is of interest for atomic physics experiments and, in particular, for portable cold-atom systems where dynamic control of alkali vapor density can enable advances in science and technology.


title: Study of the weak annihilation contributions in charmless $B_s\to VV$ decays

abstract: In this paper, in order to probe the spectator-scattering and weak annihilation contributions in charmless $B_s\to VV$ (where $V$ stands for a light vector meson) decays, we perform the $\chi^2$-analyses for the end-point parameters within the QCD factorization framework, under the constraints from the measured $\bar B_{s}\to$$\rho^0\phi$, $\phi K^{*0}$, $\phi \phi$ and $K^{*0}\bar K^{*0}$ decays. The fitted results indicate that the end-point parameters in the factorizable and nonfactorizable annihilation topologies are non-universal, which is also favored by the charmless $B\to PP$ and $PV$ (where $P$ stands for a light pseudo-scalar meson) decays observed in the previous work. Moreover, the abnormal polarization fractions $f_{L,\bot}(\bar B_{s}\to K^{*0}\bar K^{*0})=(20.1\pm7.0)\%\,,(58.4\pm8.5)\%$ measured by the LHCb collaboration can be reconciled through the weak annihilation corrections. However, the branching ratio of $\bar B_{s}\to\phi K^{*0}$ decay exhibits a tension between the data and theoretical result, which dominates the contributions to $\chi_{\rm min}^2$ in the fits. Using the fitted end-point parameters, we update the theoretical results for the charmless $B_s\to VV$ decays, which will be further tested by the LHCb and Belle-II experiments in the near future.


title: Faster Algorithms for Mean-Payoff Parity Games

abstract: Graph games provide the foundation for modeling and synthesis of reactive processes. Such games are played over graphs where the vertices are controlled by two adversarial players. We consider graph games where the objective of the first player is the conjunction of a qualitative objective (specified as a parity condition) and a quantitative objective (specified as a mean-payoff condition). There are two variants of the problem, namely, the threshold problem where the quantitative goal is to ensure that the mean-payoff value is above a threshold, and the value problem where the quantitative goal is to ensure the optimal mean-payoff value; in both cases ensuring the qualitative parity objective. The previous best-known algorithms for game graphs with $n$ vertices, $m$ edges, parity objectives with $d$ priorities, and maximal absolute reward value $W$ for mean-payoff objectives, are as follows: $O(n^{d+1} \cdot m \cdot W)$ for the threshold problem, and $O(n^{d+2} \cdot m \cdot W)$ for the value problem. Our main contributions are faster algorithms, and the running times of our algorithms are as follows: $O(n^{d-1} \cdot m \cdot W)$ for the threshold problem, and $O(n^{d} \cdot m \cdot W \cdot \log (n\cdot W))$ for the value problem. For mean-payoff parity objectives with two priorities, our algorithms match the best-known bounds of the algorithms for mean-payoff games (without conjunction with parity objectives). Our results are relevant in synthesis of reactive systems with both functional requirement (given as a qualitative objective) and performance requirement (given as a quantitative objective).


title: Beyond the "faster is slower" effect

abstract: The "faster is slower" effect raises when crowded people push each other to escape through an exit during an emergency situation. As individuals push harder, a statistical slowing down in the evacuation time can be achieved. The slowing down is caused by the presence of small groups of pedestrians (say, a small human cluster) that temporary blocks the way out when trying to leave the room. The pressure on the pedestrians belonging to this blocking cluster raises for increasing anxiety levels and/or larger number of individuals trying to leave the room through the same door. Our investigation shows, however, that very high pressures alters the dynamics in the blocking cluster, and thus, changes the statistics of the time delays along the escaping process. It can be acknowledged a reduction in the long lasting delays, while the overall evacuation performance improves. We present results on this novel phenomenon taking place beyond the "faster is slower" regime.


title: A fast algorithm for regularized focused 3-D inversion of gravity data using the randomized SVD

abstract: A fast algorithm for solving the under-determined 3-D linear gravity inverse problem based on the randomized singular value decomposition (RSVD) is developed. The algorithm combines an iteratively reweighted approach for $L_1$-norm regularization with the RSVD methodology in which the large scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low rank approximation of the system matrix with m rows is q=m, acceptable results are achievable with q<<m. In contrast to the use of the LSQR algorithm for the solution of the linear systems at each iteration, the singular values generated using the RSVD yield a good approximation of the dominant singular values of the large scale system matrix. The regularization parameter found for the small system at each iteration is thus dependent on the dominant singular values of the large scale system matrix and appropriately regularizes the dominant singular space of the large scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but are obtained at reduced computational cost. The method has been tested on synthetic models along with the real gravity data from the Morro do Engenho complex from central Brazil.


title: Performance Analysis of Inband FD-D2D Communications with Imperfect SI Cancellation for Wireless Video Distribution

abstract: Tremendous growing demand for high data rate services is the main driver for increasing traffic in wireless cellular networks. Device-to-Device (D2D) communications have recently been proposed to offload data via direct communications by bypassing cellular base stations (BSs). Such an offloading schemes increase capacity and reduce end-to-end delay in cellular networks and help to serve the dramatically increasing demand for high data rate. In this paper, we aim to analyze inband fullduplex (FD) D2D performance for the wireless video distribution by considering imperfect self-interference (SI) cancellation. Using tools from stochastic geometry, we analyze outage probability and spectral efficiency. Analytic and simulation results are used to demonstrate achievable gain against its half-duplex (HD) counterpart.


title: Matter growth in extended $\Lambda$CDM cosmology

abstract: On the basis of a previously established scalar-tensor extension of the $\Lambda$CDM model we develop an effective fluid approach for the matter growth function. This extended $\Lambda$CDM (henceforth $e_{\Phi}\Lambda$CDM) cosmology takes into account deviations from the standard model both via a modified background expansion and by the inclusion of geometric anisotropic stresses as well as of perturbations of the geometric dark-energy equivalent. The background dynamics is governed by an explicit analytic expression for the Hubble rate in which modifications of the standard model are given in terms of a single constant parameter [1]. To close the system of fluid-dynamical perturbation equations we introduce two phenomenological parameters through which the anisotropic stress is related both to the total energy density perturbation of the cosmic substratum and to relative perturbations in the effective two-component system. We quantify the impact of deviations from the standard background, of anisotropic stresses and of non-vanishing perturbations of the effective dark-energy component on the matter growth rate function $f \sigma_8$ and confront the results with recent redshift-space distortion (RSD) measurements.


title: Non-optimality of the Greedy Algorithm for subspace orderings in the method of alternating projections

abstract: The method of alternating projections involves projecting an element of a Hilbert space cyclically onto a collection of closed subspaces. It is known that the resulting sequence always converges in norm and that one can obtain estimates for the rate of convergence in terms of quantities describing the geometric relationship between the subspaces in question, namely their pairwise Friedrichs numbers. We consider the question of how best to order a given collection of subspaces so as to obtain the best estimate on the rate of convergence. We prove, by relating the ordering problem to a variant of the famous Travelling Salesman Problem, that correctness of a natural form of the Greedy Algorithm would imply that $\mathrm{P}=\mathrm{NP}$, before presenting a simple example which shows that, contrary to a claim made in the influential paper [Kayalar-Weinert, Math. Control Signals Systems, vol. 1(1), 1988], the result of the Greedy Algorithm is not in general optimal. We go on to establish sharp estimates on the degree to which the result of the Greedy Algorithm can differ from the optimal result. Underlying all of these results is a construction which shows that for any matrix whose entries satisfy certain natural assumptions it is possible to construct a Hilbert space and a collection of closed subspaces such that the pairwise Friedrichs numbers between the subspaces are given precisely by the entries of that matrix.


title: Recognizing and testing isomorphism of Cayley graphs over an abelian group of order $4p$ in polynomial time

abstract: We construct a polynomial-time algorithm that given a graph $X$ with $4p$ vertices ($p$ is prime), finds (if any) a Cayley representation of $X$ over the group $C_2\times C_2\times C_p$. This result, together with the known similar result for circulant graphs, shows that recognising and testing isomorphism of Cayley graphs over an abelian group of order $4p$ can be done in polynomial time.


title: A dynamic model for the two-parameter Dirichlet process

abstract: Let $\alpha=1/2$, $\theta>-1/2$, and $\nu_0$ be a probability measure on a type space $S$. In this paper, we investigate the stochastic dynamic model for the two-parameter Dirichlet process $\Pi_{\alpha,\theta,\nu_0}$. If $S=\mathbb{N}$, we show that the bilinear form \begin{eqnarray*} \left\{ \begin{array}{l} {\cal E}(F,G)=\frac{1}{2}\int_{{\cal P}_1(\mathbb{N})}\langle \nabla F(\mu),\nabla G(\mu)\rangle_{\mu} \Pi_{\alpha,\theta,\nu_0}(d\mu),\ \ F,G\in {\cal F},\\ {\cal F}=\{F(\mu)=f(\mu(1),\dots,\mu(d)):f\in C^{\infty}(\mathbb{R}^d), d\ge 1\} \end{array} \right. \end{eqnarray*} is closable on $L^2({\cal P}_1(\mathbb{N});\Pi_{\alpha,\theta,\nu_0})$ and its closure $({\cal E}, D({\cal E}))$ is a quasi-regular Dirichlet form. Hence $({\cal E}, D({\cal E}))$ is associated with a diffusion process in ${\cal P}_1(\mathbb{N})$ which is time-reversible with the stationary distribution $\Pi_{\alpha,\theta,\nu_0}$. If $S$ is a general locally compact, separable metric space, we discuss properties of the model \begin{eqnarray*} \left\{ \begin{array}{l} {\cal E}(F,G)=\frac{1}{2}\int_{{\cal P}_1(S)}\langle \nabla F(\mu),\nabla G(\mu)\rangle_{\mu} \Pi_{\alpha,\theta,\nu_0}(d\mu),\ \ F,G\in {\cal F},\\ {\cal F}=\{F(\mu)=f(\langle \phi_1,\mu\rangle,\dots,\langle \phi_d,\mu\rangle): \phi_i\in B_b(S),1\le i\le d,f\in C^{\infty}(\mathbb{R}^d),d\ge 1\}. \end{array} \right. \end{eqnarray*} In particular, we prove the Mosco convergence of its projection forms.


title: Transformation of Pan-STARRS1 gri to Stetson BVRI magnitudes. Photometry of small bodies observations

abstract: The UBVRI broad band photometric system is widely used in CCD astronomy. There are a lot of sets of standard stars for this photometric system, the Landolt's and Stetson's catalogues being the most precise and reliable. Another photometric system, recently considerably spread in CCD observations is ugriz, which originates from the Sloan Digital Sky Survey (SDSS) and has now many variations based on its 5 broad-band filters. One of the photometric systems based on it is The Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). In this paper we compare the BVRI magnitudes in the Stetson catalogue of standard stars with the magnitudes of the corresponding stars in the Pan-STARRS1 (PS1) grizyw catalogue. Transformations between these two systems are presented and discussed. An algorithm for data reduction and calibration is developed and its functionality is demonstrated in the magnitude determination of an asteroid.


title: On the spectra of geometric operators evolving with geometric flows

abstract: In this work we generalise various recent results on the evolution and monotonicity of the eigenvalues of certain geometric operators under specified geometric flows. Given a closed, compact Riemannian manifold $\big(M^n,g(t)\big)$ and a smooth function $\eta\in C^{\infty}(M)$ we consider the family of operators $\mathbb{L}=\Delta - g(\nabla\eta,\nabla\cdot)+cR$, where $R$ is the scalar curvature and $c$ is some real constant. We define a geometric flow on $M$ which encompasses the Ricci, the Ricci - Bourguignon and the Yamabe flows. Supposing that the metric $g(t)$ evolves along this general geometric flow we derive a formula for the evolution of the eigenvalues of $-\mathbb{L}$ and prove monotonicity results for the eigenvalues of both $-\Delta + g(\nabla\eta,\nabla\cdot)$ and $-\mathbb{L}$. We then prove Reilly-type formula for the operator $\mathbb{L}$ and employ it to establish an upper bound for the first variation of the eigenvalues of $-\mathbb{L}$. Finally, in the pursuit of a theoretical explanation of our generalisations, we formulate two conjectures on the monotonicity of the eigenvalues of Schr\"{o}dinger operators.


title: Proximity effects in bilayer graphene on monolayer WSe$_2$: Field-effect spin-valley locking, spin-orbit valve, and spin transistor

abstract: Proximity orbital and spin-orbit effects of bilayer graphene on monolayer WSe$_2$ are investigated from first-principles. We find that the built-in electric field induces an orbital band gap of about 10 meV in bilayer graphene. Remarkably, the proximity spin-orbit splitting for holes is two orders of magnitude---the spin-orbit splitting of the valence band at K is about 2 meV---more than for electrons. Effectively, holes experience spin-valley locking due to the strong proximity of the lower graphene layer to WSe$_2$. However, applying an external transverse electric field of some 1 V/nm, countering the built-in field of the heterostructure, completely reverses this effect and allows, instead for holes, electrons to be spin-valley locked with 2 meV spin-orbit splitting. Such a behavior constitutes a highly efficient field-effect spin-orbit valve, making bilayer graphene on WSe$_2$ a potential platform for a field-effect spin transistor.


title: A Comparison of Resampling and Recursive Partitioning Methods in Random Forest for Estimating the Asymptotic Variance Using the Infinitesimal Jackknife

abstract: The infinitesimal jackknife (IJ) has recently been applied to the random forest to estimate its prediction variance. These theorems were verified under a traditional random forest framework which uses classification and regression trees (CART) and bootstrap resampling. However, random forests using conditional inference (CI) trees and subsampling have been found to be not prone to variable selection bias. Here, we conduct simulation experiments using a novel approach to explore the applicability of the IJ to random forests using variations on the resampling method and base learner. Test data points were simulated and each trained using random forest on one hundred simulated training data sets using different combinations of resampling and base learners. Using CI trees instead of traditional CART trees as well as using subsampling instead of bootstrap sampling resulted in a much more accurate estimation of prediction variance when using the IJ. The random forest variations here have been incorporated into an open source software package for the R programming language.


title: Measurement of Minimum Bias Observables with ATLAS

abstract: The modelling of minimum bias interactions is a crucial ingredient to learn about the description of soft QCD processes. It has also a significant relevance for the simulation of the environment at the LHC with many concurrent pp interactions (pileup). The ATLAS collaboration has provided new measurements of the inclusive charged-particle multiplicity and its dependence on transverse momentum and pseudorapidity in special data sets with low LHC beam currents, recorded at center of mass energies of 8 TeV and 13 TeV. The measurements use charged-particle selections with minimum transverse momentum of both 100 MeV and 500 MeV, and various phase-space regions characterized by low and high charged-particle multiplicities.


title: The economical 3-3-1 model revisited

abstract: We show that the economical 3-3-1 model poses a very high new physics scale of the order of 1000~TeV due to the constraint on the flavor-changing neutral current. The implications of the model for neutrino masses, inflation, leptogenesis, and superheavy dark matter are newly recognized. Alternatively, we modify the model by rearranging the third quark generation differently from the first two quark generations, as well as changing the scalar sector. The resultant model now predicts a consistent new physics at TeV scale unlike the previous case and may be fully probed at the current colliders. Particularly, due to the minimal particle contents, the models under consideration manifestly accommodate dark matter candidates and neutrino masses, with novel and distinct production mechanisms. The large flavor-changing neutral currents that come from the ordinary and exotic quark mixings can be avoided due to the approximate $B-L$ symmetry.


title: Linking microscopic and macroscopic response in disordered solids

abstract: The modulus of a rigid network of harmonic springs depends on the sum of the energies in each of the bonds due to the applied distortion: compression in the case of the bulk modulus, $B$, or shear in the case of the shear modulus, $\mathcal{G}$. The distortion need not be global and we introduce a local modulus, $L_{i}$, associated with changing the equilibrium length of a single bond, $i$, in the network. We show that $L_{i}$ is useful for understanding many aspects of the mechanical response of the entire system. For example, it allows an understanding, and efficient computation, of how each bond in a network contributes to global properties such as $B$ and $\mathcal{G}$ and sheds light on how a particular bond's contribution to one modulus is, or is not, correlated with its contribution to another.


title: Analysis of the SFR - M* plane at z<3: single fitting versus multi-Gaussian decomposition

abstract: The analysis of galaxies on the star formation rate - stellar mass (SFR-M*) plane is a powerful diagnostic for galaxy evolution at different cosmic times. We consider a sample of 24463 galaxies from the CANDELS/GOODS-S survey to conduct a detailed analysis of the SFR-M* relation at redshifts 0.5$\leqslant z<$3 over more than three dex in stellar mass. To obtain SFR estimates, we utilise mid- and far-IR photometry when available, and rest-UV fluxes for all the other galaxies. We perform our analysis in different redshift bins, with two different methods: 1) a linear regression fitting of all star-forming galaxies, defined as those with specific star formation rates $\rm log_{10}(sSFR/yr^{-1}) > -9.8$, similarly to what is typically done in the literature; 2) a multi-Gaussian decomposition to identify the galaxy main sequence (MS), the starburst sequence and the quenched galaxy cloud. We find that the MS slope becomes flatter when higher stellar mass cuts are adopted, and that the apparent slope change observed at high masses depends on the SFR estimation method. In addition, the multi-Gaussian decomposition reveals the presence of a starburst population which increases towards low stellar masses and high redshifts. We find that starbursts make up ~5% of all galaxies at z=0.5-1.0, while they account for ~16% of galaxies at 2$<z<$3 with log$_{10}(M^{*})=$8.25-11.25. We conclude that the dissection of the SFR-M* in multiple components over a wide range of stellar masses is necessary to understand the importance of the different modes of star formation through cosmic time.


title: A recipe for echoes from exotic compact objects

abstract: Gravitational wave astronomy provides an unprecedented opportunity to test the nature of black holes and search for exotic, compact alternatives. Recent studies have shown that exotic compact objects (ECOs) can ring down in a manner similar to black holes, but can also produce a sequence of distinct pulses resembling the initial ringdown. These "echoes" would provide definite evidence for the existence of ECOs. In this work we study the generation of these echoes in a generic, parametrized model for the ECO, using Green's functions. We show how to reprocess radiation in the near-horizon region of a Schwarzschild black hole into the asymptotic radiation from the corresponding source in an ECO spacetime. Our methods allow us to understand the connection between distinct echoes and ringing at the resonant frequencies of the compact object. We find that the quasinormal mode ringing in the black hole spacetime plays a central role in determining the shape of the first few echoes. We use this observation to develop a simple template for echo waveforms. This template preforms well over a variety of ECO parameters, and with improvements may prove useful in the analysis of gravitational waves.