title: Conservative Discontinuous Galerkin Schemes for Nonlinear Fokker-Planck Collision Operators

abstract: We present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker-Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the `Lenard-Bernstein,' or `Dougherty,' operator. Several novel algorithmic innovations are reported. The concept of weak-equality is introduced and used to define weak-operators to compute primitive moments needed in the updates. Weak-equality is also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integration by parts are used to construct the discrete weak-form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the $v^2$ term is not included in the basis-set used in the discretization. A series of benchmark problems are presented and show that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping. We conclude with an outline of future work, in particular with some indications of how the algorithms presented here can be extended to use the Rosenbluth potentials to compute the drag and diffusion coefficients.


title: Transport coefficients from in medium quarkonium dynamics

abstract: The in medium dynamics of heavy particles are governed by transport coefficients. The heavy quark momentum diffusion coefficient, $\kappa$, is an object of special interest in the literature, but one which has proven notoriously difficult to estimate, despite the fact that it has been computed by weak-coupling methods at next-to-leading order accuracy, and by lattice simulations of the pure SU(3) gauge theory. Another coefficient, $\gamma$, has been recently identified. It can be understood as the dispersive counterpart of $\kappa$. Little is known about $\gamma$. Both $\kappa$ and $\gamma$ are, however, of foremost importance in heavy quarkonium physics as they entirely determine the in and out of equilibrium dynamics of quarkonium in a medium, if the evolution of the density matrix is Markovian, and the motion, quantum Brownian; the medium could be a strongly or weakly coupled plasma. In this paper, using the relation between $\kappa$, $\gamma$ and the quarkonium in medium width and mass shift respectively, we evaluate the two coefficients from existing 2+1 flavor lattice QCD data. The resulting range for $\kappa$ is consistent with earlier determinations, the one for $\gamma$ is the first non-perturbative determination of this quantity.


title: Thermal Heavy Quark Self-Energy from Euclidean Correlators

abstract: Brambilla, Escobedo, Soto, and Vairo have derived an effective description of quarkonium with two parameters; a momentum diffusion term and a real self-energy term. We point out that there is a similar real self-energy term for a single open heavy flavor and that it can be expressed directly in terms of Euclidean electric field correlators along a Polyakov line. This quantity can be directly studied on the lattice without the need for analytical continuation. We show that Minkowski-space calculations of this correlator correspond with the known NLO Euclidean value of the relevant electric field two-point function and that it differs from the real self-energy term for quarkonium.


title: Vanishing of the anchored isoperimetric profile in bond percolation at p c

abstract: We consider the anchored isoperimetric profile of the infinite open cluster, defined for $p > p\_c$, whose existence has been recently proved in [3]. We extend adequately the definition for $p = p\_c$, in finite boxes. We prove a partial result which implies that, if the limit defining the anchored isoperimetric profile at $p\_c$ exists, it has to vanish.


title: Trained Quantization Thresholds for Accurate and Efficient Fixed-Point Inference of Deep Neural Networks

abstract: We propose a method of training quantization thresholds (TQT) for uniform symmetric quantizers using standard backpropagation and gradient descent. Contrary to prior work, we show that a careful analysis of the straight-through estimator for threshold gradients allows for a natural range-precision trade-off leading to better optima. Our quantizers are constrained to use power-of-2 scale-factors and per-tensor scaling of weights and activations to make it amenable for hardware implementations. We present analytical support for the general robustness of our methods and empirically validate them on various CNNs for ImageNet classification. We are able to achieve near-floating-point accuracy on traditionally difficult networks such as MobileNets with less than 5 epochs of quantized (8-bit) retraining. Finally, we present Graffitist, a framework that enables automatic quantization of TensorFlow graphs for TQT (available at https://github.com/Xilinx/graffitist ).


title: Differential Interferometric Signatures of Close Binaries of Supermassive Black Holes in Active Galactic Nuclei

abstract: In the present paper, we explore opportunities of applying the GRAVITY at the Very Large Telescope Interferometry (VLTI) with unprecedented spatial resolution to identify close binaries of supermassive black holes (CB-SMBHs) in active galactic nuclei (AGNs). Each SMBH is assumed to be separately surrounded by their own broad-line regions (BLRs) composed of clouds with virialized motion. Composition of the binary orbital motion and the virial motion of clouds in each BLR determines the projected velocity fields and hence differential phase curves, which are obviously different from that of a single BLR. We calculate emission line profiles and differential phase curves of CB-SMBHs for the GRAVITY. For the simplest case where angular momentums of two BLRs and orbital motion are parallel, a phase plateau generally appears in the phase curves. For other combinations of the angular momentum, the plateau is replaced by new peaks and valleys variously depending on the situations. Given a combination, phase curves are also sensitive to changes of parameters of CB-SMBHs. All these features are easily distinguished from the well-known $S$-shaped phase curves of a single BLR so that the GRAVITY is expected to reveal signals of CB-SMBH from candidates of AGNs. With joint analysis of observations of reverberation mapping campaigns, we can reliably identify CB-SMBHs, and measure their orbital parameters in the meanwhile. This independent measurement of the orbital parameters also has implications to analysis of Pulsar Timing Array (PTA) observations for properties of low-frequency gravitational waves in future.


title: Optimization of Rate Allocation and Power Control for Rate Splitting Multiple Access (RSMA)

abstract: In this paper, the sum-rate maximization problem is studied for wireless networks that use downlink rate splitting multiple access (RSMA). In the considered model, each base station (BS) divides the messages that must be transmitted to its users into a `private' part and a `common' part. Here, the common message is a message that all users want to receive and the private message is a message that is dedicated to only a specific user. The RSMA mechanism enables a BS to adjust the split of common and private messages so as to control the interference by decoding and treating interference as noise and, thus optimizing the data rate of users. To maximize the users' sum-rate, the network can determine the rate allocation for the common message to meet the rate demand, and adjust the transmit power for the private message to reduce the interference. This problem is formulated as an optimization problem whose goal is to maximize the sum-rate of all users. To solve this nonconvex maximization problem, the optimal power used for transmitting the private message to the users is first obtained in closed form for a given rate allocation and common message power. Based on the optimal private message transmission power, the optimal rate allocation is then derived under a fixed common message transmission power. Subsequently, a one-dimensional search algorithm is proposed to obtain the optimal solution of common message transmission power. Simulation results show that the RSMA can achieve up to 15.6\% and 21.5\% gains in terms of data rate compared to non-orthogonal multiple access (NOMA) and orthogonal frequency-division multiple access (OFDMA), respectively.


title: Hankel determinant of second order for some classes of analytic functions

abstract: Let $f$ be analytic in the unit disk $\mathbb D$ and normalized so that $f(z)=z+a_2z^2+a_3z^3+\cdots$. In this paper, we give upper bounds of the Hankel determinant of second order for the classes of starlike functions of order $\alpha$, Ozaki close-to-convex functions and two other classes of analytic functions. Some of the estimates are sharp.


title: Low-Complexity Iterative Receiver for Orthogonal Chirp Division Multiplexing

abstract: This paper proposes a low-complexity iterative receiver for the recently proposed Orthogonal Chirp Division Multiplexing (OCDM) modulation scheme, where we consider a system under frequency-selective channels and constrained to channel state information availability only at the receiver. It has been shown that under these assumptions, OCDM becomes an optimal waveform in terms of performance, i.e., frame error rate (FER), when employing a receiver capable of achieving perfect feedback equalizer (PFE) performance. Thus, this work targets proposing such a receiver for OCDM with low-complexity. Our approach is based on the well accepted minimum mean squared error with parallel interference cancellation (MMSE-PIC), where we derive an approximated equalizer whose complexity is reduced to two fast Fourier transforms (FFTs) per iteration. The FER results reveal that i) the proposed low-complexity receiver perform as good as the original MMSE-PIC, ii) OCDM performs very closely to PFE, and iii) OCDM has approximately 2.5 dB improvement over OFDM.


title: The first 62 AGN observed with SDSS-IV MaNGA -- III: stellar and gas kinematics

abstract: We investigate the effects of Active Galactic Nuclei (AGN) on the gas kinematics of their host galaxies, using MaNGA data for a sample of 62 AGN hosts and 109 control galaxies (inactive galaxies). We compare orientation of the line of nodes (kinematic Position Angle - PA) measured from the gas and stellar velocity fields for the two samples. We found that AGN hosts and control galaxies display similar kinematic PA offsets between gas and stars. However, we note that AGN have larger fractional velocity dispersion $\sigma$ differences between gas and stars [$\sigma_{frac}=(\sigma_{\rm gas}-\sigma_{stars})/\sigma_{\rm stars}$] when compared to their controls, as obtained from the velocity dispersion values of the central (nuclear) pixel (2.5" diameter). The AGN have a median value of $\sigma_{\rm frac}$ of $<\sigma_{frac}>_{\rm AGN}=0.04$, while the the median value for the control galaxies is $<\sigma_{frac}>_{\rm CTR}=-0.23$. 75% of the AGN show $\sigma_{frac}>-0.13$, while 75% of the normal galaxies show $\sigma_{\rm frac}<-0.04$, thus we suggest that the parameter $\sigma_{\rm frac}$ can be used as an indicative of AGN activity. We find a correlation between the [OIII]$\lambda$5007 luminosity and $\sigma_{frac}$ for our sample. Our main conclusion is that the AGN already observed with MaNGA are not powerful enough to produce important outflows at galactic scales, but at 1-2 kpc scales, AGN feedback signatures are always present on their host galaxies.


title: Max-plus Operators Applied to Filter Selection and Model Pruning in Neural Networks

abstract: Following recent advances in morphological neural networks, we propose to study in more depth how Max-plus operators can be exploited to define morphological units and how they behave when incorporated in layers of conventional neural networks. Besides showing that they can be easily implemented with modern machine learning frameworks , we confirm and extend the observation that a Max-plus layer can be used to select important filters and reduce redundancy in its previous layer, without incurring performance loss. Experimental results demonstrate that the filter selection strategy enabled by a Max-plus is highly efficient and robust, through which we successfully performed model pruning on different neural network architectures. We also point out that there is a close connection between Maxout networks and our pruned Max-plus networks by comparing their respective characteristics. The code for reproducing our experiments is available online.

categories: math.ST cs.CV cs.LG cs.NE stat.ML stat.TH


title: Poincare' covariant light-front spectral function and transverse momentum distribution

abstract: In valence approximation the fermion correlator is simply related to the light-front spectral function. Then the leading twist time-reversal even transverse momentum distributions can be explicitly obtained from the light-front wave function of the system and the twist-three distributions are linear combinations of the transverse distributions at leading twist.


title: BotGraph: Web Bot Detection Based on Sitemap

abstract: The web bots have been blamed for consuming large amount of Internet traffic and undermining the interest of the scraped sites for years. Traditional bot detection studies focus mainly on signature-based solution, but advanced bots usually forge their identities to bypass such detection. With increasing cloud migration, cloud providers provide new opportunities for an effective bot detection based on big data to solve this issue. In this paper, we present a behavior-based bot detection scheme called BotGraph that combines sitemap and convolutional neural network (CNN) to detect inner behavior of bots. Experimental results show that BotGraph achieves ~95% recall and precision on 35-day production data traces from different customers including the Bing search engine and several sites.


title: Multi timescale bandwidth profile and its application for burst-aware fairness

abstract: We propose a resource sharing scheme that takes into account the traffic history over several predefined time scales and provides fair resource sharing considering the traffic history. Our concept builds on a simplified version of core-stateless resource sharing, where we only use a few Drop Precedences (DPs). For packet marking we introduce Multi timescale bandwidth profile. Additionally, we provide basic dimensioning concepts for the proposed schema and present its simulation based performance analysis.


title: The Changing Geopolitics in the Arab World: Implications of the 2017 Gulf Crisis for Business

abstract: The international community was caught by surprise on 5 June 2017 when Saudi Arabia, the United Arab Emirates (UAE), Bahrain and Egypt severed diplomatic ties with Qatar, accusing it of destabilizing the region. More than one year after this diplomatic rift, several questions remain unaddressed. This study focuses on the regional business costs of the year-long blockade on Qatar. We split the sample to compare the stock market performances of Qatar and its Middle Eastern neighbors before and after the Saudi-led Qatar boycott. We focus our attention on the conditional volatility process of stock market returns and risks related to financial interconnectedness. We show that the Gulf crisis had the most adverse impact on Qatar together with Saudi Arabia and the UAE. Although not to the same degree as these three countries, Bahrain and Egypt were also harmfully affected. But shocks to the volatility process tend to have short-lasting effects. Moreover, the total volatility spillovers to and from others increase but moderately after the blockade. Overall, the quartet lobbying efforts did not achieve the intended result. Our findings underscore Qatar's economic vulnerability but also the successful resilience strategy of this tiny state. The coordinated diplomatic efforts of Qatar have been able to fight the economic and political embargo.


title: Behaviour of the Stokes operators under domain perturbation

abstract: Depending of the geometry of the domain, one can define --at least-- three different Stokes operators with Dirichlet boundary conditions. We describe how the resolvents of these Stokes operators converge with respect to a converging sequence of domains.


title: Muon capture in nuclei: an ab initio approach based on quantum Monte Carlo methods

abstract: An ab initio quantum Monte Carlo method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A \leq 12$. As a first application of the method, we perform a calculation of the rate in $^4$He in a dynamical framework based on realistic two- and three-nucleon interactions and realistic nuclear charge-changing weak currents. The currents include one- and two-body terms induced by $\pi$- and $\rho$-meson exchange, and $N$-to-$\Delta$ excitation, and are constrained to reproduce the empirical value of the Gamow-Teller matrix element in tritium. We investigate the sensitivity of theoretical predictions to current parametrizations of the nucleon axial and induced pseudoscalar form factors as well as to two-body contributions in the weak currents. The large uncertainties in the measured values obtained from bubble-chamber experiments (carried out over 50 years ago) prevent us from drawing any definite conclusions.


title: Coneat Injective Modules

abstract: A module is called coneat injective if it is injective with respect to all coneat exact sequences. The class of such modues is enveloping and falls properly between injectives and pure injectives. Generalizations of coneat injectivity, like relative coneat injectivity and full invariance of a module in its coneat injective envelope, are studied. Using properties of such classes of modules, we characterize certain types of rings like von Neumann regular and right SF-rings. For instance, R is a right SF-ring if and only if every coneat injective left R-module is injective.


title: Holomorphic functions of exponential type on connected complex Lie groups

abstract: Holomorphic functions of exponential type on a complex Lie group $G$ (introduced by Akbarov) form a locally convex algebra, which is denoted by $\cO_{exp}(G)$. Our aim is to describe the structure of $\cO_{exp}(G)$ in the case when $G$ is connected. The following topics are auxiliary for the claimed purpose but of independent interest: (1) a characterization of linear complex Lie group (a~result similar to that of Luminet and Valette for real Lie groups); (2) properties of the exponential radical when $G$ is linear; (3) an asymptotic decomposition of a word length function into a sum of three summands (again for linear groups). The main result presents $\cO_{exp}(G)$ as a complete projective tensor of three factors, corresponding to the length function decomposition. As an application, it is shown that if $G$ is linear then the Arens-Michael envelope of $\cO_{exp}(G)$ is just the algebra of all holomorphic functions.


title: Kitaev quantum spin liquid - concept and materialization

abstract: A decade ago, Alexei Kitaev proposed an exactly solvable $S$ = 1/2 model on a two-dimensional honeycomb lattice, where the spins fractionalize into Majorana fermions and form a topological quantum spin liquid (QSL) in the ground state. It was soon recognized that a family of complex iridium oxides, as well as ruthenium chloride, with honeycomb structure are magnetic insulators and accommodate essential ingredients of the Kitaev model, due to the interplay of electron correlation and spin-orbit coupling. This initiated a race to materialize the Kitaev QSL and to capture the signature of Majorana fermions. In this review, we provide a wide perspective of this rapidly growing field, including theory, materials and experiment. We first summarize the theoretical background of the Kitaev QSL ground state and its materialization using spin-orbital-entangled $J_{\rm eff}$ = 1/2 moments. This is followed by an overview of candidate materials and their magnetic properties, including Na$_2$IrO$_3$, $\alpha$, $\beta$, $\gamma$-Li$_2$IrO$_3$,$\alpha$-RuCl$_3$ and H$_3$LiIr$_2$O$_6$. Finally, we review the latest exciting progress in the search for the Kitaev QSL. In particular, H$_3$LiIr$_2$O$_6$ and $\alpha$-RuCl$_3$ in applied magnetic field show signatures of the QSL state, and $\alpha$-RuCl$_3$ has unusual magnetic excitations and thermal transport properties that are consistent with spin fractionalization.


title: Learning Reciprocity in Complex Sequential Social Dilemmas

abstract: Reciprocity is an important feature of human social interaction and underpins our cooperative nature. What is more, simple forms of reciprocity have proved remarkably resilient in matrix game social dilemmas. Most famously, the tit-for-tat strategy performs very well in tournaments of Prisoner's Dilemma. Unfortunately this strategy is not readily applicable to the real world, in which options to cooperate or defect are temporally and spatially extended. Here, we present a general online reinforcement learning algorithm that displays reciprocal behavior towards its co-players. We show that it can induce pro-social outcomes for the wider group when learning alongside selfish agents, both in a $2$-player Markov game, and in $5$-player intertemporal social dilemmas. We analyse the resulting policies to show that the reciprocating agents are strongly influenced by their co-players' behavior.


title: A close look at the dwarf AGN of NGC 4395: optical and near-IR integral field spectroscopy

abstract: Intermediate mass black holes (10$^3$-10$^5$ M$_\odot$) in the center of dwarf galaxies are believed to be analogous to growing Active Galactic Nuclei (AGN) in the early Universe. Their characterization can provide insight about the early galaxies. We present optical and near-infrared integral field spectroscopy of the inner $\sim$50 pc of the dwarf galaxy NGC4395, known to harbor an AGN. NGC 4395 is an ideal candidate to investigate the nature of dwarf AGN, as it is nearby ($d\approx4.4$ Mpc) enough to allow a close look at its nucleus. The optical data were obtained with the Gemini GMOS-IFU covering the 4500 A to 7300 A spectral range at a spatial resolution of 10 pc. The J and K-band spectra were obtained with the Gemini NIFS at spatial resolutions of $\sim$5 pc. The gas kinematics show a compact, rotation disk component with a projected velocity amplitude of 25 km s$^{-1}$. We estimate a mass of $7.7\times10^5$ M$_\odot$ inside a radius of 10 pc. From the H$\alpha$ broad line component, we estimate the AGN bolometric luminosity as $L_{ bol}=(9.9\pm1.4)\times10^{40}$ erg s$^{-1}$ and a mass $M_{ BH}=(2.5^{+1.0}_{-0.8})\times10^5$ M$_\odot$ for the central black hole. The mean surface mass densities for the ionized and molecular gas are in the ranges (1-2) M$_{\odot} $pc$^{-2}$ and (1-4)$\times10^{-3}$ M${_\odot}$ pc$^{-2}$ and the average ratio between ionized and hot molecular gas masses is $\sim$500. The emission-line flux distributions reveal an elongated structure at 24 pc west of the nucleus, which is blueshifted relative to the systemic velocity of the galaxy by $\approx$30 km s$^{-1}$. We speculate that this structure is originated by the accretion of a gas-rich small satellite or by a low metallicity cosmic cloud.


title: Detailed studies of $^{100}$Mo two-neutrino double beta decay in NEMO-3

abstract: The full data set of the NEMO-3 experiment has been used to measure the half-life of the two-neutrino double beta decay of $^{100}$Mo to the ground state of $^{100}$Ru, $T_{1/2} = \left[ 6.81 \pm 0.01\,\left(\mbox{stat}\right) ^{+0.38}_{-0.40}\,\left(\mbox{syst}\right) \right] \times10^{18}$ y. The two-electron energy sum, single electron energy spectra and distribution of the angle between the electrons are presented with an unprecedented statistics of $5\times10^5$ events and a signal-to-background ratio of ~80. Clear evidence for the Single State Dominance model is found for this nuclear transition. Limits on Majoron emitting neutrinoless double beta decay modes with spectral indices of n=2,3,7, as well as constraints on Lorentz invariance violation and on the bosonic neutrino contribution to the two-neutrino double beta decay mode are obtained.


title: Experimental determination of line energies, line widths and relative transition probabilities of the Gadolinium L X-ray emission spectrum

abstract: In this work the complete L-emission spectrum of gadolinium with respect to line energies, natural line widths, and relative transition probabilities was investigated using monochromatized synchrotron radiation. The measurements were realized in the PTB laboratory at BESSY II by means of an in-house built von Hamos spectrometer based on up to two full-cylinder HAPG mosaic crystal. The von Hamos spectrometer is calibrated by means of elastically scattered photons from the employed synchrotron radiation beamline leading to a well-defined energy scale and an experimentally determined spectrometer response. A selective excitation of the gadolinium L subshells was carried out to ensure a robust deconvolution of neighboring emission lines of different L subshells. The experimental results are discussed in the context of existing data from common databases and published values since significant deviations, especially for the L$\gamma_2$ and L$\gamma_3$ emission lines, are observed. We further substantiate and discuss two satellite lines at the low-energy side of the L$\beta_{2,15}$ and L$\gamma_1$ emission lines arising from the N$_{4,5}$ subshell.


title: Regularity and Expansion for Steady Prandtl Equations

abstract: Due to degeneracy near the boundary, the question of high regularity for solutions to the steady Prandtl equations has been a longstanding open question since the celebrated work of Olenick. We settle this open question in affirmative in the absence of an external pressure. Our method is based on energy estimates for the quotient, $q = \frac{v}{\bar{u}}$, $\bar{u}$ being the classical Prandtl solution, via the linear Derivative Prandtl Equation (LDP). As a consequence, our regularity result leads to the construction of Prandtl layer expansion up to any order.


title: Proton CT -- a novel diagnostic tool in cancer therapy

abstract: Radiotherapy is one of the main methods in the successful treatment of cancer. The tumor is irradiated with photons or charged particles (e.g. protons), and in the case of massive charged particles, the treatment results in less unnecessary dose outside the tumor and therefore less side effects for the patient and a faster recovery. However, the dose planning of hadron therapy is calculated from photon CT measurements, which results in large uncertainties in the planning and therefore in a necessary enlargement of the treatment area. This uncertainty can be reduced by performing the CT scan using protons. The current contribution shows the development of a sampling calorimeter for proton CT measurements and describes the state of the project.


title: The effective field theory of dark energy diagnostic of linear Horndeski theories after GW170817 and GRB170817A

abstract: We summarise the effective field theory of dark energy construction to explore observable predictions of linear Horndeski theories. Based on \cite{Perenon:2016blf}, we review the diagnostic of these theories on the correlation of the large-scale structure phenomenological functions: the effective Newton constant, the light deflection parameter and the growth function of matter perturbations. We take this opportunity to discuss the evolution of the bounds the propagation speed of gravitational waves has undergone and use the most restrictive one to update the diagnostic.


title: Exponential mixing under controllability conditions for SDEs driven by a degenerate Poisson noise

abstract: We prove existence and uniqueness of the invariant measure and exponential mixing in the total-variation norm for a class of stochastic differential equations driven by degenerate compound Poisson processes. In addition to mild assumptions on the distribution of the jumps for the driving process, the hypotheses for our main result are that the corresponding control system is dissipative, approximately controllable and solidly controllable. The solid controllability assumption is weaker than the well-known parabolic H\"ormander condition and is only required from a single point to which the system is approximately controllable. Our analysis applies to Galerkin projections of stochastically forced parabolic partial differential equations with asymptotically polynomial nonlinearities and to networks of quasi-harmonic oscillators connected to different Poissonian baths.


title: Searching for Exosatellites Orbiting L and T Dwarfs: Connecting Planet Formation to Moon Formation and Finding New Temperate Worlds

abstract: L-type and T-type dwarfs span the boundaries between main-sequence stars, brown dwarfs, and planetary-mass objects. For these reasons, L and T dwarfs are the perfect laboratories for exploring the relationship between planet formation and moon formation, and evidence suggests they may be swarming with close-in rocky satellites, though none have been found to date. The discovery of satellites orbiting L or T dwarfs will have transformative implications for the nature of planets, moons and even life in the Universe. These transiting satellites will be prime targets for characterization with NASA's James Webb Space Telescope. In this white paper, we discuss the scientific motivations behind searching for transiting satellites orbiting L and T dwarfs and argue that robotizing current 1-to-2-meter US optical/infrared (O/IR) facilities and equipping them with recently developed low-cost infrared imagers will enable these discoveries in the next decade. Furthermore, robotizing the 1-to-2-meter O/IR fleet is highly synergistic with rapid follow-up of transient and multi-messenger events.


title: Uncountable structures are not classifiable up to bi-embeddability

abstract: Answering some of the main questions from [MR13], we show that whenever $\kappa$ is a cardinal satisfying $\kappa^{< \kappa} = \kappa > \omega$, then the embeddability relation between $\kappa$-sized structures is strongly invariantly universal, and hence complete for ($\kappa$-)analytic quasi-orders. We also prove that in the above result we can further restrict our attention to various natural classes of structures, including (generalized) trees, graphs, or groups. This fully generalizes to the uncountable case the main results of [LR05,FMR11,Wil14,CMR17].


title: Extended Falicov-Kimball model: Exact solution for finite temperatures

abstract: The extended Falicov-Kimball model is analyzed exactly for finite temperatures ($T\geq0$) in the limit of large dimensions. Onsite and intersite density-density interactions $U$ and $V$ are included in the model. Using the dynamical mean field theory formalism on the Bethe lattice we find rigorously the temperature dependent density of states (DOS) at half-filling. At $T=0$ the system is ordered to form the checkerboard pattern and the DOS has the gap $\Delta(\varepsilon_F) > 0$ at the Fermi level, if only $U\neq 0$ or $V\neq 0$. If $U <0$ or $U > 2V$, two additional subbands develop inside the principal energy gap. They become wider with increasing $T$ and at a certain $U$- and $V$-dependent temperature $T_{MI}$ they join with each other at $\varepsilon_F$. Since above $T_{MI}$ the DOS is positive at $\varepsilon_F$, we interpret $T_{MI}$ as the transformation temperature from insulator to metal. Moreover, we show that if $V\lesssim 0.54$ then $T_{MI}=0$ at two quasi-quantum critical points $U_{cr}^{\pm}$ (one positive and the other negative), whereas for $V\gtrsim 0.54$ there is only one negative $U^-_{cr}$. Having calculated the temperature dependent DOS we study thermodynamic properties of the system starting from its free energy and then we construct the phase diagrams in the variables $T$ and $U$ for a few values of $V$. Our calculations give that inclusion of the intersite coupling $V$ causes the finite temperature phase diagrams to become asymmetric with respect to a change of sign of $U$. On these phase diagrams we detected stability regions of eight different kinds of ordered phases, where both charge-order and antiferromagnetism coexists (five of them are insulating and three are conducting) and three different nonordered phases (two of them are insulating and one is conducting). Moreover, both continuous and discontinuous transitions between various phases were found.

categories: cond-mat.str-el cond-mat.other cond-mat.stat-mech physics.comp-ph


title: Phase mixing of nonlinear Alfven waves

abstract: Aims: This paper presents 2.5D numerical experiments of Alfv\'en wave phase mixing and aims to assess the effects of nonlinearities on wave behaviour and dissipation. In addition, this paper aims to quantify how effective the model presented in this work is at providing energy to the coronal volume. Methods: The model is presented and explored through the use of several numerical experiments which were carried out using the Lare2D code. The experiments study footpoint driven Alfv\'en waves in the neighbourhood of a two-dimensional x-type null point with initially uniform density and plasma pressure. A continuous sinusoidal driver with a constant frequency is used. Each experiment uses different driver amplitudes to compare weakly nonlinear experiments with linear experiments. Results: We find that the wave trains phase-mix owing to variations in the length of each field line and variations in the field strength. The nonlinearities reduce the amount of energy entering the domain, as they reduce the effectiveness of the driver, but they have relatively little effect on the damping rate (for the range of amplitudes studied). The nonlinearities produce density structures which change the natural frequencies of the field lines and hence cause the resonant locations to move. The shifting of the resonant location causes the Poynting flux associated with the driver to decrease. Reducing the magnetic diffusivity increases the energy build-up on the resonant field lines, however, it has little effect on the total amount of energy entering the system. From an order of magnitude estimate, we show that the Poynting flux in our experiments is comparable to the energy requirements of the quiet Sun corona. However a (possibly unphysically) large amount of magnetic diffusion was used however and it remains unclear if the model is able to provide enough energy under actual coronal conditions.


title: Corners for Layout: End-to-End Layout Recovery from 360 Images

abstract: The problem of 3D layout recovery in indoor scenes has been a core research topic for over a decade. However, there are still several major challenges that remain unsolved. Among the most relevant ones, a major part of the state-of-the-art methods make implicit or explicit assumptions on the scenes -- e.g. box-shaped or Manhattan layouts. Also, current methods are computationally expensive and not suitable for real-time applications like robot navigation and AR/VR. In this work we present CFL (Corners for Layout), the first end-to-end model for 3D layout recovery on 360 images. Our experimental results show that we outperform the state of the art relaxing assumptions about the scene and at a lower cost. We also show that our model generalizes better to camera position variations than conventional approaches by using EquiConvs, a type of convolution applied directly on the sphere projection and hence invariant to the equirectangular distortions. CFL Webpage: https://cfernandezlab.github.io/CFL/


title: The existence of supersingular curves of genus 4 in arbitrary characteristic

abstract: We prove that there exists a supersingular nonsingular curve of genus $4$ in arbitrary characteristic $p$. For $p>3$ we shall prove that the desingularization of a certain fiber product over $\mathbb{P}^1$ of two supersingular elliptic curves is supersingular.


title: A distance to the Large Magellanic Cloud that is precise to one per cent

abstract: In the era of precision cosmology, it is essential to empirically determine the Hubble constant with an accuracy of one per cent or better. At present, the uncertainty on this constant is dominated by the uncertainty in the calibration of the Cepheid period - luminosity relationship (also known as Leavitt Law). The Large Magellanic Cloud has traditionally served as the best galaxy with which to calibrate Cepheid period-luminosity relations, and as a result has become the best anchor point for the cosmic distance scale. Eclipsing binary systems composed of late-type stars offer the most precise and accurate way to measure the distance to the Large Magellanic Cloud. Currently the limit of the precision attainable with this technique is about two per cent, and is set by the precision of the existing calibrations of the surface brightness - colour relation. Here we report the calibration of the surface brightness-colour relation with a precision of 0.8 per cent. We use this calibration to determine the geometrical distance to the Large Magellanic Cloud that is precise to 1 per cent based on 20 eclipsing binary systems. The final distane is 49.59 +/- 0.09 (statistical) +/- 0.54 (systematic) kiloparsecs.


title: Natural Language Generation at Scale: A Case Study for Open Domain Question Answering

abstract: Current approaches to Natural Language Generation (NLG) for dialog mainly focus on domain-specific, task-oriented applications (e.g. restaurant booking) using limited ontologies (up to 20 slot types), usually without considering the previous conversation context. Furthermore, these approaches require large amounts of data for each domain, and do not benefit from examples that may be available for other domains. This work explores the feasibility of applying statistical NLG to scenarios requiring larger ontologies, such as multi-domain dialog applications or open-domain question answering (QA) based on knowledge graphs. We model NLG through an Encoder-Decoder framework using a large dataset of interactions between real-world users and a conversational agent for open-domain QA. First, we investigate the impact of increasing the number of slot types on the generation quality and experiment with different partitions of the QA data with progressively larger ontologies (up to 369 slot types). Second, we perform multi-task learning experiments between open-domain QA and task-oriented dialog, and benchmark our model on a popular NLG dataset. Moreover, we experiment with using the conversational context as an additional input to improve response generation quality. Our experiments show the feasibility of learning statistical NLG models for open-domain QA with larger ontologies.


title: Triangle singularities in $\bar{B}^0\to \chi_{c1}K^-\pi^+$ relevant to $Z_1(4050)$ and $Z_2(4250)$

abstract: $Z_1(4050)$ and $Z_2(4250)$ observed in $\bar{B}^0\to\chi_{c1}K^-\pi^+$ by the Belle Collaboration are candidates of charged charmonium-like states that minimally includes two quarks and two antiquarks. While $Z_1(4050)$ and $Z_2(4250)$ have been interpreted as tetraquark states previously, we propose a completely different scenario based on a kinematical effect called the triangle singularity. We demonstrate that the triangle singularities cause in the $\chi_{c1}\pi^+$ invariant mass distribution resonance-like bumps that fit very well the Belle data. If these bumps are simulated by the $Z_1(4050)$ and $Z_2(4250)$ resonance excitations, the spin-parity of them are predicted to be $1^-$ for $Z_1(4050)$ and $1^+$ or $1^-$ for $Z_2(4250)$. The bump corresponding to $Z_1(4050)$ has a highly asymmetric shape, which the Belle data exactly indicate. We show that the asymmetric shape originates from an interplay between the triangle singularity and the opening of the $X(3872)\pi^+$ channel near the triangle-singularity energy. This characteristic lineshape could be used to discriminate different interpretations of $Z_1(4050)$. An interesting prediction from interpretting $Z_1(4050)$ and $Z_2(4250)$ as the triangle singularities is that similar bumps caused by the same mechanisms possibly appear also in $\bar{B}^0\to J/\psi K^-\pi^+$ data; the already observed $Z_c(4200)$ corresponds to $Z_2(4250)$ of $J^P=1^+$.


title: Exploring the Interstellar Medium Using an Asymmetric X-ray Dust Scattering Halo

abstract: SWIFT J1658.2-4242 is an X-ray transient discovered recently in the Galactic plane, with severe X-ray absorption corresponding to an equivalent hydrogen column density of $N_{\rm H,abs}\sim2\times10^{23}$ cm$^{-2}$. Using new Chandra and XMM-Newton data, we discover a strong X-ray dust scattering halo around it. The halo profile can be well fitted by the scattering from at least three separated dust layers. During the persistent emission phase of SWIFT J1658.2-4242, the best-fit dust scattering $N_{\rm H,sca}$ based on the COMP-AC-S dust grain model is consistent with $N_{\rm H,abs}$. The best-fit halo models show that 85-90 percent of the intervening gas and dust along the line of sight of SWIFT J1658.2-4242 are located in the foreground ISM in the Galactic disk. The dust scattering halo also shows significant azimuthal asymmetry, which appears consistent with the inhomogeneous distribution of foreground molecular clouds. By matching the different dust layers to the distribution of molecular clouds along the line of sight, we estimate the source distance to be $\sim$10 kpc, which is also consistent with the results given by several other independent methods of distance estimation. The dust scattering opacity and the existence of a halo can introduce a significant spectral bias, the level of which depends on the shape of the instrumental point spread function and the source extraction region. We create the Xspec dscor model to correct for this spectral bias for different X-ray instruments. Our study reenforces the importance of considering the spectral effects of dust scattering in other absorbed X-ray sources.


title: Residual Deep Convolutional Neural Network for EEG Signal Classification in Epilepsy

abstract: Epilepsy is the fourth most common neurological disorder, affecting about 1% of the population at all ages. As many as 60% of people with epilepsy experience focal seizures which originate in a certain brain area and are limited to part of one cerebral hemisphere. In focal epilepsy patients, a precise surgical removal of the seizure onset zone can lead to effective seizure control or even a seizure-free outcome. Thus, correct identification of the seizure onset zone is essential. For clinical evaluation purposes, electroencephalography (EEG) recordings are commonly used. However, their interpretation is usually done manually by physicians and is time-consuming and error-prone. In this work, we propose an automated epileptic signal classification method based on modern deep learning methods. In contrast to previous approaches, the network is trained directly on the EEG recordings, avoiding hand-crafted feature extraction and selection procedures. This exploits the ability of deep neural networks to detect and extract relevant features automatically, that may be too complex or subtle to be noticed by humans. The proposed network structure is based on a convolutional neural network with residual connections. We demonstrate that our network produces state-of-the-art performance on two benchmark data sets, a data set from Bonn University and the Bern-Barcelona data set. We conclude that modern deep learning approaches can reach state-of-the-art performance on epileptic EEG classification and automated seizure onset zone identification tasks when trained on raw EEG data. This suggests that such approaches have potential for improving clinical practice.


title: Distortion analysis of crystalline and locally quasicrystalline 2D photonic structures with GISAXS

abstract: In this study, grazing incidence small-angle X-ray scattering (GISAXS) is used to collect statistical information on dimensional parameters in an area of 20 mm x 15 mm on photonic structures produced by nanoimprint lithography. The photonic structures are composed of crystalline and locally quasicrystalline two-dimensional patterns with structure sizes between about 100 nm and 10 $\mu$m to enable broadband visible light absorption for use in solar energy harvesting. These first GISAXS measurements on locally quasicrystalline samples demonstrate that GISAXS is capable of showing the locally quasicrystalline nature of the samples while at the same time revealing the long-range periodicity introduced due to the lattice design. We describe the scattering qualitatively in the framework of the distorted wave Born approximation using a hierarchical model mirroring the sample design, which consists of a rectangular and locally quasicrystalline supercell which is repeated periodically to fill the whole surface. The nanoimprinted samples are compared to a sample manufactured using electron beam lithography and the distortions of the periodic and locally quasiperiodic samples are quantified statistically. Due to the high sensitivity of GISAXS to deviations from the perfect lattice, the misalignment of the crystallographic axes was measured with a resolution of 0.015{\deg}, showing distortions up to +/- 0.15{\deg} in the investigated samples.


title: Injecting Software Vulnerabilities with Voltage Glitching

abstract: We show how voltage glitching can cause timing violations in CMOS behavior. Then we attack a real, security hardened, consumer device to gain code execution and dump the secure boot ROM.


title: Particle-Set Identification method to study multiplicity fluctuations

abstract: In this paper a new method of experimental data analysis, the Particle-Set Identification method, is presented. The method allows to reconstruct moments of multiplicity distribution of identified particles. The difficulty the method copes with is due to incomplete particle identification -- a particle mass is frequently determined with a resolution which does not allow for a unique determination of the particle type. Within this method the moments of order $k$ are calculated from mean multiplicities of $k$-particle sets of a given type. The Particle-Set Identification method remains valid even in the case of correlations between mass measurements for different particles. This distinguishes it from the Identity method introduced by us previously to solve the problem of incomplete particle identification in studies of particle fluctuations.


title: Union of holomorphically convex spaces

abstract: In this short note, we collect some results regarding the Remmert reduction of holomorphically convex space and its application to a variation of the usual union problem. Classically, the union problem asks the following question: is a complex space, which is an increasing union of Stein subspaces $X_1\Subset X_2\Subset\cdots$, a Stein space itself? The variation we are interested in is the following: is a complex space, which is an increasing union of holomorphically convex subspaces $X_1\Subset X_2\Subset\cdots$, holomorphically convex itself? The results presented here are close analogues of (some of) those alredy present in the literature for the Stein case; our aim is only to collect such material for reference, as we consider it well known.


title: Principal component analysis of nonequilibrium molecular dynamics simulations

abstract: Principal component analysis (PCA) represents a standard approach to identify collective variables $\{x_i\}\!=\!\boldsymbol{x}$, which can be used to construct the free energy landscape $\Delta G(\boldsymbol{x})$ of a molecular system. While PCA is routinely applied to equilibrium molecular dynamics (MD) simulations, it is less obvious how to extend the approach to nonequilibrium simulation techniques. This includes, e.g., the definition of the statistical averages employed in PCA, as well as the relation between the equilibrium free energy landscape $\Delta G(\boldsymbol{x})$ and energy landscapes $\Delta{\cal G} (\boldsymbol{x})$ obtained from nonequilibrium MD. As an example for a nonequilibrium method, `targeted MD' is considered which employs a moving distance constraint to enforce rare transitions along some biasing coordinate $s$. The introduced bias can be described by a weighting function $P(s)$, which provides a direct relation between equilibrium and nonequilibrium data, and thus establishes a well-defined way to perform PCA on nonequilibrium data. While the resulting distribution ${\cal P}(\boldsymbol{x})$ and energy $\Delta{\cal G} \propto \ln {\cal P}$ will not reflect the equilibrium state of the system, the nonequilibrium energy landscape $\Delta{\cal G} (\boldsymbol{x})$ may directly reveal the molecular reaction mechanism. Applied to targeted MD simulations of the unfolding of decaalanine, for example, a PCA performed on backbone dihedral angles is shown to discriminate several unfolding pathways. Although the formulation is in principle exact, its practical use depends critically on the choice of the biasing coordinate $s$, which should account for a naturally occurring motion between two well-defined end-states of the system.


title: Study of magnetization reversal in Neel and Bloch regime of Ni and Py stripes using Kerr microscopy

abstract: We present a systematic study of the magnetization reversal of nickel and permalloy micro-stripes with Neel and Bloch domain walls using Kerr microscopy. Magnetic field driven domain propagation was observed from higher width to lower width stripes for magnetic fields applied along the length of micro-stripes. Stripe like domains were observed with nucleation starting in lower width region followed by their propagation to higher width regions for magnetic fields applied along the width of micro-stripes. The comparison of magnetization reversal in Bloch and N\'eel domain wall regime showed higher domain wall density in Bloch regime for both nickel and permalloy stripes.


title: Fibers of multi-graded rational maps and orthogonal projection onto rational surfaces

abstract: We contribute a new algebraic method for computing the orthogonal projections of a point onto a rational algebraic surface embedded in the three dimensional projective space. This problem is first turned into the computation of the finite fibers of a generically finite dominant rational map: a congruence of normal lines to the rational surface. Then, an in-depth study of certain syzygy modules associated to such a congruence is presented and applied to build elimination matrices that provide universal representations of its finite fibers, under some genericity assumptions. These matrices depend linearly in the variables of the three dimensional space. They can be pre-computed so that the orthogonal projections of points are approximately computed by means of fast and robust numerical linear algebra calculations.


title: A compact high order Alternating Direction Implicit method for three-dimensional acoustic wave equation with variable coefficient

abstract: Efficient and accurate numerical simulation of seismic wave propagation is important in various Geophysical applications such as seismic full waveform inversion (FWI) problem. However, due to the large size of the physical domain and requirement on low numerical dispersion, many existing numerical methods are inefficient for numerical modelling of seismic wave propagation in a heterogeneous media. Despite the great efforts that have been devoted during the past decades, it still remains a challenging task in the development of efficient and accurate finite difference method for the multi-dimensional acoustic wave equation with variable velocity. In this paper, we proposed a Pad\'{e} approximation based finite difference scheme for solving the acoustic wave equation in three-dimensional heterogeneous media. The new method is obtained by combining the Pad\'{e} approximation and a novel algebraic manipulation. The efficiency of the new algorithm is further improved through the Alternative Directional Implicit (ADI) method. The stability of the new algorithm has been theoretically proved by the energy method. The new method is conditionally stable with a better Courant - Friedrichs - Lewy condition (CFL) condition, which has been verified numerically. Extensive numerical examples have been solved, which demonstrated that the new method is accurate, efficient and stable.


title: Polynomial analogue of Gandhi's fixed point theorem

abstract: The problem to be solved in this paper is to construct a general method of proving whether a certain set is p-computable or not. The method is based on a polynomial analogue of the classical Gandhi's fixed point theorem. The classical Gandhi theorem uses the extension of the predicate with the help of the special operator $\Gamma^{\Omega^*}_{\Phi(x)}$ whose smallest fixed point is the $\Sigma$-set. The work uses a new type of operator - $\Delta_0^p$-operator $\Gamma_{F_{P_1^{+}},...,F_{P_n^{+}}}^{\mathfrak{M}}$, which extends predicates so that the smallest fixed point remains a p-computable set. Moreover, if in the classical Gandhi's fixed point theorem the special $\Sigma$-formula $\Phi(\overline {x})$ is used in the construction of the operator, then in the new operator, instead of a single formula, special generating families of formulas $F_ {P_1 ^ {+}},...,F_{P_n^{+}}$. This work opens up broad prospects for the application of the polynomial analogue of the Gandhi theorem in the construction of new types of terms and formulas, in the construction of new data types and programs of polynomial computational complexity in Turing complete languages.


title: Online Non-Convex Learning: Following the Perturbed Leader is Optimal

abstract: We study the problem of online learning with non-convex losses, where the learner has access to an offline optimization oracle. We show that the classical Follow the Perturbed Leader (FTPL) algorithm achieves optimal regret rate of $O(T^{-1/2})$ in this setting. This improves upon the previous best-known regret rate of $O(T^{-1/3})$ for FTPL. We further show that an optimistic variant of FTPL achieves better regret bounds when the sequence of losses encountered by the learner is `predictable'.


title: Preconditioned P-ULA for Joint Deconvolution-Segmentation of Ultrasound Images -- Extended Version

abstract: Joint deconvolution and segmentation of ultrasound images is a challenging problem in medical imaging. By adopting a hierarchical Bayesian model, we propose an accelerated Markov chain Monte Carlo scheme where the tissue reflectivity function is sampled thanks to a recently introduced proximal unadjusted Langevin algorithm. This new approach is combined with a forward-backward step and a preconditioning strategy to accelerate the convergence, and with a method based on the majorization-minimization principle to solve the inner nonconvex minimization problems. As demonstrated in numerical experiments conducted on both simulated and in vivo ultrasound images, the proposed method provides high-quality restoration and segmentation results and is up to six times faster than an existing Hamiltonian Monte Carlo method.


title: Ab initio Investigation of Structural Stability and Exfoliation Energies in Transition Metal Dichalcogenides based on Ti-, V-, and Mo-Group Elements

abstract: In this work, we report an ab initio investigation based on density functional theory of the structural, energetic and electronic properties of 2D layered chalcogenides compounds based in the combination of the transition-metals (Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W and chalcogenides (S, Se, Te) in three polymorphic phases: trigonal prismatic (2H), octahedral (1T) and distorted octahedral (1T$_{\text{d}}$). We determined the most stable phases for each compound, verifying the existence of the 1T$_{\text{d}}$ phase for a small number of the compounds and we have also identified the magnetic compounds. In addition, with the determination of the exfoliation energies, we indicated the potential candidates to form one layer material and we have also found a relation between the exfoliation energy and the effective Bader charge in the metal, suggesting that when the materials present small exfoliation energy, it is due to the Coulomb repulsion between the chalcogen planes. Finally, we analyzed the electronic properties, identifying the semiconductor, semimetal and metal materials and predicting the band gap of the semiconductors. In our results, the dependence of the band gap on the $d$-orbital is explicit. In conclusion, we have investigated the properties of stable and metastable phases for a large set of TMD materials, and our findings may be auxiliary in the synthesis of metastable phases and in the development of new TMDs applications.


title: Identifying Experts in Software Libraries and Frameworks among GitHub Users

abstract: Software development increasingly depends on libraries and frameworks to increase productivity and reduce time-to-market. Despite this fact, we still lack techniques to assess developers expertise in widely popular libraries and frameworks. In this paper, we evaluate the performance of unsupervised (based on clustering) and supervised machine learning classifiers (Random Forest and SVM) to identify experts in three popular JavaScript libraries: facebook/react, mongodb/node-mongodb, and socketio/socket.io. First, we collect 13 features about developers activity on GitHub projects, including commits on source code files that depend on these libraries. We also build a ground truth including the expertise of 575 developers on the studied libraries, as self-reported by them in a survey. Based on our findings, we document the challenges of using machine learning classifiers to predict expertise in software libraries, using features extracted from GitHub. Then, we propose a method to identify library experts based on clustering feature data from GitHub; by triangulating the results of this method with information available on Linkedin profiles, we show that it is able to recommend dozens of GitHub users with evidences of being experts in the studied JavaScript libraries. We also provide a public dataset with the expertise of 575 developers on the studied libraries.


title: Exact Gaussian Processes on a Million Data Points

abstract: Gaussian processes (GPs) are flexible non-parametric models, with a capacity that grows with the available data. However, computational constraints with standard inference procedures have limited exact GPs to problems with fewer than about ten thousand training points, necessitating approximations for larger datasets. In this paper, we develop a scalable approach for exact GPs that leverages multi-GPU parallelization and methods like linear conjugate gradients, accessing the kernel matrix only through matrix multiplication. By partitioning and distributing kernel matrix multiplies, we demonstrate that an exact GP can be trained on over a million points, a task previously thought to be impossible with current computing hardware, in less than 2 hours. Moreover, our approach is generally applicable, without constraints to grid data or specific kernel classes. Enabled by this scalability, we perform the first-ever comparison of exact GPs against scalable GP approximations on datasets with $10^4 \!-\! 10^6$ data points, showing dramatic performance improvements.


title: Integral Fluctuation Relations for Entropy Production at Stopping Times

abstract: A stopping time $T$ is the first time when a trajectory of a stochastic process satisfies a specific criterion. In this paper, we use martingale theory to derive the integral fluctuation relation $\langle e^{-S_{\rm tot}(T)}\rangle=1$ for the stochastic entropy production $S_{\rm tot}$ in a stationary physical system at stochastic stopping times $T$. This fluctuation relation implies the law $\langle S_{\rm tot}(T)\rangle\geq 0$, which states that it is not possible to reduce entropy on average, even by stopping a stochastic process at a stopping time, and which we call the second law of thermodynamics at stopping times. This law implies bounds on the average amount of heat and work a system can extract from its environment when stopped at a random time. Furthermore, the integral fluctuation relation implies that certain fluctuations of entropy production are universal or are bounded by universal functions. These universal properties descend from the integral fluctuation relation by selecting appropriate stopping times: for example, when $T$ is a first-passage time for entropy production, then we obtain a bound on the statistics of negative records of entropy production. We illustrate these results on simple models of nonequilibrium systems described by Langevin equations and reveal two interesting phenomena. First, we demonstrate that isothermal mesoscopic systems can extract on average heat from their environment when stopped at a cleverly chosen moment and the second law at stopping times provides a bound on the average extracted heat. Second, we demonstrate that the average efficiency at stopping times of an autonomous stochastic heat engines, such as Feymann's ratchet, can be larger than the Carnot efficiency and the second law of thermodynamics at stopping times provides a bound on the average efficiency at stopping times.


title: Low Mass Stars as Tracers of Star Formation in Diverse Environments

abstract: Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.


title: Black hole in a superconducting plasma

abstract: The generalized vortical formalism provides an electrodynamic description for superconducting states---in the generalized vortical formalism, a superconducting state may be defined by the vanishing of an appropriate generalized vorticity and characterized by zero generalized helicity for incompressible fluids. In this article, we investigate these states for compressible plasmas in black hole spacetime geometries using the curved spacetime generalization of the grand generalized vortical formalism. If the magnetic field is axisymmetric and the thermodynamic properties are symmetric about the equatorial plane, the resulting states are characterized by a vanishing skin depth and a complete expulsion of the magnetic field at the equator of the black hole horizon. Moreover, if the thermodynamic properties of the plasma are uniform at the horizon, we find that the magnetic field is completely expelled from the horizon, and the plasma behaves as a perfect superconductor near the horizon. This result is independent of the spin of a black hole, holding even for a (nonrotating) Schwarzschild black hole, and demonstrates that the geometry near black hole horizons can have a significant effect on the electrodynamics of surrounding plasmas.


title: Recovery of non-smooth coefficients appearing in anisotropic wave equations

abstract: We study the problem of unique recovery of a non-smooth one-form $\mathcal A$ and a scalar function $q$ from the Dirichlet to Neumann map, $\Lambda_{\mathcal A,q}$, of a hyperbolic equation on a Riemannian manifold $(M,g)$. We prove uniqueness of the one-form $\mathcal A$ up to the natural gauge, under weak regularity conditions on $\mathcal A,q$ and under the assumption that $(M,g)$ is simple. Under an additional regularity assumption, we also derive uniqueness of the scalar function $q$. The proof is based on the geometric optic construction and inversion of the light ray transform extended as a Fourier Integral Operator to non-smooth parameters and functions.


title: Gauss-Bonnet black holes with a massive scalar field

abstract: In the present paper we consider the extended scalar-tensor-Gauss-Bonnet gravity with a massive scalar field. We prove numerically the existence of Gauss-Bonnet black holes for three different forms of the coupling function including the case of spontaneous scalarization. We have performed a systematic study of the black hole characteristics such as the area of the horizon, the entropy and the temperature for these coupling functions and compared them to the Schwarzschild solutions. The introduction of scalar field mass leads to a suppression of the scalar field and the increase of this mass brings the black holes closer to the Schwarzschild case. For linear and exponential coupling, a nonzero scalar field mass expands the domain of existence of black holes solutions. Larger deviations from the GR case are observed only for small masses and these differences decrease with the increase of the scalar field mass. In the case of a coupling function which leads to scalarization the scalar field mass has a significant influence on the bifurcation points where the scalarized black holes branch out of the Schwarzschild solution. The largest deviation from the case with a massless scalar field are observed for black hole masses close to the bifurcation point.


title: Actuator Placement for Optimizing Network Performance under Controllability Constraints

abstract: With the rising importance of large-scale network control, the problem of actuator placement has received increasing attention. Our goal in this paper is to find a set of actuators minimizing the metric that measures the average energy consumption of the control inputs while ensuring structural controllability of the network. As this problem is intractable, greedy algorithm can be used to obtain an approximate solution. To provide a performance guarantee for this approach, we first define the submodularity ratio for the metric under consideration and then reformulate the structural controllability constraint as a matroid constraint. This shows that the problem under study can be characterized by a matroid optimization involving a weakly submodular objective function. Then, we derive a novel performance guarantee for the greedy algorithm applied to this class of optimization problems. Finally, we show that the matroid feasibility check for the greedy algorithm can be cast as a maximum matching problem in a certain auxiliary bipartite graph related to the network graph.


title: Achievement of FCC specification in critical current density for Nb3Sn superconductors with artificial pinning centers

abstract: In this letter we demonstrate achievement of record non-Cu critical current density (Jc,non-Cu) in ternary, multifilamentary Nb3Sn conductors by the introduction of artificial pinning centers (APC). In the past two years, we have made great progress in the development of APC Nb3Sn wires. Recent resistivity vs magnetic field measurements confirmed the high upper critical field (Bc2) of ternary APC wires, which at 4.2 K was ~28 T, about 1-2 T higher than present state-of-the-art conductors. In addition to high Bc2, it was found that APC wires have noticeably higher Sn content in the Nb3Sn layers as compared to standard wires. The Jc,non-Cu values of the most-recent APC wires have met the Jc,non-Cu-B specification required by the Future Circular Collider (FCC), with the best heat treatment leading to a Jc,nonCu 29% higher than the FCC specification at 21 T. Microscopy analysis shows that the APC wires still have overly high residual Nb fractions due to too low of a Sn/Nb ratio, indicating that there is still great potential for further Jc,non-Cu improvement. The development of APC wires is ongoing; this letter details some of the steps forward in the optimization and lays out a roadmap to push the APC wires towards practical, magnet-grade conductors.


title: Gemini Imaging of the Host Galaxies of Changing-Look Quasars

abstract: Changing-look quasars are a newly-discovered class of luminous active galactic nuclei that undergo rapid ($\lesssim$10 year) transitions between Type 1 and Type 1.9/2, with an associated change in their continuum emission. We characterize the host galaxies of four faded changing-look quasars using broadband optical imaging. We use \textit{gri} images obtained with the Gemini Multi Object Spectrograph (GMOS) on Gemini North to characterize the surface brightness profiles of the quasar hosts and search for [O III] $\lambda4959,\lambda5007$ emission from spatially extended regions, or voorwerpjes, with the goal of using them to examine past luminosity history. Although we do not detect, voorwerpjes surrounding the four quasar host galaxies, we take advantage of the dim nuclear emission to characterize the colors and morphologies of the host galaxies. Three of the four galaxies show morphological evidence of merger activity or tidal features in their residuals. The three galaxies which are not highly distorted are fit with a single S\'ersic profile to characterize their overall surface brightness profiles. The single-S\'ersic fits give intermediate S\'ersic indices between the $n=1$ of disk galaxies and the $n=4$ of ellipticals. On a color-magnitude diagram, our changing-look quasar host galaxies reside in the blue cloud, with other AGN host galaxies and star-forming galaxies. On a color-S\'ersic index diagram the changing-look quasar hosts reside with other AGN hosts in the "green valley". Our analysis suggests that the hosts of changing-look quasars are predominantly disrupted or merging galaxies that resemble AGN hosts, rather than inactive galaxies.


title: Residual finiteness and strict distortion of cyclic subgroups of solvable groups

abstract: We provide polynomial lower bounds for residual finiteness of residually finite, finitely generated solvable groups that admit infinite order elements in the Fitting subgroup of strict distortion at least exponential. For this class of solvable groups which include polycyclic groups with a nontrivial exponential radical and the metabelian Baumslag-Solitar groups, we improve the lower bounds found in the literature. Additionally, for the class of residually finite, finitely generated solvable groups of infinite Pr\"{u}fer rank that satisfy the conditions of our theorem, we provide the first nontrivial lower bounds.


title: Regularity for higher order quasiconvex problems with linear growth from below

abstract: We announce new existence and $\varepsilon$-regularity results for minimisers of the relaxation of strongly quasiconvex integrals that on smooth maps $u\colon\Omega\subset\mathbb{R}^{n}\to\mathbb{R}^{N}$ are defined by $$u\mapsto \int_{\Omega}F(\nabla^{k}u)dx.$$ The results cover the case of integrands $F$ with $(1,q)$-growth in the full range of exponents $1<q<\frac{n}{n-1}$ for which a measure representation of the relaxed functional is possible and the minimizers belong to the space $BV^k$ of maps whose $k$-th order derivatives are measures.


title: Predictive clustering

abstract: We show how to convert any clustering into a prediction set. This has the effect of converting the clustering into a (possibly overlapping) union of spheres or ellipsoids. The tuning parameters can be chosen to minimize the size of the prediction set. When applied to k-means clustering, this method solves several problems: the method tells us how to choose k, how to merge clusters and how to replace the Voronoi partition with more natural shapes. We show that the same reasoning can be applied to other clustering methods.


title: Observation of $e^{+}e^{-}\rightarrow \pi^{+}\pi^{-}\psi(3770)$ and $D_{1}(2420)^{0}\bar{D}^{0} + c.c.$

abstract: Several intermediate states of the reaction channels $e^{+}e^{-} \rightarrow \pi^{+} \pi^{-} D^{0} \bar{D}^{0}$ and $e^{+}e^{-} \rightarrow \pi^{+} \pi^{-} D^{+}D^{-}$ are studied using the data samples collected with the BESIII detector at center-of-mass energies above 4.08 GeV. For the first time in this final state, a $\psi(3770)$ signal is seen in the $D\bar{D}$ invariant mass spectrum, with a statistical significance of $5.2\sigma$ at $\sqrt{s} = 4.42$ GeV. There is also evidence for this resonance at $\sqrt{s}$ = 4.26 and 4.36 GeV with statistical significance of 3.2$\sigma$ and 3.3$\sigma$, respectively. In addition, the Born cross section of $e^{+}e^{-}\to \pi^{+}\pi^{-}\psi(3770)$ is measured. The proposed heavy-quark-spin-symmetry partner of the $X(3872)$, the state $X_{2}(4013)$, is also searched for in the $D\bar{D}$ invariant mass spectra. No obvious signal is found. The upper limit of the Born cross section of the process $e^{+}e^{-}\to \rho^{0}X_{2}(4013)$ combined with the branching fraction is measured. Also, the processes $e^{+}e^{-}\to D_{1}(2420)\bar{D} + c.c.$ are investigated. The neutral mode with $D_{1}(2420)^{0}\to D^{0}\pi^{+}\pi^{-}$ is reported with statistical significance of 7.4$\sigma$ at $\sqrt{s} = 4.42$ GeV for the first time, and evidence with statistical significance of 3.2$\sigma$ and 3.3$\sigma$ at $\sqrt{s}$ = 4.36 and 4.60 GeV is seen, respectively. No evident signal for the process $e^{+}e^{-}\to D_{1}(2420)^{0}\bar{D}^{0} + c.c., D_{1}(2420)^{0}\to D^{*+}\pi^{-}$ is reported. Evidence for $e^{+}e^{-}\to D_{1}(2420)^{+}D^{-} +~c.c., D_{1}(2420)^{+}\to D^{+}\pi^{+}\pi^{-}$ is reported with statistical significance of 3.1$\sigma$ and 3.0$\sigma$ at $\sqrt{s}$ = 4.36 and 4.42 GeV, respectively.


title: Asymptotic Max-Min SINR Analysis of Reconfigurable Intelligent Surface Assisted MISO Systems

abstract: This work focuses on the downlink of a single-cell multi-user system in which a base station (BS) equipped with $M$ antennas communicates with $K$ single-antenna users through a reconfigurable intelligent surface (RIS) installed in the line-of-sight (LoS) of the BS. RIS is envisioned to offer unprecedented spectral efficiency gains by utilizing $N$ passive reflecting elements that induce phase shifts on the impinging electromagnetic waves to smartly reconfigure the signal propagation environment. We study the minimum signal-to-interference-plus-noise ratio (SINR) achieved by the optimal linear precoder (OLP), that maximizes the minimum SINR subject to a given power constraint for any given RIS phase matrix, for the cases where the LoS channel matrix between the BS and the RIS is of rank-one and of full-rank. In the former scenario, the minimum SINR achieved by the RIS-assisted link is bounded by a quantity that goes to zero with $K$. For the high-rank scenario, we develop accurate deterministic approximations for the parameters of the asymptotically OLP, which are then utilized to optimize the RIS phase matrix. Simulation results show that RISs can outperform half-duplex relays with a small number of passive reflecting elements while large RISs are needed to outperform full-duplex relays.


title: The Dynamic Infrared Sky

abstract: Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade.


title: Hyper-Parameter Sweep on AlphaZero General

abstract: Since AlphaGo and AlphaGo Zero have achieved breakground successes in the game of Go, the programs have been generalized to solve other tasks. Subsequently, AlphaZero was developed to play Go, Chess and Shogi. In the literature, the algorithms are explained well. However, AlphaZero contains many parameters, and for neither AlphaGo, AlphaGo Zero nor AlphaZero, there is sufficient discussion about how to set parameter values in these algorithms. Therefore, in this paper, we choose 12 parameters in AlphaZero and evaluate how these parameters contribute to training. We focus on three objectives~(training loss, time cost and playing strength). For each parameter, we train 3 models using 3 different values~(minimum value, default value, maximum value). We use the game of play 6$\times$6 Othello, on the AlphaZeroGeneral open source re-implementation of AlphaZero. Overall, experimental results show that different values can lead to different training results, proving the importance of such a parameter sweep. We categorize these 12 parameters into time-sensitive parameters and time-friendly parameters. Moreover, through multi-objective analysis, this paper provides an insightful basis for further hyper-parameter optimization.


title: Spin-polarized Correlated Insulator and Superconductor in Twisted Double Bilayer Graphene

abstract: Ferromagnetism and superconductivity typically compete with each other since the internal magnetic field generated in a magnet suppresses the formation of spin-singlet Cooper pairs in conventional superconductors. Only a handful of ferromagnetic superconductors are known in heavy fermion systems, where many-body electron interactions promoted by the narrow energy bands play a key role in stabilizing these emergent states. Recently, interaction-driven superconductivity and ferromagnetism have been demonstrated as separate phenomena in different density regimes of flat bands enabled by graphene moire superlattices. Combining superconductivity and magnetism in a single ground state may lead to more exotic quantum phases. Here, employing van der Waals heterostructures of twisted double bilayer graphene (TDBG), we realize a flat electron band that is tunable by perpendicular electric fields. Similar to the magic angle twisted bilayer graphene, TDBG exhibits energy gaps at the half and quarter filled flat bands, indicating the emergence of correlated insulating states. We find that the gaps of these insulating states increase with in-plane magnetic field, suggesting a ferromagnetic order. Upon doping the ferromagnetic half-filled insulator, superconductivity emerges with a critical temperature controlled by both density and electric fields. We observe that the in-plane magnetic field enhances the superconductivity in the low field regime, suggesting spin-polarized electron pairing. Spin-polarized superconducting states discovered in TDBG provide a new route to engineering interaction-driven topological superconductivity.


title: Kernel-based Translations of Convolutional Networks

abstract: Convolutional Neural Networks, as most artificial neural networks, are commonly viewed as methods different in essence from kernel-based methods. We provide a systematic translation of Convolutional Neural Networks (ConvNets) into their kernel-based counterparts, Convolutional Kernel Networks (CKNs), and demonstrate that this perception is unfounded both formally and empirically. We show that, given a Convolutional Neural Network, we can design a corresponding Convolutional Kernel Network, easily trainable using a new stochastic gradient algorithm based on an accurate gradient computation, that performs on par with its Convolutional Neural Network counterpart. We present experimental results supporting our claims on landmark ConvNet architectures comparing each ConvNet to its CKN counterpart over several parameter settings.


title: ExplainIt! -- A declarative root-cause analysis engine for time series data (extended version)

abstract: We present ExplainIt!, a declarative, unsupervised root-cause analysis engine that uses time series monitoring data from large complex systems such as data centres. ExplainIt! empowers operators to succinctly specify a large number of causal hypotheses to search for causes of interesting events. ExplainIt! then ranks these hypotheses, reducing the number of causal dependencies from hundreds of thousands to a handful for human understanding. We show how a declarative language, such as SQL, can be effective in declaratively enumerating hypotheses that probe the structure of an unknown probabilistic graphical causal model of the underlying system. Our thesis is that databases are in a unique position to enable users to rapidly explore the possible causal mechanisms in data collected from diverse sources. We empirically demonstrate how ExplainIt! had helped us resolve over 30 performance issues in a commercial product since late 2014, of which we discuss a few cases in detail.


title: Layer-specific connectivity revealed by diffusion-weighted functional MRI in the rat thalamocortical pathway

abstract: Investigating neural activity from a global brain perspective in-vivo has been in the domain of functional Magnetic Resonance Imaging (fMRI) over the past few decades. The intricate neurovascular couplings that govern fMRI's blood-oxygenation-level-dependent (BOLD) functional contrast are invaluable in mapping active brain regions, but they also entail significant limitations, such as non-specificity of the signal to active foci. Diffusion-weighted functional MRI (dfMRI) with relatively high diffusion-weighting strives to ameliorate this shortcoming as it offers functional contrasts more intimately linked with the underlying activity. Insofar, apart from somewhat smaller activation foci, dfMRI's contrasts have not been convincingly shown to offer significant advantages over BOLD, and its contrasts relied on significant modelling. Here, we study whether dfMRI could offer a better representation of neural activity in the thalamocortical pathway compared to its (spin-echo (SE)) BOLD counterpart. Using high-end forepaw stimulation experiments in the rat at 9.4 T, and with significant sensitivity enhancements due to the use of cryocoils, we show for the first time that dfMRI signals exhibit layer specificity, and, additionally, display signals in areas devoid of SE-BOLD responses. We find that dfMRI signals in the thalamocortical pathway cohere with each other, namely, dfMRI signals in the ventral posterolateral (VPL) thalamic nucleus cohere specifically with layers IV and V in the somatosensory cortex. These activity patterns are much better correlated (compared with SE-BOLD signals) with literature-based electrophysiological recordings in the cortex as well as thalamus. All these findings suggest that dfMRI signals better represent the underlying neural activity in the pathway. In turn, this may entail significant implications towards a much more specific and accurate (...)


title: Astro2020 Science White Paper: A Shocking Shift in Paradigm for Classical Novae

abstract: The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories.


title: A Local Limit Theorem and Delocalization of Eigenvectors for Polynomials in Two Matrices

abstract: We propose a boundary regularity condition for the $M_n(\mathbb{C})$-valued subordination functions in free probability to prove the local limit theorem and delocalization of eigenvectors for polynomials in two random matrices. We prove this through estimating the pair of $M_n(\mathbb{C})$-valued approximate subordination functions for the sum of two $M_n(\mathbb{C})$-valued random matrices $\gamma_1\otimes C_N+\gamma_2\otimes U_N^*D_NU_N$, where $C_N$, $D_N$ are deterministic diagonal matrices, and $U_N$ is Haar unitary.


title: Debiasing Community Detection: The Importance of Lowly-Connected Nodes

abstract: Community detection is an important task in social network analysis, allowing us to identify and understand the communities within the social structures. However, many community detection approaches either fail to assign low degree (or lowly-connected) users to communities, or assign them to trivially small communities that prevent them from being included in analysis. In this work, we investigate how excluding these users can bias analysis results. We then introduce an approach that is more inclusive for lowly-connected users by incorporating them into larger groups. Experiments show that our approach outperforms the existing state-of-the-art in terms of F1 and Jaccard similarity scores while reducing the bias towards low-degree users.


title: Probing the trilinear Higgs boson coupling in di-Higgs production at NLO QCD including parton shower effects

abstract: We present results for Higgs boson pair production with variations of the trilinear Higgs boson self-coupling at next-to-leading order (NLO) in QCD including the full top quark mass dependence. Differential results at 14 TeV are presented, and we discuss the implications of anomalous trilinear couplings as well as differences between the PYTHIA 8.2 and HERWIG 7.1 parton showers in combination with POWHEG. The implementation of the NLO QCD calculation with variable Higgs boson self-coupling is made publicly available in the POWHEG-BOX-V2 Monte Carlo framework. A simple method for using the new implementation to study also variations of the top quark Yukawa coupling is described.


title: On the Sprague-Grundy function of compound games

abstract: The classical game of {\sc Nim} can be naturally extended and played on an arbitrary hypergraph $\cH \subseteq 2^V \setminus \{\emptyset\}$ whose vertices $V = \{1, \ldots, n\}$ correspond to piles of stones. By one move a player chooses an edge $H$ of $\cH$ and reduces arbitrarily all piles $i \in H$. In 1901 Bouton solved the classical {\sc Nim} for which $\cH = \{\{1\}, \ldots, \{n\}\}$. In 1910 Moore introduced and solved a more general game $k$-{\sc Nim}, for which $\cH = \{H \subseteq V \mid |H| \leq k\}$, where $1 \leq k < n$. In 1980 Jenkyns and Mayberry obtained an explicit formula for the Sprague-Grundy function of Moore's {\sc Nim} for the case $k+1 = n$. Recently it was shown that the same formula works for a large class of hypergraphs. In this paper we study combinatorial properties of these hypergraphs and obtain explicit formulas for the Sprague-Grundy functions of the conjunctive and selective compounds of the corresponding hypergraph {\sc Nim} games.


title: Regular domains and surfaces of constant Gaussian curvature in three-dimensional affine space

abstract: Generalizing the notion of domains of dependence in the Minkowski space, we define and study regular domains in the affine space with respect to a proper convex cone. In dimension three, we show that every proper regular domain is uniquely foliated by a particular kind of surfaces with constant affine Gaussian curvature. The result is based on the analysis of a Monge-Amp\`ere equation with extended-real-valued lower semicontinuous boundary condition.


title: Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry

abstract: PURPOSE: Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures. METHODS: PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared. RESULTS: Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of 2 to 4x, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution. CONCLUSION: Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.


title: Butterfly in a Cocoon, Understanding the origin and morphology of Globular Cluster Streams: The case of GD-1

abstract: Tidally disrupted globular cluster streams are usually observed, and therefore perceived, as narrow, linear and one-dimensional structures in the 6D phase-space. Here we show that the GD-1 stellar stream ($\approx$ 30 pc wide), which is the tidal debris of a disrupted globular cluster, possesses a secondary diffuse and extended stellar component ($\approx$ 100 pc wide) around it, detected at >5$\sigma$ confidence level. Similar morphological properties are seen in synthetic streams that are produced from star clusters that are formed within dark matter sub-halos and then accrete onto a massive host galaxy. This lends credence to the idea that the progenitor of the highly retrograde GD-1 stream was originally formed outside of the Milky Way in a now defunct dark satellite galaxy. We deem that in future studies, this newly found $cocoon$ component may serve as a structural hallmark to distinguish between the in-situ and ex-situ (accreted) formed globular cluster streams.


title: Device-to-Device Coded Caching with Distinct Cache Sizes

abstract: This paper considers a cache-aided device-to-device (D2D) system where the users are equipped with cache memories of different size. During low traffic hours, a server places content in the users' cache memories, knowing that the files requested by the users during peak traffic hours will have to be delivered by D2D transmissions only. The worst-case D2D delivery load is minimized by jointly designing the uncoded cache placement and linear coded D2D delivery. Next, a novel lower bound on the D2D delivery load with uncoded placement is proposed and used in explicitly characterizing the minimum D2D delivery load (MD2DDL) with uncoded placement for several cases of interest. In particular, having characterized the MD2DDL for equal cache sizes, it is shown that the same delivery load can be achieved in the network with users of unequal cache sizes, provided that the smallest cache size is greater than a certain threshold. The MD2DDL is also characterized in the small cache size regime, the large cache size regime, and the three-user case. Comparisons of the server-based delivery load with the D2D delivery load are provided. Finally, connections and mathematical parallels between cache-aided D2D systems and coded distributed computing (CDC) systems are discussed.


title: Stimulated binding of polymer chains by narrow tube confinement

abstract: In biology, there are several processes in which unfolded protein chains are transported along narrow-tube channels. Normally, without such a severe configurational constraint, unfolded polypeptides would not bind to each other. However, when chain entropy is much reduced in the narrow channel, we find that polypeptide chains have a propensity to bind, even if there is no great potential energy gain in doing so. We find the average length of binding m* (the number of monomers at the chain ends that form bonds) and the critical tube diameter at which such constrained binding occurs. We carry out Brownian dynamics simulations of tightly confined chains, demonstrating their binding over the characteristic length m*, changing in tubes of different diameter.


title: Where are the Intermediate Mass Black Holes?

abstract: Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well as wandering in more massive galaxy halos. However, these environments are not conducive to the accretion events or dynamical signatures that allow us to detect IMBHs. The Laser Interferometer Space Antenna (LISA) will demystify IMBHs by detecting the mergers of these objects out to extremely high redshifts, while measuring their masses with extremely high precision. These observations of merging IMBHs will allow us to constrain the formation mechanism and subsequent evolution of massive black holes, from the 'dark ages' to the present day, and reveal the role that IMBHs play in hierarchical galaxy evolution.


title: BiHom-Novikov algebras and infinitesimal BiHom-bialgebras

abstract: We introduce and study infinitesimal BiHom-bialgebras, BiHom-Novikov algebras, BiHom-Novikov-Poisson algebras, and find some relations among these concepts. Our main result is to show how to obtain a left BiHom-pre-Lie algebra from an infinitesimal BiHom-bialgebra


title: Factorization method versus migration imaging in a waveguide

abstract: We present a comparative study of two qualitative imaging methods in an acoustic waveguide with sound hard walls. The waveguide terminates at one end and contains unknown obstacles of compact support, to be determined from data gathered by an array of sensors that probe the obstacles with waves and measure the scattered response. The first imaging method, known as the factorization method, is based on the factorization of the far field operator. It is designed to image at single frequency and estimates the support of the obstacles by a Picard range criterion. The second imaging method, known as migration, works either with one or multiple frequencies. It forms an image by backpropagating the measured scattered wave to the search points, using the Green's function in the empty waveguide. We study the connection between these methods with analysis and numerical simulations.


title: Classification of $(1{,}2)$-reflective anisotropic hyperbolic lattices of rank $4$

abstract: A hyperbolic lattice is called \textit{$(1{,}2)$-reflective} if its automorphism group is generated by $1$- and $2$-reflections up to finite index. In this paper we prove that the fundamental polyhedron of a $\mathbb{Q}$-arithmetic cocompact reflection group in the three-dimensional Lobachevsky space contains an edge such that the distance between its framing faces is small enough. Using this fact we obtain a classification of $(1{,}2)$-reflective anisotropic hyperbolic lattices of rank $4$.


title: Lifetimes of Doubly Heavy Baryons ${\cal B}_{bb}$ and ${\cal B}_{bc}$

abstract: Lifetimes of the doubly heavy baryons ${\cal B}_{bb}$ and ${\cal B}_{bc}$ are analyzed within the framework of the heavy quark expansion (HQE). Lifetime differences arise from the spectator effects such as $W$-exchange and Pauli interference. For doubly bottom baryons, the lifetime pattern is $\tau(\Omega_{bb}^-)\sim \tau(\Xi_{bb}^{-})>\tau(\Xi_{bb}^0)$. The $\Xi_{bb}^{0}$ baryon is shortest-lived owing to the $W$-exchange contribution, while $\Xi_{bb}^{-}$ and $\Omega_{bb}^{-}$ have similar lifetimes as they both receive contributions from destructive Pauli interference. We find the lifetime ratio $\tau(\Xi_{bb}^{-})/\tau(\Xi_{bb}^0)=1.26$\,. The large $W$-exchange contribution to $\Xi_{bc}^0$ through the subprocess $cd\to us\to cd$ and the sizable destructive Pauli interference contribution to $\Xi_{bc}^+$ imply a substantial lifetime difference between $\Xi_{bc}^+$ and $\Xi_{bc}^0$. In the presence of subleading $1/m_c$ and $1/m_b$ corrections to the spectator effects, we find that $\tau(\Omega_{bc}^0)$ becomes longest-lived. This is because $\Gamma^{\rm int}_+$ and $\Gamma^{\rm semi}$ for $\Omega_{bc}^0$ are subject to large cancellation between dimension-6 and -7 operators. This implies that the subleading corrections are too large to justify the validity of the HQE. Demanding that $\Gamma^{cs}_{{\rm int+}}(\Omega_{bc}^0)$, $\Gamma^{{\rm SL},cs}_{\rm int}(\Omega_{bc}^0)$ be positive and $\Gamma^{cu}_{{\rm int-}}(\Xi^+_{bc})$ be negative, we conjecture that $1.68\times 10^{-13}s<\tau(\Omega_{bc}^0)< 3.70\times 10^{-13}s$ , $4.09\times 10^{-13}s<\tau(\Xi_{bc}^+)< 6.07\times 10^{-13}s$ and $0.93\times 10^{-13}s<\tau(\Xi_{bc}^0)< 1.18\times 10^{-13}s$. Hence, the lifetime hierarchy of ${\cal B}_{bc}$ baryons is expected to be the pattern $\tau(\Xi_{bc}^{+})>\tau(\Omega_{bc}^0)>\tau(\Xi_{bc}^0)$.