title: A High-Performance Algorithm for Identifying Frequent Items in Data Streams

abstract: Estimating frequencies of items over data streams is a common building block in streaming data measurement and analysis. Misra and Gries introduced their seminal algorithm for the problem in 1982, and the problem has since been revisited many times due its practicality and applicability. We describe a highly optimized version of Misra and Gries' algorithm that is suitable for deployment in industrial settings. Our code is made public via an open source library called DataSketches that is already used by several companies and production systems. Our algorithm improves on two theoretical and practical aspects of prior work. First, it handles weighted updates in amortized constant time, a common requirement in practice. Second, it uses a simple and fast method for merging summaries that asymptotically improves on prior work even for unweighted streams. We describe experiments confirming that our algorithms are more efficient than prior proposals.


title: Role of scalar mesons in the beam asymmetry of $p \bar p$ and $\Lambda \bar \Lambda$ photoproduction at JLab

abstract: We suggest a description of the beam asymmetry in $p \bar p$ and $\Lambda \bar \Lambda$ photoproduction off the proton $\vec\gamma + p \to p \bar p + p$ and $\vec\gamma + p \to \Lambda \bar\Lambda + p$, takes into account the contribution of the scalar mesons $f_0(1370)$, $f_0(1500)$, and $f_0(1710)$. These scalars are considered as mixed states of a glueball and nonstrange and strange quarkonia in the framework based on the use of effective hadronic Lagrangians. Present results can be used to guide the possible search for this reaction by the GlueX Collaboration at JLab. Also, we did an estimate of contribution of heavier scalar meson states $f_0(2020)$, $f_0(2100)$, and $f_0(2200)$.


title: The relationship between undularity and L dependence of the nuclear optical model potential

abstract: The contribution of collective or reaction channels to a local optical model potential, OMP, can be readily calculated as a dynamical polarization potential, DPP. The resulting local DPPs commonly have undulatory (`wavy') features, often including local regions of emissivity in the imaginary component. We show here that this undularity arises from $l$-dependence of the underlying formal non-local and $l$-dependent DPP. The $l$-independent proton OMPs, that have the same $S$-matrix $S_{lj}$ as phenomenological $l$-dependent potentials, exhibit undulations that are qualitatively similar to undulations of local DPPs generated by channel coupling. The $l$-dependent phenomenological potentials studied are the potentials that give the best existing fits to the relevant elastic scattering data and the undulatory potentials presented here, being $S$-matrix equivalent (i.e.\ having the same $S$-matrix, $S_{lj}$) give exactly the same scattering. In addition, we present calculations strongly suggesting that undularity (`waviness') is a generic property of $l$-independent potentials that are $S$-matrix equivalent to $l$-dependent potentials. Implications for the validity of folding models based on a local density model are noted.


title: The Cryogenic AntiCoincidence detector for ATHENA X-IFU: a scientific assessment of the observational capabilities in the hard X-ray band

abstract: ATHENA is a large X-ray observatory, planned to be launched by ESA in 2028 towards an L2 orbit. One of the two instruments of the payload is the X-IFU: a cryogenic spectrometer based on a large array of TES microcalorimeters, able to perform integral field spectrography in the 0.2-12 keV band (2.5 eV FWHM at 6 keV). The X-IFU sensitivity is highly degraded by the particle background expected in the L2 orbit, which is induced by primary protons of both galactic and solar origin, and mostly by secondary electrons. To reduce the particle background level and enable the mission science goals, the instrument incorporates a Cryogenic AntiCoincidence detector (CryoAC). It is a 4 pixel TES based detector, placed <1 mm below the main array. In this paper we report a scientific assessment of the CryoAC observational capabilities in the hard X-ray band (E>10 keV). The aim of the study has been to understand if the present detector design can be improved in order to enlarge the X-IFU scientific capability on an energy band wider than the TES array. This is beyond the CryoAC baseline, being this instrument aimed to operate as anticoincidence particle detector and not conceived to perform X-ray observations.


title: Optimization of edge state velocity in the integer quantum Hall regime

abstract: Observation of interference in the quantum Hall regime may be hampered by a small edge state velocity due to finite phase coherence time. Therefore designing two quantum point contact (QPCs) interferometers having a high edge state velocity is desirable. Here, we present a new simulation method for realistically modeling edge states near QPCs in the integer quantum Hall effect (IQHE) regime. We calculate the filling fraction in the center of the QPC and the velocity of the edge states, and predict structures with high edge state velocity. The 3D Schr\"odinger equation is split into 1D and 2D parts. Quasi-1D Schr\"odinger and Poisson equations are solved self-consistently in the IQHE regime to obtain the potential profile near the edges, and quantum transport is used to solve for the edge state wavefunctions. The velocity of edge states is found to be $\left< E \right> / B$, where $\left< E \right>$ is the expectation value of the electric field for the edge state. Anisotropically etched trench gated heterostructures with double sided delta doping have the highest edge state velocity among the structures considered.


title: Bayesian Nonparametric Poisson-Process Allocation for Time-Sequence Modeling

abstract: Analyzing the underlying structure of multiple time-sequences provides insights into the understanding of social networks and human activities. In this work, we present the \emph{Bayesian nonparametric Poisson process allocation} (BaNPPA), a latent-function model for time-sequences, which automatically infers the number of latent functions. We model the intensity of each sequence as an infinite mixture of latent functions, each of which is obtained using a function drawn from a Gaussian process. We show that a technical challenge for the inference of such mixture models is the unidentifiability of the weights of the latent functions. We propose to cope with the issue by regulating the volume of each latent function within a variational inference algorithm. Our algorithm is computationally efficient and scales well to large data sets. We demonstrate the usefulness of our proposed model through experiments on both synthetic and real-world data sets.


title: Vector-like Leptons: Muon g-2 Anomaly, Lepton Flavor Violation, Higgs Decays, and Lepton Non-Universality

abstract: In this paper, we consider the Standard Model (SM) with one family of vector-like (VL) leptons, which couple to all three families of the SM leptons. We study the constraints on this model coming from the heavy charged lepton mass bound, electroweak precision data, the muon anomalous magnetic moment, lepton flavor violation, Higgs decay constraints and a recently measured lepton non-universality observable, $R_{K^{*0}}$. We find that the strongest constraints are coming from the muon $g-2, \ R_{\mu\mu}=\Gamma(h\to\mu\mu)/\Gamma(h\to\mu\mu)_\text{SM}, \ R_{\gamma\gamma}$ and $\text{BR}(\mu\to e\gamma)$. Although VL leptons couple to all three families of the SM leptons, the ratio of electron-VL to muon-VL coupling is constrained to be $\langle\lambda_e/\lambda_\mu\rangle\lesssim10^{-4}$. We also find that our model cannot fit the observed value of $R_{K^{*0}}$.


title: A Lightweight Regression Method to Infer Psycholinguistic Properties for Brazilian Portuguese

abstract: Psycholinguistic properties of words have been used in various approaches to Natural Language Processing tasks, such as text simplification and readability assessment. Most of these properties are subjective, involving costly and time-consuming surveys to be gathered. Recent approaches use the limited datasets of psycholinguistic properties to extend them automatically to large lexicons. However, some of the resources used by such approaches are not available to most languages. This study presents a method to infer psycholinguistic properties for Brazilian Portuguese (BP) using regressors built with a light set of features usually available for less resourced languages: word length, frequency lists, lexical databases composed of school dictionaries and word embedding models. The correlations between the properties inferred are close to those obtained by related works. The resulting resource contains 26,874 words in BP annotated with concreteness, age of acquisition, imageability and subjective frequency.


title: Bianchi I solutions of the Einstein-Boltzmann system with a positive cosmological constant

abstract: In this paper we study the future global existence and late-time behaviour of the Einstein-Boltzmann system with Bianchi I symmetry and a positive cosmological constant $\Lambda>0$. For the Boltzmann equation we consider the scattering kernel of Israel particles which are the relativistic counterpart of Maxwellian particles. Under a smallness assumption on initial data in a suitable norm we show that solutions exist globally in time and isotropize at late times.


title: Fundamental mode exact schemes for unsteady problems

abstract: The problem of increasing the accuracy of an approximate solution is considered for boundary value problems for parabolic equations. For ordinary differential equations (ODEs), nonstandard finite difference schemes are in common use for this problem. They are based on a modification of standard discretizations of time derivatives and, in some cases, allow to obtain the exact solution of problems. For multidimensional problems, we can consider the problem of increasing the accuracy only for the most important components of the approximate solution. In the present work, new unconditionally stable schemes for parabolic problems are constructed, which are exact for the fundamental mode. Such two-level schemes are designed via a modification of standard schemes with weights using Pad\'{e} approximations. Numerical results obtained for a model problem demonstrate advantages of the proposed fundamental mode exact schemes.


title: Achievable Information Rates for Coded Modulation with Hard Decision Decoding for Coherent Fiber-Optic Systems

abstract: We analyze the achievable information rates (AIRs) for coded modulation schemes with QAM constellations with both bit-wise and symbol-wise decoders, corresponding to the case where a binary code is used in combination with a higher-order modulation using the bit-interleaved coded modulation (BICM) paradigm and to the case where a nonbinary code over a field matched to the constellation size is used, respectively. In particular, we consider hard decision decoding, which is the preferable option for fiber-optic communication systems where decoding complexity is a concern. Recently, Liga \emph{et al.} analyzed the AIRs for bit-wise and symbol-wise decoders considering what the authors called \emph{hard decision decoder} which, however, exploits \emph{soft information} of the transition probabilities of discrete-input discrete-output channel resulting from the hard detection. As such, the complexity of the decoder is essentially the same as the complexity of a soft decision decoder. In this paper, we analyze instead the AIRs for the standard hard decision decoder, commonly used in practice, where the decoding is based on the Hamming distance metric. We show that if standard hard decision decoding is used, bit-wise decoders yield significantly higher AIRs than symbol-wise decoders. As a result, contrary to the conclusion by Liga \emph{et al.}, binary decoders together with the BICM paradigm are preferable for spectrally-efficient fiber-optic systems. We also design binary and nonbinary staircase codes and show that, in agreement with the AIRs, binary codes yield better performance.


title: A note on the metastability in three modifications of the standard Ising model

abstract: We consider three extensions of the standard 2D Ising model with Glauber dynamics on a finite torus at low temperature. The first model is an anisotropic version, where the interaction energy takes different values on vertical and on horizontal bonds. The second model adds next-nearest-neighbor attraction to the standard Ising model. And the third model associates different alternating signs for the magnetic fields on even and odd rows. All these models have already been studied, and results concerning metastability have been established using the so-called pathwise approach. In this text, we extend these earlier results, and apply the potential-theoretic approach to metastability to obtain more precise asymptotic information on the transition time from the metastable phase to the stable phase.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: Limits on quantum deletion from no signaling principle

abstract: One of the fundamental restrictions that quantum mechanics imposes is the "No deletion Theorem" which tells us that given two identical unknown quantum states, it is impossible to delete one of them. But nevertheless if not perfect, people have tried to delete it approximately. In these approximate deleting processes our basic target is to delete one of the two identical copies as much as possible while preserving the other copy. In this brief report, by using the No communication theorem (NCT) (impossibility of sending signal faster than light using a quantum resource) as a guiding principle, we obtain a bound on the sum of the fidelity of deletion and the fidelity of preservation. Our result not only brings out the complementary relation between these two fidelities but also predicts the optimal value of the fidelity of deletion achievable for a given fidelity of preservation under no signaling constraint. This work eventually saturates the quest for finding out the optimal value of deletion within the NCT framework.


title: Typical Relaxation of Isolated Many-Body Systems Which Do Not Thermalize

abstract: We consider isolated many-body quantum systems which do not thermalize, i.e., expectation values approach an (approximately) steady longtime limit which disagrees with the microcanonical prediction of equilibrium statistical mechanics. A general analytical theory is worked out for the typical temporal relaxation behavior in such cases. The main prerequisites are initial conditions which appreciably populate many energy levels and do not give rise to significant spatial inhomogeneities on macroscopic scales. The theory explains very well the experimental and numerical findings in a trapped-ion quantum simulator exhibiting many-body localization, in ultracold atomic gases, and in integrable hard-core boson and XXZ models.


title: Segmentation of 3D High-frequency Ultrasound Images of Human Lymph Nodes Using Graph Cut with Energy Functional Adapted to Local Intensity Distribution

abstract: Previous studies by our group have shown that three-dimensional high-frequency quantitative ultrasound methods have the potential to differentiate metastatic lymph nodes from cancer-free lymph nodes dissected from human cancer patients. To successfully perform these methods inside the lymph node parenchyma, an automatic segmentation method is highly desired to exclude the surrounding thin layer of fat from quantitative ultrasound processing and accurately correct for ultrasound attenuation. In high-frequency ultrasound images of lymph nodes, the intensity distribution of lymph node parenchyma and fat varies spatially because of acoustic attenuation and focusing effects. Thus, the intensity contrast between two object regions (e.g., lymph node parenchyma and fat) is also spatially varying. In our previous work, nested graph cut demonstrated its ability to simultaneously segment lymph node parenchyma, fat, and the outer phosphate-buffered saline bath even when some boundaries are lost because of acoustic attenuation and focusing effects. This paper describes a novel approach called graph cut with locally adaptive energy to further deal with spatially varying distributions of lymph node parenchyma and fat caused by inhomogeneous acoustic attenuation. The proposed method achieved Dice similarity coefficients of 0.937+-0.035 when compared to expert manual segmentation on a representative dataset consisting of 115 three-dimensional lymph node images obtained from colorectal cancer patients.


title: Analytic Few Photon Scattering in Waveguide QED

abstract: We develop an approach to light-matter coupling in waveguide QED based upon scattering amplitudes evaluated via Dyson series. For optical states containing more than single photons, terms in this series become increasingly complex and we provide a diagrammatic recipe for their evaluation, which is capable of yielding analytic results. Our method fully specifies a combined emitter-optical state that permits investigation of light-matter entanglement generation protocols. We use our expressions to study two-photon scattering from a $\Lambda$-system and find that the pole structure of the transition amplitude is dramatically altered as the two ground states are tuned from degeneracy.


title: Boundary amenability of $Out(F_N)$

abstract: We prove that $Out(F_N)$ is boundary amenable. This also holds more generally for $Out(G)$, where $G$ is either a toral relatively hyperbolic group or a finitely generated right-angled Artin group. As a consequence, all these groups satisfy the Novikov conjecture on higher signatures.


title: On Packet Scheduling with Adversarial Jamming and Speedup

abstract: In Packet Scheduling with Adversarial Jamming packets of arbitrary sizes arrive over time to be transmitted over a channel in which instantaneous jamming errors occur at times chosen by the adversary and not known to the algorithm. The transmission taking place at the time of jamming is corrupt, and the algorithm learns this fact immediately. An online algorithm maximizes the total size of packets it successfully transmits and the goal is to develop an algorithm with the lowest possible asymptotic competitive ratio, where the additive constant may depend on packet sizes. Our main contribution is a universal algorithm that works for any speedup and packet sizes and, unlike previous algorithms for the problem, it does not need to know these properties in advance. We show that this algorithm guarantees 1-competitiveness with speedup 4, making it the first known algorithm to maintain 1-competitiveness with a moderate speedup in the general setting of arbitrary packet sizes. We also prove a lower bound of $\phi+1\approx 2.618$ on the speedup of any 1-competitive deterministic algorithm, showing that our algorithm is close to the optimum. Additionally, we formulate a general framework for analyzing our algorithm locally and use it to show upper bounds on its competitive ratio for speedups in $[1,4)$ and for several special cases, recovering some previously known results, each of which had a dedicated proof. In particular, our algorithm is 3-competitive without speedup, matching both the (worst-case) performance of the algorithm by Jurdzinski et al. and the lower bound by Anta et al.


title: Model-Robust Counterfactual Prediction Method

abstract: We develop a novel method for counterfactual analysis based on observational data using prediction intervals for units under different exposures. Unlike methods that target heterogeneous or conditional average treatment effects of an exposure, the proposed approach aims to take into account the irreducible dispersions of counterfactual outcomes so as to quantify the relative impact of different exposures. The prediction intervals are constructed in a distribution-free and model-robust manner based on the conformal prediction approach. The computational obstacles to this approach are circumvented by leveraging properties of a tuning-free method that learns sparse additive predictor models for counterfactual outcomes. The method is illustrated using both real and synthetic data.


title: Controlling supercurrents and their spatial distribution in ferromagnets

abstract: Spin-triplet Cooper pairs induced in ferromagnets form the centrepiece of the emerging field of superconducting spintronics [1,2]. Usually the focus of research is on the spin polarization of the triplets, potentially enabling low-dissipation magnetization switching and domain wall motion. However, the fundamental mechanism for generating triplet pairs [3,4] also permits control over a parameter which has not been addressed before, namely the spatial distribution of the supercurrent. Here we demonstrate this control by tailoring distinct supercurrent pathways in the ferromagnetic weak link of a Josephson junction. Combining micromagnetic simulations with three-dimensional critical current calculations, based on the Usadel description of mesoscopic superconductivity [5], we designed a disk-shaped structure with a magnetic vortex, which induces two distinct supercurrent channels across the junction. The design was successfully tested with superconducting quantum interferometry (SQI). Moreover, we show how the position of the pathways can be controlled by moving the vortex with a magnetic field. This novel approach allows adaptable supercurrent paths to be dynamically reconfigured to switch between different functionalities in the same device.


title: Automorphisms of Toeplitz $\mathscr{B}$-free systems

abstract: For each $\mathscr{B}$-free subshift given by $\mathscr{B}=\{2^ib_i\}_{i\in\mathbb{N}}$, where $\{b_i\}_{i\in\mathbb{N}}$ is a set of pairwise coprime odd numbers greater than one, it is shown that its automorphism group consists solely of powers of the shift.


title: A rigorous derivation of the stationary compressible Reynolds equation via the Navier-Stokes equations

abstract: We provide a rigorous derivation of the compressible Reynolds system as a singular limit of the compressible (barotropic) Navier-Stokes system on a thin domain. In particular, the existence of solutions to the Navier-Stokes system with non-homogeneous boundary conditions is shown that may be of independent interest. Our approach is based on new a priori bounds available for the pressure law of hard sphere type. Finally, uniqueness for the limit problem is established in the 1D case.


title: On a free boundary problem for polymeric fluids: Global existence of weak solutions

abstract: We investigate the stability and global existence of weak solutions to a free boundary problem governing the evolution of polymeric fluids. We construct weak solutions of the two-phase model by performing the asymptotic limit of a macroscopic model governing the suspensions of rod-like molecules (known as Doi-Model) in compressible fluids as the adiabatic exponent $\gamma$ goes to $\infty.$ The convergence of these solutions, up to a subsequence, to the free-boundary problem is established using techniques in the spirit of Lions and Masmoudi \cite{LionsMasmoudi-1999}.


title: Optimal prevention with possibilistic and mixed background risk

abstract: In this paper the effect of posibilistic or mixed background risk on the level of optimal prevention is studied. In the framework of five purely possibilistic or mixed models, necessary and sufficient conditions are found such that the level of optimal saving decreases or increases as a result of the actions of various types of background risk. This way our results complete those obtained by Courbage and Rey for some prevention models with probabilistic background risk.


title: Effective Representations of Clinical Notes

abstract: Clinical notes are a rich source of information about patient state. However, using them to predict clinical events with machine learning models is challenging. They are very high dimensional, sparse and have complex structure. Furthermore, training data is often scarce because it is expensive to obtain reliable labels for many clinical events. These difficulties have traditionally been addressed by manual feature engineering encoding task specific domain knowledge. We explored the use of neural networks and transfer learning to learn representations of clinical notes that are useful for predicting future clinical events of interest, such as all causes mortality, inpatient admissions, and emergency room visits. Our data comprised 2.7 million notes and 115 thousand patients at Stanford Hospital. We used the learned representations, along with commonly used bag of words and topic model representations, as features for predictive models of clinical events. We evaluated the effectiveness of these representations with respect to the performance of the models trained on small datasets. Models using the neural network derived representations performed significantly better than models using the baseline representations with small ($N < 1000$) training datasets. The learned representations offer significant performance gains over commonly used baseline representations for a range of predictive modeling tasks and cohort sizes, offering an effective alternative to task specific feature engineering when plentiful labeled training data is not available.


title: A new complete Calabi-Yau metric on $\mathbb{C}^3$

abstract: Motivated by the study of collapsing Calabi-Yau threefolds with a Lefschetz K3 fibration, we construct a complete Calabi-Yau metric on $\mathbb{C}^3$ with maximal volume growth, which in the appropriate scale is expected to model the collapsing metric near the nodal point. This new Calabi-Yau metric has singular tangent cone at infinity, and its Riemannian geometry has certain non-standard features near the singularity of the tangent cone $\mathbb{C}^2/\mathbb{Z}_2 \times \mathbb{C}$, which are more typical of adiabatic limit problems. The proof uses an existence result in H-J. Hein's PhD thesis to perturb an asymptotic approximate solution into an actual solution, and the main difficulty lies in correcting the slowly decaying error terms.


title: CaloCube: a novel calorimeter for high-energy cosmic rays in space

abstract: In order to extend the direct observation of high-energy cosmic rays up to the PeV region, highly performing calorimeters with large geometrical acceptance and high energy resolution are required. Within the constraint of the total mass of the apparatus, crucial for a space mission, the calorimeters must be optimized with respect to their geometrical acceptance, granularity and absorption depth. CaloCube is a homogeneous calorimeter with cubic geometry, to maximise the acceptance being sensitive to particles from every direction in space; granularity is obtained by relying on small cubic scintillating crystals as active elements. Different scintillating materials have been studied. The crystal sizes and spacing among them have been optimized with respect to the energy resolution. A prototype, based on CsI(Tl) cubic crystals, has been constructed and tested with particle beams. Some results of tests with different beams at CERN are presented.


title: Cryogenic Characterization of FBK RGB-HD SiPMs

abstract: We report on the cryogenic characterization of Red Green Blue - High Density (RGB-HD) SiPMs developed at Fondazione Bruno Kessler (FBK) as part of the DarkSide program of dark matter searches with liquid argon time projection chambers. A dedicated setup was used to measure the primary dark noise, the correlated noise, and the gain of the SiPMs at varying temperatures. A custom-made data acquisition system and analysis software were used to precisely characterize these parameters. We demonstrate that FBK RGB-HD SiPMs with low quenching resistance (RGB-HD-LR$_q$) can be operated from 40 K to 300 K with gains in the range $10^5$ to $10^6$ and noise rates on the order of a few Hz/mm$^2$.


title: Giant Interatomic Energy-Transport Amplification with Nonreciprocal Photonic Topological Insulators

abstract: We show that the energy-transport efficiency in a chain of two-level emitters can be drastically enhanced by the presence of a photonic topological insulator (PTI). This is obtained by exploiting the peculiar properties of its nonreciprocal surface plasmon polariton (SPP), which is unidirectional, and immune to backscattering, and propagates in the bulk band gap. This amplification of transport efficiency can be as much as 2 orders of magnitude with respect to reciprocal SPPs. Moreover, we demonstrate that despite the presence of considerable imperfections at the interface of the PTI, the efficiency of the SPP-assisted energy transport is almost unaffected by discontinuities. We also show that the SPP properties allow energy transport over considerably much larger distances than in the reciprocal case, and we point out a particularly simple way to tune the transport. Finally, we analyze the specific case of a two-emitter chain and unveil the origin of the efficiency amplification. The efficiency amplification and the practical advantages highlighted in this work might be particularly useful in the development of new devices intended to manage energy at the atomic scale.


title: Three Loop Massive Operator Matrix Elements and Asymptotic Wilson Coefficients with Two Different Masses

abstract: Starting at 3-loop order, the massive Wilson coefficients for deep-inelastic scattering and the massive operator matrix elements describing the variable flavor number scheme receive contributions of Feynman diagrams carrying quark lines with two different masses. In the case of the charm and bottom quarks, the usual decoupling of one heavy mass at a time no longer holds, since the ratio of the respective masses, $\eta = m_c^2/m_b^2 \sim 1/10$, is not small enough. Therefore, the usual variable flavor number scheme (VFNS) has to be generalized. The renormalization procedure in the two--mass case is different from the single mass case derived in \cite{Bierenbaum:2009mv}. We present the moments $N=2,4$ and $6$ for all contributing operator matrix elements, expanding in the ratio $\eta$. We calculate the analytic results for general values of the Mellin variable $N$ in the flavor non-singlet case, as well as for transversity and the matrix element $A_{gq}^{(3)}$. We also calculate the two-mass scalar integrals of all topologies contributing to the gluonic operator matrix element $A_{gg}$. As it turns out, the expansion in $\eta$ is usually inapplicable for general values of $N$. We therefore derive the result for general values of the mass ratio. From the single pole terms we derive, now in a two-mass calculation, the corresponding contributions to the 3-loop anomalous dimensions. We introduce a new general class of iterated integrals and study their relations and present special values. The corresponding functions are implemented in computer-algebraic form.


title: Hamilton cycles in infinite cubic graphs

abstract: Investigating a problem of B. Mohar, we show that every one-ended Hamiltonian cubic graph with end degree 3 contains a second Hamilton cycle. We also construct two examples showing that this result does not extend to give a third Hamilton cycle, nor that it extends to the two-ended case.


title: Characterizations of smooth spaces by $\rho_*$-orthogonality

abstract: The aim of this paper is to present some results concerning the $\rho_*$-orthogonality in real normed spaces and its preservation by linear operators. Among other things, we prove that if $T\,: X \longrightarrow Y$ is a nonzero linear $(I, \rho_*)$-orthogonality preserving mapping between real normed spaces, then $$\frac{1}{3}\|T\|\|x\|\leq\|Tx\|\leq 3[T]\|x\|, \qquad (x\in X)$$ where $[T]:=\inf\{\|Tx\|: \,x\in X, \|x\|=1\}$. We also show that the pair $(X,\perp_{\rho_*})$ is an orthogonality space in the sense of R\"{a}tz. Some characterizations of smooth spaces are given based on the $\rho_*$-orthogonality.


title: Maximal monotone operator theory and its applications to thin film equation in epitaxial growth on vicinal surface

abstract: In this work we consider $$ w_t=[(w_{hh}+c_0)^{-3}]_{hh},\qquad w(0)=w^0, $$ which is derived from a thin film equation for epitaxial growth on vicinal surface. We formulate the problem as the gradient flow of a suitably-defined convex functional in a non-reflexive space. Then by restricting it to a Hilbert space and proving the uniqueness of its sub-differential, we can apply the classical maximal monotone operator theory. The mathematical difficulty is due to the fact that $w_{hh}$ can appear as a positive Radon measure. We prove the existence of a global strong solution. In particular, the equation holds almost everywhere when $w_{hh}$ is replaced by its absolutely continuous part.


title: The nonspinning binary black hole merger scenario revisited

abstract: We present the results of 14 simulations of nonspinning black hole binaries with mass ratios $q=m_1/m_2$ in the range $1/100\leq q\leq1$. For each of these simulations we perform three runs at increasing resolution to assess the finite difference errors and to extrapolate the results to infinite resolution. For $q\geq 1/6$, we follow the evolution of the binary typically for the last ten orbits prior to merger. By fitting the results of these simulations, we accurately model the peak luminosity, peak waveform frequency and amplitude, and the recoil of the remnant hole for unequal mass nonspinning binaries. We verify the accuracy of these new models and compare them to previously existing empirical formulas. These new fits provide a basis for a hierarchical approach to produce more accurate remnant formulas in the generic precessing case. They also provide input to gravitational waveform modeling.

categories: gr-qc astro-ph.CO astro-ph.GA astro-ph.HE


title: Quantum spin Hall effect and topological phase transitions in honeycomb antiferromagnets

abstract: While the quantum spin Hall (QSH) effect and antiferromagnetic order constitute two of the most promising phenomena for embedding basic spintronic concepts into future technologies, almost all of the QSH insulators known to date are non-magnetic. Here, based on tight-binding arguments and first-principles theory, we predict two-dimensional antiferromagnets with honeycomb lattice structure to exhibit the QSH effect due to the combined symmetry of time reversal and spatial inversion. We identify functionalized Sn films as experimentally feasible examples which reveal large band gaps rendering these systems ideal for energy efficient spintronics applications. Remarkably, we discover that tensile strain can tune the magnetic order in these materials, accompanied by a topological phase transition from the QSH to the quantum anomalous Hall phase.


title: Gross-Hopkins Duals of Higher Real K-theory Spectra

abstract: We determine the Gross-Hopkins duals of certain higher real $K$-theory spectra. More specifically, let $p$ be an odd prime, and consider the Morava $E$-theory spectrum of height $n=p-1$. It is known, in the expert circles, that for certain finite subgroups $G$ of the Morava stabilizer group, the homotopy fixed point spectra $E_n^{hG}$ are Gross-Hopkins self-dual up to a shift. In this paper, we determine the shift for those finite subgroups $G$ which contain $p$-torsion. This generalizes previous results for $n=2$ and $p=3$.


title: Solutions of the system of operator equations $BXA=B=AXB$ via $*$-order

abstract: In this paper, we establish some necessary and sufficient conditions for the existence of solutions to the system of operator equations $ BXA=B=AXB $ in the setting of bounded linear operators on a Hilbert space, where the unknown operator $X$ is called the inverse of $A$ along $B$. After that, under some mild conditions we prove that an operator $X$ is a solution of $ BXA=B=AXB $ if and only if $B \stackrel{*}{ \leq} AXA$, where the $*$-order $C\stackrel{*}{ \leq} D$ means $CC^*=DC^*, C^*C=C^*D$. Moreover we present the general solution of the equation above. Finally, we present some characterizations of $C \stackrel{*}{ \leq} D$ via other operator equations.


title: The Landscape of Deep Learning Algorithms

abstract: This paper studies the landscape of empirical risk of deep neural networks by theoretically analyzing its convergence behavior to the population risk as well as its stationary points and properties. For an $l$-layer linear neural network, we prove its empirical risk uniformly converges to its population risk at the rate of $\mathcal{O}(r^{2l}\sqrt{d\log(l)}/\sqrt{n})$ with training sample size of $n$, the total weight dimension of $d$ and the magnitude bound $r$ of weight of each layer. We then derive the stability and generalization bounds for the empirical risk based on this result. Besides, we establish the uniform convergence of gradient of the empirical risk to its population counterpart. We prove the one-to-one correspondence of the non-degenerate stationary points between the empirical and population risks with convergence guarantees, which describes the landscape of deep neural networks. In addition, we analyze these properties for deep nonlinear neural networks with sigmoid activation functions. We prove similar results for convergence behavior of their empirical risks as well as the gradients and analyze properties of their non-degenerate stationary points. To our best knowledge, this work is the first one theoretically characterizing landscapes of deep learning algorithms. Besides, our results provide the sample complexity of training a good deep neural network. We also provide theoretical understanding on how the neural network depth $l$, the layer width, the network size $d$ and parameter magnitude determine the neural network landscapes.


title: Geometric aspects of $p$-angular and skew $p$-angular distances

abstract: Corresponding to the concept of $p$-angular distance $\alpha_p[x,y]:=\left\lVert\lVert x\rVert^{p-1}x-\lVert y\rVert^{p-1}y\right\rVert$, we first introduce the notion of skew $p$-angular distance $\beta_p[x,y]:=\left\lVert \lVert y\rVert^{p-1}x-\lVert x\rVert^{p-1}y\right\rVert$ for non-zero elements of $x, y$ in a real normed linear space and study some of significant geometric properties of the $p$-angular and the skew $p$-angular distances. We then give some results comparing two different $p$-angular distances with each other. Finally, we present some characterizations of inner product spaces related to the $p$-angular and the skew $p$-angular distances. In particular, we show that if $p>1$ is a real number, then a real normed space $\mathcal{X}$ is an inner product space, if and only if for any $x,y\in \mathcal{X}\smallsetminus{\lbrace 0\rbrace}$, it holds that $\alpha_p[x,y]\geq\beta_p[x,y]$.


title: Neutrino Mass, Dark Matter and Baryon Asymmetry without Lepton Number Violation

abstract: We propose a model to explain tiny masses of neutrinos with the lepton number conservation, where neither too heavy particles beyond the TeV-scale nor tiny coupling constants are required. Assignments of conserving lepton numbers to new fields result in an unbroken $Z_2$ symmetry that stabilizes the dark matter candidate (the lightest $Z_2$-odd particle). In this model, $Z_2$-odd particles play an important role to generate the mass of neutrinos. The scalar dark matter in our model can satisfy constraints on the dark matter abundance and those from direct searches. It is also shown that the strong first-order phase transition, which is required for the electroweak baryogenesis, can be realized in our model. In addition, the scalar potential can in principle contain CP-violating phases, which can also be utilized for the baryogenesis. Therefore, three problems in the standard model, namely absence of neutrino masses, the dark matter candidate, and the mechanism to generate baryon asymmetry of the Universe, may be simultaneously resolved at the TeV-scale. Phenomenology of this model is also discussed briefly.


title: Posterior sampling for reinforcement learning: worst-case regret bounds

abstract: We present an algorithm based on posterior sampling (aka Thompson sampling) that achieves near-optimal worst-case regret bounds when the underlying Markov Decision Process (MDP) is communicating with a finite, though unknown, diameter. Our main result is a high probability regret upper bound of $\tilde{O}(DS\sqrt{AT})$ for any communicating MDP with $S$ states, $A$ actions and diameter $D$. Here, regret compares the total reward achieved by the algorithm to the total expected reward of an optimal infinite-horizon undiscounted average reward policy, in time horizon $T$. This result closely matches the known lower bound of $\Omega(\sqrt{DSAT})$. Our techniques involve proving some novel results about the anti-concentration of Dirichlet distribution, which may be of independent interest.


title: Relative Entropy and Tsallis Entropy of two Accretive Operators

abstract: Let $A$ and $B$ be two accretive operators. We first introduce the weighted geometric mean of $A$ and $B$ together with some related properties. Afterwards, we define the relative entropy as well as the Tsallis entropy of $A$ and $B$. The present definitions and their related results extend those already introduced in the literature for positive invertible operators.


title: Standardizing densities on Gaussian spaces

abstract: In the present note we investigate the problem of standardizing random variables taking values on infinite dimensional Gaussian spaces. In particular, we focus on the transformations induced on densities by the selected standardization procedure. We discover that, under certain conditions, the Wick exponentials are the key ingredients for treating this kind of problems.


title: Scaling maps of $s$-ordered quasiprobabilities are either nonpositive or completely positive

abstract: Continuous-variable systems in quantum theory can be fully described through any one of the ${\rm s}$-ordered family of quasiprobabilities $\Lambda_{\rm s}(\alpha)$, ${\rm s} \in [-1,1]$. We ask for what values of $({\rm s}, a)$ is the scaling map $\Lambda_{\rm s}(\alpha) \rightarrow a^{-2} \Lambda_{\rm s}(a^{-1}\alpha)$ a positive map? Our analysis based on a duality we establish settles this issue (i) the scaling map generically fails to be positive, showing that there is no useful entanglement witness of the scaling type beyond the transpose map, and (ii) in the two particular cases $({\rm s}=1, |a| \leq 1)$ and $({\rm s}=-1, |a| \geq 1)$, and only in these two non-trivial cases, the map is not only positive but also completely positive as seen through the noiseless attenuator and amplifier channels. We also present a `phase diagram' for the behaviour of the scaling maps in the ${\rm s}-a$ parameter space with regard to its positivity, obtained from the viewpoint of symmetric-ordered characteristic functions. This also sheds light on similar diagrams for the practically relevant attenuation and amplification maps with respect to the noise parameter, especially in the range below the complete-positivity (or quantum-limited) threshold.


title: Doping of Graphene Nanoribbons via Functional Group Edge Modification

abstract: We report on the on-surface synthesis of 7 armchair graphene nanoribbons (7-AGNRs) substituted with nitrile (CN) functional groups. The CN groups are attached to the GNR backbone by modifying the 7-AGNR precursor. While many of these groups survive the on-surface synthesis, the reaction process causes the cleavage of some CN from the ribbon backbone and the on-surface cycloisomerization of few nitriles onto pyridine rings. Scanning Tunneling Spectroscopy and Density Functional Theory reveal that CN groups behave as very efficient n-dopants, significantly downshifting the bands of the ribbon, and introducing deep impurity levels associated to the nitrogen electron lone pairs.


title: Nonlinear energy transfer and current sheet development in localized Alfven wavepacket collisions in the strong turbulence limit

abstract: In space and astrophysical plasmas, turbulence is responsible for transferring energy from large scales driven by violent events or instabilities, to smaller scales where turbulent energy is ultimately converted into plasma heat by dissipative mechanisms. The nonlinear interaction between counterpropagating Alfven waves, denoted Alfven wave collisions, drives this turbulent energy cascade, as recognized by early work with incompressible magnetohydrodynamic (MHD) equations. Recent work employing analytical calculations and nonlinear gyrokinetic simulations of Alfven wave collisions in an idealized periodic initial state have demonstrated the key properties that strong Alfven wave collisions mediate effectively the transfer of energy to smaller perpendicular scales and self-consistently generate current sheets. For the more realistic case of the collision between two initially separated Alfven wavepackets, we use a nonlinear gyrokinetic simulation to show here that these key properties persist: strong Alfven wavepacket collisions indeed facilitate the perpendicular cascade of energy and give rise to current sheets. Furthermore, the evolution shows that nonlinear interactions occur only while the wavepackets overlap, followed by a clean separation of the wavepackets with straight uniform magnetic fields and the cessation of nonlinear evolution in between collisions, even in the gyrokinetic simulation presented here which resolves dispersive and kinetic effects beyond the reach of the MHD theory.


title: Faceted classification: management and use

abstract: The paper discusses issues related to the use of faceted classifications in an online environment. The author argues that knowledge organization systems can be fully utilized in information retrieval only if they are exposed and made available for machine processing. The experience with classification automation to date may be used to speed up and ease the conversion of existing faceted schemes or the creation of management tools for new systems. The author suggests that it is possible to agree on a set of functional requirements for supporting faceted classifications online that are equally relevant for the maintenance of classifications, the creation of classification indexing tools, or the management of classifications in an authority file. It is suggested that a set of requirements for analytico-synthetic classifications may be put forward to improve standards for the use and exchange of knowledge organization systems.


title: Linear regression without correspondence

abstract: This article considers algorithmic and statistical aspects of linear regression when the correspondence between the covariates and the responses is unknown. First, a fully polynomial-time approximation scheme is given for the natural least squares optimization problem in any constant dimension. Next, in an average-case and noise-free setting where the responses exactly correspond to a linear function of i.i.d. draws from a standard multivariate normal distribution, an efficient algorithm based on lattice basis reduction is shown to exactly recover the unknown linear function in arbitrary dimension. Finally, lower bounds on the signal-to-noise ratio are established for approximate recovery of the unknown linear function by any estimator.


title: What are the Receptive, Effective Receptive, and Projective Fields of Neurons in Convolutional Neural Networks?

abstract: In this work, we explain in detail how receptive fields, effective receptive fields, and projective fields of neurons in different layers, convolution or pooling, of a Convolutional Neural Network (CNN) are calculated. While our focus here is on CNNs, the same operations, but in the reverse order, can be used to calculate these quantities for deconvolutional neural networks. These are important concepts, not only for better understanding and analyzing convolutional and deconvolutional networks, but also for optimizing their performance in real-world applications.


title: Thoma type results for discrete quantum groups

abstract: Thoma's theorem states that a group algebra $C^*(\Gamma)$ is of type I if and only if $\Gamma$ is virtually abelian. We discuss here some similar questions for the quantum groups, our main result stating that, under suitable virtually abelianity conditions on a discrete quantum group $\Gamma$, we have a stationary model of type $\pi:C^*(\Gamma)\to M_F(C(L))$, with $F$ being a finite quantum group, and with $L$ being a compact group. We discuss then some refinements of these results in the quantum permutation group case, $\widehat{\Gamma}\subset S_N^+$, by restricting the attention to the matrix models which are quasi-flat, in the sense that the images of the standard coordinates, known to be projections, have rank $\leq1$.


title: Speeding up Memory-based Collaborative Filtering with Landmarks

abstract: Recommender systems play an important role in many scenarios where users are overwhelmed with too many choices to make. In this context, Collaborative Filtering (CF) arises by providing a simple and widely used approach for personalized recommendation. Memory-based CF algorithms mostly rely on similarities between pairs of users or items, which are posteriorly employed in classifiers like k-Nearest Neighbor (kNN) to generalize for unknown ratings. A major issue regarding this approach is to build the similarity matrix. Depending on the dimensionality of the rating matrix, the similarity computations may become computationally intractable. To overcome this issue, we propose to represent users by their distances to preselected users, namely landmarks. This procedure allows to drastically reduce the computational cost associated with the similarity matrix. We evaluated our proposal on two distinct distinguishing databases, and the results showed our method has consistently and considerably outperformed eight CF algorithms (including both memory-based and model-based) in terms of computational performance.


title: Magnetic monopole mass bounds from heavy ion collisions and neutron stars

abstract: Magnetic monopoles, if they exist, would be produced amply in strong magnetic fields and high temperatures via the thermal Schwinger process. Such circumstances arise in heavy ion collisions and in neutron stars, both of which imply lower bounds on the mass of possible magnetic monopoles. In showing this, we construct the cross section for pair production of magnetic monopoles in heavy ion collisions, which indicates that they are particularly promising for experimental searches such as MoEDAL.


title: Macroscopic superpositions require tremendous measurement devices

abstract: We consider fundamental limits on the detectable size of macroscopic quantum superpositions. We argue that a full quantum mechanical treatment of system plus measurement device is required, and that a (classical) reference frame for phase or direction needs to be established to certify the quantum state. When taking the size of such a classical reference frame into account, we show that to reliably distinguish a quantum superposition state from an incoherent mixture requires a measurement device that is quadratically bigger than the superposition state. Whereas for moderate system sizes such as generated in previous experiments this is not a stringent restriction, for macroscopic superpositions of the size of a cat the required effort quickly becomes intractable, requiring measurement devices of the size of the Earth. We illustrate our results using macroscopic superposition states of photons, spins, and position. Finally, we also show how this limitation can be circumvented by dealing with superpositions in relative degrees of freedom.


title: Solar cycle indices from the photosphere to the corona: measurements and underlying physics

abstract: A variety of indices have been proposed in order to represent the many different observables modulated by the solar cycle. Most of these indices are highly correlated with each other owing to their intrinsic link with the solar magnetism and the dominant eleven year cycle, but their variations may differ in fine details, as well as on short- and long-term trends. In this paper we present an overview of the indices that are often employed to describe the many features of the solar cycle, moving from the ones referring to direct observations of the inner solar atmosphere, the photosphere and chromosphere, to those deriving from measurements of the transition region and solar corona. For each index, we summarize existing measurements {\bf and typical use}, and for those that quantify physical observables, we describe the underlying physics.


title: Excitations and relaxation dynamics in multiferroic GeV4S8 studied by THz and dielectric spectroscopy

abstract: We report on THz time-domain spectroscopy on multiferroic GeV4S8, which undergoes orbital ordering at a Jahn-Teller transition at 30.5 K and exhibits antiferromagnetic order below 14.6 K. The THz experiments are complemented by dielectric experiments at audio and radio frequencies. We identify a low-lying excitation close to 15 cm-1, which is only weakly temperature dependent and probably corresponds to a molecular excitation within the electronic level scheme of the V4 clusters. In addition, we detect complex temperature-dependent behavior of a low-lying phononic excitation, closely linked to the onset of orbitally-driven ferroelectricity. In the high-temperature cubic phase, which is paramagnetic and orbitally disordered, this excitation is of relaxational character, becomes an overdamped Lorentzian mode in the orbitally ordered phase below the Jahn-Teller transition, and finally appears as well-defined phonon excitation in the antiferromagnetic state. Abrupt changes in real and imaginary parts of the complex dielectric permittivity show that orbital ordering appears via a structural phase transition with strong first-order character and that the onset of antiferromagnetic order is accompanied by significant structural changes, which are of first-order character, too. Dielectric spectroscopy documents that, at low frequencies, significant dipolar relaxations are present in the orbitally ordered, paramagnetic phase only. In contrast to the closely related GaV4S8, this relaxation dynamics that most likely mirrors coupled orbital and polar fluctuations does not seem to be related to the dynamic processes detected in the THz regime.


title: Thermodynamic signatures of the field-induced states of graphite

abstract: When a magnetic field confines the carriers of a Fermi sea to their lowest Landau level, electron-electron interactions are expected to play a significant role in determining the electronic ground state. Graphite is known to host a sequence of magnetic field-induced states driven by such interactions. Three decades after their discovery, thermodynamic signatures of these instabilities are still elusive. Here, we report the first detection of these transitions with sound velocity. We find that the phase transition occurs in the vicinity of a magnetic field at which at least one of the Landau levels cross the Fermi energy. The evolution of elastic constant anomalies with temperature and magnetic field draws a detailed phase diagram which shows that the ground state evolves in a sequence of thermodynamic phase transitions.Our analysis indicates that electron-electron interaction is not the sole driving force of these transitions and that lattice degrees of freedom play an important role.


title: Masked Autoregressive Flow for Density Estimation

abstract: Autoregressive models are among the best performing neural density estimators. We describe an approach for increasing the flexibility of an autoregressive model, based on modelling the random numbers that the model uses internally when generating data. By constructing a stack of autoregressive models, each modelling the random numbers of the next model in the stack, we obtain a type of normalizing flow suitable for density estimation, which we call Masked Autoregressive Flow. This type of flow is closely related to Inverse Autoregressive Flow and is a generalization of Real NVP. Masked Autoregressive Flow achieves state-of-the-art performance in a range of general-purpose density estimation tasks.


title: Classification revisited: a web of knowledge

abstract: The vision of the Semantic Web (SW) is gradually unfolding and taking shape through a web of linked data, a part of which is built by capturing semantics stored in existing knowledge organization systems (KOS), subject metadata and resource metadata. The content of vast bibliographic collections is currently categorized by some widely used bibliographic classification and we may soon see them being mined for information and linked in a meaningful way across the Web. Bibliographic classifications are designed for knowledge mediation which offers both a rich terminology and different ways in which concepts can be categorized and related to each other in the universe of knowledge. From 1990-2010 they have been used in various resource discovery services on the Web and continue to be used to support information integration in a number of international digital library projects. In this chapter we will revisit some of the ways in which universal classifications, as language independent concept schemes, can assist humans and computers in structuring and presenting information and formulating queries. Most importantly, we highlight issues important to understanding bibliographic classifications, both in terms of their unused potential and technical limitations.


title: Pressure-tuning of bond-directional exchange interactions and magnetic frustration in hyperhoneycomb iridate $\beta$-$\mathrm{Li_2IrO_3}$

abstract: We explore the response of Ir $5d$ orbitals to pressure in $\beta$-$\mathrm{Li_2IrO_3}$, a hyperhoneycomb iridate in proximity to a Kitaev quantum spin liquid (QSL) ground state. X-ray absorption spectroscopy reveals a reconstruction of the electronic ground state below 2 GPa, the same pressure range where x-ray magnetic circular dichroism shows an apparent collapse of magnetic order. The electronic reconstruction, which manifests a reduction in the effective spin-orbit (SO) interaction in $5d$ orbitals, pushes $\beta$-$\mathrm{Li_2IrO_3}$ further away from the pure $J_{\rm eff}=1/2$ limit. Although lattice symmetry is preserved across the electronic transition, x-ray diffraction shows a highly anisotropic compression of the hyperhoneycomb lattice which affects the balance of bond-directional Ir-Ir exchange interactions driven by spin-orbit coupling at Ir sites. An enhancement of symmetric anisotropic exchange over Kitaev and Heisenberg exchange interactions seen in theoretical calculations that use precisely this anisotropic Ir-Ir bond compression provides one possible route to realization of a QSL state in this hyperhoneycomb iridate at high pressures.


title: H4-Alkanes: A new class of hydrogen storage material?

abstract: The methane-based material (H$_2$)$_4$CH$_4$, also called H4M for short, is in essence a methane molecule with 4 physisorbed H$_2$ molecules. While H4M has exceptionally high hydrogen storage densities when it forms a molecular solid, unfortunately, this solid is only stable at impractically high pressures and/or low temperatures. To overcome this limitation, we show through simulations that longer alkanes (methane is the shortest alkane) also form stable structures that still physisorb 4 H$_2$ molecules per carbon atom; we call those structures H4-alkanes. We further show via molecular dynamics simulations that the stability field of molecular solids formed from H4-alkanes increases remarkably with chain length compared to H4M, just as it does for regular alkanes. From our simulations of H4-alkanes with lengths 1, 4, 10, and 20, we see that e.g. for the 20-carbon the stability field is doubled at higher pressures. While even longer chains show only insignificant improvements, we discuss various other options to stabilize H4-alkanes more. Our proof-of-principle results lay the groundwork to show that H4-alkanes can become viable hydrogen storage materials.


title: Black Hole Shadows and Invariant Phase Space Structures

abstract: Utilizing concepts from dynamical systems theory, we demonstrate how the existence of light rings, or fixed points, in a spacetime will give rise to families of periodic orbits and invariant manifolds in phase space. It is shown that these structures define the shape of the black hole shadow as well as a number of salient features of the spacetime lensing. We illustrate this through the analysis of lensing by a hairy black hole.


title: MRI-PET Registration with Automated Algorithm in Pre-clinical Studies

abstract: Magnetic Resonance Imaging (MRI) and Positron Emission Tomography (PET) automatic 3-D registration is implemented and validated for small animal image volumes so that the high-resolution anatomical MRI information can be fused with the low spatial resolution of functional PET information for the localization of lesion that is currently in high demand in the study of tumor of cancer (oncology) and its corresponding preparation of pharmaceutical drugs. Though many registration algorithms are developed and applied on human brain volumes, these methods may not be as efficient on small animal datasets due to lack of intensity information and often the high anisotropy in voxel dimensions. Therefore, a fully automatic registration algorithm which can register not only assumably rigid small animal volumes such as brain but also deformable organs such as kidney, cardiac and chest is developed using a combination of global affine and local B-spline transformation models in which mutual information is used as a similarity criterion. The global affine registration uses a multi-resolution pyramid on image volumes of 3 levels whereas in local B-spline registration, a multi-resolution scheme is applied on the B-spline grid of 2 levels on the finest resolution of the image volumes in which only the transform itself is affected rather than the image volumes. Since mutual information lacks sufficient spatial information, PCA is used to inject it by estimating initial translation and rotation parameters. It is computationally efficient since it is implemented using C++ and ITK library, and is qualitatively and quantitatively shown that this PCA-initialized global registration followed by local registration is in close agreement with expert manual registration and outperforms the one without PCA initialization tested on small animal brain and kidney.


title: Double diffusivity model under stochastic forcing

abstract: The "double diffusivity" model was proposed in the late 1970s, and reworked in the early 1980s, as a continuum counterpart to existing discrete models of diffusion corresponding to high diffusivity paths, such as grain boundaries and dislocation lines. Technically, the model pans out as a system of coupled {\it Fick type} diffusion equations to represent "regular" and "high" diffusivity paths with "source terms" accounting for the mass exchange between the two paths. The model remit was extended by analogy to describe flow in porous media with double porosity, as well as to model heat conduction in media with two non-equilibrium local temperature baths e.g. ion and electron baths. Uncoupling of the two partial differential equations leads to a higher-ordered diffusion equation, solutions of which could be obtained in terms of classical diffusion equation solutions. Similar equations could also be derived within an "internal length" gradient (ILG) mechanics formulation applied to diffusion problems, {\it i.e.}, by introducing nonlocal effects, together with inertia and viscosity, in a mechanics based formulation of diffusion theory. This issue becomes particularly important in the case of diffusion in nanopolycrystals whose deterministic ILG based theoretical calculations predict a relaxation time that is only about one-tenth of the actual experimentally verified timescale. This article provides the "missing link" in this estimation by adding a vital element in the ILG structure, that of stochasticity, that takes into account all boundary layer fluctuations. Our stochastic-ILG diffusion calculation confirms rapprochement between theory and experiment, thereby benchmarking a new generation of gradient-based continuum models that conform closer to real life fluctuating environments.


title: Improving two - qubit state teleportation affected by amplitude damping noise based on choosing appropriate quantum channel

abstract: We consider two qubit teleportation via quantum channel affected by amplitude damping noise. Addressing the same problem, X. Hu, Y. Gu, Q. Gong and G. Guo [Phys. Rev. A 81, 054302, (2010)] recently showed that in presence of noise, subjecting more qubits in quantum channel to amplitude damping can increase the fidelity of teleportation protocol. However, in this paper, by making some adjustments on quantum channel, we obtain teleportation fidelity which is even higher than one in the case of X. Hu et al. Moreover, our strategy is simpler than quantum distillation and compared to using weak measurement, it is deterministic. Furthermore, explicit analysis of fidelity is provided, we show that in general, choosing appropriate quantum channel enhances the ability of teleportation better and negates the fact that more amplitude damping noise more quality.


title: Analyzing Privacy Breaches in the Aircraft Communications Addressing and Reporting System (ACARS)

abstract: The manner in which Aircraft Communications, Addressing and Reporting System (ACARS) is being used has significantly changed over time. Whilst originally used by commercial airliners to track their flights and provide automated timekeeping on crew, today it serves as a multi-purpose air-ground data link for many aviation stakeholders including private jet owners, state actors and military. Since ACARS messages are still mostly sent in the clear over a wireless channel, any sensitive information sent with ACARS can potentially lead to a privacy breach for users. Naturally, different stakeholders consider different types of data sensitive. In this paper we propose a privacy framework matching aviation stakeholders to a range of sensitive information types and assess the impact for each. Based on more than one million ACARS messages, collected over several months, we then demonstrate that current ACARS usage systematically breaches privacy for all stakeholder groups. We further support our findings with a number of cases of significant privacy issues for each group and analyze the impact of such leaks. While it is well-known that ACARS messages are susceptible to eavesdropping attacks, this work is the first to quantify the extent and impact of privacy leakage in the real world for the relevant aviation stakeholders.


title: Cycle decompositions of pathwidth-6 graphs

abstract: Haj\'os conjecture asserts that a simple Eulerian graph on n vertices can be decomposed into at most (n - 1)/2 cycles. The conjecture is only proved for graph classes in which every element contains vertices of degree 2 or 4. We develop new techniques to construct cycle decompositions. They work on the common neighbourhood of two degree-6 vertices. With these techniques we find structures that cannot occur in a minimal counterexample to Haj\'os conjecture and verify the conjecture for Eulerian graphs of pathwidth at most 6. This implies that these graphs satisfy the small cycle double cover conjecture.


title: Electrically tunable valley polarization in Weyl semimetals with tilted energy dispersion

abstract: Tunneling transport across the p-n-p junction of Weyl semimetal with tilted energy dispersion is investigated. We report that the electrons around different valleys experience opposite direction refractions at the barrier interface when the energy dispersion is tilted along one of the transverse directions. Chirality dependent refractions at the barrier interface polarize the Weyl fermions in angle-space according to their valley index. A real magnetic barrier configuration is used to select allowed transmission angles, which results in electrically controllable and switchable valley polarization. Our findings may pave the way for experimental investigation of valley polarization, as well as valleytronic and electron optic applications in Weyl semimetals.


title: From refined estimates for spherical harmonics to a sharp multiplier theorem on the Grushin sphere

abstract: We prove a sharp multiplier theorem of Mihlin-H\"ormander type for the Grushin operator on the unit sphere in $\mathbb{R}^3$, and a corresponding boundedness result for the associated Bochner-Riesz means. The proof hinges on precise pointwise bounds for spherical harmonics.


title: CacheShuffle: An Oblivious Shuffle Algorithm Using Caches

abstract: We consider Oblivious Shuffling and K-Oblivious Shuffling, a refinement thereof. We provide efficient algorithms for both and discuss their application to the design of Oblivious RAM. The task of K-Oblivious Shuffling is to obliviously shuffle N encrypted blocks that have been randomly allocated on the server in such a way that an adversary learns nothing about the new allocation of blocks. The security guarantee should hold also with respect to an adversary that has learned the initial position of K touched blocks out of the N blocks. The classical notion of Oblivious Shuffling is obtained for K = N. We present a family of algorithms for Oblivious Shuffling. Our first construction, CacheShuffleRoot, is tailored for clients with $O(\sqrt{N})$ blocks of memory and uses $(4+\epsilon)N$ blocks of bandwidth, for every $\epsilon > 0$. CacheShuffleRoot is a 4.5x improvement over previous best known results on practical sizes of N. We also present CacheShuffle that obliviously shuffles using O(S) blocks of client memory with $O(N\log_S N)$ blocks of bandwidth. We then turn to K-Oblivious Shuffling and give algorithms that require 2N + f(K) blocks of bandwidth, for some function f. That is, any extra bandwidth above the 2N lower bound depends solely on K. We present KCacheShuffleBasic that uses O(K) client storage and exactly 2N blocks of bandwidth. For smaller client storage requirements, we show KCacheShuffle, which uses O(S) client storage and requires $2N+(1+\epsilon)O(K\log_S K)$ blocks of bandwidth. Finally, we consider the case in which, in addition to the N blocks, the server stores D dummy blocks whose content is is irrelevant but still their positions must be hidden by the shuffling. For this case, we design algorithm KCacheShuffleDummy that, for N + D blocks and K touched blocks, uses O(K) client storage and $D+(2+\epsilon)N$ blocks of bandwidth.


title: EE-Grad: Exploration and Exploitation for Cost-Efficient Mini-Batch SGD

abstract: We present a generic framework for trading off fidelity and cost in computing stochastic gradients when the costs of acquiring stochastic gradients of different quality are not known a priori. We consider a mini-batch oracle that distributes a limited query budget over a number of stochastic gradients and aggregates them to estimate the true gradient. Since the optimal mini-batch size depends on the unknown cost-fidelity function, we propose an algorithm, {\it EE-Grad}, that sequentially explores the performance of mini-batch oracles and exploits the accumulated knowledge to estimate the one achieving the best performance in terms of cost-efficiency. We provide performance guarantees for EE-Grad with respect to the optimal mini-batch oracle, and illustrate these results in the case of strongly convex objectives. We also provide a simple numerical example that corroborates our theoretical findings.


title: Sub-picosecond spin dynamics of excited states in the topological insulator Bi$_2$Te$_3$

abstract: Using time-, spin- and angle-resolved photoemission, we investigate the ultrafast spin dynamics of hot electrons on the surface of the topological insulator Bi$_2$Te$_3$ following optical excitation by fs-infrared pulses. We observe two surface-resonance states above the Fermi level coexisting with a transient population of Dirac fermions that relax in about $\sim$2 ps. One state is located below $\sim$0.4 eV just above the bulk continuum, the other one at $\sim$0.8 eV inside a projected bulk band gap. At the onset of the excitation, both states exhibit a reversed spin texture with respect to that of the transient Dirac bands, in agreement with our one-step photoemission calculations. Our data reveal that the high-energy state undergoes spin relaxation within $\sim$0.5 ps, a process that triggers the subsequent spin dynamics of both the Dirac cone and the low-energy state, which behave as two dynamically-locked electron populations. We discuss the origin of this behavior by comparing the relaxation times observed for electrons with opposite spins to the ones obtained from a microscopic Boltzmann model of ultrafast band cooling introduced into the photoemission calculations. Our results demonstrate that the nonequilibrium surface dynamics is governed by electron-electron rather than electron-phonon scattering, with a characteristic time scale unambiguously determined by the complex spin texture of excited states above the Fermi level. Our findings reveal the critical importance of detecting momentum and energy-resolved spin textures with fs resolution to fully understand the sub-ps dynamics of transient electrons on the surface of topological insulators.


title: Lecture Notes on the Statistical Mechanics of Disordered Systems

abstract: This material complements David Chandler's Introduction to Modern Statistical Mechanics (Oxford University Press, 1987) in a graduate-level, one-semester course I teach in the Department of Chemistry at Duke University. Students enter this course with some knowledge of statistical thermodynamics and quantum mechanics, usually acquired from undergraduate physical chemistry at the level of D. A. McQuarrie & J. D. Simon's Physical Chemistry: A Molecular Approach (University Science Books, 1997). These notes, which introduce students to a modern treatment of glassiness and to the replica method, build on the material and problems contained in the eight chapters of Chandler's textbook.


title: Unified Fock space representation of fractional quantum Hall states

abstract: Many bosonic (fermionic) fractional quantum Hall states, such as Laughlin, Moore-Read and Read-Rezayi wavefunctions, belong to a special class of orthogonal polynomials: the Jack polynomials (times a Vandermonde determinant). This fundamental observation allows to point out two different recurrence relations for the coefficients of the permanent (Slater) decomposition of the bosonic (fermionic) states. Here we provide an explicit Fock space representation for these wavefunctions by introducing a two-body squeezing operator which represents them as a Jastrow operator applied to reference states, which are in general simple periodic one dimensional patterns. Remarkably, this operator representation is the same for bosons and fermions, and the different nature of the two recurrence relations is an outcome of particle statistics.


title: On the number of faces of Gelfand-Zetlin polytopes

abstract: We obtain a recurrence relation for the f-polynomial of Gelfand-Zetlin polytopes by analyzing geometric properties of a linear projection of the Gelfand-Zetlin polytope onto a cube. We apply this recurrence relation to find explicit formulas for the f-polynomials and h-polynomials of Gelfand-Zetlin polytopes that belong to several simplest one-parameter families.


title: Ultra-brilliance isolated attosecond gamma-ray light source from nonlinear Compton scattering

abstract: The explosion in attosecond technology has opened the gate to investigating many unexplored areas which require ultrahigh spatial and temporal resolution. In the area of nuclear physics, using gamma-rays with ultrahigh resolution in time and space will help to investigate intra-nuclear dynamics in an unprecedentedly explicit way. However, the generation of ultrahigh brilliance attosecond gamma-ray pulses with current-generation laser facilities has not been reported. In this letter, we propose a novel method to generate high charge (~1nC) attosecond (<200 attosecond) electron bunch by the near-threshold self-injection in a wakefield accelerator. We demonstrate the ability to generate an ultrahigh-brilliance (> 2*1024 photons s-1mm-2mrad-2 per 0.1%BW) attosecond (<200 attosecond) gamma-ray (Emax > 3 MeV) pulse via nonlinear Compton scattering. To the best of our knowledge, this is the first method reported to generate attosecond gamma-ray photon source using current-generation laser. This is the shortest gamma-ray photon and the highest brilliance photon source in MeV range (orders higher than the results reported). This method can be widely applied for experimental generation of 100 keV to several MeV high brilliance attosecond gamma-ray sources with current ~100 TW laser facilities, which will benefit basic science such as application in ultra-high resolution radiography.


title: Ferromagnetic Dirac-Nodal Semimetal Phase of Stretched Chromium Dioxide

abstract: We show by first-principles calculations that the Dirac nodal-line semimetal phase can co-exist with the ferromagnetic order at room temperature in chromium dioxide, a widely used material in magnetic tape applications, under small tensile hydrostatic strains. An ideally flat Dirac nodal ring close to the Fermi energy is placed in the reflection-invariant boundary of the Brillouin zone perpendicular to the magnetic order, and is topologically protected by the unitary mirror symmetry of the magnetic group $D_{4h}(C_{4h})$, which quantizes the corresponding Berry phase into integer multiples of $\pi$. The symmetry-dependent topological stability is demonstrated through showing that only the topologically protected nodal ring can persistently exist under small anisotropic stains preservingthe symmetry $D_{4h}(C_{4h})$, while the other seeming band touching points are generically gapped. Our work provides a practical platform for the investigation of novel physics and potential applications of the Dirac nodal-line and drumhead fermions, in particular those related to ferromagnetic properties


title: What do We Learn by Semantic Scene Understanding for Remote Sensing imagery in CNN framework?

abstract: Recently, deep convolutional neural network (DCNN) achieved increasingly remarkable success and rapidly developed in the field of natural image recognition. Compared with the natural image, the scale of remote sensing image is larger and the scene and the object it represents are more macroscopic. This study inquires whether remote sensing scene and natural scene recognitions differ and raises the following questions: What are the key factors in remote sensing scene recognition? Is the DCNN recognition mechanism centered on object recognition still applicable to the scenarios of remote sensing scene understanding? We performed several experiments to explore the influence of the DCNN structure and the scale of remote sensing scene understanding from the perspective of scene complexity. Our experiment shows that understanding a complex scene depends on an in-depth network and multiple-scale perception. Using a visualization method, we qualitatively and quantitatively analyze the recognition mechanism in a complex remote sensing scene and demonstrate the importance of multi-objective joint semantic support.


title: Laser-induced persistent photovoltage on the surface of a ternary topological insulator at room temperature

abstract: Using time- and angle-resolved photoemission, we investigate the ultrafast response of excited electrons in the ternary topological insulator (Bi$_{1 x}$Sb$_{x}$)$_2$Te$_3$ to fs-infrared pulses. We demonstrate that at the critical concentration $x$=0.55, where the system becomes bulk insulating, a surface voltage can be driven at room temperature through the topological surface state solely by optical means. We further show that such a photovoltage persists over a time scale that exceeds $\sim$6 $\mu$s, i.e, much longer than the characteristic relaxation times of bulk states. We attribute the origin of the photovoltage to a laser-induced band-bending effect which emerges near the surface region on ultrafast time scales. The photovoltage is also accompanied by a remarkable increase in the relaxation times of excited states as compared to undoped topological insulators. Our findings are relevant in the context of applications of topological surface states in future optical devices.


title: Scalable Variational Inference for Dynamical Systems

abstract: Gradient matching is a promising tool for learning parameters and state dynamics of ordinary differential equations. It is a grid free inference approach, which, for fully observable systems is at times competitive with numerical integration. However, for many real-world applications, only sparse observations are available or even unobserved variables are included in the model description. In these cases most gradient matching methods are difficult to apply or simply do not provide satisfactory results. That is why, despite the high computational cost, numerical integration is still the gold standard in many applications. Using an existing gradient matching approach, we propose a scalable variational inference framework which can infer states and parameters simultaneously, offers computational speedups, improved accuracy and works well even under model misspecifications in a partially observable system.


title: Bitwise Operations of Cellular Automaton on Gray-scale Images

abstract: Cellular Automata (CA) theory is a discrete model that represents the state of each of its cells from a finite set of possible values which evolve in time according to a pre-defined set of transition rules. CA have been applied to a number of image processing tasks such as Convex Hull Detection, Image Denoising etc. but mostly under the limitation of restricting the input to binary images. In general, a gray-scale image may be converted to a number of different binary images which are finally recombined after CA operations on each of them individually. We have developed a multinomial regression based weighed summation method to recombine binary images for better performance of CA based Image Processing algorithms. The recombination algorithm is tested for the specific case of denoising Salt and Pepper Noise to test against standard benchmark algorithms such as the Median Filter for various images and noise levels. The results indicate several interesting invariances in the application of the CA, such as the particular noise realization and the choice of sub-sampling of pixels to determine recombination weights. Additionally, it appears that simpler algorithms for weight optimization which seek local minima work as effectively as those that seek global minima such as Simulated Annealing.


title: Secure Computation of Randomized Functions: Further Results

abstract: We consider secure computation of randomized functions between two users, where both the users (Alice and Bob) have inputs, Alice sends a message to Bob over a rate-limited, noise-free link, and then Bob produces the output. We study two cases: (i) when privacy condition is required only against Bob, who tries to learn more about Alice's input from the message than what can be inferred by his own input and output, and (ii) when there is no privacy requirement. For both the problems, we give single-letter expressions for the optimal rates. For the first problem, we also explicitly characterize securely computable randomized functions when input has full support, which leads to a much simpler expression for the optimal rate. Recently, Data (ISIT 2016) studied the remaining two cases (first, when privacy conditions are against both the users; and second, when privacy condition is only against Alice) and obtained single-letter expressions for optimal rates in both the scenarios.


title: Gravity and handedness of photons

abstract: Vacuum fluctuations of quantum fields are altered in presence of a strong gravitational background, with important physical consequences. We argue that a non-trivial spacetime geometry can act as an optically active medium for quantum electromagnetic radiation, in such a way that the state of polarization of radiation changes in time, even in the absence of electromagnetic sources. This is a quantum effect, and is a consequence of an anomaly related to the classical invariance under electric-magnetic duality rotations in Maxwell theory.


title: Generalized bilinear forms graphs and MDR codes

abstract: We investigate the generalized bilinear forms graph $\Gamma_d$ over a residue class ring $\mathbb{Z}_{p^s}$. We show that $\Gamma_d$ is a connected vertex transitive graph, and completely determine its independence number, clique number, chromatic number and maximum cliques. We also prove that cores of both $\Gamma_d$ and its complement are maximum cliques. The graph $\Gamma_d$ is useful for error-correcting codes. We show that every largest independent set of $\Gamma_d$ is both an MRD code over $\mathbb{Z}_{p^s}$ and a usual MDS code. Moreover, there is a largest independent set of $\Gamma_d$ to be a linear code over $\mathbb{Z}_{p^s}$.


title: Impact of ejecta morphology and composition on the electromagnetic signatures of neutron star mergers

abstract: The electromagnetic transients accompanying compact binary mergers ($\gamma$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer calculations. We employ the radiative transfer code SuperNu and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts ($z$). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the FRDM model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed nIR excess after GRB130603B, a more optimistic heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models.


title: A special class of congruences on $\kappa$-frames

abstract: Madden has shown that in contrast to the situation with frames, the smallest dense quotient of a $\kappa$-frame need not be Boolean. We characterise these so-called d-reduced $\kappa$-frames as those which may be embedded as a generating sub-$\kappa$-frame of a Boolean frame. We introduce the notion of the closure of a $\kappa$-frame congruence and call a congruence clear if it is the largest congruence with a given closure. These ideas are used to prove $\kappa$-frame analogues of known results concerning Boolean frame quotients. In particular, we show that d-reduced $\kappa$-frames are precisely the quotients of $\kappa$-frames by clear congruences and that every $\kappa$-frame congruence is the meet of clear congruences.


title: Estimating Accuracy from Unlabeled Data: A Probabilistic Logic Approach

abstract: We propose an efficient method to estimate the accuracy of classifiers using only unlabeled data. We consider a setting with multiple classification problems where the target classes may be tied together through logical constraints. For example, a set of classes may be mutually exclusive, meaning that a data instance can belong to at most one of them. The proposed method is based on the intuition that: (i) when classifiers agree, they are more likely to be correct, and (ii) when the classifiers make a prediction that violates the constraints, at least one classifier must be making an error. Experiments on four real-world data sets produce accuracy estimates within a few percent of the true accuracy, using solely unlabeled data. Our models also outperform existing state-of-the-art solutions in both estimating accuracies, and combining multiple classifier outputs. The results emphasize the utility of logical constraints in estimating accuracy, thus validating our intuition.


title: From 4d Ambitwistor Strings to On Shell Diagrams and Back

abstract: We investigate the relation between 4d ambitwistor string theory and on-shell diagrams for planar N=4 super-Yang-Mills and N=8 supergravity, and deduce several new results about their scattering amplitudes at tree-level and 1-loop. In particular, we derive new Grassmannian integral formulae for tree-level amplitudes and obtain new worldsheet formulae for 1-loop amplitudes which are manifestly supersymmetric and supported on scattering equations refined by MHV degree.


title: Semantics of higher inductive types

abstract: Higher inductive types are a class of type-forming rules, introduced to provide basic (and not-so-basic) homotopy-theoretic constructions in a type-theoretic style. They have proven very fruitful for the "synthetic" development of homotopy theory within type theory, as well as in formalizing ordinary set-level mathematics in type theory. In this article, we construct models of a wide range of higher inductive types in a fairly wide range of settings. We introduce the notion of cell monad with parameters: a semantically-defined scheme for specifying homotopically well-behaved notions of structure. We then show that any suitable model category has *weakly stable typal initial algebras* for any cell monad with parameters. When combined with the local universes construction to obtain strict stability, this specializes to give models of specific higher inductive types, including spheres, the torus, pushout types, truncations, the James construction, and general localisations. Our results apply in any sufficiently nice Quillen model category, including any right proper, simplicially locally cartesian closed, simplicial Cisinski model category (such as simplicial sets) and any locally presentable locally cartesian closed category (such as sets) with its trivial model structure. In particular, any locally presentable locally cartesian closed $(\infty,1)$-category is presented by some model category to which our results apply.


title: A Numerical Relativity Waveform Surrogate Model for Generically Precessing Binary Black Hole Mergers

abstract: A generic, non-eccentric binary black hole (BBH) system emits gravitational waves (GWs) that are completely described by 7 intrinsic parameters: the black hole spin vectors and the ratio of their masses. Simulating a BBH coalescence by solving Einstein's equations numerically is computationally expensive, requiring days to months of computing resources for a single set of parameter values. Since theoretical predictions of the GWs are often needed for many different source parameters, a fast and accurate model is essential. We present the first surrogate model for GWs from the coalescence of BBHs including all $7$ dimensions of the intrinsic non-eccentric parameter space. The surrogate model, which we call NRSur7dq2, is built from the results of $744$ numerical relativity simulations. NRSur7dq2 covers spin magnitudes up to $0.8$ and mass ratios up to $2$, includes all $\ell \leq 4$ modes, begins about $20$ orbits before merger, and can be evaluated in $\sim~50\,\mathrm{ms}$. We find the largest NRSur7dq2 errors to be comparable to the largest errors in the numerical relativity simulations, and more than an order of magnitude smaller than the errors of other waveform models. Our model, and more broadly the methods developed here, will enable studies that would otherwise require millions of numerical relativity waveforms, such as parameter inference and tests of general relativity with GW observations.


title: The VLA-COSMOS 3~GHz Large Project: Cosmic evolution of radio AGN and implications for radio-mode feedback since z~5

abstract: Based on a sample of over 1,800 radio AGN at redshifts out to z~5, which have typical stellar masses within ~3x(10^{10}-10^{11}) Msol, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L_1.4GHz ~ 10^{22}-10^{27} W/Hz) out to z~5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Phi*~(1+z)^{(2.00+/-0.18)-(0.60+/-0.14)z}, and L*~(1+z)^{(2.88+/-0.82)-(0.84+/-0.34)z}, respectively, with a turnover in number and luminosity densities of the population at z~1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z~5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.


title: CNN-Based Joint Clustering and Representation Learning with Feature Drift Compensation for Large-Scale Image Data

abstract: Given a large unlabeled set of images, how to efficiently and effectively group them into clusters based on extracted visual representations remains a challenging problem. To address this problem, we propose a convolutional neural network (CNN) to jointly solve clustering and representation learning in an iterative manner. In the proposed method, given an input image set, we first randomly pick k samples and extract their features as initial cluster centroids using the proposed CNN with an initial model pre-trained from the ImageNet dataset. Mini-batch k-means is then performed to assign cluster labels to individual input samples for a mini-batch of images randomly sampled from the input image set until all images are processed. Subsequently, the proposed CNN simultaneously updates the parameters of the proposed CNN and the centroids of image clusters iteratively based on stochastic gradient descent. We also proposed a feature drift compensation scheme to mitigate the drift error caused by feature mismatch in representation learning. Experimental results demonstrate the proposed method outperforms start-of-the-art clustering schemes in terms of accuracy and storage complexity on large-scale image sets containing millions of images.


title: Wealth dynamics in a sentiment-driven market

abstract: We study dynamics of a simulated world with stock and money, driven by the externally given processes which we refer to as sentiments. The considered sentiments influence the buy/sell stock trading attitude, the perceived price uncertainty, and the trading intensity of all or a part of the market participants. We study how the wealth of market participants evolves in time in such an environment. We discuss the opposite perspective in which the parameters of the sentiment processes can be inferred a posteriori from the observed market behavior.


title: Incompressible wind accretion

abstract: We present a simple, analytic model for the accretion flow of an incompressible wind onto a gravitating object. This solution corresponds to the Newtonian limit of a previously known relativistic model for a fluid obeying a stiff equation of state for which the sound speed is constant everywhere and equal to the speed of light. The new solution should be useful as a benchmark test for numerical hydrodynamics codes and, moreover, it can be used as an illustrative example in a gas dynamics course.


title: Criticality for charged black branes

abstract: We show that the inclusion of higher curvature terms in the gravitational action can lead to phase transitions and critical behaviour for charged black branes. The higher curvature terms considered here belong to the recently constructed generalized quasi-topological class [ArXiv: 1703.01631], which possess a number of interesting properties, such as being ghost-free on constant curvature backgrounds and non-trivial in four dimensions. We show that critical behaviour is a generic feature of the black branes in all dimensions $d \ge 4$, and contextualize the results with a review of the properties of black branes in Lovelock and quasi-topological gravity, where critical behaviour is not possible. These results may have interesting implications for the CFTs dual to this class of theories.


title: Induction of Interpretable Possibilistic Logic Theories from Relational Data

abstract: The field of Statistical Relational Learning (SRL) is concerned with learning probabilistic models from relational data. Learned SRL models are typically represented using some kind of weighted logical formulas, which make them considerably more interpretable than those obtained by e.g. neural networks. In practice, however, these models are often still difficult to interpret correctly, as they can contain many formulas that interact in non-trivial ways and weights do not always have an intuitive meaning. To address this, we propose a new SRL method which uses possibilistic logic to encode relational models. Learned models are then essentially stratified classical theories, which explicitly encode what can be derived with a given level of certainty. Compared to Markov Logic Networks (MLNs), our method is faster and produces considerably more interpretable models.