title: Deep Neural Networks for Surface Segmentation Meet Conditional Random Fields

abstract: Automated surface segmentation is important and challenging in many medical image analysis applications. Recent deep learning based methods have been developed for various object segmentation tasks. Most of them are a classification based approach (e.g., U-net), which predicts the probability of being target object or background for each voxel. One problem of those methods is lacking of topology guarantee for segmented objects, and usually post processing is needed to infer the boundary surface of the object. In this paper, a novel model based on 3-D convolutional neural networks (CNNs) and Conditional Random Fields (CRFs) is proposed to tackle the surface segmentation problem with end-to-end training. To the best of our knowledge, this is the first study to apply a 3-D neural network with a CRFs model for direct surface segmentation. Experiments carried out on NCI-ISBI 2013 MR prostate dataset and Medical Segmentation Decathlon Spleen dataset demonstrated promising segmentation results.


title: Asymptotic analysis of exit time for dynamical systems with a single well potential

abstract: We study the exit time from a bounded multi-dimensional domain $\Omega$ of the stochastic process $\mathbf{Y}_\varepsilon=\mathbf{Y}_\varepsilon(t,a)$, $t\geqslant 0$, $a\in \mathcal{A}$, governed by the overdamped Langevin dynamics \begin{equation*} d\mathbf{Y}_\varepsilon =-\nabla V(\mathbf{Y}_\varepsilon) dt +\sqrt{2}\varepsilon\, d\mathbf{W}, \qquad \mathbf{Y}_\varepsilon(0,a)\equiv x\in\Omega \end{equation*} where $\varepsilon$ is a small positive parameter, $\mathcal{A}$ is a sample space, $\mathbf{W}$ is a $n$-dimensional Wiener process. The exit time corresponds to the first hitting of $\partial\Omega$ by the trajectories of the above dynamical system and the expectation value of this exit time solves the boundary value problem \begin{equation*} (-\varepsilon^2\Delta +\nabla V\cdot \nabla)u_\varepsilon=1\quad\text{in}\quad\Omega,\qquad u_\varepsilon=0\quad\text{on}\quad\partial\Omega. \end{equation*} We assume that the function $V$ is smooth enough and has the only minimum at the origin (contained in $\Omega$); the minimum can be degenerate. At other points of $\Omega$, the gradient of $V$ is non-zero and the normal derivative of $V$ at the boundary $\partial\Omega$ does not vanish as well. Our main result is a complete asymptotic expansion for $u_\varepsilon$ as well as for the lowest eigenvalue of the considered problem and for the associated eigenfunction. The asymptotics for $u_\varepsilon$ involves a term exponentially large $\varepsilon$; we find this term in a closed form. Apart of this term, we also construct a power in $\varepsilon$ asymptotic expansion such that this expansion and a mentioned exponentially large term approximate $u_\varepsilon$ up to arbitrarily power of $\varepsilon$. We also discuss some probabilistic aspects of our results.


title: Learning the Graphical Structure of Electronic Health Records with Graph Convolutional Transformer

abstract: Effective modeling of electronic health records (EHR) is rapidly becoming an important topic in both academia and industry. A recent study showed that using the graphical structure underlying EHR data (e.g. relationship between diagnoses and treatments) improves the performance of prediction tasks such as heart failure prediction. However, EHR data do not always contain complete structure information. Moreover, when it comes to claims data, structure information is completely unavailable to begin with. Under such circumstances, can we still do better than just treating EHR data as a flat-structured bag-of-features? In this paper, we study the possibility of jointly learning the hidden structure of EHR while performing supervised prediction tasks on EHR data. Specifically, we discuss that Transformer is a suitable basis model to learn the hidden EHR structure, and propose Graph Convolutional Transformer, which uses data statistics to guide the structure learning process. The proposed model consistently outperformed previous approaches empirically, on both synthetic data and publicly available EHR data, for various prediction tasks such as graph reconstruction and readmission prediction, indicating that it can serve as an effective general-purpose representation learning algorithm for EHR data.


title: XENON1T Dark Matter Data Analysis: Signal Reconstruction, Calibration and Event Selection

abstract: The XENON1T experiment at the Laboratori Nazionali del Gran Sasso is the most sensitive direct detection experiment for dark matter in the form of weakly interacting particles (WIMPs) with masses above $6\,$GeV/$c^2$ scattering off nuclei. The detector employs a dual-phase time projection chamber with 2.0 metric tons of liquid xenon in the target. A one metric $\mathrm{ton}\times\mathrm{year}$ exposure of science data was collected between October 2016 and February 2018. This article reports on the performance of the detector during this period and describes details of the data analysis that led to the most stringent exclusion limits on various WIMP-nucleon interaction models to date. In particular, signal reconstruction, event selection and calibration of the detector response to nuclear and electronic recoils in XENON1T are discussed.


title: Skyrmion-electronics: Writing, deleting, reading and processing magnetic skyrmions toward spintronic applications

abstract: The field of magnetic skyrmions has been actively investigated across a wide range of topics during the last decades. In this topical review, we mainly review and discuss key results and findings in skyrmion research since the first experimental observation of magnetic skyrmions in 2009. We particularly focus on the theoretical, computational and experimental findings and advances that are directly relevant to the spintronic applications based on magnetic skyrmions, i.e. their writing, deleting, reading and processing driven by magnetic field, electric current and thermal energy. We then review several potential applications including information storage, logic computing gates and non-conventional devices such as neuromorphic computing devices. Finally, we discuss possible future research directions on magnetic skyrmions, which also cover rich topics on other topological textures such as antiskyrmions and bimerons in antiferromagnets and frustrated magnets.


title: The $h^*$-polynomials of locally anti-blocking lattice polytopes and their $\gamma$-positivity

abstract: A lattice polytope $\mathcal{P} \subset \mathbb{R}^d$ is called a locally anti-blocking polytope if for any closed orthant $\mathbb{R}^d_{\varepsilon}$ in $\mathbb{R}^d$, $\mathcal{P} \cap \mathbb{R}^d_{\varepsilon}$ is unimodularly equivalent to an anti-blocking polytope by reflections of coordinate hyperplanes. In the present paper, we give a formula for the $h^*$-polynomials of locally anti-blocking lattice polytopes. In particular, we discuss the $\gamma$-positivity of the $h^*$-polynomials of locally anti-blocking reflexive polytopes.


title: Snapshot projection optical tomography

abstract: We present a plenoptic microscopy configuration for 3D snapshot imaging, which is dual telecentric and can directly record true projection images corresponding with different viewing angles. It also allows blocking high-angle stray rays without sacrificing the light collection efficiency. This configuration named as snapshot projection optical tomography (SPOT) arranges an objective lens and a microlens array (MLA) in a 4-f telecentric configuration and places an aperture stop at the back focal plane of a relay lens. We develop a forward imaging model for SPOT, which can also be applied to existing light field microscopy techniques using an MLA as tube lens. Using the developed system, we demonstrate snapshot 3D imaging of various fluorescent beads and a biological cell, which confirms the capability of SPOT for imaging specimens with an extended fluorophore distribution as well as isolated fluorochromes. The transverse and vertical resolutions are measured to be 0.8 $\mu$m and 1.6 $\mu$m, respectively.


title: Data-Free Quantization Through Weight Equalization and Bias Correction

abstract: We introduce a data-free quantization method for deep neural networks that does not require fine-tuning or hyperparameter selection. It achieves near-original model performance on common computer vision architectures and tasks. 8-bit fixed-point quantization is essential for efficient inference on modern deep learning hardware. However, quantizing models to run in 8-bit is a non-trivial task, frequently leading to either significant performance reduction or engineering time spent on training a network to be amenable to quantization. Our approach relies on equalizing the weight ranges in the network by making use of a scale-equivariance property of activation functions. In addition the method corrects biases in the error that are introduced during quantization. This improves quantization accuracy performance, and can be applied to many common computer vision architectures with a straight forward API call. For common architectures, such as the MobileNet family, we achieve state-of-the-art quantized model performance. We further show that the method also extends to other computer vision architectures and tasks such as semantic segmentation and object detection.


title: A Hybrid multiphase model based on lattice Boltzmann method direct simulations

abstract: By means of the multicomponent Shan-Chen lattice Boltzmann method (LBM), we investigate the multiphase flow through porous media. Despite the excellent accuracy of the LBM, large domains result in unaffordable computational expenses. The Hybrid model developed in this study is based on a pore-network (PN) approach that enhances a decomposition of the granular assembly into small subsets (pore throats). Lattice Boltzmann simulations are performed for each pore throat to determine the hydrodynamic properties (entry capillary pressure, primary drainage curve, liquid morphology, etc) at the microscale. The local properties obtained with LBM are incorporated at the network to solve the larger-scale problem. This strategy leads to a significant decrease of the computation time at the sample-scale compared to a fully resolved method. Fluid morphology and phase distribution are evaluated during the drainage of a small granular assembly using the Hybrid model (PN-LBM). Results are contrasted with those obtained in a fully resolved simulation (LBM). The agreement between the two models illustrates the capability of the Hybrid method, which combines the efficiency of the PN approach and the accuracy of the LBM at the pore scale.


title: Double-peaked OI profile: a likely signature of the gaseous ring around KH 15D

abstract: KH 15D is a well-known spectroscopic binary because of its unique and dramatic photometric variability. The variability is explained by a circumbinary dust ring but the ring itself was never directly detected. We present a new interpretation of the double-peaked [O~{\scriptsize I}]\,$\lambda$6300 profiles as originating from the hot disk surface of KH 15D. By modeling these profiles, we measure emitting radii between $\sim$0.5--5 au, basically a gaseous ring very similar in radial extent to the dust ring inferred from modeling the system's photometric variability. We also discuss the possibility that external photoevaporation driven by UV photons from the nearby massive star HD~47887 has truncated the outer edge of the disk to the observed value.


title: Large Scale Structure of Neural Network Loss Landscapes

abstract: There are many surprising and perhaps counter-intuitive properties of optimization of deep neural networks. We propose and experimentally verify a unified phenomenological model of the loss landscape that incorporates many of them. High dimensionality plays a key role in our model. Our core idea is to model the loss landscape as a set of high dimensional \emph{wedges} that together form a large-scale, inter-connected structure and towards which optimization is drawn. We first show that hyperparameter choices such as learning rate, network width and $L_2$ regularization, affect the path optimizer takes through the landscape in a similar ways, influencing the large scale curvature of the regions the optimizer explores. Finally, we predict and demonstrate new counter-intuitive properties of the loss-landscape. We show an existence of low loss subspaces connecting a set (not only a pair) of solutions, and verify it experimentally. Finally, we analyze recently popular ensembling techniques for deep networks in the light of our model.


title: Clouds of Oriented Gradients for 3D Detection of Objects, Surfaces, and Indoor Scene Layouts

abstract: We develop new representations and algorithms for three-dimensional (3D) object detection and spatial layout prediction in cluttered indoor scenes. We first propose a clouds of oriented gradient (COG) descriptor that links the 2D appearance and 3D pose of object categories, and thus accurately models how perspective projection affects perceived image gradients. To better represent the 3D visual styles of large objects and provide contextual cues to improve the detection of small objects, we introduce latent support surfaces. We then propose a "Manhattan voxel" representation which better captures the 3D room layout geometry of common indoor environments. Effective classification rules are learned via a latent structured prediction framework. Contextual relationships among categories and layout are captured via a cascade of classifiers, leading to holistic scene hypotheses that exceed the state-of-the-art on the SUN RGB-D database.


title: What Kind of Language Is Hard to Language-Model?

abstract: How language-agnostic are current state-of-the-art NLP tools? Are there some types of language that are easier to model with current methods? In prior work (Cotterell et al., 2018) we attempted to address this question for language modeling, and observed that recurrent neural network language models do not perform equally well over all the high-resource European languages found in the Europarl corpus. We speculated that inflectional morphology may be the primary culprit for the discrepancy. In this paper, we extend these earlier experiments to cover 69 languages from 13 language families using a multilingual Bible corpus. Methodologically, we introduce a new paired-sample multiplicative mixed-effects model to obtain language difficulty coefficients from at-least-pairwise parallel corpora. In other words, the model is aware of inter-sentence variation and can handle missing data. Exploiting this model, we show that "translationese" is not any easier to model than natively written language in a fair comparison. Trying to answer the question of what features difficult languages have in common, we try and fail to reproduce our earlier (Cotterell et al., 2018) observation about morphological complexity and instead reveal far simpler statistics of the data that seem to drive complexity in a much larger sample.


title: A Massive Magnetic Helium Atmosphere White Dwarf Binary in a Young Star Cluster

abstract: We have searched the Gaia DR2 catalogue for previously unknown hot white dwarfs in the direction of young open star clusters. The aim of this experiment was to try and extend the initial-final mass relation (IFMR) to somewhat higher masses, potentially providing a tension with the Chandrasekhar limit currently thought to be around 1.38 M$_{\odot}$. We discovered a particularly interesting white dwarf in the direction of the young $\sim$150 Myr old cluster Messier 47 (NGC 2422). All Gaia indicators (proper motion, parallax, location in the Gaia colour-magnitude diagram) suggest that it is a cluster member. Its spectrum, obtained from Gemini South, yields a number of anomalies: it is a DB (helium-rich atmosphere) white dwarf, it has a large magnetic field (2.5 MG), is of high mass ($\sim$1.06 M$_\odot$) and its colours are very peculiar --- particularly the redder ones ($r$, $i$, $z$ and $y$), which suggest that it has a late-type companion. This is the only magnetized, detached binary white dwarf with a non-degenerate companion of any spectral type known in or out of a star cluster. If the white dwarf is a cluster member, as all indicators suggest, its progenitor had a mass just over 6 M$_\odot$. It may, however, be telling an even more interesting story than the one related to the IFMR, one about the origin of stellar magnetic fields, Type I supernovae and gravitational waves from low mass stellar systems.


title: Shapes and Context: In-the-Wild Image Synthesis & Manipulation

abstract: We introduce a data-driven approach for interactively synthesizing in-the-wild images from semantic label maps. Our approach is dramatically different from recent work in this space, in that we make use of no learning. Instead, our approach uses simple but classic tools for matching scene context, shapes, and parts to a stored library of exemplars. Though simple, this approach has several notable advantages over recent work: (1) because nothing is learned, it is not limited to specific training data distributions (such as cityscapes, facades, or faces); (2) it can synthesize arbitrarily high-resolution images, limited only by the resolution of the exemplar library; (3) by appropriately composing shapes and parts, it can generate an exponentially large set of viable candidate output images (that can say, be interactively searched by a user). We present results on the diverse COCO dataset, significantly outperforming learning-based approaches on standard image synthesis metrics. Finally, we explore user-interaction and user-controllability, demonstrating that our system can be used as a platform for user-driven content creation.


title: B-mode Targets

abstract: The CMB-S4 Science Book, Astro2020 Science White Paper on gravitational waves, and PICO report contain an extensive discussion of many ongoing and planned efforts in the search for gravitational waved produced by inflation. Here we give a short executive summary of the results obtained in our recent papers, which specify the simplest available inflationary models providing physically motivated targets for these searches. Our conclusions are specific for the $10^{-3} \lesssim r \lesssim 10^{-2}$ range, where we present the B-mode benchmarks of the U-duality symmetric class of $\alpha$-attractors, and for $ r \lesssim 10^{-3}$, where we present B-mode targets, for which the future precision measurements of $n_s$ will be decisive. We show that a combination of the simplest $\alpha$-attractors and KKLTI models of D-brane inflation covers most of the area favored by Planck 2018.


title: Constraining primordial non-Gaussianity using two galaxy surveys and CMB lensing

abstract: Next-generation galaxy surveys will be able to measure perturbations on scales beyond the equality scale. On these ultra-large scales, primordial non-Gaussianity leaves signatures that can shed light on the mechanism by which perturbations in the early Universe are generated. We perform a forecast analysis for constraining local type non-Gaussianity and its two-parameter extension with a simple scale-dependence. We combine different clustering measurements from future galaxy surveys -- a 21cm intensity mapping survey and two photometric galaxy surveys -- via the multi-tracer approach. Furthermore we then include CMB lensing from a CMB Stage 4 experiment in the multi-tracer, which can improve the constraints on bias parameters. We forecast $\sigma(f_{\rm NL}) \simeq 0.9$ (1.4) by combining SKA1, a Euclid-like (LSST-like) survey, and CMB-S4 lensing. With CMB lensing, the precision on $f_{\rm NL}$ improves by up to a factor of 2, showing that a joint analysis is important. In the case with running of $f_{\rm NL}$, our results show that the combination of upcoming cosmological surveys could achieve $\sigma(n_{\rm NL}) \simeq 0.12$ (0.22) on the running index.


title: Two toy models for the motion of a leaky tank car

abstract: I present two toy models for the motion of a tank car as fluid drains out of an off-center opening, based on replacing the fluid with particles. No knowledge beyond simple Lagrangian mechanics is required. The first toy model is solved analytically and the second can be simulated numerically for a variety of initial conditions and mass ratios. Both toy models show the car turning around once as is required by conservation of momentum.


title: Convergence analysis of a Crank-Nicolson Galerkin method for an inverse source problem for parabolic equations with boundary observations

abstract: This work is devoted to an inverse problem of identifying a source term depending on both spatial and time variables in a parabolic equation from single Cauchy data on a part of the boundary. A Crank-Nicolson Galerkin method is applied to the least squares functional with an quadratic stabilizing penalty term. The convergence of finite dimensional regularized approximations to the sought source as measurement noise levels and mesh sizes approach to zero with an appropriate regularization parameter is proved. Moreover, under a suitable source condition, an error bound and corresponding convergence rates are proved. Finally, several numerical experiments are presented to illustrate the theoretical findings.


title: DualDICE: Behavior-Agnostic Estimation of Discounted Stationary Distribution Corrections

abstract: In many real-world reinforcement learning applications, access to the environment is limited to a fixed dataset, instead of direct (online) interaction with the environment. When using this data for either evaluation or training of a new policy, accurate estimates of discounted stationary distribution ratios -- correction terms which quantify the likelihood that the new policy will experience a certain state-action pair normalized by the probability with which the state-action pair appears in the dataset -- can improve accuracy and performance. In this work, we propose an algorithm, DualDICE, for estimating these quantities. In contrast to previous approaches, our algorithm is agnostic to knowledge of the behavior policy (or policies) used to generate the dataset. Furthermore, it eschews any direct use of importance weights, thus avoiding potential optimization instabilities endemic of previous methods. In addition to providing theoretical guarantees, we present an empirical study of our algorithm applied to off-policy policy evaluation and find that our algorithm significantly improves accuracy compared to existing techniques.


title: Incremental Classifier Learning Based on PEDCC-Loss and Cosine Distance

abstract: The main purpose of incremental learning is to learn new knowledge while not forgetting the knowledge which have been learned before. At present, the main challenge in this area is the catastrophe forgetting, namely the network will lose their performance in the old tasks after training for new tasks. In this paper, we introduce an ensemble method of incremental classifier to alleviate this problem, which is based on the cosine distance between the output feature and the pre-defined center, and can let each task to be preserved in different networks. During training, we make use of PEDCC-Loss to train the CNN network. In the stage of testing, the prediction is determined by the cosine distance between the network latent features and pre-defined center. The experimental results on EMINST and CIFAR100 show that our method outperforms the recent LwF method, which use the knowledge distillation, and iCaRL method, which keep some old samples while training for new task. The method can achieve the goal of not forgetting old knowledge while training new classes, and solve the problem of catastrophic forgetting better.


title: On the Universality of Noiseless Linear Estimation with Respect to the Measurement Matrix

abstract: In a noiseless linear estimation problem, one aims to reconstruct a vector x* from the knowledge of its linear projections y=Phi x*. There have been many theoretical works concentrating on the case where the matrix Phi is a random i.i.d. one, but a number of heuristic evidence suggests that many of these results are universal and extend well beyond this restricted case. Here we revisit this problematic through the prism of development of message passing methods, and consider not only the universality of the l1 transition, as previously addressed, but also the one of the optimal Bayesian reconstruction. We observed that the universality extends to the Bayes-optimal minimum mean-squared (MMSE) error, and to a range of structured matrices.

categories: stat.ML cs.IT cs.LG eess.SP math.IT math.ST stat.TH


title: Improving Reproducible Deep Learning Workflows with DeepDIVA

abstract: The field of deep learning is experiencing a trend towards producing reproducible research. Nevertheless, it is still often a frustrating experience to reproduce scientific results. This is especially true in the machine learning community, where it is considered acceptable to have black boxes in your experiments. We present DeepDIVA, a framework designed to facilitate easy experimentation and their reproduction. This framework allows researchers to share their experiments with others, while providing functionality that allows for easy experimentation, such as: boilerplate code, experiment management, hyper-parameter optimization, verification of data integrity and visualization of data and results. Additionally, the code of DeepDIVA is well-documented and supported by several tutorials that allow a new user to quickly familiarize themselves with the framework.


title: Dealing with Non-Stationarity in Multi-Agent Deep Reinforcement Learning

abstract: Recent developments in deep reinforcement learning are concerned with creating decision-making agents which can perform well in various complex domains. A particular approach which has received increasing attention is multi-agent reinforcement learning, in which multiple agents learn concurrently to coordinate their actions. In such multi-agent environments, additional learning problems arise due to the continually changing decision-making policies of agents. This paper surveys recent works that address the non-stationarity problem in multi-agent deep reinforcement learning. The surveyed methods range from modifications in the training procedure, such as centralized training, to learning representations of the opponent's policy, meta-learning, communication, and decentralized learning. The survey concludes with a list of open problems and possible lines of future research.


title: Normal curves on a smooth immersed surface

abstract: The aim of this paper is to investigate the sufficient condition for the invariance of a normal curve on a smooth immersed surface under isometry. We also find the the deviations of the tangential and normal components of the curve with respect to the given isometry.


title: Trip Table Estimation and Prediction for Dynamic Traffic Assignment Applications

abstract: The study focuses on estimating and predicting time-varying origin to destination (OD) trip tables for a dynamic traffic assignment (DTA) model. A bi-level optimisation problem is formulated and solved to estimate OD flows from pre-existent demand matrix and historical traffic flow counts. The estimated demand is then considered as an input for a time series OD demand prediction model to support the DTA model for short-term traffic condition forecasting. Results show a high capability of the proposed OD demand estimation method to reduce the DTA model error through an iterative solution algorithm. Moreover, the applicability of the OD demand prediction approach is investigated for an incident analysis application for a major corridor in Sydney, Australia.


title: Phase-field material point method for dynamic brittle fracture with isotropic and anisotropic surface energy

abstract: A novel phase field material point method is introduced for robust simulation of dynamic fracture in elastic media considering the most general case of anisotropic surface energy. Anisotropy is explicitly introduced through a properly defined crack density functional. The particular case of impact driven fracture is treated by employing a discrete field approach within the material point method setting. In this, the equations of motion and phase field governing equations are solved independently for each discrete field using a predictor-corrector algorithm. Contact at the interface is resolved through frictional contact conditions. The proposed method is verified using analytical predictions. The influence of surface energy anisotropy and loading conditions on the resulting crack paths is assessed through a set of benchmark problems. Comparisons are made with the standard Phase Field Finite Element Method and experimental observations.


title: Measuring clock jumps using pulsar timing

abstract: In this paper, we investigate the statistical signal-processing algorithm to measure the instant local clock jump from the timing data of multiple pulsars. Our algorithm is based on the framework of Bayesian statistics. In order to make the Bayesian algorithm applicable with limited computational resources, we dedicated our efforts to the analytic marginalization of irrelevant parameters. We found that the widely used parameter for pulsar timing systematics, the `Efac' parameter, can be analytically marginalized. This reduces the Gaussian likelihood to a function very similar to the Student's $t$-distribution. Our iterative method to solve the maximum likelihood estimator is also explained in the paper. Using pulsar timing data from the Yunnan Kunming 40m radio telescope, we demonstrate the application of the method, where 80-ns level precision for the clock jump can be achieved. Such a precision is comparable to that of current commercial time transferring service using satellites. We expect that the current method could help developing the autonomous pulsar time scale.


title: Mass Loss from the Exoplanet WASP-12b Inferred from $\textit{Spitzer}$ Phase Curves

abstract: The exoplanet WASP-12b is the prototype for the emerging class of ultra-hot, Jupiter-mass exoplanets. Past models have predicted---and near ultra-violet observations have shown---that this planet is losing mass. We present an analysis of two sets of 3.6 $\mu$m and 4.5 $\mu$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu$m, not at 3.6 $\mu$m, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, $\lesssim$ 600 K gas. It is unclear why WASP-12b is the only ultra-hot Jupiter to exhibit this mass loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of $6.4\sigma$ ($46.2^{\circ}$) at 3.6 $\mu$m.


title: R-Parity Violating Supersymmetry and the 125 GeV Higgs signals

abstract: We study the impact of R-parity violating Supersymmetry (RPV SUSY) on the 125 GeV Higgs production and decay modes at the LHC. We assume a heavy SUSY spectrum with multi-TeV squarks and SU(2) scalar singlets as well as the decoupling limit in the SUSY Higgs sector. In this case the lightest CP-even Higgs is SM-like when R-parity is conserved. In contrast, we show that when R-parity violating interactions are added to the SUSY framework, significant deviations may occur in some production and decay channels of the 125 GeV Higgs-like state. Indeed, we assume a single-flavor (mostly third generation) Bilinear RPV (BRPV) interactions, which generate Higgs-sneutrino mixing, lepton-chargino mixing and neutrino-neutralino mixing, and find that notable deviations of ${\cal O}(20-30\%)$ may be expected in the Higgs signal strength observables in some channels, e.g., in $p p \to h \to \mu^+ \mu^-, \tau^+ \tau^-$. Moreover, we find that new and detectable signals associated with BRPV Higgs decays to gauginos, $h\rightarrow\nu_{\tau}\tilde{\chi}_{2}^{0}$ and $h\rightarrow\tau^{\pm}\chi_{2}^{\mp}$, may also arise in this scenario. These decays yield a typical signature of $h \to \tau^\pm \ell^\mp + {{\not\!\! E_T}}$ ($\ell =e,\mu,\tau$) that can be much larger than in the SM, and may also be accompanied by an ${\cal O}(20-30\%)$ enhancement in the di-photon signal $pp \to h \to \gamma \gamma$. We also examine potential interesting signals of Trilinear R-parity violation (TRPV) interactions in the production and decays of the Higgs-sneutrino BRPV mixed state (assuming it is the 125 GeV scalar) and show that, in this case also, large deviations up to ${\cal O}(100\%)$ are expected in e.g., $p p \to h \to \mu^+ \mu^-, \tau^+ \tau^-$, which are sensitive to the BRPV$\times$TRPV coupling product.


title: Inquiring into the nature of the Abell 2667 Brightest Cluster Galaxy: physical properties from MUSE

abstract: Based on HST and MUSE data, we probe the stellar and gas properties (i.e. kinematics, stellar mass, star formation rate) of the radio-loud brightest cluster galaxy (BCG) located at the centre of the X-ray luminous cool core cluster Abell 2667 (z = 0.2343). The bi-dimensional modelling of the BCG surface brightness profile reveals the presence of a complex system of substructures extending all around the galaxy. Clumps of different size and shape plunged into a more diffuse component constitute these substructures, whose intense 'blue' optical colour hints to the presence of a young stellar population. Our results depict the BCG as a massive (M_star ~ 1.38 x 10^11 M_sun) dispersion-supported spheroid (v_star < 150 km/s, sigma_0 ~ 216 km/s) hosting an active supermassive black hole (M_SMBH ~ 3.8 x 10^9 M_sun) whose optical features are typical of low ionisation nuclear emission line regions. Although the velocity pattern of the stars in the BCG is irregular, the stellar kinematics in the regions of the clumps show a positive velocity of ~ 100 km/s, similarly to the gas component. An analysis of the mechanism giving rise to the observed lines in the clumps through empirical diagnostic diagrams points out that the emission is composite, suggesting the contribution from both star formation and AGN. We conclude our analysis describing how scenarios of both chaotic cold accretion and merging with a gas-rich disc galaxy can efficaciously explain the phenomena the BCG is undergoing.


title: Holographic Bulk Reconstruction And Cosmological Singularities

abstract: We study the structure of entanglement wedges in the Kasner-AdS geometry, which provides an example of AdS/CFT engineered cosmological singularity. We investigate the specific limitations of causal reconstruction methods, imposed by the presence of the cosmological singularities, and we show the supremacy of modular reconstruction. This model provides an example where modular reconstruction based on a proper operator subalgebra is more powerful than the strongest possible causal reconstruction, based on the complete operator algebra.


title: The GIST Pipeline: A Multi-Purpose Tool for the Analysis and Visualisation of (Integral-field) Spectroscopic Data

abstract: We present a convenient, all-in-one framework for the scientific analysis of fully reduced, (integral-field) spectroscopic data. The GIST pipeline (Galaxy IFU Spectroscopy Tool) is entirely written in Python3 and conducts all steps from the preparation of input data, over the scientific analysis to the production of publication-quality plots. In its basic setup, it extracts stellar kinematics, performs an emission-line analysis and derives stellar population properties from full spectral fitting as well as via the measurement of absorption line-strength indices by exploiting the well-known pPXF and GandALF routines, where the latter has now been implemented in Python. The pipeline is not specific to any instrument or analysis technique and provides easy means of modification and further development, as of its modular code architecture. An elaborate, Python-native parallelisation is implemented and tested on various machines. The software further features a dedicated visualization routine with a sophisticated graphical user interface. This allows an easy, fully-interactive plotting of all measurements, spectra, fits, and residuals, as well as star formation histories and the weight distribution of the models. The pipeline has successfully been applied to both low and high-redshift data from MUSE, PPAK (CALIFA), and SINFONI, as well as to simulated data for HARMONI@ELT and WEAVE and is currently being used by the TIMER, Fornax3D, and PHANGS collaborations. We demonstrate its capabilities by applying it to MUSE TIMER observations of NGC 1433.


title: Cooper pairing of incoherent electrons: an electron-phonon version of the Sachdev-Ye-Kitaev model

abstract: We introduce and solve a model of interacting electrons and phonons that is a natural generalization of the Sachdev-Ye-Kitaev-model and that becomes superconducting at low temperatures. In the normal state two Non-Fermi liquid fixed points with distinct universal exponents emerge. At weak coupling superconductivity prevents the onset of low-temperature quantum criticality, reminiscent of the behavior in several heavy-electron and iron-based materials. At strong coupling, pairing of highly incoherent fermions sets in deep in the Non-Fermi liquid regime, a behavior qualitatively similar to that in underdoped cuprate superconductors. The pairing of incoherent time-reversal partners is protected by a mechanism similar to Anderson's theorem for disordered superconductors. The superconducting ground state is characterized by coherent quasiparticle excitations and higher-order bound states thereof, revealing that it is no longer an ideal gas of Cooper pairs, but a strongly coupled pair fluid. The normal-state incoherency primarily acts to suppress the weight of the superconducting coherence peak and reduce the condensation energy. Based on this we expect strong superconducting fluctuations, in particular at strong coupling.


title: Distinguishing cosmologies using the turn-around radius near galaxy clusters

abstract: Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $\Lambda$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $\Lambda$CDM universe due to the large error-bars arising from the non-trivial cluster environments.


title: Joint 3D Localization and Classification of Space Debris using a Multispectral Rotating Point Spread Function

abstract: We consider the problem of joint three-dimensional (3D) localization and material classification of unresolved space debris using a multispectral rotating point spread function (RPSF). The use of RPSF allows one to estimate the 3D locations of point sources from their rotated images acquired by a single 2D sensor array, since the amount of rotation of each source image about its x, y location depends on its axial distance z. Using multi-spectral images, with one RPSF per spectral band, we are able not only to localize the 3D positions of the space debris but also classify their material composition. We propose a three-stage method for achieving joint localization and classification. In Stage 1, we adopt an optimization scheme for localization in which the spectral signature of each material is assumed to be uniform, which significantly improves efficiency and yields better localization results than possible with a single spectral band. In Stage 2, we estimate the spectral signature and refine the localization result via an alternating approach. We process classification in the final stage. Both Poisson noise and Gaussian noise models are considered, and the implementation of each is discussed. Numerical tests using multispectral data from NASA show the efficiency of our three-stage approach and illustrate the improvement of point source localization and spectral classification from using multiple bands over a single band.


title: Constraining primordial black hole masses with the isotropic gamma ray background

abstract: Primordial black holes can represent all or most of the dark matter in the window $10^{17}-10^{22}\,$g. Here we present an extension of the constraints on PBHs of masses $10^{13}-10^{18}\,$g arising from the isotropic diffuse gamma ray background. Primordial black holes evaporate by emitting Hawking radiation that should not exceed the observed background. Generalizing from monochromatic distributions of Schwarzschild black holes to extended mass functions of Kerr rotating black holes, we show that the lower part of this mass window can be closed for near-extremal black holes.


title: Star-forming clumps in the Lyman Alpha Reference Sample of galaxies -- I. Photometric analysis and clumpiness

abstract: We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $\alpha = -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $\Sigma_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $\Sigma_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/\sigma_0$). We find that in galaxies with higher $\Sigma_{SFR}$ or lower $v_s/\sigma_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$\alpha$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution.


title: Hall Viscosity in a Strongly Coupled Magnetized Plasma

abstract: We show how a Hall viscosity induced by a magnetic field can be generated in strongly coupled theories with a holographic dual. This is achieved by considering parity-breaking higher derivative terms in the gravity dual. These terms couple the Riemann curvature tensor to the field strength of a gauge field dual to the charge current, and have an analog in the field theory side as a coupling between the "Euler current" and the electromagnetic field. As a concrete example, we study the effect of the new terms in the thermodynamic and transport properties of a strongly coupled magnetized plasma dual to a dyonic black hole in $AdS_4$. As a new property of the holographic model, we find that for a state that is initially neutral at zero magnetic field, a charge density and non-dissipative Hall transport are present when the magnetic field is turned on. Remarkably, we also observe that the results from the holographic model are consistent with hydrodynamics even at magnetic fields much larger than temperature.


title: Measuring and exploiting the cloud consolidation of the Web

abstract: We present measurements showing that the top one million most popular Web domains are reachable within 13ms (in the median) from a collective of just 12 cloud data centers. We explore the consequences of this Web "consolidation", focusing on its potential for speeding up the evolution of the Web. That most popular services reside in or near a small number of data centers implies that new application and transport technologies can be rapidly deployed for these Web services, without the involvement of their operators. We show how this may be achieved by orchestrating a handful of reverse proxies deployed in the same data centers, with new technologies deployed at these proxies being nearly as effective as deploying them directly to the Web servers. We present early measurements of this approach, demonstrating a >50% reduction in Web page load times for users with high latencies to Web servers. We also show that this model, using a small number of proxies, can be surprisingly competitive with extensive CDN deployments, especially in geographies with high last-mile latencies.


title: Improved regularity for the porous medium equation along zero level-sets

abstract: In the present work we establish sharp regularity estimates for the solutions of the porous medium equation, along their zero level-sets. We work under a proximity regime on the exponent governing the nonlinearity of the problem. Then, we prove that solutions are locally of class $\mathcal{C}^{1-,\frac{1}{2}-}$ along free boundary points $x_0\in \partial\left\lbrace u>0\right\rbrace$, both in time and space. Our argument consists of importing information from the heat equation, through approximation and localization methods.


title: New Physics Probes: Atomic Parity Violation, Polarized Electron Scattering and Neutrino-Nucleus Coherent Scattering

abstract: Atomic Parity Violation (APV) is usually quantified in terms of the weak nuclear charge $Q_W$ of a nucleus, which depends on the coupling strength between the atomic electrons and quarks. In this work, we review the importance of APV to probing new physics using effective field theory. Furthermore, using $SU(2)$ invariance, we correlate our findings with those from neutrino-nucleus coherent scattering. Moreover, we investigate signs of parity violation in polarized electron scattering and show how precise measurements on the Weinberg angle, $\sin \theta_W$, will give rise to competitive bounds on light mediators over a wide range of masses and interactions strength. Lastly, apply our bounds to several models namely, Dark Z, Two Higgs Doublet Model-$U(1)_X$ and 3-3-1, considering both light and heavy mediator regimes.


title: Dynamic triad interactions and evolving turbulence spectra

abstract: We investigate the effect of a four-dimensional Fourier transform on the formulation of Navier-Stokes equation in Fourier space and the way the energy is transferred between Fourier components (modal interactions, commonly referred to as triad interactions in the classical 3-dimensional analysis). We specifically consider the effect of a spatially and temporally finite, digitally sampled velocity record on the modal interactions and find that Fourier components may interact within a broadened frequency window as compared to the usual integrals over infinite ranges. We also see how these finite velocity records have a significant effect on the efficiency of the different modal interactions and thereby on the shape and development of velocity power spectra. The observation that mismatches in the wavevector triadic interactions may be compensated by a corresponding mismatch in the frequencies supports the empirically deduced delayed interactions reported in [Josserand \textit{et al.}, \textit{J. Stat. Phys.} (2017)]. Collectively, these results explain the occurrence and time development of the so-called Richardson cascade and also why deviations from the classical Richardson cascade may occur. Finally, we quote results from companion papers that deal with measurements and computer simulations of the time development of velocity power spectra in a turbulent jet flow into which a single Fourier mode (narrow band oscillation) is injected.


title: Universal scaling of active nematic turbulence

abstract: A landmark of turbulence is the emergence of universal scaling laws, such as Kolmogorov's $E(q)\sim q^{-5/3}$ scaling of the kinetic energy spectrum of inertial turbulence with the wave vector $q$. In recent years, active fluids have been shown to exhibit turbulent-like flows at low Reynolds number. However, the existence of universal scaling properties in these flows has remained unclear. To address this issue, here we propose a minimal defect-free hydrodynamic theory for two-dimensional active nematic fluids at vanishing Reynolds number. By means of large-scale simulations and analytical arguments, we show that the kinetic energy spectrum exhibits a universal scaling $E(q)\sim q^{-1}$ at long wavelengths. We find that the energy injection due to activity has a peak at a characteristic length scale, which is selected by a nonlinear mechanism. In contrast to inertial turbulence, energy is entirely dissipated at the scale where it is injected, thus precluding energy cascades. Nevertheless, the non-local character of the Stokes flow establishes long-ranged velocity correlations, which lead to the scaling behavior. We conclude that active nematic fluids define a distinct universality class of turbulence at low Reynolds number.


title: Bridging the Physics and Chemistry of graphene. From Huckel aromaticity to Dirac cones and topological insulators

abstract: By bridging graphene and benzene through a well-defined sequence of polycyclic aromatic hydrocarbons and their inherent shell structure, J. Phys. Chem. C, 2018, 122, 17526, it is shown that graphene is actually a coherent arrangement of interwoven benzene molecules, which are coordinated by aromaticity, shell structure and topology, all interrelated. The exotic properties of graphene are in fact macroscopic manifestations of aromaticity. At the molecular level this is revealed as a dynamical flipping of the atomic pz-orbitals belonging to different sublattices, leading to a dynamical interchange between aromatic and non-aromatic rings, associated with frontier occupied-unoccupied orbitals or valence-conduction bands interchange with opposite parities, in full agreement with many-body and topological insulators theory. Such dynamical interchange or coupling is driven by inversion symmetry-breaking, activated through geometrical frustration between sublattice and space group symmetry. This, for rectangular nanographenes, in contrast to hexagonal, leads to topological gapless edge states, which, contrary to opposite claims, are not spin polarized.


title: Globular clusters as probes of dark matter cusp-core transformations

abstract: Bursty star formation in dwarf galaxies can slowly transform a steep dark matter cusp into a constant density core. We explore the possibility that globular clusters (GCs) retain a dynamical memory of this transformation. To test this, we use the nbody6df code to simulate the dynamical evolution of GCs, including stellar evolution, orbiting in static and time-varying potentials for a Hubble time. We find that GCs orbiting within a cored dark matter halo, or within a halo that has undergone a cusp-core transformation, grow to a size that is substantially larger ($R_{\rm eff} > 10$ pc) than those in a static cusped dark matter halo. They also produce much less tidal debris. We find that the cleanest signal of an historic cusp-core transformation is the presence of large GCs with tidal debris. However, the effect is small and will be challenging to observe in real galaxies. Finally, we qualitatively compare our simulated GCs with the observed GC populations in the Fornax, NGC 6822, IKN and Sagittarius dwarf galaxies. We find that the GCs in these dwarf galaxies are systematically larger ($\langle R_{\rm eff}\rangle \simeq 7.8$ pc), and have substantially more scatter in their sizes, than in-situ metal rich GCs in the Milky Way and young massive star clusters forming in M83 ($\langle R_{\rm eff} \rangle \simeq 2.5$ pc). We show that the size, scatter and survival of GCs in dwarf galaxies are all consistent with them having evolved in a constant density core, or a potential that has undergone a cusp-core transformation, but not in a dark matter cusp.


title: Unsupervised Discovery of Gendered Language through Latent-Variable Modeling

abstract: Studying the ways in which language is gendered has long been an area of interest in sociolinguistics. Studies have explored, for example, the speech of male and female characters in film and the language used to describe male and female politicians. In this paper, we aim not to merely study this phenomenon qualitatively, but instead to quantify the degree to which the language used to describe men and women is different and, moreover, different in a positive or negative way. To that end, we introduce a generative latent-variable model that jointly represents adjective (or verb) choice, with its sentiment, given the natural gender of a head (or dependent) noun. We find that there are significant differences between descriptions of male and female nouns and that these differences align with common gender stereotypes: Positive adjectives used to describe women are more often related to their bodies than adjectives used to describe men.


title: PerspectroScope: A Window to the World of Diverse Perspectives

abstract: This work presents PerspectroScope, a web-based system which lets users query a discussion-worthy natural language claim, and extract and visualize various perspectives in support or against the claim, along with evidence supporting each perspective. The system thus lets users explore various perspectives that could touch upon aspects of the issue at hand.The system is built as a combination of retrieval engines and learned textual-entailment-like classifiers built using a few recent developments in natural language understanding. To make the system more adaptive, expand its coverage, and improve its decisions over time, our platform employs various mechanisms to get corrections from the users. PerspectroScope is available at github.com/CogComp/perspectroscope.


title: Deep 2FBSDEs For Systems With Control Multiplicative Noise

abstract: We present a deep recurrent neural network architecture to solve a class of stochastic optimal control problems described by fully nonlinear Hamilton Jacobi Bellmanpartial differential equations. Such PDEs arise when one considers stochastic dynamics characterized by uncertainties that are additive and control multiplicative. Stochastic models with the aforementioned characteristics have been used in computational neuroscience, biology, finance and aerospace systems and provide a more accurate representation of actuation than models with additive uncertainty. Previous literature has established the inadequacy of the linear HJB theory and instead rely on a non-linear Feynman-Kac lemma resulting in a second order forward-backward stochastic differential equations representation. However, the proposed solutions that use this representation suffer from compounding errors and computational complexity leading to lack of scalability. In this paper, we propose a deep learning based algorithm that leverages the second order Forward-Backward SDE representation and LSTM based recurrent neural networks to not only solve such Stochastic Optimal Control problems but also overcome the problems faced by previous approaches and scales well to high dimensional systems. The resulting control algorithm is tested on non-linear systems in robotics and biomechanics to demonstrate feasibility and out-performance against previous methods.


title: Estimating the Number of Fatal Victims of the Peruvian Internal Armed Conflict, 1980-2000: an application of modern multi-list Capture-Recapture techniques

abstract: We estimate the number of fatal victims of the Peruvian internal armed conflict between 1980-2000 using stratified seven-list Capture-Recapture methods based on Dirichlet process mixtures, which we extend to accommodate incomplete stratification information. We use matched data from six sources, originally analyzed by the Peruvian Truth and Reconciliation Commission in 2003, together with a new large dataset, originally published in 2006 by the Peruvian government. We deal with missing stratification labels by developing a general framework and estimation methods based on MCMC sampling for jointly fitting generic Bayesian Capture-Recapture models and the missing labels. Through a detailed exploration driven by domain-knowledge, modeling and refining, with special precautions to avoid cherry-picking of results, we arrive to a conservative posterior estimate of 58,234 (CI95% = [56,741, 61,289]), and a more liberal estimate of 65,958 (CI95% = [61,462, 75,387]) fatal victims. We also determine that the Shining Path guerrillas killed more people than the Peruvian armed forces. We additionally explore and discuss estimates based on log-linear modeling and multiple-imputation. We finish by discussing several lessons learned about the use of Capture-Recapture methods for estimating casualties in conflicts.


title: Formation of structure for gas-liquid, non-equilibrium polymer media

abstract: This paper examines the mechanisms of destruction and synthesis for macromolecules, which may be propelled by external acoustic fields in the flows of polymeric liquids containing a large number of gas bubbles. The dynamics of these bubbles is assumed to govern by changing the flow geometry and exciting sound oscillations in the flow. Mechanically-induced kinetic changes in macromolecules (destruction and synthesis of polymer chains) will occur when the bubbles collapse.


title: The Prolog debugger and declarative programming

abstract: Logic programming is a declarative programming paradigm. Programming language Prolog makes logic programming possible, at least to a substantial extent. However the Prolog debugger works solely in terms of the operational semantics. So it is incompatible with declarative programming. This report discusses this issue and tries to find how the debugger may be used from the declarative point of view. The results are rather not encouraging. Also, the box model of Byrd, used by the debugger, is explained in terms of SLD-resolution.


title: Evolution speed of open quantum dynamics

abstract: The space of density matrices is embedded in a Euclidean space to deduce the dynamical equation satisfied by the state of an open quantum system. The Euclidean norm is used to obtain an explicit expression for the speed of the evolution of the state. The unitary contribution to the evolution speed is given by the modified skew information of the Hamiltonian, while the radial component of the evolution speed, connected to the rate at which the purity of the state changes, is shown to be determined by the modified skew information of the Lindblad operators. An open-system analogue of the quantum navigation problem is posed, and a perturbative analysis is presented to identify the amount of change on the speed. Properties of the evolution speed are examined further through example systems, showing that the evolution speed need not be a decreasing function of time.


title: Background field method and nonlinear gauges

abstract: We present a reformulation of the background field method for Yang-Mills type theories, based on using a superalgebra of generators of BRST and background field transformations. The new approach enables one to implement and consistently use non-linear gauges in a natural way, by using the requirement of invariance of the fermion gauge-fixing functional under the background field transformations.


title: Enumerating linear systems on graphs

abstract: The divisor theory of graphs views a finite connected graph $G$ as a discrete version of a Riemann surface. Divisors on $G$ are formal integral combinations of the vertices of $G$, and linear equivalence of divisors is determined by the discrete Laplacian operator for $G$. As in the case of Riemann surfaces, we are interested in the complete linear system $|D|$ of a divisor $D$---the collection of nonnegative divisors linearly equivalent to $D$. Unlike the case of Riemann surfaces, the complete linear system of a divisor on a graph is always finite. We compute generating functions encoding the sizes of all complete linear systems on $G$ and interpret our results in terms of polyhedra associated with divisors and in terms of the invariant theory of the (dual of the) Jacobian group of $G$. If $G$ is a cycle graph, our results lead to a bijection between complete linear systems and binary necklaces. The final section generalizes our results to a model based on integral $M$-matrices.


title: The radiation field on product cones

abstract: We consider the wave equation on a product cone and find a joint asymptotic expansion for solutions near null and future infinities. The rates of decay seen in the expansion at future infinity are the resonances of a hyperbolic cone and were computed by the authors in a previous paper. The expansion treats an asymptotic regime not considered in the influential work of Cheeger and Taylor. The main result adapts the blueprint laid out in the asymptotically Minkowski setting with key new elements including propagation estimates near the conic singularities. The proof of the propagation estimates extends prior work of Melrose-Vasy-Wunsch and Gannot-Wunsch.


title: Dissecting a disk-instability outburst in a symbiotic star: $NuSTAR$, and $Swift$ observations of T Coronae Borealis during the rise to the "super-active" state

abstract: The current $super-active$ state of the recurrent nova T CrB has been observed with unprecedented detail. Previously published observations provide strong evidence that this state is due to an enhancement of the flow of material through the accretion disk, which increased the optical depth of its most internal region, the boundary layer. $NuSTAR$ and $Swift$ observed T CrB in 2015 September, roughly halfway through the rise to optical maximum. In our analysis of these data, we have found that: $i$) the UV emission, as observed with $Swift$/UVOT in 2015, was already as bright as it became in 2017, after the optical peak; $ii$) the soft X-ray emission (E $\lesssim$ 0.6 keV) observed in 2017 after the optical peak, on the other hand, had not yet developed during the rising phase in 2015; $iii$) the hard X-ray emitting plasma (E $\gtrsim$ 2 keV) had the same temperature and about half the flux of that observed during quiescence in 2006. This phenomenology is akin to that observed during dwarf novae in outburst, but with the changes in the spectral energy distribution happening on a far longer time scale.


title: Deep Forward-Backward SDEs for Min-max Control

abstract: This paper presents a novel approach to numerically solve stochastic differential games for nonlinear systems. The proposed approach relies on the nonlinear Feynman-Kac theorem that establishes a connection between parabolic deterministic partial differential equations and forward-backward stochastic differential equations. Using this theorem the Hamilton-Jacobi-Isaacs partial differential equation associated with differential games is represented by a system of forward-backward stochastic differential equations. Numerical solution of the aforementioned system of stochastic differential equations is performed using importance sampling and a Long-Short Term Memory recurrent neural network, which is trained in an offline fashion. The resulting algorithm is tested on two example systems in simulation and compared against the standard risk neutral stochastic optimal control formulations.


title: A Systematic Comparison of English Noun Compound Representations

abstract: Building meaningful representations of noun compounds is not trivial since many of them scarcely appear in the corpus. To that end, composition functions approximate the distributional representation of a noun compound by combining its constituent distributional vectors. In the more general case, phrase embeddings have been trained by minimizing the distance between the vectors representing paraphrases. We compare various types of noun compound representations, including distributional, compositional, and paraphrase-based representations, through a series of tasks and analyses, and with an extensive number of underlying word embeddings. We find that indeed, in most cases, composition functions produce higher quality representations than distributional ones, and they improve with computational power. No single function performs best in all scenarios, suggesting that a joint training objective may produce improved representations.


title: Parametrized spectra, a low-tech approach

abstract: We give an alternative treatment of the foundations of parametrized spectra, with an eye toward applications in fixed-point theory. We cover most of the central results from the book of May and Sigurdsson, sometimes with weaker hypotheses, and give a new construction of the bicategory $\mathcal Ex$ of parametrized spectra. We also give a careful account of coherence results at the level of homotopy categories. This is an expository work, akin to a set of lecture notes, encompassing and extending the more novel material in the paper "A convenient category of parametrized spectra." Our primary goal is to give careful and explicit explanations for the standard facts surrounding the bicategory of parametrized spectra and the Reidemeister trace. In particular, our proof of Costenoble-Waner duality gives a more explicit account of how the concrete index formulas of Dold are related to the abstract theory of bicategorical duality from May and Sigurdsson, Ponto, and Shulman.


title: Issues with post-hoc counterfactual explanations: a discussion

abstract: Counterfactual post-hoc interpretability approaches have been proven to be useful tools to generate explanations for the predictions of a trained blackbox classifier. However, the assumptions they make about the data and the classifier make them unreliable in many contexts. In this paper, we discuss three desirable properties and approaches to quantify them: proximity, connectedness and stability. In addition, we illustrate that there is a risk for post-hoc counterfactual approaches to not satisfy these properties.


title: The Performance Of Convolutional Coding Based Cooperative Communication: Relay Position And Power Allocation Analysis

abstract: Wireless communication faces adversities due to noise, fading, and path loss. Multiple-Input Multiple-Output (MIMO) systems are used to overcome individual fading effect by employing transmit diversity. Duo to user single-antenna, Cooperation between at least two users is able to provide spatial diversity. This paper presents the evaluation of the performances of the Amplify and Forward (AF) cooperative system for different relay positions using several network topologies over Rayleigh and Rician fading channel. Furthermore, we present the performances of AF cooperative system with various power allocation. The results show that cooperative communication with convolutional coding shows an outperformance compared to the non-convolutional, which is a promising solution for high data-rate networks such as (WSN), Ad hoc, (IoT), and even mobile networks. When topologies are compared, the simulation shows that, linear topology offers the best BER performance, in contrast when the relay acts as source and the source take the relay place, the analysis result shows that, equilateral triangle topology has the best BER performance and stability, and the system performance with inter-user Rician fading channel is better than the performance of the system with inter-user Rayleigh fading channel.


title: Distribution-Free Multisample Test Based on Optimal Matching with Applications to Single Cell Genomics

abstract: In this paper we propose a nonparametric graphical test based on optimal matching, for assessing the equality of multiple unknown multivariate probability distributions. Our procedure pools the data from the different classes to create a graph based on the minimum non-bipartite matching, and then utilizes the number of edges connecting data points from different classes to examine the closeness between the distributions. The proposed test is exactly distribution-free (the null distribution does not depend on the distribution of the data) and can be efficiently applied to multivariate as well as non-Euclidean data, whenever the inter-point distances are well-defined. We show that the test is universally consistent, and prove a distributional limit theorem for the test statistic under general alternatives. Through simulation studies, we demonstrate its superior performance against other common and well-known multisample tests. In scenarios where our test suggests distributional differences across classes, we also propose an approach for identifying which class or group contributes to this overall difference. The method is applied to single cell transcriptomics data obtained from the peripheral blood, cancer tissue, and tumor-adjacent normal tissue of human subjects with hepatocellular carcinoma and non-small-cell lung cancer. Our method unveils patterns in how biochemical metabolic pathways are altered across immune cells in a cancer setting, depending on the tissue location. All of the methods described herein are implemented in the R package multicross.


title: Estimating Homogeneous Data-driven BRDF Parameters from a Reflectance Map under Known Natural Lighting

abstract: In this paper we demonstrate robust estimation of the model parameters of a fully-linear data-driven BRDF model from a reflectance map under known natural lighting. To regularize the estimation of the model parameters, we leverage the reflectance similarities within a material class. We approximate the space of homogeneous BRDFs using a Gaussian mixture model, and assign a material class to each Gaussian in the mixture model. We formulate the estimation of the model parameters as a non-linear maximum a-posteriori optimization, and introduce a linear approximation that estimates a solution per material class from which the best solution is selected. We demonstrate the efficacy and robustness of our method using the MERL BRDF database under a variety of natural lighting conditions, and we provide a proof-of-concept real-world experiment.


title: Generating Pareto optimal dose distributions for radiation therapy treatment planning

abstract: Radiotherapy treatment planning currently requires many trail-and-error iterations between the planner and treatment planning system, as well as between the planner and physician for discussion/consultation. The physician's preferences for a particular patient cannot be easily quantified and precisely conveyed to the planner. In this study we present a real-time volumetric Pareto surface dose generation deep learning neural network that can be used after segmentation by the physician, adding a tangible and quantifiable end-point to portray to the planner. From 70 prostate patients, we first generated 84,000 intensity modulated radiation therapy plans (1,200 plans per patient) sampling the Pareto surface, representing various tradeoffs between the planning target volume (PTV) and the organs-at-risk (OAR), including bladder, rectum, left femur, right femur, and body. We divided the data to 10 test patients and 60 training/validation patients. We then trained a hierarchically densely connected convolutional U-net (HD U-net), to take the PTV and avoidance map representing OARs masks and weights, and predict the optimized plan. The HD U-net is capable of accurately predicting the 3D Pareto optimal dose distributions, with average [mean, max] dose errors of [3.4%, 7.7%](PTV), [1.6%, 5.6%](bladder), [3.7%, 4.2%](rectum), [3.2%, 8.0%](left femur), [2.9%, 7.7%](right femur), and [0.04%, 5.4%](body) of the prescription dose. The PTV dose coverage prediction was also very similar, with errors of 1.3% (D98) and 2.0% (D99). Homogeneity was also similar, differing by 0.06 on average. The neural network can predict the dose within 1.7 seconds. Clinically, the optimization and dose calculation is much slower, taking 5-10 minutes.


title: Positivity of the fundamental solution for fractional diffusion and wave equations

abstract: We study the question of positivity of the fundamental solution for fractional diffusion and wave equations of the form, which may be of fractional order both in space and time. We give a complete characterization for the positivity of the fundamental solution in terms of the order of the time derivative $\alpha\in(0,2)$, the order of the spatial derivative $\beta\in (0,2]$ and the spatial dimension $d$. It turns out that the fundamental solution fails to be positive for all $\alpha\in (1,2)$, and either $\beta\in (0,2]$ and $d\ge 2$ or $\beta<\alpha$ and $d=1$, whereas in the other cases it remains positive. The proof is based on delicate properties of the Fox H-functions and the Mittag-Leffler functions.


title: Superheat conductivity and electric activity of superfluid systems

abstract: In this article it is shown that in He\nobreakdash-II placed in a magnetic field $\mathbf{H}$, the heat flow under the influence of a temperature gradient $\nabla T$ causes appearance of an electric field $\mathbf{E}\sim\mathbf{H}\times\nabla T$. This effect occurs in superfluid dielectric systems due to their specific property named superheat conductivity. The magnitude of the electric field significantly depends on the shape of the helium sample and the magnetic field direction relative to the sample. The effect exists both at static and non-stationary temperature gradient (e. g. during propagation of second sound).


title: Path Cohomology of Locally Finite Digraphs,Hodge's Theorem and the $p$-Lazy Random Walk

abstract: The study of Markov chains on discrete spaces, such as digraphs, has captivated mathematicians in recent decades due to its interconnectedness with topology, geometry, dynamics, spectral theory, and differential equations. Furthermore, extensive exploration of these multifaceted relationships has been pursued for their practical utility in diverse fields, including machine learning and image segmentation. In recent times, these interrelations have been generalized to higher dimensions within the framework of finite-dimensional simplicial complexes. In this paper, we embark on a further extension of these concepts. Initially, we introduce a cohomology of infinite (though locally finite) digraphs in arbitrary dimensions. Subsequently, in the latter portion of this manuscript, we define a fresh family of Laplace operators and conduct an examination of their spectrum, culminating in the proof of the Hodge Decomposition Theorem within this framework. Finally, we conclude by presenting a Markov chain, the $p$-Lazy Random Walk, whose asymptotic behavior is intrinsically linked to these cohomologies, while its mixing time is related to the the spectrum of our Laplace operators. This development opens doors to numerous unexplored questions, particularly regarding potential generalizations of the Ollivier-Ricci curvature to this topology and these Laplacians.


title: Second-best Beam-Alignment via Bayesian Multi-Armed Bandits

abstract: Millimeter-wave (mm-wave) systems rely on narrow-beams to cope with the severe signal attenuation in the mm-wave frequency band. However, susceptibility to beam mis-alignment due to mobility or blockage requires the use of beam-alignment schemes, with huge cost in terms of overhead and use of system resources. In this paper, a beam-alignment scheme is proposed based on Bayesian multi-armed bandits, with the goal to maximize the alignment probability and the data-communication throughput. A Bayesian approach is proposed, by considering the state as a posterior distribution over angles of arrival (AoA) and of departure (AoD), given the history of feedback signaling and of beam pairs scanned by the base-station (BS) and the user-end (UE). A simplified sufficient statistics for optimal control is derived, in the form of preference of BS-UE beam pairs. By bounding a value function, the second-best preference policy is derived, which strikes an optimal balance between exploration and exploitation by selecting the beam pair with the current second-best preference. Through Monte-Carlo simulation with analog beamforming, the superior performance of the second-best preference policy is demonstrated in comparison to existing schemes based on first-best preference, linear Thompson sampling, and upper confidence bounds, with up to 7%, 10% and 30% improvements in alignment probability, respectively.


title: Instability of flux flow and production of vortex-antivortex pairs by current-driven Josephson vortices in layered superconductors

abstract: We report numerical simulations of the nonlinear dynamics of Josephson vortices driven by strong dc currents in layered superconductors. Dynamic equations for interlayer phase differences in a stack of coupled superconducting layers were solved to calculate a drag coefficient $\eta(J)$ of the vortex as a function of the perpendicular dc current density $J$. It is shown that Cherenkov radiation produced by a moving vortex causes significant radiation drag increasing $\eta(v)$ at high vortex velocities $v$ and striking instabilities of driven Josephson vortices moving faster than a terminal velocity $v_c$. The steady-state flux flow breaks down at $v>v_c$ as the vortex starts producing a cascade of expanding vortex-antivortex pairs evolving into either planar macrovortex structures or branching flux patterns propagating both along and across the layers. The pair production triggered by a rapidly moving vortex is most pronounced in a stack of underdamped planar junctions where it can occur at $J>J_s$ well below the interlayer Josephson critical current density. Both $v_c$ and $J_s$ were calculated as functions of the quasiparticle damping parameter, and the dc magnetic field applied parallel to the layers. The effects of vortex interaction on the Cherenkov instability of moving vortex chains and lattices in annular stacks of Josephson junctions were considered. It is shown that a vortex driven by a current density $J>J_s$ in a multilayer of finite length excites self-sustained large-amplitude standing waves of magnetic flux, resulting in temporal oscillations of the total magnetic moment. We evaluated a contribution of this effect to the power $W$ radiated by the sample and showed that $W$ increases strongly as the number of layers increases.


title: Stability of Graph Scattering Transforms

abstract: Scattering transforms are non-trainable deep convolutional architectures that exploit the multi-scale resolution of a wavelet filter bank to obtain an appropriate representation of data. More importantly, they are proven invariant to translations, and stable to perturbations that are close to translations. This stability property dons the scattering transform with a robustness to small changes in the metric domain of the data. When considering network data, regular convolutions do not hold since the data domain presents an irregular structure given by the network topology. In this work, we extend scattering transforms to network data by using multiresolution graph wavelets, whose computation can be obtained by means of graph convolutions. Furthermore, we prove that the resulting graph scattering transforms are stable to metric perturbations of the underlying network. This renders graph scattering transforms robust to changes on the network topology, making it particularly useful for cases of transfer learning, topology estimation or time-varying graphs.


title: Stein's method and the distribution of the product of zero mean correlated normal random variables

abstract: Over the last 80 years there has been much interest in the problem of finding an explicit formula for the probability density function of two zero mean correlated normal random variables. Motivated by this historical interest, we use a recent technique from the Stein's method literature to obtain a simple new proof, which also serves as an exposition of a general method that may be useful in related problems.


title: Identifying resonances of the Galactic bar in Gaia DR2: I. Clues from action space

abstract: Action space synthesizes the orbital information of stars and is well-suited to analyse the rich kinematic substructure of the disc in the \emph{Gaia} DR2 radial velocity sample (RVS). We revisit the strong perturbation induced in the Milky Way (MW) disc by an $m=2$ bar, using test particle simulations and the actions $(J_R,L_z,J_z)$ estimated in an axisymmetric potential. These make three useful diagnostics cleanly visible. (1.) We use the well-known characteristic flip from outward to inward motion at the Outer Lindblad Resonance (OLR, $l=+1,m=2$), which occurs along the axisymmetric resonance line (ARL) in $(L_z,J_R)$, to identify in the \emph{Gaia} action data three candidates for the bar's OLR and pattern speed $\Omega_\text{bar}$: $1.85\Omega_0$, $1.20\Omega_0$, and $1.63\Omega_0$ (with $\sim0.1\Omega_0$ systematic uncertainty). The \emph{Gaia} data is therefore consistent with both slow and fast bar models in the literature, but disagrees with recent measurements of $\sim1.45\Omega_0$. (2.) For the first time, we demonstrate that bar resonances -- especially the OLR -- cause a gradient in vertical action $\langle J_z \rangle$ with $L_z$ around the ARL via "$J_z$-sorting" of stars. This could contribute to the observed coupling of $\langle v_R \rangle$ and $\langle | v_z | \rangle$ in the Galactic disc. (3.) We confirm prior results that the behaviour of resonant orbits is well approximated by scattering and oscillation in $(L_z,J_R)$ along a slope $\Delta J_R/\Delta L_z = l/m$ centered on the $l$:$m$ ARL. Overall, we demonstrate that axisymmetrically estimated actions are a powerful diagnostic tool even in non-axisymmetric systems.


title: Power Gradient Descent

abstract: The development of machine learning is promoting the search for fast and stable minimization algorithms. To this end, we suggest a change in the current gradient descent methods that should speed up the motion in flat regions and slow it down in steep directions of the function to minimize. It is based on a "power gradient", in which each component of the gradient is replaced by its versus-preserving $H$-th power, with $0<H<1$. We test three modern gradient descent methods fed by such variant and by standard gradients, finding the new version to achieve significantly better performances for the Nesterov accelerated gradient and AMSGrad. We also propose an effective new take on the ADAM algorithm, which includes power gradients with varying $H$.


title: CMB foreground measurements through broad-band radio spectro-polarimetry: prospects of the SKA-MPG telescope

abstract: Precise measurement of the foreground synchrotron emission, which contaminates the faint polarized cosmic microwave background radiation (CMB), is a major challenge for the next-generation of CMB experiments. To address this, dedicated foreground measurement experiments are being undertaken at radio frequencies between 2 and 40 GHz. Foreground polarized synchrotron emission measurements are particularly challenging, primarily due to the complicated frequency dependence in the presence of Faraday rotation, and are best recovered through broad fractional-bandwidth polarization measurements at frequencies $\lesssim 5$ GHz. A unique opportunity for measuring the foreground polarized synchrotron emission will be provided by the 15-m SKA-MPG telescope operating in the frequency range 1.7 to 3.5~GHz (S-Band). Here, we present the scope of a Southern sky survey in S-Band at 1 degree angular resolution and explore its added advantage for application of powerful techniques, such as, Stokes $Q$, $U$ fitting and RM-synthesis. A full Southern-sky polarization survey with this telescope, when combined with other on-going efforts at slightly higher frequencies, will provide an excellent frequency coverage for modeling and extrapolating the foreground polarized synchrotron emission to CMB frequencies ($\gtrsim80$~GHz) with rms brightness temperature better than 10~nK per 1 degree$^2$. We find that this survey will be crucial for understanding the effects of Faraday depolarization, especially in low Galactic latitude regions. This will allow better foreground cleaning and thus will contribute significantly in further improving component separation analyses and increase usable sky area for cosmological analysis of the \textit{Planck} data, and the \textit{LiteBIRD} mission in the future.


title: On pro-$p$ groups with quadratic cohomology

abstract: The main purpose of this article is to study pro-$p$ groups with quadratic $\mathbb{F}_p$-cohomology algebra, i.e. $H^\bullet$-quadratic pro-$p$ groups. Prime examples of such groups are the maximal Galois pro-$p$ groups of fields containing a primitive root of unity of order $p$. We show that the amalgamated free product and HNN-extension of $H^\bullet$-quadratic pro-$p$ groups is $H^\bullet$-quadratic, under certain necessary conditions. Moreover, we introduce and investigate a new family of pro-$p$ groups that yields many new examples of $H^\bullet$-quadratic groups: $p$-RAAGs. These examples generalise right angled Artin groups in the category of pro-$p$ groups. Finally, we explore "Tits alternative behaviour" of $H^\bullet$-quadratic pro-$p$ groups.


title: Mathematical and numerical analysis of a nonlocal Drude model in nanoplasmonics

abstract: In this paper, we consider the frequency-domain Maxwell's equations coupled to a nonlocal Drude model which describes the nonlocal optical response in metallic nanostructures. We prove the existence and uniqueness of weak solutions to the coupled equations. A Galerkin finite element method based on the Raviart--Thomas and N\'{e}d\'{e}lec elements is proposed to solve the equations and the associated error estimates are given. This is the first work on the mathematical and numerical analysis of this model. Numerical examples are presented to verify our theoretical analysis.


title: Discovery of a 2.8 s pulsar in a 2 d orbit High-Mass X-ray Binary powering the Ultraluminous X-ray source ULX-7 in M51

abstract: We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).


title: Quantization of dynamical symplectic reduction

abstract: A long-standing problem in quantum gravity and cosmology is the quantization of systems in which time evolution is generated by a constraint that must vanish on solutions. Here, an algebraic formulation of this problem is presented, together with new structures and results, which show that specific conditions need to be satisfied in order for well-defined evolution to be possible.


title: Phase Crystals

abstract: Superconductivity owes its properties to the phase of the electron pair condensate that breaks the $U(1)$ symmetry. In the most traditional ground state, the phase is uniform and rigid. The normal state can be unstable towards special inhomogeneous superconducting states: the Abrikosov vortex state, and the Fulde-Ferrell-Larkin-Ovchinnikov state. Here we show that the phase-uniform superconducting state can go into a fundamentally different and more ordered non-uniform ground state, that we denote as a phase crystal. The new state breaks translational invariance through formation of a spatially periodic modulation of the phase, manifested by unusual superflow patterns and circulating currents, that also break time-reversal symmetry. We list the general conditions needed for realization of phase crystals. Using microscopic theory we then derive an analytic expression for the superfluid density tensor for the case of a non-uniform environment in a semi-infinite superconductor. We demonstrate how the surface quasiparticle states enter the superfluid density and identify phase crystallization as the main player in several previous numerical observations in unconventional superconductors, and predict existence of a similar phenomenon in superconductor-ferromagnetic structures. This analytic approach provides a new unifying aspect for the exploration of boundary-induced quasiparticles and collective excitations in superconductors. More generally, we trace the origin of phase crystallization to non-local properties of the gradient energy, which implies existence of similar pattern-forming instabilities in many other contexts.


title: A comparative study of $2d$ Ising model at different boundary conditions using non-deterministic Hexagonal Cellular Automata

abstract: The spin system of the $2d$ Ising model having a hexagonal-lattice is simulated using non-deterministic Cellular Automata. The method to implement this program is outlined and our results show a good approximation to the exact analytic solution. The phase transition in $2d$ Ising model is studied with a $40\times40$ hexagonal-lattice with five different boundary conditions (bcs) i.e., adiabatic, periodic, reflexive, fixed $+1$ and fixed $-1$ with random orientation of spins as initial conditions in the absence of an external applied magnetic field. The critical temperature below which the spontaneous magnetization appears as well as other physical quantities such as the magnetization, energy, specific heat, susceptibility and entropy with each of the bcs are calculated. The phase transition occurs around $T^H_c$ = 1.5 which approximates well with the result obtained from exact analytic solution by Wannier and Houtappel. We compare the behavior of magnetisation per cell for five different types of bcs by calculating the number of points close to the line of zero magnetization for $T>T^H_c$. We find that the periodic, adiabatic and reflexive bcs give closer approximation to the value of $T^H_c$ than fixed $+1$ and fixed $-1$ bcs with all three initial conditions for lattice size less than $50\times50$. However, for lattice size between $50\times50$ and $200\times200$, fixed $+1$ bc and fixed $-1$ bc give closer approximation to the $T^H_c$ with initial conditions in which all spins are in down configuration and all spins are in up configuration respectively.


title: Relative quantum cohomology

abstract: We establish a system of PDE, called open WDVV, that constrains the bulk-deformed superpotential and associated open Gromov-Witten invariants of a Lagrangian submanifold $L \subset X$ with a bounding chain. Simultaneously, we define the quantum cohomology algebra of $X$ relative to $L$ and prove its associativity. We also define the relative quantum connection and prove it is flat. A wall-crossing formula is derived that allows the interchange of point-like boundary constraints and certain interior constraints in open Gromov-Witten invariants. Another result is a vanishing theorem for open Gromov-Witten invariants of homologically non-trivial Lagrangians with more than one point-like boundary constraint. In this case, the open Gromov-Witten invariants with one point-like boundary constraint are shown to recover certain closed invariants. From open WDVV and the wall-crossing formula, a system of recursive relations is derived that entirely determines the open Gromov-Witten invariants of $(X,L) = (\mathbb{C}P^n, \mathbb{R}P^n)$ with $n$ odd, defined in previous work of the authors. Thus, we obtain explicit formulas for enumerative invariants defined using the Fukaya-Oh-Ohta-Ono theory of bounding chains.


title: High pressure melt locus of iron from atom-in-jellium calculations

abstract: Although usually considered as a technique for predicting electron states in dense plasmas, atom-in-jellium calculations can be used to predict the mean displacement of the ion from its equilibrium position in colder matter, as a function of compression and temperature. The Lindemann criterion of a critical displacement for melting can then be employed to predict the melt locus, normalizing for instance to the observed melt temperature or to more direct simulations such as molecular dynamics (MD). This approach reproduces the high pressure melting behavior of Al as calculated using the Lindemann model and thermal vibrations in the solid. Applied to Fe, we find that it reproduces the limited-range melt locus of a multiphase equation of state (EOS) and the results of ab initio MD simulations, and agrees less well with a Lindemann construction using an older EOS. The resulting melt locus lies significantly above the older melt locus for pressures above 1.5\,TPa, but is closer to recent ab initio MD results and extrapolations of an analytic fit to them. This study confirms the importance of core freezing in massive exoplanets, predicting that a slightly smaller range of exoplanets than previously assessed would be likely to exhibit dynamo generation of magnetic fields by convection in the liquid portion of the core.


title: Displacements representations for the problems with spherical and circular material surfaces with surface tension

abstract: The displacements representations of the type used by Christensen and Lo (1979) are modified to allow for analytical treatment of problems involving spherical and circular material surfaces that possess constant surface tension. The modified representations are used to derive closed-form expressions for the local elastic fields and effective moduli of a macroscopically isotropic composite materials containing spherical and circular inhomogeneities with the interfaces described by the complete Gurtin-Murdoch and Steigmann-Ogden models.


title: Table-Based Neural Units: Fully Quantizing Networks for Multiply-Free Inference

abstract: In this work, we propose to quantize all parts of standard classification networks and replace the activation-weight--multiply step with a simple table-based lookup. This approach results in networks that are free of floating-point operations and free of multiplications, suitable for direct FPGA and ASIC implementations. It also provides us with two simple measures of per-layer and network-wide compactness as well as insight into the distribution characteristics of activationoutput and weight values. We run controlled studies across different quantization schemes, both fixed and adaptive and, within the set of adaptive approaches, both parametric and model-free. We implement our approach to quantization with minimal, localized changes to the training process, allowing us to benefit from advances in training continuous-valued network architectures. We apply our approach successfully to AlexNet, ResNet, and MobileNet. We show results that are within 1.6% of the reported, non-quantized performance on MobileNet using only 40 entries in our table. This performance gap narrows to zero when we allow tables with 320 entries. Our results give the best accuracies among multiply-free networks.


title: Algebraic cycles and Verra fourfolds

abstract: This note is about the Chow ring of Verra fourfolds. For a general Verra fourfold, we show that the Chow group of homologically trivial $1$-cycles is generated by conics. We also show that Verra fourfolds admit a multiplicative Chow-K\"unneth decomposition, and draw some consequences for the intersection product in the Chow ring of Verra fourfolds.