title: Exploring cosmic origins with CORE: effects of observer peculiar motion

abstract: We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.


title: Formation of Vacancies in Si- and Ge-based Clathrates: Role of Electron Localization and Symmetry Breaking

abstract: The formation of framework vacancies in Si- and Ge-based type-I clathrates is studied as function of filling the cages with K and Ba atoms using density-functional theory. Our analysis reveals the relevance of structural disorder, geometric relaxation, electronic saturation, as well as vibrational and configurational entropy. In the Si clathrates we find that vacancies are unstable, but very differently, in Ge clathrates up to three vacancies per unit cell can be stabilized. This contrasting behavior is largely driven by the different energy gain on populating the electronic vacancy states, which originates from the different degree of localization of the valence orbitals of Si and Ge. This also actuates a qualitatively different atomic relaxation of the framework.


title: A Concept For Cancelling The Leakage Field Inside The Stored Beam Chamber Of A Septum Magnet

abstract: The Advanced Photon Source is in the process of up-grading its storage ring from a double-bend to a multi-bend lattice as part of the APS Upgrade Project (APS-U). A swap-out injection scheme is planned for the APS-U to keep a constant beam current and to enable a small dynamic aperture. A novel concept that cancels out the effect of leakage field inside the stored beam chamber was introduced in the design of the septum magnet. As a result, the horizontal deflecting angle of the stored beam was reduced to below 1 micro-rad with a 2 mm septum thick-ness and 1.06T normal injection field. The concept helped to minimize the integrated skew quadrupole field and normal sextupole fields inside stored beam chamber as well. The designed septum magnet deflects the injected electron beam by 89 mrad with a ring energy of 6 GeV. The stored beam chamber has an 8 mm x 6 mm super-ellipsoidal aperture. The magnet is straight; however, it is tilted in yaw, roll, and pitch from the stored beam cham-ber to meet the on axis swap out injection requirements for the APS-U lattice.


title: Spatio-temporal analysis of regional unemployment rates: A comparison of model based approaches

abstract: This study aims to analyze the methodologies that can be used to estimate the total number of unemployed, as well as the unemployment rates for 28 regions of Portugal, designated as NUTS III regions, using model based approaches as compared to the direct estimation methods currently employed by INE (National Statistical Institute of Portugal). Model based methods, often known as small area estimation methods (Rao, 2003), "borrow strength" from neighbouring regions and in doing so, aim to compensate for the small sample sizes often observed in these areas. Consequently, it is generally accepted that model based methods tend to produce estimates which have lesser variation. Other benefit in employing model based methods is the possibility of including auxiliary information in the form of variables of interest and latent random structures. This study focuses on the application of Bayesian hierarchical models to the Portuguese Labor Force Survey data from the 1st quarter of 2011 to the 4th quarter of 2013. Three different data modeling strategies are considered and compared: Modeling of the total unemployed through Poisson, Binomial and Negative Binomial models; modeling of rates using a Beta model; and modeling of the three states of the labor market (employed, unemployed and inactive) by a Multinomial model. The implementation of these models is based on the \textit{Integrated Nested Laplace Approximation} (INLA) approach, except for the Multinomial model which is implemented based on the method of Monte Carlo Markov Chain (MCMC). Finally, a comparison of the performance of these models, as well as the comparison of the results with those obtained by direct estimation methods at NUTS III level are given.


title: The Luminous Blue Variable RMC127 as seen with ALMA and ATCA

abstract: We present ALMA and ATCA observations of the luminous blue variable \rmc. The radio maps show for the first time the core of the nebula and evidence that the nebula is strongly asymmetric with a Z-pattern shape. Hints of this morphology are also visible in the archival \emph{HST} $\rm H\alpha$ image, which overall resembles the radio emission. The emission mechanism in the outer nebula is optically thin free-free in the radio. At high frequencies, a component of point-source emission appears at the position of the star, up to the ALMA frequencies. The rising flux density distribution ($S_{\nu}\sim \nu^{0.78\pm0.05}$) of this object suggests thermal emission from the ionized stellar wind and indicates a departure from spherical symmetry with $n_{e}(r)\propto r^{-2}$. We examine different scenarios to explain this excess of thermal emission from the wind and show that this can arise from a bipolar outflow, supporting the suggestion by other authors that the stellar wind of \rmc is aspherical. We fit the data with two collimated ionized wind models and we find that the mass-loss rate can be a factor of two or more smaller than in the spherical case. We also fit the photometry obtained by IR space telescopes and deduce that the mid- to far-IR emission must arise from extended, cool ($\sim80\,\rm K$) dust within the outer ionized nebula. Finally we discuss two possible scenarios for the nebular morphology: the canonical single star expanding shell geometry, and a precessing jet model assuming presence of a companion star.


title: Probability of loops formation in star polymers in long range correlated disorder

abstract: We analyze the statistics of loops formation in $f$-branched star polymers in an environment with structural defects, correlated at large distances $r$ according to a power law $\sim r^{-a}$. Applying the direct polymer renormalization approach, we found the values of the set of universal exponents, governing the scaling of probabilities of various types of loops in macromolecules.


title: The real projective spaces in homotopy type theory

abstract: Homotopy type theory is a version of Martin-L\"of type theory taking advantage of its homotopical models. In particular, we can use and construct objects of homotopy theory and reason about them using higher inductive types. In this article, we construct the real projective spaces, key players in homotopy theory, as certain higher inductive types in homotopy type theory. The classical definition of RP(n), as the quotient space identifying antipodal points of the n-sphere, does not translate directly to homotopy type theory. Instead, we define RP(n) by induction on n simultaneously with its tautological bundle of 2-element sets. As the base case, we take RP(-1) to be the empty type. In the inductive step, we take RP(n+1) to be the mapping cone of the projection map of the tautological bundle of RP(n), and we use its universal property and the univalence axiom to define the tautological bundle on RP(n+1). By showing that the total space of the tautological bundle of RP(n) is the n-sphere, we retrieve the classical description of RP(n+1) as RP(n) with an (n+1)-cell attached to it. The infinite dimensional real projective space, defined as the sequential colimit of the RP(n) with the canonical inclusion maps, is equivalent to the Eilenberg-MacLane space K(Z/2Z,1), which here arises as the subtype of the universe consisting of 2-element types. Indeed, the infinite dimensional projective space classifies the 0-sphere bundles, which one can think of as synthetic line bundles. These constructions in homotopy type theory further illustrate the utility of homotopy type theory, including the interplay of type theoretic and homotopy theoretic ideas.


title: Polar factorization of conformal and projective maps of the sphere in the sense of optimal mass transport

abstract: Let M be a compact Riemannian manifold and let $\mu$,d be the associated measure and distance on M. Robert McCann obtained, generalizing results for the Euclidean case by Yann Brenier, the polar factorization of Borel maps S : M -> M pushing forward $\mu$ to a measure $\nu$: each S factors uniquely a.e. into the composition S = T \circ U, where U : M -> M is volume preserving and T : M -> M is the optimal map transporting $\mu$ to $\nu$ with respect to the cost function d^2/2. In this article we study the polar factorization of conformal and projective maps of the sphere S^n. For conformal maps, which may be identified with elements of the identity component of O(1,n+1), we prove that the polar factorization in the sense of optimal mass transport coincides with the algebraic polar factorization (Cartan decomposition) of this Lie group. For the projective case, where the group GL_+(n+1) is involved, we find necessary and sufficient conditions for these two factorizations to agree.


title: A power-law coupled three-form dark energy model

abstract: We consider a field theory model of coupled dark energy which treats dark energy as a three-form field and dark matter as a spinor field. By assuming the effective mass of dark matter as a power-law function of the three-form field and neglecting the potential term of dark energy, we obtain three solutions of the autonomous system of evolution equations, including a de Sitter attractor, a tracking solution and an approximate solution. To understand the strength of the coupling, we confront the model with the latest Type Ia Supernova (SN Ia), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave Backround (CMB) radiation observations, with the conclusion that the combination of these three databases marginalized over the present dark matter density parameter $\Omega_{m0}$ and the present three-form field $\kappa X_{0}$ gives stringent constraints on the coupling constant, $-0.017< \lambda <0.047$ ($2\sigma$ confidence level), by which we give out the model applicable parameter range.


title: Derivation of the Asymptotic Eigenvalue Distribution for Causal 2D-AR Models under Upscaling

abstract: This technical report describes the derivation of the asymptotic eigenvalue distribution for causal 2D-AR models under an upscaling scenario. Specifically, it tackles the analytical derivation of the asymptotic eigenvalue distribution of the sample autocorrelation matrix corresponding to genuine and upscaled images. It also includes the pseudocode of the derived approaches for resampling detection and resampling factor estimation that are based on this analysis.


title: On the large-Q^2 behavior of the pion transition form factor

abstract: We study the transition of non-perturbative to perturbative QCD in situations with possible violations of scaling limits. To this end we consider the singly- and doubly-virtual pion transition form factor $\pi^0\to\gamma\gamma$ at all momentum scales of symmetric and asymmetric photon momenta within the Dyson-Schwinger/Bethe-Salpeter approach. For the doubly virtual form factor we find good agreement with perturbative asymptotic scaling laws. For the singly-virtual form factor our results agree with the Belle data. At very large off-shell photon momenta we identify a mechanism that introduces quantitative modifications to Efremov-Radyushkin-Brodsky-Lepage scaling.


title: Deep Occlusion Reasoning for Multi-Camera Multi-Target Detection

abstract: People detection in single 2D images has improved greatly in recent years. However, comparatively little of this progress has percolated into multi-camera multi-people tracking algorithms, whose performance still degrades severely when scenes become very crowded. In this work, we introduce a new architecture that combines Convolutional Neural Nets and Conditional Random Fields to explicitly model those ambiguities. One of its key ingredients are high-order CRF terms that model potential occlusions and give our approach its robustness even when many people are present. Our model is trained end-to-end and we show that it outperforms several state-of-art algorithms on challenging scenes.


title: Accurate Single Stage Detector Using Recurrent Rolling Convolution

abstract: Most of the recent successful methods in accurate object detection and localization used some variants of R-CNN style two stage Convolutional Neural Networks (CNN) where plausible regions were proposed in the first stage then followed by a second stage for decision refinement. Despite the simplicity of training and the efficiency in deployment, the single stage detection methods have not been as competitive when evaluated in benchmarks consider mAP for high IoU thresholds. In this paper, we proposed a novel single stage end-to-end trainable object detection network to overcome this limitation. We achieved this by introducing Recurrent Rolling Convolution (RRC) architecture over multi-scale feature maps to construct object classifiers and bounding box regressors which are "deep in context". We evaluated our method in the challenging KITTI dataset which measures methods under IoU threshold of 0.7. We showed that with RRC, a single reduced VGG-16 based model already significantly outperformed all the previously published results. At the time this paper was written our models ranked the first in KITTI car detection (the hard level), the first in cyclist detection and the second in pedestrian detection. These results were not reached by the previous single stage methods. The code is publicly available.


title: The 2015-2016 outburst of the classical EXor V1118 Ori

abstract: After a quiescence period of about 10 years, the classical EXor source V1118 Ori has undergone an accretion outburst in 2015 September. The maximum brightness (DV > 4 mag) was reached in 2015 December and was maintained for several months. Since 2016 September, the source is in a declining phase. Photometry and low/ high-resolution spectroscopy were obtained with MODS and LUCI2 at the {\it Large Binocular Telescope}, with the facilities at the Asiago 1.22 and 1.82 m telescopes, and with GIANO at the {\it Telescopio Nazionale Galileo}. The spectra are dominated by emission lines of \hi\ and neutral metallic species. From line and continuum analysis we derive the mass accretion rate and its evolution during the outburst. Considering that extinction may vary between 1.5 and 2.9 mag, we obtain m_acc= 0.3$-$2.0 10$^{-8}$ m_sun/yr, in quiescence and m_acc= 0.2$-$1.9 10$^{-6}$ m_sun/yr, at the outburst peak. The Balmer decrement shape has been interpreted by means of line excitation models, finding that from quiescence to outburst peak, the electron density has increased from $\sim$ 2 10$^9$ cm$^{-3}$ to $\sim$ 4 10$^{11}$ cm$^{-3}$. The profiles of the metallic lines are symmetric and narrower than 100 km s$^{-1}$, while \hi\, and \hei\,\,lines show prominent wings extending up to $\pm$ 500 km s$^{-1}$. The metallic lines likely originate at the base of the accretion columns, where neutrals are efficiently shielded against the ionizing photons, while faster ionized gas is closer to the star. Outflowing activity is testified by the detection of a variable P Cyg-like profile of the H$\alpha$ and \hei\, 1.08\,$\mu$m lines.


title: Two-mode generalization of the Jaynes-Cummings and Anti-Jaynes-Cummings models

abstract: We introduce two generalizations of the Jaynes-Cummings (JC) model for two modes of oscillation. The first model is formed by two Jaynes-Cummings interactions, while the second model is written as a simultaneous Jaynes-Cummings and Anti-Jaynes-Cummings (AJC) interactions. We study some of its properties and obtain the energy spectrum and eigenfunctions of these models by using the tilting transformation and the Perelomov number coherent states of the two-dimensional harmonic oscillator. Moreover, as physical applications, we connect these new models with two important and novelty problems: The relativistic non-degenerate parametric amplifier and the relativistic problem of two coupled oscillators.


title: A Catalan Subset of Descending Plane Partitions

abstract: Descending plane partitions, alternating sign matrices, and totally symmetric self-complementary plane partitions are equinumerous combinatorial sets for which no explicit bijection is known. In this paper, we isolate a subset of descending plane partitions counted by the Catalan numbers. The proof follows by constructing a generating tree on these descending plane partitions that has the same structure as the generating tree for 231-avoiding permutations. We hope this result will provide insight on the search for a bijection with alternating sign matrices and/or totally symmetric self-complementary plane partitions, since these also contain Catalan subsets.


title: A remark on projective limits of function spaces

abstract: Let Omega subset of C^d be an open set and Km, m = 1, 2, . . . an exhaustion of Omega by compact subsets of Omega. We set Omega_m = int(Km) and let Xm(Omega_m) be a topological space of holomorphic functions on Omega_m between A^ infinity (Omega_m) and H(Omega_m). Then we show that the projective limit of the family Xm(Omega_m), m = 1, 2, . . . , under the restriction maps is homeomorphic and linearly isomorphic to the Frechet space H(Omega), idependently of the choice of the spaces Xm(Omega_m).


title: Redefining Context Windows for Word Embedding Models: An Experimental Study

abstract: Distributional semantic models learn vector representations of words through the contexts they occur in. Although the choice of context (which often takes the form of a sliding window) has a direct influence on the resulting embeddings, the exact role of this model component is still not fully understood. This paper presents a systematic analysis of context windows based on a set of four distinct hyper-parameters. We train continuous Skip-Gram models on two English-language corpora for various combinations of these hyper-parameters, and evaluate them on both lexical similarity and analogy tasks. Notable experimental results are the positive impact of cross-sentential contexts and the surprisingly good performance of right-context windows.


title: Positive Semidefiniteness and Positive Definiteness of a Linear Parametric Interval Matrix

abstract: We consider a symmetric matrix, the entries of which depend linearly on some parameters. The domains of the parameters are compact real intervals. We investigate the problem of checking whether for each (or some) setting of the parameters, the matrix is positive definite (or positive semidefinite). We state a characterization in the form of equivalent conditions, and also propose some computationally cheap sufficient\,/\,necessary conditions. Our results extend the classical results on positive (semi-)definiteness of interval matrices. They may be useful for checking convexity or non-convexity in global optimization methods based on branch and bound framework and using interval techniques.


title: Automation of NLO QCD and EW corrections with Sherpa and Recola

abstract: This publication presents the combination of the one-loop matrix-element generator Recola with the multipurpose Monte Carlo program Sherpa. Since both programs are highly automated, the resulting Sherpa+Recola framework allows for the computation of -in principle- any Standard Model process at both NLO QCD and EW accuracy. To illustrate this, three representative LHC processes have been computed at NLO QCD and EW: vector-boson production in association with jets, off-shell Z-boson pair production, and the production of a top-quark pair in association with a Higgs boson. In addition to fixed-order computations, when considering QCD corrections, all functionalities of Sherpa, i.e. particle decays, QCD parton showers, hadronisation, underlying events, etc. can be used in combination with Recola. This is demonstrated by the merging and matching of one-loop QCD matrix elements for Drell-Yan production in association with jets to the parton shower. The implementation is fully automatised, thus making it a perfect tool for both experimentalists and theorists who want to use state-of-the-art predictions at NLO accuracy.


title: Continuous gravitational wave searches with pulsar timing arrays: Maximization versus marginalization over pulsar phase parameters

abstract: Resolvable Supermassive Black Hole Binaries are promising sources for Pulsar Timing Array based gravitational wave searches. Search algorithms for such targets must contend with the large number of so-called pulsar phase parameters in the joint log-likelihood function of the data. We compare the localization accuracy for two approaches: Maximization over the pulsar phase parameters (MaxPhase) against marginalization over them (AvPhase). Using simulated data from a pulsar timing array with 17 pulsars, we find that for weak and moderately strong signals, AvPhase outperforms MaxPhase significantly, while they perform comparably for strong signals.


title: Thermodynamic and Holographic Information Dual to Volume

abstract: In this paper, we will analyze the connection between the fidelity susceptibility, the holographic complexity and the thermodynamic volume. We will regularize the fidelity susceptibility and the holographic complexity by subtracting the contribution of the background AdS spacetime from the deformation of the AdS spacetime. It will be demonstrated that this regularized fidelity susceptibility has the same behavior as the thermodynamic volume and that the regularized complexity has a very different behavior. As the information dual to different volumes in the bulk would be measured by the fidelity susceptibility and the holographic complexity, this paper will establish a connection between thermodynamics and information dual to a volume.


title: Importance Sampled Stochastic Optimization for Variational Inference

abstract: Variational inference approximates the posterior distribution of a probabilistic model with a parameterized density by maximizing a lower bound for the model evidence. Modern solutions fit a flexible approximation with stochastic gradient descent, using Monte Carlo approximation for the gradients. This enables variational inference for arbitrary differentiable probabilistic models, and consequently makes variational inference feasible for probabilistic programming languages. In this work we develop more efficient inference algorithms for the task by considering importance sampling estimates for the gradients. We show how the gradient with respect to the approximation parameters can often be evaluated efficiently without needing to re-compute gradients of the model itself, and then proceed to derive practical algorithms that use importance sampled estimates to speed up computation.We present importance sampled stochastic gradient descent that outperforms standard stochastic gradient descent by a clear margin for a range of models, and provide a justifiable variant of stochastic average gradients for variational inference.


title: On the Statistical Foundations of Kaluza's Magnetohydrodynamics

abstract: The introduction of electromagnetic fields into the Boltzmann equation following a 5D general relativistic approach is considered in order to establish the transport equations for dilute charged fluids in the presence of a weak electromagnetic field. The conserved 5D stress-energy tensor is evaluated using the J\"uttner function for non-degenerate relativistic gases in local equilibrium, and the evolution equations for the local thermodynamic variables are established by means of relativistic kinetic theory. An outline of the possibilities offered by the Kaluza-type approach to MHD is also included.


title: Tests for CPT sum rule and U-spin violation in Time-dependent CP violation of $B^0_s \to K^+ K^-$ and $B^0_d \to \pi^+ \pi^-$

abstract: Recent LHCb data for time-dependent CP violation in $B_d^0 \to \pi^+\pi^-$ and $B^0_s\to K^+K^-$ show deviations from theoretical predictions. Besides their central values for $\mathcal C_{K^+K^-}$, $\mathcal S_{K^+K^-}$ and $\mathcal A^{\Delta \Gamma}_{K^+K^-}$ violate quantum mechanic CPT invariant sum rule (CPT sum rule) prediction of $|\mathcal C_{K^+K^-}|^2 + |\mathcal S_{K^+K^-}| ^2 + |\mathcal A^{\Delta \Gamma}_{K^+K^-}|^2 = 1$ (LHCb data imply the sum to be $0.67\pm 0.20$.), their values for $\mathcal C_{K^+ K^-}= 0.24\pm 0.06\pm {0.02}$ and $\mathcal C_{\pi^+ \pi^-} = - 0.24\pm 0.07\pm 0.01$ also show large violation of SU(3) or its U-spin sub-group symmetry (SU(3)/U) relation $\mathcal C_{K^+ K^-} /\mathcal C_{\pi^- \pi^+} = - \mathcal B(B_d^0 \to \pi^- \pi^+)\tau_{B^0_s}/\mathcal B(B^0_s \to K^+ K^-)\tau_{B_d^0}$ (LHCb data imply the ratio of left-side to right-side to be $4.67\pm 1.88$.) . The LHCb results need to be further confirmed to be taken seriously. We suggest to use time-dependent CP violation in $B_s\to K^0\bar K^0, \pi^+\pi^-, \pi^0\pi^0$ to further test the CPT sum rule. Assuming that the sum rule holds, we propose that violation of the SU(3)/U relation may indicate a large FSI phase difference in the $\pi^+\pi^-$ and $K^+K^-$ re-scattering. We suggest several other U-spin pairs of $B\to PP$ decays to further test SU(3)/U relations.


title: When black holes collide: Probing the interior composition by the spectrum of ringdown modes and emitted gravitational waves

abstract: The merger of colliding black holes (BHs) should lead to the production of ringdown or quasinormal modes (QNMs), which may very well be sensitive to the state of the interior. We put this idea to the test with a recent proposal that the interior of a BH consists of a bound state of highly excited, long, closed, interacting strings; figuratively, a collapsed polymer. We show that such BHs do indeed have a distinct signature in their QNM spectrum: A new class of modes whose frequencies are parametrically lower than the lowest-frequency mode of a classical BH and whose damping times are parametrically longer. The reason for the appearance of the new modes is that our model contains another scale, the string length, which is parametrically larger than the Planck length. This distinction between the collapsed-polymer model and general-relativistic BHs could be made with gravitational-wave observations and offers a means for potentially measuring the strength of the coupling in string theory. For example, GW150914 already allows us to probe the strength of the string coupling near the regime which is predicted by the unification of the gravitational and gauge-theory couplings. We also derive bounds on the amplitude of the collapsed-polymer QNMs that can be placed by current and future gravitational-wave observations.


title: Stability and instability of the sub-extremal Reissner-Nordstr\"om black hole interior for the Einstein-Maxwell-Klein-Gordon equations in spherical symmetry

abstract: We show non-linear stability and instability results in spherical symmetry for the interior of a charged black hole -approaching a sub-extremal Reissner-Nordstr\"om background fast enough at infinity- in presence of a massive and charged scalar field, motivated by the strong cosmic censorship conjecture in that setting : 1. Stability : We prove that spherically symmetric characteristic initial data to the Einstein-Maxwell- Klein-Gordon equations approaching a Reissner-Nordstr\"om background with a sufficiently decaying polynomial decay rate on the event horizon gives rise to a space-time possessing a Cauchy horizon in a neighbourhood of time-like infinity. Moreover if the decay is even stronger, we prove that the spacetime metric admits a continuous extension to the Cauchy horizon. This generalizes the celebrated stability result of Dafermos for Einstein-Maxwell-real-scalar-field in spherical symmetry. 2. Instability : We prove that for the class of space-times considered in the stability part, whose scalar field in addition obeys a polynomial averaged-L^2 (consistent) lower bound on the event horizon, the scalar field obeys an integrated lower bound transversally to the Cauchy horizon. As a consequence we prove that the non-degenerate energy is infinite on any null surface crossing the Cauchy horizon and the curvature of a geodesic vector field blows up at the Cauchy horizon near time-like infinity. This generalizes an instability result due to Luk and Oh for Einstein-Maxwell-real-scalar-field in spherical symmetry. This instability of the black hole interior can also be viewed as a step towards the resolution of the C^2 strong cosmic censorship conjecture for one-ended asymptotically initial data.


title: On the universality of the nonstationary ideal

abstract: Burke \cite{MR1472122} proved that the generalized nonstationary ideal, denoted NS, is universal in the following sense: every normal ideal, and every tower of normal ideals of inaccessible height, is a canonical Rudin-Keisler projection of the restriction of $\text{NS}$ to some stationary set. We investigate how far Burke's theorem can be pushed, by analyzing the universality properties of NS with respect to the wider class of \emph{$\mathcal{C}$-systems of filters} introduced by Audrito-Steila \cite{AudritoSteila}. First we answer a question of \cite{AudritoSteila}, by proving that $\mathcal{C}$-systems of filters do not capture all kinds of set-generic embeddings. We provide a characterization of supercompactness in terms of short extenders and canonical projections of NS, without any reference to the strength of the extenders; as a corollary, NS can consistently fail to canonically project to arbitrarily strong short extenders. We prove that $\omega$-cofinal towers of normal ultrafilters---e.g.\ the kind used to characterize I2 and I3 embeddings---are well-founded if and only if they are canonical projections of NS. Finally, we provide a characterization of "$\aleph_\omega$ is Jonsson" in terms of canonical projections of NS.


title: d-Complete posets: local structural axioms, properties, and equivalent definitions

abstract: Although d-complete posets arose along the interface between algebraic combinatorics and Lie theory, they are defined using only requirements on their local structure. These posets are a mutual generalization of rooted trees, shapes, and shifted shapes. They possess Stanley's hook product property for their P-partition generating functions and Schutzenberger's well defined jeu de taquin rectification property. The original definition of d-complete poset was lengthy, but more succinct definitions were later developed. Here several definitions are shown to be equivalent. The basic properties of d-complete posets are summarized. Background and a partial bibliography for these posets is given.


title: Hilbert-Poincare series for spaces of commuting elements in Lie groups

abstract: In this article we study the homology of spaces ${\rm Hom}(\mathbb{Z}^n,G)$ of ordered pairwise commuting $n$-tuples in a Lie group $G$. We give an explicit formula for the Poincare series of these spaces in terms of invariants of the Weyl group of $G$. By work of Bergeron and Silberman, our results also apply to ${\rm Hom}(F_n/\Gamma_n^m,G)$, where the subgroups $\Gamma_n^m$ are the terms in the descending central series of the free group $F_n$. Finally, we show that there is a stable equivalence between the space ${\rm Comm}(G)$ studied by Cohen-Stafa and its nilpotent analogues.


title: Gravitational convergence, shear deformation and rotation of magnetic forcelines

abstract: We consider the "kinematics" of spacelike congruences and apply them to a family of self-gravitating magnetic forcelines. Our aim is to investigate the convergence and the possible focusing of these lines, as well as their rotation and shear deformation. In so going, we introduce a covariant 1+2 splitting of the 3-dimensional space, parallel and orthogonal to the direction of the field lines. The convergence, or the divergence, of the latter is monitored by the associated Raychaudhuri equation, which is obtained by propagating the spatial divergence of a unitary spacelike vector along its own direction. Applied to a magnetic vector, the resulting expression shows that, although the convergence of the magnetic forcelines is affected by the gravitational pull of all the other sources, it is unaffected by the field's own gravity, irrespective of how strong the latter is. This rather counterintuitive result is entirely due to the magnetic tension, namely to the negative pressure the field exerts parallel to its lines of force. In particular, the magnetic tension always cancels out the field's energy-density input to the Raychaudhuri equation, leaving the latter free of any direct magnetic-energy contribution. Similarly, the rotation and the shear deformation of the aforementioned forcelines are also unaffected by the magnetic input to the total gravitational energy. In a sense, the magnetic lines do not seem to "feel" their own gravity no matter how strong the latter may be.


title: Sorting sums of binary decision summands

abstract: A sum where each of the $N$ summands can be independently chosen from two choices yields $2^N$ possible summation outcomes. There is an $\mathcal{O}(K^2)$-algorithm that finds the $K$ smallest/largest of these sums by evading the enumeration of all sums.


title: Network Dissection: Quantifying Interpretability of Deep Visual Representations

abstract: We propose a general framework called Network Dissection for quantifying the interpretability of latent representations of CNNs by evaluating the alignment between individual hidden units and a set of semantic concepts. Given any CNN model, the proposed method draws on a broad data set of visual concepts to score the semantics of hidden units at each intermediate convolutional layer. The units with semantics are given labels across a range of objects, parts, scenes, textures, materials, and colors. We use the proposed method to test the hypothesis that interpretability of units is equivalent to random linear combinations of units, then we apply our method to compare the latent representations of various networks when trained to solve different supervised and self-supervised training tasks. We further analyze the effect of training iterations, compare networks trained with different initializations, examine the impact of network depth and width, and measure the effect of dropout and batch normalization on the interpretability of deep visual representations. We demonstrate that the proposed method can shed light on characteristics of CNN models and training methods that go beyond measurements of their discriminative power.


abstract: We consider Tikhonov regularization of control-constrained optimal control problems. We present new a-priori estimates for the regularization error assuming measure and source-measure conditions. In the special case of bang-bang solutions, we introduce another assumption to obtain the same convergence rates. This new condition turns out to be useful in the derivation of error estimates for the discretized problem. The necessity of the just mentioned assumptions to obtain certain convergence rates is analyzed. Finally, a numerical example confirms the analytical findings.


title: A complete dichotomy for complex-valued Holant^c

abstract: Holant problems are a family of counting problems on graphs, parametrised by sets of complex-valued functions of Boolean inputs. Holant^c denotes a subfamily of those problems, where any function set considered must contain the two unary functions pinning inputs to values 0 or 1. The complexity classification of Holant problems usually takes the form of dichotomy theorems, showing that for any set of functions in the family, the problem is either #P-hard or it can be solved in polynomial time. Previous such results include a dichotomy for real-valued Holant^c and one for Holant^c with complex symmetric functions. Here, we derive a dichotomy theorem for Holant^c with complex-valued, not necessarily symmetric functions. The tractable cases are the complex-valued generalisations of the tractable cases of the real-valued Holant^c dichotomy. The proof uses results from quantum information theory, particularly about entanglement.


title: Kinetic Monte Carlo simulations of the grain-surface back-diffusion effect

abstract: Rate-equation models are a widely-used and inexpensive tool for the simulation of interstellar chemistry under a range of physical conditions. However, their application to grain-surface chemical systems necessitates a number of simplifying assumptions, due to the requirement to treat only the total population of each species, using averaged rates, rather than treating each surface particle as an independent entity. While the outputs from rate-equation models are strictly limited to such population information, the inputs -- in the form of the averaged rates that control the time-evolution of chemical populations -- can be guided by the results from more exact simulation methods. Here, we examine the effects of back-diffusion, wherein particles diffusing on a surface revisit binding sites on the lattice, slowing the total reaction rate. While this effect has been studied for two-particle systems, its influence at greater surface coverage of reactants has not been explored. Results from two Monte Carlo kinetics models (one a 2-D periodic lattice, the other the surface of a three dimensionally-realized grain) were used to develop a means to incorporate the grain-surface back-diffusion effect into rate-equation methods. The effects of grain size, grain morphology, and surface coverage on the magnitude of the back-diffusion effect were studied for the simple H+H reaction system. The results were fit with expressions that can be easily incorporated into astrochemical rate-equation models to reproduce accurately the effects of back-diffusion on grain-surface reaction rates. Back-diffusion reduces reaction rates by a maximum factor of around 5 for the canonical grain of $\sim$10$^6$ surface sites, but this falls to unity at close to full surface coverage.


title: Design and Fabrication of Waveguide Optics for Imaging Applications

abstract: Due to the non-ionizing property, researchers have chosen to investigate terahertz radiation (THz) Imaging instrumentation for Bio-Sensing applications. The present work is to design and fabricate a near field lens that can focus guided terahertz radiation to a microscopic region for the detection of cancer-affected cells in Biological tissue. Operational characteristics such as field of view, optical loss factor, and hydrophobicity must be included to achieve an effective design of the lens.


title: Gender Disparities in Science? Dropout, Productivity, Collaborations and Success of Male and Female Computer Scientists

abstract: Scientific collaborations shape ideas as well as innovations and are both the substrate for, and the outcome of, academic careers. Recent studies show that gender inequality is still present in many scientific practices ranging from hiring to peer-review processes and grant applications. In this work, we investigate gender-specific differences in collaboration patterns of more than one million computer scientists over the course of 47 years. We explore how these patterns change over years and career ages and how they impact scientific success. Our results highlight that successful male and female scientists reveal the same collaboration patterns: compared to scientists in the same career age, they tend to collaborate with more colleagues than other scientists, seek innovations as brokers and establish longer-lasting and more repetitive collaborations. However, women are on average less likely to adapt the collaboration patterns that are related with success, more likely to embed into ego networks devoid of structural holes, and they exhibit stronger gender homophily as well as a consistently higher dropout rate than men in all career ages.


title: Small sets in dense pairs

abstract: Let $\widetilde{\mathcal M}=\langle \mathcal M, P\rangle$ be an expansion of an o-minimal structure $\mathcal M$ by a dense set $P\subseteq M$, such that three tameness conditions hold. We prove that the induced structure on $P$ by $\mathcal M$ eliminates imaginaries. As a corollary, we obtain that every small set $X$ definable in $\widetilde{\mathcal M}$ can be definably embedded into some $P^l$, uniformly in parameters, settling a question from [10]. We verify the tameness conditions in three examples: dense pairs of real closed fields, expansions of $\mathcal M$ by a dense independent set, and expansions by a dense divisible multiplicative group with the Mann property. Along the way, we point out a gap in the proof of a relevant elimination of imaginaries result in Wencel [17]. The above results are in contrast to recent literature, as it is known in general that $\widetilde{\mathcal M}$ does not eliminate imaginaries, and neither it nor the induced structure on $P$ admits definable Skolem functions.


title: Monitoring and manipulating Higgs and Goldstone modes in a supersolid quantum gas

abstract: Access to collective excitations lies at the heart of our understanding of quantum many-body systems. We study the Higgs and Goldstone modes in a supersolid quantum gas that is created by coupling a Bose-Einstein condensate symmetrically to two optical cavities. The cavity fields form a U(1)-symmetric order parameter that can be modulated and monitored along both quadratures in real time. This enables us to measure the excitation energies across the superfluid-supersolid phase transition, establish their amplitude and phase nature, as well as characterize their dynamics from an impulse response. Furthermore, we can give a tunable mass to the Goldstone mode at the crossover between continuous and discrete symmetry by changing the coupling of the quantum gas with either cavity.


title: Search for magnetic inelastic dark matter with XENON100

abstract: We present the first search for dark matter-induced delayed coincidence signals in a dual-phase xenon time projection chamber, using the 224.6 live days of the XENON100 science run II. This very distinct signature is predicted in the framework of magnetic inelastic dark matter which has been proposed to reconcile the modulation signal reported by the DAMA/LIBRA collaboration with the null results from other direct detection experiments. No candidate event has been found in the region of interest and upper limits on the WIMP's magnetic dipole moment are derived. The scenarios proposed to explain the DAMA/LIBRA modulation signal by magnetic inelastic dark matter interactions of WIMPs with masses of 58.0 GeV/c$^2$ and 122.7 GeV/c$^2$ are excluded at 3.3 $\sigma$ and 9.3 $\sigma$, respectively.


title: Darboux Integrability of Trapezoidal $H^{4}$ and $H^{6}$ Families of Lattice Equations II: General Solutions

abstract: In this paper we construct the general solutions of two families of quad-equations, namely the trapezoidal $H^{4}$ equations and the $H^{6}$ equations. These solutions are obtained exploiting the properties of the first integrals in the Darboux sense, which were derived in [Gubbiotti G., Yamilov R.I., J. Phys. A: Math. Theor. 50 (2017), 345205, 26 pages, arXiv:1608.03506]. These first integrals are used to reduce the problem to the solution of some linear or linearizable non-autonomous ordinary difference equations which can be formally solved.


title: Global Stability for a HIV/AIDS Model

abstract: We investigate global stability properties of a HIV/AIDS population model with constant recruitment rate, mass action incidence, and variable population size. Existence and uniqueness results for disease-free and endemic equilibrium points are proved. Global stability of the equilibria is obtained through Lyapunov's direct method and LaSalle's invariance principle.


title: Integrable Open Spin Chains from Flavored ABJM Theory

abstract: We compute the two-loop anomalous dimension matrix in the scalar sector of planar ${\cal N}=3$ flavored ABJM theory. Using coordinate Bethe ansatz, we obtain the reflection matrix and confirm that the boundary Yang-Baxter equations are satisfied. This establishes the integrability of this theory in the scalar sector at the two-loop order.


title: Reliable Probabilistic Gossip over Large-Scale Random Topologies

abstract: This paper studies reliability of probabilistic neighbor-aware gossip algorithms over three well- known large-scale random topologies, namely Bernoulli (or Erd\H{o}s-R\'enyi) graph, the random geometric graph, and the scale-free graph. We propose a new and simple algorithm which ensures higher reliability at lower message complexity than the three families of gossip algorithms over every topology in our study. We also present a uniform approach to model the reliability of probabilistic gossip algorithms in the different random graphs, whose properties, in fact, are quite different. In our model a forwarding probability is derived with consideration of parameters in gossip algorithms and graph properties. Our simulations show that our model gives a reasonable prediction of the trade-off between reliability and message complexity for all probabilistic neighbor-aware gossip algorithms in various random networks. Therefore, it allows to fine-tune the input parameters in the gossip protocols to achieve a desirable reliability with tolerable message complexity.


title: The free boundary Schur process and applications I

abstract: We investigate the free boundary Schur process, a variant of the Schur process introduced by Okounkov and Reshetikhin, where we allow the first and the last partitions to be arbitrary (instead of empty in the original setting). The pfaffian Schur process, previously studied by several authors, is recovered when just one of the boundary partitions is left free. We compute the correlation functions of the process in all generality via the free fermion formalism, which we extend with the thorough treatment of "free boundary states". For the case of one free boundary, our approach yields a new proof that the process is pfaffian. For the case of two free boundaries, we find that the process is not pfaffian, but a closely related process is. We also study three different applications of the Schur process with one free boundary: fluctuations of symmetrized last passage percolation models, limit shapes and processes for symmetric plane partitions, and for plane overpartitions.


title: The spectrum, radiation conditions and the Fredholm property for the Dirichlet Laplacian in a perforated plane with semi-infinite inclusions

abstract: We consider the spectral Dirichlet problem for the Laplace operator in the plane $\Omega^{\circ}$ with double-periodic perforation but also in the domain $\Omega^{\bullet}$ with a semi-infinite foreign inclusion so that the Floquet-Bloch technique and the Gelfand transform do not apply directly. We describe waves which are localized near the inclusion and propagate along it. We give a formulation of the problem with radiation conditions that provides a Fredholm operator of index zero. The main conclusion concerns the spectra $\sigma^{\circ}$ and $\sigma^{\bullet}$ of the problems in $\Omega^{\circ}$ and $\Omega^{\bullet},$ namely we present a concrete geometry which supports the relation $\sigma^{\circ}\varsubsetneqq\sigma^{\bullet}$ due to a new non-empty spectral band caused by the semi-infinite inclusion called an open waveguide in the double-periodic medium.


title: Online Weighted Degree-Bounded Steiner Networks via Novel Online Mixed Packing/Covering

abstract: We design the first online algorithm with poly-logarithmic competitive ratio for the edge-weighted degree-bounded Steiner forest(EW-DB-SF) problem and its generalized variant. We obtain our result by demonstrating a new generic approach for solving mixed packing/covering integer programs in the online paradigm. In EW-DB-SF we are given an edge-weighted graph with a degree bound for every vertex. Given a root vertex in advance we receive a sequence of terminal vertices in an online manner. Upon the arrival of a terminal we need to augment our solution subgraph to connect the new terminal to the root. The goal is to minimize the total weight of the solution while respecting the degree bounds on the vertices. In the offline setting edge-weighted degree-bounded Steiner tree (EW-DB-ST) and its many variations have been extensively studied since early eighties. Unfortunately the recent advancements in the online network design problems are inherently difficult to adapt for degree-bounded problems. In contrast in this paper we obtain our result by using structural properties of the optimal solution, and reducing the EW-DB-SF problem to an exponential-size mixed packing/covering integer program in which every variable appears only once in covering constraints. We then design a generic integral algorithm for solving this restricted family of IPs. We demonstrate a new technique for solving mixed packing/covering integer programs. Define the covering frequency k of a program as the maximum number of covering constraints in which a variable can participate. Let m denote the number of packing constraints. We design an online deterministic integral algorithm with competitive ratio of O(k log m) for the mixed packing/covering integer programs. We believe this technique can be used as an interesting alternative for the standard primal-dual techniques in solving online problems.


title: Gaussian and Weibull noncommutative charged black holes

abstract: We derive and investigate the physical properties of asymptotically flat noncommutative regular charged black holes with a Gaussian mass density distribution and a Weibull electric charge density distribution. Both distributions replace the Dirac one and introduce a set of substitution rules constituting new ways of counting. The solutions have a de Sitter behavior in the vicinity of the origin provided the electric charge density is Weibull of the form $r^{n/2}{\rm e}^{-r^2/(4\theta^2)}$ with $n\geq 1$. The electric field and temperature are finite for all values of the radial coordinate, mass, and charge. The charge is bounded from above for stability reasons and stable charged quantum particles have mass and charge bounded from below and from above within the simplified semi-classical model we present in this work.


title: Electron-phonon scattering effects on electronic and optical properties of orthorhombic GeS

abstract: Group-VI monochalcogenides are attracting a great deal of attention due to their peculiar anisotropic properties. Very recently, it has been suggested that GeS could act as a promissory absorbing material with high input-output ratios, relevant features for designing prospective optoelectronic devices. In this work, we use the \emph{ab-initio} many body perturbation theory to study the role of the electron-phonon coupling on orthorhombic GeS. We identify the vibrational modes that efficiently couple with the electronic states responsible for giving rise to the first and second excitonic state. We also study the finite-temperature optical absorption and show that even at $T\to0K$, the role of the electron-phonon interaction is crucial to properly describe the main experimental excitation peaks position and width. Our results suggest that the electron-phonon coupling is essential to properly describe the optical properties of the monochalcogenides family.


title: Multiplicative quiver varieties and generalised Ruijsenaars-Schneider models

abstract: We study some classical integrable systems naturally associated with multiplicative quiver varieties for the (extended) cyclic quiver with $m$ vertices. The phase space of our integrable systems is obtained by quasi-Hamiltonian reduction from the space of representations of the quiver. Three families of Poisson-commuting functions are constructed and written explicitly in suitable Darboux coordinates. The case $m=1$ corresponds to the tadpole quiver and the Ruijsenaars-Schneider system and its variants, while for $m>1$ we obtain new integrable systems that generalise the Ruijsenaars-Schneider system. These systems and their quantum versions also appeared recently in the context of supersymmetric gauge theory and cyclotomic DAHAs, as well as in the context of the Macdonald theory.


title: Headphones on the wire

abstract: We analyze a dataset providing the complete information on the effective plays of thousands of music listeners during several months. Our analysis confirms a number of properties previously highlighted by research based on interviews and questionnaires, but also uncover new statistical patterns, both at the individual and collective levels. In particular, we show that individuals follow common listening rhythms characterized by the same fluctuations, alternating heavy and light listening periods, and can be classified in four groups of similar sizes according to their temporal habits --- 'early birds', 'working hours listeners', 'evening listeners' and 'night owls'. We provide a detailed radioscopy of the listeners' interplay between repeated listening and discovery of new content. We show that different genres encourage different listening habits, from Classical or Jazz music with a more balanced listening among different songs, to Hip Hop and Dance with a more heterogeneous distribution of plays. Finally, we provide measures of how distant people are from each other in terms of common songs. In particular, we show that the number of songs $S$ a DJ should play to a random audience of size $N$ such that everyone hears at least one song he/she currently listens to, is of the form $S\sim N^\alpha$ where the exponent depends on the music genre and is in the range $[0.5,0.8]$. More generally, our results show that the recent access to virtually infinite catalogs of songs does not promote exploration for novelty, but that most users favor repetition of the same songs.


title: Analytical study of the "master-worker" framework scalability on multiprocessors with distributed memory

abstract: The paper is devoted to an analytical study of the "master-worker" framework scalability on multiprocessors with distributed memory. A new model of parallel computations called BSF is proposed. The BSF model is based on BSP and SPMD models. The scope of BSF model is the compute-intensive applications. The architecture of BSF-computer is defined. The structure of BSF-program is described. The Using this metric, the upper scalability bounds of BSF programs on distributed memory multiprocessors are evaluated. The formulas for estimating the parallel efficiency of BSF programs also proposed.


title: Learn to Model Motion from Blurry Footages

abstract: It is difficult to recover the motion field from a real-world footage given a mixture of camera shake and other photometric effects. In this paper we propose a hybrid framework by interleaving a Convolutional Neural Network (CNN) and a traditional optical flow energy. We first conduct a CNN architecture using a novel learnable directional filtering layer. Such layer encodes the angle and distance similarity matrix between blur and camera motion, which is able to enhance the blur features of the camera-shake footages. The proposed CNNs are then integrated into an iterative optical flow framework, which enable the capability of modelling and solving both the blind deconvolution and the optical flow estimation problems simultaneously. Our framework is trained end-to-end on a synthetic dataset and yields competitive precision and performance against the state-of-the-art approaches.


title: Anomalous Scaling of Stochastic Processes and the Moses Effect

abstract: The state of a stochastic process evolving over a time $t$ is typically assumed to lie on a normal distribution whose width scales like $t^{1/2}$. However, processes where the probability distribution is not normal and the scaling exponent differs from $\frac{1}{2}$ are known. The search for possible origins of such "anomalous" scaling and approaches to quantify them are the motivations for the work reported here. In processes with stationary increments, where the stochastic process is time-independent, auto-correlations between increments and infinite variance of increments can cause anomalous scaling. These sources have been referred to as the $\it{Joseph}$ $\it{effect}$ the $\it{Noah}$ $\it{effect}$, respectively. If the increments are non-stationary, then scaling of increments with $t$ can also lead to anomalous scaling, a mechanism we refer to as the $\it{Moses}$ $\it{effect}$. Scaling exponents quantifying the three effects are defined and related to the Hurst exponent that characterizes the overall scaling of the stochastic process. Methods of time series analysis that enable accurate independent measurement of each exponent are presented. Simple stochastic processes are used to illustrate each effect. Intraday Financial time series data is analyzed, revealing that its anomalous scaling is due only to the Moses effect. In the context of financial market data, we reiterate that the Joseph exponent, not the Hurst exponent, is the appropriate measure to test the efficient market hypothesis.

categories: q-fin.ST cond-mat.stat-mech math.DS physics.data-an


title: Coupling of COAMPS and WAVEWATCH with Improved Wave Physics

abstract: The Model Coupling Executable Library (MCEL), developed at the University of Southern Mississippi's Center of Higher Learning, has been successfully used to couple the Coupled Ocean/Atmospheric Mesoscale Prediction System (COAMPS) and the ocean wave model WAVEWATCH. An example of its application is shown for Hurricane Gordon, showing that two-way coupling results affects boundary layer physics differently than one-way coupling --- in this case, resulting in larger z_o and, consequently, larger surface fluxes and a more intense hurricane. However, since analyzing MCEL is difficult because the wave physics is inaccurate, improvements to the wave algorithms are also part of the deliverables. A new analytical expression for the wind/wave growth factor has been derived based on normal modes analysis and rapid distortion theory valid for all wave regimes except for tropical cyclone conditions. This new algorithm is validated against a numerical simulation of the Reynolds-stress transport equations and matches well. In contrast, other wave growth expression used in ocean models like the WAve Model (WAM) and WAVEWATCH do not produce the same results, with larger wave growth values peaking at smaller wave age values. These differences are attributed to the application of curve fitting by the other algorithms, while the new formulation is an analytical expression derived from first principles and includes factors missing in previous schemes such as turbulent interaction. If the Reynolds-stress transport equations solutions are reasonably accurate, it indicates that all the previous wave growth schemes, including WAVEWATCH, have serious deficiencies. Another unique result from this work includes a second analytical wave growth formulation valid for tropical cyclone conditions.


title: Noise-Tolerant Interactive Learning from Pairwise Comparisons

abstract: We study the problem of interactively learning a binary classifier using noisy labeling and pairwise comparison oracles, where the comparison oracle answers which one in the given two instances is more likely to be positive. Learning from such oracles has multiple applications where obtaining direct labels is harder but pairwise comparisons are easier, and the algorithm can leverage both types of oracles. In this paper, we attempt to characterize how the access to an easier comparison oracle helps in improving the label and total query complexity. We show that the comparison oracle reduces the learning problem to that of learning a threshold function. We then present an algorithm that interactively queries the label and comparison oracles and we characterize its query complexity under Tsybakov and adversarial noise conditions for the comparison and labeling oracles. Our lower bounds show that our label and total query complexity is almost optimal.


title: Geometric modeling of aerosol surface roughness and numerical simulation of its impact on aerosol optical properties

abstract: A stochastic method is introduced for geometric modeling aerosol surface roughness with random field and discrete differential geometry theory. Optical scattering properties are computed for randomly oriented spheroidal particles with uniformly random surface roughness. Invariant imbedding T-Matrix and geometric optics method are applied to compute light scattering properties of aerosol particles covered from Rayleigh to geometric optics region. We simulatied scattering phase matrix of model particles, better performace is obtained than smooth particles. We aslo introduce introduce entropy and relative entropy as measures to distinguish the optical scattering similarities between different particles.


title: Deterministic Quantum Annealing Expectation-Maximization Algorithm

abstract: Maximum likelihood estimation (MLE) is one of the most important methods in machine learning, and the expectation-maximization (EM) algorithm is often used to obtain maximum likelihood estimates. However, EM heavily depends on initial configurations and fails to find the global optimum. On the other hand, in the field of physics, quantum annealing (QA) was proposed as a novel optimization approach. Motivated by QA, we propose a quantum annealing extension of EM, which we call the deterministic quantum annealing expectation-maximization (DQAEM) algorithm. We also discuss its advantage in terms of the path integral formulation. Furthermore, by employing numerical simulations, we illustrate how it works in MLE and show that DQAEM outperforms EM.

categories: stat.ML cond-mat.stat-mech physics.comp-ph quant-ph


title: Interlacement of double curves of immersed spheres

abstract: We characterize those unions of embedded disjoint circles in the 2-sphere which can be the multiple point set of a generic immersion of the 2-sphere into 3-dimensional space in terms of the interlacement of the given circles. Our result is the one higher dimensional analogue of Rosenstiehl's characterization of words being Gauss codes of self-crossing plane curves. Our proof uses a result of Lippner and we further generalize the ideas of Fraysseix and Ossona de Mendez, which leads us to directed interlacement graphs of paired trees and their local complementation.


title: Anisotropic features in the electronic structure of the two-dimensional transition metal trichalcogenide TiS$_3$: electron doping and plasmons

abstract: Analysis of the band structure of TiS$_3$ single-layers suggests the possibility of changing their physical behaviour by injecting electron carriers. The anisotropy of the valence and conduction bands is explained in terms of their complex orbital composition. The nature of the Fermi surface and Lindhard response function for different doping concentrations is studied by means of first-principles DFT calculations. It is suggested that for electron doping levels $x$ (number of electrons per unit cell) $\sim$ 0.18-0.30$e^-$ the system could exhibit incommensurate charge or spin modulations which, however, would keep the metallic state whereas systems doped with smaller $x$ would be 2D metals without any electronic instability. The effect of spin-orbit coupling in the band dispersion is analysed. The DFT effective masses are used to study the plasmon spectrum from an effective low energy model. We find that this material supports highly anisotropic plasmons, with opposite anisotropy for the electron and hole bands.


title: Thermal Rayleigh-Marangoni convection in a three-layer liquid-metal-battery model

abstract: The combined effects of buoyancy-driven Rayleigh-B\'{e}nard convection (RC) and surface tension-driven Marangoni convection (MC) are studied in a triple-layer configuration which serves as a simplified model for a liquid metal battery (LMB). The three-layer model consists of a liquid metal alloy cathode, a molten salt separation layer, and a liquid metal anode at the top. Convection is triggered by the temperature gradient between the hot electrolyte and the colder electrodes, which is a consequence of the release of resistive heat during operation. We present a linear stability analysis of the state of pure thermal conduction in combination with three-dimensional direct numerical simulations of the nonlinear turbulent evolution on the basis of a pseudospectral method. Five different modes of convection are identified in the configuration, which are partly coupled to each other: RC in the upper electrode, RC with internal heating in the molten salt layer, MC at both interfaces between molten salt and electrode as well as anti-convectionin the middle layer and lower electrode. The linear stability analysis confirms that the additional Marangoni effect in the present setup increases the growth rates of the linearly unstable modes, i.e. Marangoni and Rayleigh-B\'{e}nard instability act together in the molten salt layer. The critical Grashof and Marangoni numbers decrease with increasing middle layer thickness. The calculated thresholds for the onset of convection are found for realistic current densities of laboratory-sized LMBs. The global turbulent heat transfer follows scaling predictions for internally heated RC. The global turbulent momentum transfer is comparable with turbulent convection in the classical Rayleigh-B\'{e}nard case.


title: Multi-scale detection of hierarchical community architecture in structural and functional brain networks

abstract: Community detection algorithms have been widely used to study the organization of complex systems like the brain. A principal appeal of these techniques is their ability to identify a partition of brain regions (or nodes) into communities, where nodes within a community are densely interconnected. In their simplest application, community detection algorithms are agnostic to the presence of community hierarchies, but a common characteristic of many neural systems is a nested hierarchy. To address this limitation, we exercise a multi-scale extension of a community detection technique known as modularity maximization, and we apply the tool to both synthetic graphs and graphs derived from human structural and functional imaging data. Our multi-scale community detection algorithm links a graph to copies of itself across neighboring topological scales, thereby becoming sensitive to conserved community organization across neighboring levels of the hierarchy. We demonstrate that this method allows for a better characterization of topological inhomogeneities of the graph's hierarchy by providing a local (node) measure of community stability and inter-scale reliability across topological scales. We compare the brain's structural and functional network architectures and demonstrate that structural graphs display a wider range of topological scales than functional graphs. Finally, we build a multimodal multiplex graph that combines structural and functional connectivity in a single model, and we identify the topological scales where resting state functional connectivity and underlying structural connectivity show similar versus unique hierarchical community architecture. Together, our results showcase the advantages of the multi-scale community detection algorithm in studying hierarchical community structure in brain graphs, and they illustrate its utility in modeling multimodal neuroimaging data.


title: Givental's non-linear Maslov index on lens spaces

abstract: Givental's non-linear Maslov index, constructed in 1990, is a quasimorphism on the universal cover of the identity component of the contactomorphism group of real projective space. This invariant was used by several authors to prove contact rigidity phenomena such as orderability, unboundedness of the discriminant and oscillation metrics, and a contact geometric version of the Arnold conjecture. In this article we give an analogue for lens spaces of Givental's construction and its applications.


title: GIER: A Danish computer from 1961 with a role in the modern revolution of astronomy

abstract: A Danish computer, GIER, from 1961 played a vital role in the development of a new method for astrometric measurement. This method, photon counting astrometry, ultimately led to two satellites with a significant role in the modern revolution of astronomy. A GIER was installed at the Hamburg Observatory in 1964 where it was used to implement the entirely new method for the measurement of stellar positions by means of a meridian circle, then the fundamental instrument of astrometry. An expedition to Perth in Western Australia with the instrument and the computer was a success. This method was also implemented in space in the first ever astrometric satellite Hipparcos launched by ESA in 1989. The Hipparcos results published in 1997 revolutionized astrometry with an impact in all branches of astronomy from the solar system and stellar structure to cosmic distances and the dynamics of the Milky Way. In turn, the results paved the way for a successor, the one million times more powerful Gaia astrometry satellite launched by ESA in 2013. Preparations for a Gaia successor in twenty years are making progress.


title: Kesten's bound for sub-exponential densities on the real line and its multi-dimensional analogues

abstract: We study the tail asymptotic of sub-exponential probability densities on the real line. Namely, we show that the n-fold convolution of a sub-exponential probability density on the real line is asymptotically equivalent to this density times n. We prove Kesten's bound, which gives a uniform in n estimate of the n-fold convolution by the tail of the density. We also introduce a class of regular sub-exponential functions and use it to find an analogue of Kesten's bound for functions on $\mathbb{R}^d$. The results are applied for the study of the fundamental solution to a nonlocal heat-equation.


title: Effects of geometric confinement on a droplet between two parallel planes

abstract: When a droplet (which size is characterized by R) is confined between two parallel planes, its morphology will change accordingly to either varying the volume of the droplet or the separation (characterized by h) between the planes. We are aiming at investigating how such a geometric confinement affects the wetting behaviours of a droplet. Our focus lies on two distinguished regimes: (1) a pancake shape in a Hele-Shaw cell when the droplet is highly compressed (i.e. h/R << 1), in which particular attention is paid on nonwetting and wetting cases, respectively; and (2) a liquid Hertzian contact rendered by a slight confinement (i.e. h/(2R) --> 1) in a nonwetting case. To realize this aim, we first develop strict analytical expressions of the shape of the droplet which are available for arbitrary contact angles between the liquid and the solid planes, but in which the elliptic integrals indicate that these solutions are implicit. By employing asymptotic methods, we are able to give expressions of relevant geometrical and physical parameters (the Laplace pressure, droplet volume, surface energy and external force) in terms of sole functions of R and h in an explicit manner. Comparisons suggest that over a large range of h/R, our asymptotic results quantitatively agree well with the numerical solutions of the analytical expressions. This systematic study of the parameter space allows a comprehensive understanding of the geometric confinement on wetting, especially a wide existence of logarithmic behaviours in a liquid Hertzian contact, to be identified.


title: Learning to Generate Long-term Future via Hierarchical Prediction

abstract: We propose a hierarchical approach for making long-term predictions of future frames. To avoid inherent compounding errors in recursive pixel-level prediction, we propose to first estimate high-level structure in the input frames, then predict how that structure evolves in the future, and finally by observing a single frame from the past and the predicted high-level structure, we construct the future frames without having to observe any of the pixel-level predictions. Long-term video prediction is difficult to perform by recurrently observing the predicted frames because the small errors in pixel space exponentially amplify as predictions are made deeper into the future. Our approach prevents pixel-level error propagation from happening by removing the need to observe the predicted frames. Our model is built with a combination of LSTM and analogy based encoder-decoder convolutional neural networks, which independently predict the video structure and generate the future frames, respectively. In experiments, our model is evaluated on the Human3.6M and Penn Action datasets on the task of long-term pixel-level video prediction of humans performing actions and demonstrate significantly better results than the state-of-the-art.


title: SkiMap: An Efficient Mapping Framework for Robot Navigation

abstract: We present a novel mapping framework for robot navigation which features a multi-level querying system capable to obtain rapidly representations as diverse as a 3D voxel grid, a 2.5D height map and a 2D occupancy grid. These are inherently embedded into a memory and time efficient core data structure organized as a Tree of SkipLists. Compared to the well-known Octree representation, our approach exhibits a better time efficiency, thanks to its simple and highly parallelizable computational structure, and a similar memory footprint when mapping large workspaces. Peculiarly within the realm of mapping for robot navigation, our framework supports realtime erosion and re-integration of measurements upon reception of optimized poses from the sensor tracker, so as to improve continuously the accuracy of the map.


title: Finite temperature QCD four-point function in the presence of a weak magnetic field within the hard thermal loop approximation

abstract: We compute the weak magnetic field correction to the four-point quark-gluon vertex at finite temperature and chemical potential within the hard thermal loop approximation. The resulting four-point quark-gluon vertex satisfies the Ward identity with the three-point quark-gluon vertex.


title: On large gaps between zeros of $L$-functions from branches

abstract: It is commonly believed that the normalized gaps between consecutive ordinates $t_n$ of the zeros of the Riemann zeta function on the critical line can be arbitrarily large. In particular, drawing on analogies with random matrix theory, it has been conjectured that $$\lambda' ={ lim ~ sup} ~( t_{n+1} - t_n ) \frac{ \log( t_n /2 \pi e)}{2\pi}$$ equals $\infty$. In this article we provide arguments, although not a rigorous proof, that $\lambda'$ is finite. Conditional on the Riemann Hypothesis, we show that if there are no changes of branch between consecutive zeros then $\lambda' \leq 3$, otherwise $\lambda'$ is allowed to be greater than $3$. Additional arguments lead us to propose $\lambda'\leq 5$. This proposal is consistent with numerous calculations that place lower bounds on $\lambda'$. We present the generalization of this result to all Dirichlet $L$-functions and those based on cusp forms.


title: A leptoquark model to accommodate $R_K^\mathrm{exp} <R_K^\mathrm{SM} $ and $R_{K^\ast}^\mathrm{exp}< R_{K^\ast}^\mathrm{SM}$

abstract: We show that a modification of the model with a low energy scalar leptoquark state carrying hypercharge $Y=7/6$ allows to accommodate both $R_K<1$ and $R_{K^\ast}<1$, through loop effects, consistent with recent observations made at LHCb. We describe details of the model, compute the relevant Wilson coefficient and, after discussing a number of constraints, we examine the phenomenological consequences of the model. The bounds on the lepton flavor violating decay rates, induced by this model, include $\mathcal{B}(Z\to\mu\tau) \lesssim \mathcal{O}(10^{-7})$, and $\mathcal{B}(B\to K\mu\tau) \lesssim \mathcal{O}(10^{-9})$. We also comment on the interpretation of the bounds on the leptoquark mass obtained from the direct searches at LHC.


title: Beating 1-1/e for Ordered Prophets

abstract: Hill and Kertz studied the prophet inequality on iid distributions [The Annals of Probability 1982]. They proved a theoretical bound of $1-\frac{1}{e}$ on the approximation factor of their algorithm. They conjectured that the best approximation factor for arbitrarily large n is $\frac{1}{1+1/e} \approx 0.731$. This conjecture remained open prior to this paper for over 30 years. In this paper we present a threshold-based algorithm for the prophet inequality with n iid distributions. Using a nontrivial and novel approach we show that our algorithm is a 0.738-approximation algorithm. By beating the bound of $\frac{1}{1+1/e}$, this refutes the conjecture of Hill and Kertz. Moreover, we generalize our results to non-iid distributions and discuss its applications in mechanism design.


title: Scalable computation of Jordan chains

abstract: We present an algorithm to compute the Jordan chain of a nearly defective matrix with a $2\times2$ Jordan block. The algorithm is based on an inverse-iteration procedure and only needs information about the invariant subspace corresponding to the Jordan chain, making it suitable for use with large matrices arising in applications, in contrast with existing algorithms which rely on an SVD. The algorithm produces the eigenvector and Jordan vector with $O(\varepsilon)$ error, with $\varepsilon$ being the distance of the given matrix to an exactly defective matrix. As an example, we demonstrate the use of this algorithm in a problem arising from electromagnetism, in which the matrix has size $212^2\times 212^2$. An extension of this algorithm is also presented which can achieve higher order convergence [$O(\varepsilon^2)$] when the matrix derivative is known.


title: Generative Face Completion

abstract: In this paper, we propose an effective face completion algorithm using a deep generative model. Different from well-studied background completion, the face completion task is more challenging as it often requires to generate semantically new pixels for the missing key components (e.g., eyes and mouths) that contain large appearance variations. Unlike existing nonparametric algorithms that search for patches to synthesize, our algorithm directly generates contents for missing regions based on a neural network. The model is trained with a combination of a reconstruction loss, two adversarial losses and a semantic parsing loss, which ensures pixel faithfulness and local-global contents consistency. With extensive experimental results, we demonstrate qualitatively and quantitatively that our model is able to deal with a large area of missing pixels in arbitrary shapes and generate realistic face completion results.


title: Entanglement Wedge Reconstruction via Universal Recovery Channels

abstract: We apply and extend the theory of universal recovery channels from quantum information theory to address the problem of entanglement wedge reconstruction in AdS/CFT. It has recently been proposed that any low-energy local bulk operators in a CFT boundary region's entanglement wedge can be reconstructed on that boundary region itself. Existing work arguing for this proposal relies on algebraic consequences of the exact equivalence between bulk and boundary relative entropies, namely the theory of operator algebra quantum error correction. However, bulk and boundary relative entropies are only approximately equal in bulk effective field theory, and in similar situations it is known that predictions from exact entropic equalities can be qualitatively incorrect. The framework of universal recovery channels provides a robust demonstration of the entanglement wedge reconstruction conjecture in addition to new physical insights. Most notably, we find that a bulk operator acting in a given boundary region's entanglement wedge can be expressed as the response of the boundary region's modular Hamiltonian to a perturbation of the bulk state in the direction of the bulk operator. This formula can be interpreted as a noncommutative version of Bayes' rule that attempts to undo the noise induced by restricting to only a portion of the boundary, and has an integral representation in terms of modular flows. To reach these conclusions, we extend the theory of universal recovery channels to finite-dimensional operator algebras and demonstrate that recovery channels approximately preserve the multiplicative structure of the operator algebra.


title: Squeezed Fourier Meets Toeplitz Algebras

abstract: We look for new steps on the dynamical operations that may squeeze simultaneously some families of quantum states, independently of their initial shape, induced by softly acting external fields which might produce the squeezing of the canonical observables $q,p$ of charged particles. Also, we present some exactly solvable cases of the problem which appear in the symmetric evolution intervals permitting to find explicitly the time dependence of the external fields needed to generate the required evolution operators. Curiously, our results are interrelated with a simple, non--trivial, anti--commuting algebra of Toeplitz which describes the problem more easily than the frequently used Ermakov--Milne invariants.


title: The Magic Barrier Revisited: Accessing Natural Limitations of Recommender Assessment

abstract: Recommender systems nowadays have many applications and are of great economic benefit. Hence, it is imperative for success-oriented companies to compare different of such systems and select the better one for their purposes. To this end, various metrics of predictive accuracy are commonly used, such as the Root Mean Square Error (RMSE), or precision and recall. All these metrics more or less measure how well a recommender system can predict human behaviour. Unfortunately, human behaviour is always associated with some degree of uncertainty, making the evaluation difficult, since it is not clear whether a deviation is system-induced or just originates from the natural variability of human decision making. At this point, some authors speculated that we may be reaching some Magic Barrier where this variability prevents us from getting much more accurate. In this article, we will extend the existing theory of the Magic Barrier into a new probabilistic but a yet pragmatic model. In particular, we will use methods from metrology and physics to develop easy-to-handle quantities for computation to describe the Magic Barrier for different accuracy metrics and provide suggestions for common application. This discussion is substantiated by comprehensive experiments with real users and large-scale simulations on a high-performance cluster.


title: Jet Substructure Studies with CMS Open Data

abstract: We use public data from the CMS experiment to study the 2-prong substructure of jets. The CMS Open Data is based on 31.8/pb of 7 TeV proton-proton collisions recorded at the Large Hadron Collider in 2010, yielding a sample of 768,687 events containing a high-quality central jet with transverse momentum larger than 85 GeV. Using CMS's particle flow reconstruction algorithm to obtain jet constituents, we extract the 2-prong substructure of the leading jet using soft drop declustering. We find good agreement between results obtained from the CMS Open Data and those obtained from parton shower generators, and we also compare to analytic jet substructure calculations performed to modified leading-logarithmic accuracy. Although the 2010 CMS Open Data does not include simulated data to help estimate systematic uncertainties, we use track-only observables to validate these substructure studies.


title: Properties of the cosmological filament between two clusters: A possible detection of a large-scale accretion shock by $Suzaku$

abstract: We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected value from a typical cluster temperature profile) to $kT\sim$6.5 keV. Our analysis also shows that filament plasma is present out to a radial distance of 15' (1.3 Mpc) from a line connecting the two clusters. The temperature profile is characterized by an almost flat radial shape with $kT\sim$6-7 keV within 10' or $\sim$0.8 Mpc. Across $r$=8'~from the axis, the temperature of the filament plasma shows a drop from 6.3 keV to 5.1 keV, indicating the presence of a shock front. The Mach number based on the temperature drop is estimated to be ${\cal M}\sim$1.3. We also successfully determined the abundance profile up to 15' (1.3 Mpc), showing an almost constant value ($Z$=0.3 solar) at the cluster outskirt. We estimated the Compton $y$-parameter to be $\sim$14.5$\pm1.3\times10^{-6}$, which is in agreement with $Planck$'s results (14-17$\times10^{-6}$ on the filament). The line of sight depth of the filament is $l\sim$1.1 Mpc, indicating that the geometry of filament is likely a pancake shape rather than cylindrical. The total mass of the filamentary structure is $\sim$7.7$\times10^{13}~\rm M_{\odot}$. We discuss a possible interpretation of the drop of X-ray emission at the rim of the filament, which was pushed out by the merging activity and formed by the accretion flow induced by the gravitational force of the filament.


title: Non-resonant Collider Signatures of a Singlet-Driven Electroweak Phase Transition

abstract: We analyze the collider signatures of the real singlet extension of the Standard Model in regions consistent with a strong first-order electroweak phase transition and a singlet-like scalar heavier than the Standard Model-like Higgs. A definitive correlation exists between the strength of the phase transition and the trilinear coupling of the Higgs to two singlet-like scalars, and hence between the phase transition and non-resonant scalar pair production involving the singlet at colliders. We study the prospects for observing these processes at the LHC and a future 100 TeV $pp$ collider, focusing particularly on double singlet production. We also discuss correlations between the strength of the electroweak phase transition and other observables at hadron and future lepton colliders. Searches for non-resonant singlet-like scalar pair production at 100 TeV would provide a sensitive probe of the electroweak phase transition in this model, complementing resonant di-Higgs searches and precision measurements. Our study illustrates a strategy for systematically exploring the phenomenologically viable parameter space of this model, which we hope will be useful for future work.


title: A small-scale dynamo in feedback-dominated galaxies - II. The saturation phase and the final magnetic configuration

abstract: Magnetic fields in galaxies are believed to be the result of dynamo amplification of initially weak seed fields, reaching equipartition strength inside the interstellar medium. The small-scale dynamo appears to be a viable mechanism to explain observations of strong magnetic fields in present-day and high-redshift galaxies, considering the extreme weakness of seed fields predicted by battery mechanisms or primordial fields. Performing high-resolution adaptive mesh magneto-hydrodynamic simulations of a small mass, isolated cooling halo with an initial magnetic seed field strength well below equipartition, we follow the small-scale dynamo amplification from supernova-induced turbulence up to saturation of the field. We find that saturation occurs when the average magnetic pressure reaches only 3 % to 5 % of the turbulent pressure. The magnetic energy growth transitions from exponential to linear, and finally comes to halt. The saturation level increases slightly with grid resolution. These results are in good agreement with theoretical predictions for magnetic Prandtl numbers of order $\mathrm{Pr_M} \sim 1$ and turbulent Mach numbers of order $\mathrm{M} \sim 10$. When we suppress supernova feedback after our simulation has reached saturation, we find that turbulence decays and that the gas falls back onto a thin disk with the magnetic field in local equipartition. We propose a scenario in which galactic magnetic fields are amplified from weak seed fields in the early stages of the Universe to sub-equipartition fields, owing to the turbulent environment of feedback-dominated galaxies at high redshift, and are evolved further in a later stage up to equipartition, as galaxies transformed into more quiescent, large spiral disks.


title: Integrated optics prototype beam combiner for long baseline interferometry in the L and M bands

abstract: In the last few years, integrated optics (IO) beam combiners have facilitated the emergence of 4-telescope interferometers such as PIONIER or GRAVITY, boosting the imaging capabilities of the VLTI. However, the spectral range beyond 2.2microns is not ideally covered by the conventional silica based IO. Here, we propose to consider new laser-written IO prototypes made of GLS glasses, a material that permits access to the mid-infrared spectral regime. Our goal is to conduct a full characterization of our mid-IR IO 2-telescope coupler in order to measure the performance levels directly relevant for long-baseline interferometry. We focus in particular on the exploitation of the L and M astronomical bands. We use a dedicated Michelson-interferometer setup to perform Fourier Transform spectroscopy on the coupler and measure its broadband interferometric performance. We also analyze the polarization properties of the coupler, the differential dispersion and phase degradation as well as the modal behavior and the total throughput. We measure broadband interferometric contrasts of 94.9% and 92.1% for unpolarized light in the L and M bands. Spectrally integrated splitting ratios are close to 50% but show chromatic dependence over the considered bandwidths. Additionally, the phase variation due to the combiner is measured and does not exceed 0.04rad and 0.07rad across the band L and M band, respectively. The total throughput of the coupler including Fresnel and injection losses from free-space is 25.4%. The laser-written IO GLS prototype combiners prove to be a reliable technological solution with promising performance for mid-infrared long-baseline interferometry. In the next steps, we will consider more advanced optical functions as well as a fiber-fed input and revise the optical design parameters in order the further enhance the total throughput and achromatic behavior.


title: The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

abstract: The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.


title: Kernel methods for interpretable machine learning of order parameters

abstract: Machine learning is capable of discriminating phases of matter, and finding associated phase transitions, directly from large data sets of raw state configurations. In the context of condensed matter physics, most progress in the field of supervised learning has come from employing neural networks as classifiers. Although very powerful, such algorithms suffer from a lack of interpretability, which is usually desired in scientific applications in order to associate learned features with physical phenomena. In this paper, we explore support vector machines (SVMs) which are a class of supervised kernel methods that provide interpretable decision functions. We find that SVMs can learn the mathematical form of physical discriminators, such as order parameters and Hamiltonian constraints, for a set of two-dimensional spin models: the ferromagnetic Ising model, a conserved-order-parameter Ising model, and the Ising gauge theory. The ability of SVMs to provide interpretable classification highlights their potential for automating feature detection in both synthetic and experimental data sets for condensed matter and other many-body systems.


title: Reconsidering the One Leptoquark solution: flavor anomalies and neutrino mass

abstract: We reconsider a model introducing a scalar leptoquark $\phi \sim (\mathbf{3}, \mathbf{1}, -1/3)$ to explain recent deviations from the standard model in semileptonic $B$ decays. The leptoquark can accommodate the persistent tension in the decays $\bar{B}\rightarrow D^{(*)}\tau \bar{\nu}$ as long as its mass is lower than approximately $10 \text{ TeV}$, and we show that a sizeable Yukawa coupling to the right-chiral tau lepton is necessary for an acceptable explanation. Agreement with the measured $\bar{B}\rightarrow D^{(*)}\tau \bar{\nu}$ rates is mildly compromised for parameter choices addressing the tensions in $b \to s \mu \mu$, where the model can significantly reduce the discrepancies in angular observables, branching ratios and the lepton-flavor-universality observables $R_K$ and $R_{K^*}$. The leptoquark can also reconcile the predicted and measured value of the anomalous magnetic moment of the muon and appears naturally in models of radiative neutrino mass derived from lepton-number violating effective operators. As a representative example, we incorporate the particle into an existing two-loop neutrino mass scenario derived from a dimension-nine operator. In this specific model, the structure of the neutrino mass matrix provides enough freedom to explain the small masses of the neutrinos in the region of parameter space dictated by agreement with the anomalies in $\bar{B}\rightarrow D^{(*)}\tau \bar{\nu}$, but not the $b \to s$ transition. This is achieved without excessive fine-tuning in the parameters important for neutrino mass.


title: Supernova Kicks and Dynamics of Compact Remnants in the Galactic Centre

abstract: The Galactic Centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we investigate the role of supernova (SN) explosions occurring in massive binary systems lying in a disc-like structure within the innermost parsec. We use a regularized algorithm to simulate 30,000 isolated three-body systems composed of a stellar binary orbiting the SMBH. We start the integration when the primary member undergoes a SN explosion, and we analyze the impact of SN kicks on the orbits of stars and compact remnants. We find that SN explosions scatter the lighter stars in the pair on completely different orbits, with higher eccentricity and inclination. In contrast, stellar-mass black holes (BHs) and massive stars retain memory of the orbit of their progenitor star. Our results suggest that SN kicks are not sufficient to eject BHs from the GC. We thus predict that all BHs that form in situ in the central parsec of our Galaxy remain in the GC, building up a cluster of dark remnants. In addition, the change of NS orbits induced by SNe may partially account for the observed dearth of NSs in the GC. About 40 per cent of remnants stay bound to the stellar companion after the kick; we expect up to 70 per cent of them might become X-ray binaries through Roche-lobe filling. Finally, the eccentricity of some light stars becomes > 0.7 as an effect of the SN kick, producing orbits similar to those of the G1 and G2 dusty objects.