title: 3D Radiation Non-ideal Magnetohydrodynamical Simulations Of The Inner Rim In Protoplanetary Disks

abstract: Many planets orbit within an AU of their stars, raising questions about their origins. Particularly puzzling are the planets found near the silicate sublimation front. We investigate conditions near the front in the protostellar disk around a young intermediate-mass star, using the first global 3-D radiation non-ideal MHD simulations in this context. We treat the starlight heating; the silicate grains sublimation and deposition at the local, time-varying temperature and density; temperature-dependent Ohmic dissipation; and various initial magnetic fields. The results show magnetorotational turbulence around the sublimation front at 0.5 AU. The disk interior to 0.8 AU is turbulent, with velocities exceeding 10% of the sound speed. Beyond 0.8 AU is the dead zone, cooler than 1000 K and with turbulence orders of magnitude weaker. A local pressure maximum just inside the dead zone concentrates solid particles, favoring their growth. Over many orbits, a vortex develops at the dead zones inner edge, increasing the disks thickness locally by around 10%. We synthetically observe the results using Monte Carlo transfer calculations, finding the sublimation front is near-infrared bright. The models with net vertical magnetic fields develop extended, magnetically-supported atmospheres that reprocess extra starlight, raising the near-infrared flux 20%.The vortex throws a nonaxisymmetric shadow on the outer disk. At wavelengths > 2 micron, the flux varies several percent on monthly timescales. The variations are more regular when the vortex is present. The vortex is directly visible as an arc at ultraviolet through near-infrared wavelengths, given sub-AU spatial resolution.


title: Coupling Distributed and Symbolic Execution for Natural Language Queries

abstract: Building neural networks to query a knowledge base (a table) with natural language is an emerging research topic in deep learning. An executor for table querying typically requires multiple steps of execution because queries may have complicated structures. In previous studies, researchers have developed either fully distributed executors or symbolic executors for table querying. A distributed executor can be trained in an end-to-end fashion, but is weak in terms of execution efficiency and explicit interpretability. A symbolic executor is efficient in execution, but is very difficult to train especially at initial stages. In this paper, we propose to couple distributed and symbolic execution for natural language queries, where the symbolic executor is pretrained with the distributed executor's intermediate execution results in a step-by-step fashion. Experiments show that our approach significantly outperforms both distributed and symbolic executors, exhibiting high accuracy, high learning efficiency, high execution efficiency, and high interpretability.


title: Joint Hand Detection and Rotation Estimation by Using CNN

abstract: Hand detection is essential for many hand related tasks, e.g. parsing hand pose, understanding gesture, which are extremely useful for robotics and human-computer interaction. However, hand detection in uncontrolled environments is challenging due to the flexibility of wrist joint and cluttered background. We propose a deep learning based approach which detects hands and calibrates in-plane rotation under supervision at the same time. To guarantee the recall, we propose a context aware proposal generation algorithm which significantly outperforms the selective search. We then design a convolutional neural network(CNN) which handles object rotation explicitly to jointly solve the object detection and rotation estimation tasks. Experiments show that our method achieves better results than state-of-the-art detection models on widely-used benchmarks such as Oxford and Egohands database. We further show that rotation estimation and classification can mutually benefit each other.


title: Update of the Hadronic Vacuum Polarisation Contribution to the muon g-2

abstract: Precise data on e^+e^- to hadrons have recently become available and are used to compute the lowest-order hadronic vacuum polarisation contribution to the muon magnetic anomaly through dispersion relations. This is the case for the dominant pi+ pi- channel, but the most significant progress comes from the near completion of the BABAR program of measuring exclusive processes below 2 GeV with the initial-state radiation method which allows an efficient coverage of a large range of energies.. In this paper we briefly review the data treatment, the achieved improvements, and the result obtained for the full Standard Model prediction of the muon magnetic anomaly. The value obtained, a_mu (had~LO)=(692.6 +- 3.3)x 10^{-10} is 20% more precise than our last estimate in 2010. It deviates from the direct experimental determination by (27.4 +- 7.6)x 10^{-10} (3.6 sigma). Perpectives for further improvement are discussed.


title: Ionic Coulomb blockade and anomalous mole fraction effect in NaChBac bacterial ion channels

abstract: We report an experimental study of the influences of the fixed charge and bulk ionic concentrations on the conduction of biological ion channels, and we consider the results within the framework of the ionic Coulomb blockade model of permeation and selectivity. Voltage clamp recordings were used to investigate the Na$^+$/Ca$^{2+}$ anomalous mole fraction effect (AMFE) exhibited by the bacterial sodium channel NaChBac and its mutants. Site-directed mutagenesis was used to study the effect of either increasing or decreasing the fixed charge in their selectivity filters for comparison with the predictions of the Coulomb blockade model. The model was found to describe well some aspects of the experimental (divalent blockade and AMFE) and simulated (discrete multi-ion conduction and occupancy band) phenomena, including a concentration-dependent shift of the Coulomb staircase. These results substantially extend the understanding of ion channel selectivity and may also be applicable to biomimetic nanopores with charged walls.


title: Instantaneously complete Yamabe flow on hyperbolic space

abstract: We prove global existence of instantaneously complete Yamabe flows on hyperbolic space of arbitrary dimension $m\geq3$. The initial metric is assumed to be conformally hyperbolic with conformal factor and scalar curvature bounded from above. We do not require initial completeness or bounds on the Ricci curvature. If the initial data are rotationally symmetric, the solution is proven to be unique in the class of instantaneously complete, rotationally symmetric Yamabe flows.


title: Potential reduction interior point algorithm unfolding of neutron energy spectrum measured by recoil proton method

abstract: A reconstruction algorithm for unfolding the neutron energy spectrum has been developed based on the potential reduction interior point. This algorithm can be easily applied to the neutron energy spectrum reconstruction in the recoil proton method. We transform neutron energy spectrum unfolding problem into a typical nonnegative linear complementarity problem. The recoil proton energy spectrum and response matrix at angles of 0{\deg} and 30{\deg} are generated by the Geant4 simulation toolkit. Several different neutron energy test spectrums are also employed. It is found that this unfolding algorithm is stable and provide efficiency accurate results.


title: An approach to the topological complexity of the Klein bottle

abstract: Recently, Cohen and Vandembroucq proved that the reduced topological complexity of the Klein bottle is 4. Simultaneously and independently, we announced a proof of the same result. Mistakes were found in our argument, which was quite different than theirs. After correcting these, we found that our description of the obstruction class agreed with theirs. Our approach to showing that this obstruction is nonzero failed to do so, while theirs did not fail. Here we discuss our approach, which deals more directly with the simplicial structure of the Klein bottle.


title: Two-dimensional ReN$_2$ materials from first principles

abstract: The discovery of graphene makes it highly desirable to seek new two-dimensional materials. Through first-principles investigation, we predict two-dimensional materials of ReN$_{2}$: honeycomb and tetragonal structures. The phonon spectra establish the dynamical stability for both of the two structures, and the calculated in-plane stiffness constants proves their mechanical stability. The energy bands near the Fermi level consist of N-p and Re-d orbitals for the honeycomb structure, and are mainly from Re d orbitals for the tetragonal structure. While the tetragonal structure is non-magnetic, the honeycomb structure has N-based ferromagnetism, which will transit to anti-ferromagnetism under 14$\%$ biaxial strain. The calculated electron localization function and spin density indicate that direct N-N bond can occur only in the honeycomb structure. The ferromagnetism allows us to distinguish the two 2D phases easily. The tetragonal phase has lower energy than the honeycomb one, which means that the tetragonal phase is more stable, but the hexagonal phase has much larger bulk, shear, and Young's muduli than the tetragonal phase. The tetragonal phase is a three-bands metal, and the hexagonal phase is a ferromagnetic semi-metal. The special structural, electronic, magnetic, and optical properties in the honeycomb and tetragonal structures make them promising for novel applications.


title: Adaptive mesh point selection for the efficient solution of scalar IVPs

abstract: We discuss adaptive mesh point selection for the solution of scalar IVPs. We consider a method that is optimal in the sense of the speed of convergence, and aim at minimizing the local errors. Although the speed of convergence cannot be improved by using the adaptive mesh points compared to the equidistant points, we show that the factor in the error expression can be significantly reduced. We obtain formulas specifying the gain achieved in terms of the number of discretization subintervals, as well as in terms of the prescribed level of the local error. Both nonconstructive and constructive versions of the adaptive mesh selection are shown, and a numerical example is given.


title: QEX: a framework for lattice field theories

abstract: We present a new software framework for simulating lattice field theories. It features an intuitive programming interface, while simultaneously achieving high performance supercomputing, all in one programming language, Nim. With a macro system based on its abstract syntax tree, the language enables us to check and optimize our code at compile time. It also allows us to code intrinsics that map directly to machine instructions, and generates efficient native code. We show how we use Nim's metaprogramming features in our code, and present the current status of the code and future plans.


title: Protein-Ligand Scoring with Convolutional Neural Networks

abstract: Computational approaches to drug discovery can reduce the time and cost associated with experimental assays and enable the screening of novel chemotypes. Structure-based drug design methods rely on scoring functions to rank and predict binding affinities and poses. The ever-expanding amount of protein-ligand binding and structural data enables the use of deep machine learning techniques for protein-ligand scoring. We describe convolutional neural network (CNN) scoring functions that take as input a comprehensive 3D representation of a protein-ligand interaction. A CNN scoring function automatically learns the key features of protein-ligand interactions that correlate with binding. We train and optimize our CNN scoring functions to discriminate between correct and incorrect binding poses and known binders and non-binders. We find that our CNN scoring function outperforms the AutoDock Vina scoring function when ranking poses both for pose prediction and virtual screening.


title: Phase segregation facilitates exfoliation of franckeite crystals to a single unit cell thickness

abstract: Weak interlayer van der Waals interactions in bulk crystals facilitate their mechanical exfoliation to monolayer and few-layer two-dimensional (2D) materials, which exhibit striking physical phenomena absent in their bulk form. Here we study a 2D form of a mineral franckeite and show that phase segregation into discrete layers at the sub-nanometre scale facilitates its layered structure and basal cleavage. This behaviour is likely to be common in a wider family of complex crystals and could be exploited for a single-step synthesis of van der Waals heterostructures, as an alternative to stacking of 2D materials. Mechanical exfoliation allowed us to produce crystals down to a single unit cell thickness and rationalise its basal cleavage by atomic-resolution scanning transmission electron microscopy (STEM). We demonstrate p-type electrical conductivity and remarkable electrochemical properties in exfoliated crystals, which shows promise for energy storage applications.


title: Integrability via geometry: dispersionless differential equations in three and four dimensions

abstract: We prove that the existence of a dispersionless Lax pair with spectral parameter for a nondegenerate hyperbolic second order partial differential equation (PDE) is equivalent to the canonical conformal structure defined by the symbol being Einstein-Weyl on any solution in 3D, and self-dual on any solution in 4D. The first main ingredient in the proof is a characteristic property for dispersionless Lax pairs. The second is the projective behaviour of the Lax pair with respect to the spectral parameter. Both are established for nondegenerate determined systems of PDEs of any order. Thus our main result applies more generally to any such PDE system whose characteristic variety is a quadric hypersurface.


title: On the non-equivalence of the Bernoulli and K properties in dimension four

abstract: We study skew products where the base is a hyperbolic automorphism of $\mathbb{T}^2$, the fiber is a smooth area preserving flow on $\mathbb{T}^2$ with one fixed point (of high degeneracy) and the skewing function is a smooth non coboundary with non-zero integral. The fiber dynamics can be represented as a special flow over an irrational rotation and a roof function with one power singularity. We show that for a full measure set of rotations the corresponding skew product is $K$ and not Bernoulli. As a consequence we get a natural class of volume-preserving diffeomorphisms of $\mathbb{T}^4$ which are $K$ and not Bernoulli.


title: Negative thermal expansion and local lattice distortion in the (Sc1-xTix)F3 and related solid solutions

abstract: Negative thermal expansion (NTE) is unusual but important property for control of thermal expansion. In the present study, the chemical modification is utilized to engineer controllable thermal expansion in cubic NTE ScF3. A broad window of the coefficient of thermal expansion (CTE = -1.51 ~ -3.4*10^-6 K-1, 300-800 K) has been achieved in Sc1-xMxF3 (M = Ti, Al and Ga). The long-range crystallographic structure of (Sc1-xTix)F3 adheres to the cubic Pm-3m symmetry according to the analysis of high-energy synchrotron X-ray powder diffraction. Pair distribution function (PDF) analysis of synchrotron X-ray total scattering was performed to investigate the local lattice distortion. It was found that the weakness of NTE has a close correlation with the local lattice distortion. Based on the coupled rotation model, it is presumably that this local distortion might dampen the transverse vibration of F atoms and thus reduce NTE. The present work provides a possible reference for the design of controllable NTE in open framework solids.


title: Smoothed Projections over Weakly Lipschitz Domains

abstract: We develop finite element exterior calculus over weakly Lipschitz domains. Specifically, we construct commuting projections from $L^p$ de~Rham complexes over weakly Lipschitz domains onto finite element de~Rham complexes. These projections satisfy uniform bounds for finite element spaces with bounded polynomial degree over shape-regular families of triangulations. Thus we extend the theory of finite element differential forms to polyhedral domains that are weakly Lipschitz but not strongly Lipschitz. As new mathematical tools, we use the collar theorem in the Lipschitz category, and we show that the degrees of freedom in finite element exterior calculus are flat chains in the sense of geometric measure theory.


title: Hierarchy through Composition with Linearly Solvable Markov Decision Processes

abstract: Hierarchical architectures are critical to the scalability of reinforcement learning methods. Current hierarchical frameworks execute actions serially, with macro-actions comprising sequences of primitive actions. We propose a novel alternative to these control hierarchies based on concurrent execution of many actions in parallel. Our scheme uses the concurrent compositionality provided by the linearly solvable Markov decision process (LMDP) framework, which naturally enables a learning agent to draw on several macro-actions simultaneously to solve new tasks. We introduce the Multitask LMDP module, which maintains a parallel distributed representation of tasks and may be stacked to form deep hierarchies abstracted in space and time.


title: Including heavy spin effects in the prediction of a $\bar{b} \bar{b} u d$ tetraquark with lattice QCD potentials

abstract: We investigate spin effects in four-quark systems consisting of two heavy anti-bottom quarks and two light up/down quarks. To this end we use the Born-Oppenheimer approximation. We utilize potentials of two static antiquarks in the presence of two quarks of finite mass computed via lattice QCD and solve a coupled-channel Schr\"odinger equation for the anti-bottom-anti-bottom separation. Without taking heavy quark spins into account this approach predicted a $u d \bar b \bar b$ tetraquark bound state with quantum numbers $I(J^P) = 0(1^+)$. We now extend this Born-Oppenheimer approach with coupled channel Schr\"odinger equations allowing us to incorporate effects due to the heavy $\bar b$ spins. We confirm the existence of the $u d \bar b \bar b$ tetraquark.


title: The Galilean Superstring

abstract: The action for a Galilean superstring is found from a non-relativistic limit of the closed Green-Schwarz (GS) superstring; it has zero tension and provides an example of a massless super-Galilean system A Wess-Zumino term leads to a topological central charge in the Galilean supersymmetry algebra, such that unitarity requires a upper bound on the total momentum. This Galilean-invariant bound, which is also implied by the classical phase-space constraints, is saturated by solutions of the superstring equations of motion that half-preserve supersymmetry. We discuss briefly the extension to the Galilean supermembrane.


title: Misiurewicz parameters and dynamical stability of polynomial-like maps of large topological degree

abstract: Given a family of polynomial-like maps of large topological degree, we relate the presence of Misiurewicz parameters to a growth condition of the postcritical volume. This allows us to generalize to this setting the theory of stability and bifurcation developed by Berteloot, Dupont and the author for endomorphisms of $\mathbb{P}^k$.


title: A Maximum A Posteriori Estimation Framework for Robust High Dynamic Range Video Synthesis

abstract: High dynamic range (HDR) image synthesis from multiple low dynamic range (LDR) exposures continues to be actively researched. The extension to HDR video synthesis is a topic of significant current interest due to potential cost benefits. For HDR video, a stiff practical challenge presents itself in the form of accurate correspondence estimation of objects between video frames. In particular, loss of data resulting from poor exposures and varying intensity make conventional optical flow methods highly inaccurate. We avoid exact correspondence estimation by proposing a statistical approach via maximum a posterior (MAP) estimation, and under appropriate statistical assumptions and choice of priors and models, we reduce it to an optimization problem of solving for the foreground and background of the target frame. We obtain the background through rank minimization and estimate the foreground via a novel multiscale adaptive kernel regression technique, which implicitly captures local structure and temporal motion by solving an unconstrained optimization problem. Extensive experimental results on both real and synthetic datasets demonstrate that our algorithm is more capable of delivering high-quality HDR videos than current state-of-the-art methods, under both subjective and objective assessments. Furthermore, a thorough complexity analysis reveals that our algorithm achieves better complexity-performance trade-off than conventional methods.


title: String Diagrams For Double Categories and Equipments

abstract: A popular graphical calculus for monoidal categories makes computations tactile and intuitive. Complicated diagram chases can be expressed in a few pictures and discovered by playing with a shoelace. Joyal and Street's proof of the soundness of this calculus says that any deformation of a diagram, any bending of the strings, describes the same morphism. In this paper, we extend the graphical calculus to double categories and proarrow equipments in order to bring the string diagrammatic method to formal category theory. Our main theorem proves this calculus sound with the help of Dawson and Pare's results on composition in double categories.


title: Computing the Bouligand derivative of a class of piecewise-differentiable flows

abstract: Event--selected $C^r$ vector fields yield piecewise-differentiable flows, which possess a continuous and piecewise-linear Bouligand (or B-)derivative; here we provide an algorithm for computing this B-derivative. The number of "pieces" of the piecewise-linear B-derivative is factorial ($d!$) in the dimension ($d$) of the space, precluding a polynomial-time algorithm. We show how an exponential number ($2^d$) of points can be used to represent the B-derivative as a piecewise-linear homeomorphism in such a way that evaluating the derivative reduces to linear algebra computations involving a matrix constructed from $d$ of these points.


title: Charged current neutrino interactions in hot and dense matter

abstract: We derive the charged current absorption rate of electron and anti-electron neutrinos in dense matter using a fully relativistic approach valid at arbitrary matter degeneracy. We include mean field energy shifts due to nuclear interactions and the corrections due to weak magnetism. The rates are derived both from the familiar Fermi's Golden Rule, and from the techniques of finite temperature field theory, and shown to be equivalent. In various limits, these results can also be used to calculate neutral current opacities. We find that some pieces of the response have been left out in previous derivations and their contribution at high density can be significant. Useful formulae and detailed derivations are presented and we provide a new open source implementation of these opacities for use in radiation hydrodynamic simulations of core-collapse supernovae and neutron star mergers.


title: In-vacuo-dispersion features for GRB neutrinos and photons

abstract: Over the last 15 years there has been considerable interest in the possibility of quantum-gravity-induced in-vacuo dispersion, the possibility that spacetime itself might behave essentially like a dispersive medium for particle propagation. Two very recent studies have exposed what might be in-vacuo dispersion features for GRB (gamma-ray-burst) neutrinos of energy in the range of 100 TeV and for GRB photons with energy in the range of 10 GeV. We here show that these two features are roughly compatible with a description such that the same effects apply over 4 orders of magnitude in energy. We also characterize quantitatively how rare it would be for such features to arise accidentally, as a result of (still unknown) aspects of the mechanisms producing photons at GRBs or as a result of background neutrinos accidentally fitting the profile of a GRB neutrino affected by in-vacuo dispersion.


title: Feedback Neural Network for Weakly Supervised Geo-Semantic Segmentation

abstract: Learning from weakly-supervised data is one of the main challenges in machine learning and computer vision, especially for tasks such as image semantic segmentation where labeling is extremely expensive and subjective. In this paper, we propose a novel neural network architecture to perform weakly-supervised learning by suppressing irrelevant neuron activations. It localizes objects of interest by learning from image-level categorical labels in an end-to-end manner. We apply this algorithm to a practical challenge of transforming satellite images into a map of settlements and individual buildings. Experimental results show that the proposed algorithm achieves superior performance and efficiency when compared with various baseline models.


title: Medium effects in ${\Lambda}K^+$ pair production by 2.83 GeV protons on nuclei

abstract: We study ${\Lambda}K^+$ pair production in the interaction of protons of 2.83 GeV kinetic energy with C, Cu, Ag, and Au target nuclei in the framework of the nuclear spectral function approach for incoherent primary proton--nucleon and secondary pion--nucleon production processes, and processes associated with the creation of intermediate ${\Sigma^0}K^+$ pairs. The approach accounts for the initial proton and final $\Lambda$ hyperon absorption, final $K^+$ meson distortion in nuclei, target nucleon binding, and Fermi motion, as well as nuclear mean-field potential effects on these processes. We calculate the $\Lambda$ momentum dependence of the absolute ${\Lambda}K^+$ yield from the target nuclei considered, in the kinematical conditions of the ANKE experiment, performed at COSY, within the different scenarios for the $\Lambda$-nucleus effective scalar potential. We show that the above observable is appreciably sensitive to this potential in the low-momentum region. Therefore, direct comparison of the results of our calculations with the data from the ANKE-at-COSY experiment can help to determine the above potential at finite momenta. We also demonstrate that the two-step pion--nucleon production channels dominate in the low-momentum ${\Lambda}K^+$ production in the chosen kinematics and, therefore, they have to be taken into account in the analysis of these data.


title: X-ray flux variability of active galactic nuclei observed using NuSTAR

abstract: We present results on a systematic study of flux variability on hourly time-scales in a large sample of active galactic nuclei (AGN) in the 3-79 keV band using data from Nuclear Spectroscopic Telescope Array. Our sample consists of 4 BL Lac objects (BL Lacs), 3 flat spectrum radio quasars (FSRQs) 24 Seyfert 1, 42 Seyfert 2 and 8 narrow line Seyfert 1 (NLSy1) galaxies. We find that in the 3-79 keV band, about 65% of the sources in our sample show significant variations on hourly time scales. Using Mann-Whitney U-test and Kolmogorov-Smirnov test, we find no difference in the variability behaviour between Seyfert 1 and 2 galaxies. The blazar sources (FSRQs and BL Lacs) in our sample, are more variable than Seyfert galaxies that include Seyfert 1 and Seyfert 2 in the soft (3-10 keV), hard (10-79 keV) and total (3-79 keV) bands. NLSy1 galaxies show the highest duty cycle of variability (87%), followed by BL Lacs (82%), Seyfert galaxies (56%) and FSRQs (23%). We obtained flux doubling/halving time in the hard X-ray band less than 10 min in 11 sources. The flux variations between the hard and soft bands in all the sources in our sample are consistent with zero lag.


title: Lorentz invariance violation and simultaneous emission of electromagnetic and gravitational waves

abstract: In this work, we compute some phenomenological bounds for the electromagnetic and massive gravitational high-derivative extensions supposing that it is possible to have an astrophysical process that generates simultaneously gravitational and electromagnetic waves. We present Lorentz invariance violating (LIV) higher-order derivative models, following the Myers-Pospelov approach, to electrodynamics and massive gravitational waves. We compute the corrected equation of motion of these models, their dispersion relations and the velocities. The LIV parameters for the gravitational and electromagnetic sectors, $\xi_{g}$ and $\xi_{\gamma}$, respectively, were also obtained for three different approaches: luminal photons, time delay of flight and the difference of graviton and photon velocities. These LIV parameters depend on the mass scales where the LIV-terms become relevant, $M$ for the electromagnetic sector and $M_{1}$ for the gravitational one. We obtain, using the values for $M$ and $M_{1}$ found in the literature, that $\xi_{g}\sim10^{-2}$, which is expected to be phenomenologically relevant and $\xi_{\gamma}\sim10^{3}$, which cannot be suitable for an effective LIV theory. However, we show that $\xi_{\gamma}$ can be interesting in a phenomenological point of view if $M\gg M_{1}$. Finally the relation between the variation of the velocities of the photon and the graviton in relation to the speed of light was calculated and resulted in $\Delta v_{g}/\Delta v_{\gamma}\lesssim1.82\times 10^{-3}$.


title: Autocatalytic Sets and RNA Secondary Structure

abstract: The dominant paradigm in origin of life research is that of an RNA world. However, despite experimental progress towards the spontaneous formation of RNA, the RNA world hypothesis still has its problems. Here, we introduce a novel computational model of chemical reaction networks based on RNA secondary structure and analyze the existence of autocatalytic sub-networks in random instances of this model, by combining two well-established computational tools. Our main results are that (i) autocatalytic sets are highly likely to exist, even for very small reaction networks and short RNA sequences, and (ii) sequence diversity seems to be a more important factor in the formation of autocatalytic sets than sequence length. These findings could shed new light on the probability of the spontaneous emergence of an RNA world as a network of mutually collaborative ribozymes.


title: Probing the BFKL dynamics in inclusive three jet production at the LHC

abstract: We propose the study of new observables in LHC inclusive events with three tagged jets, one in the forward direction, one in the backward direction and both well-separated in rapidity from the each other (Mueller-Navelet jets), together with a third jet tagged in central regions of rapidity. Since non-tagged associated mini-jet multiplicity is allowed, we argue that projecting the cross sections on azimuthal-angle components can provide several distinct tests of the BFKL dynamics. Realistic LHC kinematical cuts are introduced.


title: Chiral Tensors of Mixed Young Symmetry

abstract: Chiral tensors of mixed Young symmetry, which exist in the same spacetime dimensions $2 + 4n$ where chiral $p$-forms can be defined, are investigated. Such chiral tensors have been argued to play a central role in exotic formulations of gravity in 6 dimensions and possess intriguing properties. A variational principle that yields the chiral equations of motion is explicitly constructed and related to the action for a non-chiral tensor. The use of prepotentials turns out to be essential in our analysis. We also comment on dimensional reduction.


title: Cohomology of Polychromatic Configuration Spaces of Euclidean Space

abstract: Recently, the homology and cohomology of non-k-overlapping discs, or, equivalently, no k-equal subspaces of Euclidean space, were calculated by Dobrinskaya and Turchin. We calculate the homology and cohomology of two classes of more general spaces: decreasing polychromatic configuration spaces and bicolored configuration spaces. Instead of all points behaving similarly, we allow for varying behavior between points.


title: The distribution of radioactive $^{44}$Ti in Cassiopeia A

abstract: The distribution of elements produced in the inner-most layers of a supernova explosion is a key diagnostic for studying the collapse of massive stars. Here we present the results of a 2.4 Ms \textit{NuSTAR} observing campaign aimed at studying the supernova remnant Cassiopeia A (Cas A). We perform spatially-resolved spectroscopic analyses of the $^{44}$Ti ejecta which we use to determine the Doppler shift and thus the three-dimensional (3D) velocities of the $^{44}$Ti ejecta. We find an initial $^{44}$Ti mass of 1.54 $\pm$ 0.21 $\times 10^{-4}$ M$_{\odot}$ which has a present day average momentum direction of 340$^{\circ}$ $\pm$ 15$^{\circ}$ projected on to the plane of the sky (measured clockwise from Celestial North) and tilted by 58$^{\circ}$ $\pm$ 20$^{\circ}$ into the plane of the sky away from the observer, roughly opposite to the inferred direction of motion of the central compact object. We find some $^{44}$Ti ejecta that are clearly interior to the reverse shock and some that are clearly exterior to the reverse shock. Where we observe $^{44}$Ti ejecta exterior to the reverse shock we also see shock-heated iron; however, there are regions where we see iron but do not observe $^{44}$Ti. This suggests that the local conditions of the supernova shock during explosive nucleosynthesis varied enough to suppress the production of $^{44}$Ti in some regions by at least a factor of two, even in regions that are assumed to be the result of processes like $\alpha$-rich freezeout that should produce both iron and titanium.


title: Continuum limit and stochastic homogenization of discrete ferromagnetic thin films

abstract: We study the discrete-to-continuum limit of ferromagnetic spin systems when the lattice spacing tends to zero. We assume that the atoms are part of a (maybe) non-periodic lattice close to a flat set in a lower dimensional space, typically a plate in three dimensions. Scaling the particle positions by a small parameter $\varepsilon>0$ we perform a $\Gamma$-convergence analysis of properly rescaled interfacial-type energies. We show that, up to subsequences, the energies converge to a surface integral defined on partitions of the flat space. In the second part of the paper we address the issue of stochastic homogenization in the case of random stationary lattices. A finer dependence of the homogenized energy on the average thickness of the random lattice is analyzed for an example of magnetic thin system obtained by a random deposition mechanism.


title: Discovery of XO-6b: a hot Jupiter transiting a fast rotating F5 star on an oblique orbit

abstract: Only a few hot Jupiters are known to orbit around fast rotating stars. These exoplanets are harder to detect and characterize and may be less common than around slow rotators. Here, we report the discovery of the transiting hot Jupiter XO-6b, which orbits a bright, hot, and fast rotating star: V = 10.25, Teff = 6720 +/- 100 K, v sin i = 48 +/- 3 km/s. We detected the planet from its transits using the XO instruments and conducted a follow-up campaign. Because of the fast stellar rotation, radial velocities taken along the orbit do not yield the planet's mass with a high confidence level, but we secure a 3-sigma upper limit Mp < 4.4 MJup. We also obtain high resolution spectroscopic observations of the transit with the SOPHIE spectrograph at the 193-cm telescope of the Observatoire de Haute-Provence and analyze the stellar lines profile by Doppler tomography. The transit is clearly detected in the spectra. The radii measured independently from the tomographic analysis and from the photometric lightcurves are consistent, showing that the object detected by both methods is the same and indeed transits in front of XO-6. We find that XO-6b lies on a prograde and misaligned orbit with a sky-projected obliquity lambda = -20.7 +/- 2.3 deg. The rotation period of the star is shorter than the orbital period of the planet: Prot < 2.12 days, Porb = 3.77 days. Thus, this system stands in a largely unexplored regime of dynamical interactions between close-in giant planets and their host stars.


title: Inverse Electromagnetic Diffraction by Biperiodic Dielectric Gratings

abstract: Consider the incidence of a time-harmonic electromagnetic plane wave onto a biperiodic dielectric grating, where the surface is assumed to be a small and smooth perturbation of a plane. The diffraction is modeled as a transmission problem for Maxwell's equations in three dimensions. This paper concerns the inverse diffraction problem which is to reconstruct the grating surface from either the diffracted field or the transmitted field. A novel approach is developed to solve the challenging nonlinear and ill-posed inverse problem. The method requires only a single incident field and is realized via the fast Fourier transform. Numerical results show that it is simple, fast, and stable to reconstruct biperiodic dielectric grating surfaces with super-resolved resolution.


title: Continued Fractions for Square Series Generating Functions

abstract: We consider new series expansions for variants of the so-termed ordinary geometric square series generating functions originally defined in the recent article titled "Square Series Generating Function Transformations" (arXiv: 1609.02803). Whereas the original square series transformations article adapts known generating function transformations to construct integral representations for these square series functions enumerating the square powers of $q^{n^2}$ for some fixed non-zero $q$ with $|q| < 1$, we study the expansions of these special series through power series generated by Jacobi-type continued fractions, or J-fractions. We prove new exact expansions of the $h^{th}$ convergents to these continued fraction series and show that the limiting case of these convergent generating functions exists. We also prove new infinite $q$-series representations of special square series expansions involving square-power terms of the series parameter $q$, the $q$-Pochhammer symbol, and double sums over the $q$-binomial coefficients. Applications of the new results we prove within the article include new $q$-series representations for the ordinary generating functions of the special sequences, $r_p(n)$, and $\sigma_1(n)$, as well as parallels to the examples of the new integral representations for theta functions, series expansions of infinite products and partition function generating functions, and related unilateral special function series cited in the first square series transformations article.


title: Triclinic nematic colloidal crystals from competing elastic and electrostatic interactions

abstract: Self-assembly of nanoparticles can enable composites with pre-designed properties but remains challenged by reproducing structural diversity of atomic and molecular crystals. We combine anisotropic elastic and weakly screened electrostatic interactions to guide both orientational and triclinic positional self-ordering of inorganic nanocrystals in a nematic fluid host. The lattice periodicity of these low-symmetry colloidal crystals is more than an order of magnitude larger than the nanoparticle size. Orientations of nanocrystals, as well as crystallographic axes of ensuing triclinic colloidal crystals, are coupled to the uniform alignment direction of the nematic host, which can be readily controlled on large scales. We probe colloidal pair and many-body interactions and demonstrate how triclinic crystals with orientational ordering of the semiconductor nanorods emerge from competing long-range elastic and electrostatic forces.


title: Improved generator objectives for GANs

abstract: We present a framework to understand GAN training as alternating density ratio estimation and approximate divergence minimization. This provides an interpretation for the mismatched GAN generator and discriminator objectives often used in practice, and explains the problem of poor sample diversity. We also derive a family of generator objectives that target arbitrary $f$-divergences without minimizing a lower bound, and use them to train generative image models that target either improved sample quality or greater sample diversity.


title: Cooperative effect of doping and temperature on the polaronic band structure in strongly correlated electron systems with strong electron-phonon interaction

abstract: In this work we investigate doping and temperature dependences of electronic structure of system with strong electronic correlations and strong electron-phonon interaction modeling cuprates in the frameworks of the three-band p-d-Holstein model by a polaronic version of the generalized tight binding (GTB) method. Within this approach the electronic structure is formed by polaronic quasiparticles constructed as excitations between initial and final polaronic multielectron states. Doping and temperature effects are taken into account by occupation numbers of local excited polaronic states and variations in the magnitude of spin-spin correlation functions. Both effects are manifested in the reconstruction of band structure, Fermi contours, density of states and redistribution of the spectral weight over the Hubbard polaron subbands. Doping leads to transformation of Fermi contour from small hole pockets around k=(pi/2,pi/2) with inhomogeneous spectral weight distribution at small hole concentration to large contour around k=(pi,pi) in the overdoped compound as a result of two quantum phase transitions. In the system with phonon subsystem and EPI doping results in the top of the valence band splitting off and new polaron subbands appearance. Temperature increasing in the system with doped holes and moderate EPI leads to formation of the flatband around k=(pi,pi) and transfer of the spectral weight to the splitted off top of the valence band.


title: Ergodic Theory and the Structure of Noncommutative Space-Time

abstract: We develop further our fibre bundle construct of non-commutative space-time on a Minkowski base space. We assume space-time is non-commutative due to the existence of additional non-commutative algebraic structure at each point x of space-time, forming a quantum operator 'fibre algebra' A(x). This structure then corresponds to the single fibre of a fibre bundle. A gauge group acts on each fibre algebra locally, while a 'section' through this bundle is then a quantum field of the form {A(x); x in M} with M the underlying space-time manifold. In addition, we assume a local algebra O(D) corresponding to the algebra of sections of such a principal fibre bundle with base space a finite and bounded subset of space-time, D. The algebraic operations of addition and multiplication are assumed defined fibrewise for this algebra of sections. We characterise 'ergodic' extremal quantum states of the fibre algebra invariant under the subgroup T of local translations of space-time of the Poincare group P in terms of a non-commutative extension of entropy applied to the subgroup T. We also characterise the existence of T- invariant states by generalizing to the non-commutative case Kakutani's work on wandering projections. This leads on to a classification of the structure of the local algebra O(D) by using a 'T-Twisted' equivalence relation, including a full analysis of the T-type III case. In particular we show that O(D) is T-type III if and only if the crossed product algebra O(D)xT is type III in the sense of Murray-von Neumann.


title: Inequalities for means of two mean variables

abstract: Motivated by the work of Anderson, Vamanamurthy and Vuorinen \cite{avv}, in this paper authors study the log-convexity and log-concavity of Power mean, Identric mean, weighted Power mean, Lehmer mean, Modified Alzer mean, and establish the relation of these means with each other.


title: Jonckheere-Terpstra test for nonclassical error versus log-sensitivity relationship of quantum spin network controllers

abstract: Selective information transfer in spin ring networks by energy landscape shaping control has the property that the error 1-prob, where prob is the transfer success probability, and the sensitivity of the error to spin coupling uncertainties are statistically increasing across a family of controllers of increasing error. The need for a statistical Hypothesis Testing of a concordant trend is made necessary by the noisy behavior of the sensitivity versus the error as a consequence of the optimization of the controllers in a challenging error landscape. Here, we examine the concordant trend between the error and another measure of performance - the logarithmic sensitivity - used in robust control to formulate a well known fundamental limitation. Contrary to error versus sensitivity, the error versus logarithmic sensitivity trend is less obvious, because of the amplification of the noise due to the logarithmic normalization. This results in the Kendall {\tau} test for rank correlation between the error and the log sensitivity to be somewhat pessimistic with marginal significance level. Here it is shown that the Jonckheere-Terpstra test, because it tests the Alternative Hypothesis of an ordering of the medians of some groups of log sensitivity data, alleviates this statistical problem. This identifies cases of concordant trend between the error and the logarithmic sensitivity, a highly anti-classical features that goes against the well know sensitivity versus complementary sensitivity limitation.


title: Incorporating Quotation and Evaluation Into Church's Type Theory

abstract: ${\rm CTT}_{\rm qe}$ is a version of Church's type theory that includes quotation and evaluation operators that are similar to quote and eval in the Lisp programming language. With quotation and evaluation it is possible to reason in ${\rm CTT}_{\rm qe}$ about the interplay of the syntax and semantics of expressions and, as a result, to formalize syntax-based mathematical algorithms. We present the syntax and semantics of ${\rm CTT}_{\rm qe}$ as well as a proof system for ${\rm CTT}_{\rm qe}$. The proof system is shown to be sound for all formulas and complete for formulas that do not contain evaluations. We give several examples that illustrate the usefulness of having quotation and evaluation in ${\rm CTT}_{\rm qe}$.


title: On the muon transfer from protium to neon at low collision energies

abstract: The direct muon transfer from protium to neon at low collision energies is considered. The question is discussed how the energy dependence of the transfer rate can be extracted from available experimental data. A model is suggested which is based on two assumptions. First, the muon transfer occurs within an interaction sphere out of which the entrance and transfer channels are not coupled. Secondly, for the s- and p-waves at low collision energies, the complex logarithmic derivatives of the radial wave functions of the entrance channel on the interaction sphere are energy-independent. With only the s-wave taken into account, two values of the logarithmic derivative were found, which agree with the transfer rates measured at the temperatures of 20 and 300 K. Within this temperature interval, the rates evaluated for these two values differ only slightly. Noticeable differences arise at temperatures of a few kelvins. The choice of the preferred value of the logarithmic derivative was made by the comparison of the energy dependence of the S-matrix with results of previous calculations. Thereby, the proposed procedure may be a way to correct calculations at low collision energies.


title: The metallic winds in dwarf galaxies

abstract: We present results from models of galactic winds driven by energy injected from nuclear (at the galactic center) and non-nuclear starbursts. The total energy of the starburst is provided by very massive young stellar clusters,which can push the galactic interstellar medium and produce an important outflow. Such outflow can be a well, or partially mixed wind, or a highly metallic wind. We have performed adiabatic 3D N-Body/Smooth Particle Hydrodynamics simulations of galactic winds using the GADGET-2 code. The numerical models cover a wide range of parameters, varying the galaxy concentration index, gas fraction of the galactic disk, and radial distance of the starburst. We show that an off-center starburst in dwarf galaxies is the most effective mechanism to produce a significant loss of metals (material from the starburst itself). At the same time a non-nuclear starburst produce a high efficiency of metal loss, in spite of having a moderate to low mass loss rate.


title: Faster Space-Efficient Algorithms for Subset Sum, k-Sum and Related Problems

abstract: We present space efficient Monte Carlo algorithms that solve Subset Sum and Knapsack instances with $n$ items using $O^*(2^{0.86n})$ time and polynomial space, where the $O^*(\cdot)$ notation suppresses factors polynomial in the input size. Both algorithms assume random read-only access to random bits. Modulo this mild assumption, this resolves a long-standing open problem in exact algorithms for NP-hard problems. These results can be extended to solve Binary Linear Programming on $n$ variables with few constraints in a similar running time. We also show that for any constant $k\geq 2$, random instances of $k$-Sum can be solved using $O(n^{k-0.5}polylog(n))$ time and $O(\log n)$ space, without the assumption of random access to random bits. Underlying these results is an algorithm that determines whether two given lists of length $n$ with integers bounded by a polynomial in $n$ share a common value. Assuming random read-only access to random bits, we show that this problem can be solved using $O(\log n)$ space significantly faster than the trivial $O(n^2)$ time algorithm if no value occurs too often in the same list.


title: In the transmission of information, the great potential of model-based coding with the SP theory of intelligence

abstract: Model-based coding, described by John Pierce in 1961, has great potential to reduce the volume of information that needs to be transmitted in moving big data, without loss of information, from one place to another, or in lossless communications via the internet. Compared with ordinary compression methods, this potential advantage of model-based coding in the transmission of data arises from the fact that both the transmitter ("Alice") and the receiver ("Bob") are equipped with a grammar for the kind of data that is to be transmitted, which means that, to achieve lossless transmission of a body of data from Alice and Bob, a relatively small amount of information needs to be sent. Preliminary trials indicate that, with model-based coding, the volume of information to be sent from Alice to Bob to achieve lossless transmission of a given body of data may be less than $6\%$ of the volume of information that needs to be sent when ordinary compression methods are used. Until recently, it has not been feasible to convert John Pierce's vision into something that may be applied in practice. Now, with the development of the "SP theory of intelligence" and its realisation in the "SP computer model", there is clear potential to realise the three main functions that will be needed: unsupervised learning of a grammar for the kind of data that is to be transmitted using a relatively powerful computer that is independent of Alice and Bob; the encoding by Alice of any one example of such data in terms of the grammar; and, with the grammar, decoding of the encoding by Bob to retrieve the given example. It appears now to be feasible, within reasonable timescales, to bring these capabilities to a level where they may be applied to the transmission of realistically large bodies of data.


title: Six winters of photometry from Dome C, Antarctica: challenges, improvements, and results from the ASTEP experiment

abstract: ASTEP (Antarctica Search for Transiting ExoPlanets) is a pilot project that aims at searching and characterizing transiting exoplanets from Dome C in Antarctica and to qualify this site for photometry in the visible. Two instruments were installed at Dome C and ran for six winters in total. The analysis of the collected data is nearly complete. We present the operation of the instruments, and the technical challenges, limitations, and possible solutions in light of the data quality. The instruments performed continuous observations during the winters. Human interventions are required mainly for regular inspection and ice dust removal. A defrosting system is efficient at preventing and removing ice on the mirrors. The PSF FWHM is 4.5 arcsec on average which is 2.5 times larger than the specification, and is highly variable; the causes are the poor ground-level seeing, the turbulent plumes generated by the heating system, and to a lower extent the imperfect optical alignment and focusing, and some astigmatism. We propose solutions for each of these aspects that would largely increase the PSF stability. The astrometric and guiding precisions are satisfactory and would deserve only minor improvements. Major issues are encountered with the camera shutter which did not close properly after two winters; we minimized this issue by heating the shutter and by developing specific image calibration algorithms. Finally, we summarize the site testing and science results obtained with ASTEP. Overall, the ASTEP experiment will serve as a basis to design and operate future optical and near-infrared telescopes in Antarctica.


title: Unitary Correlation Sets

abstract: The unitary correlation sets defined by the first author in conjunction with tensor products of $\mathcal{U}_{nc}(n)$ are further studied. We show that Connes' embedding problem is equivalent to deciding whether or not two smaller versions of the unitary correlation sets are equal. Moreover, we obtain the result that Connes' embedding problem is equivalent to deciding whether or not two cross norms on $M_n \otimes M_n$ are equal for all $n \geq 2$.


title: Quantum propagation across cosmological singularities

abstract: The initial singularity is the most troubling feature of the standard cosmology, which quantum effects are hoped to resolve. In this paper, we study quantum cosmology with conformal (Weyl) invariant matter. We show it is natural to extend the scale factor to negative values, allowing a large, collapsing Universe to evolve across a quantum "bounce" into an expanding Universe like ours. We compute the Feynman propagator for Friedmann-Robertson-Walker backgrounds exactly, identifying curious pathologies in the case of curved (open or closed) universes. We then include anisotropies, fixing the operator ordering of the quantum Hamiltonian by imposing covariance under field redefinitions and again finding exact solutions. We show how complex classical solutions allow one to circumvent the singularity while maintaining the validity of the semiclassical approximation. The simplest isotropic universes sit on a critical boundary, beyond which there is qualitatively different behavior, with potential for instability. Additional scalars improve the theory's stability. Finally, we study the semiclassical propagation of inhomogeneous perturbations about the flat, isotropic case, at linear and nonlinear order, showing that, at least at this level, there is no particle production across the bounce. These results form the basis for a promising new approach to quantum cosmology and the resolution of the big bang singularity.


title: Calculation of Momentum Distribution Function of a Non-thermal Fermionic Dark Matter

abstract: The most widely studied scenario in dark matter phenomenology is the thermal WIMP scenario. Inspite of numerous efforts to detect WIMP, till now we have no direct evidence for it. A possible explanation for this non-observation of dark matter could be because of its very feeble interaction strength and hence, failing to thermalise with the rest of the cosmic soup. In other words, the dark matter might be of non-thermal origin where the relic density is obtained by the so-called freeze-in mechanism. Furthermore, if this non-thermal dark matter is itself produced substantially from the decay of another non-thermal mother particle, then their distribution functions may differ in both size and shape from the usual equilibrium distribution function. In this work, we have studied such a non-thermal (fermionic) dark matter scenario in the light of a new type of $U(1)_{\rm B-L}$ model. The $U(1)_{\rm B-L}$ model is interesting, since, besides being anomaly free, it can give rise to neutrino mass by Type II see-saw mechanism. Moreover, as we will show, it can accommodate a non-thermal fermionic dark matter as well. Starting from the collision terms, we have calculated the momentum distribution function for the dark matter by solving a coupled system of Boltzmann equations. We then used it to calculate the final relic abundance, as well as other relevant physical quantities. We have also compared our result with that obtained from solving the usual Boltzmann (or rate) equations directly in terms of comoving number density, $Y$. Our findings suggest that the latter approximation is valid only in cases where the system under study is close to equilibrium, and hence should be used with caution.


title: Change point detection in heteroscedastic time series

abstract: Many time series exhibit changes both in level and in variability. Generally, it is more important to detect a change in the level, and changing or smoothly evolving variability can confound existing tests. This paper develops a framework for testing for shifts in the level of a series which accommodates the possibility of changing variability. The resulting tests are robust both to heteroskedasticity and serial dependence. They rely on a new functional central limit theorem for dependent random variables whose variance can change or trend in a substantial way. This new result is of independent interest as it can be applied in many inferential contexts applicable to time series. Its application to change point tests relies on a new approach which utilizes Karhunen--Lo{\'e}ve expansions of the limit Gaussian processes. After presenting the theory in the most commonly encountered setting of the detection of a change point in the mean, we show how it can be extended to linear and nonlinear regression. Finite sample performance is examined by means of a simulation study and an application to yields on US treasury bonds.


title: Safety Verification and Control for Collision Avoidance at Road Intersections

abstract: This paper presents the design of a supervisory algorithm that monitors safety at road intersections and overrides drivers with a safe input when necessary. The design of the supervisor consists of two parts: safety verification and control design. Safety verification is the problem to determine if vehicles will be able to cross the intersection without colliding with current drivers' inputs. We translate this safety verification problem into a jobshop scheduling problem, which minimizes the maximum lateness and evaluates if the optimal cost is zero. The zero optimal cost corresponds to the case in which all vehicles can cross each conflict area without collisions. Computing the optimal cost requires solving a Mixed Integer Nonlinear Programming (MINLP) problem due to the nonlinear second-order dynamics of the vehicles. We therefore estimate this optimal cost by formulating two related Mixed Integer Linear Programming (MILP) problems that assume simpler vehicle dynamics. We prove that these two MILP problems yield lower and upper bounds of the optimal cost. We also quantify the worst case approximation errors of these MILP problems. We design the supervisor to override the vehicles with a safe control input if the MILP problem that computes the upper bound yields a positive optimal cost. We theoretically demonstrate that the supervisor keeps the intersection safe and is non-blocking. Computer simulations further validate that the algorithms can run in real time for problems of realistic size.


title: Characterizing a high spin magnetic impurity via Andreev reflection spectroscopy

abstract: The ground state properties of a high spin magnetic impurity and its interaction with an electronic spin are probed via Andreev reflection. We see that through the charge and spin conductance one can effectively estimate the interaction strength, the ground state spin and magnetic moment of any high spin magnetic impurity. We show how a high spin magnetic impurity at the junction between a normal metal and superconductor can contribute to superconducting spintronics applications. Particularly, while spin conductance is absent below the gap for Ferromagnet-Insulator-Superconductor junctions we show that in the case of a Normal metal-High spin magnetic impurity-Normal Metal-Insulator-Superconductor (NMNIS) junction it is present. Further, it is seen that pure spin conduction can exist without any accompanying charge conduction in the NMNIS junction.


title: Locality and Unitarity from Singularities and Gauge Invariance

abstract: We conjecture that the leading two-derivative tree-level amplitudes for gluons and gravitons can be derived from gauge invariance together with mild assumptions on their singularity structure. Assuming locality (that the singularities are associated with the poles of cubic graphs), we prove that gauge-invariance in just (n-1) particles together with minimal power-counting uniquely fixes the amplitude. Unitarity in the form of factorization then follows from locality and gauge invariance. We also give evidence for a stronger conjecture: assuming only that singularities occur when the sum of a subset of external momenta go on-shell, we show in non-trivial examples that gauge-invariance and power-counting demand a graph structure for singularities. Thus both locality and unitarity emerge from singularities and gauge invariance. Similar statements hold for theories of Goldstone bosons like the non-linear sigma model and Dirac-Born-Infeld, by replacing the condition of gauge invariance with an appropriate degree of vanishing in soft limits.


title: Octonionic two-qubit separability probability conjectures

abstract: We study, further, a conjectured formula for generalized two-qubit Hilbert-Schmidt separability probabilities that has recently been proven by Lovas and Andai (https://arxiv.org/pdf/1610.01410.pdf) for its real (two-rebit) asserted value ($\frac{29}{64}$), and that has also been very strongly supported numerically for its complex ($\frac{8}{33}$), and quaternionic ($\frac{26}{323}$) counterparts. Now, we seek to test the presumptive octonionic value of $\frac{44482}{4091349} \approx 0.0108722$. We are somewhat encouraged by certain numerical computations, indicating that this (51-dimensional) instance of the conjecture might be fulfilled by setting a certain determinantal-power parameter $a$, introduced by Forrester (https://arxiv.org/pdf/1610.08081.pdf), to 0 (or possibly near to 0). Hilbert-Schmidt measure being the case $k=0$ of random induced measure, for $k=1$, the corresponding octonionic separability probability conjecture is $\frac{7612846}{293213345} \approx 0.0259635$, while for $k=2$, it is $\frac{4893392}{95041567} \approx 0.0514869, \ldots$. The relation between the parameters $a$ and $k$ is explored.


title: On the volume of double twist link cone-manifolds

abstract: We consider the double twist link $J(2m+1, 2n+1)$ which is the two-bridge link corresponding to the continued fraction $(2m+1)-1/(2n+1)$. It is known that $J(2m+1, 2n+1)$ has reducible nonabelian $SL_2(\mathbb{C})$-character variety if and only if $m=n$. In this paper we give a formula for the volume of hyperbolic cone-manifolds of $J(2m+1,2m+1)$. We also give a formula for the A-polynomial 2-tuple corresponding to the canonical component of the character variety of $J(2m+1,2m+1)$.


title: Strong convergence of tamed $\theta$-EM scheme for neutral SDDEs with one-sided Lipschitz drift

abstract: This paper is concerned with strong convergence of a tamed $\theta$-Euler-Maruyama scheme for neutral stochastic differential delay equations with superlinearly growing coefficients. We not only prove the strong convergence of implicit schemes, but also reveal the convergence rate for these equations driven by Brownian motion and pure jumps, respectively.


title: Discovering Conversational Dependencies between Messages in Dialogs

abstract: We investigate the task of inferring conversational dependencies between messages in one-on-one online chat, which has become one of the most popular forms of customer service. We propose a novel probabilistic classifier that leverages conversational, lexical and semantic information. The approach is evaluated empirically on a set of customer service chat logs from a Chinese e-commerce website. It outperforms heuristic baselines.


title: Interactive Prior Elicitation of Feature Similarities for Small Sample Size Prediction

abstract: Regression under the "small $n$, large $p$" conditions, of small sample size $n$ and large number of features $p$ in the learning data set, is a recurring setting in which learning from data is difficult. With prior knowledge about relationships of the features, $p$ can effectively be reduced, but explicating such prior knowledge is difficult for experts. In this paper we introduce a new method for eliciting expert prior knowledge about the similarity of the roles of features in the prediction task. The key idea is to use an interactive multidimensional-scaling (MDS) type scatterplot display of the features to elicit the similarity relationships, and then use the elicited relationships in the prior distribution of prediction parameters. Specifically, for learning to predict a target variable with Bayesian linear regression, the feature relationships are used to construct a Gaussian prior with a full covariance matrix for the regression coefficients. Evaluation of our method in experiments with simulated and real users on text data confirm that prior elicitation of feature similarities improves prediction accuracy. Furthermore, elicitation with an interactive scatterplot display outperforms straightforward elicitation where the users choose feature pairs from a feature list.


title: The Physical Systems Behind Optimization Algorithms

abstract: We use differential equations based approaches to provide some {\it \textbf{physics}} insights into analyzing the dynamics of popular optimization algorithms in machine learning. In particular, we study gradient descent, proximal gradient descent, coordinate gradient descent, proximal coordinate gradient, and Newton's methods as well as their Nesterov's accelerated variants in a unified framework motivated by a natural connection of optimization algorithms to physical systems. Our analysis is applicable to more general algorithms and optimization problems {\it \textbf{beyond}} convexity and strong convexity, e.g. Polyak-\L ojasiewicz and error bound conditions (possibly nonconvex).


title: Study of Dark-Matter Admixed Neutron Stars using the Equation of State from the Rotational Curves of Galaxies

abstract: In this work, we employ the dark matter equations of state (DMEOSs) obtained from the rotational curves of galaxies as well as the fermionic DMEOS with m = 1.0 GeV to study the structure of darkmatter admixed neutron stars (DMANSs). Applying the equation of state in the Skyrme framework for the neutron matter (NM), we calculate the mass-radius relation for different DMANSs with various DMEOSs and central pressure of dark matter (DM) to NM ratios. Our results show that for some DMEOSs, the mass-radius relations are in agreement with new observations, e.g. EXO 1745-248, 4U 1608-52, and 4U 1820-30, which are inconsistent with the normal neutron stars. We conclude that both DMEOSs and central pressure ratios of DM to NM affect the slope of the mass-radius relation of DMANSs. This is because of the interaction between DM and NM, which leads to gravitationally or self-bound DMANSs. We study the radius of the NM sphere as well as the radius of the DM halo for different DMANSs. The results confirm that, in some cases, a NM sphere with a small radius is surrounded by a halo of DM with a larger radius. Our calculations verify that, due to the different degrees of DM domination in DMANSs, with a value of the visible radius of a star two possible DMANSs with different masses can be exist. The gravitational redshift is also calculated for DMANSs with different DMEOSs and central pressure ratios. The results explain that the existence of DM in a DMANS leads to higher values of gravitational redshift of the star.


title: Exploring the efficacy and limitations of shock-cooling models: new analysis of Type II supernovae observed by the Kepler mission

abstract: Modern transient surveys have begun discovering and following supernovae (SNe) shortly after first light---providing systematic measurements of the rise of Type II SNe. We explore how analytic models of early shock-cooling emission from core-collapse SNe can constrain the progenitor's radius, explosion velocity, and local host extinction. We simulate synthetic photometry in several realistic observing scenarios and, assuming the models describe the typical explosions well, find that ultraviolet observations can constrain the progenitor's radius to a statistical uncertainty of $\pm10-15\%$, with a systematic uncertainty of $\pm 20\%$. With these observations the local host extinction ($A_V$) can be constrained to a factor of two and the shock velocity to $\pm 5\%$ with a systematic uncertainty of $\pm 10\%$. We also re-analyze the SN light curves presented in Garnavich et al. (2016) and find that KSN 2011a can be fit by a blue supergiant model with a progenitor radius of $R_s < 7.7 + 8.8 {\rm (stat)} + 1.9 {\rm (sys)}R_\odot$, while KSN 2011d can be fit with a red supergiant model with a progenitor radius of $R_s = 111_{-21 {\rm (stat)}}^{+89 {\rm (stat)}} \phantom{}_{-1{\rm (sys)}}^{+49 {\rm (sys)}}$. Our results do not agree with those of Garnavich et al. (2016). Moreover, we re-evaluate their claims and find that there is no statistically significant evidence for a shock breakout flare in the light curve of KSN 2011d.


title: Quantum autoencoders for efficient compression of quantum data

abstract: Classical autoencoders are neural networks that can learn efficient low dimensional representations of data in higher dimensional space. The task of an autoencoder is, given an input $x$, is to map $x$ to a lower dimensional point $y$ such that $x$ can likely be recovered from $y$. The structure of the underlying autoencoder network can be chosen to represent the data on a smaller dimension, effectively compressing the input. Inspired by this idea, we introduce the model of a quantum autoencoder to perform similar tasks on quantum data. The quantum autoencoder is trained to compress a particular dataset of quantum states, where a classical compression algorithm cannot be employed. The parameters of the quantum autoencoder are trained using classical optimization algorithms. We show an example of a simple programmable circuit that can be trained as an efficient autoencoder. We apply our model in the context of quantum simulation to compress ground states of the Hubbard model and molecular Hamiltonians.


title: Random versus maximum entropy models of neural population activity

abstract: The principle of maximum entropy provides a useful method for inferring statistical mechanics models from observations in correlated systems, and is widely used in a variety of fields where accurate data are available. While the assumptions underlying maximum entropy are intuitive and appealing, its adequacy for describing complex empirical data has been little studied in comparison to alternative approaches. Here data from the collective spiking activity of retinal neurons is reanalysed. The accuracy of the maximum entropy distribution constrained by mean firing rates and pairwise correlations is compared to a random ensemble of distributions constrained by the same observables. In general, maximum entropy approximates the true distribution better than the typical or mean distribution from that ensemble. This advantage improves with population size, with groups as small as 8 being almost always better described by maximum entropy. Failure of maximum entropy to outperform random models is found to be associated with strong correlations in the population.


title: 3D Shape Segmentation with Projective Convolutional Networks

abstract: This paper introduces a deep architecture for segmenting 3D objects into their labeled semantic parts. Our architecture combines image-based Fully Convolutional Networks (FCNs) and surface-based Conditional Random Fields (CRFs) to yield coherent segmentations of 3D shapes. The image-based FCNs are used for efficient view-based reasoning about 3D object parts. Through a special projection layer, FCN outputs are effectively aggregated across multiple views and scales, then are projected onto the 3D object surfaces. Finally, a surface-based CRF combines the projected outputs with geometric consistency cues to yield coherent segmentations. The whole architecture (multi-view FCNs and CRF) is trained end-to-end. Our approach significantly outperforms the existing state-of-the-art methods in the currently largest segmentation benchmark (ShapeNet). Finally, we demonstrate promising segmentation results on noisy 3D shapes acquired from consumer-grade depth cameras.


title: Uniform Atmospheric Retrieval Analysis of Ultracool Dwarfs II: Properties of 11 T-dwarfs

abstract: Brown dwarf spectra are rich in information revealing of the chemical and physical processes operating in their atmospheres. We apply a recently developed atmospheric retrieval tool to an ensemble of late T-dwarf (600-800K) near infrared spectra. With these spectra we are able to place direct constraints the molecular abundances of H$_2$O, CH$_4$, CO, CO$_2$, NH$_3$, H$_2$S, and Na+K, gravity, thermal structure (and effective temperature), photometric radius, and cloud optical depths. We find that ammonia, water, methane, and the alkali metals are present and well constrained in all 11 objects. From the abundance constraints we find no significant trend in the water, methane, or ammonia abundances with temperature, but find a very strong ($>$25$\sigma$) increasing trend in the alkali metal abundances with effective temperature, indicative of alkali rainout. We also find little evidence for optically thick clouds. With the methane and water abundances, we derive the intrinsic atmospheric metallicity and carbon-to-oxygen ratios. We find in our sample, that metallicities are typically sub solar and carbon-to-oxygen ratios are somewhat super solar, different than expectations from the local stellar population. We also find that the retrieved vertical thermal profiles are consistent with radiative equilibrium over the photospheric regions. Finally, we find that our retrieved effective temperatures are lower than previous inferences for some objects and that our radii are larger than expectations from evolutionary models, possibly indicative of un-resolved binaries. This investigation and methodology represents a paradigm in linking spectra to the determination of the fundamental chemical and physical processes governing cool brown dwarf atmospheres.


title: The BRAVE Program - I: Improved Bulge Stellar Velocity Dispersion Estimates for a Sample of Active Galaxies

abstract: We present new bulge stellar velocity dispersion measurements for 10 active galaxies with secure $M_{BH}$ determinations from reverberation-mapping. These new velocity dispersion measurements are based on spatially resolved kinematics from integral-field (IFU) spectroscopy. In all but one case, the field of view of the IFU extends beyond the effective radius of the galaxy, and in the case of Mrk 79 the field of view extends to almost one half the effective radius. This combination of spatial resolution and field of view allows for secure determinations of stellar velocity dispersion within the effective radius for all 10 target galaxies. Spatially resolved maps of the first (V) and second ($\sigma_{\star}$) moments of the line-of-sight velocity distribution (LOSVD) indicate the presence of kinematic substructure in most cases. In future projects we plan to explore methods of correcting for the effects of kinematic substructure in the derived bulge stellar velocity dispersion measurements.


title: Analysis of Multi-Index Monte Carlo Estimators for a Zakai SPDE

abstract: In this article, we propose a space-time Multi-Index Monte Carlo (MIMC) estimator for a one-dimensional parabolic stochastic partial differential equation (SPDE) of Zakai type. We compare the complexity with the Multilevel Monte Carlo (MLMC) method of Giles and Reisinger (2012), and find, by means of Fourier analysis, that the MIMC method: (i) has suboptimal complexity of $O(\varepsilon^{-2}|\log\varepsilon|^3)$ for a root mean square error (RMSE) $\varepsilon$ if the same spatial discretisation as in the MLMC method is used; (ii) has a better complexity of $O(\varepsilon^{-2}|\log\varepsilon|)$ if a carefully adapted discretisation is used; (iii) has to be adapted for non-smooth functionals. Numerical tests confirm these findings empirically.


title: Advances in using Internet searches to track dengue

abstract: Dengue is a mosquito-borne disease that threatens more than half of the world's population. Despite being endemic to over 100 countries, government-led efforts and mechanisms to timely identify and track the emergence of new infections are still lacking in many affected areas. Multiple methodologies that leverage the use of Internet-based data sources have been proposed as a way to complement dengue surveillance efforts. Among these, the trends in dengue-related Google searches have been shown to correlate with dengue activity. We extend a methodological framework, initially proposed and validated for flu surveillance, to produce near real-time estimates of dengue cases in five countries/regions: Mexico, Brazil, Thailand, Singapore and Taiwan. Our result shows that our modeling framework can be used to improve the tracking of dengue activity in multiple locations around the world.


title: Prospects of models with light sgoldstino in electron beam dump experiment at CERN SPS

abstract: We discuss phenomenology of light scalar sgoldstino in context of CERN electron beam dump experiment NA64. We calculate sgoldstino production rate for this experiment taking into account sgoldstino mixing with the Higgs boson and find a region in the model parameter space which can be tested in NA64.


title: Task-Guided and Path-Augmented Heterogeneous Network Embedding for Author Identification

abstract: In this paper, we study the problem of author identification under double-blind review setting, which is to identify potential authors given information of an anonymized paper. Different from existing approaches that rely heavily on feature engineering, we propose to use network embedding approach to address the problem, which can automatically represent nodes into lower dimensional feature vectors. However, there are two major limitations in recent studies on network embedding: (1) they are usually general-purpose embedding methods, which are independent of the specific tasks; and (2) most of these approaches can only deal with homogeneous networks, where the heterogeneity of the network is ignored. Hence, challenges faced here are two folds: (1) how to embed the network under the guidance of the author identification task, and (2) how to select the best type of information due to the heterogeneity of the network. To address the challenges, we propose a task-guided and path-augmented heterogeneous network embedding model. In our model, nodes are first embedded as vectors in latent feature space. Embeddings are then shared and jointly trained according to task-specific and network-general objectives. We extend the existing unsupervised network embedding to incorporate meta paths in heterogeneous networks, and select paths according to the specific task. The guidance from author identification task for network embedding is provided both explicitly in joint training and implicitly during meta path selection. Our experiments demonstrate that by using path-augmented network embedding with task guidance, our model can obtain significantly better accuracy at identifying the true authors comparing to existing methods.


title: Recalibration of the $M_{BH}-\sigma_{\star}$ Relation for AGN

abstract: We present a re-calibration of the $M_{BH}-\sigma_{\star}$ relation, based on a sample of 16 reverberation-mapped galaxies with newly determined bulge stellar velocity dispersions ($\sigma_{\star}$) from integral-field spectroscopy (IFS), and a sample of 32 quiescent galaxies with publicly available IFS. For both samples, $\sigma_{\star}$ is determined via two different methods that are popular in the literature, and we provide fits for each sample based on both sets of $\sigma_{\star}$. We find the fit to the AGN sample is shallower than the fit to the quiescent galaxy sample, and that the slopes for each sample are in agreement with previous investigations. However, the intercepts to the quiescent galaxy relations are notably higher than those found in previous studies, due to the systematically lower $\sigma_{\star}$ measurements that we obtain from IFS. We find that this may be driven, in part, by poorly constrained measurements of bulge effective radius ($r_{e}$) for the quiescent galaxy sample, which may bias the $\sigma_{\star}$ measurements low. We use these quiescent galaxy parameterizations, as well as one from the literature, to recalculate the virial scaling factor $f$. We assess the potential biases in each measurement, and suggest $f=4.82\pm1.67$ as the best currently available estimate. However, we caution that the details of how $\sigma_{\star}$ is measured can significantly affect $f$, and there is still much room for improvement.


title: A Generalization of the Curry-Howard Correspondence

abstract: We present a variant of the calculus of deductive systems developed in (Lambek 1972, 1974), and give a generalization of the Curry-Howard-Lambek theorem giving an equivalence between the category of typed lambda-calculi and the category of cartesian closed categories and exponential-preserving morphisms that leverages the theory of generalized categories (Schoenbaum 2016). We discuss potential applications and extensions.


title: Constraints On The Dynamical Environments Of Supermassive Black-hole Binaries Using Pulsar-timing Arrays

abstract: We introduce a technique for gravitational-wave analysis, where Gaussian process regression is used to emulate the strain spectrum of a stochastic background using population-synthesis simulations. This leads to direct Bayesian inference on astrophysical parameters. For PTAs specifically, we interpolate over the parameter space of supermassive black-hole binary environments, including 3-body stellar scattering, and evolving orbital eccentricity. We illustrate our approach on mock data, and assess the prospects for inference with data similar to the NANOGrav 9-yr data release.


title: Mapping the Extinction Curve in 3D: Structure on Kiloparsec Scales

abstract: Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made possible precise measurements of the shape of the extinction curve for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dimensional maps of dust column density with R(V) measurements to constrain the three-dimensional distribution of R(V) in the Milky Way. We find that variations in R(V) are correlated on kiloparsec scales. In particular, most of the dust within one kiloparsec in the outer Galaxy, including many local molecular clouds (Orion, Taurus, Perseus, California, Cepheus), has a significantly lower R(V) than more distant dust in the Milky Way. These results provide new input to models of dust evolution and processing, and complicate application of locally derived extinction curves to more distant regions of the Milky Way and to other galaxies.


title: The spinning Kerr-black-hole-mirror bomb: A lower bound on the radius of the reflecting mirror

abstract: The intriguing superradiant amplification phenomenon allows an orbiting scalar field to extract rotational energy from a spinning Kerr black hole. Interestingly, the energy extraction rate can grow exponentially in time if the black-hole-field system is placed inside a reflecting mirror which prevents the field from radiating its energy to infinity. This composed Kerr-black-hole-scalar-field-mirror system, first designed by Press and Teukolsky, has attracted the attention of physicists over the last four decades. Previous numerical studies of this spinning {\it black-hole bomb} have revealed the interesting fact that the superradiant instability shuts down if the reflecting mirror is placed too close to the black-hole horizon. In the present study we use analytical techniques to explore the superradiant instability regime of this composed Kerr-black-hole-linearized-scalar-field-mirror system. In particular, it is proved that the lower bound ${{r_{\text{m}}}\over{r_+}}>{1\over 2}\Big(\sqrt{1+{{8M}\over{r_-}}}-1\Big)$ provides a necessary condition for the development of the exponentially growing superradiant instabilities in this composed physical system, where $r_{\text{m}}$ is the radius of the confining mirror and $r_{\pm}$ are the horizon radii of the spinning Kerr black hole. We further show that, in the linearized regime, this {\it analytically} derived lower bound on the radius of the confining mirror agrees with direct {\it numerical} computations of the superradiant instability spectrum which characterizes the spinning black-hole-mirror bomb.


title: Anomalous screening of quantum impurities by a neutral environment

abstract: It is a common knowledge that an effective interaction of a quantum impurity with an electromagnetic field can be screened by surrounding charge carriers, whether mobile or static. Here we demonstrate that very strong, `anomalous' screening can take place in the presence of a neutral, weakly-polarizable environment, due to an exchange of orbital angular momentum between the impurity and the bath. Furthermore, we show that it is possible to generalize all phenomena related to isolated impurities in an external field to the case when a many-body environment is present, by casting the problem in terms of the angulon quasiparticle. As a result, the relevant observables such as the effective Rabi frequency, geometric phase, and impurity spatial alignment are straightforward to evaluate in terms of a single parameter: the angular-momentum-dependent screening factor.

categories: cond-mat.mes-hall cond-mat.quant-gas physics.atom-ph physics.chem-ph


title: Evolution of the real-space correlation function from next generation cluster surveys

abstract: We investigate to which accuracy it is possible to recover the real-space two-point correlation function of galaxy clusters from cluster catalogues based on photometric redshifts, and test our ability to measure the redshift and mass evolution of the correlation length and the bias parameter as a function of the redshift uncertainty. We calculate the correlation function for cluster sub-samples covering various mass and redshift bins selected from a light-cone catalogue. To simulate the distribution of clusters in photometric redshift space, we assign to each cluster a redshift randomly extracted from a Gaussian distribution. The dispersion is varied in the range $\sigma_{(z=0)} = 0.001$ to $0.050$. The correlation function in real-space is computed through estimation and deprojection of $w_{p}(r_{p})$. Four mass ranges (from $M_{halo}> 2 \times 10^{13}$ to $M_{halo}> 2 \times 10^{14}$) and six redshift slices covering the redshift range [0,2] are investigated, using cosmological redshifts and photo-z configurations. We find a clear increase of the correlation amplitude as a function of redshift and mass for the $z_{c}$ samples. The evolution of the derived bias parameter is in agreement with theoretical expectations. From our pilot sample limited to $M_{halo}> 5 \times 10^{13} (0.4 < z < 0.7)$, we find that the real-space correlation function can be recovered by deprojection of $w_{p}(r_{p})$ within an accuracy of 5% for $\sigma_{z} = 0.001 \times (1 + z_{c})$ and within 10% for $\sigma_{z} = 0.03 \times (1 + z_{c})$. The evolution of the correlation in redshift and mass is clearly detected for all $\sigma_{z}$ tested. The best-fit parameters $(r_{0}$ and ${\gamma})$ as well as the bias obtained from the deprojection method for all $\sigma_{z}$ are within the $1 \sigma$ uncertainty of the $z_{c}$ sample.


title: The 5 Gyr evolution of sub-M* galaxies

abstract: We gathered two complete samples of $M_{AB}(r)<-18$ ($M_{star} > 10^{9} M_{\odot}$) galaxies, which are representative of the present-day galaxies and their counterparts at 5 Gyr ago. We analysed their 2D luminosity profiles and carefully decomposed them into bulges, bars and discs. This was done in a very consistent way at the two epochs, by using the same image quality and same (red) filters at rest. We classified them into elliptical, lenticular, spiral, and peculiar galaxies on the basis of a morphological decision tree. We found that at $z=0$, sub-M* ($10^{9} M_{\odot} < M_{star} < 1.5 \times 10^{10} M_{\odot}$) galaxies follow a similar Hubble Sequence compared to their massive counterparts, though with a considerable larger number of (1) peculiar galaxies and (2) low surface brightness galaxies. These trends persist in the $z\sim0.5$ sample, suggesting that sub-M* galaxies have not reached yet a virialised state, conversely to their more massive counterparts. The fraction of peculiar galaxies is always high, consistent with a hierarchical scenario in which minor mergers could have played a more important role for sub-M* galaxies than for more massive galaxies. Interestingly, we also discovered that more than 10\% of the sub-M* galaxies at z=0.5 are low surface brightness galaxies with clumpy and perturbed features that suggest merging. Even more enigmatic is the fact that their disc scale-lengths are comparable to that of M31, while their stellar masses are similar to that of the LMC.


title: Relic galaxies: where are they?

abstract: The finding that massive galaxies grow with cosmic time fired the starting gun for the search of objects which could have survived up to the present day without suffering substantial changes (neither in their structures, neither in their stellar populations). Nevertheless, and despite the community efforts, up to now only one firm candidate to be considered one of these relics is known: NGC 1277. Curiously, this galaxy is located at the centre of one of the most rich near galaxy clusters: Perseus. Is its location a matter of chance? Should relic hunters focus their search on galaxy clusters? In order to reply this question, we have performed a simultaneous and analogous analysis using simulations (Millennium I-WMAP7) and observations (New York University Value-Added Galaxy Catalogue). Our results in both frameworks agree: it is more probable to find relics in high density environments.


title: A Search for L/T Transition Dwarfs With Pan-STARRS1 and WISE. II. L/T Transition Atmospheres and Young Discoveries

abstract: The evolution of brown dwarfs from L to T spectral types is one of the least understood aspects of the ultracool population, partly for lack of a large, well-defined, and well-characterized sample in the L/T transition. To improve the existing census, we have searched $\approx$28,000 deg$^2$ using the Pan-STARRS1 and WISE surveys for L/T transition dwarfs within 25 pc. We present 130 ultracool dwarf discoveries with estimated distances $\approx9-130$ pc, including 21 that were independently discovered by other authors and 3 that were previously identified as photometric candidates. Seventy-nine of our objects have near-IR spectral types of L6-T4.5, the most L/T transition dwarfs from any search to date, and we have increased the census of L9-T1.5 objects within 25 pc by over 50%. The color distribution of our discoveries provides further evidence for the "L/T gap," a deficit of objects with $(J-K)_{\rm MKO}\approx0.0-0.5$ mag in the L/T transition, and thus reinforces the idea that the transition from cloudy to clear photospheres occurs rapidly. Among our discoveries are 31 candidate binaries based on their low-resolution spectral features. Two of these candidates are common proper motion companions to nearby main sequence stars; if confirmed as binaries, these would be rare benchmark systems with the potential to stringently test ultracool evolutionary models. Our search also serendipitously identified 23 late-M and L dwarfs with spectroscopic signs of low gravity implying youth. Finally, we identify 10 candidate members of nearby young moving groups (YMG) with spectral types L7-T4.5, including three showing spectroscopic signs of low gravity. If confirmed, any of these would be among the coolest known YMG members and would help to determine the effective temperature at which young brown dwarfs cross the L/T transition. (Abridged)


title: A Sommerfeld Toolbox for Colored Dark Sectors

abstract: We present analytical formulae for the Sommerfeld corrections to the annihilation of massive colored particles into quarks and gluons through the strong interaction. These corrections are essential to accurately compute the dark matter relic density for coannihilation with colored partners. Our formulae allow to compute the Sommerfeld effect, not only for the lowest term in the angular momentum expansion of the amplitude, but for all orders in the partial wave expansion. In particular, we carefully account for the effects of the spin of the annihilating particle on the symmetry of the two-particle wave function. This work focuses on strongly interacting particles of arbitrary spin in the triplet, sextet and octet color representations. For typical velocities during freeze-out, we find that including Sommerfeld corrections on the next-to-leading order partial wave leads to modifications of up to ten to twenty percent on the total annihilation cross section. Complementary to QCD, we generalize our results to particles charged under an arbitrary unbroken SU(N) gauge group, as encountered in dark glueball models. In connection with this paper a Mathematica notebook is provided to compute the Sommerfeld corrections for colored particles up to arbitrary order in the angular momentum expansion.


title: Probing the assembly history and dynamical evolution of massive black hole binaries with pulsar timing arrays

abstract: We consider the inverse problem in pulsar timing array (PTA) analysis, investigating what astrophysical information about the underlying massive black hole binary (MBHB) population can be recovered from the detection of a stochastic gravitational wave background (GWB). We employ a physically motivated model that connects the GWB spectrum to a series of parameters describing the underlying redshift evolution of the MBHB mass function and to the typical eccentricity they acquire while interacting with the dense environment of post merger galactic nuclei. This allows the folding in of information about the spectral shape of the GWB into the analysis. The priors on the model parameters are assumed to be uninformative and consistent with the current lack of secure observations of sub-parsec MBHBs. We explore the implications of current upper limits as well as of future detections with a variety of PTA configurations. We confirm our previous finding that current upper limits can only place an upper bound on the overall MBHB merger rate. Depending on the properties of the array, future detections can also constrain several MBHB population models at different degrees of fidelity. In particular, a simultaneous detection of a steepening of the spectrum at high frequency and a bending at low frequency will place strong constraints on both the MBHB mass function and on the typical eccentricity of inspiralling MBHBs, providing insights on MBHB astrophysics unlikely to be achievable by any other means.


title: Velocity Segregation and Systematic Biases In Velocity Dispersion Estimates With the SPT-GMOS Spectroscopic Survey

abstract: The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplemented by additional published spectroscopy, resulting in a final spectroscopic sample of 4148 galaxy spectra---2868 cluster members. The velocity dispersion of star-forming cluster galaxies is $17\pm4$% greater than that of passive cluster galaxies, and the velocity dispersion of bright ($m < m^{*}-0.5$) cluster galaxies is $11\pm4$% lower than the velocity dispersion of our total member population. We find good agreement with simulations regarding the shape of the relationship between the measured velocity dispersion and the fraction of passive vs. star-forming galaxies used to measure it, but we find a small offset between this relationship as measured in data and simulations in which suggests that our dispersions are systematically low by as much as 3\% relative to simulations. We argue that this offset could be interpreted as a measurement of the effective velocity bias that describes the ratio of our observed velocity dispersions and the intrinsic velocity dispersion of dark matter particles in a published simulation result. Measuring velocity bias in this way suggests that large spectroscopic surveys can improve dispersion-based mass-observable scaling relations for cosmology even in the face of velocity biases, by quantifying and ultimately calibrating them out.


title: Ray-tracing and polarized radiative transfer in General Relativity

abstract: We discuss the problem of polarized radiative transfer in general relativity. We present a set of equations suitable for solving the problem numerically for the case of an arbitrary space-time metric, and show numerical solutions to example problems. The solutions are computed with a new ray-tracing code, Arcmancer, developed by the authors.