title: Length of Perverse Sheaves on Hyperplane Arrangements

abstract: In this article we address the length of perverse sheaves arising as direct images of rank one local systems on complements of hyperplane arrangements. In the case of a cone over an essential line arrangement with at most triple points, we provide combinatorial formulas for these lengths. As by-products, we also obtain in this case combinatorial formulas for the intersection cohomology Betti numbers of rank one local systems on the complement with same monodromy around the planes.


title: Non-Gaussian and Gottesman-Kitaev-Preskill state preparation by photon catalysis

abstract: Continuous-variable quantum-computing (CVQC) is the most scalable implementation of QC to date but requires non-Gaussian resources to allow exponential speedup and quantum correction, using error encoding such as Gottesman-Kitaev-Preskill (GKP) states. However, GKP state generation is still an experimental challenge. We show theoretically that photon catalysis, the interference of coherent states with single-photon states followed by photon-number-resolved detection, is a powerful enabler for non-Gaussian quantum state engineering such as exactly displaced single-photon states and $M$-symmetric superpositions of squeezed vacuum (SSV), including squeezed cat states ($M=2$). By including photon-counting based state breeding, we demonstrate the potential to enlarge SSV states and produce GKP states.


title: Topology of Hybrid Analytifications

abstract: We investigate the topological properties of Berkovich analytifications over hybrid fields, that is a field equipped with the maximum of its native norm and the trivial norm. We prove that the analytification of the affine line or of a smooth projective curve over a countable Archimedean hybrid field is contractible, and show that it can be non-contractible when the field is uncountable. Further, we prove that the analytification of affine space over a non-Archimedean hybrid field or over a discrete valuation ring is contractible. As an application, we show that the Berkovich affine line over the ring of integers of a number field is contractible.


title: Local differential privacy: Elbow effect in optimal density estimation and adaptation over Besov ellipsoids

abstract: We address the problem of non-parametric density estimation under the additional constraint that only privatised data are allowed to be published and available for inference. For this purpose, we adopt a recent generalisation of classical minimax theory to the framework of local $\alpha$-differential privacy and provide a lower bound on the rate of convergence over Besov spaces $B^s_{pq}$ under mean integrated $\mathbb L^r$-risk. This lower bound is deteriorated compared to the standard setup without privacy, and reveals a twofold elbow effect. In order to fulfil the privacy requirement, we suggest adding suitably scaled Laplace noise to empirical wavelet coefficients. Upper bounds within (at most) a logarithmic factor are derived under the assumption that $\alpha$ stays bounded as $n$ increases: A linear but non-adaptive wavelet estimator is shown to attain the lower bound whenever $p \geq r$ but provides a slower rate of convergence otherwise. An adaptive non-linear wavelet estimator with appropriately chosen smoothing parameters and thresholding is shown to attain the lower bound within a logarithmic factor for all cases.


title: Breathing chromium spinels: a showcase for a variety of pyrochlore Heisenberg Hamiltonians

abstract: We address the long-standing problem of the microscopic origin of the richly diverse phenomena in the chromium breathing pyrochlore material family. Combining electronic structure and renormalization group techniques we resolve the magnetic interactions and analyze their reciprocal-space susceptibility. We show that the physics of these materials is principally governed by long-range Heisenberg Hamiltonian interactions, a hitherto unappreciated fact. Our calculations uncover that in these isostructural compounds, the choice of chalcogen triggers a proximity of the materials to classical spin liquids featuring degenerate manifolds of wave-vectors of different dimensions: A Coulomb phase with three-dimensional degeneracy for LiInCr4O8 and LiGaCr4O8, a spiral spin liquid with two-dimensional degeneracy for CuInCr4Se8 and one-dimensional line degeneracies characteristic of the face-centered cubic antiferromagnet for LiInCr4S8, LiGaCr4S8 and CuInCr4S8. The surprisingly complex array of prototypical pyrochlore behaviors we discovered in chromium spinels may inspire studies of transition paths between different semi-classical spin liquids by doping or pressure.


title: Uniform boundedness for Brauer groups of forms in positive characteristic

abstract: Let $k$ be a finitely generated field of characteristic $p>0$ and $X$ a smooth and proper scheme over $k$. Recent works of Cadoret, Hui and Tamagawa show that, if $X$ satisfies the $\ell$-adic Tate conjecture for divisors for every prime $\ell\neq p$, the Galois invariant subgroup $Br(X_{\overline k})[p']^{\pi_1(k)}$ of the prime-to-$p$ torsion of the geometric Brauer group of $X$ is finite. The main result of this note is that, for every integer $d\geq 1$, there exists a constant $C:=C(X,d)$ such that for every finite field extension $k \subseteq k'$ with $[k':k]\leq d$ and every $(\overline k/k')$-form $Y$ of $X$ one has $|(Br(Y\times_{k'}\overline k)[p']^{\pi_1(k')}|\leq C$. The theorem is a consequence of general results on forms of compatible systems of $\pi_1(k)$-representations and it extends to positive characteristic a recent result of Orr and Skorobogatov in characteristic zero.


title: A Deep Learning based approach to VM behavior identification in cloud systems

abstract: Cloud computing data centers are growing in size and complexity to the point where monitoring and management of the infrastructure become a challenge due to scalability issues. A possible approach to cope with the size of such data centers is to identify VMs exhibiting a similar behavior. Existing literature demonstrated that clustering together VMs that show a similar behavior may improve the scalability of both monitoring andmanagement of a data center. However, available techniques suffer from a trade-off between accuracy and time to achieve this result. Throughout this paper we propose a different approach where, instead of an unsupervised clustering, we rely on classifiers based on deep learning techniques to assigna newly deployed VMs to a cluster of already-known VMs. The two proposed classifiers, namely DeepConv and DeepFFT use a convolution neural network and (in the latter model) exploits Fast Fourier Transformation to classify the VMs. Our proposal is validated using a set of traces describing the behavior of VMs from a realcloud data center. The experiments compare our proposal with state-of-the-art solutions available in literature, demonstrating that our proposal achieve better performance. Furthermore, we show that our solution issignificantly faster than the alternatives as it can produce a perfect classification even with just a few samples of data, making our proposal viable also toclassify on-demand VMs that are characterized by a short life span.


title: O-GAN: Extremely Concise Approach for Auto-Encoding Generative Adversarial Networks

abstract: In this paper, we propose Orthogonal Generative Adversarial Networks (O-GANs). We decompose the network of discriminator orthogonally and add an extra loss into the objective of common GANs, which can enforce discriminator become an effective encoder. The same extra loss can be embedded into any kind of GANs and there is almost no increase in computation. Furthermore, we discuss the principle of our method, which is relative to the fully-exploiting of the remaining degrees of freedom of discriminator. As we know, our solution is the simplest approach to train a generative adversarial network with auto-encoding ability.


title: Escaping Saddle Points with the Successive Convex Approximation Algorithm

abstract: Optimizing non-convex functions is of primary importance in the vast majority of machine learning algorithms. Even though many gradient descent based algorithms have been studied, successive convex approximation based algorithms have been recently empirically shown to converge faster. However, such successive convex approximation based algorithms can get stuck in a first-order stationary point. To avoid that, we propose an algorithm that perturbs the optimization variable slightly at the appropriate iteration. In addition to achieving the same convergence rate results as the non-perturbed version, we show that the proposed algorithm converges to a second order stationary point. Thus, the proposed algorithm escapes the saddle point efficiently and does not get stuck at the first order saddle points.


title: Sparkling saddle loops of vector fields on surfaces

abstract: An orientation-preserving non-contractible separatrix loop of a hyperbolic saddle of a vector field on a two-dimensional surface may be accumulated by a separatrix of the same saddle. We study the unfolding of such loops in generic one-parameter families of vector fields as a semi-local bifurcation. As a byproduct, we construct a countable family of pairwise non-equivalent germs of bifurcation diagrams that appear in locally generic one-parameter families.


title: Coherent Multi-Transducer Ultrasound Imaging

abstract: An extended aperture has the potential to greatly improve ultrasound imaging performance. This work extends the effective aperture size by coherently compounding the received radio frequency data from multiple transducers. A framework is developed in which an ultrasound imaging system consisting of $N$ synchronized matrix arrays, each with partly shared field of view, take turns to transmit plane waves. Only one individual transducer transmits at each time while all $N$ transducers simultaneously receive. The subwavelength localization accuracy required to combine information from multiple transducers is achieved without the use of any external tracking device. The method developed in this study is based on the study of the backscattered echoes received by the same transducer and resulting from a targeted scatterer point in the medium insonated by the multiple ultrasound probes of the system. The current transducer locations along with the speed of sound in the medium are deduced by optimizing the cross-correlation between these echoes. The method is demonstrated experimentally in 2-D using ultrasound point and anechoic lesion phantoms and a first demonstration of a free-hand experiment is also shown. Results demonstrate that the coherent multi-transducer imaging has the potential to improve ultrasound image quality, improving resolution and target detectability. Lateral resolution, contrast and contrast-to-noise ratio improved from 0.67 mm, -6.708 dB and 0.702, respectively, when using a single probe, to 0.18 mm, -7.251 dB and 0.721 in the coherent multi-transducer imaging case.


title: Do BL Lac Objects and FR I radiogalaxies inhabit the same galaxy environment?

abstract: We investigate the environments of galaxies around BL Lac objects and FR I radiogalaxies, the alleged parent populations of misaligned sources. We compare the environment of a sample of 50 BL Lac objects at 0.1<z<0.33 with that of a sample of 90 FR I galaxies at 0.1<z< 0.15. The galaxy environment is estimated using the SDSS images in the i-band. We find that the galaxy excess density within 0.5 Mpc around of FR I radiogalaxies is a factor ~ 2 larger than that around BL Lacs. This implies a reconsideration of the parent population of BL Lac objects.


title: Temperature-dependent stability of polytypes and stacking faults in SiC: reconciling theory and experiments

abstract: The relative stability of SiC polytypes, changing with temperature, has been considered a paradox for about thirty years, due to discrepancies between theory and experiments. Based on ab-initio calculations including van der Waals corrections, a temperature-dependent polytypic diagram consistent with the experimental observations is obtained. Results are easily interpreted based on the influence of the hexagonality on both cohesive energy and entropy. Temperature-dependent stability of stacking faults is also analyzed and found to be in agreement with experimental evidences. Our results suggest that lower temperatures during SiC crystal deposition are advantageous in order to reduce ubiquitous stacking faults in SiC-based power devices.


title: Returning functions with closed graph are continuous

abstract: A function $f:X\to \mathbb R$ defined on a topological space $X$ is called returning if for any point $x\in X$ there exists a positive real number $M_x$ such that for every path-connected subset $C_x\subset X$ containing the point $x$ and any $y\in C_x\setminus\{x\}$ there exists a point $z\in C_x\setminus\{x,y\}$ such that $|f(z)|\le \max\{M_x,|f(y)|\}$. A topological space $X$ is called path-inductive if a subset $U\subset X$ is open if and only if for any path $\gamma:[0,1]\to X$ the preimage $\gamma^{-1}(U)$ is open in $[0,1]$. The class of path-inductive spaces includes all first-countable locally path-connected spaces and all sequential locally contractible space. We prove that a function $f:X\to \mathbb R$ defined on a path-inductive space $X$ is continuous if and only of it is returning and has closed graph. This implies that a (weakly) \'Swi\c atkowski function $f:\mathbb R\to\mathbb R$ is continuous if and only if it has closed graph, which answers a problem of Maliszewski, inscibed to Lviv Scottish Book.


title: Goos-Hanchen shift at the reflection of light from the complex structures composed of superconducting and dielectric layers

abstract: The Goos-Hanchen effect of light reflected from sandwich (three-layered) structures composed of a superconducting YBa2Cu3O7 film and two different dielectric films is investigated theoretically. It has been shown that optical anisotropy of YBa2Cu3O7 film, as well as its positions in the three-layer specimen, strongly effects on the lateral shift values. We have shown that, for all positions of the superconducting film in the three-layered structure, variation of temperature makes possible to control the values of the lateral shift of TE-polarized light at the incidence angles close to pseudo-Brewster angles, whereas for TM-polarized light the lateral shift is only significant at grazing incidence.


title: Universal approximations of permutation invariant/equivariant functions by deep neural networks

abstract: In this paper, we develop a theory about the relationship between $G$-invariant/equivariant functions and deep neural networks for finite group $G$. Especially, for a given $G$-invariant/equivariant function, we construct its universal approximator by deep neural network whose layers equip $G$-actions and each affine transformations are $G$-equivariant/invariant. Due to representation theory, we can show that this approximator has exponentially fewer free parameters than usual models.


title: Entanglement types for two-qubit states with real amplitudes

abstract: We study the set of two-qubit pure states with real amplitudes and their geometrical representation in the three-dimensional sphere. In this representation, we show that the maximally entangled states --those locally equivalent to the Bell States --form two disjoint circles perpendicular to each other. We also show that taking the natural Riemannian metric on the sphere, the set of states connected by local gates are equidistant to this pair of circles. Moreover, the unentangled, or so-called product states, are $\pi/4$ units away to the maximally entangled states. This is, the unentangled states are the farthest away to the maximally entangled states. In this way, if we define two states to be equivalent if they are connected by local gates, we have that there are as many equivalent classes as points in the interval $[0,\pi/4]$ with the point $0$ corresponding to the maximally entangled states. The point $\pi/4$ corresponds to the unentangled states which geometrically are described by a torus. Finally, for every $0< d < \pi/4$ the point $d$ corresponds to a disjoint pair of torus. We also show that if a state is $d$ units away from the maximally entangled states, then its entanglement entropy is $S(d) = 1- \log_2 \sqrt{\frac{(1+\sin 2 d)^{1+\sin 2 d}}{(1-\sin 2 d)^{-1+\sin 2 d}}}$. Finally, we also show how this geometrical interpretation allows us to clearly see that any pair of two-qubit states with real amplitudes can be connected with a circuit that only has single-qubit gates and one controlled-Z gate.


title: Efficient all-optical helicity dependent switching of spins in a Pt/Co/Pt film by a dual-pulse excitation

abstract: All-optical helicity dependent switching (AO-HDS), deterministic control of magnetization by circularly polarized laser pulses, allows to efficiently manipulate spins without the need of a magnetic field. However, AO-HDS in ferromagnetic metals so far requires many laser pulses for fully switching their magnetic states. Using a combination of a short, 90-fs linearly polarized pulse and a subsequent longer, 3-ps circularly polarized pulse, we demonstrate that the number of pulses for full magnetization reversal can be reduced to 4 pulse pairs in a single stack of Pt/Co/Pt. The obtained results suggest that the dual-pulse approach is a potential route towards realizing efficient AO-HDS in ferromagnetic metals.


title: A new wavelength calibration for echelle spectrographs using Fabry-Perot etalons

abstract: The study of Earth-mass extrasolar planets via the radial-velocity technique and the measurement of the potential cosmological variability of fundamental constants call for very-high-precision spectroscopy at the level of $\updelta\lambda/\lambda<10^{-9}$. Wavelength accuracy is obtained by providing two fundamental ingredients: 1) an absolute and information-rich wavelength source and 2) the ability of the spectrograph and its data reduction of transferring the reference scale (wavelengths) to a measurement scale (detector pixels) in a repeatable manner. The goal of this work is to improve the wavelength calibration accuracy of the HARPS spectrograph by combining the absolute spectral reference provided by the emission lines of a thorium-argon hollow-cathode lamp (HCL) with the spectrally rich and precise spectral information of a Fabry-P\'erot-based calibration source. On the basis of calibration frames acquired each night since the Fabry-P\'erot etalon was installed on HARPS in 2011, we construct a combined wavelength solution which fits simultaneously the thorium emission lines and the Fabry-P\'erot lines. The combined fit is anchored to the absolute thorium wavelengths, which provide the `zero-point' of the spectrograph, while the Fabry-P\'erot lines are used to improve the (spectrally) local precision. The obtained wavelength solution is verified for auto-consistency and tested against a solution obtained using the HARPS Laser-Frequency Comb (LFC). The combined thorium+Fabry-P\'erot wavelength solution shows significantly better performances compared to the thorium-only calibration. The presented techniques will therefore be used in the new HARPS and HARPS-N pipeline, and will be exported to the ESPRESSO spectrograph.


title: Invariance of immersed Floer cohomology under Lagrangian surgery

abstract: We show that Floer cohomology of an immersed Lagrangian brane is invariant under smoothing of a self-intersection point if the quantum valuation of the weakly bounding cochain vanishes and the Lagrangian has dimension at least two. The chain-level map replaces the two orderings of the self-intersection point with meridianal and longitudinal cells on the handle created by the surgery, and uses a bijection between holomorphic disks developed by Fukaya-Oh-Ohta-Ono. Our result generalizes invariance of potentials for certain Lagrangian surfaces in Dimitroglou-Rizell--Ekholm--Tonkonog, and implies the invariance of Floer cohomology under mean curvature flow with this type of surgery, as conjectured by Joyce.


title: Generative Adversarial Nets for Robust Scatter Estimation: A Proper Scoring Rule Perspective

abstract: Robust scatter estimation is a fundamental task in statistics. The recent discovery on the connection between robust estimation and generative adversarial nets (GANs) by Gao et al. (2018) suggests that it is possible to compute depth-like robust estimators using similar techniques that optimize GANs. In this paper, we introduce a general learning via classification framework based on the notion of proper scoring rules. This framework allows us to understand both matrix depth function and various GANs through the lens of variational approximations of $f$-divergences induced by proper scoring rules. We then propose a new class of robust scatter estimators in this framework by carefully constructing discriminators with appropriate neural network structures. These estimators are proved to achieve the minimax rate of scatter estimation under Huber's contamination model. Our numerical results demonstrate its good performance under various settings against competitors in the literature.


title: MS-TCN: Multi-Stage Temporal Convolutional Network for Action Segmentation

abstract: Temporally locating and classifying action segments in long untrimmed videos is of particular interest to many applications like surveillance and robotics. While traditional approaches follow a two-step pipeline, by generating frame-wise probabilities and then feeding them to high-level temporal models, recent approaches use temporal convolutions to directly classify the video frames. In this paper, we introduce a multi-stage architecture for the temporal action segmentation task. Each stage features a set of dilated temporal convolutions to generate an initial prediction that is refined by the next one. This architecture is trained using a combination of a classification loss and a proposed smoothing loss that penalizes over-segmentation errors. Extensive evaluation shows the effectiveness of the proposed model in capturing long-range dependencies and recognizing action segments. Our model achieves state-of-the-art results on three challenging datasets: 50Salads, Georgia Tech Egocentric Activities (GTEA), and the Breakfast dataset.


title: Performance Analysis of NOMA-based Cooperative Relaying in {\alpha} - {\mu} Fading Channels

abstract: Non-orthogonal multiple access (NOMA) is widely recognized as a potential multiple access technology for efficient radio spectrum utilization in the fifth-generation (5G) wireless communications standard. In this paper, we study the average achievable rate and outage probability of a cooperative relaying system (CRS) based on NOMA (CRS-NOMA) over wireless links governed by the $\alpha$-$\mu$ generalized fading model; here $\alpha$ and $\mu$ designate the nonlinearity and clustering parameters, respectively, of each link. The average achievable rate is represented in closed-form using Meijer's G-function and the extended generalized bivariate Fox's H-function (EGBFHF), and the outage probability is represented using the lower incomplete Gamma function. Our results confirm that the CRS-NOMA outperforms the CRS with conventional orthogonal multiple access (CRS-OMA) in terms of spectral efficiency at high transmit signal-to-noise ratio (SNR). It is also evident from our results that with an increase in the value of the nonlinearity/clustering parameter, the SNR at which the CRS-NOMA outperforms its OMA based counterpart becomes higher. Furthermore, the asymptotic analysis of the outage probability reveals the dependency of the diversity order of each symbol in the CRS-NOMA system on the $\alpha$ and $\mu$ parameters of the fading links.


title: Bulky Hamiltonian isotopies of Lagrangian tori with applications

abstract: We exhibit monotone Lagrangian tori inside the standard symplectic four-dimensional unit ball that become Hamiltonian isotopic to the Clifford torus, i.e.~the standard product torus, only when considered inside a strictly larger ball (they are not even symplectomorphic to a standard torus inside the unit ball). These tori are then used to construct new examples of symplectic embeddings of toric domains into the unit ball which are symplectically knotted in the sense of J.~Gutt and M.~Usher. We also give a characterisation of the Clifford torus inside the ball as well as the projective plane in terms of quantitative considerations; more specifically, we show that a torus is Hamiltonian isotopic to the Clifford torus whenever one can find a symplectic embedding of a sufficiently large ball in its complement.


title: Influence of misfit strain on the Goos-Hanchen shift upon reflection from a magnetic film on a nonmagnetic substrate

abstract: The influence of the misfit strain on the lateral shift (Goos-Hanchen effect) experienced by a near-infrared electromagnetic wave upon reflection from the surface of a bilayer consisting of a magnetic, gyrotropic (i.e., whose permittivity tensor elements depend upon magnetization) yttrium-iron garnet film deposited on a nonmagnetic gadolinium-gallium garnet substrate is investigated theoretically. In the geometry of the transverse magneto-optical Kerr effect, it is shown that the mechanical strain near the geometrical film/substrate interface can induce a significant lateral shift of the beam for incidence angles close to normal incidence, where no shift appears in the absence of strain. Our calculations demonstrate positive as well as negative values of the lateral shift, depending on the incident light polarization and on the film thickness. In contrast to that of the misfit strain, the influence of the magnetization of the gyrotropic film on the lateral shift is more noticeable for a TM- than for a TE-polarized wave.


title: TableBank: A Benchmark Dataset for Table Detection and Recognition

abstract: We present TableBank, a new image-based table detection and recognition dataset built with novel weak supervision from Word and Latex documents on the internet. Existing research for image-based table detection and recognition usually fine-tunes pre-trained models on out-of-domain data with a few thousand human-labeled examples, which is difficult to generalize on real-world applications. With TableBank that contains 417K high quality labeled tables, we build several strong baselines using state-of-the-art models with deep neural networks. We make TableBank publicly available and hope it will empower more deep learning approaches in the table detection and recognition task. The dataset and models are available at \url{https://github.com/doc-analysis/TableBank}.


title: Pyronia: Intra-Process Access Control for IoT Applications

abstract: Third-party code plays a critical role in IoT applications, which generate and analyze highly privacy-sensitive data. Unlike traditional desktop and server settings, IoT devices mostly run a dedicated, single application. As a result, vulnerabilities in third-party libraries within a process pose a much bigger threat than on traditional platforms. We present Pyronia, a fine-grained access control system for IoT applications written in high-level languages. Pyronia exploits developers' coarse-grained expectations about how imported third-party code operates to restrict access to files, devices, and specific network destinations, at the granularity of individual functions. To efficiently protect such sensitive OS resources, Pyronia combines three techniques: system call interposition, stack inspection, and memory domains. This design avoids the need for application refactoring, or unintuitive data flow analysis, while enforcing the developer's access policy at run time. Our Pyronia prototype for Python runs on a custom Linux kernel, and incurs moderate performance overhead on unmodified Python applications.


title: Light Element Discontinuities Suggest an Early Termination of Star Formation in the Globular Cluster NGC 6402 (M14)

abstract: NGC 6402 is among the most massive globular clusters in the Galaxy, but little is known about its detailed chemical composition. Therefore, we obtained radial velocities and/or chemical abundances of 11 elements for 41 red giant branch cluster members using high resolution spectra obtained with the Magellan-M2FS instrument. We find NGC 6402 to be only moderately metal-poor with $\langle$[Fe/H]$\rangle$ = $-$1.13 dex ($\sigma$ = 0.05 dex) and to have a mean heliocentric radial velocity of $-$61.1 km s$^{\rm -1}$ ($\sigma$ = 8.5 km s$^{\rm -1}$). In general, NGC 6402 exhibits mean composition properties that are similar to other inner Galaxy clusters, such as [$\alpha$/Fe] $\sim$ $+$0.3 dex, [Cr,Ni/Fe] $\sim$ 0.0 dex, and $\langle$[La/Eu]$\rangle$ = $-$0.08 dex. Similarly, we find large star-to-star abundance variations for O, Na, Mg, Al, and Si that are indicative of gas that experienced high temperature proton-capture burning. Interestingly, we detect three distinct populations, but also find large gaps in the [O/Fe], [Na/Fe], and [Al/Fe] distributions that may provide the first direct evidence of delayed formation for intermediate composition stars. A qualitative enrichment model is discussed where clusters form stars through an early ($<$ 5$-$10 Myr) phase, which results in first generation and "extreme" composition stars, and a delayed phase ($>$ 40 Myr), which results in the dilution of processed and pristine gas and the formation of intermediate composition stars. For NGC 6402, the missing intermediate composition stars suggest the delayed phase terminated prematurely, and as a result the cluster may uniquely preserve details of the chemical enrichment process.


title: Weak frames in Hilbert C*-modules with application in Gabor analysis

abstract: In the first part of the paper we describe the dual \ell^2(A)^{\prime} of the standard Hilbert C*-module \ell^2(A) over an arbitrary (not necessarily unital) C*-algebra A. When A is a von Neumann algebra, this enables us to construct explicitly a self-dual Hilbert A-module \ell^2_{\text{strong}}(A) that is isometrically isomorphic to \ell^2(A)^{\prime}, which contains \ell^2(A), and whose A-valued inner product extends the original inner product on \ell^2(A). This serves as a concrete realization of a general construction for Hilbert C*-modules over von Neumann algebras introduced by W. Paschke. Then we introduce a concept of a weak Bessel sequence and a weak frame in Hilbert C*-modules over von Neumann algebras. The dual \ell^2(A)^{\prime} is recognized as a suitable target space for the analysis operator. We describe fundamental properties of weak frames such as the correspondence with surjective adjointable operators, the canonical dual, the reconstruction formula, etc; first for self-dual modules and then, working in the dual, for general modules. In the last part of the paper we describe a class of Hilbert C*-modules over L^{\infty}(I), where I is a bounded interval on the real line, that appear naturally in connection with Gabor (i.e. Weyl-Heisenberg) systems. We then demonstrate that Gabor Bessel systems and Gabor frames in L^2(\Bbb R) are in a bijective correspondence with weak Bessel systems and weak frames of translates by a in these modules over L^{\infty}[0,1/b], where a,b>0 are the lattice parameters. In this setting some well known results on Gabor systems are discussed and some new are obtained.


abstract: We apply our abstract gradient inequalities developed by the authors in arXiv:1510.03817 to prove Lojasiewicz--Simon gradient inequalities for the harmonic map energy function using Sobolev spaces which impose minimal regularity requirements on maps between closed, Riemannian manifolds. Our Lojasiewicz--Simon gradient inequalities for the harmonic map energy function generalize those of Kwon (2002), Liu and Yang (2010), Simon (1983, 1985), and Topping (1997).


title: Elusive Longer-Run Impacts of Head Start: Replications Within and Across Cohorts

abstract: Using an additional decade of CNLSY data, this study replicated and extended Deming's (2009) evaluation of Head Start's life-cycle skill formation impacts in three ways. Extending the measurement interval for Deming's adulthood outcomes, we found no statistically significant impacts on earnings and mixed evidence of impacts on other adult outcomes. Applying Deming's sibling comparison framework to more recent birth cohorts born to CNLSY mothers revealed mostly negative Head Start impacts. Combining all cohorts shows generally null impacts on school-age and early adulthood outcomes.


title: The Lustre Storage Architecture

abstract: This lengthy document often referred to as the "Lustre Book", contains a detailed outline of Lustre file system architecture, as it was created between 2001 and 2005, in accordance with the requirements from various users. Now, in 2019, most features have been implemented, but some only recently, and some along different lines of thought.


title: Probing the Emission Mechanism and Magnetic Field of Neutrino Blazars with Multi-Wavelength Polarization Signatures

abstract: The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $\gamma$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $\gamma$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.


title: Maxwell's (D, H ) excitation fields: lessons from permanent magnets

abstract: Macroscopic Maxwellian electrodynamics consists of four field quantities along with electric charges and electric currents. The fields occur in pairs, the primary ones being the electric and magnetic fields (E,B), and the other the excitation fields (D,H). The link between the two pairs of field is provided by constitutive relations, which specify (D,H) in terms of (E,B); this last connection enabling Maxwell's (differential) equations to be combined in a way that supports waves. In this paper we examine the role played by the excitation fields (D,H), showing that they can be regarded as not having a physical existence, and are merely playing a mathematically convenient role. This point of view is made particularly relevant when we consider competing constitutive models of permanent magnets, which although having the same measurable magnetic properties, have startlingly different behaviours for the magnetic excitation field H.


title: Bondi Accretion in the Spherically Symmetric Johannsen-Psaltis Spacetime

abstract: The Johannsen-Psaltis spacetime explicitly violates the no hair theorem. It describes rotating black holes with scalar hair in the form of parametric deviations from the Kerr metric. In principle, black hole solutions in any modified theory of gravity could be written in terms of the Johannsen-Psaltis metric. We study the accretion of gas onto a static limit of this spacetime. We utilise a recently proposed pseudo-Newtonian formulation of the dynamics around arbitrary static, spherically symmetric spacetimes. We obtain a potential that generalises the Paczy\'nski-Wiita potential to the static Johannsen-Psaltis metric. We also perform a fully relativistic analysis of the geodesic equations in the static Johannsen-Psaltis spacetime. We find that positive values of the scalar hair parameter, $\epsilon_{3}$, lower the accretion rate and vice versa. Similarly, positive (negative) values of $\epsilon_{3}$ reduce (increase) the gravitational acceleration of radially infalling massive particles.


title: Learning Exploration Policies for Navigation

abstract: Numerous past works have tackled the problem of task-driven navigation. But, how to effectively explore a new environment to enable a variety of down-stream tasks has received much less attention. In this work, we study how agents can autonomously explore realistic and complex 3D environments without the context of task-rewards. We propose a learning-based approach and investigate different policy architectures, reward functions, and training paradigms. We find that the use of policies with spatial memory that are bootstrapped with imitation learning and finally finetuned with coverage rewards derived purely from on-board sensors can be effective at exploring novel environments. We show that our learned exploration policies can explore better than classical approaches based on geometry alone and generic learning-based exploration techniques. Finally, we also show how such task-agnostic exploration can be used for down-stream tasks. Code and Videos are available at: https://sites.google.com/view/exploration-for-nav.


title: Fermi gas throughout the BCS-BEC crossover: Comparative study of t-matrix approaches with various degrees of self-consistency

abstract: The diagrammatic t-matrix approximation has often been adopted to describe a dilute Fermi gas. This approximation, originally considered by Galitskii for a repulsive inter-particle interaction (Galitskii-1958), has later been widely utilized for an attractive Fermi gas to describe the BCS-BEC crossover from strongly overlapping Cooper pairs in weak coupling to non-overlapping composite bosons in strong coupling. Several variants of the $t$-matrix approximation have been considered in the literature, which are distinguished by the degree of self-consistency allowed in the building blocks of the diagrammatic structure. Here, we perform a systematic and comparative study of all possible variants on the degree of self-consistency for the $t$-matrix approximation in an attractive Fermi gas, which enables us to confront their outcomes for thermodynamic and dynamical quantities on the same footing in an unbiased way. For definiteness, only the normal phase above the superfluid critical temperature is considered. The dispute that can be raised in this context, about the adequateness of introducing progressive degrees of self-consistency over and above the non-self-consistent $t$-matrix approximation for an attractive Fermi gas, parallels the recent interest in the literature on assessing the importance of various degrees of self-consistency in the context of semiconductors and insulators.


title: Optimal Battery Energy Storage Placement for Transient Voltage Stability Enhancement

abstract: A placement problem for multiple Battery Energy Storage System (BESS) units is formulated towards power system transient voltage stability enhancement in this paper. The problem is solved by the Cross-Entropy (CE) optimization method. A simulation-based approach is adopted to incorporate higher-order dynamics and nonlinearities of generators and loads. The objective is to maximize the voltage stability index, which is set up based on certain grid-codes. Formulations of the optimization problem are then discussed. Finally, the proposed approach is implemented in MATLAB/DIgSILENT and tested on the New England 39-Bus system. Results indicate that installing BESS units at the optimized location can alleviate transient voltage instability issue compared with the original system with no BESS. The CE placement algorithm is also compared with the classic PSO (Particle Swarm Optimization) method, and its superiority is demonstrated in terms of fewer iterations for convergence with better solution qualities.


title: Cyclotomic Coincidences

abstract: In this paper, we show that if $m$ and $n$ are distinct positive integers and $x$ is a nonzero real number with $\Phi_m(x)=\Phi_n(x)$, then $\frac{1}{2}<|x|<2$ except when $\{m,n\}=\{2,6\}$ and $x=2$. We also observe that 2 appears to be the largest limit point of the set of values of $x$ for which $\Phi_m(x)=\Phi_n(x)$ for some $m\neq n$.


title: First Results on 3D Pixel Sensors Interconnected to the RD53A Readout Chip after Irradiation to $1\times$$10^{16}$neq cm$^{-2}$

abstract: Results obtained with 3D columnar pixel sensors bump-bonded to the RD53A prototype readout chip are reported. The interconnected modules have been tested in a hadron beam before and after irradiation to a fluence of about $1\times$$10^{16}$neq cm$^{-2}$ (1MeV equivalent neutrons). All presented results are part of the CMS R&D activities in view of the pixel detector upgrade for the High Luminosity phase of the LHC at CERN (HL-LHC). A preliminary analysis of the collected data shows hit detection efficiencies around 97% measured after proton irradiation.


title: Complexity Results and Algorithms for Bipolar Argumentation

abstract: Bipolar Argumentation Frameworks (BAFs) admit several interpretations of the support relation and diverging definitions of semantics. Recently, several classes of BAFs have been captured as instances of bipolar Assumption-Based Argumentation, a class of Assumption-Based Argumentation (ABA). In this paper, we establish the complexity of bipolar ABA, and consequently of several classes of BAFs. In addition to the standard five complexity problems, we analyse the rarely-addressed extension enumeration problem too. We also advance backtracking-driven algorithms for enumerating extensions of bipolar ABA frameworks, and consequently of BAFs under several interpretations. We prove soundness and completeness of our algorithms, describe their implementation and provide a scalability evaluation. We thus contribute to the study of the as yet uninvestigated complexity problems of (variously interpreted) BAFs as well as of bipolar ABA, and provide the lacking implementations thereof.


title: Structurally stable families of periodic solutions in sweeping processes of networks of elastoplastic springs

abstract: Networks of elastoplastic springs (elastoplastic systems) have been linked to differential equations with polyhedral constraints in the pioneering paper by Moreau (1974). Periodic loading of an elastoplastic system, therefore, corresponds to a periodic motion of the polyhedral constraint. According to Krejci (1996), every solution of a sweeping process with a periodically moving constraint asymptotically converges to a periodic orbit. Understanding whether such an asymptotic periodic orbit is unique or there can be an entire family of asymptotic periodic orbits (that form a periodic attractor) has been an open problem since then. Since suitable small perturbation of a polyhedral constraint seems to be always capable to destroy a potential family of periodic orbits, it is expected that none of potential periodic attractor is structurally stable. In the present paper we give a simple example to prove that even though the periodic attractor (of non-stationary periodic solutions) can be destroyed by little perturbation of the moving constraint, the periodic attractor resists perturbations of the physical parameters of the mechanical model (i.e. the parameters of the network of elastoplastic springs).


title: Dealing with Qualitative and Quantitative Features in Legal Domains

abstract: In this work, we enrich a formalism for argumentation by including a formal characterization of features related to the knowledge, in order to capture proper reasoning in legal domains. We add meta-data information to the arguments in the form of labels representing quantitative and qualitative data about them. These labels are propagated through an argumentative graph according to the relations of support, conflict, and aggregation between arguments.


title: Flavor Quantum Dots and Artificial Quark Model in Transition Metal Dichalcogenides

abstract: We show that the triply degenerate Q valleys in few-layer transition metal dichalcogenides provide a unique platform for exploring the rare flavor SU(3) symmetry in quantum dot geometry. The single and double dots are reminiscent of the quark model and eightfold way, and their many-body triplets and octets may be regarded as artificial quarks and hadrons. For the artificial quark transistor, each level hosts one central and two side Coulomb peaks of irrational height ratios, and flavor Kondo effects occur at 1/3 and 2/3 fillings with fractional conductance quantization in the unitary limit.


title: Augmented Reality Prosthesis Training Setup for Motor Skill Enhancement

abstract: Adjusting to amputation can often time be difficult for the body. Post-surgery, amputees have to wait for up to several months before receiving a properly fitted prosthesis. In recent years, there has been a trend toward quantitative outcome measures. In this paper, we developed the augmented reality (AR) version of one such measure, the Prosthetic Hand Assessment Measure (PHAM). The AR version of the PHAM - HoloPHAM, offers amputees the advantage to train with pattern recognition, at their own time and convenience, pre- and post-prosthesis fitting. We provide a rigorous analysis of our system, focusing on its ability to simulate reach, grasp, and touch in AR. Similarity of motion joint dynamics for reach in physical and AR space were compared, with experiments conducted to illustrate how depth in AR is perceived. To show the effectiveness and validity of our system for prosthesis training, we conducted a 10-day study with able-bodied subjects (N = 3) to see the effect that training on the HoloPHAM had on other established functional outcome measures. A washout phase of 5 days was incorporated to observe the effect without training. Comparisons were made with standardized outcome metrics, along with the progression of kinematic variability over time. Statistically significant (p<0.05) improvements were observed between pre- and post-training stages. Our results show that AR can be an effective tool for prosthesis training with pattern recognition systems, fostering motor learning for reaching movement tasks, and paving the possibility of replacing physical training.


title: PDP: A General Neural Framework for Learning Constraint Satisfaction Solvers

abstract: There have been recent efforts for incorporating Graph Neural Network models for learning full-stack solvers for constraint satisfaction problems (CSP) and particularly Boolean satisfiability (SAT). Despite the unique representational power of these neural embedding models, it is not clear how the search strategy in the learned models actually works. On the other hand, by fixing the search strategy (e.g. greedy search), we would effectively deprive the neural models of learning better strategies than those given. In this paper, we propose a generic neural framework for learning CSP solvers that can be described in terms of probabilistic inference and yet learn search strategies beyond greedy search. Our framework is based on the idea of propagation, decimation and prediction (and hence the name PDP) in graphical models, and can be trained directly toward solving CSP in a fully unsupervised manner via energy minimization, as shown in the paper. Our experimental results demonstrate the effectiveness of our framework for SAT solving compared to both neural and the state-of-the-art baselines.


title: X-rays from Galaxies

abstract: This is a paper to appear as a book chapter in The Chandra X-ray Observatory: Exploring the high energy universe, Ed B Wilkes and W Tucker (Bristol: IOP Publishing Ltd) AAS-IOP ebooks It reviews the results of the observations of galaxies with the Chandra X-ray Observatory to-date, including: populations of X-ray binaries and their evolution; the hot gaseous component; hidden AGNs, and AGN-ISM interaction in nearby galaxies.


title: The knots that lie above all shadows

abstract: We show that for each even integer $m\ge 2$, every reduced shadow with sufficiently many crossings is a shadow of a torus knot T(2,m+1), or of a twist knot $T_m$, or of a connected sum of $m$ trefoil knots.


title: Thermodynamic uncertainty for run-and-tumble type processes

abstract: Thermodynamic uncertainty relations have emerged as universal bounds on current fluctuations in non-equilibrium systems. Here we derive a new bound for a particular class of run-and-tumble type processes using the mathematical framework of renewal-reward theory which can be applied to both Markovian and non-Markovian systems. We demonstrate the results for selected single-particle models as well as a variant of the asymmetric simple exclusion process with collective tumbles. Our bound is relatively tight for a broad parameter regime and only requires knowledge of the statistics of run lengths and the mean entropy production rate of tumbles.


title: Learning Latent Plans from Play

abstract: Acquiring a diverse repertoire of general-purpose skills remains an open challenge for robotics. In this work, we propose self-supervising control on top of human teleoperated play data as a way to scale up skill learning. Play has two properties that make it attractive compared to conventional task demonstrations. Play is cheap, as it can be collected in large quantities quickly without task segmenting, labeling, or resetting to an initial state. Play is naturally rich, covering ~4x more interaction space than task demonstrations for the same amount of collection time. To learn control from play, we introduce Play-LMP, a self-supervised method that learns to organize play behaviors in a latent space, then reuse them at test time to achieve specific goals. Combining self-supervised control with a diverse play dataset shifts the focus of skill learning from a narrow and discrete set of tasks to the full continuum of behaviors available in an environment. We find that this combination generalizes well empirically---after self-supervising on unlabeled play, our method substantially outperforms individual expert-trained policies on 18 difficult user-specified visual manipulation tasks in a simulated robotic tabletop environment. We additionally find that play-supervised models, unlike their expert-trained counterparts, are more robust to perturbations and exhibit retrying-till-success behaviors. Finally, we find that our agent organizes its latent plan space around functional tasks, despite never being trained with task labels. Videos, code and data are available at learning-from-play.github.io


title: Gradient Coding with Clustering and Multi-message Communication

abstract: Gradient descent (GD) methods are commonly employed in machine learning problems to optimize the parameters of the model in an iterative fashion. For problems with massive datasets, computations are distributed to many parallel computing servers (i.e., workers) to speed up GD iterations. While distributed computing can increase the computation speed significantly, the per-iteration completion time is limited by the slowest straggling workers. Coded distributed computing can mitigate straggling workers by introducing redundant computations; however, existing coded computing schemes are mainly designed against persistent stragglers, and partial computations at straggling workers are discarded, leading to wasted computational capacity. In this paper, we propose a novel gradient coding (GC) scheme which allows multiple coded computations to be conveyed from each worker to the master per iteration. We numerically show that the proposed GC with multi-message communication (MMC) together with clustering provides significant improvements in the average completion time (of each iteration), with minimal or no increase in the communication load.


title: MARTINI-based Coarse-grained Model for Poly(alpha-peptoid)s

abstract: In this paper, we present a new coarse-grained (CG) model for poly (alpha-peptoid)s that is compatible with the MARTINI CG FF. In the proposed model, CG poly (alpha-peptoid) is composed by a CG backbone (here we select polysarcosine as the backbone) and side chain beads. The CG model of the backbone (polysarcosine) in a solvent is first developed and then extended to poly (alpha-peptoid)s with different side groups that can be obtained from MARTINI FF. We demonstrate that our CG model has good transferability. For example, the CG potentials for polysarcosine can be transferred to predict hydration free energy of other peptoids. Also, the CG polypeptoid model accurately predicts the radius of gyration over a wide range of chain lengths and the solvation free energy for relatively short peptoid molecules in good solvents. We use the CG model to study sequenced diblock polypeptoid in binary solvent mixtures and compare the results with the experimentally observed coil-globule transition.


title: Spectral Visibility Graphs: Application to Similarity of Harmonic Signals

abstract: Graph theory is emerging as a new source of tools for time series analysis. One promising method is to transform a signal into its visibility graph, a representation which captures many interesting aspects of the signal. Here we introduce the visibility graph for audio spectra and propose a novel representation for audio analysis: the spectral visibility graph degree. Such representation inherently captures the harmonic content of the signal whilst being resilient to broadband noise. We present experiments demonstrating its utility to measure robust similarity between harmonic signals in real and synthesised audio data. The source code is available online.


title: Crowdsourced Behavior-Driven Development: Implementing Microservices through Microtasks

abstract: Key to the effectiveness of crowdsourcing approaches for software engineering is workflow design, describing how complex work is organized into small, relatively independent microtasks. In this paper, we introduce a Behavior-Driven Development (BDD) workflow for accomplishing programming work through self-contained microtasks, implemented as a preconfigured environment called Crowd Microservices. In our approach, a client, acting on behalf of a software team, describes a microservice as a set of endpoints with paths, requests, and responses. A crowd then implements the endpoints, identifying individual endpoint behaviors which they test, implement, and debug, creating new functions and interacting with persistence APIs as needed. To evaluate our approach, we conducted a feasibility study in which a small crowd worked to implement a small ToDo microservice. The crowd created an implementation with only four defects, completing 350 microtasks and implementing 13 functions. We discuss the implications of these findings for incorporating crowdsourced programming contributions into traditional software projects.


title: Multi-tribrackets

abstract: We introduce multi-tribrackets, algebraic structures for region coloring of diagrams of knots and links with different operations at different kinds of crossings. In particular we consider the case of component multi-tribrackets which have different tribracket operations at single-component crossings and multi-component crossings. We provide examples to show that the resulting counting invariants can distinguish links which are not distinguished by the counting invariants associated to the standard tribracket coloring. We reinterpret the results of [11] in terms of multi-tribrackets and consider futuredirections for multi-tribracket theory.


title: Spike-and-Slab Group Lassos for Grouped Regression and Sparse Generalized Additive Models

abstract: We introduce the spike-and-slab group lasso (SSGL) for Bayesian estimation and variable selection in linear regression with grouped variables. We further extend the SSGL to sparse generalized additive models (GAMs), thereby introducing the first nonparametric variant of the spike-and-slab lasso methodology. Our model simultaneously performs group selection and estimation, while our fully Bayes treatment of the mixture proportion allows for model complexity control and automatic self-adaptivity to different levels of sparsity. We develop theory to uniquely characterize the global posterior mode under the SSGL and introduce a highly efficient block coordinate ascent algorithm for maximum a posteriori (MAP) estimation. We further employ de-biasing methods to provide uncertainty quantification of our estimates. Thus, implementation of our model avoids the computational intensiveness of Markov chain Monte Carlo (MCMC) in high dimensions. We derive posterior concentration rates for both grouped linear regression and sparse GAMs when the number of covariates grows at nearly exponential rate with sample size. Finally, we illustrate our methodology through extensive simulations and data analysis.


title: Statistical Guarantees for the Robustness of Bayesian Neural Networks

abstract: We introduce a probabilistic robustness measure for Bayesian Neural Networks (BNNs), defined as the probability that, given a test point, there exists a point within a bounded set such that the BNN prediction differs between the two. Such a measure can be used, for instance, to quantify the probability of the existence of adversarial examples. Building on statistical verification techniques for probabilistic models, we develop a framework that allows us to estimate probabilistic robustness for a BNN with statistical guarantees, i.e., with a priori error and confidence bounds. We provide experimental comparison for several approximate BNN inference techniques on image classification tasks associated to MNIST and a two-class subset of the GTSRB dataset. Our results enable quantification of uncertainty of BNN predictions in adversarial settings.


title: Flame-like Ellerman Bombs and Their Connection to Solar UV Bursts

abstract: Ellerman bombs (EBs) are small-scale intense brightenings in H$\alpha$ wing images, which are generally believed to be signatures of magnetic reconnection events around the temperature minimum region of the solar atmosphere. They have a flame-like morphology when observed near the solar limb. Recent observations from the Interface Region Imaging Spectrograph (IRIS) reveal another type of small-scale reconnection events, termed UV bursts, in the lower atmosphere. Though previous observations have shown a clear coincidence of some UV bursts and EBs, the exact relationship between these two phenomena is still under debate. We investigate the spatial and temporal relationship between flame-like EBs and UV bursts using joint near-limb observations between the 1.6--meter Goode Solar Telescope (GST) and IRIS. In total 161 EBs have been identified from the GST observations, and 20 of them reveal signatures of UV bursts in the IRIS images. Interestingly, we find that these UV bursts have a tendency to appear at the upper parts of their associated flame-like EBs. The intensity variations of most EB-related UV bursts and their corresponding EBs match well. Our results suggest that these UV bursts and EBs are likely formed at different heights during a common reconnection process.


title: Lessons Learned from a Decade of Providing Interactive, On-Demand High Performance Computing to Scientists and Engineers

abstract: For decades, the use of HPC systems was limited to those in the physical sciences who had mastered their domain in conjunction with a deep understanding of HPC architectures and algorithms. During these same decades, consumer computing device advances produced tablets and smartphones that allow millions of children to interactively develop and share code projects across the globe. As the HPC community faces the challenges associated with guiding researchers from disciplines using high productivity interactive tools to effective use of HPC systems, it seems appropriate to revisit the assumptions surrounding the necessary skills required for access to large computational systems. For over a decade, MIT Lincoln Laboratory has been supporting interactive, on-demand high performance computing by seamlessly integrating familiar high productivity tools to provide users with an increased number of design turns, rapid prototyping capability, and faster time to insight. In this paper, we discuss the lessons learned while supporting interactive, on-demand high performance computing from the perspectives of the users and the team supporting the users and the system. Building on these lessons, we present an overview of current needs and the technical solutions we are building to lower the barrier to entry for new users from the humanities, social, and biological sciences.


title: A General size-biased distribution

abstract: We generalize a size-biased distribution related to the Riemann xi function using the work of Ferrar. Some analysis and properties of this more general distribution are offered as well.


title: An Accurate Circuit Model for Photoconductive Antennas

abstract: THz photoconductive antennas are among the most popular devices for THz wave generation and detection. A novel and accurate circuit is proposed to model the performance of such antennas, especially those which are frequency dependant. The suggested equivalent circuit consists of an external source voltage, a source input admittance, a loss resistance related to the electrode loss resistance, a frequency-dependant capacitance associated with capacitive behavior of the gap between the electrodes, and a controlled voltage dependant source representing the screening of photo-carriers. Based on this model, the radiation from a simple frequency dependent antenna, i.e. a dipole antenna, is investigated. Excellent agreement is observed between the obtained results by the proposed equivalent circuit model and those by measurement.


title: Detailed Performance Loss Analysis of Silicon Solar Cells using High-Throughput Metrology Methods

abstract: In this work, novel, high-throughput metrology methods are used to perform a detailed performance loss analysis of approximately 400 industrial crystalline silicon solar cells, all coming from the same production line. The characterization sequence includes a non-destructive transfer length method (TLM) measurement technique featuring circular TLM structures hidden within the busbar region of the cells. It also includes a very fast external quantum efficiency and reflectance measurement technique. More traditional measurements, like illuminated current-voltage, Suns-VOC, and photoluminescence imaging are also used to carry out the loss analysis. The variance of the individual loss parameters and their impact on cell performance are investigated and quantified for this large group of industrial solar cells. Some important correlations between the measured loss parameters are found. The nature of these distributions and correlations provide important insights about loss mechanisms in a cell and help prioritize efforts to optimize the performance of the production line.


title: On the period mod $m$ of polynomially-recursive sequences: a case study

abstract: Polynomially-recursive sequences generally have a periodic behavior mod $m$. In this paper, we analyze the period mod $m$ of a second order polynomially-recursive sequence. The problem originally comes from an enumeration of avoiding pattern permutations and appears to be linked with nice number theory notions (the Carmichael function, Wieferich primes, algebraic integers). We give the mod $a^k$ supercongruences, and generalize these results to a class of recurrences.


title: Multi-epoch Ultraviolet HST Observations of Accreting Low-mass Stars

abstract: Variability is a defining characteristic of young low-mass stars that are still accreting material from their primordial protoplanetary disk. Here we present the largest \textit{HST} variability study of Classical T Tauri stars (CTTS) to date. For 5 of these objects, we obtained a total of 25 spectra with the Space Telescope Imaging Spectrograph (STIS). Mass accretion rates and the fraction of the star covered by accretion columns (i.e., filling factors) were inferred using 1D NLTE physical models whose parameters were fit within a Bayesian framework. On week long timescales, typical changes in the mass accretion rates range up to a factor of $\sim2$, while changes of up to a factor of $\sim5$ are inferred for the filling factors. In addition to this, we observed a possible accretion burst in the transitional disk system GM Aur, and an incident we interpret as a chance alignment of an accretion column and the undisturbed photosphere along our line of sight in the full disk system VW Cha. We also measure correlations between mass accretion rate and line luminosities for use as secondary tracers of accretion. We place our objects in context with recent high-cadence photometric surveys of low-mass star formation regions and highlight the need for more broad-wavelength, contemporaneous data to better understand the physical mechanisms behind accretion variability in CTTS.


title: A New Method to Quantify Environment and Model Ram-Pressure Stripping in N-Body Simulations

abstract: We introduce a Local Background Environment (LBE) estimator that can be measured in and around every galaxy or its dark matter subhalo in high-resolution cosmological simulations. The LBE is designed to capture the influence of environmental effects such as ram-pressure stripping on the formation and evolution of galaxies in semi-analytical models. We define the LBE directly from the particle data within an adaptive spherical shell, and devise a Gaussian mixture estimator to separate background particles from previously unidentified subhalo particles. Analyzing the LBE properties, we find that the LBE of satellite galaxies is not at rest with respect to their host halo, in contrast to typical assumptions. The orientations of the velocities of a subhalo and its LBE are well aligned in the outer infall regions of haloes, but decorrelated near halo center. Significantly, there is no abrupt change in LBE velocity or density at the halo virial radius. This suggests that stripping should also happen beyond this radius. Therefore, we use the time-evolving LBE of galaxies to develop a method to better account for ram-pressure stripping within the Munich semi-analytical model, L-Galaxies. Overall, our new approach results in a significant increase in gas stripping across cosmic time. Central galaxies, as well as satellites beyond the virial radius, can lose a significant fraction of their hot halo gas. As a result, the gas fractions and star formation rates of satellite galaxies are suppressed relative to the fiducial model, although the stellar masses and global stellar mass functions are largely unchanged.


title: Periodic Q-mode modulation in PSR J1825$-$0935 (PSR B1822$-$09)

abstract: PSR J1825$-$0935 (PSR B1822$-$09) switches between radio-quiet (Q-mode) and radio-bright (B-mode) modes. The Q-mode is known to have a periodic fluctuation that modulates both the interpulse and the main pulse with the same period. Earlier investigators argued that the periodic Q-mode modulation is associated with drifting subpulses. We report on single-pulse observations of PSR J1825$-$0935 that were made using the Parkes 64-m radio telescope with a central frequency of 1369 MHz. The high-sensitivity observations revealed that the periodic Q-mode modulation is in fact a periodic longitude-stationary intensity modulation occurring in the interpulse and the main pulse. The fluctuation spectral analysis showed that the modulation period is about $43 P_1$, where $P_1$ is the rotation period of the pulsar. Furthermore, we confirm that the modulation patterns in the interpulse and the main pulse are phase-locked. Specifically, the intensities of the interpulse and the immediately following main pulse are more highly correlated than for the main pulse and interpulse at any other lag. Polarization properties of the strong and weak Q-mode states are different, even for the trailing part of the main pulse which does not show the periodic intensity modulation.


title: Modelling annual and orbital variations in the scintillation of the relativistic binary PSR J1141$-$6545

abstract: We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pulsar timing alone. These include: the orbital inclination $i$; the longitude of the ascending node $\Omega$; and the pulsar system transverse velocity. We use the annual variations to resolve the previous ambiguity in the sense of the inclination angle. Using the correct sense, and a prior probability distribution given by a constraint from pulsar timing ($i=73\pm3^\circ$), we find $\Omega=24.8\pm1.8^\circ$ and we estimate the pulsar distance to be $D=10^{+4}_{-3}$ kpc. This then gives us an estimate of this pulsar's proper motion of $\mu_{\alpha}\cos{\delta}=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $\mu_{\delta}=1.8\pm0.6$ mas yr$^{-1}$ in declination. Finally, we obtain measurements of the spatial structure of the interstellar electron density fluctuations, including: the spatial scale and anisotropy of the diffraction pattern; the distribution of scattering material along the line of sight; and spatial variation in the strength of turbulence from epoch to epoch. We find that the scattering is dominated by a thin screen at a distance of $(0.724\pm0.008)D$, with an anisotropy axial ratio $A_{\rm r} = 2.14\pm0.11$.


title: Hoeffding-Type and Bernstein-Type Inequalities for Right Censored Data

abstract: We present Hoeffding-type and Bernstein-type inequalities for right-censored data. The inequalities bound the difference between an inverse of the probability of censoring weighting (IPCW) estimator and its expectation. We first discuss the asymptotic properties of the estimator and provide conditions for its efficiency. We present standard, data dependent, and uniform Hoeffding-type inequalities. We then present a Bernstein-type inequality. Finally, we show how to apply these inequalities in an empirical risk minimization setting.


title: Enhanced thermal Hall effect in the square-lattice N\'eel state

abstract: Recent experiments on several cuprate compounds have identified an enhanced thermal Hall response in the pseudogap phase. Most strikingly, this enhancement persists even in the undoped system, which challenges our understanding of the insulating parent compounds. To explain these surprising observations, we study the quantum phase transition of a square-lattice antiferromagnet from a confining N\'eel state to a state with coexisting N\'eel and semion topological order. The transition is driven by an applied magnetic field and involves no change in the symmetry of the state. The critical point is described by a strongly-coupled conformal field theory with an emergent global $SO(3)$ symmetry. The field theory has four different formulations in terms of $SU(2)$ or $U(1)$ gauge theories, which are all related by dualities; we relate all four theories to the lattice degrees of freedom. We show how proximity of the confining N\'eel state to the critical point can explain the enhanced thermal Hall effect seen in experiment.


title: Factorization of jet cross sections in heavy-ion collisions

abstract: We propose a new phenomenological approach to establish QCD factorization of jet cross sections in the heavy-ion environment. Starting from a factorization formalism in proton-proton collisions, we introduce medium modified jet functions to capture the leading interaction of jets with the hot and dense QCD medium. A global analysis using a Monte Carlo sampling approach is performed in order to reliably determine the new jet functions from the nuclear modification factor of inclusive jets at the LHC. We find that gluon jets are significantly more suppressed due to the presence of the medium than quark jets. In addition, we observe that the jet radius dependence is directly related to the relative suppression of quark and gluon jets. Our approach may help to improve the extraction of medium properties from data.


title: Is patchy reionization an obstacle in detecting the primordial gravitational wave signal?

abstract: The large-scale CMB B-mode polarization is the direct probe to the low frequency primordial gravitational wave signal. However, unambiguous measurement of this signal requires a precise understanding of the possible contamination. One such potential contamination arises from the patchiness in the spatial distribution of free electrons during the epoch of reionization. We estimate the B-mode power spectrum due to patchy reionization using a combination of \emph{photon-conserving} semi-numerical simulation and analytical calculation, and compare its amplitude with the primordial B-mode signal. For a reionization history which is in agreement with several latest observations, we find that a stronger secondary B-mode polarization signal is produced when the reionization is driven by the sources in massive halos and its amplitude can be comparable to the recombination bump for tensor to scalar ratio $(r) \lesssim 5 \times 10^{-4}$. If contamination from patchy reionization is neglected in the analysis of B-mode polarization data, then for the models of reionization considered in this analysis, we find a maximum bias of about $30\%$ in the value of $r=\,10^{-3}$ when spatial modes between $\ell \in [50, 200]$ are used with a delensing efficiency of $50\%$. The inferred bias from patchy reionization is not a severe issue for the upcoming ground-based CMB experiment Simons Observatory, but can be a potential source of confusion for proposed CMB experiments which target to detect the value of $r< 10^{-3}$. However, this obstacle can be removed by utilizing the difference in the shape of the power spectrum from the primordial signal.


title: Enlarging habitable zones around binary stars in hostile environments

abstract: Habitable zones are regions around stars where large bodies of liquid water can be sustained on a planet or satellite. As many stars form in binary systems with non-zero eccentricity, the habitable zones around the component stars of the binary can overlap and be enlarged when the two stars are at periastron (and less often when the stars are at apastron). We perform N-body simulations of the evolution of dense star-forming regions and show that binary systems where the component stars originally have distinct habitable zones can undergo interactions that push the stars closer together, causing the habitable zones to merge and become enlarged. Occasionally, overlapping habitable zones can occur if the component stars move further apart, but the binary becomes more eccentric. Enlargement of habitable zones happens to 1-2 binaries from an average initial total of 352 in each simulated star-forming region, and demonstrates that dense star-forming regions are not always hostile environments for planet formation and evolution.


title: Black hole mass estimation for Active Galactic Nuclei from a new angle

abstract: The scaling relations between supermassive black holes and their host galaxy properties are of fundamental importance in the context black hole-host galaxy co-evolution throughout cosmic time. In this work, we use a novel algorithm that identifies smooth trends in complex datasets and apply it to a sample of 2,000 type I active galactic nuclei (AGN) spectra. We detect a sequence in emission line shapes and strengths which reveals a correlation between the narrow L([OIII])/L(H$\beta$) line ratio and the width of the broad H$\alpha$. This scaling relation ties the kinematics of the gas clouds in the broad line region to the ionisation state of the narrow line region, connecting the properties of gas clouds kiloparsecs away from the black hole to material gravitationally bound to it on sub-parsec scales. This relation can be used to estimate black hole masses from narrow emission lines only. It therefore enables black hole mass estimation for obscured type II AGN and allows us to explore the connection between black holes and host galaxy properties for thousands of objects, well beyond the local Universe. Using this technique, we present the $M_{\mathrm{BH}} - \sigma$ and $M_{\mathrm{BH}} - M_{*}$ scaling relations for a sample of about 10,000 type II AGN from SDSS. These relations are remarkably consistent with those observed for type I AGN, suggesting that this new method may perform as reliably as the classical estimate used in non-obscured type I AGN. These findings open a new window for studies of black hole-host galaxy co-evolution throughout cosmic time.


title: Implicit Regularization in Over-parameterized Neural Networks

abstract: Over-parameterized neural networks generalize well in practice without any explicit regularization. Although it has not been proven yet, empirical evidence suggests that implicit regularization plays a crucial role in deep learning and prevents the network from overfitting. In this work, we introduce the gradient gap deviation and the gradient deflection as statistical measures corresponding to the network curvature and the Hessian matrix to analyze variations of network derivatives with respect to input parameters, and investigate how implicit regularization works in ReLU neural networks from both theoretical and empirical perspectives. Our result reveals that the network output between each pair of input samples is properly controlled by random initialization and stochastic gradient descent to keep interpolating between samples almost straight, which results in low complexity of over-parameterized neural networks.


title: Statistically-informed deep learning for gravitational wave parameter estimation

abstract: We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell=m=2$ bar mode, $(\omega_R, \omega_I)$. Our neural networks combine a modified $\texttt{WaveNet}$ architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_f, \omega_R, \omega_I)$ of five binary black holes: $\texttt{GW150914}, \texttt{GW170104}, \texttt{GW170814}, \texttt{GW190521}$ and $\texttt{GW190630}$. We use $\texttt{PyCBC Inference}$ to directly compare traditional Bayesian methodologies for parameter estimation with our deep-learning-based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90$\%$ confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 $\texttt{NVIDIA}$ GPU to produce median values and posterior distributions within two milliseconds for each event. This neural network, and a tutorial for its use, are available at the $\texttt{Data and Learning Hub for Science}$.


title: Observation of Chiral Surface Excitons in a Topological Insulator Bi$_2$Se$_3$

abstract: The protected electron states at the boundaries or on the surfaces of topological insulators (TIs) have been the subject of intense theoretical and experimental investigations. Such states are enforced by very strong spin-orbit interaction in solids composed of heavy elements. Here, we study the composite particles -- chiral excitons -- formed by the Coulomb attraction between electrons and holes residing on the surface of an archetypical three-dimensional topological insulator (TI), Bi$_2$Se$_3$. Photoluminescence (PL) emission arising due to recombination of excitons in conventional semiconductors is usually unpolarized because of scattering by phonons and other degrees of freedom during exciton thermalization. On the contrary, we observe almost perfectly polarization-preserving PL emission from chiral excitons. We demonstrate that the chiral excitons can be optically oriented with circularly polarized light in a broad range of excitation energies, even when the latter deviate from the (apparent) optical band gap by hundreds of meVs, and that the orientation remains preserved even at room temperature. Based on the dependences of the PL spectra on the energy and polarization of incident photons, we propose that chiral excitons are made from massive holes and massless (Dirac) electrons, both with chiral spin textures enforced by strong spin-orbit coupling. A theoretical model based on such proposal describes quantitatively the experimental observations. The optical orientation of composite particles, the chiral excitons, emerges as a general result of strong spin-orbit coupling in a 2D electron system. Our findings can potentially expand applications of TIs in photonics and optoelectronics.


title: Vanishing Hall Response of Charged Fermions in a Transverse Magnetic Field

abstract: We study the Hall response of quasi-two-dimensional lattice systems of charged fermions under a weak transverse magnetic field, in the ballistic coherent limit. We identify a setup in which this response vanishes over a wide range of parameters: the paradigmatic "Landauer-B\"uttiker" setup commonly studied for coherent quantum transport, consisting of a strip contacted to biased ideal reservoirs of charges. We show that the effect does not rely on particle-hole symmetry, and is robust to a variety of perturbations including variations of the transverse magnetic field, chemical potential, and temperature. We trace this robustness back to a topological property of the Fermi surface: the number of Fermi points (central charge) of the system. We argue that the mechanism responsible for the vanishing Hall response can operate both in noninteracting and interacting systems, which we verify in concrete examples using density-matrix renormalization group (DMRG) simulations.


title: Renormalization-group study of the many-body localization transition in one dimension

abstract: Using a new approximate strong-randomness renormalization group (RG), we study the many-body localized (MBL) phase and phase transition in one-dimensional quantum systems with short-range interactions and quenched disorder. Our RG is built on those of Zhang $\textit{et al.}$ [1] and Goremykina $\textit{et al.}$ [2], which are based on thermal and insulating blocks. Our main addition is to characterize each insulating block with two lengths: a physical length, and an internal decay length $\zeta$ for its effective interactions. In this approach, the MBL phase is governed by a RG fixed line that is parametrized by a global decay length $\tilde{\zeta}$, and the rare large thermal inclusions within the MBL phase have a fractal geometry. As the phase transition is approached from within the MBL phase, $\tilde{\zeta}$ approaches the finite critical value corresponding to the avalanche instability, and the fractal dimension of large thermal inclusions approaches zero. Our analysis is consistent with a Kosterlitz-Thouless-like RG flow, with no intermediate critical MBL phase.


title: RELICS: Reionization Lensing Cluster Survey

abstract: Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.


title: The Milky Way bar/bulge in proper motions: a 3D view from VIRAC & Gaia

abstract: We have derived absolute proper motions of the entire Galactic bulge region from VIRAC and Gaia. We present these as both integrated on-sky maps and, after isolating standard candle red clump (RC) stars, as a function of distance using RC magnitude as a proxy. These data provide a new global, 3-dimensional view of the Milky Way barred bulge kinematics. We find a gradient in the mean longitudinal proper motion, $<\mu_l^\star>$, between the different sides of the bar, which is sensitive to the bar pattern speed. The split RC has distinct proper motions and is colder than other stars at similar distance. The proper motion correlation map has a quadrupole pattern in all magnitude slices showing no evidence for a separate, more axisymmetric inner bulge component. The line-of-sight integrated kinematic maps show a high central velocity dispersion surrounded by a more asymmetric dispersion profile. $\sigma_{\mu_l} / \sigma_{\mu_b}$ is smallest, $\sim1.1$, near the minor axis and reaches $\sim1.4$ near the disc plane. The integrated $<\mu_b>$ pattern signals a superposition of bar rotation and internal streaming motion, with the near part shrinking in latitude and the far part expanding. To understand and interpret these remarkable data, we compare to a made-to-measure barred dynamical model, folding in the VIRAC selection function to construct mock maps. We find that our model of the barred bulge, with a pattern speed of 37.5 $\mathrm{km \, s^{-1} \, kpc^{-1}}$, is able to reproduce all observed features impressively well. Dynamical models like this will be key to unlocking the full potential of these data.


title: Formation of short-period planets by disk migration

abstract: Protoplanetary disks are thought to be truncated at orbital periods of around 10 days. Therefore, origin of rocky short period planets with $P < 10$ days is a puzzle. We propose that many of these planets may form through the Type-I migration of planets locked into a chain of mutual mean motion resonances. We ran N-body simulations of planetary embryos embedded in a protoplanetary disk. The embryos experienced gravitational scatterings, collisions, disk torques, and dampening of orbital eccentricity and inclination. We then modelled Kepler observations of these planets using a forward model of both the transit probability and the detection efficiency of the Kepler pipeline. We found that planets become locked into long chains of mean motion resonances that migrate in unison. When the chain reaches the edge of the disk, the inner planets are pushed past the edge due to the disk torques acting on the planets farther out in the chain. Our simulated systems successfully reproduce the observed period distribution of short period Kepler planets between 1 and 2 $R_\oplus$. However, we obtain fewer closely packed short period planets than in the Kepler sample. Our results provide valuable insight into the planet formation process, and suggests that resonance locks, migration, and dynamical instabilities play important roles the the formation and evolution of close-in small exoplanets.


title: Sub-second variability in black-hole X-ray binary jets

abstract: In the last 10 years multi-wavelength fast variability studies of low mass X-ray binaries have shown a dramatic development. A key discovery was the detection of O-IR sub-second fluctuations in two black-hole transients, lagging the X-rays by $\approx$0.1 s. This demonstrated how the fluctuations observed in the inflow could be transferred to the jet, allowing therefore also to study in a completely new way the physical processes which take place at the base of the jet. In this paper we review the latest developments of the study of jets with this new approach, focusing on the results obtained with cross-spectral analysis techniques.


title: High-Energy Neutrinos from Blazar Flares and Implications of TXS 0506+056

abstract: Motivated by the observation of a $>290$ TeV muon neutrino by IceCube, coincident with a $\sim$6 month-long $\gamma$-ray flare of the blazar TXS 0506+056, and an archival search which revealed $13 \pm 5$ further, lower-energy neutrinos in the direction of the source in 2014-2015, we discuss the likely contribution of blazars to the diffuse high-energy neutrino intensity, the implications for neutrino emission from TXS 0506+056 based on multi-wavelength observations of the source, and a multi-zone model that allows for sufficient neutrino emission so as to reconcile the multi-wavelength cascade constraints with the neutrino emission seen by IceCube in the direction of TXS 0506+056.


title: Cosmic web anisotropy is the primary indicator of halo assembly bias

abstract: The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $\alpha$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $\alpha$, and the other between $\alpha$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $\beta$) and angular momentum (dimensionless spin $\lambda$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $\alpha$, $b_1$ and any of the internal properties $c\in\{\beta,c_v/a_v,c/a,c_{\rm vir},\lambda\}$ is consistent with $p(\alpha,b_1,c)\simeq p(\alpha)p(b_1|\alpha)p(c|\alpha)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-\alpha$ and $b_1-\alpha$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $\alpha$ does not explain their assembly bias trends. The overarching importance of $\alpha$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias.


title: Transverse kinematics of the Galactic bar-bulge from VVV and Gaia

abstract: We analyse the kinematics of the Galactic bar-bulge using proper motions from the ESO public survey Vista Variables in the Via Lactea (VVV) and the second Gaia data release. Gaia has provided some of the first absolute proper motions within the bulge and the near-infrared VVV multi-epoch catalogue complements Gaia in highly-extincted low-latitude regions. We discuss the relative-to-absolute calibration of the VVV proper motions using Gaia. Along lines of sight spanning $-10<\ell/\,\mathrm{deg}<10$ and $-10<b/\,\mathrm{deg}<5$, we probabilistically model the density and velocity distributions as a function of distance of $\sim45$ million stars. The transverse velocities confirm the rotation signature of the bar seen in spectroscopic surveys. The differential rotation between the double peaks of the magnitude distribution confirms the X-shaped nature of the bar-bulge. Both transverse velocity components increase smoothly along the near-side of the bar towards the Galactic centre, peak at the Galactic centre and decline on the far-side. The anisotropy is $\sigma_\ell/\sigma_b\approx1.1-1.3$ within the bulk of the bar, reducing to $0.9-1.1$ when rotational broadening is accounted for, and exhibits a clear X-shaped signature. The vertex deviation in $\ell$ and $b$ is significant $|\rho_{\ell b}|\lesssim0.2$, greater on the near-side of the bar and produces a quadrupole signature across the bulge indicating approximate radial alignment. We have re-constructed the 3D kinematics from the assumption of triaxiality, finding good agreement with spectroscopic survey results. In the co-rotating frame, we find evidence of bar-supporting x1 orbits and tangential bias in the in-plane dispersion field.


title: The pattern speed of the Milky Way bar from transverse velocities

abstract: We use the continuity equation to derive a method for measuring the pattern speed of the Milky Way's bar/bulge from proper motion data. The method has minimal assumptions but requires complete coverage of the non-axisymmetric component in two of the three Galactic coordinates. We apply our method to the proper motion data from a combination of Gaia DR2 and VISTA Variables in the Via Lactea (VVV) to measure the pattern speed of the bar as $\Omega_\mathrm{p}=(41\pm 3)\,\mathrm{km\,s^{-1}\,kpc^{-1}}$ (where the error is statistical). This puts the corotation radius at $(5.7\pm0.4)\,\mathrm{kpc}$, under the assumptions of the standard peculiar motion of the Sun and the absence of non-axisymmetric streaming in the Solar neighbourhood. The obtained result uses only data on the near-side of the bar which produces consistent measurements of the distance and velocity of the centre of the Galaxy. Addition of the data on the far-side of the bar pulls the pattern speed down to $\Omega_\mathrm{p}=(31\pm 1)\,\mathrm{km\,s^{-1}\,kpc^{-1}}$ but requires a lower transverse velocity for the Galactic centre than observed. This suggests systematics of $5-10\,\mathrm{km\,s^{-1}kpc^{-1}}$ dominate the uncertainty. We demonstrate using a dynamically-formed bar/bulge simulation that even with the limited field of view of the VVV survey our method robustly recovers the pattern speed.


title: No Run Gravity

abstract: Considering the dark energy/gravity landscape if next generation surveys of galaxies, cosmic microwave background radiation, and gravitational waves do not find clear modification of gravity, we develop No Run Gravity as a counterexample to the conclusion that this would imply general relativity with an expansion history described by an equation of state $w(z)$. No Run Gravity is a cubic Horndeski theory with a constant Planck mass, no gravitational slip, and no modification of gravitational waves, but a rich phenomenology beyond $w(z)$. We calculate the evolution of gravitational strength, sound speed, and cosmic growth within the theory and project sensitivities for upcoming DESI redshift space distortion data.


title: Experimentally reducing the quantum measurement back-action in work distributions by a collective measurement

abstract: In quantum thermodynamics, the standard approach to estimate work fluctuations in unitary processes is based on two projective measurements, one performed at the beginning of the process and one at the end. The first measurement destroys any initial coherence in the energy basis, thus preventing later interference effects. In order to decrease this back-action, a scheme based on collective measurements has been proposed in~[PRL 118, 070601 (2017)]. Here, we report its experimental implementation in an optical system. The experiment consists of a deterministic collective measurement on identically prepared two qubits, encoded in the polarisation and path degree of a single photon. The standard two projective measurement approach is also experimentally realized for comparison. Our results show the potential of collective schemes to decrease the back-action of projective measurements, and capture subtle effects arising from quantum coherence.


title: Universal skein theory for group actions

abstract: Given a group action on a finite set, we define the group-action model which consists of tensor network diagrams which are invariant under the group symmetry. In particular, group-action models can be realized as the even part of group-subgroup subfactor planar algebras. Moreover, all group-subgroup subfactor planar algebras arise in this way from transitive actions. In this paper, we provide a universal skein theory for those planar algebras. With the help of this skein theory, we give a positive answer to a question asked by Vaughan Jones in the late nineties.