title: Gravitoviscous protoplanetary disks with a dust component. I. The importance of the inner sub-au region

abstract: The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of prestellar cores of 0.5--1.0~$M_\odot$ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central "smart" cell (CSC) -- a simplified model that simulates physical processes that may occur in this region. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks and they are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates reaching the limit when streaming instability becomes operational. We argue that gravitational instability, together with streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits (Abridged).


title: Quantum correlations from dynamically modulated optical nonlinear interactions

abstract: We investigate optical nonlinear interactions in a dynamic environment by studying generation of photons in spontaneous parametric down conversion inside a nonlinear cavity where the optical path length is periodically modulated in time. We show that the temporal dynamics of the cavity modify the nonlinear interaction and the generated continuous variable time-frequency entangled bi-photon state evolves into a tunable discrete higher dimensional state in the non-adiabatic modulation regime where the modulation time scales are much faster than the photon lifetime. In this regime, the system mimics effects of a quantum random walk in a photonic lattice with many associated effects including localized and delocalized wavefunctions of the generated photons. We also propose generation of time-frequency hyper-entangled states in the adiabatic limit. Our analysis shows that the proposed system is promising for applications in quantum simulation and information processing in the time-frequency domain.


title: Entanglement subvolume law for 2D frustration-free spin systems

abstract: Let $H$ be a frustration-free Hamiltonian describing a 2D grid of qudits with local interactions, a unique ground state, and local spectral gap lower bounded by a positive constant. For any bipartition defined by a vertical cut of length $L$ running from top to bottom of the grid, we prove that the corresponding entanglement entropy of the ground state of $H$ is upper bounded by $\tilde{O}(L^{5/3})$. For the special case of a 1D chain, our result provides a new area law which improves upon prior work, in terms of the scaling with qudit dimension and spectral gap. In addition, for any bipartition of the grid into a rectangular region $A$ and its complement, we show that the entanglement entropy is upper bounded as $\tilde{O}(|\partial A|^{5/3})$ where $\partial A$ is the boundary of $A$. This represents the first subvolume bound on entanglement in frustration-free 2D systems. In contrast with previous work, our bounds depend on the local (rather than global) spectral gap of the Hamiltonian. We prove our results using a known method which bounds the entanglement entropy of the ground state in terms of certain properties of an approximate ground state projector (AGSP). To this end, we construct a new AGSP which is based on a robust polynomial approximation of the AND function and we show that it achieves an improved trade-off between approximation error and entanglement.


title: Estimates for the SVD of the truncated Fourier transform on L2(exp(b|$\times$|)) and stable analytic continuation

abstract: The Fourier transform truncated on [-c,c] is usually analyzed when acting on L^2(-1/b,1/b) and its right-singular vectors are the prolate spheroidal wave functions. This paper considers the operator acting on the larger space L^2(exp(b|.|)) on which it remains injective. We give nonasymptotic upper and lower bounds on the singular values with similar qualitative behavior in m (the index), b, and c. The lower bounds are used to obtain rates of convergence for stable analytic continuation of possibly nonbandlimited functions whose Fourier transform belongs to L^2(exp(b|.|)). We also derive bounds on the sup-norm of the singular functions. Finally, we propose a numerical method to compute the SVD and apply it to stable analytic continuation when the function is observed with error on an interval.


title: Magnetic Fields and Afterglows of BdHNe: Inferences from GRB 130427A, GRB 160509A, GRB 160625B, GRB 180728A and GRB 190114C

abstract: GRB 190114C is the first binary-driven hypernova (BdHN) fully observed from the initial supernova appearance to the final emergence of the optical SN signal. It offers an unprecedented testing ground for the BdHN theory and it is here determined and further extended to additional gamma-ray bursts (GRBs). BdHNe comprise two subclasses of long GRBs with progenitors a binary system composed of a carbon-oxygen star (CO$_\textrm{core}$) and a neutron star (NS) companion. The CO$_\textrm{core}$ explodes as a SN leaving at its center a newborn NS ($\nu$NS). The SN ejecta hypercritically accretes both on the $\nu$NS and the NS companion. BdHNe I are the tightest binaries where the accretion leads the companion NS to gravitational collapse into a black hole (BH). In BdHN II the accretion onto the NS is lower, so there is no BH formation. We observe the same structure of the afterglow for GRB 190114C and other selected examples of BdHNe I (GRB 130427A, GRB 160509A, GRB 160625B) and for BdHN II (GRB 180728A). In all the cases the explanation of the afterglow is reached via the synchrotron emission powered by the $\nu$NS: their magnetic fields structures and their spin are determined. For BdHNe I, we discuss the properties of the magnetic field embedding the newborn BH, inherited from the collapsed NS and amplified during the gravitational collapse process, and surrounded by the SN ejecta.


title: Parallel and Communication Avoiding Least Angle Regression

abstract: We are interested in parallelizing the Least Angle Regression (LARS) algorithm for fitting linear regression models to high-dimensional data. We consider two parallel and communication avoiding versions of the basic LARS algorithm. The two algorithms have different asymptotic costs and practical performance. One offers more speedup and the other produces more accurate output. The first is bLARS, a block version of LARS algorithm, where we update b columns at each iteration. Assuming that the data are row-partitioned, bLARS reduces the number of arithmetic operations, latency, and bandwidth by a factor of b. The second is Tournament-bLARS (T-bLARS), a tournament version of LARS where processors compete by running several LARS computations in parallel to choose b new columns to be added in the solution. Assuming that the data are column-partitioned, T-bLARS reduces latency by a factor of b. Similarly to LARS, our proposed methods generate a sequence of linear models. We present extensive numerical experiments that illustrate speedups up to 4x compared to LARS without any compromise in solution quality.


title: Dynamical stability analysis of tumor growth and invasion: A reaction-diffusion model

abstract: The acid-mediated tumor invasion hypothesis proposes that altered glucose metabolism exhibited by the vast majority of tumors leads to increased acid (H+ ion) production which subsequently facilitates tumor invasion [1-3]. The reaction-diffusion model [2] that captures the key elements of the hypothesis shows how the densities of normal cells, tumor cells, and excess H+ ions change with time due to both chemical reactions between these three populations and density-dependent diffusion by which they spread out in three-dimensional space. Moreover, it proposes that each cell has an optimal pH for survival; that is, if the local pH deviates from the optimal value in either an acidic or alkaline direction, the cells begin to die, and that the death rate saturates at some maximum value when the microenvironment is extremely acidic or alkaline. We have previously studied in detail how the death-rate functions of the normal and tumor populations depend upon the H+ ion density [4]. Here, we extend previous work by investigating how the equilibrium densities (at which the time rates of change of the cellular densities are equal to zero) reached by the normal and tumor populations in three-dimensional space are affected by the presence of the H+ ions, and we present detailed analytical and computational techniques to analyze the dynamical stability of these equilibrium densities. For a sample set of biological input parameters and within the acid-mediation hypothesis, our model predicts the transformation to a malignant behavior, as indicated by the presence of unstable sets of equilibrium densities.


title: Social Cards Probably Provide For Better Understanding Of Web Archive Collections

abstract: Used by a variety of researchers, web archive collections have become invaluable sources of evidence. If a researcher is presented with a web archive collection that they did not create, how do they know what is inside so that they can use it for their own research? Search engine results and social media links are represented as surrogates, small easily digestible summaries of the underlying page. Search engines and social media have a different focus, and hence produce different surrogates than web archives. Search engine surrogates help a user answer the question "Will this link meet my information need?" Social media surrogates help a user decide "Should I click on this?" Our use case is subtly different. We hypothesize that groups of surrogates together are useful for summarizing a collection. We want to help users answer the question of "What does the underlying collection contain?" But which surrogate should we use? With Mechanical Turk participants, we evaluate six different surrogate types against each other. We find that the type of surrogate does not influence the time to complete the task we presented the participants. Of particular interest are social cards, surrogates typically found on social media, and browser thumbnails, screen captures of web pages rendered in a browser. At $p=0.0569$, and $p=0.0770$, respectively, we find that social cards and social cards paired side-by-side with browser thumbnails probably provide better collection understanding than the surrogates currently used by the popular Archive-It web archiving platform. We measure user interactions with each surrogate and find that users interact with social cards less than other types. The results of this study have implications for our web archive summarization work, live web curation platforms, social media, and more.


title: Markov jump processes and collision-like models in the kinetic description of multi-agent systems

abstract: Multi-agent systems can be successfully described by kinetic models, which allow one to explore the large scale aggregate trends resulting from elementary microscopic interactions. The latter may be formalised as collision-like rules, in the spirit of the classical kinetic approach in gas dynamics, but also as Markov jump processes, which assume that every agent is stimulated by the other agents to change state according to a certain transition probability distribution. In this paper we establish a parallelism between these two descriptions, whereby we show how the understanding of the kinetic jump process models may be improved taking advantage of techniques typical of the collisional approach.


title: Modulation-resonance mechanism for surface waves in a two-layer fluid system

abstract: We propose a Boussinesq-type model to study the surface/interfacial wave manifestation of an underlying, slowly-varying, long-wavelength, baroclinic flow in a two-layer, density-stratified system. The results of our model show numerically that, under strong nonlinearity, surface waves, with their typical wavenumber being the resonant $k_{\mathrm{res}}$, can be generated locally at the leading edge of the underlying slowly-varying, long-wavelength baroclinic flow. Here, the resonant $k_{\mathrm{res}}$ satisfies the class 3 triad resonance condition among two short-mode waves and one long-mode wave in which all waves propagate in the same direction. Moreover, when the slope of the baroclinic flow is sufficiently small, only one spatially-localized large-amplitude surface wave packet can be generated at the leading edge. This localized surface wave packet becomes high in amplitude and large in group velocity after the interaction with its surrounding waves. These results are qualitatively consistent with various experimental observations including resonant surface waves at the leading edge of an internal wave. Subsequently, we propose a mechanism, referred to as the modulation-resonance mechanism, underlying these surface phenomena, based on our numerical simulations. The proposed modulation-resonance mechanism combines the linear modulation (ray-based) theory for the spatiotemporal asymmetric behavior of surface waves and the nonlinear class 3 triad resonance theory for the energy focusing of surface waves around the resonant wavenumber $k_{\mathrm{res}}$ in Fourier space.


title: Expression of Fractals Through Neural Network Functions

abstract: To help understand the underlying mechanisms of neural networks (NNs), several groups have, in recent years, studied the number of linear regions $\ell$ of piecewise linear functions generated by deep neural networks (DNN). In particular, they showed that $\ell$ can grow exponentially with the number of network parameters $p$, a property often used to explain the advantages of DNNs over shallow NNs in approximating complicated functions. Nonetheless, a simple dimension argument shows that DNNs cannot generate all piecewise linear functions with $\ell$ linear regions as soon as $\ell > p$. It is thus natural to seek to characterize specific families of functions with $\ell$ linear regions that can be constructed by DNNs. Iterated Function Systems (IFS) generate sequences of piecewise linear functions $F_k$ with a number of linear regions exponential in $k$. We show that, under mild assumptions, $F_k$ can be generated by a NN using only $\mathcal{O}(k)$ parameters. IFS are used extensively to generate, at low computational cost, natural-looking landscape textures in artificial images. They have also been proposed for compression of natural images, albeit with less commercial success. The surprisingly good performance of this fractal-based compression suggests that our visual system may lock in, to some extent, on self-similarities in images. The combination of this phenomenon with the capacity, demonstrated here, of DNNs to efficiently approximate IFS may contribute to the success of DNNs, particularly striking for image processing tasks, as well as suggest new algorithms for representing self similarities in images based on the DNN mechanism.


title: Error Analysis and Correction for Weighted A*'s Suboptimality (Extended Version)

abstract: Weighted A* (wA*) is a widely used algorithm for rapidly, but suboptimally, solving planning and search problems. The cost of the solution it produces is guaranteed to be at most W times the optimal solution cost, where W is the weight wA* uses in prioritizing open nodes. W is therefore a suboptimality bound for the solution produced by wA*. There is broad consensus that this bound is not very accurate, that the actual suboptimality of wA*'s solution is often much less than W times optimal. However, there is very little published evidence supporting that view, and no existing explanation of why W is a poor bound. This paper fills in these gaps in the literature. We begin with a large-scale experiment demonstrating that, across a wide variety of domains and heuristics for those domains, W is indeed very often far from the true suboptimality of wA*'s solution. We then analytically identify the potential sources of error. Finally, we present a practical method for correcting for two of these sources of error and experimentally show that the correction frequently eliminates much of the error.


title: Spatial realization of a Lie algebra and Bar construction of a group

abstract: We prove that the spatial realization of a rational complete Lie algebra $L$, concentrated in degree 0, is isomorphic to the simplicial bar construction on the group, obtained from the Baker-Campbell-Hausdorff product on $L$.


title: Ultracold collisions of Cs in excited hyperfine and Zeeman states

abstract: We investigate Cs+Cs scattering in excited Zeeman and hyperfine states. We calculate the real and imaginary parts of the s-wave scattering length; the imaginary part directly provides the rate coefficient for 2-body inelastic loss, while the real part allows us to identify regions of magnetic field where 3-body recombination will be slow. We identify field regions where Cs in its $(f,m_f)=(3,+2)$ and $(3,+1)$ states may be stable enough to allow Bose-Einstein condensation, and additional regions for these and the $(3,0)$ and $(3,-3)$ states where high-density clouds should be long-lived.


title: Full-Stack, Real-System Quantum Computer Studies: Architectural Comparisons and Design Insights

abstract: In recent years, Quantum Computing (QC) has progressed to the point where small working prototypes are available for use. Termed Noisy Intermediate-Scale Quantum (NISQ) computers, these prototypes are too small for large benchmarks or even for Quantum Error Correction, but they do have sufficient resources to run small benchmarks, particularly if compiled with optimizations to make use of scarce qubits and limited operation counts and coherence times. QC has not yet, however, settled on a particular preferred device implementation technology, and indeed different NISQ prototypes implement qubits with very different physical approaches and therefore widely-varying device and machine characteristics. Our work performs a full-stack, benchmark-driven hardware-software analysis of QC systems. We evaluate QC architectural possibilities, software-visible gates, and software optimizations to tackle fundamental design questions about gate set choices, communication topology, the factors affecting benchmark performance and compiler optimizations. In order to answer key cross-technology and cross-platform design questions, our work has built the first top-to-bottom toolflow to target different qubit device technologies, including superconducting and trapped ion qubits which are the current QC front-runners. We use our toolflow, TriQ, to conduct {\em real-system} measurements on 7 running QC prototypes from 3 different groups, IBM, Rigetti, and University of Maryland. From these real-system experiences at QC's hardware-software interface, we make observations about native and software-visible gates for different QC technologies, communication topologies, and the value of noise-aware compilation even on lower-noise platforms. This is the largest cross-platform real-system QC study performed thus far; its results have the potential to inform both QC device and compiler design going forward.


title: On Mappings on the Hypercube with Small Average Stretch

abstract: Let $A \subseteq \{0,1\}^n$ be a set of size $2^{n-1}$, and let $\phi \colon \{0,1\}^{n-1} \to A$ be a bijection. We define the average stretch of $\phi$ as ${\sf avgStretch}(\phi) = {\mathbb E}[{\sf dist}(\phi(x),\phi(x'))]$, where the expectation is taken over uniformly random $x,x' \in \{0,1\}^{n-1}$ that differ in exactly one coordinate. In this paper we continue the line of research studying mappings on the discrete hypercube with small average stretch. We prove the following results. (1) For any set $A \subseteq \{0,1\}^n$ of density $1/2$ there exists a bijection $\phi_A \colon \{0,1\}^{n-1} \to A$ such that ${\sf avgstretch}(\phi_A) = O(\sqrt{n})$. (2) For $n = 3^k$ let $A_{{\sf rec\text{-}maj}} = \{x \in \{0,1\}^n : {\sf rec\text{-}maj}(x) = 1\}$, where ${\sf rec\text{-}maj} : \{0,1\}^n \to \{0,1\}$ is the function recursive majority of 3's. There exists a bijection $\phi_{{\sf rec\text{-}maj}} \colon \{0,1\}^{n-1} \to A_{\sf rec\text{-}maj}$ such that ${\sf avgstretch}(\phi_{\sf rec\text{-}maj}) = O(1)$. (3) Let $A_{\sf tribes} = \{x \in \{0,1\}^n : {\sf tribes}(x) = 1\}$. There exists a bijection $\phi_{{\sf tribes}} \colon \{0,1\}^{n-1} \to A_{\sf tribes}$ such that ${\sf avgstretch}(\phi_{{\sf tribes}}) = O(\log(n))$. These results answer the questions raised by Benjamini et al.\ (FOCS 2014).


title: On the descriptive power of Neural-Networks as constrained Tensor Networks with exponentially large bond dimension

abstract: In many cases, Neural networks can be mapped into tensor networks with an exponentially large bond dimension. Here, we compare different sub-classes of neural network states, with their mapped tensor network counterpart for studying the ground state of short-range Hamiltonians. We show that when mapping a neural network, the resulting tensor network is highly constrained and thus the neural network states do in general not deliver the naive expected drastic improvement against the state-of-the-art tensor network methods. We explicitly show this result in two paradigmatic examples, the 1D ferromagnetic Ising model and the 2D antiferromagnetic Heisenberg model, addressing the lack of a detailed comparison of the expressiveness of these increasingly popular, variational ans\"atze.

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


title: Detection and timing of gamma-ray pulsations from the $707$ Hz pulsar J0952$-$0607

abstract: The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.


title: CoRide: Joint Order Dispatching and Fleet Management for Multi-Scale Ride-Hailing Platforms

abstract: How to optimally dispatch orders to vehicles and how to tradeoff between immediate and future returns are fundamental questions for a typical ride-hailing platform. We model ride-hailing as a large-scale parallel ranking problem and study the joint decision-making task of order dispatching and fleet management in online ride-hailing platforms. This task brings unique challenges in the following four aspects. First, to facilitate a huge number of vehicles to act and learn efficiently and robustly, we treat each region cell as an agent and build a multi-agent reinforcement learning framework. Second, to coordinate the agents from different regions to achieve long-term benefits, we leverage the geographical hierarchy of the region grids to perform hierarchical reinforcement learning. Third, to deal with the heterogeneous and variant action space for joint order dispatching and fleet management, we design the action as the ranking weight vector to rank and select the specific order or the fleet management destination in a unified formulation. Fourth, to achieve the multi-scale ride-hailing platform, we conduct the decision-making process in a hierarchical way where a multi-head attention mechanism is utilized to incorporate the impacts of neighbor agents and capture the key agent in each scale. The whole novel framework is named as CoRide. Extensive experiments based on multiple cities real-world data as well as analytic synthetic data demonstrate that CoRide provides superior performance in terms of platform revenue and user experience in the task of city-wide hybrid order dispatching and fleet management over strong baselines.


title: Stacking, Strain, & Twist in 2D Materials Quantified by 3D Electron Diffraction

abstract: The field of two-dimensional (2D) materials has expanded to multilayered systems where electronic, optical, and mechanical properties change-often dramatically-with stacking order, thickness, twist, and interlayer spacing [1-5]. For transition metal dichalcogenides (TMDs), bond coordination within a single van der Waals layer changes the out-of-plane symmetry that can cause metal-insulator transitions [1, 6] or emergent quantum behavior [7]. Discerning these structural order parameters is often difficult using real-space measurements, however, we show 2D materials have distinct, conspicuous three-dimensional (3D) structure in reciprocal space described by near infinite oscillating Bragg rods. Combining electron diffraction and specimen tilt we probe Bragg rods in all three dimensions to identify multilayer structure with sub-Angstrom precision across several 2D materials-including TMDs (MoS2, TaSe2, TaS2) and multilayer graphene. We demonstrate quantitative determination of key structural parameters such as surface roughness, inter- & intra-layer spacings, stacking order, and interlayer twist using a rudimentary transmission electron microscope (TEM). We accurately characterize the full interlayer stacking order of multilayer graphene (1-, 2-, 6-, 12-layers) as well the intralayer structure of MoS2 and extract a chalcogen-chalcogen layer spacing of 3.07 +/- 0.11 Angstrom. Furthermore, we demonstrate quick identification of multilayer rhombohedral graphene.


title: Binary Neutron Star (BNS) merger: What we learned from relativistic ejecta of GW/GRB~170817A

abstract: Gravitational waves from coalescence of a Binary Neutron Star (BNS) and its accompagning short Gamma-Ray Burst GW/GRB~170817A confirmed the presumed origin of these puzzeling transients and opened up the way for relating properties of short GRBs to those of their progenitor stars and their surroundings. Here we review an extensive analysis of the prompt gamma-ray and late afterglows of this event. We show that a fraction of polar ejecta from the merger had been accelerated to ultra-relativistic speeds. This structured jet had an initial Lorentz factor of about $260$ in our direction - $\mathcal{O}(10^\circ)$ from the jet's axis - and was a few orders of magnitude less dense than in typical short GRBs. At the time of arrival to circum-burst material the ultra-relativistic jet had a close to Gaussian profile and a Lorentz factor $\gtrsim 130$ in its core. It had retained in some extent its internal collimation and coherence, but had extended laterally to create mildly relativistic lobes - a {\it cocoon}. External shocks on the far from center inhomogeneous circum-burst material and low density of colliding shells generated slow rising afterglows. The circum-burst material was somehow correlated with the merger and it is possible that it contained recently ejected material from glitching, which had resumed due to the deformation of neutron stars crust by tidal forces in the latest stages of inspiral but well before their merger. By comparing these findings with the results of relativistic MHD simulations and observed gravitational waves we conclude that progenitor neutron stars were old, had close masses and highly reduced magnetic fields. In addition, they probably had oppositely directed spins due to the encounter and gravitational interaction with other stars.


title: The molecular gas content of shell galaxies

abstract: Shells are fine stellar structures identified by their arc-like shapes present around a galaxy and currently thought to be vestiges of galaxy interactions and/or mergers. The study of their number, geometry, stellar populations and gas content can help to derive the interaction/merger history of a galaxy. Numerical simulations have proposed a mechanism of shell formation through phase wrapping during a radial minor merger. Alternatively, there could be barely a space wrapping, when particles have not made any radial oscillation yet, but are bound by their radial expansion, or produce an edge-brightened feature. These can be distinguished, because they are expected to keep a high radial velocity. While shells are first a stellar phenomenon, HI and CO observations have revealed neutral gas associated with shells. Some of the gas, the most diffuse and dissipative, is expected to be driven quickly to the center if it is travelling on nearly radial orbits. Molecular gas, distributed in dense clumps, is less dissipative, and may be associated to shells, and determine their velocity, too difficult to obtain from stars. We present here a search for molecular gas in nine shell galaxies with the IRAM-30m telescope. Six of them are detected in their galaxy center, and in three galaxies, we clearly detect molecular gas in shells. The derived amount of molecular gas varies from 1.5 10$^8$ to 3.4 10$^9$ M$_\odot$ in the shells. For two of them (Arp 10 and NGC 3656), the shells are characteristic of an oblate system. Their velocity is nearly systemic, and we conclude that these shells are phase-wrapped. For the third one (NGCB3934) the shells appear to participate to the rotation, and follow up with higher spatial resolution is required to conclude.


title: Cubic interaction vertices for N=1 arbitrary spin massless supermultiplets in flat space

abstract: In the framework of light-cone gauge formulation, massless arbitrary spin N=1 supermultiplets in four-dimensional flat space are considered. We study both the integer (super)spin and half-integer (super)spin supermultiplets. For such supermultiplets, formulation in terms of unconstrained light-cone gauge superfields defined on chiral momentum superspace is used. Superfield representation for all cubic interaction vertices of the supermultiplets is obtained. Representation of the cubic vertices in terms of component fields is derived. Realization of relativistic symmetries of N=1 Poincare superalgebra on space of interacting superfields is also found.


title: Straight to Shapes++: Real-time Instance Segmentation Made More Accurate

abstract: Instance segmentation is an important problem in computer vision, with applications in autonomous driving, drone navigation and robotic manipulation. However, most existing methods are not real-time, complicating their deployment in time-sensitive contexts. In this work, we extend an existing approach to real-time instance segmentation, called `Straight to Shapes' (STS), which makes use of low-dimensional shape embedding spaces to directly regress to object shape masks. The STS model can run at 35 FPS on a high-end desktop, but its accuracy is significantly worse than that of offline state-of-the-art methods. We leverage recent advances in the design and training of deep instance segmentation models to improve the performance accuracy of the STS model whilst keeping its real-time capabilities intact. In particular, we find that parameter sharing, more aggressive data augmentation and the use of structured loss for shape mask prediction all provide a useful boost to the network performance. Our proposed approach, `Straight to Shapes++', achieves a remarkable 19.7 point improvement in mAP (at IOU of 0.5) over the original method as evaluated on the PASCAL VOC dataset, thus redefining the accuracy frontier at real-time speeds. Since the accuracy of instance segmentation is closely tied to that of object bounding box prediction, we also study the error profile of the latter and examine the failure modes of our method for future improvements.


title: Large-Scale Statistical Survey of Magnetopause Reconnection

abstract: The Magnetospheric Multiscale Mission (MMS) seeks to study the micro-physics of reconnection, which occurs at the magnetopause boundary layer between the magnetosphere of Earth and the interplanetary magnetic field originating from the sun. Identifying this region of space automatically will allow for statistical analysis of reconnection events. The magnetopause region is difficult to identify automatically using simple models, and time consuming for scientists to classify by hand. We introduced a hierarchical Bayesian mixture model with linear and auto regressive components to identify the magnetopause. Using data from the MMS mission with the programming languages R and Stan, we modeled and predicted possible regions and evaluated our performance against a boosted regression tree model. Our model selects twice as many magnetopause regions as the comparison model, without significant over selection, achieving a 31\% true positive rate and 93\% true negative rate. Our method will allow scientists to study the micro-physics of reconnection events in the magnetopause using the large body of MMS data without manual classification.


title: Brick: Asynchronous Payment Channels

abstract: Off-chain protocols (channels) are a promising solution to the scalability and privacy challenges of blockchain payments. Current proposals, however, require synchrony assumptions to preserve the safety of a channel, leaking to an adversary the exact amount of time needed to control the network for a successful attack. In this paper, we introduce Brick, the first payment channel that remains secure under network asynchrony and concurrently provides correct incentives. The core idea is to incorporate the conflict resolution process within the channel by introducing a rational committee of external parties, called Wardens. Hence, if a party wants to close a channel unilaterally, it can only get the committee's approval for the last valid state. Brick provides sub-second latency because it does not employ heavy-weight consensus. Instead, Brick uses consistent broadcast to announce updates and close the channel, a light-weight abstraction that is powerful enough to preserve safety and liveness to any rational parties. Furthermore, we consider permissioned blockchains, where the additional property of auditability might be desired for regulatory purposes. We introduce Brick+, an off-chain construction that provides auditability on top of Brick without conflicting with its privacy guarantees. We formally define the properties our payment channel construction should fulfill, and prove that both Brick and Brick+ satisfy them. We also design incentives for Brick such that honest and rational behavior aligns. Finally, we provide a reference implementation of the smart contracts in Solidity.


title: Efficient candidate screening under multiple tests and implications for fairness

abstract: When recruiting job candidates, employers rarely observe their underlying skill level directly. Instead, they must administer a series of interviews and/or collate other noisy signals in order to estimate the worker's skill. Traditional economics papers address screening models where employers access worker skill via a single noisy signal. In this paper, we extend this theoretical analysis to a multi-test setting, considering both Bernoulli and Gaussian models. We analyze the optimal employer policy both when the employer sets a fixed number of tests per candidate and when the employer can set a dynamic policy, assigning further tests adaptively based on results from the previous tests. To start, we characterize the optimal policy when employees constitute a single group, demonstrating some interesting trade-offs. Subsequently, we address the multi-group setting, demonstrating that when the noise levels vary across groups, a fundamental impossibility emerges whereby we cannot administer the same number of tests, subject candidates to the same decision rule, and yet realize the same outcomes in both groups.


title: Aspects of the Levi form

abstract: We discuss various analytical and geometrical aspects of the Levi form, which is associated with a CR manifold having any CR dimension and any CR codimension.


title: Complementation, Local Complementation, and Switching in Binary Matroids

abstract: In 2004, Ehrenfeucht, Harju, and Rozenberg showed that any graph on a vertex set $V$ can be obtained from a complete graph on $V$ via a sequence of the operations of complementation, switching edges and non-edges at a vertex, and local complementation. The last operation involves taking the complement in the neighbourhood of a vertex. In this paper, we consider natural generalizations of these operations for binary matroids and explore their behaviour. We characterize all binary matroids obtainable from the binary projective geometry of rank $r$ under the operations of complementation and switching. Moreover, we show that not all binary matroids of rank at most $r$ can be obtained from a projective geometry of rank $r$ via a sequence of the three generalized operations. We introduce a fourth operation and show that, with this additional operation, we are able to obtain all binary matroids.


title: Cosmological constraints and phenomenology of a beyond-Horndeski model

abstract: We study observational constraints on a specific dark energy model in the framework of Gleyzes-Langlois-Piazza-Vernizzi theories, which extends the Galileon ghost condensate (GGC) to the domain of beyond Horndeski theories. In this model, we show that the Planck cosmic microwave background (CMB) data, combined with datasets of baryon acoustic oscillations, supernovae type Ia, and redshift-space distortions, give the tight upper bound $|\alpha_{\rm H}^{(0)}| \le {\cal O}(10^{-6})$ on today's beyond-Horndeski (BH) parameter $\alpha_{\rm H}$. This is mostly attributed to the shift of CMB acoustic peaks induced by the early-time changes of cosmological background and perturbations arising from the dominance of $\alpha_{\rm H}$ in the dark energy density. In comparison to the $\Lambda$-cold-dark-matter ($\Lambda$CDM) model, our BH model suppresses the large-scale integrated-Sachs-Wolfe (ISW) tail of CMB temperature anisotropies due to the existence of cubic Galileons, and it modifies the small-scale CMB power spectrum because of the different background evolution. We find that the BH model considered fits the data better than $\Lambda$CDM according to the $\chi^2$ statistics, yet the deviance information criterion (DIC) slightly favors the latter. Given the fact that our BH model with $\alpha_{\rm H}=0$ (i.e., the GGC model) is favored over $\Lambda$CDM even by the DIC, there are no particular signatures for the departure from Horndeski theories in current observations.


title: Context-Aware Middleware: A Review

abstract: During previous years several studies have introduced the concept of "context-aware middleware" and also proposed solutions under this title; however, these systems are different in functionality. In this chapter, context-aware middleware is investigated from the standpoint of functional and non-functional requirements. Afterward, some well-known middleware systems are reviewed and, finally, open research directions as well as conclusion remarks are presented.


title: Supporting Software Engineering Research and Education by Annotating Public Videos of Developers Programming

abstract: Software engineering has long studied how software developers work, building a body of work which forms the foundation of many software engineering best practices, tools, and theories. Recently, some developers have begun recording videos of themselves engaged in programming tasks contributing to open source projects, enabling them to share knowledge and socialize with other developers. We believe that these videos offer an important opportunity for both software engineering research and education. In this paper, we discuss the potential use of these videos as well as open questions for how to best enable this envisioned use. We propose creating a central repository of programming videos, enabling analyzing and annotating videos to illustrate specific behaviors of interest such as asking and answering questions, employing strategies, and software engineering theories. Such a repository would offer an important new way in which both software engineering researchers and students can understand how software developers work.


title: Crosscap numbers of alternating knots via unknotting splices

abstract: Ito-Takimura recently defined a splice-unknotting number $u^-(D)$ for knot diagrams. They proved that this number provides an upper bound for the crosscap number of any prime knot, asking whether equality holds in the alternating case. We answer their question in the affirmative. (Ito has independently proven the same result.) As an application, we compute the crosscap numbers of all prime alternating knots through at least 13 crossings, using Gauss codes.


title: Simple and Effective Regularization Methods for Training on Noisily Labeled Data with Generalization Guarantee

abstract: Over-parameterized deep neural networks trained by simple first-order methods are known to be able to fit any labeling of data. Such over-fitting ability hinders generalization when mislabeled training examples are present. On the other hand, simple regularization methods like early-stopping can often achieve highly nontrivial performance on clean test data in these scenarios, a phenomenon not theoretically understood. This paper proposes and analyzes two simple and intuitive regularization methods: (i) regularization by the distance between the network parameters to initialization, and (ii) adding a trainable auxiliary variable to the network output for each training example. Theoretically, we prove that gradient descent training with either of these two methods leads to a generalization guarantee on the clean data distribution despite being trained using noisy labels. Our generalization analysis relies on the connection between wide neural network and neural tangent kernel (NTK). The generalization bound is independent of the network size, and is comparable to the bound one can get when there is no label noise. Experimental results verify the effectiveness of these methods on noisily labeled datasets.


title: Object Discovery with a Copy-Pasting GAN

abstract: We tackle the problem of object discovery, where objects are segmented for a given input image, and the system is trained without using any direct supervision whatsoever. A novel copy-pasting GAN framework is proposed, where the generator learns to discover an object in one image by compositing it into another image such that the discriminator cannot tell that the resulting image is fake. After carefully addressing subtle issues, such as preventing the generator from `cheating', this game results in the generator learning to select objects, as copy-pasting objects is most likely to fool the discriminator. The system is shown to work well on four very different datasets, including large object appearance variations in challenging cluttered backgrounds.


title: Coexistence of absolute negative mobility and anomalous diffusion

abstract: Using extensive numerical studies we demonstrate that absolute negative mobility of a Brownian particle (i.e. the net motion into the direction opposite to a constant biasing force acting around zero bias) does coexist with anomalous diffusion. The latter is characterized in terms of a nonlinear scaling with time of the mean-square deviation of the particle position. Such anomalous diffusion covers "coherent" motion (i.e. the position dynamics x(t) approaches in evolving time a constant dispersion), ballistic diffusion, subdiffusion, superdiffusion and hyperdiffusion. In providing evidence for this coexistence we consider a paradigmatic model of an inertial Brownian particle moving in a one-dimensional symmetric periodic potential being driven by both an unbiased time-periodic force and a constant bias. This very setup allows for various sorts of different physical realizations.


title: On the OA(1536,13,2,7) and related orthogonal arrays

abstract: With a computer-aided approach based on the connection with equitable partitions, we establish the uniqueness of the orthogonal array OA$(1536,13,2,7)$, constructed in [D.G.Fon-Der-Flaass. Perfect $2$-Colorings of a Hypercube, Sib. Math. J. 48 (2007), 740-745] as an equitable partition of the $13$-cube with quotient matrix $[[0,13],[3,10]]$. By shortening the OA$(1536,13,2,7)$, we obtain $3$ inequivalent orthogonal arrays OA$(768,12,2,6)$, which is a complete classification for these parameters too. After our computing, the first parameters of unclassified binary orthogonal arrays OA$(N,n,2,t)$ attending the Friedman bound $N\ge 2^n(1-n/2(t+1))$ are OA$(2048,14,2,7)$. Such array can be obtained by puncturing any binary $1$-perfect code of length $15$. We construct orthogonal arrays with these and similar parameters OA$(N=2^{n-m+1},n=2^m-2,2,t=2^{m-1}-1)$, $m\ge 4$, that are not punctured $1$-perfect codes. Additionally, we prove that any orthogonal array OA$(N,n,2,t)$ with even $t$ attending the bound $N \ge 2^n(1-(n+1)/2(t+2))$ induces an equitable $3$-partition of the $n$-cube.


title: The Different Regimes of Axion Gauge Field Inflation

abstract: In axion gauge field inflation an axion-like particle driving cosmic inflation is coupled to the Chern-Simons density of an Abelian or non-Abelian gauge group. In the case of a non-Abelian gauge group, this can lead to the formation of a stable, homogeneous and isotropic gauge field background. We study the dynamics of the inflaton and gauge fields in terms of the two effective coupling parameters: the gauge coupling and the axion decay constant. Starting from the Bunch-Davies vacuum in the far past, we find that the non-trivial gauge field background arises only significantly after the cosmic microwave background (CMB) scales have left the horizon. At these scales, the model thus closely resembles Abelian axion inflation, thereby naturally reconciling the tension of non-Abelian axion gauge field inflation with the latest CMB observations. We further consider two exemplary UV-completions of this setup: multiple Peccei-Quinn axions and axion monodromy in string theory. In both cases we find that the majority of the parameter space is excluded by theoretical or observational constraints. The remaining parameter space can be divided into three regimes. (i) For small gauge couplings we recover natural inflation. For large gauge couplings the non-Abelian gauge theory either (ii) mimics the Abelian theory or (iii) non-linear interactions prohibit a linear analysis of the gauge field perturbations.


title: Revisiting Stochastic Extragradient

abstract: We fix a fundamental issue in the stochastic extragradient method by providing a new sampling strategy that is motivated by approximating implicit updates. Since the existing stochastic extragradient algorithm, called Mirror-Prox, of (Juditsky et al., 2011) diverges on a simple bilinear problem when the domain is not bounded, we prove guarantees for solving variational inequality that go beyond existing settings. Furthermore, we illustrate numerically that the proposed variant converges faster than many other methods on bilinear saddle-point problems. We also discuss how extragradient can be applied to training Generative Adversarial Networks (GANs) and how it compares to other methods. Our experiments on GANs demonstrate that the introduced approach may make the training faster in terms of data passes, while its higher iteration complexity makes the advantage smaller.


title: A Unifying Causal Framework for Analyzing Dataset Shift-stable Learning Algorithms

abstract: Recent interest in the external validity of prediction models (i.e., the problem of different train and test distributions, known as dataset shift) has produced many methods for finding predictive distributions that are invariant to dataset shifts and can be used for prediction in new, unseen environments. However, these methods consider different types of shifts and have been developed under disparate frameworks, making it difficult to theoretically analyze how solutions differ with respect to stability and accuracy. Taking a causal graphical view, we use a flexible graphical representation to express various types of dataset shifts. Given a known graph of the data generating process, we show that all invariant distributions correspond to a causal hierarchy of graphical operators which disable the edges in the graph that are responsible for the shifts. The hierarchy provides a common theoretical underpinning for understanding when and how stability to shifts can be achieved, and in what ways stable distributions can differ. We use it to establish conditions for minimax optimal performance across environments, and derive new algorithms that find optimal stable distributions. Using this new perspective, we empirically demonstrate that that there is a tradeoff between minimax and average performance.


title: Probes of the Standard Model effective field theory extended with a right-handed neutrino

abstract: If neutrinos are Dirac particles and, as suggested by the so far null LHC results, any new physics lies at energies well above the electroweak scale, the Standard Model effective field theory has to be extended with operators involving the right-handed neutrinos. In this paper, we study this effective field theory and set constraints on the different dimension-six interactions. To that aim, we use LHC searches for associated production of light (and tau) leptons with missing energy, monojet searches, as well as pion and tau decays. Our bounds are generally above the TeV for order one couplings. One particular exception is given by operators involving top quarks. These provide new signals in top decays not yet studied at colliders. Thus, we also design an LHC analysis to explore these signatures in the $t\overline{t}$ production. Our results are also valid if the right-handed neutrinos are Majorana and long-lived.


title: Complex planar curves homeomorphic to a line have at most four singular points

abstract: We show that a complex planar curve homeomorphic to the projective line has at most four singular points. If it has exactly four then it has degree five and is unique up to a projective equivalence.


title: FlightGoggles: A Modular Framework for Photorealistic Camera, Exteroceptive Sensor, and Dynamics Simulation

abstract: FlightGoggles is a photorealistic sensor simulator for perception-driven robotic vehicles. The key contributions of FlightGoggles are twofold. First, FlightGoggles provides photorealistic exteroceptive sensor simulation using graphics assets generated with photogrammetry. Second, it provides the ability to combine (i) synthetic exteroceptive measurements generated in silico in real time and (ii) vehicle dynamics and proprioceptive measurements generated in motio by vehicle(s) in a motion-capture facility. FlightGoggles is capable of simulating a virtual-reality environment around autonomous vehicle(s). While a vehicle is in flight in the FlightGoggles virtual reality environment, exteroceptive sensors are rendered synthetically in real time while all complex extrinsic dynamics are generated organically through the natural interactions of the vehicle. The FlightGoggles framework allows for researchers to accelerate development by circumventing the need to estimate complex and hard-to-model interactions such as aerodynamics, motor mechanics, battery electrochemistry, and behavior of other agents. The ability to perform vehicle-in-the-loop experiments with photorealistic exteroceptive sensor simulation facilitates novel research directions involving, e.g., fast and agile autonomous flight in obstacle-rich environments, safe human interaction, and flexible sensor selection. FlightGoggles has been utilized as the main test for selecting nine teams that will advance in the AlphaPilot autonomous drone racing challenge. We survey approaches and results from the top AlphaPilot teams, which may be of independent interest.


title: The empirical aspects of Eliashberg formalism for the superconducting mechanism in Iron-based superconductors

abstract: We have investigated experimentally how properties of NdFeAsO0.8F0.2 superconductor affected due to the substitution of the Ca2+/Nd3+ doping. Based on the XRD data refinement, various structural parameters such as lattice parameters, bond angles, bond length, and etc. were studied. We have determined the upper limit of the calcium solubility in the NdFeAsO0.8F0.2 phase and it is restricted to x=0.05. Also, we have found that the lattice parameters and the cell volume decreased by increasing the calcium content. According to the XRD data analysis, we have argued that these reductions are due to the variations in the bond lengths and the bond angles of (O/F)-Nd-(O/F) and As-Fe-As i.e. alpha, beta upon increasing the calcium dopant. So, we have expected that the superconducting transition temperature (TC) will be sensitive to the calcium doping values. Experimentally, the TC of our samples was reduced from 53 K (for x = 0) to 48 K (for x = 0.01) and 27 K (for x=0.025) and disappeared for our other sample. Then we have studied the dependence of TC and bond angles, bond length, the pnictogen height, and the lattice parameter to examine the available theories from an empirical point of view. The consistency of our experimental results and the theoretical reports based on the spin- and the orbital-fluctuation theories shows that these models play an important role in the pairing mechanism of the iron-based superconductors.


title: A New Non-Linear Conjugate Gradient Algorithm for Destructive Cure Rate Model and a Simulation Study: Illustration with Negative Binomial Competing Risks

abstract: In this paper, we propose a new estimation methodology based on a projected non-linear conjugate gradient (PNCG) algorithm with an efficient line search technique. We develop a general PNCG algorithm for a survival model incorporating a proportion cure under a competing risks setup, where the initial number of competing risks are exposed to elimination after an initial treatment (known as destruction). In the literature, expectation maximization (EM) algorithm has been widely used for such a model to estimate the model parameters. Through an extensive Monte Carlo simulation study, we compare the performance of our proposed PNCG with that of the EM algorithm and show the advantages of our proposed method. Through simulation, we also show the advantages of our proposed methodology over other optimization algorithms (including other conjugate gradient type methods) readily available as R software packages. To show these we assume the initial number of competing risks to follow a negative binomial distribution although our general algorithm allows one to work with any competing risks distribution. Finally, we apply our proposed algorithm to analyze a well-known melanoma data.


title: On Star-critical (K1,n,K1,m + e) Ramsey numbers

abstract: Let $G, H$ be finite graphs without loops or multiple edges and $K_n$ denote the complete graph on $n$ vertices. If for every red/blue colouring of edges of the complete graph $K_n$, there exists a red copy of $G$, or a blue copy of $H$, we will say that $K_n\rightarrow (G,H)$. The Ramsey number $r(G, H)$ is defined as the smallest positive integer $n$ such that $K_{n} \rightarrow (G, H)$. Star-critical Ramsey number $r_*(G, H)$ is defined as the largest value of $k$ such that $K_{r(G,H)-1} \sqcup K_{1,k} \rightarrow (G, H)$. In this paper, we will find $r_*(K_{1,n}, K_{1,m}+e)$ for all $n,m \geq 3$.


title: Trust but Verify: An Information-Theoretic Explanation for the Adversarial Fragility of Machine Learning Systems, and a General Defense against Adversarial Attacks

abstract: Deep-learning based classification algorithms have been shown to be susceptible to adversarial attacks: minor changes to the input of classifiers can dramatically change their outputs, while being imperceptible to humans. In this paper, we present a simple hypothesis about a feature compression property of artificial intelligence (AI) classifiers and present theoretical arguments to show that this hypothesis successfully accounts for the observed fragility of AI classifiers to small adversarial perturbations. Drawing on ideas from information and coding theory, we propose a general class of defenses for detecting classifier errors caused by abnormally small input perturbations. We further show theoretical guarantees for the performance of this detection method. We present experimental results with (a) a voice recognition system, and (b) a digit recognition system using the MNIST database, to demonstrate the effectiveness of the proposed defense methods. The ideas in this paper are motivated by a simple analogy between AI classifiers and the standard Shannon model of a communication system.


title: State-Reification Networks: Improving Generalization by Modeling the Distribution of Hidden Representations

abstract: Machine learning promises methods that generalize well from finite labeled data. However, the brittleness of existing neural net approaches is revealed by notable failures, such as the existence of adversarial examples that are misclassified despite being nearly identical to a training example, or the inability of recurrent sequence-processing nets to stay on track without teacher forcing. We introduce a method, which we refer to as \emph{state reification}, that involves modeling the distribution of hidden states over the training data and then projecting hidden states observed during testing toward this distribution. Our intuition is that if the network can remain in a familiar manifold of hidden space, subsequent layers of the net should be well trained to respond appropriately. We show that this state-reification method helps neural nets to generalize better, especially when labeled data are sparse, and also helps overcome the challenge of achieving robust generalization with adversarial training.


title: Elliptic solutions to integrable nonlinear equations and many-body systems

abstract: We review elliptic solutions to integrable nonlinear partial differential and difference equations (KP, mKP, BKP, Toda) and derive equations of motion for poles of the solutions. The pole dynamics is given by an integrable many-body system (Calogero-Moser, Ruijsenaars-Schneider). The basic tool is the auxiliary linear problems for the wave function which yield equations of motion together with their Lax representation. We also discuss integrals of motion and properties of the spectral curves.


title: Partial minimization of strict convex functions and tensor scaling

abstract: Assume that f is a strict convex function with a unique minimum in R^n. We divide the vector of n-variables to d groups of vector subvariables with d at least two. We assume that we can find the partial minimum of f with respect to each vector subvariable while other variables are fixed. We then describe an algorithm that partially minimizes each time on a specifically chosen vector subvariable. This algorithm converges geometrically to the unique minimum. The rate of convergence depends on the uniform bounds on the eigenvalues of the Hessian of f in the compact sublevel set f whose values are at most f(x_0), where x_0 is the starting point of the algorithm. In the case where f is a polynomial of degree two, with positive definite quadratic term, and d=n our method can be considered as a generalization of the classical conjugate gradient method. The main result of this paper is the observation that the celebrated Sinkhorn diagonal scaling algorithm for matrices, and the corresponding diagonal scaling of tensors, can be viewed as partial minimization of certain logconvex functions.


title: Investigate the effects of the existence of correlation between two consecutive use of the quantum channel on quantum speed limit time

abstract: Memory effects play an important role in the theory of open quantum systems. There are two completely independent insights about memory for quantum channels. In quantum information theory, the memory of the quantum channel is depicted by the correlations between consecutive uses of the channel on a set of quantum systems. In the theory of open quantum systems memory effects result from correlations which are created during the quantum evolution. Here, we study the quantum speed limit time for correlated quantum channel i.e. when there exist correlation between consecutive uses of quantum channel . Quantum speed limit time is the bound on the minimal time evolution between initial and target states. It is apply for quantifying the maximum speed of quantum evolution. In this work, we will consider correlated pure dephasing colored noise as an example of unital quantum channels and correlated squeezed generalized amplitude damping channel as an example of non-unital quantum channels. We evaluate the quantum speed limit time for these two channel. We find that regardless of whether the channel is unital or not, the quantum speed limit time is increased by increasing correlation between two consecutive uses of the channel.


title: Large Sample Properties of Matching for Balance

abstract: Matching methods are widely used for causal inference in observational studies. Among them, nearest neighbor matching is arguably the most popular. However, nearest neighbor matching does not generally yield an average treatment effect estimator that is $\sqrt{n}$-consistent (Abadie and Imbens, 2006). Are matching methods not $\sqrt{n}$-consistent in general? In this paper, we study a recent class of matching methods that use integer programming to directly target aggregate covariate balance as opposed to finding close neighbor matches. We show that under suitable conditions these methods can yield simple estimators that are $\sqrt{n}$-consistent and asymptotically optimal.


title: Automatic Delineation of Kidney Region in DCE-MRI

abstract: Delineation of the kidney region in dynamic contrast-enhanced magnetic resonance Imaging (DCE-MRI) is required during post-acquisition analysis in order to quantify various aspects of renal function, such as filtration and perfusion or blood flow. However, this can be obfuscated by the Partial Volume Effect (PVE), caused due to the mixing of any single voxel with two or more signal intensities from adjacent regions such as liver region and other tissues. To avoid this problem, firstly, a kidney region of interest (ROI) needs to be defined for the analysis. A clinician may choose to select a region avoiding edges where PV mixing is likely to be significant. However, this approach is time-consuming and labour intensive. To address this issue, we present Dynamic Mode Decomposition (DMD) coupled with thresholding and blob analysis as a framework for automatic delineation of the kidney region. This method is first validated on synthetically generated data with ground-truth available and then applied to ten healthy volunteers' kidney DCE-MRI datasets. We found that the result obtained from our proposed framework is comparable to that of a human expert. For example, while our result gives an average Root Mean Square Error (RMSE) of 0.0097, the baseline achieves an average RMSE of 0.1196 across the 10 datasets. As a result, we conclude automatic modelling via DMD framework is a promising approach.


title: Locating fixed points in the phase plane

abstract: The critical point is a fixed point in finite-size scaling. To quantify the behaviour of such a fixed point, we define, at a given temperature and scaling exponent ratio, the width of scaled observables for different sizes. The minimum of the width reveals the position of fixed point, its corresponding phase transition temperature, and scaling exponent ratio. The value of this ratio tells the nature of fixed point, which can be a critical point, a point of the first order phase transition line, or a point of the crossover region. To demonstrate the effectiveness of this method, we apply it to three typical samples produced by the 3D three-state Potts model. Results show the method to be more precise and effective than conventional methods. Finally, we discuss a possible application at the beam energy scan program of relativistic heavy-Ion collision.


title: Direct and inverse theorems of approximation of functions in weighted Orlicz type spaces with variable exponent

abstract: In weighted Orlicz type spaces ${\mathcal S}_{_{\scriptstyle \mathbf p,\,\mu}}$ with a variable summation exponent, the direct and inverse approximation theorems are proved in terms of best approximations of functions and moduli of smoothness of fractional order. It is shown that the constant obtained in the inverse approximation theorem is in a certain sense the best. Some applications of the results are also proposed. In particular, the constructive characteristics of functional classes defined by such moduli of smoothness are given. Equivalence between moduli of smoothness and certain Peetre $K$-functionals is shown in the spaces ${\mathcal S}_{_{\scriptstyle \mathbf p,\,\mu}}$.


title: Towards the full classification of exceptional scattered polynomials

abstract: Let $f(X) \in \mathbb{F}_{q^r}[X]$ be a $q$-polynomial. If the $\mathbb{F}_q$-subspace $U=\{(x^{q^t},f(x)) \mid x \in \mathbb{F}_{q^n}\}$ defines a maximum scattered linear set, then we call $f(X)$ a scattered polynomial of index $t$. The asymptotic behaviour of scattered polynomials of index $t$ is an interesting open problem. In this sense, exceptional scattered polynomials of index $t$ are those for which $U$ is a maximum scattered linear set in ${\rm PG}(1,q^{mr})$ for infinitely many $m$. The complete classifications of exceptional scattered monic polynomials of index $0$ (for $q>5$) and of index 1 were obtained by Bartoli and Zhou. In this paper we complete the classifications of exceptional scattered monic polynomials of index $0$ for $q \leq 4$. Also, some partial classifications are obtained for arbitrary $t$. As a consequence, the complete classification of exceptional scattered monic polynomials of index $2$ is given.


title: Using Neural Networks for Relation Extraction from Biomedical Literature

abstract: Using different sources of information to support automated extracting of relations between biomedical concepts contributes to the development of our understanding of biological systems. The primary comprehensive source of these relations is biomedical literature. Several relation extraction approaches have been proposed to identify relations between concepts in biomedical literature, namely, using neural networks algorithms. The use of multichannel architectures composed of multiple data representations, as in deep neural networks, is leading to state-of-the-art results. The right combination of data representations can eventually lead us to even higher evaluation scores in relation extraction tasks. Thus, biomedical ontologies play a fundamental role by providing semantic and ancestry information about an entity. The incorporation of biomedical ontologies has already been proved to enhance previous state-of-the-art results.


title: Successive Cancellation List Decoding of Semi-random Unit Memory Convolutional Codes

abstract: We present in this paper a special class of unit memory convolutional codes (UMCCs), called semi-random UMCCs (SRUMCCs), where the information block is first encoded by a short block code and then transmitted in a block Markov (random) superposition manner. We propose a successive cancellation list decoding algorithm, by which a list of candidate codewords are generated serially until one passes an empirical divergence test instead of the conventional cyclic redundancy check (CRC). The threshold for testing the correctness of candidate codewords can be learned off-line based on the statistical behavior of the introduced empirical divergence function (EDF). The performance-complexity tradeoff and the performance-delay tradeoff can be achieved by adjusting the statistical threshold and the decoding window size. To analyze the performance, a closed-form upper bound and a simulated lower bound are derived. Simulation results verify our analysis and show that: 1) The proposed list decoding algorithm with empirical divergence test outperforms the sequential decoding in high signal-to-noise ratio (SNR) region; 2) Taking the tail-biting convolutional codes (TBCC) as the basic codes, the proposed list decoding of SRUMCCs have comparable performance with the polar codes under the constraint of equivalent decoding delay.


title: A Self-Attention Joint Model for Spoken Language Understanding in Situational Dialog Applications

abstract: Spoken language understanding (SLU) acts as a critical component in goal-oriented dialog systems. It typically involves identifying the speakers intent and extracting semantic slots from user utterances, which are known as intent detection (ID) and slot filling (SF). SLU problem has been intensively investigated in recent years. However, these methods just constrain SF results grammatically, solve ID and SF independently, or do not fully utilize the mutual impact of the two tasks. This paper proposes a multi-head self-attention joint model with a conditional random field (CRF) layer and a prior mask. The experiments show the effectiveness of our model, as compared with state-of-the-art models. Meanwhile, online education in China has made great progress in the last few years. But there are few intelligent educational dialog applications for students to learn foreign languages. Hence, we design an intelligent dialog robot equipped with different scenario settings to help students learn communication skills.


title: An Accelerated Decentralized Stochastic Proximal Algorithm for Finite Sums

abstract: Modern large-scale finite-sum optimization relies on two key aspects: distribution and stochastic updates. For smooth and strongly convex problems, existing decentralized algorithms are slower than modern accelerated variance-reduced stochastic algorithms when run on a single machine, and are therefore not efficient. Centralized algorithms are fast, but their scaling is limited by global aggregation steps that result in communication bottlenecks. In this work, we propose an efficient \textbf{A}ccelerated \textbf{D}ecentralized stochastic algorithm for \textbf{F}inite \textbf{S}ums named ADFS, which uses local stochastic proximal updates and randomized pairwise communications between nodes. On $n$ machines, ADFS learns from $nm$ samples in the same time it takes optimal algorithms to learn from $m$ samples on one machine. This scaling holds until a critical network size is reached, which depends on communication delays, on the number of samples $m$, and on the network topology. We provide a theoretical analysis based on a novel augmented graph approach combined with a precise evaluation of synchronization times and an extension of the accelerated proximal coordinate gradient algorithm to arbitrary sampling. We illustrate the improvement of ADFS over state-of-the-art decentralized approaches with experiments.


title: Multi-Modal Graph Interaction for Multi-Graph Convolution Network in Urban Spatiotemporal Forecasting

abstract: Graph convolution network based approaches have been recently used to model region-wise relationships in region-level prediction problems in urban computing. Each relationship represents a kind of spatial dependency, like region-wise distance or functional similarity. To incorporate multiple relationships into spatial feature extraction, we define the problem as a multi-modal machine learning problem on multi-graph convolution networks. Leveraging the advantage of multi-modal machine learning, we propose to develop modality interaction mechanisms for this problem, in order to reduce generalization error by reinforcing the learning of multimodal coordinated representations. In this work, we propose two interaction techniques for handling features in lower layers and higher layers respectively. In lower layers, we propose grouped GCN to combine the graph connectivity from different modalities for more complete spatial feature extraction. In higher layers, we adapt multi-linear relationship networks to GCN by exploring the dimension transformation and freezing part of the covariance structure. The adapted approach, called multi-linear relationship GCN, learns more generalized features to overcome the train-test divergence induced by time shifting. We evaluated our model on ridehailing demand forecasting problem using two real-world datasets. The proposed technique outperforms state-of-the art baselines in terms of prediction accuracy, training efficiency, interpretability and model robustness.


title: A State Observer Design for Simultaneous Estimation of Charge State and Crossover in Self-Discharging Disproportionation Redox Flow Batteries

abstract: This paper presents an augmented state observer design for the simultaneous estimation of charge state and crossover flux in disproportionation redox flow batteries, which exhibits exponential estimation error convergence to a bounded residual set. The crossover flux of vanadium through the porous separator is considered as an unknown function of the battery states, model-approximated as the output of a persistently excited linear system. This parametric model and the simple isothermal lumped parameter model of the battery are combined to form an augmented space state representation suitable for the observer design, which is carried out via Lyapunov stability theory including the error-uncertainty involved in the approximation of the crossover flux. The observer gain is calculated by solving a polytopic linear matrix inequality problem via convex optimization. The performance of this design is evaluated with a laboratory flow battery prototype undergoing self-discharge.


title: Are sample means in multi-armed bandits positively or negatively biased?

abstract: It is well known that in stochastic multi-armed bandits (MAB), the sample mean of an arm is typically not an unbiased estimator of its true mean. In this paper, we decouple three different sources of this selection bias: adaptive \emph{sampling} of arms, adaptive \emph{stopping} of the experiment, and adaptively \emph{choosing} which arm to study. Through a new notion called ``optimism'' that captures certain natural monotonic behaviors of algorithms, we provide a clean and unified analysis of how optimistic rules affect the sign of the bias. The main takeaway message is that optimistic sampling induces a negative bias, but optimistic stopping and optimistic choosing both induce a positive bias. These results are derived in a general stochastic MAB setup that is entirely agnostic to the final aim of the experiment (regret minimization or best-arm identification or anything else). We provide examples of optimistic rules of each type, demonstrate that simulations confirm our theoretical predictions, and pose some natural but hard open problems.


title: New decomposition formulas associated with the Lauricella multivariable hypergeometric functions

abstract: Decomposition formulas associated with the Lauricella multivariable hypergeometric functions were known, however, due to the recurrence of those formulas, additional difficulties may arise in the applications. Further study of the properties of the famous expansion formulas showed that it can be reduced to a more convenient form. In addition, this paper contains applications of new expansion formulas to the solving of boundary value problems for a multidimensional elliptic equation with several singular coefficients.


title: Nuclear-matter distribution in the proton-rich nuclei $^7$Be and $^8$B from intermediate energy proton elastic scattering in inverse kinematics

abstract: Absolute differential cross sections for elastic $p^7$Be and $p^8$B small-angle scattering were measured in inverse kinematics at an energy of 0.7 GeV/u at GSI Darmstadt. The hydrogen-filled ionization chamber IKAR was used as an active target to detect the recoil protons. The projectile tracking and isotope identification were performed with multi-wire proportional chambers and scintillation detectors. The measured cross sections were analysed using the Glauber multiple-scattering theory. The root-mean-square (rms) nuclear matter radii $R_{\rm m} = 2.42 (4)$ fm for $^7$Be and $R_{\rm m} = 2.58 (6)$ fm for $^8$B were obtained. The radial density distribution deduced for $^8$B exhibits a proton halo structure with the rms halo radius $R_{\rm h} = 4.24 (25)$ fm. A comparison of the deduced experimental radii is displayed with existing experimental and theoretical data.


title: Vibrations of an elastic bar, isospectral deformations, and modified Camassa-Holm equations

abstract: Vibrations of an elastic rod are described by a Sturm-Liouville system. We present a general discussion of isospectral (spectrum preserving) deformations of such a system. We interpret one family of such deformations in terms of a two-component modified Camassa-Holm equation (2-mCH) and solve completely its dynamics for the case of discrete measures (multipeakons). We show that the underlying system is Hamiltonian and prove its Liouville integrability. The present paper generalizes our previous work on interlacing multipeakons of the 2-mCH and multipeakons of the 1-mCH. We give a unified approach to both equations, emphasizing certain natural family of interpolation problems germane to the solution of the inverse problem for 2-mCH as well as to this type of a Sturm-Liouville system with singular coefficients.


title: Flexible and inflexible $CR$ submanifolds

abstract: In this paper we prove new embedding results for compactly supported deformations of $CR$ submanifolds of $\mathbb{C}^{n+d}$: We show that if $M$ is a $2$-pseudoconcave $CR$ submanifold of type $(n,d)$ in $\mathbb{C}^{n+d}$, then any compactly supported $CR$ deformation stays in the space of globally $CR$ embeddable in $\mathbb{C}^{n+d}$ manifolds. This improves an earlier result, where $M$ was assumed to be a quadratic $2$-pseudoconcave $CR$ submanifold of $\mathbb{C}^{n+d}$. We also give examples of weakly $2$-pseudoconcave $CR$ manifolds admitting compactly supported $CR$ deformations that are not even locally $CR$ embeddable.


title: Lyman continuum observations across cosmic time: recent developments, future requirements

abstract: Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.


title: The Curious Case of NiRh$_2$O$_4$: A Spin-Orbit Entangled Diamond Lattice Paramagnet

abstract: Motivated by the interest in topological quantum paramagnets in candidate spin-$1$ magnets, we investigate the diamond lattice compound NiRh$_2$O$_4$ using {\it ab initio} theory and model Hamiltonian approaches. Our density functional study, taking into account the unquenched orbital degrees of freedom, shows stabilization of $S=1$ and $L=1$ state. We highlight the importance of spin-orbit coupling, in addition to Coulomb correlations, in driving the insulating gap, and uncover frustrating large second-neighbor exchange mediated by Ni-Rh covalency. A single-site model Hamiltonian incorporating the large tetragonal distortion is shown to give rise to a spin-orbit entangled non-magnetic ground state, largely accounting for the entropy, magnetic susceptibility, and inelastic neutron scattering results. Incorporating inter-site exchange within a slave-boson theory, we show that exchange frustration can suppress exciton condensation. We capture the dispersive gapped magnetic modes, uncover `dark states' invisible to neutrons, and make predictions for future experiments.


title: Life in the fast lane: a direct view of the dynamics, formation, and evolution of the Milky Way's bar

abstract: Studies of the ages, abundances, and motions of individual stars in the Milky Way provide one of the best ways to study the evolution of disk galaxies over cosmic time. The formation of the Milky Way's barred inner region in particular is a crucial piece of the puzzle of disk galaxy evolution. Using data from APOGEE and Gaia, we present maps of the kinematics, elemental abundances, and age of the Milky Way bulge and disk that show the barred structure of the inner Milky Way in unprecedented detail. The kinematic maps allow a direct, purely kinematic determination of the bar's pattern speed of 41+/-3 km/s/kpc and of its shape and radial profile. We find the bar's age, metallicity, and abundance ratios to be the same as those of the oldest stars in the disk that are formed in its turbulent beginnings, while stars in the bulge outside of the bar are younger and more metal-rich. This implies that the bar likely formed ~8 Gyr ago, when the decrease in turbulence in the gas disk allowed a thin disk to form that quickly became bar-unstable. The bar's formation therefore stands as a crucial epoch in the evolution of the Milky Way, a picture that is in line with the evolutionary path that emerges from observations of the gas kinematics in external disk galaxies over the last ~10 Gyr.


title: Quantized Fermi-arc-mediated transport in Weyl semimetal nanowires

abstract: We study longitudinal transport in Weyl semimetal nanowires, both in the absence and in the presence of a magnetic flux threading the nanowires. We identify two qualitatively different regimes of transport with respect to the chemical potential in the nanowires. In the "surface regime", for low doping, most of the conductance occurs through the Fermi-arc surface states, and it rises in steps of one quantum of conductance as a function of the chemical potential; furthermore, with varying flux the conductance changes in steps of one quantum of conductance with characteristic Fabry-P\'erot interference oscillations. In the "bulk-surface regime", for highly-doped samples, the dominant contribution to the conductance is quadratic in the chemical potential, and mostly conditioned by the bulk states; the flux dependence shows clearly that both the surface and the bulk states contribute. The two aforementioned regimes prove that the contribution of the Fermi-arc surface states is salient and, therefore, crucial for understanding transport properties of finite-size Weyl semimetal systems. Last but not least, we demonstrate that both regimes are robust to disorder.


title: Particle motion near high-spin black holes

abstract: General relativity predicts that the Kerr black hole develops qualitatively new and surprising features in the limit of maximal spin. Most strikingly, the region of spacetime near the event horizon stretches into an infinitely long throat and displays an emergent conformal symmetry. Understanding dynamics in this NHEK (Near-Horizon Extreme Kerr) geometry is necessary for connecting theory to upcoming astronomical observations of high-spin black holes. We review essential properties of NHEK and its relationship to the rapidly rotating Kerr black hole. We then completely solve the geodesic equation in the NHEK region and describe how the resulting trajectories transform under the action of its enhanced symmetries. In the process, we derive explicit expressions for the angular integrals appearing in the Kerr geodesic equation and obtain a useful formula, valid at arbitrary spin, for a particle's polar angle in terms of its radial motion. These results will aid in the analytic computation of astrophysical observables relevant to ongoing and future experiments.


title: Dark matter seeping through dynamic gauge kinetic mixing

abstract: We show for the first time that the loop-driven kinetic mixing between visible and dark Abelian gauge bosons can facilitate dark matter production in the early Universe by creating a 'dynamic' portal, which depends on the energy of the process. The required smallness of the strength of the portal interaction, suited for freeze-in, is justified by a suppression arising from the mass of a heavy vector-like fermion. The strong temperature sensitivity associated with the interaction is responsible for most of the dark matter production during the early stages of reheating.


title: Magnetic Field Induced Competing Phases in Spin-Orbital Entangled Kitaev Magnets

abstract: There has been a great interest in magnetic field induced quantum spin liquids in Kitaev magnets after the discovery of neutron scattering continuum and half quantized thermal Hall conductivity in the material $\alpha$-RuCl$_3$. In this work, we provide a semiclassical analysis of the relevant theoretical models on large system sizes, and compare the results to previous studies on quantum models with small system sizes. We find a series of competing magnetic orders with fairly large unit cells at intermediate magnetic fields, which are most likely missed by previous approaches. We show that quantum fluctuations are typically strong in these large unit cell orders, while their magnetic excitations may resemble a scattering continuum and give rise to a large thermal Hall conductivity. Our work provides an important basis for a thorough investigation of emergent spin liquids and competing phases in Kitaev magnets.


title: Chiral soliton lattice in QCD-like theories

abstract: Recently, it has been shown that the ground state of quantum chromodynamics (QCD) in sufficiently strong magnetic fields and at moderate baryon number chemical potential carries a crystalline condensate of neutral pions: the chiral soliton lattice (CSL). While the result was obtained in a model-independent manner using effective field theory techniques, its realization from first principles using lattice Monte Carlo simulation is hampered by the infamous sign problem. Here we show that CSL, or a similar inhomogeneous phase, also appears in the phase diagram of a class of vector-like gauge theories that do not suffer from the sign problem even in the presence of a baryon chemical potential and external magnetic field. We also show that the onset of nonuniform order manifests itself already in the adjacent homogeneous Bose-Einstein-condensation phase through a characteristic roton-like minimum in the dispersion relation of the lowest-lying quasiparticle mode. Last but not least, our work gives a class of explicit counterexamples to the long-standing conjecture that positivity of the determinant of the Dirac operator (that is, absence of the sign problem) in a vector-like gauge theory precludes spontaneous breaking of translational invariance, and thus implies the absence of inhomogeneous phases in the phase diagram of the theory.


title: Brown Dwarf Atmospheres as the Potentially Most Detectable and Abundant Sites for Life

abstract: We show that the total habitable volume in the atmospheres of cool brown dwarfs with effective temperatures of $\sim 250$-$350$ K is possibly larger by two orders of magnitude than that of Earth-like planets. We also study the role of aerosols, nutrients and photosynthesis in facilitating life in brown dwarf atmospheres. Our predictions might be testable through searches for spectral edges in the near-infrared and chemical disequilibrium in the atmospheres of nearby brown dwarfs that are either free-floating or within several AU of stars. For the latter category, we find that the James Webb Space Telescope (JWST) may be able to achieve a signal-to-noise ratio of $\sim 5$ after a few hours of integration time per source for the detection of biogenic spectral features in $\sim 10^3$ cool brown dwarfs.


title: Higher-order topological odd-parity superconductors

abstract: The topological property of a gapped odd-parity superconductor is jointly determined by its pairing nodes and Fermi surfaces in normal state. We reveal that the contractibility of Fermi surfaces without crossing any time-reversal invariant momentum and the presence of nontrivial Berry phase on Fermi surfaces are two key conditions for the realization of higher-order topological odd-parity superconductors (HOTOPSCs). When the normal state is a normal metal, we reveal the necessity of removable Dirac pairing nodes and provide a general and simple principle to realize HOTOPSCs. Our findings can not only be applied to analyze the topological property of odd-parity superconductors, but also be used as a guiding principle to design new platforms of higher-order topological superconductors, as well as higher-order topological insulators owing to their direct analogy in Hamiltonian description.


title: The Hubbard model on a triangular lattice -- pseudogap due to spin-density-wave fluctuations

abstract: We calculate the fermionic spectral function $A_k (\omega)$ in the spiral spin-density-wave (SDW) state of the Hubbard model on a quasi-2D triangular lattice at small but finite temperature $T$. The spiral SDW order $\Delta (T)$ develops below $T = T_N$ and has momentum ${ \bf K} = (4\pi/3,0)$. We pay special attention to fermions with momenta ${\bf k}$, for which ${\bf k}$ and ${\bf k} + {\bf K}$ are close to Fermi surface in the absence of SDW. At the mean field level, $A_k (\omega)$ for such fermions has peaks at $\omega = \pm \Delta (T)$ at $T < T_N$ and displays a conventional Fermi liquid behavior at $T > T_N$. We show that this behavior changes qualitatively beyond mean-field due to singular self-energy contributions from thermal fluctuations in a quasi-2D system. We use a non-perturbative eikonal approach and sum up infinite series of thermal self-energy terms. We show that $A_k (\omega)$ shows peak/dip/hump features at $T < T_N$, with the peak position at $\Delta (T)$ and hump position at $\Delta (T=0)$. Above $T_N$, the hump survives up to $T = T_p > T_N$, and in between $T_N$ and $T_p$ the spectral function displays the pseudogap behavior. We show that the difference between $T_p$ and $T_N$ is controlled by the ratio of in-plane and out-of-plane static spin susceptibilities, which determines the combinatoric factors in the diagrammatic series for the self-energy. For certain values of this ratio, $T_p = T_N$, i.e., the pseudogap region collapses. In this last case, thermal fluctuations are logarithmically singular, yet they do not give rise to pseudogap behavior. Our computational method can be used to study pseudogap physics due to thermal fluctuations in other systems.


title: Kinematics of $z\geq 6$ galaxies from [CII] line emission

abstract: We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$\mu$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5\sigma$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects.


title: How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

abstract: We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...