title: Basker: A Threaded Sparse LU Factorization Utilizing Hierarchical Parallelism and Data Layouts

abstract: Scalable sparse LU factorization is critical for efficient numerical simulation of circuits and electrical power grids. In this work, we present a new scalable sparse direct solver called Basker. Basker introduces a new algorithm to parallelize the Gilbert-Peierls algorithm for sparse LU factorization. As architectures evolve, there exists a need for algorithms that are hierarchical in nature to match the hierarchy in thread teams, individual threads, and vector level parallelism. Basker is designed to map well to this hierarchy in architectures. There is also a need for data layouts to match multiple levels of hierarchy in memory. Basker uses a two-dimensional hierarchical structure of sparse matrices that maps to the hierarchy in the memory architectures and to the hierarchy in parallelism. We present performance evaluations of Basker on the Intel SandyBridge and Xeon Phi platforms using circuit and power grid matrices taken from the University of Florida sparse matrix collection and from Xyce circuit simulations. Basker achieves a geometric mean speedup of 5.91x on CPU (16 cores) and 7.4x on Xeon Phi (32 cores) relative to KLU. Basker outperforms Intel MKL Pardiso (PMKL) by as much as 53x on CPU (16 cores) and 13.3x on Xeon Phi (32 cores) for low fill-in circuit matrices. Furthermore, Basker provides 5.4x speedup on a challenging matrix sequence taken from an actual Xyce simulation.


title: Comet formation in collapsing pebble clouds. What cometary bulk density implies for the cloud mass and dust-to-ice ratio

abstract: Comets are remnants of the icy planetesimals that formed beyond the ice line in the Solar Nebula. Growing from micrometre-sized dust and ice particles to km-sized objects is, however, difficult because of growth barriers and time scale constraints. The gravitational collapse of pebble clouds that formed through the streaming instability may provide a suitable mechanism for comet formation. We study the collisional compression of cm-sized porous ice/dust-mixed pebbles in collapsing pebble clouds. For this, we developed a collision model for pebbles consisting of a mixture of ice and dust, characterised by their dust-to-ice mass ratio. Using the final compression of the pebbles, we constrain combinations of initial cloud mass, initial pepple porosity, and dust-to-ice ratio that lead to cometesimals which are consistent with observed bulk properties of cometary nuclei. We find that observed high porosity and low density of ~0.5 g/cc of comet nuclei can only be explained if comets formed in clouds with mass approximately M>1e18 g. Lower mass clouds would only work if the pebbles were initially very compact. Furthermore, the dust-to-ice ratio must be in the range of between 3 and 9 to match the observed bulk properties of comet nuclei. (abridged version)


title: A Gr\"atzer-Schmidt theorem for the Lindenbaum-Tarski algebra of IPC

abstract: We prove a version of the Gr\"atzer-Schmidt theorem for the Lindenbaum-Tarski algebra associated to the Implicational Propositional Calculus.


title: Rattling in spatially discrete diffusion equations with hysteresis

abstract: The paper treats a reaction-diffusion equation with hysteretic nonlinearity on a one-dimensional lattice. It arises as a result of the spatial discretization of the corresponding continuous model with so-called nontransverse initial data and exhibits a propagating microstructure --- which we call {\em rattling} --- in the hysteretic component of the solution. We analyze this microstructure and determine the speed of its propagation depending on the parameters of hysteresis and the nontransversality coefficient in the initial data.


title: Diphoton excess in phenomenological spin-2 resonance scenarios

abstract: We provide a possible explanation of a 750 GeV diphoton excess recently reported by both the ATLAS and CMS collaborations in the context of phenomenological spin-2 resonance scenarios, where the independent effective couplings of the resonance with gluons, quarks and photons are considered. We find a parameter region where the excess can be accounted for without conflicting with dijet constraints. We also show that the kinematical distributions might help to determine the couplings to gluons and quarks.


title: Non-collinear order and gapless superconductivity in s-wave magnetic superconductors

abstract: We study the behavior of magnetic superconductors which involve a local attractive interaction between electrons, and a coupling between local moments and the electrons. We solve this `Hubbard-Kondo' model through a variational minimization at zero temperature and validate the results via a Monte Carlo based on static auxiliary field decomposition of the Hubbard interaction. Over a magnetic coupling window that widens with increasing attractive interaction the ground state supports simultaneous magnetic and superconducting order. The pairing amplitude remains s-wave like, without significant spatial modulation, while the magnetic phase evolves from a ferromagnet, through non-collinear `spiral' states, to a Neel state with increasing density and magnetic coupling. We find that at intermediate magnetic coupling the antiferromagnetic-superconducting state is gapless, except for the regime of Neel order. We map out the phase diagram in terms of density, magnetic coupling and attractive interaction, establish the electron dispersion and effective `Fermi surface' in the ground state, provide an estimate of the magnetic and superconducting temperature scales via Monte Carlo, and compare our results to available data on the borocarbides.


title: Phase diagram and nucleation in the Polyakov-loop-extended Quark-Meson truncation of QCD with the unquenched Polyakov-loop potential

abstract: Unquenching of the Polyakov-loop potential showed to be an important improvement for the description of the phase structure and thermodynamics of strongly-interacting matter at zero quark chemical potentials with Polyakov-loop extended chiral models. This work constitutes the first application of the quark backreaction on the Polyakov-loop potential at nonzero density. The observation is that it links the chiral and deconfinement phase transition also at small temperatures and large quark chemical potentials. The build up of the surface tension in the Polyakov-loop extended Quark-Meson model is explored by investigating the two and 2+1-flavour Quark-Meson model and analysing the impact of the Polyakov-loop extension. In general, the order of magnitude of the surface tension is given by the chiral phase transition. The coupling of the chiral and deconfinement transition with the unquenched Polyakov-loop potential leads to the fact that the Polyakov-loop contributes at all temperatures.


title: Occupation probabilities and current densities of bulk and edge states of a Floquet topological insulator

abstract: Results are presented for the occupation probabilities and current densities of bulk and edge states of half-filled graphene in a cylindrical geometry, and irradiated by a circularly polarized laser. It is assumed that the system is closed, and that the laser has been switched on as a quench. Laser parameters corresponding to some representative topological phases are studied: one where the Chern number of the Floquet bands equals the number of chiral edge modes, a second where anomalous edge states appear in the Floquet Brillouin zone boundaries, and a third where the Chern number is zero, yet topological edge states appear at the center and boundaries of the Floquet Brillouin zone. Qualitative differences are found for the high frequency off-resonant and low frequency on-resonant laser with edge states arising due to resonant processes occupied with a high effective temperature on the one hand, while edge states arising due to off-resonant processes occupied with a low effective temperature on the other. For an ideal half-filled system where only one of the bands in the Floquet Brillouin zone is occupied and the other empty, particle-hole and inversion symmetry of the Floquet Hamiltonian implies zero current density. However the laser switch-on protocol breaks the inversion symmetry, resulting in a net cylindrical sheet of current density at steady-state. Due to the underlying chirality of the system, this current density profile is associated with a net charge imbalance between the top and bottom of the cylinders.


title: Hot ammonia around young O-type stars. III. High-mass star formation and hot core activity in W51~Main

abstract: This paper is the third in a series of ammonia multilevel imaging studies in well-known high-mass star forming regions. Using the JVLA, we have mapped the hot and dense molecular gas in W51 Main, with about 0.2 arcsec angular resolution, in five highly-excited metastable inversion transitions of ammonia (ammonia): (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13). We have identified and characterised two main centers of high-mass star formation in W51-Main: the W51e2 complex and the W51e8 core (6 arcsec southward of W51e2). The former breaks down into three further sub-cores: W51e2-W, which surrounds the well known HC HII region, where hot ammonia is observed in absorption, and two additional dusty cores, W51e2-E (~0.8" to the East) and W51e2-NW (~1" to the North), where hot ammonia is observed in emission. The velocity maps towards the HC HII region show a clear velocity gradient that may indicate rotation, though we do not directly observe a Keplerian velocity profile. The absence of outflow and/or maser activity and the low amount of molecular gas available for accretion (~5 solar masses) with respect to the mass of the central YSO (>20 solar masses), both indicate that the central YSO has already accreted most of its final mass. On the other hand, in the nearby W51e2-E object, the relatively large amount of hot molecular gas available for accretion (~20 solar masses, within about half an arcsecond or 2500 AU), along with strong outflow and maser activity, indicates that the main accretion center in the W51e2 complex is W51e2-E rather than W51e2-W. Finally, W51e2-NW and W51e8, although less dense, are also hot cores and contain a significant amount of molecular gas (~30 and 70 solar masses, respectively). We speculate that they may host high-mass YSOs either at a previous evolutionary stage or with lower mass than W51e2-E and W51e2-W.


title: Pressure-induced phase transformation in zircon-type orthovanadate SmVO4 from experiment and theory

abstract: The compression behavior of zircon-type samarium orthovanadate, SmVO4, has been investigated using synchrotron-based powder x-ray diffraction and ab-initio calculations up to 21 GPa. The results indicate the instability of ambient zircon phase at around 6 GPa, which transforms to a high-density scheelite-type phase. The high-pressure phase remains stable up to 21 GPa, the highest pressure reached in the present investigations. On pressure release, the scheelite phase is recovered. Crystal structure of high-pressure phase and equations of state (EOS) for the zircon- and scheelite-type phases have been determined. Various compressibilities such as bulk, axial and bond, estimated from the experimental data are found to be in good agreement with the results obtained from theoretical calculations. Calculated elastic constants show that the zircon structure becomes mechanically unstable beyond the transition pressure. Overall there is good agreement between experimental and theoretical findings.


title: Addressing spin transitions on 209Bi donors in silicon using circularly-polarized microwaves

abstract: Over the past decade donor spin qubits in isotopically enriched $^{28}$Si have been intensely studied due to their exceptionally long coherence times. More recently bismuth donor electron spins have become popular because Bi has a large nuclear spin which gives rise to clock transitions (first-order insensitive to magnetic field noise). At every clock transition there are two nearly degenerate transitions between four distinct states which can be used as a pair of qubits. Here it is experimentally demonstrated that these transitions are excited by microwaves of opposite helicity such that they can be selectively driven by varying microwave polarization. This work uses a combination of a superconducting coplanar waveguide (CPW) microresonator and a dielectric resonator to flexibly generate arbitrary elliptical polarizations while retaining the high sensitivity of the CPW.


title: High resolution spectroscopic analysis of seven giants in the bulge globular cluster NGC 6723

abstract: Globular clusters associated with the Galactic bulge are important tracers of stellar populations in the inner Galaxy. High resolution analysis of stars in these clusters allows us to characterize them in terms of kinematics, metallicity, and individual abundances, and to compare these fingerprints with those characterizing field populations. We present iron and element ratios for seven red giant stars in the globular cluster NGC~6723, based on high resolution spectroscopy. High resolution spectra ($R\sim48~000$) of seven K giants belonging to NGC 6723 were obtained with the FEROS spectrograph at the MPG/ESO 2.2m telescope. Photospheric parameters were derived from $\sim130$ FeI and FeII transitions. Abundance ratios were obtained from line-to-line spectrum synthesis calculations on clean selected features. An intermediate metallicity of [Fe/H]$=-0.98\pm0.08$ dex and a heliocentric radial velocity of $v_{hel}=-96.6\pm1.3~km s^{-1}$ were found for NGC 6723. Alpha-element abundances present enhancements of $[O/Fe]=0.29\pm0.18$ dex, $[Mg/Fe]=0.23\pm0.10$ dex, $[Si/Fe]=0.36\pm0.05$ dex, and $[Ca/Fe]=0.30\pm0.07$ dex. Similar overabundance is found for the iron-peak Ti with $[Ti/Fe]=0.24\pm0.09$ dex. Odd-Z elements Na and Al present abundances of $[Na/Fe]=0.00\pm0.21$ dex and $[Al/Fe]=0.31\pm0.21$ dex, respectively. Finally, the s-element Ba is also enhanced by $[Ba/Fe]=0.22\pm0.21$ dex. The enhancement levels of NGC 6723 are comparable to those of other metal-intermediate bulge globular clusters. In turn, these enhancement levels are compatible with the abundance profiles displayed by bulge field stars at that metallicity. This hints at a possible similar chemical evolution with globular clusters and the metal-poor of the bulge going through an early prompt chemical enrichment.


title: Baby Steps in Quantum Ring Theory: towards a background independent framework for Quantum Gravity

abstract: We investigate gravity as a gauge theory in the language of fiber bundles with tools from algebraic geometry. Compelled by the construction of the Eilenberg-MacLane classifying space via Fox derivations in an integral group ring, the origin of locally curved spacetime as a an abelian group extension is sought algebraically by means of a suitable derivation rule that conforms with the algebraic content of Einstein's equation. Accordingly, a peculiar flavor of background independence is advanced, curved spacetime should be thought of as an abelian proxy for the Lorentz group, causal and spin throughout. Thus, diffeomorphism invariance amounts to "choosing" an appropriate algebro-geometric--differential tool as to make algebraic sense of Einstein's equation or of a sensible evolution of such an equation. In this regard, it is suggested that a left-$\alpha$ derivation rule may facilitate the construction. On the other hand, we are led to interpret "spacetime curvature" in the broad sense of the exterior covariant derivative of a Cartan connection which includes both the torsion of the vierbein as well as the standard curvature of the spin connection. Making manifest use of the index structure in the standard Hilbert action--in particular, the distinction between internal and external indices--inclusion of the torsion 2-form is found to be inconsistent with contraction of the spacetime curvature as a scalar density. Aiming to erect a geometrically justified gravitational action that incorporates torsion, we lay out several fundamental geometric, differential and algebraic building blocks while highlighting descent relations among the different pieces. The groundwork is thus laid out for the search of the sought after action functional.


title: Diversifying Software Architecture for Sustainability: A Value-based Perspective

abstract: Although the concept of software diversity has been thoroughly adopted by software architects for many years, yet the advent of using diversity to achieve sustainability is overlooked. We argue that option thinking is an effective decision making tool to evaluate the trade-offs between architectural strategies and their long-term values under uncertainty. Our method extends cost-benefit analysis method CBAM. Unlike CBAM, our focus is on valuing the options which diversification can embed in the architecture and their corresponding value using real options pricing theory. The intuitive assumption is that the value of these options can provide the architect with insights on the long-term performance of these decisions in relation to some scenarios of interest and use them as the basis for reasoning about sustainability. The method aims to answer the following: (1) Is diversification of architectural decisions beneficial and when they can help in sustaining the software, (2) When, where and to what extent. The proposed model is illustrated and evaluated using a case study from the literature referred to as GridStix.


title: A heuristic derivation of radiative power loss and radiation reaction from the kinetic power of electric inertial mass of a charge

abstract: It is shown that formulas for the radiative power loss and radiation reaction from a charge can be derived in a heuristic manner from the kinetic power (rate of change of the kinetic energy) of its electric inertial mass. The derivation assumes a non-relativistic but otherwise an arbitrary motion of the charge. We exploit the fact that as the charge velocity changes because of a constant acceleration, there are accompanying modifications in its electromagnetic fields which can remain concurrent with the charge motion because the velocity as well as acceleration information enters into the field expression. However, if the acceleration of the charge is varying, information about that being not present in the field expressions, the electromagnetic fields get `out of step' with the actual charge motion. Accordingly we arrive at a radiation reaction formula for an arbitrarily moving charge, obtained hitherto in literature from the self-force, derived in a rather cumbersome way from the detailed mutual interaction between various constituents of a small charged sphere. This way we demonstrate that a power loss from a charge occurs only when there is a change in its acceleration and the derived instantaneous power loss is directly proportional to the scalar product of the velocity and the rate of change of the acceleration of the charge.


title: Local universality for real roots of random trigonometric polynomials

abstract: Consider a random trigonometric polynomial $X_n: \mathbb R \to \mathbb R$ of the form $$ X_n(t) = \sum_{k=1}^n \left( \xi_k \sin (kt) + \eta_k \cos (kt)\right), $$ where $(\xi_1,\eta_1),(\xi_2,\eta_2),\ldots$ are independent identically distributed bivariate real random vectors with zero mean and unit covariance matrix. Let $(s_n)_{n\in\mathbb N}$ be any sequence of real numbers. We prove that as $n\to\infty$, the number of real zeros of $X_n$ in the interval $[s_n+a/n, s_n+ b/n]$ converges in distribution to the number of zeros in the interval $[a,b]$ of a stationary, zero-mean Gaussian process with correlation function $(\sin t)/t$. We also establish similar local universality results for the centered random vectors $(\xi_k,\eta_k)$ having an arbitrary covariance matrix or belonging to the domain of attraction of a two-dimensional $\alpha$-stable law.


title: Efficient magnetic fields for supporting toroidal plasmas

abstract: The magnetic field that supports tokamak and stellarator plasmas must be produced by coils well separated from the plasma. However the larger the separation, the more difficult it is to produce a given magnetic field in the plasma region, so plasma configurations should be chosen that can be supported as efficiently as possible by distant coils. The efficiency of an externally-generated magnetic field is a measure of the field's shaping component magnitude at the plasma compared to the magnitude near the coils; the efficiency of a plasma equilibrium can be measured using the efficiency of the required external shaping field. Counterintuitively, plasma shapes with low curvature and spectral width may have low efficiency, whereas plasma shapes with sharp edges may have high efficiency. Two precise measures of magnetic field efficiency, which correctly identify such differences in difficulty, will be examined. These measures, which can be expressed as matrices, relate the externally-produced normal magnetic field on the plasma surface to the either the normal field or current on a distant control surface. A singular value decomposition (SVD) of either matrix yields an efficiency ordered basis for the magnetic field distributions. Calculations are carried out for both tokamak and stellarator cases. For axisymmetric surfaces with circular cross-section, the SVD is calculated analytically, and the range of poloidal and toroidal mode numbers that can be controlled to a given desired level is determined. If formulated properly, these efficiency measures are independent of the coordinates used to parameterize the surfaces.


title: Generalised unitarity for dimensionally regulated amplitudes within FDF

abstract: We review the Four-Dimensional-Formulation variant of the Four-Dimensional-Helicity scheme, by showing two applications of this regularisation scheme. The first one is the computation of one-loop helicity amplitudes, for which we present preliminary results for the analytic expressions of the one-loop Higgs plus five- gluon amplitudes. In the second part, we study the Colour-Kinematics duality for off-shell diagrams in gauge theories coupled to matter, showing in a diagrammatic way that the Jacobi relations for the kinematic numerators of off-shell diagrams, built with Feynman rules in axial gauge, reduce to definite set of violating terms due to the contributions of sub-graphs only.


title: Linear independence of indefinite iterated Eisenstein integrals

abstract: We prove linear independence of indefinite iterated Eisenstein integrals over the fraction field of the ring of formal power series $\mathbb{Z}[[q]]$. Our proof relies on a general criterium for linear independence of iterated integrals, which has been established by Deneufch{\^a}tel, Duchamp, Minh and Solomon. As a corollary, we obtain $\mathbb{C}$-linear independence of indefinite iterated Eisenstein integrals, which has applications to the study of elliptic multiple zeta values, as defined by Enriquez.


title: Sub-Optimal Multi-Phase Path Planning: A Method for Solving Rubik's Revenge

abstract: Rubik's Revenge, a 4x4x4 variant of the Rubik's puzzles, remains to date as an unsolved puzzle. That is to say, we do not have a method or successful categorization to optimally solve every one of its approximately $7.401 \times 10^{45}$ possible configurations. Rubik's Cube, Rubik's Revenge's predecessor (3x3x3), with its approximately $4.33 \times 10^{19}$ possible configurations, has only recently been completely solved by Rokicki et. al, further finding that any configuration requires no more than 20 moves. With the sheer dimension of Rubik's Revenge and its total configuration space, a brute-force method of finding all optimal solutions would be in vain. Similar to the methods used by Rokicki et. al on Rubik's Cube, in this paper we develop a method for solving arbitrary configurations of Rubik's Revenge in phases, using a combination of a powerful algorithm known as IDA* and a useful definition of distance in the cube space. While time-series results were not successfully gathered, it will be shown that this method far outweighs current human-solving methods and can be used to determine loose upper bounds for the cube space. Discussion will suggest that this method can also be applied to other puzzles with the proper transformations.


title: Some reflections on variational methods for partial differential inclusions

abstract: We discuss the application of variational methods, based on non-smooth critical point theory, to a general class of partial differential inclusions.


title: Effects of Sfermion Mixing induced by RGE in the CMSSM

abstract: Even within the Constrained Minimal Supersymmetric Standard Model (CMSSM) it is possible to induce sfermion flavor mixing through the Renormalization Group Equations (RGE) when the full structure of the Yukawa couplings is considered. We analyse the impact of including those effects on the accurate computation of $B$-physics observables, electroweak precision observables (EWPO) and the Higgs boson mass predictions.


title: Folding Flat Crease Patterns with Thick Materials

abstract: Modeling folding surfaces with nonzero thickness is of practical interest for mechanical engineering. There are many existing approaches that account for material thickness in folding applications. We propose a new systematic and broadly applicable algorithm to transform certain flat-foldable crease patterns into new crease patterns with similar folded structure but with a facet-separated folded state. We provide conditions on input crease patterns for the algorithm to produce a thickened crease pattern avoiding local self intersection, and provide bounds for the maximum thickness that the algorithm can produce for a given input. We demonstrate these results in parameterized numerical simulations and physical models.


title: Sampling-Based Query Re-Optimization

abstract: Despite of decades of work, query optimizers still make mistakes on "difficult" queries because of bad cardinality estimates, often due to the interaction of multiple predicates and correlations in the data. In this paper, we propose a low-cost post-processing step that can take a plan produced by the optimizer, detect when it is likely to have made such a mistake, and take steps to fix it. Specifically, our solution is a sampling-based iterative procedure that requires almost no changes to the original query optimizer or query evaluation mechanism of the system. We show that this indeed imposes low overhead and catches cases where three widely used optimizers (PostgreSQL and two commercial systems) make large errors.


title: Gravitational lensing of charged Ayon-Beato-Garcia black holes and non-linear effects of Maxwell fields

abstract: Non-singular Ayon-Beato-Garcia (ABG) spherically symmetric static black hole (BH) with charge to mass ratio $q$ is metric solution of Born-Infeld nonlinear Maxwell-Einstein theory. Central region of the BH behaves as (anti-) de Sitter for $(|q|>1)~|q|<1 .$ In case of $|q|=1$ the BH central region behaves as Minkowski flat metric. Nonlinear Electromagnetic (NEM) fields counterpart causes to deviate light geodesics and so light rays will forced to move on effective metric. In this paper we study weak and strong gravitational lensing of light rays by seeking affects of NEM fields counterpart on image locations and corresponding magnification. We set our calculations to experimentally observed Sgr A$^*$ BH. In short we obtained: For large distances the NEM counterpart is negligible reaching to linear Maxwell fields. The NEM makes enlarge the BH photon sphere radius as linearly by raising $|q|>1$ but deceases by raising $|q|\leq1.$ Sign of deflection angle of bending light rays is changed in presence of NEM effects with respect to ones obtained in absence of NEM fields. Absolute value of deflection angle raises by increasing $|q|\to1.$ Weak image locations decreases (increases) by raising $0<|q|<1$ in presence (absence) of NEM fields. By raising the closest distance of the bending light rays weak image locations changes from left (right) to right (left) in absence (presence) of NEM fields. Einstein rings radius and corresponding Magnification centroid become larger (smaller) in presence (absence) of NEM fields in case of weak lensing. Angular separation $s$ between the innermost and outermost relativistic images increases (decreases) by increasing $0<|q|<1$ in absence (presence) of NEM fields. Corresponding magnification $r$ decreases (increases) by raising $0<|q|<1$ in absence (presence) of NEM fields. $s (r)$ raises (decreases) by increasing $|q|>>1.$


title: The trilinear restriction estimate with sharp dependence on the transversality

abstract: We improve the Bennett--Carbery--Tao trilinear restriction estimate for subsets of the paraboloid in three dimensions, giving the sharp factor depending on the transversality.


title: On the Motion of a Free Particle in the de Sitter Manifold

abstract: Let $M=SO(1,4)/SO(1,3)\simeq S^{3}\times\mathbb{R}$ (a parallelizable manifold) be a submanifold in the structure $(\mathring{M}% ,\boldsymbol{\mathring{g}})$ (hereafter called the bulk) where $\mathring {M}\simeq\mathbb{R}^{5}$ and $\boldsymbol{\mathring{g}}$ is a pseudo Euclidian metric of signature $(1,4)$. Let $\boldsymbol{i}:M\rightarrow\mathbb{R}^{5}$ be the inclusion map and let \ $\boldsymbol{g}=\boldsymbol{i}^{\ast }\boldsymbol{\mathring{g}}$ be the pullback metric on $M$. It has signature $(1,3)$ Let $\boldsymbol{D}$ be the Levi-Civita connection of $\boldsymbol{g}% $. We call the structure $(M,\boldsymbol{g})$ a de Sitter manifold and $M^{dSL}=(M=\mathbb{R\times}S^{3},\boldsymbol{g},\boldsymbol{D},\tau _{\boldsymbol{g}},\uparrow)$ a de Sitter spacetime structure, which is \ of course orientable by $\tau_{\boldsymbol{g}}\in\sec% %TCIMACRO{\tbigwedge \nolimits^{4}}% %BeginExpansion {\textstyle\bigwedge\nolimits^{4}} %EndExpansion T^{\ast}M$ and time orientable (by $\uparrow$).\ Under these conditions we prove that if the motion of a free particle moving on $M$ happens with constant \emph{bulk} angular momentum then its motion in the structure $M^{dSL}$ is a timelike geodesic. Also any geodesic motion in the structure $M^{dSL}$ implies that the particle has constant angular momentum in the bulk.


title: New spiral state and skyrmion lattice in 3D model of chiral magnets

abstract: We present the phase diagram of magnetic states for films of isotropic chiral magnets calculated as function of applied magnetic field and thickness of the film. We have found a novel magnetic state driven by the natural confinement of the crystal, localized at the surface and stacked on top of the conical bulk phase. This magnetic surface state has a three-dimensional (3D) chiral spin-texture described by the superposition of helical and cycloidal spin spirals. This surface state exists for a large range of applied magnetic fields and for any film thickness beyond a critical one. We also identified the whole thickness and field range for which the skyrmion lattice becomes the ground state of the system. Below a certain critical thickness the surface state and bulk conical phase are suppressed in favor of the skyrmion lattice. Unraveling of those phases and the construction of the phase diagram became possible using advanced computational techniques for direct energy minimization applied to a basic 3D model for chiral magnets. Presented results provide a comprehensive theoretical description for those effects already observed in experiments on thin films of chiral magnets, predict new effects important for applications and open perspectives for experimental studies of such systems.


title: Ultra-confined acoustic THz graphene plasmons revealed by photocurrent nanoscopy

abstract: Terahertz (THz) fields are widely applied for sensing, communication and quality control. In future applications, they could be efficiently confined, enhanced and manipulated - well below the classical diffraction limit - through the excitation of graphene plasmons (GPs). These possibilities emerge from the strongly reduced GP wavelength, lp, compared to the photon wavelength, l0, which can be controlled by modulating the carrier density of graphene via electrical gating. Recently, GPs in a graphene-insulator-metal configuration have been predicted to exhibit a linear dispersion (thus called acoustic plasmons) and a further reduced wavelength, implying an improved field confinement, analogous to plasmons in two-dimensional electron gases (2DEGs) near conductive substrates. While infrared GPs have been visualised by scattering-type scanning near-field optical microscopy (s-SNOM), the real-space imaging of strongly confined THz plasmons in graphene and 2DEGs has been elusive so far - only GPs with nearly free-space wavelength have been observed. Here we demonstrate real-space imaging of acoustic THz plasmons in a graphene photodetector with split-gate architecture. To that end, we introduce nanoscale-resolved THz photocurrent near-field microscopy, where near-field excited GPs are detected thermoelectrically rather than optically. The on-chip GP detection simplifies GP imaging, as sophisticated s-SNOM detection schemes can be avoided. The photocurrent images reveal strongly reduced GP wavelengths (lp = l0/66), a linear dispersion resulting from the coupling of GPs with the metal gate below the graphene, and that plasmon damping at positive carrier densities is dominated by Coulomb impurity scattering. Acoustic GPs could thus strongly benefit the development of deep subwavelength-scale THz devices.


title: Applying a Differential Evolutionary Algorithm to a Constraint-based System to support Separation of OTDR Superimposed Signal after Passive Optical Network Splitters

abstract: The FTTH (Fiber To The Home) market currently needs new network maintenance technologies that can, economically and effectively, cope with massive fiber plants. However, operating these networks requires adequate means for an effective monitoring cost. Especially for troubleshooting faults that are associated with the possibility of remote identification of fiber breaks, which may exist in the network. Optical Time Domain Reflectometry (OTDR) techniques are widely used in point-to-point optical network topologies. Nevertheless, it has major limitations in tree-structured PONs (Passive Optical Networks), where all different branches backscatter the light in just one conventional OTDR trace with combined signals arriving on the OLT (Optical Line Terminal) side. Furthermore, passive power splitters used in FTTH networks input large attenuation, impoverishing the reflected signal. This makes the identification of the very branch affected by the problem practically impossible, when considering conventional analyses. The use of constraint-based techniques have been applied in a large amount of applications for Engineering Design, where the duties imposed for graphics and equations constraints result in valued features to CAD/CAE software capabilities. Currently, it provides a faster decision making capacity for engineers. This work applies the constraint based approach along with a Differential Evolutionary Algorithm to separate the superimposed OTDR signals, after the splitters of a FTTH Passive Optical Networks. This research introduces a new set of algorithms performing a coupling to an Optical Network (ON) CAD Design with its correspondent OTDR measurement signal, considering its geographical distribution branches of different lengths after the splitter. Results of this work are presented in a FTTN (Fiber To The Node) prototype arrangement, using a 1:8 passive power splitter.


title: The Anderson--Holstein Model in Two Flavors of the Non--Crossing Approximation

abstract: The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for non-equilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary non-crossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.


title: Strong convergence rates for nonlinearity-truncated Euler-type approximations of stochastic Ginzburg-Landau equations

abstract: This article proposes and analyzes explicit and easily implementable temporal numerical approximation schemes for additive noise-driven stochastic partial differential equations (SPDEs) with polynomial nonlinearities such as, e.g., stochastic Ginzburg-Landau equations. We prove essentially sharp strong convergence rates for the considered approximation schemes. Our analysis is carried out for abstract stochastic evolution equations on separable Banach and Hilbert spaces including the above mentioned SPDEs as special cases. We also illustrate our strong convergence rate results by means of a numerical simulation in Matlab.


title: Interference and dynamics of light from a distance-controlled atom pair in an optical cavity

abstract: Interference is central to quantum physics and occurs when indistinguishable paths exist, like in a double-slit experiment. Replacing the two slits with two single atoms introduces optical non-linearities for which nontrivial interference phenomena are predicted. Their observation, however, has been hampered by difficulties in preparing the required atomic distribution, controlling the optical phases and detecting the faint light. Here we overcome all of these experimental challenges by combining an optical lattice for atom localisation, an imaging system with single-site resolution, and an optical resonator for light steering. We observe resonator-induced saturation of resonance fluorescence for constructive interference of the scattered light and nonzero emission with huge photon bunching for destructive interference. The latter is explained by atomic saturation and photon pair generation. Our experimental setting is scalable and allows one to realize the Tavis-Cummings model for any number of atoms and photons, explore fundamental aspects of light-matter interaction, and implement new quantum information processing protocols.


title: A New Pivot Selection Algorithm for Symmetric Indefinite Factorization Arising in Quadratic Programming with Block Constraint Matrices

abstract: Quadratic programmingis a class of constrained optimization problem with quadratic objective functions and linear constraints. It has applications in many areas and is also used to solve nonlinear optimization problems. This article focuses on the equality constrained quadratic programs whose constraint matrices are block diagonal. Using the direct solution method, we propose a new pivot selection algorithm for the factorization of the Karush-Kuhn-Tucker(KKT) matrix for this problem that maintains the sparsity and stability of the problem. Our experiments show that our pivot selection algorithm appears to produce no fill-ins in the factorizationof such matrices. In addition, we compare our method with MA57 and find that the factors produced by our algorithm are sparser.


title: Fermionic scenario for the destruction of superconductivity in ultrathin MoC films evidenced by STM measurements

abstract: We use sub-Kelvin scanning tunneling spectroscopy to investigate the suppression of superconductivity in homogeneously disordered ultrathin MoC films. We observe that the superconducting state remains spatially homogeneous even on the films of 3 nm thickness. The vortex imaging suggests the global phase coherence in our films. Upon decreasing thickness, when the superconducting transition drops from 8.5 to 1.2 K, the superconducting energy gap follows perfectly Tc. All this is pointing to a two-stage fermionic scenario of the superconductor-insulator transition (SIT) via a metallic state as an alternative to the direct bosonic SIT scenario with a Cooper-pair insulating state evidenced by the last decade STM experiments.


title: On the accuracy of surface hopping dynamics in condensed phase non-adiabatic problems

abstract: We perform extensive benchmark comparisons of surface hopping dynamics with numerically exact calculations for the spin-boson model over a wide range of energetic and coupling parameters as well as temperature. We find that deviations from golden-rule scaling in the Marcus regime are generally small and depend sensitively on the energetic bias between electronic states. Fewest switched surface hopping (FSSH) is found to be surprisingly accurate over a large swath of parameter space. The inclusion of decoherence corrections via the augmented FSSH (A-FSSH) algorithm improves the accuracy of dynamical behavior compared to exact simulations, but the effects are generally not dramatic, at least for the case of an environment modeled with the commonly used Debye spectral density.


title: Super-resolution by means of Beurling minimal extrapolation

abstract: Let $M(\mathbb{T}^d)$ be the space of complex bounded Radon measures defined on the $d$-dimensional torus group $(\mathbb{R}/\mathbb{Z})^d=\mathbb{T}^d$, equipped with the total variation norm $\|\cdot\|$; and let $\hat\mu$ denote the Fourier transform of $\mu\in M(\mathbb{T}^d)$. We address the super-resolution problem: For given spectral (Fourier transform) data defined on a finite set $\Lambda\subset\mathbb{Z}^d$, determine if there is a unique $\mu\in M(\mathbb{T}^d)$ of minimal norm for which $\hat\mu$ equals this data on $\Lambda$. Without additional assumptions on $\mu$ and $\Lambda$, our main theorem shows that the solutions to the super-resolution problem, which we call minimal extrapolations, depend crucially on the set $\Gamma\subset\Lambda$, defined in terms of $\mu$ and $\Lambda$. For example, when $\Gamma=0$, the minimal extrapolations are singular measures supported in the zero set of an analytic function, and when $\Gamma\geq 2$, the minimal extrapolations are singular measures supported in the intersection of $\Gamma\choose 2$ hyperplanes. By theory and example, we show that the case $\Gamma=1$ is different from other cases and is deeply connected with the existence of positive minimal extrapolations. This theorem has implications to the possibility and impossibility of uniquely recovering $\mu$ from $\Lambda$. We illustrate how to apply our theory to both directions, by computing pertinent analytical examples. These examples are of interest in both super-resolution and deterministic compressed sensing. Our concept of an admissibility range fundamentally connects Beurling's theory of minimal extrapolation with Candes and Fernandez-Granda's work on super-resolution. This connection is exploited to address situations where current algorithms fail to compute a numerical solution to the super-resolution problem.


title: Cores, joins and the Fano-flow conjectures

abstract: The Fan-Raspaud Conjecture states that every bridgeless cubic graph has three 1-factors with empty intersection. A weaker one than this conjecture is that every bridgeless cubic graph has two 1-factors and one join with empty intersection. Both of these two conjectures can be related to conjectures on Fano-flows. In this paper, we show that these two conjectures are equivalent to some statements on cores and weak cores of a bridgeless cubic graph. In particular, we prove that the Fan-Raspaud Conjecture is equivalent to a conjecture proposed in [E. Steffen, 1-factor and cycle covers of cubic graphs, J. Graph Theory 78 (2015) 195-206]. Furthermore, we disprove a conjecture proposed in [G. Mazzuoccolo, New conjectures on perfect matchings in cubic graphs, Electron. Notes Discrete Math. 40 (2013) 235-238] and we propose a new version of it under a stronger connectivity assumption. The weak oddness of a cubic graph $G$ is the minimum number of odd components in the complement of a join of $G$. We obtain an upper bound of weak oddness in terms of weak cores, and thus an upper bound of oddness in terms of cores as a by-product.


title: Codeword Stabilized Quantum Codes for Asymmetric Channels

abstract: We discuss a method to adapt the codeword stabilized (CWS) quantum code framework to the problem of finding asymmetric quantum codes. We focus on the corresponding Pauli error models for amplitude damping noise and phase damping noise. In particular, we look at codes for Pauli error models that correct one or two amplitude damping errors. Applying local Clifford operations on graph states, we are able to exhaustively search for all possible codes up to length $9$. With a similar method, we also look at codes for the Pauli error model that detect a single amplitude error and detect multiple phase damping errors. Many new codes with good parameters are found, including nonadditive codes and degenerate codes.


title: A Confidence-Based Approach for Balancing Fairness and Accuracy

abstract: We study three classical machine learning algorithms in the context of algorithmic fairness: adaptive boosting, support vector machines, and logistic regression. Our goal is to maintain the high accuracy of these learning algorithms while reducing the degree to which they discriminate against individuals because of their membership in a protected group. Our first contribution is a method for achieving fairness by shifting the decision boundary for the protected group. The method is based on the theory of margins for boosting. Our method performs comparably to or outperforms previous algorithms in the fairness literature in terms of accuracy and low discrimination, while simultaneously allowing for a fast and transparent quantification of the trade-off between bias and error. Our second contribution addresses the shortcomings of the bias-error trade-off studied in most of the algorithmic fairness literature. We demonstrate that even hopelessly naive modifications of a biased algorithm, which cannot be reasonably said to be fair, can still achieve low bias and high accuracy. To help to distinguish between these naive algorithms and more sensible algorithms we propose a new measure of fairness, called resilience to random bias (RRB). We demonstrate that RRB distinguishes well between our naive and sensible fairness algorithms. RRB together with bias and accuracy provides a more complete picture of the fairness of an algorithm.


title: Model uncertainty and reference value of the Planck constant

abstract: Statistical parametric models are proposed to explain the values of the Planck constant obtained by comparing electrical and mechanical powers and by counting atoms in Si 28 enriched crystals. They assume that uncertainty contributions -- having heterogeneous, datum-specific, variances -- might not be included in the error budgets of some of the measured values. Model selection and model averaging are used to investigate data consistency, to identify a reference value of the Planck constant, and to include the model uncertainty in the error budget.


title: Adaptive confidence sets in shape restricted regression

abstract: A simple construction of adaptive confidence sets is proposed in isotonic, convex and unimodal regression. In univariate isotonic regression, the proposed confidence set enjoys uniform coverage over all non-decreasing regression functions. Furthermore, the diameter of the proposed confidence set automatically adapts to the unknown number of pieces of the true parameter, in the sense that the diameter is bounded from above by the minimax risk over the class of $k$-piecewise constant functions. The diameter of the confidence set is a simple increasing function of the number of jumps of the isotonic least-squares estimate. A similar construction is proposed in convex regression where the true regression function is convex and piecewise affine. Here, the confidence set enjoys uniform coverage and its diameter automatically adapt to the number of affine pieces of the true regression function. The diameter of the confidence set is an increasing function of the number of affine pieces of the convex least-squares estimate. We explain how to extend this technique to a non-convex set by proposing a similar adaptive confidence set in unimodal regression. The confidence set automatically adapts to the number of jumps of the true unimodal regression function and its diameter is an increasing function of the number of jumps of the unimodal least-squares estimate.


title: Spatial Scaling of Satellite Soil Moisture using Temporal Correlations and Ensemble Learning

abstract: A novel algorithm is developed to downscale soil moisture (SM), obtained at satellite scales of 10-40 km by utilizing its temporal correlations to historical auxiliary data at finer scales. Including such correlations drastically reduces the size of the training set needed, accounts for time-lagged relationships, and enables downscaling even in the presence of short gaps in the auxiliary data. The algorithm is based upon bagged regression trees (BRT) and uses correlations between high-resolution remote sensing products and SM observations. The algorithm trains multiple regression trees and automatically chooses the trees that generate the best downscaled estimates. The algorithm was evaluated using a multi-scale synthetic dataset in north central Florida for two years, including two growing seasons of corn and one growing season of cotton per year. The time-averaged error across the region was found to be 0.01 $\mathrm{m}^3/\mathrm{m}^3$, with a standard deviation of 0.012 $\mathrm{m}^3/\mathrm{m}^3$ when 0.02% of the data were used for training in addition to temporal correlations from the past seven days, and all available data from the past year. The maximum spatially averaged errors obtained using this algorithm in downscaled SM were 0.005 $\mathrm{m}^3/\mathrm{m}^3$, for pixels with cotton land-cover. When land surface temperature~(LST) on the day of downscaling was not included in the algorithm to simulate "data gaps", the spatially averaged error increased minimally by 0.015 $\mathrm{m}^3/\mathrm{m}^3$ when LST is unavailable on the day of downscaling. The results indicate that the BRT-based algorithm provides high accuracy for downscaling SM using complex non-linear spatio-temporal correlations, under heterogeneous micro meteorological conditions.


title: Syntax-Semantics Interaction Parsing Strategies. Inside SYNTAGMA

abstract: This paper discusses SYNTAGMA, a rule based NLP system addressing the tricky issues of syntactic ambiguity reduction and word sense disambiguation as well as providing innovative and original solutions for constituent generation and constraints management. To provide an insight into how it operates, the system's general architecture and components, as well as its lexical, syntactic and semantic resources are described. After that, the paper addresses the mechanism that performs selective parsing through an interaction between syntactic and semantic information, leading the parser to a coherent and accurate interpretation of the input text.


title: Can Negligible Cooperation Increase Network Reliability?

abstract: In network cooperation strategies, nodes work together with the aim of increasing transmission rates or reliability. This paper demonstrates that enabling cooperation between the transmitters of a two-user multiple access channel, via a cooperation facilitator that has access to both messages, always results in a network whose maximal- and average-error sum-capacities are the same---even when those capacities differ in the absence of cooperation and the information shared with the encoders is negligible. From this result, it follows that if a multiple access channel with no transmitter cooperation has different maximal- and average-error sum-capacities, then the maximal-error sum-capacity of the network consisting of this channel and a cooperation facilitator is not continuous with respect to the output edge capacities of the facilitator. This shows that there exist networks where sharing even a negligible number of bits per channel use with the encoders yields a non-negligible benefit.


title: Higher spin six vertex model and symmetric rational functions

abstract: We consider a fully inhomogeneous stochastic higher spin six vertex model in a quadrant. For this model we derive concise integral representations for multi-point q-moments of the height function and for the q-correlation functions. At least in the case of the step initial condition, our formulas degenerate in appropriate limits to many known formulas of such type for integrable probabilistic systems in the (1+1)d KPZ universality class, including the stochastic six vertex model, ASEP, various q-TASEPs, and associated zero range processes. Our arguments are largely based on properties of a family of symmetric rational functions which can be defined as partition functions of the inhomogeneous higher spin six vertex model for suitable domains. In the homogeneous case, such functions were previously studied in http://arxiv.org/abs/1410.0976; they also generalize classical Hall-Littlewood and Schur polynomials. A key role is played by Cauchy-like summation identities for these functions, which are obtained as a direct corollary of the Yang-Baxter equation for the higher spin six vertex model.


title: Nonperturbative True Muonium on the Light Front with TMSWIFT

abstract: The true muonium $(\mu\bar{\mu})$ bound state presents an interesting test of light-cone quantization techniques. In addition to exhibiting the standard problems of handling non-perturbative calculations, true muonium requires correct treatment of $e\bar{e}$ Fock-state contributions. Having previously produced a crude model of true muonium using the method of iterated resolvents, our current work has focused on the inclusion of the box diagrams to improve the cutoff-dependent issues of the model. Further, a parallel computer code, TMSWIFT, allowing for smaller numerical uncertainties, has been developed. This work focuses on the current state of these efforts to develop a model of true muonium that is testable at near-term experiments.


title: Phase Transitions and Duality in Adiabatic Memristive Networks

abstract: The development of neuromorphic systems based on memristive elements - resistors with memory - requires a fundamental understanding of their collective dynamics when organized in networks. Here, we study an experimentally inspired model of two-dimensional disordered memristive networks subject to a slowly ramped voltage and show that they undergo a first-order phase transition in the conductivity for sufficiently high values of memory, as quantified by the memristive ON/OFF ratio. We investigate the consequences of this transition for the memristive current-voltage characteristics both through simulation and theory, and uncover a duality between forward and reverse switching processes that has also been observed in several experimental systems of this sort. Our work sheds considerable light on the statistical properties of memristive networks that are presently studied both for unconventional computing and as models of neural networks.


title: Galaxy Cluster Mass Estimation from Stacked Spectroscopic Analysis

abstract: We use simulated galaxy surveys to study: i) how galaxy membership in redMaPPer clusters maps to the underlying halo population, and ii) the accuracy of a mean dynamical cluster mass, $M_\sigma(\lambda)$, derived from stacked pairwise spectroscopy of clusters with richness $\lambda$. Using $\sim\! 130,000$ galaxy pairs patterned after the SDSS redMaPPer cluster sample study of Rozo et al. (2015 RMIV), we show that the pairwise velocity PDF of central--satellite pairs with $m_i < 19$ in the simulation matches the form seen in RMIV. Through joint membership matching, we deconstruct the main Gaussian velocity component into its halo contributions, finding that the top-ranked halo contributes $\sim 60\%$ of the stacked signal. The halo mass scale inferred by applying the virial scaling of Evrard et al. (2008) to the velocity normalization matches, to within a few percent, the log-mean halo mass derived through galaxy membership matching. We apply this approach, along with mis-centering and galaxy velocity bias corrections, to estimate the log-mean matched halo mass at $z=0.2$ of SDSS redMaPPer clusters. Employing the velocity bias constraints of Guo et al. (2015), we find $\langle \ln(M_{200c})|\lambda \rangle = \ln(M_{30}) + \alpha_m \ln(\lambda/30)$ with $M_{30} = 1.56 \pm 0.35 \times 10^{14} M_\odot$ and $\alpha_m = 1.31 \pm 0.06_{stat} \pm 0.13_{sys}$. Systematic uncertainty in the velocity bias of satellite galaxies overwhelmingly dominates the error budget.


title: Factorization of stochastic maps using the Stinespring representations

abstract: In this work, we investigate the existence of a factorization for a unital completely positive map, between non-commutative probability space which do not change the expectation values of the events. These maps are called in literature stochastic maps. Using the Stinespring representations of completely positive map and assuming the existence of anti-unitary operator on Hilbert space related to these representations which satisfying some modular relations, we prove that stochastic maps with adjoint, admits a factorization.


title: Local Network Community Detection with Continuous Optimization of Conductance and Weighted Kernel K-Means

abstract: Local network community detection is the task of finding a single community of nodes concentrated around few given seed nodes in a localized way. Conductance is a popular objective function used in many algorithms for local community detection. This paper studies a continuous relaxation of conductance. We show that continuous optimization of this objective still leads to discrete communities. We investigate the relation of conductance with weighted kernel k-means for a single community, which leads to the introduction of a new objective function, $\sigma$-conductance. Conductance is obtained by setting $\sigma$ to $0$. Two algorithms, EMc and PGDc, are proposed to locally optimize $\sigma$-conductance and automatically tune the parameter $\sigma$. They are based on expectation maximization and projected gradient descent, respectively. We prove locality and give performance guarantees for EMc and PGDc for a class of dense and well separated communities centered around the seeds. Experiments are conducted on networks with ground-truth communities, comparing to state-of-the-art graph diffusion algorithms for conductance optimization. On large graphs, results indicate that EMc and PGDc stay localized and produce communities most similar to the ground, while graph diffusion algorithms generate large communities of lower quality.


title: Wireless Network Simplification : Beyond Diamond Networks

abstract: We consider an arbitrary layered Gaussian relay network with $L$ layers of $N$ relays each, from which we select subnetworks with $K$ relays per layer. We prove that: (i) For arbitrary $L, N$ and $K = 1$, there always exists a subnetwork that approximately achieves $\frac{2}{(L-1)N + 4}$ $\left(\mbox{resp.}\frac{2}{LN+2}\right)$ of the network capacity for odd $L$ (resp. even $L$), (ii) For $L = 2, N = 3, K = 2$, there always exists a subnetwork that approximately achieves $\frac{1}{2}$ of the network capacity. We also provide example networks where even the best subnetworks achieve exactly these fractions (up to additive gaps). Along the way, we derive some results on MIMO antenna selection and capacity decomposition that may also be of independent interest.


title: Testing ATLAS Z+MET Excess with LHC Run 2

abstract: The ATLAS collaboration reported a 3sigma excess in the search of events containing on-Z dilepton, jets, and large missing momentum (MET) in the 8 TeV LHC run. Motivated by this excess, many models of new physics have been proposed. Recently, the ATLAS and CMS collaborations reported new results for similar Z+MET channels in the 13 TeV run. In this paper, we comprehensively discuss the consistency between the proposed models and the LHC results of Run 1 and Run 2. We find that in models with heavy gluino production, there is generically some tension between the 8 TeV and 13 TeV results. On the other hand, models with light squark production provide relatively better fitting to both results.


title: The Maillet--Malgrange type theorem for generalized power series

abstract: There is proposed the Maillet--Malgrange type theorem for a generalized power series (having complex power exponents) formally satisfying an algebraic ordinary differential equation. The theorem describes the growth of the series coefficients.


title: CosmoLike - Cosmological Likelihood Analyses for Photometric Galaxy Surveys

abstract: We explore strategies to extract cosmological constraints from a joint analysis of cosmic shear, galaxy-galaxy lensing, galaxy clustering, cluster number counts and cluster weak lensing. We utilize the CosmoLike software to simulate results from an LSST like data set, specifically, we 1) compare individual and joint analyses of the different probes, 2) vary the selection criteria for lens and source galaxies, 3) investigate the impact of blending, 4) investigate the impact of the assumed cosmological model in multi-probe covariances, 6) quantify information content as a function of scales, and 7) explore the impact of intrinsic galaxy alignment in a multi-probe context. Our analyses account for all cross correlations within and across probes and include the higher-order (non-Gaussian) terms in the multi-probe covariance matrix. We simultaneously model cosmological parameters and a variety of systematics, e.g. uncertainties arising from shear and photo-z calibration, cluster mass-observable relation, galaxy intrinsic alignment, and galaxy bias (up to 54 parameters altogether). We highlight two results: First, increasing the number density of source galaxies by ~30%, which corresponds to solving blending for LSST, only gains little information. Second, including small scales in clustering and galaxy-galaxy lensing, by utilizing HODs, can substantially boost cosmological constraining power. The CosmoLike modules used to compute the results in this paper will be made publicly available at https://github.com/elikrause/CosmoLike_Forecasts.


title: Majorana-time-reversal symmetries: a fundamental principle for sign-problem-free quantum Monte Carlo simulations

abstract: A fundamental open issue in physics is whether and how the fermion sign problem in quantum Monte Carlo (QMC) simulations can be solved generically. Here, we show that Majorana-time-reversal (MTR) symmetries can provide a unifying principle to solve the fermion sign problem in interacting fermionic models. By systematically classifying Majorana-bilinear operators according to the anti-commuting MTR symmetries they respect, we rigorously proved that there are two and only two fundamental symmetry classes which are sign-problem-free and which we call the "Majorana class" and "Kramers class", respectively. Novel sign-problem-free models in the Majorana class include interacting topological superconductors and interacting models of charge-4e superconductors. We believe that our MTR unifying principle could shed new light on sign-problem-free QMC simulation on strongly correlated systems and interacting topological matters.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.supr-con hep-lat


title: Transformation properties for Dyson's rank function

abstract: At the 1987 Ramanujan Centenary meeting Dyson asked for a coherent group-theoretical structure for Ramanujan's mock theta functions analogous to Hecke's theory of modular forms. Many of Ramanujan's mock theta functions can be written in terms of $R(\zeta,q)$, where $R(z,q)$ is the two-variable generating function of Dyson's rank function and $\zeta$ is a root of unity. Building on earlier work of Watson, Zwegers, Gordon and McIntosh, and motivated by Dyson's question, Bringmann, Ono and Rhoades studied transformation properties of $R(\zeta,q)$. In this paper we strengthen and extend the results of Bringmann, Rhoades and Ono, and the later work of Ahlgren and Treneer. As an application we give a new proof of Dyson's rank conjecture and show that Ramanujan's Dyson rank identity modulo $5$ from the Lost Notebook has an analogue for all primes greater than $3$. The proof of this analogue was inspired by recent work of Jennings-Shaffer on overpartition rank differences mod $7$.


title: First-Principles Equation of State and Electronic Properties of Warm Dense Oxygen

abstract: We perform all-electron path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of oxygen. Our simulations cover a wide density-temperature range of $1-100$~g$\,$cm$^{-3}$ and $10^4-10^9$~K. By combining results from PIMC and DFT-MD, we are able to compute pressures and internal energies from first-principles at all temperatures and provide a coherent equation of state. We compare our first-principles calculations with analytic equations of state, which tend to agree for temperatures above 8$\times$10$^6$~K. Pair-correlation functions and the electronic density of states reveal an evolving plasma structure and ionization process that is driven by temperature and density. As we increase the density at constant temperature, we find that the ionization fraction of the 1s state decreases while the other electronic states move towards the continuum. Finally, the computed shock Hugoniot curves show an increase in compression as the first and second shells are ionized.


title: Uniform observability estimates for linear waves

abstract: In this article, we give a completely constructive proof of the observability/controllability of the wave equation on a compact manifold under optimal geometric conditions. This contrasts with the original proof of Bardos-Lebeau-Rauch, which contains two non-constructive arguments. Our method is based on the Dehman-Lebeau Egorov approach to treat the high-frequencies, and the optimal unique continuation stability result of the authors for the low-frequencies. As an application, we first give estimates of the blowup of the observability constant when the time tends to the limit geometric control time (for wave equations with possibly lower order terms). Second, we provide (on manifolds with or without boundary) with an explicit dependence of the observability constant with respect to the addition of a bounded potential to the equation.


title: A Note on Antenna Selection in Gaussian MIMO Channels: Capacity Guarantees and Bounds

abstract: We consider the problem of selecting $k_t \times k_r$ antennas from a Gaussian MIMO channel with $n_t \times n_r$ antennas, where $k_t \leq n_t$ and $k_r \leq n_r$. We prove the following two results that hold universally, in the sense that they do not depend on the channel coefficients: (i) The capacity of the best $k_t \times k_r$ subchannel is always lower bounded by a fraction $\frac{k_t k_r}{n_t n_r}$ of the full capacity (with $n_t \times n_r$ antennas). This bound is tight as the channel coefficients diminish in magnitude. (ii) There always exists a selection of $k_t \times k_r$ antennas (including the best) that achieves a fraction greater than $\frac{\min(k_t ,k_r)}{\min(n_t,n_r)}$ of the full capacity within an additive constant that is independent of the coefficients in the channel matrix. The key mathematical idea that allows us to derive these universal bounds is to directly relate the determinants of principle sub-matrices of a Hermitian matrix to the determinant of the entire matrix.


title: Atom-optics simulator of lattice transport phenomena

abstract: We experimentally investigate a scheme for studying lattice transport phenomena, based on the controlled momentum-space dynamics of ultracold atomic matter waves. In the effective tight-binding models that can be simulated, we demonstrate that this technique allows for a local and time-dependent control over all system parameters, and additionally allows for single-site resolved detection of atomic populations. We demonstrate full control over site-to-site off-diagonal tunneling elements (amplitude and phase) and diagonal site-energies, through the observation of continuous-time quantum walks, Bloch oscillations, and negative tunneling. These capabilities open up new prospects in the experimental study of disordered and topological systems.


title: CFHTLenS revisited: assessing concordance with Planck including astrophysical systematics

abstract: We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited


title: First-Principles Equation of State Calculations of Warm Dense Nitrogen

abstract: Using path integral Monte Carlo (PIMC) and density functional molecular dynamics (DFT-MD) simulation methods, we compute a coherent equation of state (EOS) of nitrogen that spans the liquid, warm dense matter (WDM), and plasma regimes. Simulations cover a wide range of density-temperature space, $1.5-13.9$~g$~$cm$^{-3}$ and $10^3-10^9$~K. In the molecular dissociation regime, we extend the pressure-temperature phase diagram beyond previous studies, providing dissociation and Hugoniot curves in good agreement with experiments and previous DFT-MD work. Analysis of pair-correlation functions and the electronic density of states in the WDM regime reveals an evolving plasma structure and ionization process that is driven by temperature and pressure. Our Hugoniot curves display a sharp change in slope in the dissociation regime and feature two compression maxima as the K and L shells are ionized in the WDM regime, which have some significant differences from the predictions of plasma models.


title: Automatic Matching of Bullet Land Impressions

abstract: In 2009, the National Academy of Sciences published a report questioning the scientific validity of many forensic methods including firearm examination. Firearm examination is a forensic tool used to help the court determine whether two bullets were fired from the same gun barrel. During the firing process, rifling, manufacturing defects, and impurities in the barrel create striation marks on the bullet. Identifying these striation markings in an attempt to match two bullets is one of the primary goals of firearm examination. We propose an automated framework for the analysis of the 3D surface measurements of bullet land impressions which transcribes the individual characteristics into a set of features that quantify their similarities. This makes identification of matches easier and allows for a quantification of both matches and matchability of barrels. The automatic matching routine we propose manages to (a) correctly identify land impressions (the surface between two bullet groove impressions) with too much damage to be suitable for comparison, and (b) correctly identify all 10,384 land-to-land matches of the James Hamby study.


title: Weak Antilocalization and Universal Conductance Fluctuations in Bismuth Telluro-Sulfide Topological Insulators

abstract: We report on van der Waals epitaxial growth, materials characterization and magnetotransport experiments in crystalline nanosheets of Bismuth Telluro-Sulfide (BTS). Highly layered, good-quality crystalline nanosheets of BTS are obtained on SiO$_2$ and muscovite mica. Weak-antilocalization (WAL), electron-electron interaction-driven insulating ground state and universal conductance fluctuations are observed in magnetotransport experiments on BTS devices. Temperature, thickness and magnetic field dependence of the transport data indicate the presence of two-dimensional surface states along with bulk conduction, in agreement with theoretical models. An extended-WAL model is proposed and utilized in conjunction with a two-channel conduction model to analyze the data, revealing a surface component and evidence of multiple conducting channels. A facile growth method and detailed magnetotransport results indicating BTS as an alternative topological insulator material system are presented.


title: How quickly can anyons be braided? Or: How I learned to stop worrying about diabatic errors and love the anyon

abstract: Topological phases of matter are a potential platform for the storage and processing of quantum information with intrinsic error rates that decrease exponentially with inverse temperature and with the length scales of the system, such as the distance between quasiparticles. However, it is less well-understood how error rates depend on the speed with which non-Abelian quasiparticles are braided. In general, diabatic corrections to the holonomy or Berry's matrix vanish at least inversely with the length of time for the braid, with faster decay occurring as the time-dependence is made smoother. We show that such corrections will not affect quantum information encoded in topological degrees of freedom, unless they involve the creation of topologically nontrivial quasiparticles. Moreover, we show how measurements that detect unintentionally created quasiparticles can be used to control this source of error.


title: Unveiling the submerged part of the iceberg: radio-loud narrow-line Seyfert 1s with SKA

abstract: Narrow-line Seyfert 1 galaxies (NLS1) are active galactic nuclei (AGN) known to have small masses of the central black hole and high accretion rates. NLS1s are generally radio-quiet, but a small part of them (about 7\%) are radio-loud. The recent discovery of powerful relativistic jets in radio-loud NLS1s (RLNLS1s), emitting at high-energy $\gamma$-rays, opened intriguing questions. The observed luminosity of the jet is generally weak, smaller than blazars, although when rescaled for the mass of the central black hole, it becomes of the same order of magnitude of the latter. The weak luminosity, and hence observed flux, resulted in a small number of known RLNLS1. From a recent survey of RLNLS1s, it was found that only 8 out of 42 sources had radio flux density at 1.4 GHz greater than 100 mJy, while 21 out of 42 had flux density smaller than 10 mJy. In addition, given the strong variability at all wavelengths, with present-day facilities RLNLS1s can often only be detected during high activity periods. The Square Kilometer Array (SKA), with its superior sensitivity, will break this limit, allowing us to unveil a relatively unknown population of jetted AGN. We present the results of a study aimed at evaluating the scenario that could emerge after the advent of SKA.


title: A Survey on Security Metrics

abstract: The importance of security metrics can hardly be overstated. Despite the attention that has been paid by the academia, government and industry in the past decades, this important problem stubbornly remains open. In this survey, we present a survey of knowledge on security metrics. The survey is centered on a novel taxonomy, which classifies security metrics into four categories: metrics for measuring the system vulnerabilities, metrics for measuring the defenses, metrics for measuring the threats, and metrics for measuring the situations. The insight underlying the taxonomy is that situations (or outcomes of cyber attack-defense interactions) are caused by certain threats (or attacks) against systems that have certain vulnerabilities (including human factors) and employ certain defenses. In addition to systematically reviewing the security metrics that have been proposed in the literature, we discuss the gaps between the state of the art and the ultimate goals.


title: Shift-Invariant and Sampling Spaces Associated with the Special Affine Fourier Transform

abstract: The Special Affine Fourier Transformation or the SAFT generalizes a number of well known unitary transformations as well as signal processing and optics related mathematical operations. Shift-invariant spaces also play an important role in sampling theory, multiresolution analysis, and many other areas of signal and image processing. Shannon's sampling theorem, which is at the heart of modern digital communications, is a special case of sampling in shift-invariant spaces. Furthermore, it is well known that the Poisson summation formula is equivalent to the sampling theorem and that the Zak transform is closely connected to the sampling theorem and the Poisson summation formula. These results have been known to hold in the Fourier transform domain for decades and were recently shown to hold in the Fractional Fourier transform domain by A. Bhandari and A. Zayed. The main goal of this article is to show that these results also hold true in the SAFT domain. We provide a short, self-contained proof of Shannon's theorem for functions bandlimited in the SAFT domain and then show that sampling in the SAFT domain is equivalent to orthogonal projection of functions onto a subspace of bandlimited basis associated with the SAFT domain. This interpretation of sampling leads to least-squares optimal sampling theorem. Furthermore, we show that this approximation procedure is linked with convolution and semi-discrete convolution operators that are associated with the SAFT domain. We conclude the article with an application of fractional delay filtering of SAFT bandlimited functions.


title: Proof of bijection for combinatorial number system

abstract: Combinatorial number system represents a non-negative natural numbers as sum of binomial coefficients. This paper presents an induction proof that there exists unique representation of every non-negative natural number $m$ as sum of $r$ binomial coefficients.


title: A study of phantom scalar field cosmology using Lie and Noether symmetries

abstract: The paper deals with phantom scalar field cosmology in Einstein gravity. At first using Lie symmetry, the coupling function to the kinetic term and the potential function of the scalar field and the equation of state parameter of the matter field are determined and a simple solution is obtained. Subsequently, Noether symmetry is imposed on the Lagrangian of the system. The symmetry vector is obtained and the potential takes a very general form from which potential using Lie Symmetry can be obtained as a particular case. Then we choose a point transformation $(a,\phi)\rightarrow(u,v)$ such that one of the transformed variables (say u) is a cyclic for the Lagrangian. Using conserved charge (corresponding to the cyclic coordinate) and the constant of motion, solutions are obtained.


title: Gamow Vectors Explain the Shock "Batman" Profile

abstract: The description of shock waves beyond the shock point is a challenge in nonlinear physics. Finding solutions to the global dynamics of dispersive shock waves is not always possible due to the lack of integrability. Here we propose a new method based on the eigenstates (Gamow vectors) of a reversed harmonic oscillator in a rigged Hilbert space. These vectors allow analytical formulation for the development of undular bores of shock waves in a nonlinear nonlocal medium. Experiments by a photothermal induced nonlinearity confirm theoretical predictions: as the undulation period as a function of power and the characteristic quantized decays of Gamow vectors. Our results demonstrate that Gamow vector are a novel and effective paradigm for describing extreme nonlinear phenomena.


title: Young and Millisecond Pulsar GeV Gamma-ray Fluxes from the Galactic Center and Beyond

abstract: Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.


title: Neutrino propagation in the galactic dark matter halo

abstract: Neutrino oscillations are a widely observed and well established phenomenon. It is also well known that deviations with respect to flavor conversion probabilities in vacuum arise due to neutrino interactions with matter. In this work, we analyze the impact of new interactions between neutrinos and the dark matter present in the Milky Way on the neutrino oscillation pattern. The dark matter-neutrino interaction is modeled by using an effective coupling proportional to the Fermi constant $G_F$ with no further restrictions on its flavor structure. For the galactic dark matter profile we consider an homogeneous distribution as well as several density profiles, estimating in all cases the size of the interaction required to get an observable effect at different neutrino energies. Our discussion is mainly focused in the PeV neutrino energy range, to be explored in observatories like IceCube and KM3NeT. The obtained results may be interpreted in terms of a light $\mathcal{O}$(sub-eV--keV) or WIMP-like dark matter particle or as a new interaction with a mediator of $\mathcal{O}$(sub-eV--keV) mass.


title: The second law of quantum thermodynamics as an equality

abstract: We investigate the connection between recent results in quantum thermodynamics and fluctuation relations by adopting a fully quantum mechanical description of thermodynamics. By including a work system whose energy is allowed to fluctuate, we derive a set of equalities which all thermodynamical transitions have to satisfy. This extends the condition for maps to be Gibbs-preserving to the case of fluctuating work, providing a more general characterisation of maps commonly used in the information theoretic approach to thermodynamics. For final states, block diagonal in the energy basis, this set of equalities are necessary and sufficient conditions for a thermodynamical state transition to be possible. The conditions serves as a parent equation which can be used to derive a number of results. These include writing the second law of thermodynamics as an equality featuring a fine-grained notion of the free energy. It also yields a generalisation of the Jarzynski fluctuation theorem which holds for arbitrary initial states, and under the most general manipulations allowed by the laws of quantum mechanics. Furthermore, we show that each of these relations can be seen as the quasi-classical limit of three fully quantum identities. This allows us to consider the free energy as an operator, and allows one to obtain more general and fully quantum fluctuation relations from the information theoretic approach to quantum thermodynamics.


title: NuSTAR unveils a heavily obscured low-luminosity Active Galactic Nucleus in the Luminous Infrared Galaxy NGC 6286

abstract: We report the detection of a heavily obscured Active Galactic Nucleus (AGN) in the luminous infrared galaxy (LIRG) NGC 6286, identified in a 17.5 ks NuSTAR observation. The source is in an early merging stage, and was targeted as part of our ongoing NuSTAR campaign observing local luminous and ultra-luminous infrared galaxies in different merger stages. NGC 6286 is clearly detected above 10 keV and, by including the quasi-simultaneous Swift/XRT and archival XMM-Newton and Chandra data, we find that the source is heavily obscured [$N_{\rm\,H}\simeq (0.95-1.32)\times 10^{24}\rm\,cm^{-2}$], with a column density consistent with being Compton-thick [CT, $\log (N_{\rm\,H}/\rm cm^{-2})\geq 24$]. The AGN in NGC 6286 has a low absorption-corrected luminosity ($L_{2-10\rm\,keV}\sim 3-20\times 10^{41}\rm\,erg\,s^{-1}$) and contributes $\lesssim$1\% to the energetics of the system. Because of its low-luminosity, previous observations carried out in the soft X-ray band ($<10$ keV) and in the infrared did not notice the presence of a buried AGN. NGC 6286 has multi-wavelength characteristics typical of objects with the same infrared luminosity and in the same merger stage, which might imply that there is a significant population of obscured low-luminosity AGN in LIRGs that can only be detected by sensitive hard X-ray observations.


title: Resolving Microlensing Events with Triggered VLBI

abstract: Microlensing events provide a unique capacity to study the stellar remnant population of the Galaxy. Optical microlensing suffers from a near complete degeneracy between the mass, the velocity and the distance. However, a subpopulation of lensed stars, Mira variable stars, are also radio bright, exhibiting strong SiO masers. These are sufficiently bright and compact to permit direct imaging using existing very long baseline interferometers such as the Very Long Baseline Array (VLBA). We show that these events are relatively common, occurring at a rate of $\approx 2~{\rm yr^{\ -1}}$ of which $0.1~{\rm yr^{-1}}$ are associated with Galactic black holes. Features in the associated images, e.g., the Einstein ring, are sufficiently well resolved to fully reconstruct the lens properties, enabling the measurement of mass, distance, and tangential velocity of the lensing object to a precision better than 15%. Future radio microlensing surveys conducted with upcoming radio telescopes combined with modest improvements in the VLBA could increase the rate of Galactic black hole events to roughly 10~${\rm yr}^{-1}$, sufficient to double the number of known stellar mass black holes in a couple years, and permitting the construction of distribution functions of stellar mass black hole properties.


title: Cosmic Reionization On Computers: Numerical and Physical Convergence

abstract: In this paper I show that simulations of reionization performed under the Cosmic Reionization On Computers (CROC) project do converge in space and mass, albeit rather slowly. A fully converged solution (for a given star formation and feedback model) can be determined at a level of precision of about 20%, but such a solution is useless in practice, since achieving it in production-grade simulations would require a large set of runs at various mass and spatial resolutions, and computational resources for such an undertaking are not yet readily available. In order to make progress in the interim, I introduce a weak convergence correction factor in the star formation recipe, which allows one to approximate the fully converged solution with finite resolution simulations. The accuracy of weakly converged simulations approaches a comparable, ~20% level of precision for star formation histories of individual galactic halos and other galactic properties that are directly related to star formation rates, like stellar masses and metallicities. Yet other properties of model galaxies, for example, their HI masses, are recovered in the weakly converged runs only within a factor of two.


title: A physical model for the [CII]-FIR deficit in luminous galaxies

abstract: Observations of ionised carbon at 158 micron ([CII]) from luminous star-forming galaxies at z~0 show that their ratios of [CII] to far infrared (FIR) luminosity are systematically lower than those of more modestly star-forming galaxies. In this paper, we provide a theory for the origin of this so called "[CII] deficit" in galaxies. Our model treats the interstellar medium as a collection of clouds with radially-stratified chemical and thermal properties, which are dictated by the clouds' volume and surface densities, as well as the interstellar radiation and cosmic ray fields to which they are exposed. [CII] emission arises from the outer, HI dominated layers of clouds, and from regions where the hydrogen is H2 but the carbon is predominantly C+. In contrast, the most shielded regions of clouds are dominated by CO and produce little [CII] emission. This provides a natural mechanism to explain the observed [CII]-star formation relation: galaxies' star formation rates are largely driven by the surface densities of their clouds. As this rises, so does the fraction of gas in the CO-dominated phase that produces little [CII] emission. Our model further suggests that the apparent offset in the [CII]-FIR relation for high-z sources compared to those at present epoch may arise from systematically larger gas masses at early times: a galaxy with a large gas mass can sustain a high star formation rate even with relatively modest surface density, allowing copious [CII] emission to coexist with rapid star formation.


title: Resonance locking as the source of rapid tidal migration in the Jupiter and Saturn moon systems

abstract: The inner moons of Jupiter and Saturn migrate outwards due to tidal energy dissipation within the planets, the details of which remain poorly understood. We demonstrate that resonance locking between moons and internal oscillation modes of the planet can produce rapid tidal migration. Resonance locking arises due to the internal structural evolution of the planet and typically produces an outward migration rate comparable to the age of the solar system. Resonance locking predicts a similar migration timescale but a different effective tidal quality factor $Q$ governing the migration of each moon. It also predicts nearly constant migration timescales a function of semi-major axis, such that effective $Q$ values were larger in the past. Recent measurements of Jupiter and Saturn's moon systems find effective $Q$ values that are smaller than expected (and are different between moons), and which correspond to migration timescales of $\sim$10 Gyr. If confirmed, the measurements are broadly consistent with resonance locking as the dominant source of tidal dissipation in Jupiter and Saturn. Resonance locking also provides solutions to several problems posed by current measurements: it naturally explains the exceptionally small $Q$ governing Rhea's migration, it allows the large heating rate of Enceladus to be achieved in an equilibrium eccentricity configuration, and it resolves evolutionary problems arising from present-day migration/heating rates.


title: Momentum Fractionation on Superstrata

abstract: Superstrata are bound states in string theory that carry D1, D5, and momentum charges, and whose supergravity descriptions are parameterized by arbitrary functions of (at least) two variables. In the D1-D5 CFT, typical three-charge states reside in high-degree twisted sectors, and their momentum charge is carried by modes that individually have fractional momentum. Understanding this momentum fractionation holographically is crucial for understanding typical black-hole microstates in this system. We use solution-generating techniques to add momentum to a multi-wound supertube and thereby construct the first examples of asymptotically-flat superstrata. The resulting supergravity solutions are horizonless and smooth up to well-understood orbifold singularities. Upon taking the AdS3 decoupling limit, our solutions are dual to CFT states with momentum fractionation. We give a precise proposal for these dual CFT states. Our construction establishes the very nontrivial fact that large classes of CFT states with momentum fractionation can be realized in the bulk as smooth horizonless supergravity solutions.


title: RadioAstron Observations of the Quasar 3C273: a Challenge to the Brightness Temperature Limit

abstract: Inverse Compton cooling limits the brightness temperature of the radiating plasma to a maximum of $10^{11.5}$ K. Relativistic boosting can increase its observed value, but apparent brightness temperatures much in excess of $10^{13}$ K are inaccessible using ground-based very long baseline interferometry (VLBI) at any wavelength. We present observations of the quasar 3C273, made with the space VLBI mission RadioAstron on baselines up to 171,000 km, which directly reveal the presence of angular structure as small as 26 $\mu$as (2.7 light months) and brightness temperature in excess of $10^{13}$ K. These measurements challenge our understanding of the non-thermal continuum emission in the vicinity of supermassive black holes and require a much higher Doppler factor than what is determined from jet apparent kinematics.


title: Inevitable Power-law Behavior of Isolated Many-Body Quantum Systems and How It Anticipates Thermalization

abstract: Despite being ubiquitous, out-of-equilibrium quantum systems are much less understood than systems at equilibrium. Progress in the field has benefited from a symbiotic relationship between theoretical studies and new experiments on coherent dynamics. The present work strengthens this connection by providing a general picture of the relaxation process of isolated lattice many-body quantum systems that are routinely studied in experiments with cold atoms, ions traps, and nuclear magnetic resonance. We show numerically and analytically that the long-time decay of the probability for finding the system in its initial state necessarily shows a power-law behavior $\propto t^{ - \gamma }$. This happens independently of the details of the system, such as integrability, level repulsion, and the presence or absence of disorder. Information about the spectrum, the structure of the initial state, and the number of particles that interact simultaneously is contained in the value of $\gamma$. From it, we can anticipate whether the initial state will or will not thermalize.


title: Transport and mixing of r-process elements in neutron star binary merger blast waves

abstract: The r-process nuclei are robustly synthesized in the material ejected during a neutron star binary merger (NSBM), as tidal torques transport angular momentum and energy through the outer Lagrange point in the form of a vast tidal tail. If NSBM are indeed solely responsible for the solar system r- process abundances, a galaxy like our own would require to host a few NSBM per million years, with each event ejecting, on average, about 5x10^{-2} M_sun of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than in ordinary supernovae, NSBM deposit a comparable amount of energy into the interstellar medium (ISM). In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution in a cosmological context, we perform a suite of three-dimensional hydrodynamic simulations with optically-thin radiative cooling of isolated NSBM ejecta expanding in environments with conditions adopted from Milky Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution that results from thermal instability in atomic gas is remarkably similar to that of a standard supernova blast wave. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBM blast waves. Among other quantities, we examine the radius, time, mass and kinetic energy content of the NSBM remnant at shell formation as well as the momentum injected to the ISM. We find that the shell formation epoch is attained when the swept-up mass is about 10^3 M_sun, at this point the mass fraction of r-process material is drastically enhanced up to two orders of magnitude in relation to a solar metallicity ISM.


title: On the Possibility of Non-Geodesic Motion of the Massless Spinning Top

abstract: The motion of spinning massless particles in gravitationally curved backgrounds is revisited by considering new types of constraints. Those constraints guarantee zero mass ($P_\mu P^\mu=0$) and they allow for the possibility of trajectories which are not simply null geodesics. To exemplify this previously unknown possibility, the equations of motion are solved for radial motion in Schwarzschild background. It is found that the particle experiences a spin-induced energy shift, which is proportional to the Hawking temperature of the black hole background.