title: Connecting and dissecting galaxies' angular momenta and neutral gas in a hierarchical universe: cue DARK SAGE

abstract: We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional disc structure, making both j_disc and q predicted quantities. We show that Dark Sage produces a clear correlation between gas fraction and j_disc at fixed disc mass, in line with recent results from observations and hydrodynamic simulations. This translates to a tight q--f_atm sequence for star-forming central galaxies, which closely tracks the analytic prediction of Obreschkow et al. The scatter in this sequence is driven by the probability distribution function of mass as a function of j (PDF of j) within discs, specifically where it peaks. We find that halo mass is primarily responsible for the peak location of the PDF of j, at least for low values of q. Two main mechanisms of equal significance are then identified for disconnecting f_atm from q. Mergers in the model can trigger quasar winds, with the potential to blow out most of the gas disc, while leaving the stellar disc relatively unharmed. Ram-pressure stripping of satellite galaxies has a similar effect, where f_atm can drop drastically with only a minimal effect to q. We highlight challenges associated with following these predictions up with observations.


title: Combined Explanations of the $b \to s \mu^+ \mu^-$ and $b \to c \tau^- {\bar\nu}$ Anomalies: a General Model Analysis

abstract: There are four models of tree-level new physics (NP) that can potentially simultaneously explain the $b \to s \mu^+ \mu^-$ and $b \to c \tau^- {\bar\nu}$ anomalies. They are the S3, U3 and U1 leptoquarks (LQs), and a triplet of SM-like vector bosons (VBs). Under the theoretical assumption that the NP couples predominantly to the third generation, previous analyses found that, when constraints from other processes are taken into account, the S3, U3 and VB models cannot explain the B anomalies, but U1 is viable. In this paper, we reanalyze these models, but without any assumption about their couplings. We find that, even in this most general case, S3 and U3 are excluded. For the U1 model, constraints from the semileptonic lepton-flavour-violating (LFV) processes $B \to K^{(*)} \mu^\pm \tau^\mp$, $\tau \to \mu \phi$ and $\Upsilon \to \mu \tau$, which have been largely ignored previously, are found to be very important. Because of the LFV constraints, the pattern of couplings of the U1 LQ is similar to that obtained with the above theoretical assumption. Also, the LFV constraints render unimportant those constraints obtained using the renormalization group equations. As for the VB model, it is excluded if the above theoretical assumption is made due to the additional constraints from $B^0_s$-${\bar B}^0_s$ mixing, $\tau\to 3\mu$ and $\tau \to \mu \nu {\bar\nu}$. By contrast, we find a different set of NP couplings that both explains the $b \to s \mu^+ \mu^-$ anomaly and is compatible with all constraints. However, it does not reproduce the measured values of the $b \to c \tau^- {\bar\nu}$ anomalies -- it would be viable only if future measurements find that the central values of these anomalies are reduced. Even so, this VB model is excluded by the LHC bounds on high-mass resonant dimuon pairs. This conclusion is reached without any assumptions about the NP couplings.


title: Approximating real-rooted and stable polynomials, with combinatorial applications

abstract: Let $p(x)=a_0 + a_1 x + \ldots + a_n x^n$ be a polynomial with all roots real and satisfying $x \leq -\delta$ for some $0<\delta <1$. We show that for any $0 < \epsilon <1$, the value of $p(1)$ is determined within relative error $\epsilon$ by the coefficients $a_k$ with $k \leq {c \over \sqrt{\delta}} \ln {n \over \epsilon \sqrt{ \delta}}$ for some absolute constant $c > 0$. Consequently, if $m_k(G)$ is the number of matchings with $k$ edges in a graph $G$, then for any $0 < \epsilon < 1$, the total number $M(G)=m_0(G)+m_1(G) + \ldots $ of matchings is determined within relative error $\epsilon$ by the numbers $m_k(G)$ with $k \leq c \sqrt{\Delta} \ln (v /\epsilon)$, where $\Delta$ is the largest degree of a vertex, $v$ is the number of vertices of $G$ and $c >0$ is an absolute constant. We prove a similar result for polynomials with complex roots satisfying $\Re\thinspace z \leq -\delta$ and apply it to estimate the number of unbranched subgraphs of $G$.


title: Shell-model coupled-cluster method for open-shell nuclei

abstract: We present an approach to derive effective shell-model interactions from microscopic nuclear forces. The similarity-transformed coupled-cluster Hamiltonian decouples the single-reference state of a closed-shell nucleus and provides us with a core for the shell model. We use a second similarity transformation to decouple a shell-model space from the excluded space. We show that the three-body terms induced by both similarity transformations are crucial for an accurate computation of ground and excited states. As a proof of principle we use a nucleon-nucleon interaction from chiral effective field theory, employ a $^4$He core, and compute low-lying states of $^{6-8}$He and $^{6-8}$Li in $p$-shell model spaces. Our results agree with benchmarks from full configuration interaction.


title: Contrastive Hebbian Learning with Random Feedback Weights

abstract: Neural networks are commonly trained to make predictions through learning algorithms. Contrastive Hebbian learning, which is a powerful rule inspired by gradient backpropagation, is based on Hebb's rule and the contrastive divergence algorithm. It operates in two phases, the forward (or free) phase, where the data are fed to the network, and a backward (or clamped) phase, where the target signals are clamped to the output layer of the network and the feedback signals are transformed through the transpose synaptic weight matrices. This implies symmetries at the synaptic level, for which there is no evidence in the brain. In this work, we propose a new variant of the algorithm, called random contrastive Hebbian learning, which does not rely on any synaptic weights symmetries. Instead, it uses random matrices to transform the feedback signals during the clamped phase, and the neural dynamics are described by first order non-linear differential equations. The algorithm is experimentally verified by solving a Boolean logic task, classification tasks (handwritten digits and letters), and an autoencoding task. This article also shows how the parameters affect learning, especially the random matrices. We use the pseudospectra analysis to investigate further how random matrices impact the learning process. Finally, we discuss the biological plausibility of the proposed algorithm, and how it can give rise to better computational models for learning.


title: Speaker Adapted Beamforming for Multi-Channel Automatic Speech Recognition

abstract: This paper presents, in the context of multi-channel ASR, a method to adapt a mask based, statistically optimal beamforming approach to a speaker of interest. The beamforming vector of the statistically optimal beamformer is computed by utilizing speech and noise masks, which are estimated by a neural network. The proposed adaptation approach is based on the integration of the beamformer, which includes the mask estimation network, and the acoustic model of the ASR system. This allows for the propagation of the training error, from the acoustic modeling cost function, all the way through the beamforming operation and through the mask estimation network. By using the results of a first pass recognition and by keeping all other parameters fixed, the mask estimation network can therefore be fine tuned by retraining. Utterances of a speaker of interest can thus be used in a two pass approach, to optimize the beamforming for the speech characteristics of that specific speaker. It is shown that this approach improves the ASR performance of a state-of-the-art multi-channel ASR system on the CHiME-4 data. Furthermore the effect of the adaptation on the estimated speech masks is discussed.


title: Segre Class Computation and Practical Applications

abstract: Let $X \subset Y$ be closed (possibly singular) subschemes of a smooth projective toric variety $T$. We show how to compute the Segre class $s(X,Y)$ as a class in the Chow group of $T$. Building on this, we give effective methods to compute intersection products in projective varieties, to determine algebraic multiplicity without working in local rings, and to test pairwise containment of subvarieties of $T$. Our methods may be implemented without using Groebner bases; in particular any algorithm to compute the number of solutions of a zero-dimensional polynomial system may be used.


title: Built-in Vulnerabilities to Imperceptible Adversarial Perturbations

abstract: Designing models that are robust to small adversarial perturbations of their inputs has proven remarkably difficult. In this work we show that the reverse problem---making models more vulnerable---is surprisingly easy. After presenting some proofs of concept on MNIST, we introduce a generic tilting attack that injects vulnerabilities into the linear layers of pre-trained networks by increasing their sensitivity to components of low variance in the training data without affecting their performance on test data. We illustrate this attack on a multilayer perceptron trained on SVHN and use it to design a stand-alone adversarial module which we call a steganogram decoder. Finally, we show on CIFAR-10 that a poisoning attack with a poisoning rate as low as 0.1% can induce vulnerabilities to chosen imperceptible backdoor signals in state-of-the-art networks. Beyond their practical implications, these different results shed new light on the nature of the adversarial example phenomenon.


title: Stochastic Geometry Analysis of Interference Coordination in NOMA-based Heterogeneous Networks

abstract: In this paper, we present a stochastic geometry-based analysis for the coverage probability of multi-tier heterogeneous networks employing interference coordination for non-orthogonal multiple access (NOMA). Two basic and conventional interference coordination alternatives; strict fractional frequency reuse (FFR) and soft frequency reuse (SFR) are analyzed, and their performance gains on different instances are illustrated. Coverage probability analysis has been extended to that of spectral efficiency for a fair performance evaluation between different FFR methods. The analysis covers both open-access networks with micro- and pico-cells, and closed-access networks with femto-cells. For heterogeneous networks, NOMA is observed to achieve better spectral efficiency than OFDMA, while OFDMA provides better coverage. Both analysis and numerical results reveal that strict FFR provides better coverage while SFR provides spectral efficiency advantage in open-access. However, strict FFR outperforms SFR for both measures in closed-access networks.


title: Desynchronization of Random Dynamical System under Perturbation by an Intrinsic Noise

abstract: In the theory of random dynamical systems (RDS), individuals with different initial states follow a same law of motion that is stochastically changing with time | called extrinsic noise. In the present work, intrin- sic noises for each individual are considered as a perturbation to an RDS. This gives rise to random Markov systems (RMS) in which the law of mo- tion is still stochastically changing with time, but individuals also exhibit statistically independent variations, with each transition having a small probability not to follow the law. As a consequence, two individuals in an RMS system go through stochastically distributed periods of synchro- nization and desynchronization, driven by extrinsic and intrinsic noises respectively. We show that in-sync time, e.g., escaping from a random attractor, has a symptotic geometric distribution.


title: Mu-tau reflection symmetry with a high scale texture-zero

abstract: The $\mu\tau$-reflection symmetric neutrino mass matrix can accommodate all known neutrino mixing angles, with maximal atmospheric angle fixed, and predicts all the unknown CP phases of the lepton sector but is unable to predict the absolute neutrino mass scale. Here we present a highly predictive scenario where $\mu\tau$-reflection is combined with a discrete abelian symmetry to enforce a texture-zero in the mass matrix of the heavy right-handed neutrinos that generate the light neutrino masses. Such a restriction reduces the free parameters of the low energy theory to zero and the absolute neutrino mass scale is restricted to few discrete regions, three in the few meV range and one extending up to around 30 meV. The heavy neutrino sector is dependent only on two free parameters which are further restricted to small regions from the requirement of successful leptogenesis. Mass degenerate heavy neutrinos are possible in one case but there is no resonant enhancement of the CP asymmetry.


title: Chaos in convolution operators on the space of entire functions of infinitely many complex variables

abstract: A classical result of Godefroy and Shapiro states that every nontrivial convolution operator on the space $\mathcal{H}(\mathbb{C}^n)$ of entire functions of several complex variables is hypercyclic. In sharp contrast with this result F\'avaro and Mujica show that no translation operator on the space $\mathcal{H}(\mathbb{C}^\mathbb{N})$ of entire functions of infinitely many complex variables is hypercyclic. In this work we study the linear dynamics of convolution operators on $\mathcal{H}(\mathbb{C}^\mathbb{N})$. First we show that no convolution operator on $\mathcal{H}(\mathbb{C}^\mathbb{N})$ is neither cyclic nor $n$-supercyclic for any positive integer $n$. After we study the notion of Li--Yorke chaos in non-metrizable topological vector spaces and we show that every nontrivial convolution operator on $\mathcal{H}(\mathbb{C}^\mathbb{N})$ is Li--Yorke chaotic.


title: Local fragmentation of thin disks in Eddington-inspired gravity

abstract: We find the generalized version of the Toomre's criterion for the stability of a rotating thin disk in the context of Eddington inspired Born-Infeld (EiBI) gravity which possesses one free parameter $\chi$. To do so we use the weak field limit of the theory and find the dispersion relation for the propagation of matter density waves on the surface of a self-gravitating and differentially rotating disk. Finally we find a new version of Toomre's stability criterion for thin disks. We show that EiBI gravity with negative $\chi$ destabilizes all the rotating thin disks. On the other hand EiBI with positive $\chi$ substantially can suppress the local fragmentation, and has stabilizing effects against axi-symmetric perturbations. More specifically, we show that only an annulus remains unstable on the surface of the disk. The width of the annulus directly depends on the magnitude of $\chi$.


title: Aggregation-Diffusion to Constrained Interaction: Minimizers & Gradient Flows in the Slow Diffusion Limit

abstract: Inspired by recent work on minimizers and gradient flows of constrained interaction energies, we prove that these energies arise as the slow diffusion limit of well-known aggregation-diffusion energies. We show that minimizers of aggregation-diffusion energies converge to a minimizer of the constrained interaction energy and gradient flows converge to a gradient flow. Our results apply to a range of interaction potentials, including singular attractive and repulsive-attractive power-law potentials. In the process of obtaining the slow diffusion limit, we also extend the well-posedness theory for aggregation-diffusion equations and Wasserstein gradient flows to admit a wide range of nonconvex interaction potentials. We conclude by applying our results to develop a numerical method for constrained interaction energies, which we use to investigate open questions on set valued minimizers.


title: Fast CapsNet for Lung Cancer Screening

abstract: Lung cancer is the leading cause of cancer-related deaths in the past several years. A major challenge in lung cancer screening is the detection of lung nodules from computed tomography (CT) scans. State-of-the-art approaches in automated lung nodule classification use deep convolutional neural networks (CNNs). However, these networks require a large number of training samples to generalize well. This paper investigates the use of capsule networks (CapsNets) as an alternative to CNNs. We show that CapsNets significantly outperforms CNNs when the number of training samples is small. To increase the computational efficiency, our paper proposes a consistent dynamic routing mechanism that results in $3\times$ speedup of CapsNet. Finally, we show that the original image reconstruction method of CapNets performs poorly on lung nodule data. We propose an efficient alternative, called convolutional decoder, that yields lower reconstruction error and higher classification accuracy.


title: Influence of Interfaces on the Transport Properties of Graphite revealed by Nanometer Thickness Reduction

abstract: We investigated the influence of thickness reduction on the transport properties of graphite microflakes. Using oxygen plasma etching we decreased the thickness of highly oriented pyrolytic graphite (HOPG) microflakes from $\sim 100$~nm to $\sim 20$~nm systematically. Keeping current and voltage electrodes intact, the electrical resistance $R(T)$, the magnetoresistance (MR) and Raman spectra were measured in every individual sample and after each etching step of a few nm. The results show that $R(T)$ and MR can increase or decrease with the sample thickness in a non-systematic way. The results indicate that HOPG samples are inhomogeneous materials, in agreement with scanning transmission electron microscopy images and X-ray diffraction data. Our results further indicate that the quantum oscillations in the MR are not an intrinsic property of the ideal graphite structure but their origin is related to internal conducting interfaces.


title: On Distributionally Robust Chance Constrained Programs with Wasserstein Distance

abstract: This paper studies a distributionally robust chance constrained program (DRCCP) with Wasserstein ambiguity set, where the uncertain constraints should be satisfied with a probability at least a given threshold for all the probability distributions of the uncertain parameters within a chosen Wasserstein distance from an empirical distribution. In this work, we investigate equivalent reformulations and approximations of such problems. We first show that a DRCCP can be reformulated as a conditional value-at-risk constrained optimization problem, and thus admits tight inner and outer approximations. We also show that a DRCCP of bounded feasible region is mixed integer representable by introducing big-M coefficients and additional binary variables. For a DRCCP with pure binary decision variables, by exploring the submodular structure, we show that it admits a big-M free formulation, which can be solved by a branch and cut algorithm. Finally, we present a numerical study to illustrate the effectiveness of the proposed formulations.


title: Automatic Design of Task-specific Robotic Arms

abstract: We present an interactive, computational design system for creating custom robotic arms given high-level task descriptions and environmental constraints. Various task requirements can be encoded as desired motion trajectories for the robot arm's end-effector. Given such end-effector trajectories, our system enables on-demand design of custom robot arms using a library of modular and reconfigurable parts such as actuators and connecting links. By searching through the combinatorial set of possible arrangements of these parts, our method generates a functional, as-simple-as-possible robot arm that is capable of tracking the desired trajectories. We demonstrate our system's capabilities by creating robot arm designs in simulation, for various trajectory following scenarios.


title: Improved Image Selection for Stack-Based HDR Imaging

abstract: Stack-based high dynamic range (HDR) imaging is a technique for achieving a larger dynamic range in an image by combining several low dynamic range images acquired at different exposures. Minimizing the set of images to combine, while ensuring that the resulting HDR image fully captures the scene's irradiance, is important to avoid long image acquisition and post-processing times. The problem of selecting the set of images has received much attention. However, existing methods either are not fully automatic, can be slow, or can fail to fully capture more challenging scenes. In this paper, we propose a fully automatic method for selecting the set of exposures to acquire that is both fast and more accurate. We show on an extensive set of benchmark scenes that our proposed method leads to improved HDR images as measured against ground truth using the mean squared error, a pixel-based metric, and a visible difference predictor and a quality score, both perception-based metrics.


title: RISE: Randomized Input Sampling for Explanation of Black-box Models

abstract: Deep neural networks are being used increasingly to automate data analysis and decision making, yet their decision-making process is largely unclear and is difficult to explain to the end users. In this paper, we address the problem of Explainable AI for deep neural networks that take images as input and output a class probability. We propose an approach called RISE that generates an importance map indicating how salient each pixel is for the model's prediction. In contrast to white-box approaches that estimate pixel importance using gradients or other internal network state, RISE works on black-box models. It estimates importance empirically by probing the model with randomly masked versions of the input image and obtaining the corresponding outputs. We compare our approach to state-of-the-art importance extraction methods using both an automatic deletion/insertion metric and a pointing metric based on human-annotated object segments. Extensive experiments on several benchmark datasets show that our approach matches or exceeds the performance of other methods, including white-box approaches. Project page: http://cs-people.bu.edu/vpetsiuk/rise/


title: Doubly Robust Estimation in Observational Studies with Partial Interference

abstract: Interference occurs when the treatment (or exposure) of one individual affects the outcomes of others. In some settings it may be reasonable to assume individuals can be partitioned into clusters such that there is no interference between individuals in different clusters, i.e., there is partial interference. In observational studies with partial interference, inverse probability weighted (IPW) estimators have been proposed of different possible treatment effects. However, the validity of IPW estimators depends on the propensity score being known or correctly modeled. Alternatively, one can estimate the treatment effect using an outcome regression model. In this paper, we propose doubly robust (DR) estimators which utilize both models and are consistent and asymptotically normal if either model, but not necessarily both, is correctly specified. Empirical results are presented to demonstrate the DR property of the proposed estimators, as well as the efficiency gain of DR over IPW estimators when both models are correctly specified. The different estimators are illustrated using data from a study examining the effects of cholera vaccination in Bangladesh.


title: "Zombie" or active? An alternative explanation to the properties of star-forming galaxies at high redshift

abstract: Star-forming galaxies at high redshift show anomalous values of infrared excess, which can be described only by extremizing the existing relations between the shape of their ultraviolet continuum emission and their infrared-to-ultraviolet luminosity ratio, or by constructing \textit{ad-hoc} models of star formation and dust distribution. We present an alternative explanation, based on unveiled AGN activity, to the existence of such galaxies. In fact, the presence of a weak AGN configures as a natural scenario in order to explain the observed spectral properties of such high-$z$ objects in terms of a continuum slope distribution rather than altered infrared excesses, due to the different shape of the AGN continuum emission with respect to quiescent galaxies. To this aim, we directly compare the infrared-to-ultraviolet properties of high-redshift galaxies to those of known categories of AGN (quasars and Seyferts). We also infer the characteristics of their possible X-ray emission. We find a strong similarity between the spectral shapes and luminosity ratios of AGN with the corresponding properties of such galaxies. In addition, we derive expected X-ray fluxes that are compatible with energetics from AGN activity. We conclude that a moderate AGN contribution to the UV emission of such high-$z$ objects is a valid alternative to explain their spectral properties. Even the presence of an active nucleus in each source would not violate the expected quasar statistics. Furthermore, we suggest that the observed similarities between anomalous star-forming galaxies and quasars may provide a benchmark for future theoretical and observational studies on the galaxy population in the early Universe.


title: COHERENT analysis of neutrino generalized interactions

abstract: Effective neutrino-quark generalized interactions are entirely determined by Lorentz invariance, so they include all possible four-fermion non derivative Lorentz structures. They contain neutrino-quark non-standard interactions as a subset, but span over a larger set that involves effective scalar, pseudoscalar, axial and tensor operators. Using recent COHERENT data, we derive constraints on the corresponding couplings by considering scalar, vector and tensor quark currents and assuming no lepton flavor dependence. We allow for mixed neutrino-quark Lorentz couplings and consider two types of scenarios in which: (i) one interaction at the nuclear level is present at a time, (ii) two interactions are simultaneously present. For scenarios (i) our findings show that scalar interactions are the most severely constrained, in particular for pseudoscalar-scalar neutrino-quark couplings. In contrast, tensor and non-standard vector interactions still enable for sizable effective parameters. We find as well that an extra vector interaction improves the data fit when compared with the result derived assuming only the standard model contribution. In scenarios (ii) the presence of two interactions relaxes the bounds and opens regions in parameter space that are otherwise closed, with the effect being more pronounced in the scalar-vector and scalar-tensor cases. We point out that barring the vector case, our results represent the most stringent bounds on effective neutrino-quark generalized interactions for mediator masses of order $\sim 1\,$GeV. They hold as well for larger mediator masses, case in which they should be compared with limits from neutrino deep-inelastic scattering data.


title: Scalar Tops and Perturbed Quadrupoles: Probing Fundamental Physics with Spin-Precessing Binaries

abstract: Parity violation in the gravitational interaction has an important impact on fundamental observables and the evolution of the universe. We here investigate for the first time our ability to probe the parity violating nature of the gravitational interaction using gravitational waves from spin-precessing binaries. Focusing on dynamical Chern-Simons gravity, we derive the spin-precession equations, calculate the gravitational waves emitted by spin-precessing, quasi-circular black hole binaries and estimate the level to which the theory could be constrained with future gravitational wave observations.


title: $T\bar{T}$ deformed partition functions

abstract: We demonstrate the presence of modular properties in partition functions of $T\bar{T}$ deformed conformal field theories. These properties are verified explicitly for the deformed free boson. The modular features facilitate a derivation of the asymptotic density of states in these theories, which turns out to interpolate between Cardy and Hagedorn behaviours. We also point out a sub-sector of the spectrum that remains undeformed under the $T\bar{T}$ flow. Finally, we comment on the deformation of the CFT vacuum character and its implications for the holographic dual.


title: Validation of Inventory models for Single-echelon Supply Chain using Discrete-event Simulation

abstract: Inventory decision frameworks proposed in the literature for single-echelon supply chain systems rely on assumptions to obtain closed form expressions. In particular, two such frameworks - one conventional and the other with a demand undershoot - determine optimal reorder point for a desired $\beta$ or Type-II service level. In this work we assess the accuracy and applicability of these frameworks with the help of a discrete event simulation framework developed in SimPy. While in several cases the closed form literature models under-predict the service level, and thus result in a higher reorder point and safety stock, we observe some situations where model predicted service levels match the simulated results with statistical significance.


title: Minimum Quadratic Helicity States

abstract: Building on previous results on the quadratic helicity in magnetohydrodynamics (MHD) we investigate particular minimum helicity states. Those are eigenfunctions of the curl operator and are shown to constitute solutions of the quasi-stationary incompressible ideal MHD equations. We then show that these states have indeed minimum quadratic helicity.


title: An affine generalization of evacuation

abstract: We establish the existence of an involution on tabloids that is analogous to Schutzenberger's evacuation map on standard Young tableaux. We find that the number of its fixed points is given by evaluating a certain Green's polynomial at $q = -1$, and satisfies a "domino-like" recurrence relation.


title: Direct measurement of the intra-pixel response function of Kepler Space Telescope's CCDs

abstract: Space missions designed for high precision photometric monitoring of stars often under-sample the point-spread function, with much of the light landing within a single pixel. Missions like MOST, Kepler, BRITE, and TESS, do this to avoid uncertainties due to pixel-to-pixel response nonuniformity. This approach has worked remarkably well. However, individual pixels also exhibit response nonuniformity. Typically, pixels are most sensitive near their centers and less sensitive near the edges, with a difference in response of as much as 50%. The exact shape of this fall-off, and its dependence on the wavelength of light, is the intra-pixel response function (IPRF). A direct measurement of the IPRF can be used to improve the photometric uncertainties, leading to improved photometry and astrometry of under-sampled systems. Using the spot-scan technique, we measured the IPRF of a flight spare e2v CCD90 imaging sensor, which is used in the Kepler focal plane. Our spot scanner generates spots with a full-width at half-maximum of $\lesssim$5 microns across the range of 400 nm - 900 nm. We find that Kepler's CCD shows similar IPRF behavior to other back-illuminated devices, with a decrease in responsivity near the edges of a pixel by $\sim$50%. The IPRF also depends on wavelength, exhibiting a large amount of diffusion at shorter wavelengths and becoming much more defined by the gate structure in the near-IR. This method can also be used to measure the IPRF of the CCDs used for TESS, which borrows much from the Kepler mission.


title: Spin-Precessing Black Hole Binaries in Dynamical Chern-Simons Gravity

abstract: Gravitational waves from spin-precessing binaries exhibit amplitude oscillations that provide an invaluable method to extract the spins of the inspiraling compact objects. The spin-spin and spin-orbit interactions that cause this effect are sensitive to the fundamental nature of gravity, which will allow us to constrain modified theories of gravity using gravitational wave observations of precessing binaries. We here consider precessing black hole binaries in dynamical Chern-Simons gravity, an effective theory of gravity that enhances parity violating interactions. We model the black holes as modified point particles using effective field theory, and derive the spin-precession equations for a binary system by working within the post-Newtonian formalism. We find that the spin-spin and quadrupole-monopole interactions of General Relativity are modified due to an interaction between the scalar dipoles of the two black holes and the modified black hole quadrupole as a result of the violation of the no hair theorems. These modifications enter the precession equations at leading post-Newtonian order. We further show that these precession equations admit seven constants of motion when neglecting radiation reaction, with only the mass-weighted effective spin being modified from General Relativity. We discuss how these may be used to reduce the precession equations to quadrature and the possibility of constructing analytic Fourier domain waveforms for generic spin-precessing binaries in dCS gravity.


title: Facial unique-maximum colorings of plane graphs with restriction on big vertices

abstract: A facial unique-maximum coloring of a plane graph is a proper coloring of the vertices using positive integers such that each face has a unique vertex that receives the maximum color in that face. Fabrici and G\"{o}ring (2016) proposed a strengthening of the Four Color Theorem conjecturing that all plane graphs have a facial unique-maximum coloring using four colors. This conjecture has been disproven for general plane graphs and it was shown that five colors suffice. In this paper we show that plane graphs, where vertices of degree at least four induce a star forest, are facially unique-maximum 4-colorable. This improves a previous result for subcubic plane graphs by Andova, Lidick\'y, Lu\v{z}ar, and \v{S}krekovski (2018). We conclude the paper by proposing some problems.


title: Law of two-sided exit by a spectrally positive strictly stable process

abstract: For a spectrally positive strictly stable process with index in (1,2), the paper obtains i) the density of the time when the process makes first exit from an interval by hitting the interval's lower end point before jumping over its upper end point, and ii) the joint distribution of the time, the undershoot, and the jump of the process when it makes first exit the other way around. For i), the density of the time of first exit is expressed as an infinite sum of functions, each the product of a polynomial and an exponential function, with all coefficients determined by the roots of a Mittag-Leffler function. For ii), conditional on the undershoot, the time and the jump of first exit are independent, and the marginal conditional densities of the time has similar features as i).


title: Griffiths Physics in an Ultracold Bose Gas

abstract: Coupled XY model systems consisting of three-dimensional (3D) systems with disordered interlayer physics are of significant theoretical interest. We realize a set of coupled quasi-2D layers of $^{87}$Rb in the presence of disordered inter-layer coupling. This is achieved with our high bandwidth arbitrary optical lattice to obviate restrictions on the dimensionality of disorder with speckle-generated optical fields. We identify phase crossover regions compatible with the existence of a pair of intermediate Griffiths phases between a thermal state and the emergence of bulk 3D superfluidity.


title: Simpler derivation of bounded pitch inequalities for set covering, and minimum knapsack sets

abstract: A valid inequality \alpha^Tx \ge \alpha_0 for a set covering problem is said to have pitch <= k ( a positive integer) if the k smallest positive \alpha_j sum to at least alpha_0. This paper presents a new, simple derivation of a relaxation for set covering problems whose solutions satisfy all valid inequalities of pitch and is of polynomial size, for each fixed . We also consider the minimum knapsack problem, and show that for each fixed integer p > 0 and 0 < \epsilon < 1 one can separate, within additive tolerance \epsilon, from the relaxation defined by the valid inequalities with coefficients in {0, 1, . . . , p} in time polynomial in the number of variables and 1/\epsilon.


title: Two Different Methods for Modelling the Likely Upper Economic Limit of the Future United Kingdom Wind Fleet

abstract: Methods for predicting the likely upper economic limit for the wind fleet in the United Kingdom should be simple to use whilst being able to cope with evolving technologies, costs and grid management strategies. This paper present two such models, both of which use data on historical wind patterns but apply different approaches to estimating the extent of wind shedding as a function of the size of the wind fleet. It is clear from the models that as the wind fleet increases in size, wind shedding will progressively increase, and as a result the overall economic efficiency of the wind fleet will be reduced. The models provide almost identical predictions of the efficiency loss and suggest that the future upper economic limit of the wind fleet will be mainly determined by the wind fleet Headroom, a concept described in some detail in the paper. The results, which should have general applicability, are presented in graphical form, and should obviate the need for further modelling using the primary data. The paper also discusses the effectiveness of the wind fleet in decarbonising the grid, and the growing competition between wind and solar fleets as sources of electrical energy for the United Kingdom.


title: Signal Processing Based Pile-up Compensation for Gated Single-Photon Avalanche Diodes

abstract: Single-photon avalanche diode (SPAD) based transient imaging suffers from an aberration called pile-up. When multiple photons arrive within a single repetition period of the illuminating laser, the SPAD records only the arrival of the first photon; this leads to a bias in the recorded light transient wherein the transient response at later time-instants are under-estimated. An unfortunate consequence of this is the need to operate the illumination at low-power levels to reduce the probability of multiple photons returning in a single period. Operating the laser at low power results in either low signal-to-noise ratio (SNR) in the measured transients or reduced frame rate due to longer exposure durations to achieve a high SNR. In this paper, we propose a signal processing-based approach to compensate pile-up in post-processing, thereby enabling high power operation of the illuminating laser. While increasing illumination does cause a fundamental information loss in the data captured by SPAD, we quantify this information loss using Cramer-Rao bound and show that the errors in our framework are only limited to this information loss. We experimentally validate our hypotheses using real data from a lab prototype.


title: Dilepton azimuthal correlations in $t \bar t$ production

abstract: The dilepton azimuthal correlation, namely the difference $\phi$ between the azimuthal angles of the positive and negative charged lepton in the laboratory frame, provides a stringent test of the spin correlation in $t \bar t$ production at the Large Hadron Collider. We introduce a parameterisation of the differential cross section $d\sigma / d\phi$ in terms of a Fourier series and show that the third-order expansion provides a sufficiently accurate approximation. This expansion can be considered as a `bridge' between theory and data, making it very simple to cast predictions in the Standard Model (SM) and beyond, and to report measurements, without the need to provide the numbers for the whole binned distribution. We show its application by giving predictions for the coefficients in the presence of (i) an anomalous top chromomagnetic dipole moment; (ii) an anomalous $tbW$ interaction. The methods presented greatly facilitate the study of this angular distribution, which is of special interest given the $3.2(3.7)\sigma$ deviation from the SM next-to-leading order prediction found by the ATLAS Collaboration in Run 2 data.


title: HybridNet: Integrating Model-based and Data-driven Learning to Predict Evolution of Dynamical Systems

abstract: The robotic systems continuously interact with complex dynamical systems in the physical world. Reliable predictions of spatiotemporal evolution of these dynamical systems, with limited knowledge of system dynamics, are crucial for autonomous operation. In this paper, we present HybridNet, a framework that integrates data-driven deep learning and model-driven computation to reliably predict spatiotemporal evolution of a dynamical systems even with in-exact knowledge of their parameters. A data-driven deep neural network (DNN) with Convolutional LSTM (ConvLSTM) as the backbone is employed to predict the time-varying evolution of the external forces/perturbations. On the other hand, the model-driven computation is performed using Cellular Neural Network (CeNN), a neuro-inspired algorithm to model dynamical systems defined by coupled partial differential equations (PDEs). CeNN converts the intricate numerical computation into a series of convolution operations, enabling a trainable PDE solver. With a feedback control loop, HybridNet can learn the physical parameters governing the system's dynamics in real-time, and accordingly adapt the computation models to enhance prediction accuracy for time-evolving dynamical systems. The experimental results on two dynamical systems, namely, heat convection-diffusion system, and fluid dynamical system, demonstrate that the HybridNet produces higher accuracy than the state-of-the-art deep learning based approach.


title: Kernel Methods for Nonlinear Connectivity Detection

abstract: In this paper, we show that the presence of nonlinear coupling between time series may be detected employing kernel feature space representations alone dispensing with the need to go back to solve the pre-image problem to gauge model adequacy. As a consequence, the canonical methodology for model construction, diagnostics, and Granger connectivity inference applies with no change other than computation using kernels in lieu of second-order moments.


title: Wall Stress Estimation of Cerebral Aneurysm based on Zernike Convolutional Neural Networks

abstract: Convolutional neural networks (ConvNets) have demonstrated an exceptional capacity to discern visual patterns from digital images and signals. Unfortunately, such powerful ConvNets do not generalize well to arbitrary-shaped manifolds, where data representation does not fit into a tensor-like grid. Hence, many fields of science and engineering, where data points possess some manifold structure, cannot enjoy the full benefits of the recent advances in ConvNets. The aneurysm wall stress estimation problem introduced in this paper is one of many such problems. The problem is well-known to be of a paramount clinical importance, but yet, traditional ConvNets cannot be applied due to the manifold structure of the data, neither does the state-of-the-art geometric ConvNets perform well. Motivated by this, we propose a new geometric ConvNet method named ZerNet, which builds upon our novel mathematical generalization of convolution and pooling operations on manifolds. Our study shows that the ZerNet outperforms the other state-of-the-art geometric ConvNets in terms of accuracy.


title: OctApps: a library of Octave functions for continuous gravitational-wave data analysis

abstract: Gravitational waves are minute ripples in spacetime, first predicted by Einstein's general theory of relativity in 1916. Gravitational waves from rapidly-rotating neutron stars, whose shape deviates from perfect axisymmetry, are a potential astrophysical source of gravitational waves, but which so far have not been detected. The search for this type of signals, also known as continuous waves, presents a significant data analysis challenge, as their weak signatures are expected to be buried deep within the instrumental noise of the LIGO and Virgo detectors. The OctApps library provides various functions, written in Octave, intended to aid research scientists who perform searches for continuous gravitational waves.


title: Time-Dependent Density Matrix Renormalization Group Algorithms for Nearly Exact Absorption and Fluorescence Spectra of Molecular Aggregates at Both Zero and Finite Temperature

abstract: We implement and apply time-dependent density matrix renormalization group (TD-DMRG) algorithms at zero and finite temperature to compute the linear absorption and fluorescence spectra of molecular aggregates. Our implementation is within a matrix product state/operator framework with an explicit treatment of the excitonic and vibrational degrees of freedom, and uses the locality of the Hamiltonian in the zero-exciton space to improve the efficiency and accuracy of the calculations. We demonstrate the power of the method by calculations on several molecular aggregate models, comparing our results against those from multi-layer multiconfiguration time- dependent Hartree and n-particle approximations. We find that TD-DMRG provides an accurate and efficient route to calculate the spectrum of molecular aggregates.


title: Entanglement entropy and $T \overline{T}$ deformation

abstract: Quantum gravity in a finite region of spacetime is conjectured to be dual to a conformal field theory deformed by the irrelevant operator $T \overline{T}$. We test this conjecture with entanglement entropy, which is sensitive to ultraviolet physics on the boundary while also probing the bulk geometry. We find that the entanglement entropy for an entangling surface consisting of two antipodal points on a sphere is finite and precisely matches the Ryu-Takayanagi formula applied to a finite region consistent with the conjecture of McGough, Mezei and Verlinde. We also consider a one-parameter family of conical entropies, which are finite and verify a conjecture due to Dong. Since ultraviolet divergences are local, we conclude that the $T \overline{T}$ deformation acts as an ultraviolet cutoff on the entanglement entropy. Our results support the conjecture that the $T \overline{T}$-deformed CFT is the holographic dual of a finite region of spacetime.


title: Confining Potentials

abstract: We use the t'Hoft-Wilson method for the generation of static fermions potential in order to derive a class of confining potentials which can describe the quark confinement. A general pattern for the construction of propagators through the localization of non-local actions is uncovered.


title: Qubit-nonlinear-oscillator systems: from the moderate-coupling limit to the ultrastrong-coupling regime

abstract: This work aims to provide an alternative approach to modeling a two-state system (qubit) coupled to a nonlinear oscillator. Within a single algebraic scheme based upon the f-deformed oscillator description, hard and soft nonlinearities are proposed to be simulated by making use of fitting algebraic models extracted from the trigonometric and modified P\"oschl-Teller potentials, respectively. In the regime where the strength of the coupling is considered to be moderate, this approach allows for an analytic, albeit approximate, diagonalization process of the proposed Hamiltonian through using the Van Vleck perturbation theory and embracing the two types of nonlinear features. In the ultrastrong-coupling regime, the effect of such nonlinearities upon the squeezing and phase space properties of the ground state of the composite system is also explored by numerical means.


title: CSI-based Outdoor Localization for Massive MIMO: Experiments with a Learning Approach

abstract: We report on experimental results on the use of a learning-based approach to infer the location of a mobile user of a cellular network within a cell, for a 5G-type Massive multiple input, multiple output (MIMO) system. We describe how the sample spatial covariance matrix computed from the CSI can be used as the input to a learning algorithm which attempts to relate it to user location. We discuss several learning approaches, and analyze in depth the application of extreme learning machines, for which theoretical approximate performance benchmarks are available, to the localization problem. We validate the proposed approach using experimental data collected on a Huawei 5G testbed, provide some performance and robustness benchmarks, and discuss practical issues related to the deployment of such a technique in 5G networks.


title: The squeezed thermal reservoir as a generalized equilibrium reservoir

abstract: We explore the perspective of considering the squeezed thermal reservoir as an equilibrium reservoir in a generalized Gibbs ensemble with two non-commuting conserved quantities. We outline the main properties of such a reservoir in terms of the exchange of energy, both heat and work, and entropy, giving some key examples to clarify its physical interpretation. This new paradigm allows for a correct and insightful interpretation of all thermodynamical features of the squeezed thermal reservoir, as well as other similar non-thermal reservoirs, including the characterization of reversibility and the first and second laws of thermodynamics.


title: Augmenting Source Code Lines with Sample Variable Values

abstract: Source code is inherently abstract, which makes it difficult to understand. Activities such as debugging can reveal concrete runtime details, including the values of variables. However, they require that a developer explicitly requests these data for a specific execution moment. We present a simple approach, RuntimeSamp, which collects sample variable values during normal executions of a program by a programmer. These values are then displayed in an ambient way at the end of each line in the source code editor. We discuss questions which should be answered for this approach to be usable in practice, such as how to efficiently record the values and when to display them. We provide partial answers to these questions and suggest future research directions.


title: Selective molecular transport in thermo-responsive polymer membranes: role of nanoscale hydration and fluctuations

abstract: For a wide range of modern soft functional materials the selective transport of sub-nanometer-sized molecules (`penetrants') through a stimuli-responsive polymeric membrane is key to the desired function. In this study, we investigate the diffusion properties of penetrants ranging from non-polar to polar molecules and ions in a matrix of collapsed Poly(N-isopropylacrylamide) (PNIPAM) polymers in water by means of extensive molecular dynamics simulations. We find that the water distributes heterogeneously in fractal-like cluster structures embedded in the nanometer-sized voids of the polymer matrix. The nano-clustered water acts as an important player in the penetrant diffusion, which proceeds via a hopping mechanism through `wet' transition states: the penetrants hop from one void to another via transient water channels opened by rare but decisive polymer fluctuations. The diffusivities of the studied penetrants extend over almost five orders of magnitude and thus enable a formulation of an analytical scaling relation with a clear non-Stokesian, exponential dependence of the diffusion coefficient on the penetrant's radius for the uncharged penetrants. Charged penetrants (ions) behave differently as they get captured in large isolated water clusters. Finally, we find large energetic activation barriers for hopping, which significantly depend on the hydration state and thereby challenge available transport theories.


title: Social media usage patterns during natural hazards

abstract: Natural hazards are becoming increasingly expensive as climate change and development are exposing communities to greater risks. Preparation and recovery are critical for climate change resilience, and social media are being used more and more to communicate before, during, and after disasters. While there is a growing body of research aimed at understanding how people use social media surrounding disaster events, most existing work has focused on a single disaster case study. In the present study, we analyze five of the costliest disasters in the last decade in the United States (Hurricanes Irene and Sandy, two sets of tornado outbreaks, and flooding in Louisiana) through the lens of Twitter. In particular, we explore the frequency of both generic and specific food-security related terms, and quantify the relationship between network size and Twitter activity during disasters. We find differences in tweet volume for keywords depending on disaster type, with people using Twitter more frequently in preparation for Hurricanes, and for real-time or recovery information for tornado and flooding events. Further, we find that people share a host of general disaster and specific preparation and recovery terms during these events. Finally, we find that among all account types, individuals with "average" sized networks are most likely to share information during these disasters, and in most cases, do so more frequently than normal. This suggests that around disasters, an ideal form of social contagion is being engaged in which average people rather than outsized influentials are key to communication. These results provide important context for the type of disaster information and target audiences that may be most useful for disaster communication during varying extreme events.


title: G51.04+0.07 and its environment: Identification of a new Galactic supernova remnant at low radio frequencies

abstract: We have identified a new supernova remnant (SNR), G51.04+0.07 (G51.04), using observations at 74 MHz from the VLA Low-Frequency Sky Survey Redux. Earlier, higher frequency radio continuum, recombination line, and infrared data had correctly inferred the presence of nonthermal radio emission within a large complex including ionised nebulae and active star formation. However, our observations have allowed us to redefine at least one SNR as a relatively small source (7'.5 x 3' in size) located at the southern periphery of the originally defined SNR candidate. The integrated flux density of G51.04 at 74 MHz is 6.1 +/- 0.8 Jy, while its radio continuum spectrum has a slope alpha=-0.52 +/- 0.05, typical of a shell-type remnant. We also measured spatial variations in the spectral index between 74 and 1400~MHz across the source, ranging from a steeper spectrum (alpha~-0.50) coincident with the brightest emission to a flatter component (alpha~-0.30) in the surrounding fainter region. To probe the interstellar medium into which the redefined SNR is likely evolving, we have analysed the surrounding atomic and molecular gas using the 21~cm neutral hydrogen (HI) and 13CO (J=1-0) emissions. We found that G51.04 is confined within an HI cavity and that its radio emission is consistent with the remains of a stellar explosion that occurred ~6300 yr ago at a distance of ~7.7 kpc. Kinematic data suggest that the newly discovered SNR lies in front of HII regions in the complex, consistent with the lack of a low frequency turnover in the continuum spectrum. The CO observations revealed molecular material that traces the central and northern parts of G51.04. The interaction between the cloud and the radio source is not conclusive. The relatively low flux density (~1.5 Jy at 1400 MHz) of G51.04 is consistent with many similar SNRs hidden along complex lines of sight towards inner Galactic emission complexes.


title: Quantitative MRI: Absolute T1, T2 and Proton Density Parameters from Deep Learning

abstract: Quantitative MRI is highly desirable in terms of intrinsic tissue parameters such as T1, T2 and proton density. This approach promises to minimize diagnostic variability and differentiate normal and pathological tissues by comparing tissue parameters to the normal ranges. Also, absolute quantification can help segment MRI tissue images with better accuracy compared to traditional qualitative segmentation methods. Currently, there are several methods proposed to quantify tissue parameters; however, all of them require excessive scan time and thus are difficult to be applied in clinical applications. In this paper, we propose a novel machine learning approach for MRI quantification, which can dramatically decrease the scan time and improve image quality.


title: Transition functions of diffusion processes with the Jack parameter on the Thoma simplex

abstract: The paper deals with a three-dimensional family of diffusion processes on an infinite-dimensional simplex. These processes were constructed by Borodin and Olshanski (arXiv:0706.1034; arXiv:0902.3395), and they include, as limit objects, the Ethier-Kurtz's infinitely-many-neutral-allels diffusion model (1981) and its extension found by Petrov (arXiv:0708.1930). Each process X from our family possesses a symmetrising measure M. Our main result is that the transition function of X has continuous density with respect to M. This is a generalization of earlier results due to Ethier (1992) and to Feng, Sun, Wang, and Xu (2011). Our proof substantially uses a special basis in the algebra of symmetric functions related to Laguerre polynomials.


title: Unravelling the stereodynamics of cold HD-H2 collisions

abstract: Measuring inelastic rates with partial wave resolution requires temperatures close to a Kelvin or below, even for the lightest molecule. In a recent experiment Perreault et al. [1] studied collisional relaxation of excited HD molecules in the v = 1, j = 2 state by para- and ortho-H2 at a temperature of about 1 K, extracting the angular distribution of scattered HD in the v = 1,j = 0 state. By state-preparation of the HD molecules, control of the angular distribution of scattered HD was demonstrated. Here, we report a first-principles simulation of that experiment which enables us to attribute the main features of the observed angular distribution to a single L = 2 partial-wave shape resonance. Our results demonstrate important stereodynamical insights that can be gained when numerically-exact quantum scattering calculations are combined with experimental results in the few-partial-wave regime.


title: Turbulence correction with artificial neural networks

abstract: We design an optical feedback network making use of machine learning techniques and demonstrate via simulations its ability to correct for the effects of turbulent propagation on optical modes. This artificial neural network scheme only relies on measuring the intensity profile of the distorted modes, making the approach simple and robust. The network results in the generation of various mode profiles at the transmitter that, after propagation through turbulence, closely resemble the desired target mode. The corrected optical mode profiles at the receiver are found to be nearly identical to the desired profiles, with near-zero mean square error indices. We are hopeful that the present results combining the fields of machine learning and optical communications will greatly enhance the robustness of free-space optical links.


title: Motion planning in high-dimensional spaces

abstract: Motion planning is a key tool that allows robots to navigate through an environment without collisions. The problem of robot motion planning has been studied in great detail over the last several decades, with researchers initially focusing on systems such as planar mobile robots and low degree-of-freedom (DOF) robotic arms. The increased use of high DOF robots that must perform tasks in real time in complex dynamic environments spurs the need for fast motion planning algorithms. In this overview, we discuss several types of strategies for motion planning in high dimensional spaces and dissect some of them, namely grid search based, sampling based and trajectory optimization based approaches. We compare them and outline their advantages and disadvantages, and finally, provide an insight into future research opportunities.


title: Photon statistics of amplified spontaneous emission

abstract: Developments in quantum technologies lead to new applications that require radiation sources with specific photon statistics. A widely used Poissonian statistics are easily produced by lasers; however, some applications require super- or sub-Poissonian statistics. Statistical properties of a light source are characterized by the second-order coherence function g^(2)(0). This function distinguishes stimulated radiation of lasers with g^(2)(0)=1 from light of other sources. For example, g^(2)(0)=2 for black-body radiation, and g^(2)(0)=0 for single-photon emission. One of the applications requiring super-Poissonian statistics (g^(2)(0)>1) is ghost imaging with thermal light. Ghost imaging also requires light with a narrow linewidth and high intensity. Currently, rather expensive and inefficient light sources are used for this purpose. In the last year, a superluminescent diode based on amplified spontaneous emission (ASE) has been considered as a new light source for ghost imaging. Even though ASE has been widely studied, its photon statistics has not been settled - there are neither reliable theoretical estimates of the second-order coherence function nor unambiguous experimental data. Our computer simulation clearly establishes that coherence properties of light produced by ASE are similar to that of a thermal source with g^(2)(0)=2 independent of pump power. This result manifests the fundamental difference between ASE and laser radiation.


title: The local principle of large deviations for compound Poisson process with catastrophes

abstract: The continuous time Markov process considered in this paper belongs to a class of population models with linear growth and catastrophes. There, the catastrophes happen at the arrival times of a Poisson process, and at each catastrophe time, a randomly selected portion of the population is eliminated. For this population process, we derive an asymptotic upper bound for the maximum value and prove the local large deviation principle.


title: Six line configurations and string dualities

abstract: We study the family of K3 surfaces of Picard rank sixteen associated with the double cover of the projective plane branched along the union of six lines, and the family of its Van Geemen-Sarti partners, i.e., K3 surfaces with special Nikulin involutions, such that quotienting by the involution and blowing up recovers the former. We prove that the family of Van Geemen-Sarti partners is a four-parameter family of K3 surfaces with $H \oplus E_7(-1) \oplus E_7(-1)$ lattice polarization. We describe explicit Weierstrass models on both families using even modular forms on the bounded symmetric domain of type $IV$. We also show that our construction provides a geometric interpretation, called geometric two-isogeny, for the F-theory/heterotic string duality in eight dimensions. As a result, we obtain novel F-theory models, dual to non-geometric heterotic string compactifications in eight dimensions with two non-vanishing Wilson line parameters.


title: Learning Treatment Regimens from Electronic Medical Records

abstract: Appropriate treatment regimens play a vital role in improving patient health status. Although some achievements have been made, few of the recent studies of learning treatment regimens have exploited different kinds of patient information due to the difficulty in adopting heterogeneous data to many data mining methods. Moreover, current studies seem too rigid with fixed intervals of treatment periods corresponding to the varying lengths of hospital stay. To this end, this work proposes a generic data-driven framework which can derive group-treatment regimens from electronic medical records by utilizing a mixed-variate restricted Boltzmann machine and incorporating medical domain knowledge. We conducted experiments on coronary artery disease as a case study. The obtained results show that the framework is promising and capable of assisting physicians in making clinical decisions.


title: Die Klassifikation der Gruppen bis zur Ordnung p^5

abstract: This paper presents a classification of all p-groups of order p^5 up to isomorphism. It contains a full list of their polycyclic presentations, a short introduction to the basic ideas of the methodes used to classify the groups, and a detailed proof of the adequacy and completeness of the classification and presentations presented in the paper.


title: The cosmological impact of future constraints on $H_0$ from gravitational-wave standard sirens

abstract: Gravitational-wave standard sirens present a novel approach for the determination of the Hubble constant. After the recent spectacular confirmation of the method thanks to GW170817 and its optical counterpart, additional standard siren measurements from future gravitational-wave sources are expected to constrain the Hubble constant to high accuracy. At the same time, improved constraints are expected from observations of cosmic microwave background (CMB) polarization and from baryon acoustic oscillations (BAO) surveys. We explore the role of future standard siren constraints on $H_0$ in light of expected CMB+BAO data. Considering a $10$-parameters cosmological model, in which curvature, the dark energy equation of state, and the Hubble constant are unbounded by CMB observations, we find that a combination of future CMB+BAO data will constrain the Hubble parameter to $\sim 1.5 \%$. Further extending the parameter space to a time-varying dark energy equation of state, we find that future CMB+BAO constraints on $H_0$ are relaxed to $\sim 3.0 \%$. These accuracies are within reach of future standard siren measurements from the Hanford-Livingston-Virgo and the Hanford-Livingston-Virgo-Japan-India networks of interferometers, showing the cosmological relevance of these sources. If future gravitational-wave standard siren measurements reach $1\%$ on $H_0$, as expected, they would significantly improve future CMB+BAO constraints on curvature and on the dark energy equation of state by up to a factor $\sim 3$. We also show that the inclusion of $H_0$ constraints from gravitational-wave standard sirens could result in a reduction of the dark energy figure-of-merit (i.e., the cosmological parameter volume) by up to a factor of $\sim 400$.


title: Exploring the Semantic Content of Unsupervised Graph Embeddings: An Empirical Study

abstract: Graph embeddings have become a key and widely used technique within the field of graph mining, proving to be successful across a broad range of domains including social, citation, transportation and biological. Graph embedding techniques aim to automatically create a low-dimensional representation of a given graph, which captures key structural elements in the resulting embedding space. However, to date, there has been little work exploring exactly which topological structures are being learned in the embeddings process. In this paper, we investigate if graph embeddings are approximating something analogous with traditional vertex level graph features. If such a relationship can be found, it could be used to provide a theoretical insight into how graph embedding approaches function. We perform this investigation by predicting known topological features, using supervised and unsupervised methods, directly from the embedding space. If a mapping between the embeddings and topological features can be found, then we argue that the structural information encapsulated by the features is represented in the embedding space. To explore this, we present extensive experimental evaluation from five state-of-the-art unsupervised graph embedding techniques, across a range of empirical graph datasets, measuring a selection of topological features. We demonstrate that several topological features are indeed being approximated by the embedding space, allowing key insight into how graph embeddings create good representations.


title: Fast, Anytime Motion Planning for Prehensile Manipulation in Clutter

abstract: Many methods have been developed for planning the motion of robotic arms for picking and placing, ranging from local optimization to global search techniques, which are effective for sparsely placed objects. Dense clutter, however, still adversely affects the success rate, computation times, and quality of solutions in many real-world setups. The current work integrates tools from existing methodologies and proposes a framework that achieves high success ratio in clutter with anytime performance. The idea is to first explore the lower dimensional end effector's task space efficiently by ignoring the arm, and build a discrete approximation of a navigation function, which guides the end effector towards the set of available grasps or object placements. This is performed online, without prior knowledge of the scene. Then, an informed sampling-based planner for the entire arm uses Jacobian-based steering to reach promising end effector poses given the task space guidance. While informed, the method is also comprehensive and allows the exploration of alternative paths over time if the task space guidance does not lead to a solution. This paper evaluates the proposed method against alternatives in picking or placing tasks among varying amounts of clutter for a variety of robotic manipulators with different end-effectors. The results suggest that the method reliably provides higher quality solution paths quicker, with a higher success rate relative to alternatives.


title: A unifying method for the design of algorithms canonizing combinatorial objects

abstract: We devise a unified framework for the design of canonization algorithms. Using hereditarily finite sets, we define a general notion of combinatorial objects that includes graphs, hypergraphs, relational structures, codes, permutation groups, tree decompositions, and so on. Our approach allows for a systematic transfer of the techniques that have been developed for isomorphism testing to canonization. We use it to design a canonization algorithm for combinatorial objects in general. This result gives new fastest canonization algorithms with an asymptotic running time matching the best known isomorphism algorithm for the following types of objects: hypergraphs, hypergraphs of bounded color class size, permutation groups (up to permutational isomorphism) and codes that are explicitly given (up to code equivalence).


title: HairNet: Single-View Hair Reconstruction using Convolutional Neural Networks

abstract: We introduce a deep learning-based method to generate full 3D hair geometry from an unconstrained image. Our method can recover local strand details and has real-time performance. State-of-the-art hair modeling techniques rely on large hairstyle collections for nearest neighbor retrieval and then perform ad-hoc refinement. Our deep learning approach, in contrast, is highly efficient in storage and can run 1000 times faster while generating hair with 30K strands. The convolutional neural network takes the 2D orientation field of a hair image as input and generates strand features that are evenly distributed on the parameterized 2D scalp. We introduce a collision loss to synthesize more plausible hairstyles, and the visibility of each strand is also used as a weight term to improve the reconstruction accuracy. The encoder-decoder architecture of our network naturally provides a compact and continuous representation for hairstyles, which allows us to interpolate naturally between hairstyles. We use a large set of rendered synthetic hair models to train our network. Our method scales to real images because an intermediate 2D orientation field, automatically calculated from the real image, factors out the difference between synthetic and real hairs. We demonstrate the effectiveness and robustness of our method on a wide range of challenging real Internet pictures and show reconstructed hair sequences from videos.


title: Degradation of Entanglement in Markovian Noise

abstract: The entanglement survival time is defined as the maximum time a system which is evolving under the action of local Markovian, homogenous in time noise, is capable to preserve the entanglement it had at the beginning of the temporal evolution. In this paper we study how this quantity is affected by the interplay between the coherent preserving and dissipative contributions of the corresponding dynamical generator. We report the presence of a counterintuitive, non-monotonic behaviour in such functional, capable of inducing sudden death of entanglement in models which, in the absence of unitary driving are capable to sustain entanglement for arbitrarily long times.


title: Micromechanical Experimental and Numerical Studies of Collagen Fibers Failure in Arterial Tissue

abstract: Arterial tissue failures lead to a number of clinical conditions that develop rapidly and unpredictably in vivo. Structural components and their interfacial mechanical strength of arterial tissue play a critical role in the process of arterail delamination. Therefore, there is a pressing need to understand the micromechanical mechanisms of arterial delamination. The objective of this study was to investigate various failure mechanisms (e.g. failure of collagen fibers) responsible for arterial interfacial delamination. In-situ tensile tests of fibers were performed on a micro-tester in the scanning electron microscope. A 3D unit cell model containing an individual fiber bridging two arterial tissue layers was constructed. An exponential cohesive zone model (CZM) was used to assess the stiffening and softening mechanical behaviors of collagen fiber bundles between the two arterial layers. An anisotropic constitutive model was implemented for characterizing the mechanical properties of the amorphous matrix which includes the fibrous cap and the underlying plaque tissue and a nonlinear elastic model was adopted for characterizing the mechanical properties of the fibers. The CZM and elastic parameter values of fibers were identified through an inverse boundary value approach that matches the load-displacement curves from simulation predictions of tensile test of collagen fibers with experimental measurements. The identified parameter values were then used as input in the 3D unit cell model, through which micromechanical factors affecting the resultant traction-separation relation for the unit cell were investigated via a parametric analysis. Results of the parametric analysis showed the applicability of the 3D unit cell model approach for evaluating the micromechanical mechanisms of arterial tissue failure processes.


title: Contrastive Explanations with Local Foil Trees

abstract: Recent advances in interpretable Machine Learning (iML) and eXplainable AI (XAI) construct explanations based on the importance of features in classification tasks. However, in a high-dimensional feature space this approach may become unfeasible without restraining the set of important features. We propose to utilize the human tendency to ask questions like "Why this output (the fact) instead of that output (the foil)?" to reduce the number of features to those that play a main role in the asked contrast. Our proposed method utilizes locally trained one-versus-all decision trees to identify the disjoint set of rules that causes the tree to classify data points as the foil and not as the fact. In this study we illustrate this approach on three benchmark classification tasks.


title: Quantum Field Theory with Classical Sources - Linearized Quantum Gravity

abstract: In a previous work and in terms of an exact quantum-mechanical framework, $\hbar$-independent causal and retarded expectation values of the second-quantized electro-magnetic fields in the Coulomb gauge were derived in the presence of a conserved classical electric current. The classical $\hbar$-independent Maxwell's equations then naturally emerged. In the present work, we extend these considerations to linear gravitational quantum deviations around a flat Minkowski space-time in a Coulomb-like gauge. The emergence of the classical causal and properly retarded linearized classical theory of general relativity with a conserved classical energy-momentum tensor is then outlined. The quantum-mechanical framework also provides for a simple approach to classical quadrupole gravitational radiation of Einstein and microscopic spontaneous graviton emission and/or absorption processes.


title: Hydroelastic wake on a thin elastic sheet floating on water

abstract: We investigate the hydroelastic wake created by a perturbation moving at constant speed along a thin elastic sheet floating at the surface of deep water. Using a high-resolution cross-correlation imaging technique, we characterize the waves as a function of the perturbation speed, for different sheet thicknesses. The general theoretical expression for the dispersion relation of hydroelastic waves includes three components: gravity, bending and stretching. The bending modulus and the tension in the sheet are independently measured. Excellent agreement is found between the experimental data and the theoretical expression.


title: An Input-Output Approach to Structured Stochastic Uncertainty

abstract: We consider linear time invariant systems with exogenous stochastic disturbances, and in feedback with structured stochastic uncertainties. This setting encompasses linear systems with both additive and multiplicative noise. Our concern is to characterize second-order properties such as mean-square stability and performance. A purely input-output treatment of these systems is given without recourse to state space models, and thus the results are applicable to certain classes of distributed systems. We derive necessary and sufficient conditions for mean-square stability in terms of the spectral radius of a linear matrix operator whose dimension is that of the number of uncertainties, rather than the dimension of any underlying state space models. Our condition is applicable to the case of correlated uncertainties, and reproduces earlier results for uncorrelated uncertainties. For cases where state space realizations are given, Linear Matrix Inequality (LMI) equivalents of the input-output conditions are given.


title: Tracing the Sources of Reionization in Cosmological Radiation Hydrodynamics Simulations

abstract: We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}<M<10^{10}{\rm M_{\odot}}h^{-1}$, stars with metallicity $10^{-3}Z_{\odot}<Z<10^{-1.5}Z_{\odot}$, and stars with age $3\,{\rm Myr} < t < 10 \, {\rm Myr}$ dominate reionization by both mass and volume. We show that the sources that reionize most of the volume of the Universe by $z=6$ are not necessarily the same sources that dominate the meta-galactic UV background at the same redshift. We further show that in our simulations, the contribution of each type of source to reionization is not uniform across different gas phases. The IGM, CGM, filaments, ISM, and rarefied supernova heated gas have all been photoionized by different classes of sources. Collisional ionisation contributes at both the lowest and highest densities. In the early stages of the formation of individual HII bubbles, reionization proceeds with the formation of concentric shells of gas ionised by different classes of sources, leading to large temperature variations as a function of galacto-centric radius. The temperature structure of individual HII bubbles may thus give insight into the star formation history of the galaxies acting as the first ionising sources. Our explorative simulations highlight how the complex nature of reionization can be better understood by using our photon tracer algorithm.


title: Learning from imperfections: constructing phase diagrams from atomic imaging of fluctuations

abstract: Materials characterization and property measurements are a cornerstone of material science, providing feedback from synthesis to applications. Traditionally, a single sample is used to derive information on a single point in composition space, and imperfections, impurities and stochastic details of material structure are deemed irrelevant or complicating factors in analysis. Here we demonstrate that atomic-scale studies of a single nominal composition can provide information on a finite area of chemical space. This information can be used to reconstruct the material properties in a finite composition and temperature range. We develop a statistical physics-based framework that incorporates chemical and structural data to infer effective atomic interactions driving segregation in a La5/8Ca3/8MnO3 thin-film. A variational autoencoder is used to determine anomalous behaviors in the composition phase diagram. This study provides a framework for creating generative models from diverse data and provides direct insight into the driving forces for cation segregation in manganites.


title: CommunityWatch: The Swiss-Army Knife of BGP Anomaly Detection

abstract: We present CommunityWatch, an open-source system that enables timely and accurate detection of BGP routing anomalies. CommunityWatch leverages meta-data encoded by AS operators on their advertised routes through the BGP Communities attribute. The BGP Communities values lack standardized semantics, offering the flexibility to attach a wide range of information, including AS relationships, location data, and route redistribution policies. Therefore, parsing and correlating Community values and their dynamics enables the detection and tracking of a variety of routing anomalies. We exhibit the efficacy of CommunityWatch through the detection of three different types of anomalies: infrastructure outages, route leaks, and traffic blackholing.


title: The Impact of Competition Between Cancer Cells and Healthy Cells on Optimal Drug Delivery

abstract: Cell competition is recognized to be instrumental to the dynamics and structure of the tumor-host interface in invasive cancers. In mild competition scenarios, the healthy tissue and cancer cells can coexist. When the competition is aggressive, competitive cells, the so called super-competitors, expand by killing other cells. Novel cytotoxic drugs and molecularly targeted drugs are commonly administered as part of cancer therapy. Both types of drugs are susceptible to various mechanisms of drug resistance, obstructing or preventing a successful outcome. In this paper, we develop a cancer growth model that accounts for the competition between cancer cells and healthy cells. The model incorporates resistance to both cytotoxic and targeted drugs. In both cases, the level of drug resistance is assumed to be a continuous variable ranging from fully-sensitive to fully-resistant. Using our model we demonstrate that when the competition is moderate, therapies using both drugs are more effective compared with single drug therapies. However, when cancer cells are highly competitive, targeted drugs become more effective. In this case, therapies that are initiated with a targeted drug and are exposed to it for a sufficiently long time are shown to have better outcomes. The results of the study stress the importance of adjusting the therapy to the pre-treatment resistance levels. We conclude with a study of the spatiotemporal propagation of drug resistance in a competitive setting, verifying that the same conclusions hold in the spatially heterogeneous case.


title: Test equations and linear stability of implicit-explicit general linear methods

abstract: Eigenvalue perturbation theory is applied to justify using complex-valued linear scalar test equations to characterize the stability of implicit-explicit general linear methods (IMEX GLMs) solving autonomous linear ordinary differential equations (ODEs) when the implicitly treated term is sufficiently stiff relative to the explicitly treated term. The stiff and non-stiff matrices are not assumed to be simultaneously diagonalizable or triangularizable and neither matrix is assumed to be symmetric or negative definite. The stability of IMEX GLMs solving complex-valued scalar linear ODEs displaying parabolic and hyperbolic stiffness is analyzed and related to the higher dimensional theory. The utility of the theoretical results is highlighted with a stability analysis of a family of IMEX Runge-Kutta methods solving IVPs of a linear 2D shallow-water model and a linear 1D advection-diffusion model.


title: Weight Thresholding on Complex Networks

abstract: Weight thresholding is a simple technique that aims at reducing the number of edges in weighted networks that are otherwise too dense for the application of standard graph theoretical methods. We show that the group structure of real weighted networks is very robust under weight thresholding, as it is maintained even when most of the edges are removed. This appears to be related to the correlation between topology and weight that characterizes real networks. On the other hand, the behavior of other properties is generally system dependent.


title: LazyFP: Leaking FPU Register State using Microarchitectural Side-Channels

abstract: Modern processors utilize an increasingly large register set to facilitate efficient floating point and SIMD computation. This large register set is a burden for operating systems, as its content needs to be saved and restored when the operating system context switches between tasks. As an optimization, the operating system can defer the context switch of the FPU and SIMD register set until the first instruction is executed that needs access to these registers. Meanwhile, the old content is left in place with the hope that the current task might not use these registers at all. This optimization is commonly called lazy FPU context switching. To make it possible, a processor offers the ability to toggle the availability of instructions utilizing floating point and SIMD registers. If the instructions are turned off, any attempt of executing them will generate a fault. In this paper, we present an attack that exploits lazy FPU context switching and allows an adversary to recover the FPU and SIMD register set of arbitrary processes or VMs. The attack works on processors that transiently execute FPU or SIMD instructions that follow an instruction generating the fault indicating the first use of FPU or SIMD instructions. On operating systems using lazy FPU context switching, the FPU and SIMD register content of other processes or virtual machines can then be reconstructed via cache side effects. With SIMD registers not only being used for cryptographic computation, but also increasingly for simple operations, such as copying memory, we argue that lazy FPU context switching is a dangerous optimization that needs to be turned off in all operating systems, if there is a chance that they run on affected processors.


title: Non-interacting fermions in hard-edge potentials

abstract: We consider the spatial quantum and thermal fluctuations of non-interacting Fermi gases of $N$ particles confined in $d$-dimensional non-smooth potentials. We first present a thorough study of the spherically symmetric pure hard-box potential, with vanishing potential inside the box, both at $T=0$ and $T>0$. We find that the correlations near the wall are described by a "hard edge" kernel, which depend both on $d$ and $T$, and which is different from the "soft edge" Airy kernel, and its higher $d$ generalizations, found for smooth potentials. We extend these results to the case where the potential is non-uniform inside the box, and find that there exists a family of kernels which interpolate between the above "hard edge" kernel and the "soft edge" kernels. Finally, we consider one-dimensional singular potentials of the form $V(x)\sim |x|^{-\gamma}$ with $\gamma>0$. We show that the correlations close to the singularity at $x=0$ are described by this "hard edge" kernel for $1\leq\gamma<2$ while they are described by a broader family of "hard edge" kernels known as the Bessel kernel for $\gamma=2$ and, finally by the Airy kernel for $\gamma>2$. These one-dimensional kernels also appear in random matrix theory, and we provide here the mapping between the $1d$ fermion models and the corresponding random matrix ensembles. Part of these results were announced in a recent Letter, EPL 120, 10006 (2017).

categories: cond-mat.stat-mech cond-mat.quant-gas math-ph math.MP


title: On the shape of invading population in oriented environments

abstract: We analyze the properties of population spreading in environments with spatial anisotropy within the frames of a lattice model of asymmetric (biased) random walkers. The expressions for the universal shape characteristics of the instantaneous configuration of population, such as asphericity $A$ and prolateness $S$ are found analytically and proved to be dependent only on the asymmetric transition probabilities in different directions. The model under consideration is shown to capture, in particular, the peculiarities of invasion in presence of an array of oriented tubes (fibers) in the environment.


title: Gradient Sensing via Cell Communication

abstract: Experimental evidence lends support to the conjecture that the ability of chains of cells to sense the gradient of an external chemical concentration could rely on cell-to-cell communication. This is the basis for the gradient sensing nature of a specific model type of the Local Excitation, Global Inhibition (LEGI) principle, wherein the strength of the external chemical field is sensed through a comparison between a local exciting species and a global inhibitor that is shared via intra-cellular reactions in the cell chain. In this study we generalize the nearest neighbor communication mechanism in the above-mentioned LEGI model in order to explore how the chemical sensing characteristics depend on the parameterization of the communication itself, cell size, and the radius of influence of neighboring cells. It was found that the radius of influence was less important than the approximating model for communication. Higher order approximations to the communication mechanism were better able to sense an external gradient. However, an analysis of the signal to noise ratio established that higher order models for communication were more prone to noise and thus have a lower signal to noise ratio. The generalization as well as the tools used in the analysis of the dynamics can be extended to more heterogeneous networks and can thus prove useful in using models and observations in the process of understanding chemical gradient via LEGI models with a communication component.


title: Adaptive Bayesian Estimation of Mixed Discrete-Continuous Distributions under Smoothness and Sparsity

abstract: We consider nonparametric estimation of a mixed discrete-continuous distribution under anisotropic smoothness conditions and possibly increasing number of support points for the discrete part of the distribution. For these settings, we derive lower bounds on the estimation rates in the total variation distance. Next, we consider a nonparametric mixture of normals model that uses continuous latent variables for the discrete part of the observations. We show that the posterior in this model contracts at rates that are equal to the derived lower bounds up to a log factor. Thus, Bayesian mixture of normals models can be used for optimal adaptive estimation of mixed discrete-continuous distributions.


title: Back and Forth Error Compensation and Correction Method for Linear Hyperbolic Systems with Application to the Maxwell's equations

abstract: We study the Back and Forth Error Compensation and Correction (BFECC) method for linear hyperbolic PDE systems. The BFECC method has been applied to schemes for advection equations to improve their stability and order of accuracy. Similar results are established in this paper for schemes for linear hyperbolic PDE systems with constant coefficients. We apply the BFECC method to central difference scheme and Lax-Friedrichs scheme for the Maxwell's equations and obtain second order accurate schemes with larger CFL number than the classical Yee scheme. The method is further applied to schemes on non-orthogonal unstructured grids. The new BFECC schemes for the Maxwell's equations operate on a single non-staggered grid and are simple to implement on unstructured grids. Numerical examples are given to demonstrate the effectiveness of the new schemes.


title: Standard Plane Detection in 3D Fetal Ultrasound Using an Iterative Transformation Network

abstract: Standard scan plane detection in fetal brain ultrasound (US) forms a crucial step in the assessment of fetal development. In clinical settings, this is done by manually manoeuvring a 2D probe to the desired scan plane. With the advent of 3D US, the entire fetal brain volume containing these standard planes can be easily acquired. However, manual standard plane identification in 3D volume is labour-intensive and requires expert knowledge of fetal anatomy. We propose a new Iterative Transformation Network (ITN) for the automatic detection of standard planes in 3D volumes. ITN uses a convolutional neural network to learn the relationship between a 2D plane image and the transformation parameters required to move that plane towards the location/orientation of the standard plane in the 3D volume. During inference, the current plane image is passed iteratively to the network until it converges to the standard plane location. We explore the effect of using different transformation representations as regression outputs of ITN. Under a multi-task learning framework, we introduce additional classification probability outputs to the network to act as confidence measures for the regressed transformation parameters in order to further improve the localisation accuracy. When evaluated on 72 US volumes of fetal brain, our method achieves an error of 3.83mm/12.7 degrees and 3.80mm/12.6 degrees for the transventricular and transcerebellar planes respectively and takes 0.46s per plane. Source code is publicly available at https://github.com/yuanwei1989/plane-detection.


title: Impact of distance determinations on Galactic structure. II. Old tracers

abstract: Here we review the efforts of a number of recent results that use old tracers to understand the build up of the Galaxy. Details that lead directly to using these old tracers to measure distances are discussed. We concentrate on the following: (1) the structure and evolution of the Galactic bulge and inner Galaxy constrained from the dynamics of individual stars residing therein; (2) the spatial structure of the old Galactic bulge through photometric observations of RR Lyrae-type stars; (3) the three\--dimensional structure, stellar density, mass, chemical composition, and age of the Milky Way bulge as traced by its old stellar populations; (4) an overview of RR Lyrae stars known in the ultra-faint dwarfs and their relation to the Galactic halo; and (5) different approaches for estimating absolute and relative cluster ages.


title: Isotropic tensor-valued polynomial functions of fourth-order tensors

abstract: Fourth-order tensor-valued functions appear in numerous fields of study. The formulation of practical models for these complex functions often requires their representation in terms of tensors of order two. In this paper, we develop an appropriate representation formula by assuming that the isotropic fourth--order tensor--valued function is a polynomial function in the components of two symmetric second-order tensors of degree $\leq$ 2. We illustrate the utility of the result by applying it to obtain a representation of the fluctuating velocity, pressure-gradient correlations of turbulence.


title: Towards the identification of Parkinson's Disease using only T1 MR Images

abstract: Parkinson's Disease (PD) is one of the most common types of neurological diseases caused by progressive degeneration of dopamin- ergic neurons in the brain. Even though there is no fixed cure for this neurodegenerative disease, earlier diagnosis followed by earlier treatment can help patients have a better quality of life. Magnetic Resonance Imag- ing (MRI) has been one of the most popular diagnostic tool in recent years because it avoids harmful radiations. In this paper, we investi- gate the plausibility of using MRIs for automatically diagnosing PD. Our proposed method has three main steps : 1) Preprocessing, 2) Fea- ture Extraction, and 3) Classification. The FreeSurfer library is used for the first and the second steps. For classification, three main types of classifiers, including Logistic Regression (LR), Random Forest (RF) and Support Vector Machine (SVM), are applied and their classification abil- ity is compared. The Parkinsons Progression Markers Initiative (PPMI) data set is used to evaluate the proposed method. The proposed system prove to be promising in assisting the diagnosis of PD.


title: Myocardial Segmentation of Contrast Echocardiograms Using Random Forests Guided by Shape Model

abstract: Myocardial Contrast Echocardiography (MCE) with micro-bubble contrast agent enables myocardial perfusion quantification which is invaluable for the early detection of coronary artery diseases. In this paper, we proposed a new segmentation method called Shape Model guided Random Forests (SMRF) for the analysis of MCE data. The proposed method utilizes a statistical shape model of the myocardium to guide the Random Forest (RF) segmentation in two ways. First, we introduce a novel Shape Model (SM) feature which captures the global structure and shape of the myocardium to produce a more accurate RF probability map. Second, the shape model is fitted to the RF probability map to further refine and constrain the final segmentation to plausible myocardial shapes. Evaluated on clinical MCE images from 15 patients, our method obtained promising results (Dice=0.81, Jaccard=0.70, MAD=1.68 mm, HD=6.53 mm) and showed a notable improvement in segmentation accuracy over the classic RF and its variants.