title: CNN-based Facial Affect Analysis on Mobile Devices

abstract: This paper focuses on the design, deployment and evaluation of Convolutional Neural Network (CNN) architectures for facial affect analysis on mobile devices. Unlike traditional CNN approaches, models deployed to mobile devices must minimise storage requirements while retaining high performance. We therefore propose three variants of established CNN architectures and comparatively evaluate them on a large, in-the-wild benchmark dataset of facial images. Our results show that the proposed architectures retain similar performance to the dataset baseline while minimising storage requirements: achieving 58% accuracy for eight-class emotion classification and average RMSE of 0.39 for valence/arousal prediction. To demonstrate the feasibility of deploying these models for real-world applications, we implement a music recommendation interface based on predicted user affect. Although the CNN models were not trained in the context of music recommendation, our case study shows that: (i) the trained models achieve similar prediction performance to the benchmark dataset, and (ii) users tend to positively rate the song recommendations provided by the interface. Average runtime of the deployed models on an iPhone 6S equates to ~45 fps, suggesting that the proposed architectures are also well suited for real-time deployment on video streams.


title: Peeking Behind Objects: Layered Depth Prediction from a Single Image

abstract: While conventional depth estimation can infer the geometry of a scene from a single RGB image, it fails to estimate scene regions that are occluded by foreground objects. This limits the use of depth prediction in augmented and virtual reality applications, that aim at scene exploration by synthesizing the scene from a different vantage point, or at diminished reality. To address this issue, we shift the focus from conventional depth map prediction to the regression of a specific data representation called Layered Depth Image (LDI), which contains information about the occluded regions in the reference frame and can fill in occlusion gaps in case of small view changes. We propose a novel approach based on Convolutional Neural Networks (CNNs) to jointly predict depth maps and foreground separation masks used to condition Generative Adversarial Networks (GANs) for hallucinating plausible color and depths in the initially occluded areas. We demonstrate the effectiveness of our approach for novel scene view synthesis from a single image.


title: Two Algorithms to Find Primes in Patterns

abstract: Let $k\ge 1$ be an integer, and let $P= (f_1(x), \ldots, f_k(x) )$ be $k$ admissible linear polynomials over the integers, or \textit{the pattern}. We present two algorithms that find all integers $x$ where $\max{ \{f_i(x) \} } \le n$ and all the $f_i(x)$ are prime. Our first algorithm takes at most $O_P(n/(\log\log n)^k)$ arithmetic operations using $O(k\sqrt{n})$ space. Our second algorithm takes slightly more time, $O_P(n/(\log \log n)^{k-1})$ arithmetic operations, but uses only $n^{1/c}$ space for a constant $c>2$. We prove correctness unconditionally, but the running time relies on two unproven but reasonable conjectures. We are unaware of any previous complexity results for this problem beyond the use of a prime sieve. We also implemented several parallel versions of our second algorithm to show it is viable in practice. In particular, we found some new Cunningham chains of length 15, and we found all quadruplet primes up to $10^{17}$.


title: Discrete-time model for a substance motion in a channel of a network. Application to a human migration channel

abstract: We discuss a discrete-time model for motion of substance in a channel of a network. For the case of stationary motion of the substance and for the case of time-independent values of the parameters of the model we obtain a new class of statistical distributions that describe the distribution of the substance along the nodes of the channel. The case of interaction between a kind of substance specific for a node of the network and another kind of substance that is leaked from the channel is studied in presence of possibility for conversion between the two substances. Several scenarios connected to the dynamics of the two kinds of substances are described. The studied models: (i) model of motion of substance through a channel of a network, and (ii) model of interaction between two kinds of substances in a network node connected to the channel, are discussed from the point of view of human migration dynamics and interaction between the population of migrants and the native population of a country.


title: A quantum Johnson-Lindenstrauss lemma via unitary t-designs

abstract: The famous Johnson-Lindenstrauss lemma states that for any set of n vectors, there is a linear transformation into a space of dimension O(log n) that approximately preserves all their lengths. In fact, a Haar random unitary transformation followed by projection onto the first O(log n) coordinates followed by a scaling works as a valid transformation with high probability. In this work, we show that the Haar random unitary can be replaced by a uniformly random unitary chosen from a finite set called an approximate unitary t-design for t = O(log n). Choosing a unitary from such a design requires only polylogarithmic random bits as opposed to exponential in dimension random bits required to choose a Haar random unitary with reasonable precision. Moreover, since such unitaries can be efficiently implemented in the superpositional setting, our result can be viewed as an efficient quantum Johnson-Lindenstrauss transform akin to efficient quantum Fourier transforms widely used in earlier work on quantum algorithms. We prove our result by leveraging a method of Low for showing concentration for approximate unitary t-designs. We discuss algorithmic advantages and limitations of our result and conclude with a toy application to private information retrieval.


title: Khovanskii-Finite Valuations, Rational Curves, and Torus Actions

abstract: We study full rank homogeneous valuations on (multi)-graded domains and ask when they have finite Khovanskii bases. We show that there is a natural reduction from multigraded to simply graded domains. As special cases, we consider projective coordinate rings of rational curves, and almost toric varieties. Our results relate to several problems posed by Kaveh and Manon, and imply that the procedure of Bossinger-Lamboglia-Mincheva-Mohammadi for producing tropical prime cones will not terminate in general.


title: Density waves and jet emission asymmetry in Bose Fireworks

abstract: A Bose condensate subject to a periodic modulation of the two-body interactions was recently observed to emit matter-wave jets resembling "fireworks" [Nature 551, 356(2017)]. In this paper, combining experiment with numerical simulation, we demonstrate that these "Bose fireworks" represent a late stage in a complex time evolution of the driven condensate. We identify a "density wave" stage which precedes jet emission and results from interference of matterwaves. The density waves self-organize and self-amplify without the breaking of long range translational symmetry. Importantly, this density wave structure deterministically establishes the template for the subsequent patterns of the emitted jets. Our simulations, in good agreement with experiment, also address the apparent asymmetry in the jet pattern and show it is fully consistent with momentum conservation.


title: The subcompleteness of diagonal Prikry forcing

abstract: Let $D$ be an infinite discrete set of measurable cardinals. It is shown that generalized Prikry forcing to add a countable sequence to each cardinal in $D$ is subcomplete. To do this it is shown that a simplified version of generalized Prikry forcing which adds a point below each cardinal in $D$, called generalized diagonal Prikry forcing, is subcomplete. Moreover, the generalized diagonal Prikry forcing associated to $D$ is subcomplete above $\mu$, where $\mu$ is any regular cardinal below the first limit point of $D$.


title: Non approximate derivability of the Takagi function

abstract: The Takagi function is a classical example of a continuous nowhere differentiable function. In this paper we prove that it is nowhere approximately derivable.


title: Fast Vessel Segmentation and Tracking in Ultra High-Frequency Ultrasound Images

abstract: Ultra High Frequency Ultrasound (UHFUS) enables the visualization of highly deformable small and medium vessels in the hand. Intricate vessel-based measurements, such as intimal wall thickness and vessel wall compliance, require sub-millimeter vessel tracking between B-scans. Our fast GPU-based approach combines the advantages of local phase analysis, a distance-regularized level set, and an Extended Kalman Filter (EKF), to rapidly segment and track the deforming vessel contour. We validated on 35 UHFUS sequences of vessels in the hand, and we show the transferability of the approach to 5 more diverse datasets acquired by a traditional High Frequency Ultrasound (HFUS) machine. To the best of our knowledge, this is the first algorithm capable of rapidly segmenting and tracking deformable vessel contours in 2D UHFUS images. It is also the fastest and most accurate system for 2D HFUS images.


title: A Study on the Strong Duality of Conic Relaxation of AC Optimal Power Flow in Radial Networks

abstract: This paper proposes a set of closed-form conditions to ensure the strong duality of second-order cone program (SOCP) formulation for AC power flow in radial power networks. In addition, numerical evaluations on IEEE 33-bus test networks and a real-world distribution network are performed to demonstrate the validity of the proposed conditions.


title: Open Gromov-Witten invariants and Boundary States

abstract: In this note we show how Kontsevich-Soibelman algebra arise naturally in Open Gromov-Witten theory for not compact geometries.


title: The large-scale geometry of right-angled Coxeter groups

abstract: This is a survey of some aspects of the large-scale geometry of right-angled Coxeter groups. The emphasis is on recent results on their negative curvature properties, boundaries, and their quasi-isometry and commensurability classification.


title: Mapping class groupoids and Thompson's groups

abstract: We introduce a groupoid ${\mathbf{\Pi MG}}}$, called the fundamental modular groupoid, which is a variant of Penner's mapping class groupoid. We study how it relates to the surface mapping class groups and Thompson's group $\mathsf T$. We also introduce larger groupoid $\mathbf{\Omega MG}$, which is related to outer automorphisms of free groups and Thompson's group $\mathsf V$ in a similar manner.


title: Dynamics of Langton's ant allowed to periodically go straight

abstract: A modified version of Langton's ant is considered. The modified automaton is allowed to go straight $N$-th step instead of turning. The cell state, however, is changed as usually. Depending on the value of $N$ the automaton exhibits different behaviors. Since the Cohen-Kung theorem is not applicable to this modified rule set, in most cases oscillating patterns are observed. For several values of $N$ the automaton leads to a creation of a highway. More interestingly, a few of the automata were found to exhibit a long-term chaotic behavior, exceeding even $10^{13}$ steps. The analysis of the dynamics of the system and emergent patterns is provided.


title: Quantum entanglement and the geometry of spacetime

abstract: This is a brief, popular-level introduction to holographic entanglement. It was published in the newsletter of the International Centre for Theoretical Sciences, Bangalore.


title: Stronger Tests of the Collapse Locality Loophole in Bell Experiments

abstract: Several versions of quantum theory assume some form of localized collapse. If measurement outcomes are indeed defined by localized collapses, then a loophole-free demonstration of Bell non-locality needs to ensure space-like separated collapses associated with the measurements of the entangled systems. This collapse locality loophole remains largely untested, with one significant exception probing Diosi's and Penrose's gravitationally induced collapse hypotheses. I describe here techniques that allow much stronger experimental tests. These apply to all the well known types of collapse postulate, including gravitationally induced collapse, spontaneous localization models and Wigner's consciousness-induced collapse.


title: PCNNA: A Photonic Convolutional Neural Network Accelerator

abstract: Convolutional Neural Networks (CNN) have been the centerpiece of many applications including but not limited to computer vision, speech processing, and Natural Language Processing (NLP). However, the computationally expensive convolution operations impose many challenges to the performance and scalability of CNNs. In parallel, photonic systems, which are traditionally employed for data communication, have enjoyed recent popularity for data processing due to their high bandwidth, low power consumption, and reconfigurability. Here we propose a Photonic Convolutional Neural Network Accelerator (PCNNA) as a proof of concept design to speedup the convolution operation for CNNs. Our design is based on the recently introduced silicon photonic microring weight banks, which use broadcast-and-weight protocol to perform Multiply And Accumulate (MAC) operation and move data through layers of a neural network. Here, we aim to exploit the synergy between the inherent parallelism of photonics in the form of Wavelength Division Multiplexing (WDM) and sparsity of connections between input feature maps and kernels in CNNs. While our full system design offers up to more than 3 orders of magnitude speedup in execution time, its optical core potentially offers more than 5 order of magnitude speedup compared to state-of-the-art electronic counterparts.


title: Exceptional points in the thermoacoustic spectrum

abstract: Exceptional points are found in the spectrum of a prototypical thermoacoustic system as the parameters of the flame transfer function are varied. At these points, two eigenvalues and the associated eigenfunctions coalesce. The system's sensitivity to changes in the parameters becomes infinite. Two eigenvalue branches collide at the exceptional point as the interaction index is increased. One branch originates from a purely acoustic mode, whereas the other branch originates from an intrinsic thermoacoustic mode. The existence of exceptional points in thermoacoustic systems has implications for physical understanding, computing, modeling and control.


title: Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert

abstract: A high-energy neutrino event detected by IceCube on 22 September 2017 was coincident in direction and time with a gamma-ray flare from the blazar TXS 0506+056. Prompted by this association, we investigated 9.5 years of IceCube neutrino observations to search for excess emission at the position of the blazar. We found an excess of high-energy neutrino events with respect to atmospheric backgrounds at that position between September 2014 and March 2015. Allowing for time-variable flux, this constitutes 3.5{\sigma} evidence for neutrino emission from the direction of TXS 0506+056, independent of and prior to the 2017 flaring episode. This suggests that blazars are the first identifiable sources of the high-energy astrophysical neutrino flux.


title: Khovanov homotopy type, periodic links and localizations

abstract: Given an $m$-periodic link $L\subset S^3$, we show that the Khovanov spectrum $\mathcal{X}_L$ constructed by Lipshitz and Sarkar admits a homology group action. We relate the Borel cohomology of $\mathcal{X}_L$ to the equivariant Khovanov homology of $L$ constructed by the second author. The action of Steenrod algebra on the cohomology of $\mathcal{X}_L$ gives an extra structure of the periodic link. Another consequence of our construction is an alternative proof of the localization formula for Khovanov homology, obtained first by Stoffregen and Zhang. By applying Dwyer-Wilkerson theorem we express Khovanov homology of the quotient link in terms of equivariant Khovanov homology of the original link.


title: Ultrahigh-density spin-polarized H and D observed via magnetization quantum beats

abstract: We measure nuclear and electron spin-polarized H and D densities of at least 10$^{19}\, cm^{-3}$ with $\sim$10 ns lifetimes, from the photodissociation of HBr and DI with circularly-polarized UV light pulses. This density is $\sim$6 orders of magnitude higher than that produced by conventional continuous-production methods, and, surprisingly, at least 100 times higher than expected densities for this photodissociation method. We observe the hyperfine quantum beating of the H and D magnetization with a pick-up coil, i.e., the respective 0.7 and 3 ns periodic transfer of polarization from the electrons to the nuclei and back. The $\rm{10^{19}\,cm^{-3}}$ spin-polarized H and D density is sufficient for laser-driven ion acceleration of spin polarized electrons, protons, or deuterons, the preparation of nuclear-spin-polarized molecules, and for the demonstration of spin-polarized D-T or D-$\rm{{^3He}}$ laser fusion, for which a reactivity enhancement of $\rm{\sim50\%}$ is expected.


title: Theta-vexillary signed permutations

abstract: Theta-vexillary signed permutations are elements in the hyperoctahedral group that index certain classes of degeneracy loci of type B and C. These permutations are described using triples of $s$-tuples of integers subject to specific conditions. The objective of this work is to present different characterizations of theta-vexillary signed permutations, describing them in terms of corners in the Rothe diagram and pattern avoidance.


title: Effective Reformulation of Query for Code Search using Crowdsourced Knowledge and Extra-Large Data Analytics

abstract: Software developers frequently issue generic natural language queries for code search while using code search engines (e.g., GitHub native search, Krugle). Such queries often do not lead to any relevant results due to vocabulary mismatch problems. In this paper, we propose a novel technique that automatically identifies relevant and specific API classes from Stack Overflow Q & A site for a programming task written as a natural language query, and then reformulates the query for improved code search. We first collect candidate API classes from Stack Overflow using pseudo-relevance feedback and two term weighting algorithms, and then rank the candidates using Borda count and semantic proximity between query keywords and the API classes. The semantic proximity has been determined by an analysis of 1.3 million questions and answers of Stack Overflow. Experiments using 310 code search queries report that our technique suggests relevant API classes with 48% precision and 58% recall which are 32% and 48% higher respectively than those of the state-of-the-art. Comparisons with two state-of-the-art studies and three popular search engines (e.g., Google, Stack Overflow, and GitHub native search) report that our reformulated queries (1) outperform the queries of the state-of-the-art, and (2) significantly improve the code search results provided by these contemporary search engines.


title: WEIRD: Wide-orbit Exoplanet search with InfraRed Direct imaging

abstract: We report results from the Wide-orbit Exoplanet search with InfraRed Direct imaging (WEIRD), a survey designed to search for Jupiter-like companions on very wide orbits (1000 to 5000 AU) around young stars ($<$120 Myr) that are known members of moving groups in the solar neighborhood ($<$70 pc). Sharing the same age, distance, and metallicity as their host while being on large enough orbits to be studied as "isolated" objects make such companions prime targets for spectroscopic observations and valuable benchmark objects for exoplanet atmosphere models. The search strategy is based on deep imaging in multiple bands across the near-infrared domain. For all 177 objects of our sample, $z_{ab}^\prime$, $J$, [3.6] and [4.5] images were obtained with CFHT/MegaCam, GEMINI/GMOS, CFHT/WIRCam, GEMINI/Flamingos-2, and $Spitzer$/IRAC. Using this set of 4 images per target, we searched for sources with red $z_{ab}^\prime$ and $[3.6]-[4.5]$ colors, typically reaching good completeness down to 2Mjup companions, while going down to 1Mjup for some targets, at separations of $1000-5000$ AU. The search yielded 4 candidate companions with the expected colors, but they were all rejected through follow-up proper motion observations. Our results constrain the occurrence of 1-13 Mjup planetary-mass companions on orbits with a semi-major axis between 1000 and 5000 AU at less than 0.03, with a 95\% confidence level.


title: Comment on "Understanding the $\gamma$-ray emission from the globular cluster 47 Tuc: evidence for dark matter?"

abstract: In a recent paper Brown et al. (2018) analyze the spectral properties of the globular cluster 47 Tucanae (47 Tuc) using 9 years of Fermi-LAT data. Brown et al. (2018) argue that the emission from 47 Tuc cannot be explained by millisecond pulsars (MSPs) alone because of a significant discrepancy between the MSP spectral properties and those of 47 Tuc. It is argued that there is a significant ($>5\sigma$) preference for a two source scenario. The second component could be from the annihilation of dark matter in a density spike surrounding the intermediate-mass black hole candidate in 47 Tuc. In this paper we argue that the claimed discrepancy arises because Brown et al. (2017) use a stacked MSP spectrum to model the emission from MSPs in 47 Tuc which is insufficient to account for the uncertainties in the spectrum of the MSPs in 47 Tuc. Contrary to the claims by Brown et al. (2018), we show that the significance of an additional dark matter component is $\lesssim 2\sigma$ when sample variance in the spectrum of a population of MSPs is accounted for. The spectrum of 47 Tuc is compatible with that of a population of MSPs similar to the disk population.


title: Invariant measures for integrable spin chains and integrable discrete NLS

abstract: We consider discrete analogues of two well-known open problems regarding invariant measures for dispersive PDE, namely, the invariance of the Gibbs measure for the continuum (classical) Heisenberg model and the invariance of white noise under focusing cubic NLS. These continuum models are completely integrable and connected by the Hasimoto transform; correspondingly, we focus our attention on discretizations that are also completely integrable and also connected by a discrete Hasimoto transform. We consider these models on the infinite lattice $\mathbb Z$. Concretely, for a completely integrable variant of the classical Heisenberg spin chain model (introduced independently by Haldane, Ishimori, and Sklyanin) we prove the existence and uniqueness of solutions for initial data following a Gibbs law (which we show is unique) and show that the Gibbs measure is preserved under these dynamics. In the setting of the focusing Ablowitz--Ladik system, we prove invariance of a measure that we will show is the appropriate discrete analogue of white noise. We also include a thorough discussion of the Poisson geometry associated to the discrete Hasimoto transform introduced by Ishimori that connects the two models studied in this article.


title: Landau diamagnetism of the free electron gas as a Fermi surface effect

abstract: The diamagnetic response of the free electron gas called the Landau diamagnetism is a complex and elusive effect requiring laborious computations. Here based on the semi-classical treatment of the problem I present a clear picture of the Landau diamagnetism at zero temperature, which offers a simple derivation of this effect and leads to important consequences: 1) the diamagnetic response is due to electron states in a very narrow Fermi surface region in the k-space, 2) small Fermi energy oscillations in an applied magnetic field are caused by redistribution (inflow or outflow) of electrons from the equatorial region of the Fermi surface. The consideration is based on a structure called magnetic tube whose electron states surround a certain Landau level in k-space. A completely filled magnetic tube does not change its energy in an applied magnetic field as if it complied with the Bohr -- van Leeuwen theorem. The intersection of tubes with the Fermi surface leads to the appearance of partially occupied tubes in the region of intersection. The reconstruction of electron states in a magnetic field in this very small narrow region gives rise to the Landau diamagnetic response. In addition to the Landau diamagnetism this approach fully describes the oscillatory de Haas - van Alphen contribution to the magnetic susceptibility from the equatorial region of the Fermi sphere.


title: Runoff on rooted trees

abstract: We introduce an idealised model for overland flow generated by rain falling on a hill-slope. Our prime motivation is to show how the coalescence of runoff streams promotes the total generation of runoff. We show that, for our model, as the rate of rainfall increases in relation to the soil infiltration rate, there is a distinct phase-change. For low rainfall (the subcritical case) only the bottom of the hill-slope contributes to the total overland runoff, while for high rainfall (the supercritical case) the whole slope contributes and the total runoff increases dramatically. We identify the critical point at which the phase-change occurs, and show how it depends on the degree of coalescence. When there is no stream coalescence the critical point occurs when the rainfall rate equals the average infiltration rate, but when we allow coalescence the critical point occurs when the rainfall rate is less than the average infiltration rate, and increasing the amount of coalescence increases the total expected runoff.


title: GPU-based Commonsense Paradigms Reasoning for Real-Time Query Answering and Multimodal Analysis

abstract: We utilize commonsense knowledge bases to address the problem of real- time multimodal analysis. In particular, we focus on the problem of multimodal sentiment analysis, which consists in the simultaneous analysis of different modali- ties, e.g., speech and video, for emotion and polarity detection. Our approach takes advantages of the massively parallel processing power of modern GPUs to enhance the performance of feature extraction from different modalities. In addition, in order to ex- tract important textual features from multimodal sources we generate domain-specific graphs based on commonsense knowledge and apply GPU-based graph traversal for fast feature detection. Then, powerful ELM classifiers are applied to build the senti- ment analysis model based on the extracted features. We conduct our experiments on the YouTube dataset and achieve an accuracy of 78% which outperforms all previous systems. In term of processing speed, our method shows improvements of several orders of magnitude for feature extraction compared to CPU-based counterparts.


title: Perfectly absorbing exceptional points and chiral absorbers

abstract: We identify a new kind of physically realizable exceptional point (EP) corresponding to degenerate coherent perfect absorption, in which two purely incoming solutions of the wave operator for electromagnetic or acoustic waves coalesce to a single state. Such non-hermitian degeneracies can occur at a real-valued frequency without any associated noise or non-linearity, in contrast to EPs in lasers. The absorption lineshape for the eigenchannel near the EP is quartic in frequency around its maximum in any dimension. In general, for the parameters at which an operator EP occurs, the associated scattering matrix does not have an EP. However, in one dimension, when the $S$-matrix does have a perfectly absorbing EP, it takes on a universal one-parameter form with degenerate values for all scattering coefficients. For absorbing disk resonators, these EPs give rise to chiral absorption: perfect absorption for only one sense of rotation of the input wave.


title: Excitation of higher-order modes in optofluidic photonic crystal fiber

abstract: Higher-order modes up to LP$_{33}$ are controllably excited in water-filled kagom\'{e}- and bandgap-style hollow-core photonic crystal fibers (HC-PCF). A spatial light modulator is used to create amplitude and phase distributions that closely match those of the fiber modes, resulting in typical launch efficiencies of 10-20% into the liquid-filled core. Modes, excited across the visible wavelength range, closely resemble those observed in air-filled kagom\'{e} HC-PCF and match numerical simulations. Mode indices are obtained by launching plane-waves at specific angles onto the fiber input-face and comparing the resulting intensity pattern to that of a particular mode. These results provide a framework for spatially-resolved sensing in HC-PCF microreactors and fiber-based optical manipulation.


title: First passage in an interval for fractional Brownian motion

abstract: Be $X_t$ a random process starting at $x \in [0,1]$ with absorbing boundary conditions at both ends of the interval. Denote $P_1(x)$ the probability to first exit at the upper boundary. For Brownian motion, $P_1(x)=x$, equivalent to $P_1'(x)=1$. For fractional Brownian motion with Hurst exponent $H$, we establish that $P_1'(x) = {\cal N} [x(1-x)]^{\frac1H -2} e^{\epsilon {\cal F}(x)+ {\cal O}(\epsilon^2)}$, where $\epsilon=H-\frac12$. The function ${\cal F}(x)$ is analytic, and well approximated by its Taylor expansion, ${\cal F}(x)\simeq 16 (C-1) (x-1/2)^2 +{\cal O}(x-1/2)^4$, where $C= 0.915...$ is the Catalan-constant. A similar result holds for moments of the exit time starting at $x$. We then consider the span of $X_t$, i.e. the size of the (compact) domain visited up to time $t$. For Brownian motion, we derive an analytic expression for the probability that the span reaches 1 for the first time, then generalized to fBm. Using large-scale numerical simulations with system sizes up to $N=2^{24}$ and a broad range of $H$, we confirm our analytic results. There are important finite-discretization corrections which we quantify. They are most severe for small $H$, necessitating to go to the large systems mentioned above.


title: Time and place of the maximum for one-dimensional diffusion bridges and meanders

abstract: For three constrained Brownian motions, the excursion, the meander, and the reflected bridge, the densities of the maximum and of the time to reach it were expressed as double series by Majumdar, Randon-Furling, Kearney, and Yor (2008). Some of these series were regularized by Abel summation. Similar results for Bessel processes were obtained by Schehr and Le Doussal (2010) using the real space renormalization group method. Here this work is reviewed, and extended from the point of view of one-dimensional diffusion theory to some other diffusion processes including skew Brownian bridges and generalized Bessel meanders. We discuss the limits of the application of this method for other diffusion processes.


title: Observation of the frozen charge of a Kondo resonance

abstract: The ability to control electronic states at the nanoscale has contributed to our modern understanding of condensed matter. In particular, quantum dot circuits represent model systems for the study of strong electronic correlations, epitomized by the Kondo effect. Here, we show that circuit Quantum Electrodynamics architectures can be used to study the internal degrees of freedom of such a many-body phenomenon. We couple a quantum dot to a high finesse microwave cavity to measure with an unprecedented sensitivity the dot electronic compressibility i.e. the ability of the dot to accommodate charges. Because it corresponds solely to the charge response of the electronic system, this quantity is not equivalent to the conductance which involves in general other degrees of freedom such as spin. By performing dual conductance/compressibility measurements in the Kondo regime, we uncover directly the charge dynamics of this peculiar mechanism of electron transfer. Strikingly, the Kondo resonance, visible in transport measurements, is transparent to microwave photons trapped in the high finesse cavity. This reveals that, in such a many body resonance, finite conduction is achieved from a charge frozen by Coulomb interaction. This previously elusive freezing of charge dynamics is in stark contrast with the physics of a free electron gas. Our setup highlights the power of circuit quantum electrodynamics architectures to study condensed matter problems. The tools of cavity quantum electrodynamics could be used in other types of mesoscopic circuits with many-body correlations and bring a promising platform to perform quantum simulation of fermion-boson problems.


title: A Proposal to Detect Dark Matter Using Axionic Topological Antiferromagnets

abstract: Antiferromagnetically doped topological insulators (A-TI) are among the candidates to host dynamical axion fields and axion-polaritons; weakly interacting quasiparticles that are analogous to the dark axion, a long sought after candidate dark matter particle. Here we demonstrate that using the axion quasiparticle antiferromagnetic resonance in A-TI's in conjunction with low-noise methods of detecting THz photons presents a viable route to detect axion dark matter with mass 0.7 to 3.5 meV, a range currently inaccessible to other dark matter detection experiments and proposals. The benefits of this method at high frequency are the tunability of the resonance with applied magnetic field, and the use of A-TI samples with volumes much larger than 1 mm$^3$.

categories: hep-ph cond-mat.mes-hall cond-mat.str-el hep-ex physics.ins-det


title: Operators in Machine Learning: Response Properties in Chemical Space

abstract: The role of response operators is well established in quantum mechanics. We investigate their use for universal quantum machine learning models of response properties in molecules. After introducing a theoretical basis, we present and discuss numerical evidence based on measuring the potential energy's response with respect to atomic displacement and to electric fields. Prediction errors for corresponding properties, atomic forces and dipole moments, improve in a systematic fashion with training set size and reach high accuracy for small training sets. Prediction of normal modes and IR-spectra of some small molecules demonstrates the usefulness of this approach for chemistry.


title: An autonomous robot for continuous tracking of millimetric-sized walkers

abstract: The precise and continuous tracking of millimetric-sized walkers --such as ants-- is quite important in behavioral studies. However, due to technical limitations, most studies concentrate on trajectories within arenas no more than 100 times bigger than the size of the walker or longer trajectories at the expense of either accuracy or continuity. Our work describes a scientific instrument designed to push the boundaries of precise and continuous motion analysis up to 1000 body lengths or more. It consists of a mobile robotic platform that uses Digital Image Processing techniques to track the targets in real time by calculating their spatial position. During the experiments, all the images are stored, and afterwards processed to estimate with higher precision the path traced by the walkers. Some preliminary results achieved using the proposed tracking system are presented.


title: Scaling laws of the Kondo problem at finite frequency

abstract: Driving a quantum system at finite frequency allows one to explore its dynamics. This has become a well mastered resource for controlling the quantum state of two level systems in the context of quantum information processing. However, this can also be of fundamental interest, especially with many-body systems which display an intricate finite frequency behavior. In condensed matter, the Kondo effect epitomizes strong electronic correlations, but the study of its dynamics and the related scaling laws has remained elusive so far. Here, we fill this gap by studying a carbon nanotube based Kondo quantum dot driven by a microwave signal. Our findings not only confirm long-standing theoretical predictions, but also allow us to establish a simple ansatz for the scaling laws on the Kondo problem at finite frequency. More generally, our technique opens a new path for understanding the dynamics of complex quantum dot circuits in the context of quantum simulation of strongly correlated electron fluids.


title: Stable symmetry-protected 3D embedded solitons in Bose-Einstein condensates

abstract: Embedded solitons are rare self-localized nonlinear structures that, counterintuitively, survive inside a continuous background of resonant states. While this topic has been widely studied in nonlinear optics, it has received almost no attention in the field of Bose-Einstein condensation. In this work, we consider experimentally realizable Bose-Einstein condensates loaded in one-dimensional optical lattices and demonstrate that they support continuous families of stable three-dimensional (3D) embedded solitons. These solitons can exist inside the resonant continuous Bloch bands because they are protected by symmetry. The analysis of the Bogoliubov excitation spectrum as well as the long-term evolution after random perturbations proves the robustness of these nonlinear structures against any weak perturbation. This may open up a way for the experimental realization of stable 3D matter-wave embedded solitons as well as for monitoring the gap-soliton to embedded-soliton transition.


title: Effect of Fluoropolymer Composition on Topochemical Synthesis of SrMnO$_{3-{\delta}}$F$_{{\gamma}}$ Oxyfluoride Films

abstract: We report the synthesis of SrMnO$_{3-{\delta}}$F$_{{\gamma}}$ perovskite oxyfluoride thin films using a vapor transport method to fluorinate as-grown SrMnO$_{2.5}$ epitaxial thin films. The influence of the fluoropolymer, which acts as a fluorine vapor source, was investigated by utilizing polyvinyl fluoride (PVF), polyvinylidene difluoride (PVDF) and polytetrafluoroethylene (PTFE) in the reaction. The same process was carried out with polyethylene (PE) to isolate the role of carbon in the vapor transport process. The F distribution was probed by X-ray photoemission spectroscopy, which confirmed the incorporation of F into the films and revealed higher F concentrations in films exposed to PVF and PVDF compared to PTFE. The c-axis parameter expands after fluorination, a result consistent with density functional theory calculations that attribute the volume expansion to elongated Mn-F bonds compared to shorter Mn-O bonds. Using X-ray absorption spectroscopy, we show that the fluorination process reduces the nominal Mn oxidation state suggesting that F substitutes on O sites in the lattice as opposed to filling anion vacancy sites, a finding further supported by calculated formation energies of different F site occupancies. These results provide new insights into topochemical fluorination of perovskite oxides, which should enable future synthesis and design efforts focused on oxyfluoride heterostructures.


title: Multi-messenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A

abstract: Individual astrophysical sources previously detected in neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017 we detected a high-energy neutrino, IceCube-170922A, with an energy of approximately 290 TeV. Its arrival direction was consistent with the location of a known gamma-ray blazar TXS 0506+056, observed to be in a flaring state. An extensive multi-wavelength campaign followed, ranging from radio frequencies to gamma-rays. These observations characterize the variability and energetics of the blazar and include the first detection of TXS 0506+056 in very-high-energy gamma-rays. This observation of a neutrino in spatial coincidence with a gamma-ray emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.


title: Equidistribution on Kuga-Sato Varieties of Torsion Points on CM Elliptic Curves

abstract: A connected Kuga-Sato variety $\mathbf{W}^r$ parameterizes tuples of $r$ points on elliptic curves (with level structure). A special point of $\mathbf{W}^r$ is a tuple of torsion points on a CM elliptic curve. A sequence of special points is strict if any CM elliptic curve appears at most finitely many times and no relation between the points in the tuple is satisfied infinitely often. The genus orbit of a special point is the $\operatorname{Gal}(\bar{\mathbb{Q}}/\mathbb{Q}^{\mathrm{ab}})$-orbit. We show that genus orbits of special points in a strict sequence equidistribute in $\mathbf{W}^r(\mathbb{C})$, assuming a congruence condition at two fixed primes. A genus orbit can be very sparse in the full Galois orbit. In particular, the number of torsion points on each elliptic curve in a genus orbit is not bounded below by the torsion order. A genus orbit corresponds to a toral packet in an extension of $\mathbf{SL}_2$ by a vector representation. These packets also arise in the study by Aka, Einsiedler and Shapira of grids orthogonal to lattice points on the $2$-sphere. As an application we establish their joint equidistribution conjecture assuming two split primes.


title: Alternative zebra-structure models in solar radio emission

abstract: In the literature, discussion continues about the nature of the zebra structure (ZS) in type IV radio bursts, and understanding even the most extended mechanism associated with double plasma resonance has been improved in series of works. Moreover, in the recent work (Ben\'a\v{c}hek, Karlick\'y, Yasnov, 2017) its ineffectiveness was shown under the usually adopted conditions in the radio source. In this case in a number of works we demonstrated the possibility of modeling with whistlers to explain many thin components of ZS stripes, taking into account the effects of scattering whistlers on fast particles. This situation stimulates the search for new mechanisms. For example, earlier we showed the importance of explosive instability, at least for large flares with the ejections of protons. In the system a weakly relativistic beam of protons, nonisothermic plasma, the slow beam mode of the space charge possesses negative energy, and in the triplet slow and fast beam modes and ion acoustic wave an explosive cascade of harmonics from ionic sound is excited. Electromagnetic waves in the form of ZS stripes appear as a result of the fast protons scattering on these harmonics. Such a mechanism can also be promising for a ZS in radio emission from the pulsar in the Crab nebula.


title: The Andromeda galaxy's most important merger about 2 billion years ago as M32's likely progenitor

abstract: Although the proximity of the Andromeda galaxy (M31) offers an opportunity to understand how mergers affect galaxies, uncertainty remains about M31's most important mergers. Previous studies focused individually on the giant stellar stream or the impact of M32 on M31's disk, thereby suggesting many substantial satellite interactions. Yet models of M31's disk heating and the similarity between the stellar populations of different tidal substructures in M31's outskirts both suggested a single large merger. M31's stellar halo (its outer low-surface-brightness regions) is built up from the tidal debris of satellites and provides information about its important mergers. Here we use cosmological models of galaxy formation to show that M31's massive and metal-rich stellar halo, containing intermediate-age stars, dramatically narrows the range of allowed interactions, requiring a single dominant merger with a large galaxy (with stellar mass about 2.5 x 10^10 solar masses, the third largest member of the Local Group) about 2 Gyr ago. This single event explains many observations that were previously considered separately: M31's compact and metal-rich satellite M32 is likely to be the stripped core of the disrupted galaxy, its rotating inner stellar halo contains most of the merger debris, and the giant stellar stream is likely to have been thrown out during the merger. This interaction may explain M31's global burst of star formation about 2 Gyr ago in which approximately a fifth of its stars were formed. Moreover, M31's disk and bulge were already in place, suggesting that mergers of this magnitude need not dramatically affect galaxy structure.


title: RAIM: Recurrent Attentive and Intensive Model of Multimodal Patient Monitoring Data

abstract: With the improvement of medical data capturing, vast amount of continuous patient monitoring data, e.g., electrocardiogram (ECG), real-time vital signs and medications, become available for clinical decision support at intensive care units (ICUs). However, it becomes increasingly challenging to model such data, due to high density of the monitoring data, heterogeneous data types and the requirement for interpretable models. Integration of these high-density monitoring data with the discrete clinical events (including diagnosis, medications, labs) is challenging but potentially rewarding since richness and granularity in such multimodal data increase the possibilities for accurate detection of complex problems and predicting outcomes (e.g., length of stay and mortality). We propose Recurrent Attentive and Intensive Model (RAIM) for jointly analyzing continuous monitoring data and discrete clinical events. RAIM introduces an efficient attention mechanism for continuous monitoring data (e.g., ECG), which is guided by discrete clinical events (e.g, medication usage). We apply RAIM in predicting physiological decompensation and length of stay in those critically ill patients at ICU. With evaluations on MIMIC- III Waveform Database Matched Subset, we obtain an AUC-ROC score of 90.18% for predicting decompensation and an accuracy of 86.82% for forecasting length of stay with our final model, which outperforms our six baseline models.


title: Magnetization curves and low-temperature thermodynamics of two spin-1/2 Heisenberg edge-shared tetrahedra

abstract: A full energy spectrum, magnetization and susceptibility of a spin-1/2 Heisenberg model on two edge-shared tetrahedra are exactly calculated by assuming two different coupling constants. It is shown that a ground state in zero field is either a singlet or a triplet state depending on a relative strength of both coupling constants. Low-temperature magnetization curves may exhibit three different sequences of intermediate plateaux at the following fractional values of the saturation magnetization: 1/3-2/3-1, 0-1/3-2/3-1 or 0-2/3-1. The inverse susceptibility displays a marked temperature dependence significantly influenced by a character of the zero-field ground state. The obtained theoretical results are confronted with recent high-field magnetization data of the mineral crystal fedotovite K2Cu3(SO4)3.


title: Inverse Mean Curvature Flow and the Stability of the Positive Mass Theorem

abstract: We study the stability of the Positive Mass Theorem (PMT) in the case where a sequence of regions of manifolds with positive scalar curvature $U_T^i\subset M_i^3$ are foliated by a smooth solution to Inverse Mean Curvature Flow (IMCF) which may not be uniformly controlled near the boundary. Then if $\partial U_T^i = \Sigma_0^i \cup \Sigma_T^i$, $m_H(\Sigma_T^i) \rightarrow 0$ and extra technical conditions are satisfied we show that $U_T^i$ converges to a flat annulus with respect to Sormani-Wenger Intrinsic Flat (SWIF) convergence.


title: Assessing Test Case Prioritization on Real Faults and Mutants

abstract: Test Case Prioritization (TCP) is an important component of regression testing, allowing for earlier detection of faults or helping to reduce testing time and cost. While several TCP approaches exist in the research literature, a growing number of studies have evaluated them against synthetic software defects, called mutants. Hence, it is currently unclear to what extent TCP performance on mutants would be representative of the performance achieved on real faults. To answer this fundamental question, we conduct the first empirical study comparing the performance of TCP techniques applied to both real-world and mutation faults. The context of our study includes eight well-studied TCP approaches, 35k+ mutation faults, and 357 real-world faults from five Java systems in the Defects4J dataset. Our results indicate that the relative performance of the studied TCP techniques on mutants may not strongly correlate with performance on real faults, depending upon attributes of the subject programs. This suggests that, in certain contexts, the best performing technique on a set of mutants may not be the best technique in practice when applied to real faults. We also illustrate that these correlations vary for mutants generated by different operators depending on whether chosen operators reflect typical faults of a subject program. This highlights the importance, particularly for TCP, of developing mutation operators tailored for specific program domains.


title: A Deterministic Distributed Algorithm for Weighted All Pairs Shortest Paths Through Pipelining

abstract: We present a new pipelined approach to compute all pairs shortest paths (APSP) in a directed graph with nonnegative integer edge weights (including zero weights) in the CONGEST model in the distributed setting. Our deterministic distributed algorithm computes shortest paths of distance at most $\Delta$ for all pairs of vertices in at most $2 n \sqrt{\Delta} + 2n$ rounds, and more generally, it computes h-hop shortest paths for k sources in $2\sqrt{nkh} + n + k$ rounds. The algorithm is simple, and it has some novel features and a nontrivial analysis.It uses only the directed edges in the graph for communication. This algorithm can be used as a base within asymptotically faster algorithms that match or improve on the current best deterministic bound of $\tilde{O}(n^{3/2})$ rounds for this problem when edge weights are $O(n)$ or shortest path distances are $\tilde{O}(n^{3/2})$.


title: Hierarchical Classification using Binary Data

abstract: In classification problems, especially those that categorize data into a large number of classes, the classes often naturally follow a hierarchical structure. That is, some classes are likely to share similar structures and features. Those characteristics can be captured by considering a hierarchical relationship among the class labels. Here, we extend a recent simple classification approach on binary data in order to efficiently classify hierarchical data. In certain settings, specifically, when some classes are significantly easier to identify than others, we showcase computational and accuracy advantages.


title: Out-of-time-ordered correlators in short-range and long-range hard-core boson models and in the Luttinger-liquid model

abstract: We study out-of-time-ordered correlators (OTOC) in hard-core boson models with short-range and long-range hopping and compare the results to the OTOC in the Luttinger liquid model. For the density operator, a related `commutator function' starts at zero and decays back to zero after the passage of the wavefront in all three models, while the wavefront broadens as $t^{1/3}$ in the short-range model and shows no broadening in the long-range model and the Luttinger liquid model. For the boson creation operator, the corresponding commutator function shows saturation inside the light cone in all three models, with similar wavefront behavior as in the density-density commutator function, despite the presence of a nonlocal string in terms of Jordan-Wigner fermions. For the long-range model and the Luttinger liquid model, the commutator function decays as power law outside the light cone in the long time regime when following different fixed-velocity rays. In all cases, the OTOCs approach their long-time values in a power-law fashion, with different exponents for different observables and short-range vs long-range cases. Our long-range model appears to capture exponents in the Luttinger liquid model (which are found to be independent of the Luttinger parameter in the model). This conclusion also bears on the OTOC calculations in conformal field theories, which we propose correspond to long-ranged models.


title: Fundamental groups and group presentations with bounded relator lengths

abstract: We study the geometry of compact geodesic spaces with trivial first Betti number admitting large finite groups of isometries. We show that if a finite group $G$ acts by isometries on a compact geodesic space $X$ whose first Betti number vanishes, then diam$(X) / $diam$(X / G ) \leq 4 \sqrt{ \vert G \vert }$. For a group $G$ and a finite symmetric generating set $S$, $P_k(\Gamma (G, S))$ denotes the 2-dimensional CW-complex whose 1-skeleton is the Cayley graph $\Gamma$ of $G$ with respect to $S$ and whose 2-cells are $m$-gons for $0 \leq m \leq k$, defined by the simple graph loops of length $m$ in $\Gamma$, up to cyclic permutations. Let $G$ be a finite abelian group with $\vert G \vert \geq 3$ and $S$ a symmetric set of generators for which $P_k(\Gamma (G,S))$ has trivial first Betti number. We show that the first nontrivial eigenvalue $-\lambda_1$ of the Laplacian on the Cayley graph satisfies $\lambda_1 \geq 2 - 2 \cos ( 2 \pi / k ) $. We also give an explicit upper bound on the diameter of the Cayley graph of $G$ with respect to $S$ of the form $O (k^2 \vert S \vert \log \vert G \vert )$. Related explicit bounds for the Cheeger constant and Kazhdan constant of the pair $(G,S)$ are also obtained.


title: SOXS: a wide band spectrograph to follow up transients

abstract: SOXS (Son Of X-Shooter) will be a spectrograph for the ESO NTT telescope capable to cover the optical and NIR bands, based on the heritage of the X-Shooter at the ESO-VLT. SOXS will be built and run by an international consortium, carrying out rapid and longer term Target of Opportunity requests on a variety of astronomical objects. SOXS will observe all kind of transient and variable sources from different surveys. These will be a mixture of fast alerts (e.g. gamma-ray bursts, gravitational waves, neutrino events), mid-term alerts (e.g. supernovae, X-ray transients), fixed time events (e.g. close-by passage of minor bodies). While the focus is on transients and variables, still there is a wide range of other astrophysical targets and science topics that will benefit from SOXS. The design foresees a spectrograph with a Resolution-Slit product ~ 4500, capable of simultaneously observing over the entire band the complete spectral range from the U- to the H-band. The limiting magnitude of R~20 (1 hr at S/N~10) is suited to study transients identified from on-going imaging surveys. Light imaging capabilities in the optical band (grizy) are also envisaged to allow for multi-band photometry of the faintest transients. This paper outlines the status of the project, now in Final Design Phase.


title: A First Look at Mobile Intelligence: Architecture, Experimentation and Challenges

abstract: Artificial intelligence (AI) technology makes mobile devices become intelligent objects which can learn and act automatically. Although AI will bring great opportunities for mobile applications, little work has focused on the architecture and the interaction with the cloud. In this article, we present three existing architectures of mobile intelligence in detail and introduce its broad application prospects. Furthermore, we conduct a series of experiments to evaluate the performance of the prevalent commercial applications and intelligent frameworks. Our results show that there is a big gap between Quality of Experience (QoE) requirements and the status quo. So far, we have seen only the tip of the iceberg. We pose issues and challenges to advance the area of mobile intelligence and hope to pave the way for the forthcoming.


title: Estimation of Time Delay caused by Point Geometry in Public Transport

abstract: Travel time prediction is a well-renowned topic of research. It is primarily influenced by traffic congestion, road conditions and route geometry. Among them route geometry at any point is not investigated enough to find a sound spatial dynamics of timing delays and quality of public transport. This study investigates the reliability of travel time to build a new key performance indicator of public transport network. We have introduced a suitable point wise clustering followed by an adapted statistical significance analysis. The outcome is an estimation of absolute delay at geometrical points independent of the delay at bus station. This outcome serves as an incremental delay towards overall delay. Our investigation suggests that this novel metric of delay time contributes around 23.56% in overall delay. The proposed methodology (here in known as Delay Time and Quality of Service (DT QS)) is capable of providing novel insights which attracts attention among policy makers of public transport operators. The novel insight is coined with assignment of high value to level-of-service attributes. Moreover, the proposed technique is independent of any assumption of normal distribution.


title: Detecting Macroscopic Indefiniteness of Cat States in Bosonic Interferometers

abstract: The paradigm of Schr\"{o}dinger's cat illustrates how quantum states preclude the assignment of definite properties to a macroscopic object (realism). In this work we develop a method to investigate the indefiniteness of cat states using currently available cold atom technology. The method we propose uses the observation of a statistical distribution to demonstrate the macroscopic distinction between dead and alive states, and uses the determination of the interferometric sensitivity (Fisher information) to detect the indefiniteness of the cat's vital status. We show how combining the two observations can provide information about the structure of the quantum state without the need for full quantum state tomography, and propose a measure of the indefiniteness based on this structure. We test this method using a cat state proposed by Gordon and Savage [Phys. Rev. A 59, 4623 (1999)] which is dynamically produced from a coherent state. As a control, we consider a set of states produced using the same dynamical procedure acting on an initial thermal distribution. Numerically simulating our proposed method, we show that as the temperature of this initial state is increased, the produced state undergoes a quantum to classical crossover where the indefiniteness of the cat's vital status is lost, while the macroscopic distinction between dead and alive states of the cat is maintained.


title: Ground Layer Adaptive Optics for the W. M. Keck Observatory: Feasibility Study

abstract: Ground-layer adaptive optics (GLAO) systems offer the possibility of improving the "seeing" of large ground-based telescopes and increasing the efficiency and sensitivity of observations over a wide field-of-view. We explore the utility and feasibility of deploying a GLAO system at the W. M. Keck Observatory in order to feed existing and future multi-object spectrographs and wide-field imagers. We also briefly summarize science cases spanning exoplanets to high-redshift galaxy evolution that would benefit from a Keck GLAO system. Initial simulations indicate that a Keck GLAO system would deliver a 1.5x and 2x improvement in FWHM at optical (500 nm) and infrared (1.5 micron), respectively. The infrared instrument, MOSFIRE, is ideally suited for a Keck GLAO feed in that it has excellent image quality and is on the telescope's optical axis. However, it lacks an atmospheric dispersion compensator, which would limit the minimum usable slit size for long-exposure science cases. Similarly, while LRIS and DEIMOS may be able to accept a GLAO feed based on their internal image quality, they lack either an atmospheric dispersion compensator (DEIMOS) or flexure compensation (LRIS) to utilize narrower slits matched to the GLAO image quality. However, some science cases needing shorter exposures may still benefit from Keck GLAO and we will investigate the possibility of installing an ADC.


title: Detecting T-cell receptors involved in immune responses from single repertoire snapshots

abstract: Hypervariable T-cell receptors (TCR) play a key role in adaptive immunity, recognising a vast diversity of pathogen-derived antigens. High throughput sequencing of TCR repertoires (RepSeq) produces huge datasets of T-cell receptor sequences from blood and tissue samples. However, our ability to extract clinically relevant information from RepSeq data is limited, mainly because little is known about TCR-disease associations. Here we present a statistical approach called ALICE (Antigen-specific Lymphocyte Identification by Clustering of Expanded sequences) that identifies TCR sequences that are actively involved in the current immune response from a single RepSeq sample, and apply it to repertoires of patients with a variety of disorders - autoimmune disease (ankylosing spondylitis), patients under cancer immunotherapy, or subject to an acute infection (live yellow fever vaccine). The method's robustness is demonstrated by the agreement of its predictions with independent assays, and is supported by its ability to selectively detect responding TCR in the memory but not in the naive subset. ALICE requires no longitudinal data collection nor large cohorts, and is thus directly applicable to most RepSeq datasets. Its results facilitate the identification of TCR variants associated with a wide variety of diseases and conditions, which can be used for diagnostics, rational vaccine design and evaluation of the adaptive immune system state.


title: Stable Multiple Time Step Simulation/Prediction from Lagged Dynamic Network Regression Models

abstract: Recent developments in computers and automated data collection strategies have greatly increased the interest in statistical modeling of dynamic networks. Many of the statistical models employed for inference on large-scale dynamic networks suffer from limited forward simulation/prediction ability. A major problem with many of the forward simulation procedures is the tendency for the model to become degenerate in only a few time steps, i.e., the simulation/prediction procedure results in either null graphs or complete graphs. Here, we describe an algorithm for simulating a sequence of networks generated from lagged dynamic network regression models DNR(V), a sub-family of TERGMs. We introduce a smoothed estimator for forward prediction based on smoothing of the change statistics obtained for a dynamic network regression model. We focus on the implementation of the algorithm, providing a series of motivating examples with comparisons to dynamic network models from the literature. We find that our algorithm significantly improves multi-step prediction/simulation over standard DNR(V) forecasting. Furthermore, we show that our method performs comparably to existing more complex dynamic network analysis frameworks (SAOM and STERGMs) for small networks over short time periods, and significantly outperforms these approaches over long time time intervals and/or large networks.


title: Valley-dependent Lorentz force and Aharonov-Bohm phase in strained graphene p-n junction

abstract: Veselago lens focusing in graphene p-n junction is promising for realizations of new generation electron optics devices. However, the effect of the strain-induced Aharonov-Bohm interference in a p-n junction has not been discussed before. We provide an experimentally feasible setup based on the Veselago lens in which the presence of strain can result in both the valley-dependent Lorentz force and Aharonov-Bohm interference. In particular, by employing the Green's function and tight binding methods, we study the strain induced by dislocations and line defects in a p-n junction and show how the resulting Aharonov-Bohm phase and interference can be detected. Furthermore, for a different strain configuration, e.g. corresponding to corrugated graphene, we find strong signatures of valley splitting induced by the fictitious magnetic field. Our proposal can be useful for mapping elastic deformations and defects, and for studying valley dependent effects in graphene.


title: A polarization-induced 2D hole gas in undoped gallium nitride quantum wells

abstract: The long-missing polarization-induced two-dimensional hole gas is finally observed in undoped gallium nitride quantum wells. Experimental results provide unambiguous proof that a 2D hole gas in GaN grown on AlN does not need acceptor doping, and can be formed entirely by the difference in the internal polarization fields across the semiconductor heterojunction. The measured 2D hole gas densities, about $4 \times 10^{13}$ cm$^{-2}$, are among the highest among all known semiconductors and remain unchanged down to cryogenic temperatures. Some of the lowest sheet resistances of all wide-bandgap semiconductors are seen. The observed results provide a new probe for studying the valence band structure and transport properties of wide-bandgap nitride interfaces, and simultaneously enable the missing component for gallium nitride-based p-channel transistors for energy-efficient electronics.


title: Invariant multi-graphs in step skew-products

abstract: We study step skew-products over a finite-state shift (base) space whose fiber maps are $C^1$ injective maps on the unit interval. We show that certain invariant sets have a multi-graph structure and can be written graphs of one, two or more functions defined on the base. In particular, this applies to any hyperbolic set and to the support of any ergodic hyperbolic measure. Moreover, within the class of step skew-products whose interval maps are 'absorbing', open and densely the phase space decomposes into attracting and repelling double-strips such that their attractors and repellers are graphs of one single-valued or bi-valued continuous function almost everywhere, respectively.


title: Fisher Information and Logarithmic Sobolev Inequality for Matrix Valued Functions

abstract: We prove a version of Talagrand's concentration inequality for subordinated sub-Laplacian on a compact Riemannian manifold using tools from noncommutative geometry. As an application, motivated by quantum information theory, we show that on a finite dimensional matrix algebra the set of self-adjoint generators satisfying a tensor stable modified logarithmic Sobolev inequality is dense.


title: Enhanced slip properties of lubricant-infused grooves

abstract: We ascertain the enhanced slip properties for a liquid flow over lubricant-infused unidirectional surfaces. This situation reflects many practical settings involving liquid flows past superhydrophobic grooves filled with gas, or past grooves infused with another, immiscible, liquid of smaller or equal viscosity, i.e. where the ratio of lubricant and liquid viscosities, $\mu \leq 1$. To maximize the slippage, we consider deep grooves aligned with the flow. The (normalized by a texture period $L$) effective slip length, $b_{\mathrm{eff}}$, is found as an expansion to first order in protrusion angle $\theta$ about a solution for a flat liquid-lubricant interface. Our results show a significant increase in $b_{\mathrm{eff}}$ with the area fraction of lubricant, $\phi$, and a strong decrease with $\mu$. By contrast, only little influence of $\theta$ on $b_{\mathrm{eff}}$ is observed. Convex meniscus slightly enhances, and concave - slightly reduces $b_{\mathrm{eff}}$ relative the case of a flat liquid-lubricant interface. The largest correction for $\theta$ is found when $\mu = 0$, it decreases with $\mu$, and disappears at $\mu = 1$. Finally, we show that lubricant-infused surfaces of small $\theta$ can be modeled as flat with patterns of local slip boundary conditions, and that the (scaled with $L$) local slip length at the liquid-lubricant interface is an universal function of $\phi$ and $\mu$ only.


title: Ergoregion instability of exotic compact objects: electromagnetic and gravitational perturbations and the role of absorption

abstract: Spinning horizonless compact objects may be unstable against an 'ergoregion instability'. We investigate this mechanism for electromagnetic perturbations of ultracompact Kerr-like objects with a reflecting surface, extending previous (numerical and analytical) work limited to the scalar case. We derive an analytical result for the frequency and the instability time scale of unstable modes which is valid at small frequencies. We argue that our analysis can be directly extended to gravitational perturbations of exotic compact objects in the black-hole limit. The instability for electromagnetic and gravitational perturbations is generically stronger than in the scalar case and it requires larger absorption to be quenched. We argue that exotic compact objects with spin $\chi\lesssim 0.7$ ($\chi\lesssim 0.9$) should have an absorption coefficient of at least $0.3\%$ ($6\%$) to remain linearly stable, and that an absorption coefficient of at least $\approx60\%$ would quench the instability for any spin. We also show that - in the static limit - the scalar, electromagnetic, and gravitatonal perturbations of the Kerr metric are related to one another through Darboux transformations.


title: Characterizing health informatics journals by subject-level dependencies: a citation network analysis

abstract: Citation network analysis has become one of methods to study how scientific knowledge flows from one domain to another. Health informatics is a multidisciplinary field that includes social science, software engineering, behavioral science, medical science and others. In this study, we perform an analysis of citation statistics from health informatics journals using data set extracted from CrossRef. For each health informatics journal, we extract the number of citations from/to studies related to computer science, medicine/clinical medicine and other fields, including the number of self-citations from the health informatics journal. With a similar number of articles used in our analysis, we show that the Journal of the American Medical Informatics Association (JAMIA) has more in-citations than the Journal of Medical Internet Research (JMIR); while JMIR has a higher number of out-citations and self-citations. We also show that JMIR cites more articles from health informatics journals and medicine related journals. In addition, the Journal of Medical Systems (JMS) cites more articles from computer science journals compared with other health informatics journals included in our analysis.


title: Character ratios, representation varieties and random generation of finite groups of Lie type

abstract: We use character theory of finite groups of Lie type to establish new results on representation varieties of Fuchsian groups, and also on probabilistic generation of groups of Lie type.


title: Pressure-induced Superconductivity in Sulfur-doped SnSe Single Crystal Using Boron-doped Diamond Electrode-prefabricated Diamond Anvil Cell

abstract: Sulfur-doped SnSe single crystal was successfully synthesized using a melt and slow-cooling method. The chemical composition and valence state of the obtained sample were analyzed by X-ray photoelectron spectroscopy. The pressure range of a diamond anvil cell with boron-doped diamond electrodes was upgraded to 104 GPa using nano-polycrystalline diamond anvil to investigate a pressure effect for the sample. Electrical resistivity measurements of sulfur-doped SnSe single crystal showed the insulator-metal-superconductor transition by applying high pressure up to 75.9 GPa.


title: Lesion segmentation using U-Net network

abstract: This paper explains the method used in the segmentation challenge (Task 1) in the International Skin Imaging Collaboration's (ISIC) Skin Lesion Analysis Towards Melanoma Detection challenge held in 2018. We have trained a U-Net network to perform the segmentation. The key elements for the training were first to adjust the loss function to incorporate unbalanced proportion of background and second to perform post-processing operation to adjust the contour of the prediction.


title: SOFIA/FORCAST Observations of the Luminous Blue Variable Candidates MN 90 and HD 168625

abstract: We present SOFIA/FORCAST imaging of the circumstellar dust shells surrounding the luminous blue variable (LBV) candidates MN 90 and HD 168625 to quantify the mineral abundances of the dust and to constrain the evolutionary state of these objects. Our image at 37.1 micron of MN 90 shows a limb-brightened, spherical dust shell. A least-squares fit to the spectral energy distribution of MN 90 yields a dust temperature of 59 \pm 10 K, with the peak of the emission at 42.7 micron. Using 2-Dust radiative transfer code, we estimate for MN 90 that mass-loss occurred at a rate of (7.3 \pm 0.4)x10^-7 M_sun/yr x (v_exp/50 km/s) to create a dust shell with a dust mass of (3.2 \pm 0.1)x10^-2 M_sun. Our images between 7.7 - 37.1 micron of HD 168625 complement previously obtained mid-IR imaging of its bipolar nebulae. The SOFIA/FORCAST imaging of HD 168625 shows evidence for the limb-brightened peaks of an equatorial torus. We estimate a dust temperature of 170 \pm 40 K for the equatorial dust surrounding HD 168625, with the peak of the emission at 18.3 micron. Our 2-Dust model for HD 168625 estimates that mass-loss occurred at a rate of (3.2 \pm 0.2)x10^-7 M_sun/yr to create a dust torus/shell with a dust mass of (2.5 \pm 0.1)x10^-3 M_sun.


title: The $g$-good neighbor conditional diagnosability of locally exchanged twisted cubes

abstract: Connectivity and diagnosability are important parameters in measuring the fault tolerance and reliability of interconnection networks. The $R^g$-vertex-connectivity of a connected graph $G$ is the minimum cardinality of a faulty set $X\subseteq V(G)$ such that $G-X$ is disconnected and every fault-free vertex has at least $g$ fault-free neighbors. The $g$-good-neighbor conditional diagnosability is defined as the maximum cardinality of a $g$-good-neighbor conditional faulty set that the system can guarantee to identify. The interconnection network considered here is the locally exchanged twisted cube $LeTQ(s,t)$. For $1\leq s\leq t$ and $0\leq g\leq s$, we first determine the $R^g$-vertex-connectivity of $LeTQ(s,t)$, then establish the $g$-good neighbor conditional diagnosability of $LeTQ(s,t)$ under the PMC model and MM$^*$ model, respectively.


title: Parametric Study of the Rossby Wave Instability in a Two-dimensional Barotropic Disk II: Non-Linear Calculations

abstract: Vortices in protoplanetary disks have attracted attention since the discovery of lopsided structures. One of the possible mechanisms for producing vortices is the Rossby Wave Instability (RWI). In our previous work, we have performed detailed linear stability analyses of the RWI with various initial conditions. In this paper, we perform numerical simulations of the vortex formation by the RWI in 2D barotropic disks using the Athena++ code. As initial conditions, we consider axisymmetric disks with a Gaussian surface density bump of various contrasts and half-widths. Perturbations grow as expected from the linear stability analyses in the linear and weakly non-linear regimes. After the saturation, multiple vortices are formed in accordance with the most unstable azimuthal mode and coalesce one after another. In the end, only one quasi-stationary vortex (the RWI vortex) remains, which migrates inward. During the RWI evolution, the axisymmetric component approaches the stable configuration. We find that the axisymmetric component reaches the marginally stable state for the most unstable azimuthal mode at the saturation and the marginally stable state for the m = 1 mode at the final vortex merger. We investigate the structure and evolution of the RWI vortices. We obtain some empirical relations between the properties of the RWI vortices and the initial conditions. Using tracer particle analyses, we find that the RWI vortex can be considered as a physical entity like a large fluid particle. Our results provide a solid theoretical ground for quantitative interpretation of the observed lopsided structures in protoplanetary disks.


title: Random Sampling and Efficient Algorithms for Multiscale PDEs

abstract: We describe a numerical framework that uses random sampling to efficiently capture low-rank local solution spaces of multiscale PDE problems arising in domain decomposition. In contrast to existing techniques, our method does not rely on detailed analytical understanding of specific multiscale PDEs, in particular, their asymptotic limits. We present the application of the framework on two examples --- a linear kinetic equation and an elliptic equation with rough media. On these two examples, this framework achieves the asymptotic preserving property for the kinetic equations and numerical homogenization for the elliptic equations.


title: Complex self-sustained oscillation patterns in modular excitable networks

abstract: We study the relationship between the modularity of scale-free excitable networks and their ability to support self-sustained oscillation patterns. We find that the probability for a network of given degree-distribution exponent to be able to support self-sustained oscillations is strongly affected by its modularity. In addition, both high- and low-modularity networks are more likely to exhibit long-period oscillation patterns than those with intermediate modularity, but the degrees of complexity and correlation in these two cases are different. The long-period oscillations cannot be explained by a minimum-length Winfree loop, but instead arise from the interplay between two or more propagating waves. Finally, we introduce a new method that can be used to analyze the structure of the self-sustained oscillation sources at different levels of detail and show that the period of the oscillation pattern is statistically correlated with the fraction of modules that are part of the oscillation source.

categories: physics.soc-ph nlin.AO nlin.PS physics.bio-ph


title: Cross-correlations in a quantum dot Cooper pair splitter with ferromagnetic leads

abstract: We investigate Andreev transport through a quantum dot attached to two external ferromagnetic leads and one superconducting electrode. The transport properties of the system are studied by means of the real-time diagrammatic technique in the sequential tunneling regime. To distinguish various contributions to Andreev current we calculate the current cross-correlations, i.e. correlations between currents flowing through two junctions with normal leads. We analyze dependence of current cross-correlations on various parameters of the considered model, both in linear and nonlinear transport regimes. The processes and mechanisms leading to enhancement, suppression or sign change of current cross-correlations are examined and discussed. Interestingly, our results show that for specific transport regimes splitted Cooper pair results in two uncorrelated electrons. However, utilizing ferromagnetic leads instead of non-magnetic electrodes can result in positive current cross-correlations.


title: A Localized Reduced-Order Modeling Approach for PDEs with Bifurcating Solutions

abstract: Reduced-order modeling (ROM) commonly refers to the construction, based on a few solutions (referred to as snapshots) of an expensive discretized partial differential equation (PDE), and the subsequent application of low-dimensional discretizations of partial differential equations (PDEs) that can be used to more efficiently treat problems in control and optimization, uncertainty quantification, and other settings that require multiple approximate PDE solutions. In this work, a ROM is developed and tested for the treatment of nonlinear PDEs whose solutions bifurcate as input parameter values change. In such cases, the parameter domain can be subdivided into subregions, each of which corresponds to a different branch of solutions. Popular ROM approaches such as proper orthogonal decomposition (POD), results in a global low-dimensional basis that does no respect not take advantage of the often large differences in the PDE solutions corresponding to different subregions. Instead, in the new method, the k-means algorithm is used to cluster snapshots so that within cluster snapshots are similar to each other and are dissimilar to those in other clusters. This is followed by the construction of local POD bases, one for each cluster. The method also can detect which cluster a new parameter point belongs to, after which the local basis corresponding to that cluster is used to determine a ROM approximation. Numerical experiments show the effectiveness of the method both for problems for which bifurcation cause continuous and discontinuous changes in the solution of the PDE.


title: Inferring interaction partners from protein sequences using mutual information

abstract: Functional protein-protein interactions are crucial in most cellular processes. They enable multi-protein complexes to assemble and to remain stable, and they allow signal transduction in various pathways. Functional interactions between proteins result in coevolution between the interacting partners, and thus in correlations between their sequences. Pairwise maximum-entropy based models have enabled successful inference of pairs of amino-acid residues that are in contact in the three-dimensional structure of multi-protein complexes, starting from the correlations in the sequence data of known interaction partners. Recently, algorithms inspired by these methods have been developed to identify which proteins are functional interaction partners among the paralogous proteins of two families, starting from sequence data alone. Here, we demonstrate that a slightly higher performance for partner identification can be reached by an approximate maximization of the mutual information between the sequence alignments of the two protein families. Our mutual information-based method also provides signatures of the existence of interactions between protein families. These results stand in contrast with structure prediction of proteins and of multi-protein complexes from sequence data, where pairwise maximum-entropy based global statistical models substantially improve performance compared to mutual information. Our findings entail that the statistical dependences allowing interaction partner prediction from sequence data are not restricted to the residue pairs that are in direct contact at the interface between the partner proteins.


title: Exact solution of some quarter plane walks with interacting boundaries

abstract: The set of random walks with different step sets (of short steps) in the quarter plane has provided a rich set of models that have profoundly different integrability properties. In particular, 23 of the 79 effectively different models can be shown to have generating functions that are algebraic or differentiably finite. Here we investigate how this integrability may change in those 23 models where in addition to length one also counts the number of sites of the walk touching either the horizontal and/or vertical boundaries of the quarter plane. This is equivalent to introducing interactions with those boundaries in a statistical mechanical context. We are able to solve for the generating function in a number of cases. For example, when counting the total number of boundary sites without differentiating whether they are horizontal or vertical, we can solve the generating function of a generalised Kreweras model. However, in many instances we are not able to solve as the kernel methodology seems to break down when including counts with the boundaries.


title: Cohomology of finite tensor categories: duality and Drinfeld centers

abstract: We consider the finite generation property for cohomology of a finite tensor category C, which requires that the self-extension algebra of the unit Ext*_C(1,1) is a finitely generated algebra and that, for each object V in C, the graded extension group Ext*_C(1,V) is a finitely generated module over the aforementioned algebra. We prove that this cohomological finiteness property is preserved under duality (with respect to exact module categories) and taking the Drinfeld center, under suitable restrictions on C. For example, the stated result holds when C is a braided tensor category of odd Frobenius-Perron dimension. By applying our general results, we obtain a number of new examples of finite tensor categories with finitely generated cohomology. In characteristic 0, we show that dynamical quantum groups at roots of unity have finitely generated cohomology. We also provide a new class of examples in finite characteristic which are constructed via infinitesimal group schemes.


title: Weak in the NEES?: Auto-tuning Kalman Filters with Bayesian Optimization

abstract: Kalman filters are routinely used for many data fusion applications including navigation, tracking, and simultaneous localization and mapping problems. However, significant time and effort is frequently required to tune various Kalman filter model parameters, e.g. process noise covariance, pre-whitening filter models for non-white noise, etc. Conventional optimization techniques for tuning can get stuck in poor local minima and can be expensive to implement with real sensor data. To address these issues, a new "black box" Bayesian optimization strategy is developed for automatically tuning Kalman filters. In this approach, performance is characterized by one of two stochastic objective functions: normalized estimation error squared (NEES) when ground truth state models are available, or the normalized innovation error squared (NIS) when only sensor data is available. By intelligently sampling the parameter space to both learn and exploit a nonparametric Gaussian process surrogate function for the NEES/NIS costs, Bayesian optimization can efficiently identify multiple local minima and provide uncertainty quantification on its results.


title: Toward a language-theoretic foundation for planning and filtering

abstract: We address problems underlying the algorithmic question of automating the co-design of robot hardware in tandem with its apposite software. Specifically, we consider the impact that degradations of a robot's sensor and actuation suites may have on the ability of that robot to complete its tasks. We introduce a new formal structure that generalizes and consolidates a variety of well-known structures including many forms of plans, planning problems, and filters, into a single data structure called a procrustean graph, and give these graph structures semantics in terms of ideas based in formal language theory. We describe a collection of operations on procrustean graphs (both semantics-preserving and semantics-mutating), and show how a family of questions about the destructiveness of a change to the robot hardware can be answered by applying these operations. We also highlight the connections between this new approach and existing threads of research, including combinatorial filtering, Erdmann's strategy complexes, and hybrid automata.


title: A Borg-Levinson theorem for magnetic Schr\"odinger operators on a Riemannian manifold

abstract: This article is concerned with uniqueness and stability issues for the inverse spectral problem of recovering the magnetic field and the electric potential in a Riemannian manifold from some asymptotic knowledge of the boundary spectral data of the corresponding Schr\"odinger operator under Dirichlet boundary conditions. The spectral data consist of some asymptotic knowledge of a subset of eigenvalues and Neumann traces of the associated eigenfunctions of the magnetic Laplacian. We also address the same question for Schr\"odinger operators under Neumann boundary conditions, in which case we measure the Dirichlet traces of eigenfunctions. In our results we characterize the uniqueness of the magnetic field from a rate of growth of the eigenvalues, combined with suitable asymptotic properties of boundary observation of eigenfunctions, of the associated magnetic Schr\"odinger operator. To our best knowledge this is the first result proving uniqueness from such general asymptotic behavior of boundary spectral data.


title: Correlated rigidity percolation and colloidal gels

abstract: Rigidity percolation (RP) occurs when mechanical stability emerges in disordered networks as constraints or components are added. Here we discuss RP with structural correlations, an effect ignored in classical theories albeit relevant to many liquid-to-amorphous-solid transitions, such as colloidal gelation, which are due to attractive interactions and aggregation. Using a lattice model, we show that structural correlations shift RP to lower volume fractions. Through molecular dynamics simulations, we show that increasing attraction in colloidal gelation increases structural correlation and thus lowers the RP transition, agreeing with experiments. Hence colloidal gelation can be understood as a RP transition, but occurs at volume fractions far below values predicted by the classical RP, due to attractive interactions which induce structural correlation.


title: Projection based embedding theory for solving Kohn-Sham density functional theory

abstract: Quantum embedding theories are playing an increasingly important role in bridging different levels of approximation to the many body Schr\"odinger equation in physics, chemistry and materials science. In this paper, we present a linear algebra perspective of the recently developed projection based embedding theory (PET) [Manby et al, J. Chem. Theory Comput. 8, 2564, 2012], restricted to the context of Kohn-Sham density functional theory. By partitioning the global degrees of freedom into a `system' part and a `bath' part, and by choosing a proper projector from the bath, PET is an in principle exact formulation to confine the calculation to the system part only, and hence can be carried out with reduced computational cost. Viewed from the perspective of the domain decomposition method, one particularly interesting feature of PET is that it does not enforce a boundary condition explicitly, and remains applicable even when the discretized Hamiltonian matrix is dense, such as in the context of the planewave discretization. In practice, the accuracy of PET depends on the accuracy of the projector for the bath. Based on the linear algebra reformulation, we develop a first order perturbation correction to the projector from the bath to improve its accuracy. Numerical results for real chemical systems indicate that with a proper choice of reference system, the perturbatively corrected PET can be sufficiently accurate even when strong perturbation is applied to very small systems, such as the computation of the ground state energy of a SiH$_3$F molecule, using a SiH$_4$ molecule as the reference system.


title: Evolution of the cold gas fraction and the star formation history: Prospects with current and future radio facilities

abstract: It has recently been shown that the abundance of cold neutral gas may follow a similar evolution as the star formation history. This is physically motivated, since stars form out of this component of the neutral gas and if the case, would resolve the longstanding issue that there is a clear disparity between the total abundance of neutral gas and star forming activity over the history of the Universe. Radio-band 21-cm absorption traces the cold gas and comparison with the Lyman-alpha absorption, which traces all of the gas, provides a measure of the cold gas fraction or the spin temperature. The recent study has shown that the spin temperature (degenerate with the ratio of the absorber/emitter extent) appears to be anti-correlated with the star formation density, undergoing a similar steep evolution as the star formation rate over redshifts of 0 < z < 3, whereas the total neutral hydrogen exhibits little evolution. Above z > 3, where the SFR shows a steep decline with redshift, there is insufficient 21-cm data to determine whether the spin temperature continues to follow the SFR. Knowing this is paramount in ascertaining whether the cold neutral gas does trace the star formation over the Universe's history. We explore the feasibility of resolving this with 21-cm observations of the largest contemporary sample of reliable damped Lyman-alpha absorption systems and conclude that, while today's largest radio interferometers can reach the required sensitivity at z < 3.5, the Square Kilometre Array is required to probe to higher redshifts.


title: Homogeneous Universal H-fields

abstract: We consider derivations $\partial$ on Conway's field $\mathbf{No}$ of surreal numbers such that the ordered differential field $(\mathbf{No},\partial)$ has constant field $\mathbb{R}$ and is a model of the model companion of the theory of $H$-fields with small derivation. We show that this determines $(\mathbf{No},\partial)$ uniquely up to isomorphism, and that this structure is absolutely homogeneous universal for models of this theory with constant field $\mathbb{R}$.


title: Enhanced Light Emission from the Ridge of Two-dimensional InSe Flakes

abstract: InSe, a newly rediscovered two-dimensional (2D) semiconductor, possesses superior electrical and optical properties as a direct bandgap semiconductor with high mobility from bulk to atomically thin layers, drastically different from transition metal dichalcogenides (TMDCs) in which the direct bandgap only exists at the single layer limit. However, absorption in InSe is mostly dominated by an out-of-plane dipole contribution which results in the limited absorption of normally incident light which can only excite the in-plane dipole at resonance. To address this challenge, we have explored a unique geometric ridge state of the 2D flake without compromising the sample quality. We observed the enhanced absorption at the ridge over a broad range of excitation frequencies from photocurrent and photoluminescence (PL) measurements. In addition, we have discovered new PL peaks at low temperature due to defect states on the ridge, which can be as much as ~ 60 times stronger than the intrinsic PL peak of InSe. Interestingly, the PL of the defects is highly tunable through an external electrical field, which can be attributed to the Stark effect of the localized defects. InSe ridges thus provide new avenues for manipulating light-matter interaction and defect-engineering which are vitally crucial for novel optoelectronic devices based on 2D semiconductors.


title: Incremental Embedding: A Density Matrix Embedding Scheme for Molecules

abstract: The idea of using fragment embedding to circumvent the high computational scaling of accurate electronic structure methods while retaining high accuracy has been a long-standing goal for quantum chemists. Traditional fragment embedding methods mainly focus on systems composed of weakly correlated parts and are insufficient when division across chemical bonds is unavoidable. Recently, density matrix embedding theory (DMET) and other methods based on the Schmidt decomposition have emerged as a fresh approach to this problem. Despite their success on model systems, these methods can prove difficult for realistic systems because they rely on either a rigid, non-overlapping partition of the system or a specification of some special sites (i.e. `edge' and `center' sites), neither of which is well-defined in general for real molecules. In this work, we present a new Schmidt decomposition-based embedding scheme called Incremental Embedding that allows the combination of arbitrary overlapping fragments without the knowledge of edge sites. This method forms a convergent hierarchy in the sense that higher accuracy can be obtained by using fragments involving more sites. The computational scaling for the first few levels is lower than that of most correlated wave function methods. We present results for several small molecules in atom-centered Gaussian basis sets and demonstrate that Incremental Embedding converges quickly with fragment size and recovers most static correlation in small basis sets even when truncated at the second lowest level.


title: How Post-Newtonian Dynamics Shape the Distribution of Stationary Binary Black Hole LISA Sources in nearby Globular Clusters

abstract: We derive the observable gravitational wave (GW) peak frequency ($f$) distribution of binary black holes (BBHs) that currently reside inside their globular clusters (GCs), with and without 2.5 Post-Newtonian (2.5PN) effects included in the dynamical evolution of the BBHs. Recent Newtonian studies have reported that a notable number of nearby non-merging BBHs, i.e. those BBHs that are expected to undergo further dynamical interactions before merger, in GCs are likely to be observable by LISA. However, our 2.5PN calculations show that the distribution of $\log f$ for the non-merging BBH population above $\sim 10^{-3.5}$ Hz scales as $f^{-34/9}$ instead of the $f^{-2/3}$ scaling found in the Newtonian case. This leads to an approximately two-orders-of-magnitude reduction in the expected number of GW sources at $\sim 10^{-3}$ Hz, which lead us to conclude that observing nearby BBHs with LISA is not as likely as has been claimed in the recent literature. In fact, our results suggest that it might be more likely that LISA detects the population of BBHs that will merge before undergoing further interactions. This interestingly suggests that the BBH merger rate derived from LIGO can be used to forecast the number of nearby LISA sources, as well as providing insight into the fraction of BBH mergers forming in GCs.