title: A Distributed Parallel Algorithm for Minimum Spanning Tree Problem

abstract: In this paper we present and evaluate a parallel algorithm for solving a minimum spanning tree (MST) problem for supercomputers with distributed memory. The algorithm relies on the relaxation of the message processing order requirement for one specific message type compared to the original GHS (Gallager, Humblet, Spira) algorithm. Our algorithm adopts hashing and message compression optimization techniques as well. To the best of our knowledge, this is the first parallel implementation of the GHS algorithm that linearly scales to more than 32 nodes (256 cores) of Infiniband cluster.


title: Multiple Emitters in a Waveguide: Non-Reciprocity and Correlated Photons at Perfect Elastic Transmission

abstract: We investigate interference and correlation effects when several detuned emitters are placed along a one-dimensional photonic waveguide. Such a setup allows multiple interactions between the photons and the strongly coupled emitters, and underlies proposed devices for quantum information processing. We show, first, that a pair of detuned two-level systems (2LS) separated by a half wavelength mimic a driven {\Lambda}-type three-level system (3LS) in both the single- and two- photon sectors. There is an interference-induced transparency peak at which the fluorescence is quenched, leaving the transmitted photons completely uncorrelated. Slightly away from this separation, we find that the inelastic scattering (fluorescence) is large, leading to nonlinear effects such as non-reciprocity (rectification). We connect this non-reciprocity to inelastic scattering caused by driving a dark pole and so derive a condition for maximum rectification. Finally, by placing a true 3LS midway between the two 2LS, we show that elastic scattering produces only transmission, but inelastic scattering nevertheless occurs (the fluorescence is not quenched) causing substantial photon correlations.


title: Deep Learning Ensembles for Melanoma Recognition in Dermoscopy Images

abstract: Melanoma is the deadliest form of skin cancer. While curable with early detection, only highly trained specialists are capable of accurately recognizing the disease. As expertise is in limited supply, automated systems capable of identifying disease could save lives, reduce unnecessary biopsies, and reduce costs. Toward this goal, we propose a system that combines recent developments in deep learning with established machine learning approaches, creating ensembles of methods that are capable of segmenting skin lesions, as well as analyzing the detected area and surrounding tissue for melanoma detection. The system is evaluated using the largest publicly available benchmark dataset of dermoscopic images, containing 900 training and 379 testing images. New state-of-the-art performance levels are demonstrated, leading to an improvement in the area under receiver operating characteristic curve of 7.5% (0.843 vs. 0.783), in average precision of 4% (0.649 vs. 0.624), and in specificity measured at the clinically relevant 95% sensitivity operating point 2.9 times higher than the previous state-of-the-art (36.8% specificity compared to 12.5%). Compared to the average of 8 expert dermatologists on a subset of 100 test images, the proposed system produces a higher accuracy (76% vs. 70.5%), and specificity (62% vs. 59%) evaluated at an equivalent sensitivity (82%).


title: Definition of fractional Laplacian for functions with polynomial growth\tha

abstract: We introduce a notion of fractional Laplacian for functions which grow more than linearly at infinity. In such case, the operator is not defined in the classical sense: nevertheless, we can give an ad-hoc definition which can be useful for applications in various fields, such as blowup and free boundary problems. In this setting, when the solution has a polynomial growth at infinity, the right hand side of the equation is not just a function, but an equivalence class of functions modulo polynomials of a fixed order. We also give a sharp version of the Schauder estimates in this framework, in which the full smooth H\"older norm of the solution is controlled in terms of the seminorm of the nonlinearity. Though the method presented is very general and potentially works for general nonlocal operators, for clarity and concreteness we focus here on the case of the fractional Laplacian.


title: Acoustic vibration problem for dissipative fluids

abstract: In this paper we analyze a finite element method for solving a quadratic eigenvalue problem derived from the acoustic vibration problem for a heterogeneous dissipative fluid. The problem is shown to be equivalent to the spectral problem for a noncompact operator and athorough spectral characterization is given. The numerical discretization of the problem is based on Raviart-Thomas finite elements. The method is proved to be free of spurious modes and to converge with optimal order. Finally, we report numerical tests which allow us to assess the performance of the method.


title: Quantum entanglement for two qubits in a nonstationary cavity

abstract: The quantum entanglement and the probability of the dynamical Lamb effect for two qubits caused by non-adiabatic fast change of the boundary conditions are studied. The conditional concurrence of the qubits for each fixed number of created photons in a nonstationary cavity is obtained as a measure of the dynamical quantum entanglement due to the dynamical Lamb effect. We discuss the physical realization of the dynamical Lamb effect, based on superconducting qubits.


title: High Power Polarized Positron Source

abstract: We discuss the basis of polarized positron production by low energy polarized electrons. Efficiency of conversion ~0.1-1% might be interesting for CEBAF facility and ILC as well.


title: Toric surfaces with equivariant Kazhdan-Lusztig atlases

abstract: A Kazhdan-Lusztig atlas, introduced by He, Knutson and Lu, on a stratified variety (V,Y) is a way of modeling the stratification Y of V locally on the stratification of Kazhdan-Lusztig varieties X^w_o \cap X_v by their intersection with opposite Schubert varieties X_u. We are interested in classifying smooth toric surfaces with Kazhdan-Lusztig atlases. This involves finding a degeneration of V to a union of Richardson varieties in the flag variety H/B_H of some Kac-Moody group H. We determine which toric surfaces have a chance at having a Kazhdan-Lusztig atlas by looking at their moment polytopes, then describe a way to find a suitable group H. More precisely, we find that (up to equivalence) there are 19 or 20 broken toric surfaces admitting simply-laced atlases, and that there are at most 7543 broken toric surfaces where H is any Kac-Moody group.


title: A Closed Form Solution to Multi-View Low-Rank Regression

abstract: Real life data often includes information from different channels. For example, in computer vision, we can describe an image using different image features, such as pixel intensity, color, HOG, GIST feature, SIFT features, etc.. These different aspects of the same objects are often called multi-view (or multi-modal) data. Low-rank regression model has been proved to be an effective learning mechanism by exploring the low-rank structure of real life data. But previous low-rank regression model only works on single view data. In this paper, we propose a multi-view low-rank regression model by imposing low-rank constraints on multi-view regression model. Most importantly, we provide a closed-form solution to the multi-view low-rank regression model. Extensive experiments on 4 multi-view datasets show that the multi-view low-rank regression model outperforms single-view regression model and reveals that multi-view low-rank structure is very helpful.


title: Szeg\H{o} kernel asymptotic expansion on CR manifolds with $S^1$ action

abstract: Let $X$ be a compact connected strongly pseudoconvex CR manifold of dimension $2n+1, n \ge 1$ with a transversal CR $S^1$ action on $X$. We establish an asymptotic expansion for the $m$-th Fourier component of the Szeg\H{o} kernel function as $m\rightarrow\infty$, where the expansion involves a contribution in terms of a distance function from lower dimensional strata of the $S^1$ action. We also obtain explicit formulas for the first three coefficients of the expansion.


title: New Hardness Results for the Permanent Using Linear Optics

abstract: In 2011, Aaronson gave a striking proof, based on quantum linear optics, showing that the problem of computing the permanent of a matrix is #P-hard. Aaronson's proof led naturally to hardness of approximation results for the permanent, and it was arguably simpler than Valiant's seminal proof of the same fact in 1979. Nevertheless, it did not prove that computing the permanent was #P-hard for any class of matrices which was not previously known. In this paper, we present a collection of new results about matrix permanents that are derived primarily via these linear optical techniques. First, we show that the problem of computing the permanent of a real orthogonal matrix is #P-hard. Much like Aaronson's original proof, this will show that even a multiplicative approximation remains #P-hard to compute. The hardness result even translates to permanents over finite fields, where the problem of computing the permanent of an orthogonal matrix is ModpP-hard in the finite field F_{p^4} for all primes p not equal to 2 or 3. Interestingly, this characterization is tight: in fields of characteristic 2, the permanent coincides with the determinant; in fields of characteristic 3, one can efficiently compute the permanent of an orthogonal matrix by a nontrivial result of Kogan. Finally, we use more elementary arguments to prove #P-hardness for the permanent of a positive semidefinite matrix, which shows that certain probabilities of boson sampling experiments with thermal states are hard to compute exactly despite the fact that they can be efficiently sampled by a classical computer.


title: Measurement of the optical dielectric function of transition metal dichalcogenide monolayers: MoS2, MoSe2, WS2 and WSe2

abstract: We report a determination of the complex in-plane dielectric function of monolayers of four transition metal dichalcogenides: MoS2, MoSe2, WS2 and WSe2, for photon energies from 1.5 - 3 eV. The results were obtained from reflection spectra using a Kramers-Kronig constrained variational analysis. From the dielectric functions, we obtain the absolute absorbance of the monolayers. We also provide a comparison of the dielectric function for the monolayers with the corresponding bulk materials.


title: A non-backtracking Polya's theorem

abstract: P\'olya's random walk theorem states that a random walk on a $d$-dimensional grid is recurrent for $d=1,2$ and transient for $d\ge3$. We prove a version of P\'olya's random walk theorem for non-backtracking random walks. Namely, we prove that a non-backtracking random walk on a $d$-dimensional grid is recurrent for $d=2$ and transient for $d=1$, $d\ge3$. Along the way, we prove several useful general facts about non-backtracking random walks on graphs. In addition, our proof includes an exact enumeration of the number of closed non-backtracking random walks on an infinite 2-dimensional grid. This enumeration suggests an interesting combinatorial link between non-backtracking random walks on grids, and trinomial coefficients.


title: Road Curb Extraction from Mobile LiDAR Point Clouds

abstract: Automatic extraction of road curbs from uneven, unorganized, noisy and massive 3D point clouds is a challenging task. Existing methods often project 3D point clouds onto 2D planes to extract curbs. However, the projection causes loss of 3D information which degrades the performance of the detection. This paper presents a robust, accurate and efficient method to extract road curbs from 3D mobile LiDAR point clouds. Our method consists of two steps: 1) extracting the candidate points of curbs based on the proposed novel energy function and 2) refining the candidate points using the proposed least cost path model. We evaluated our method on a large-scale of residential area (16.7GB, 300 million points) and an urban area (1.07GB, 20 million points) mobile LiDAR point clouds. Results indicate that the proposed method is superior to the state-of-the-art methods in terms of robustness, accuracy and efficiency. The proposed curb extraction method achieved a completeness of 78.62% and a correctness of 83.29%. These experiments demonstrate that the proposed method is a promising solution to extract road curbs from mobile LiDAR point clouds.


title: Finite-sum Composition Optimization via Variance Reduced Gradient Descent

abstract: The stochastic composition optimization proposed recently by Wang et al. [2014] minimizes the objective with the compositional expectation form: $\min_x~(\mathbb{E}_iF_i \circ \mathbb{E}_j G_j)(x).$ It summarizes many important applications in machine learning, statistics, and finance. In this paper, we consider the finite-sum scenario for composition optimization: \[\min_x f (x) := \frac{1}{n} \sum_{i = 1}^n F_i \left(\frac{1}{m} \sum_{j = 1}^m G_j (x) \right). \] We propose two algorithms to solve this problem by combining the stochastic compositional gradient descent (SCGD) and the stochastic variance reduced gradient (SVRG) technique. A constant linear convergence rate is proved for strongly convex optimization, which substantially improves the sublinear rate $O(K^{-0.8})$ of the best known algorithm.


title: On nodes of small degrees and degree profile in preferential dynamic attachment circuits

abstract: We investigate the joint distribution of nodes of small degrees and the degree profile in preferential dynamic attachment circuits. In particular, we study the joint asymptotic distribution of the number of the nodes of outdegree $0$ (terminal nodes) and outdegree $1$ in a very large circuit. The expectation and variance of the number of those two types of nodes are both asymptotically linear with respect to the age of the circuit. We show that the numbers of nodes of outdegree $0$ and $1$ asymptotically follow a two-dimensional Gaussian law via multivariate martingale methods. We also study the exact distribution of the degree of a node, as the circuit ages, via a series of P\'olya-Eggenberger urn models with "hiccups" in between. The exact expectation and variance of the degree of nodes are determined by recurrence methods. Phase transitions of these degrees are discussed briefly. This is an extension of the abstract [20].


title: Analysis and Modeling of Social Influence in High Performance Computing Workloads

abstract: Social influence among users (e.g., collaboration on a project) creates bursty behavior in the underlying high performance computing (HPC) workloads. Using representative HPC and cluster workload logs, this paper identifies, analyzes, and quantifies the level of social influence across HPC users. We show the existence of a social graph that is characterized by a pattern of dominant users and followers. This pattern also follows a power-law distribution, which is consistent with those observed in mainstream social networks. Given its potential impact on HPC workloads prediction and scheduling, we propose a fast-converging, computationally-efficient online learning algorithm for identifying social groups. Extensive evaluation shows that our online algorithm can (1) quickly identify the social relationships by using a small portion of incoming jobs and (2) can efficiently track group evolution over time.


title: Correcting Type Ia Supernova Distances for Selection Biases and Contamination in Photometrically Identified Samples

abstract: We present a new technique to create a bin-averaged Hubble Diagram (HD) from photometrically identified SN~Ia data. The resulting HD is corrected for selection biases and contamination from core collapse (CC) SNe, and can be used to infer cosmological parameters. This method, called "BBC" (BEAMS with Bias Corrections), includes two fitting stages. The first BBC fitting stage uses a posterior distribution that includes multiple SN likelihoods, a Monte Carlo simulation to bias-correct the fitted SALT-II parameters, and CC probabilities determined from a machine learning technique. The BBC fit determines 1) a bin-averaged HD (average distance vs. redshift), and 2) the nuisance parameters alpha and beta, which multiply the stretch and color (respectively) to standardize the SN brightness. In the second stage, the bin-averaged HD is fit to a cosmological model where priors can be imposed. We perform high precision tests of the BBC method by simulating large (150,000 event) data samples corresponding to the Dark Energy Survey Supernova Program. Our tests include three models of intrinsic scatter, each with two different CC rates. In the BBC fit, the SALT-II nuisance parameters alpha and beta are recovered to within 1% of their true values. In the cosmology fit, we determine the dark energy equation of state parameter w using a fixed value of Omega_matter as a prior: averaging over all six tests based on 6 x 150,000 = 900,000 SNe, there is a small w-bias of 0.006 +- 0.002. Finally, the BBC fitting code is publicly available in the SNANA package.


title: A spacetime DPG method for the Schrodinger equation

abstract: A spacetime Discontinuous Petrov Galerkin (DPG) method for the linear time-dependent Schrodinger equation is proposed. The spacetime approach is particularly attractive for capturing irregular solutions. Motivated by the fact that some irregular Schrodinger solutions cannot be solutions of certain first order reformulations, the proposed spacetime method uses the second order Schrodinger operator. Two variational formulations are proved to be well-posed: a strong formulation (with no relaxation of the original equation) and a weak formulation (also called the ultraweak formulation, that transfers all derivatives onto test functions). The convergence of the DPG method based on the ultraweak formulation is investigated using an interpolation operator. A standalone appendix analyzes the ultraweak formulation for general differential operators. Reports of numerical experiments motivated by pulse propagation in dispersive optical fibers are also included.


title: Symmetry protected line nodes in non-symmorphic magnetic space groups: Applications to UCoGe and UPd$_2$Al$_3$

abstract: We present the group-thoretical classification of gap functions in superconductors coexisting with some magnetic order in non-symmorphic magnetic space groups. Based on the weak-coupling BCS theory, we show that UCoGe-type ferromagnetic superconductors must have horizontal line nodes on either $k_z=0$ or $\pm\pi/c$ plane. Moreover, it is likely that additional Weyl point nodes exist at the axial point. On the other hand, in UPd$_2$Al$_3$-type antiferromagnetic superconductors, gap functions with $A_g$ symmetry possess horizontal line nodes in antiferromagnetic Brillouin zone boundary perpendicular to $c$-axis. In other words, the conventional fully-gapped $s$-wave superconductivity is forbidden in this type of antiferromagnetic superconductors, irrelevant to the pairing mechanism, as long as the Fermi surface crosses a zone boundary. UCoGe and UPd$_2$Al$_3$ are candidates for unconventional superconductors possessing hidden symmetry-protected line nodes, peculiar to non-symmorphic magnetic space groups.


title: How efficiently can one untangle a double-twist? Waving is believing!

abstract: It has long been known to mathematicians and physicists that while a full rotation in three-dimensional Euclidean space causes tangling, two rotations can be untangled. Formally, an untangling is a based nullhomotopy of the double-twist loop in the special orthogonal group of rotations. We study a particularly simple, geometrically defined untangling procedure, leading to new conclusions regarding the minimum possible complexity of untanglings. We animate and analyze how our untangling operates on frames in 3-space, and teach readers in a video how to wave the nullhomotopy with their hands.


title: Distributed Optimal Gas-Power Flow Using Convex Optimization and ADMM

abstract: This paper proposes a convex optimization based distributed algorithm to solve multi-period optimal gas-power flow (OGPF) in coupled energy distribution systems. At the gas distribution system side, the non-convex Weymouth gas flow equations is convexified as quadratic constraints. The optimal gas flow (OGF) subproblem is solved by an iterative second-order cone programming procedure, whose efficiency is two orders of magnitudes higher than traditional nonlinear methods. A convex quadratic program based initiation scheme is suggested, which helps to find a high-quality starting point. At the power distribution system side, convex relaxation is performed on the non-convex branch flow equations, and the optimal power flow (OPF) subproblem gives rise to a second order cone program. Tightness is guaranteed by the radial topology. In the proposed distributed algorithm, OGF and OPF are solved independently, and coordinated by the alternating direction multiplier method (ADMM). Numerical results corroborate significant enhancements on computational robustness and efficiency compared with existing centralized OGPF methods.


title: Tests of Lorentz symmetry in the gravitational sector

abstract: Lorentz symmetry is one of the pillars of both General Relativity and the Standard Model of particle physics. Motivated by ideas about quantum gravity, unification theories and violations of CPT symmetry, a significant effort has been put the last decades into testing Lorentz symmetry. This review focuses on Lorentz symmetry tests performed in the gravitational sector. We briefly review the basics of the pure gravitational sector of the Standard-Model Extension (SME) framework, a formalism developed in order to systematically parametrize hypothetical violations of the Lorentz invariance. Furthermore, we discuss the latest constraints obtained within this formalism including analyses of the following measurements: atomic gravimetry, Lunar Laser Ranging, Very Long Baseline Interferometry, planetary ephemerides, Gravity Probe B, binary pulsars, high energy cosmic rays,... In addition, we propose a combined analysis of all these results. We also discuss possible improvements on current analyses and present some sensitivity analyses for future observations.


title: Escape of mass and entropy for geodesic flows

abstract: In this paper we study the ergodic theory of the geodesic flow on negatively curved geometrically finite manifolds. We prove that the measure theoretic entropy is upper semicontinuous when there is no loss of mass. In case we are losing mass, the critical exponents of parabolic subgroups of the fundamental group have a significant meaning. More precisely, the failure of upper-semicontinuity of the entropy is determinated by the maximal parabolic critical exponent. We also study the pressure of positive H\"older continuous potentials going to zero through the cusps. We prove that the pressure map $t\mapsto P(tF)$ is differentiable until it undergoes a phase transition, after which it becomes constant. This description allows, in particular, to compute the entropy at infinity of the geodesic flow.


title: A study of adatom ripening on an Al (111) surface with machine learning force fields

abstract: Surface phenomena are increasingly becoming important in exploring nanoscale materials growth and characterization. Consequently, the need for atomistic based simulations is increasing. Nevertheless, relying entirely on quantum mechanical methods limits the length and time scales one can consider, resulting in an ever increasing dependence on alternative machine learning based force fields. Recently, we proposed a machine learning approach, known as AGNI, that allows fast and accurate atomic force predictions given the atom's neighborhood environment. Here, we make use of such force fields to study and characterize the nanoscale diffusion and growth processes occurring on an Al (111) surface. In particular we focus on the adatom ripening phenomena, confirming past experimental findings, wherein a low and high temperature growth regime were observed, using entirely molecular dynamics simulations.


title: On the Diversity of Compact Objects within Supernova Remnants II: Energy Loss Mechanisms

abstract: Energy losses from isolated neutron stars are commonly attributed to the emission of electromagnetic radiation from a rotating point-like magnetic dipole in vacuum. This emission mechanism predicts a braking index $n=3$, which is not observed in highly magnetized neutron stars. Despite this fact, the assumptions of a dipole field and rapid early rotation are often assumed a priori, typically causing a discrepancy between the characteristic age and the associated supernova remnant (SNR) age. We focus on neutron stars with `anomalous' magnetic fields that have established SNR associations and known ages. Anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) are usually described in terms of the magnetar model, which posits a large magnetic field established by dynamo action. The high magnetic field pulsars (HBPs) have extremely large magnetic fields just above QED scale (but below that of the AXPs and SGRs), and central compact objects (CCOs) may have buried fields that will emerge in the future as nascent magnetars. In the first part of this series we examined magnetic field growth as a method of uniting the CCOs with HBPs and X-ray dim isolated neutron stars (XDINSs) through evolution. In this work we constrain the characteristic age of these neutron stars using the related SNR age for a variety of energy loss mechanisms and allowing for arbitrary initial spin periods. In addition to the SNR age, we also make use of the observed braking indices and X-ray luminosities to constrain the models.


title: On the Stability of Kalman-Bucy Diffusion Processes

abstract: The Kalman-Bucy filter is the optimal state estimator for an Ornstein-Uhlenbeck diffusion given that the system is partially observed via a linear diffusion-type (noisy) sensor. Under Gaussian assumptions, it provides a finite-dimensional exact implementation of the optimal Bayes filter. It is generally the only such finite-dimensional exact instance of the Bayes filter for continuous state-space models. Consequently, this filter has been studied extensively in the literature since the seminal 1961 paper of Kalman and Bucy. The purpose of this work is to review, re-prove and refine existing results concerning the dynamical properties of the Kalman-Bucy filter so far as they pertain to filter stability and convergence. The associated differential matrix Riccati equation is a focal point of this study with a number of bounds, convergence, and eigenvalue inequalities rigorously proven. New results are also given in the form of exponential and comparison inequalities for both the filter and the Riccati flow.


title: Translation and Rotation of Transformation Media under Electromagnetic Pulse

abstract: It is well known that optical media create artificial geometry for light, and curved geometry acts as an effective optical medium. This correspondence originates from the form invariance of Maxwells equations, which recently has spawned a booming field called transformation optics. Here we investigate responses of three transformation media under electromagnetic pulses, and find that pulse radiation can induce unbalanced net force on transformation media, which will cause translation and rotation of transformation media although their final momentum can still be zero. Therefore, the transformation media do not necessarily stay the same after an electromagnetic wave passes through.


title: ExpressPass: End-to-End Credit-based Congestion Control for Datacenters

abstract: As link speeds increase in datacenter networks, existing congestion control algorithms become less effective in providing fast convergence. TCP-based algorithms that probe for bandwidth take a long time to reach the fair-share and lead to long flow completion times. An ideal congestion control algorithms for datacenter must provide 1) zero data loss, 2) fast convergence, and 3) low buffer occupancy. However, these requirements present conflicting goals. For fast convergence,flows must ramp up quickly, but this risks packet losses and large queues. Thus, even the state-of-the-art algorithms, such as TIMELY and DCQCN, rely on link layer flow control (e.g.,Priority-based Flow Control) to achieve zero loss. This paper presents a new approach, called ExpressPass, an end-to-end credit-based congestion control algorithm for datacenters. ExpressPass is inspired by credit-based flow control, but extends it to work end-to-end. The switches control the amount of credit packets by rate limiting and ensure data packets flow in the reverse direction without any loss. ExpressPass leverages this to ramp up aggressively. ExpressPass converges up to 80 times faster than DCTCP at 10Gbps link, and the gap increases as link speeds become faster. Our simulation with realistic workload shows that ExpressPass significantly reduces the flow completion time especially for small and medium size flows compared to DCTCP, HULL, and DX.


title: Discovery of RR Lyrae Stars in the Nuclear Bulge of the Milky Way

abstract: Galactic nuclei, like the one of the Milky Way, are extreme places with high stellar densities and, in most cases, hosting a supermassive black hole. One of the scenarios proposed for the formation of the Galactic nucleus is by merging of primordial globular clusters (Capuzzo-Dolcetta 1993). An implication of this model is that this region should host stars characteristically found in old Milky Way globular clusters. RR Lyrae stars are primary distance indicators, well known representatives of old and metal-poor stellar populations, and therefore regularly found in globular clusters. Here we report the discovery of a dozen RR Lyrae ab-type stars in the vicinity of the Galactic center, i.e. in the so-called nuclear stellar bulge of the Milky Way. This discovery provides the first direct observational evidence that the Galactic nuclear stellar bulge contains ancient stars (>10 Gyr old). Based on this we conclude that merging globular clusters likely contributed to building-up the high stellar density in the nuclear stellar bulge of the Milky Way.


title: Negative Circles in Signed Graphs: A Problem Collection

abstract: I propose that most problems about circles (cycles, circuits) in ordinary graphs that have odd or even length find their proper setting in the theory of signed graphs, where each edge has a sign, $+$ or $-$. Even-circle and odd-circle problems correspond to questions about positive and negative circles in signed graphs. (The sign of a circle is the product of its edge signs.) I outline questions about circles in signed graphs, that seem natural and potentially important.


title: Forbidden Induced Subgraphs

abstract: In descending generality I survey: five partial orderings of graphs, the induced-subgraph ordering, and examples like perfect, threshold, and mock threshold graphs. The emphasis is on how the induced subgraph ordering differs from other popular orderings and leads to different basic questions.


title: Bright Room-Temperature Single Photon Emission from Defects in Gallium Nitride

abstract: Single photon emitters play a central role in many photonic quantum technologies. A promising class of single photon emitters consists of atomic color centers in wide-bandgap crystals, such as diamond silicon carbide and hexagonal boron nitride. However, it is currently not possible to grow these materials as sub-micron thick films on low-refractive index substrates, which is necessary for mature photonic integrated circuit technologies. Hence, there is great interest in identifying quantum emitters in technologically mature semiconductors that are compatible with suitable heteroepitaxies. Here, we demonstrate robust single photon emitters based on defects in gallium nitride (GaN), the most established and well understood semiconductor that can emit light over the entire visible spectrum. We show that the emitters have excellent photophysical properties including a brightness in excess of 500x10^3 counts/s. We further show that the emitters can be found in a variety of GaN wafers, thus offering reliable and scalable platform for further technological development. We propose a theoretical model to explain the origin of these emitters based on cubic inclusions in hexagonal gallium nitride. Our results constitute a feasible path to scalable, integrated on-chip quantum technologies based on GaN.


title: Thick-target transmission method for excitation functions of interaction cross sections

abstract: We propose a method, called as thick-target transmission (T3) method, to obtain an excitation function of interaction cross sections. In an ordinal experiment to measure the excitation function of interaction cross sections by the transmission method, we need to change the beam energy for each cross section. In the T3 method, the excitation function is derived from the beam attenuations measured at the targets of different thicknesses without changing the beam energy. The advantage of the T3 method is the simplicity and availability for radioactive beams. To confirm the availability, we perform a simulation for the 12C+27Al system with the PHITS code instead of actual experiments. Our results have large uncertainties but well reproduce the tendency of the experimental data.


title: A model with Darwinian dynamics on a rugged landscape

abstract: We discuss the population dynamics with selection and random diffusion, keeping the total population constant, in a fitness landscape associated with Constraint Satisfaction, a paradigm for difficult optimization problems. We obtain a phase diagram in terms of the size of the population and the diffusion rate, with a glass phase inside which the dynamics keeps searching for better configurations, and outside which deleterious `mutations' spoil the performance. The phase diagram is analogous to that of dense active matter in terms of temperature and drive.


title: Non-negative matrix factorization-based subband decomposition for acoustic source localization

abstract: A novel non-negative matrix factorization (NMF) based subband decomposition in frequency spatial domain for acoustic source localization using a microphone array is introduced. The proposed method decomposes source and noise subband and emphasises source dominant frequency bins for more accurate source representation. By employing NMF, delay basis vectors and their subband information in frequency spatial domain for each frame is extracted. The proposed algorithm is evaluated in both simulated noise and real noise with a speech corpus database. Experimental results clearly indicate that the algorithm performs more accurately than other conventional algorithms under both reverberant and noisy acoustic environments.


title: The statistics of low frequency radio interference at the Murchison Radio-astronomy Observatory

abstract: We characterize the low frequency radio-frequency interference (RFI) environment at the Murchison Radio-astronomy Observatory (MRO), the location selected for the low-frequency component of the Square Kilometre Array. Data were collected from the BIGHORNS instrument, located at the MRO, which records a contiguous bandwidth between 70 and 300 MHz, between November 2014 to March 2015 inclusive. The data were processed to identify RFI, and we describe a series of statistics in both the time and frequency domain, including modeling of the RFI occupancy and signal power as a series of distribution functions, with the goal of aiding future scientists and operation staff in observation planning.


title: Dynamical Conductivity of Dirac Materials

abstract: For graphene (a Dirac material) it has been theoretically predicted and experimentally observed that DC resistivity is proportional to $ T^4$ when the temperature is much less than Bloch- Gr\"{u}neisen ($\Theta_{BG}$) temperature and T linear in opposite case ($T>>\Theta_{BG}$). Going beyond the DC case, we investigate the dynamical conductivity in graphene using the powerful method of memory function formalism. In the DC (zero frequency regime) limit, we obtained the above mention behavior which was previously obtained using the Bloch-Boltzmann kinetic equation. In the finite frequency regime, we obtained several new results: (1) the generalized Drude scattering rate, in the zero temperature limit, shows $\omega^4 $ behavior at low frequencies ($\omega << k_B \Theta_{BG}/ \hbar$) and saturates at higher frequencies. We also observed the Holstein Mechanism, however, with different power laws from that in the case of metals; (2) At higher frequencies, $\omega>>k_B \Theta_{BG}/ \hbar$, and higher temperatures $T>>\Theta_{BG}$, we observed that the generalized Drude scattering rate is linear in temperature. In addition, several other results are also obtained. With the experimental advancement of this field, these results should be experimentally tested.


title: A fast semi-discrete Kansa method to solve the two-dimensional spatiotemporal fractional diffusion

abstract: Anomalous diffusion is a common phenomenon observed in underground solute transport, soil water infiltration and sediment movement, etc. Time and space fractional derivative advection-dispersion equation (FADE) has been widely employed as the governing equation to characterize above mentioned anomalous diffusion related processes. However, a main problem in application of time and space FADE model to describe the real-world mass transport processes, is its low computation efficiency for long-time range and large irregular domain cases. This study offers a new algorithm in which the Kansa method is used for vector space fractional derivative term discretization and then analytical approach for resulted time fractional ordinary system. The influence of node distribution mode and node numbers on accuracy and convergence rate are analysed through the numerical examples in one and two dimensional cases. To test the application potentials of present methods, we offer the numerical results of two dimensional time and space FADE models in continuous and discrete cases. It shows that the solute plumes in heterogeneous and anisotropic media, can be well simulated by using present method compared with the previous time-consumed particle Monte-Carlo methods.


title: Equationally extremal semilattices

abstract: In the current paper we study extremal semilattices with respect to their equational properties. In the class $\mathbf{S}_n$ of all semilattices of order $n$ we find semilattices which have maximal (minimal) number of consistent equations. Moreover, we find a semilattice in $\mathbf{S}_n$ with maximal sum of numbers of solutions of equations.


title: Dynamics of piecewise translation maps

abstract: In this paper, we introduce a generalized piecewise translation map on the Euclidean space. We provide a special case when this map is always of finite type. For a finite type map in this case, we form conjectures on the semi-continuity of the attractor. Moreover, we also provide some conjectures and pictures on piecewise translations on a two dimensional torus.


title: Littlewood-Paley theorem, Nikolskii inequality, Besov spaces, Fourier and spectral multipliers on graded Lie groups

abstract: In this paper we investigate Besov spaces on graded Lie groups. We prove a Nikolskii type inequality (or the Reverse H\"older inequality) on graded Lie groups and as consequence we obtain embeddings of Besov spaces. We prove a version of the Littlewood-Paley theorem on graded Lie groups. The results are applied to obtain embedding properties of Besov spaces and multiplier theorems for both spectral and Fourier multipliers in Besov spaces on graded Lie groups. In particular, we give a number of sufficient conditions for the boundedness of Fourier multipliers in Besov spaces.


title: Optimal Distributed Stochastic Mirror Descent for Strongly Convex Optimization

abstract: In this paper we consider convergence rate problems for stochastic strongly-convex optimization in the non-Euclidean sense with a constraint set over a time-varying multi-agent network. We propose two efficient non-Euclidean stochastic subgradient descent algorithms based on the Bregman divergence as distance-measuring function rather than the Euclidean distances that were employed by the standard distributed stochastic projected subgradient algorithms. For distributed optimization of nonsmooth and strongly convex functions whose only stochastic subgradients are available, the first algorithm recovers the best previous known rate of O(ln(T)/T) (where T is the total number of iterations). The second algorithm is an epoch variant of the first algorithm that attains the optimal convergence rate of O(1/T), matching that of the best previously known centralized stochastic subgradient algorithm. Finally, we report some simulation results to illustrate the proposed algorithms.


title: Wire active microrheology to differentiate viscoelastic liquids from soft solids

abstract: Viscoelastic liquids are characterized by a finite static viscosity and a zero yield stress, whereas soft solids have an infinite viscosity and a non-zero yield stress. The rheological nature of viscoelastic materials has long been a challenge, and it is still a matter of debate. Here, we provide for the first time the constitutive equations of linear viscoelasticity for magnetic wires in yield stress materials, together with experimental measurements using Magnetic Rotational Spectroscopy (MRS). With MRS, the wires are submitted to a rotational magnetic field as a function of frequency and the wire motion is monitored by time-lapse microscopy. The soft solids studied are gel-forming polysaccharide aqueous dispersions (gellan gum) at concentrations above the gelification point. It is found that soft solids exhibit a clear and distinctive signature compared to viscous and viscoelastic liquids. In particular, the wire average rotation velocity equals zero over a broad frequency range. We also show the MRS technique is quantitative. From the wire oscillation amplitudes, the equilibrium elastic modulus is retrieved and agrees with polymer dynamics theory.


title: Wave-front error breakdown in LGS MOAO validated on-sky by CANARY

abstract: CANARY is the multi-object adaptive optics (MOAO) on-sky pathfinder developed in the perspective of Multi-Object Spectrograph on Extremely Large Telescopes~(ELTs). In 2013, CANARY was operated on-sky at the William Herschel telescope~(WHT), using three off-axis natural guide stars~(NGS) and four off-axis Rayleigh laser guide stars~(LGS), in open-loop, with the on-axis compensated turbulence observed with a H-band imaging camera and a Truth wave-front sensor~(TS) for diagnostic purposes. Our purpose is to establish a reliable and accurate wave-front error breakdown for LGS MOAO. This will enable a comprehensive analysis of \cana on-sky results and provide tools for validating simulations of MOAO systems for ELTs. To evaluate the MOAO performance, we compared the CANARY on-sky results running in MOAO, in Single Conjugated Adaptive Optics~(SCAO) and in Ground Layer Adaptive Optics~(GLAO) modes, over a large set of data acquired in 2013. We provide a statistical study of the seeing. We also evaluated the wave-front error breakdown from both analytic computations, one based on a MOAO system modelling and the other on the measurements from the CANARY TS. We have focussed especially on the tomographic error and we detail its vertical error decomposition~(VED). We show that CANARY obtained 30.1\%, 21.4\% and 17.1\% H-band Strehl ratios in SCAO, MOAO and GLAO respectively, for median seeing conditions with 0.66" of total seeing including 0.59" at the ground. Moreover, we get 99\% of correlation over 4,500 samples, for any AO modes, between two analytic computations of residual phase variance. Based on these variances, we obtain a reasonable Strehl-ratio~(SR) estimation when compared to the measured IR image SR. We evaluate the gain in compensation for the altitude turbulence brought by MOAO when compared to GLAO.


title: CMOS-Memristor Hybrid Integrated Pixel Sensors

abstract: Increase in image resolution require the ability of image sensors to pack an increased number of circuit components in a given area. On the the other hand a high speed processing of signals from the sensors require the ability of pixel to carry out pixel parallel operations. In the paper, we propose a modified 3T and 4T CMOS wide dynamic range pixels, which we refer as 2T-M and 3T-M configurations, comprising of MOSFETS and memristors. The low leakage currents and low area of memristors helps to achieve the objective of reducing the area, while the possibility to create arrays of memristors and MOSFETs across different layers within the chip, ensure the possibility to scale the circuit architecture.


title: Connected components of the moduli of elliptic K3 surfaces

abstract: The combinatorial type of an elliptic K3 surface with a zero section is the pair of the ADE -type of singular fibers and the torsion part of the Mordell-Weil group. We determine the set of connected components of the moduli of elliptic K3 surfaces with fixed combinatorial type. Our method relies on the theory of Miranda and Morrison on the structure of a genus of even indefinite lattices, and computer-aided calculations of p-adic quadratic forms.


title: Reasoning in the Bernays-Schoenfinkel-Ramsey Fragment of Separation Logic

abstract: Separation Logic (SL) is a well-known assertion language used in Hoare-style modular proof systems for programs with dynamically allocated data structures. In this paper we investigate the fragment of first-order SL restricted to the Bernays-Schoenfinkel-Ramsey quantifier prefix $\exists^*\forall^*$, where the quantified variables range over the set of memory locations. When this set is uninterpreted (has no associated theory) the fragment is PSPACE-complete, which matches the complexity of the quantifier-free fragment. However, SL becomes undecidable when the quantifier prefix belongs to $\exists^*\forall^*\exists^*$ instead, or when the memory locations are interpreted as integers with linear arithmetic constraints, thus setting a sharp boundary for decidability within SL. We have implemented a decision procedure for the decidable fragment of $\exists^*\forall^*$SL as a specialized solver inside a DPLL($T$) architecture, within the CVC4 SMT solver. The evaluation of our implementation was carried out using two sets of verification conditions, produced by (i) unfolding inductive predicates, and (ii) a weakest precondition-based verification condition generator. Experimental data shows that automated quantifier instantiation has little overhead, compared to manual model-based instantiation.


title: Adiabatic Quantum Computation

abstract: The quantum adiabatic theorem ensures that a slowly changing system, initially prepared in its ground state, will evolve to its final ground state with arbitrary precision. As a first result this thesis extends the original theorem to projection operators keeping the statement valid for Hamiltonians with degenerate ground spaces. Yet the main focus of this work lies in studying the efficiency of quantum circuit simulations by adabatic quantum computation. The standard Hamiltonian construction by Kitaev is based on a path graph reflecting the $L$ computation steps and influencing the scaling of the necessary evolution time by its spectral gap of $\mathcal{O}\left(\frac{1}{L^2}\right)$. Aspiring to an improved running time we generalize Kitaev's Hamiltonian to so-called standard graph Hamiltonians based on graph families with a different spectral gap. In this generalized construction the first two time derivatives of the Hamiltonian and the fraction of initial vertices appear as additional parameters of running time. In a first step the time derivatives can be proven to be constant. Expansion results from spectral graph theory however impose a trade-off between the spectral gap and the fraction of initial vertices as well as the fraction of final vertices which corresponds to the probability for obtaining the correct computational output. Graphs with spectral gap $\mathcal{O}\left(\frac{1}{L^k}\right)$, $k<2$, turn out to contradict very likely, graphs with $k<1$ even for sure at least one of the efficiency criterias for running time, output probability or Hamiltonian implementation. The above results may also be obstacles for a possible quantum PCP-theorem in complexity theory claiming the local Hamiltonian problem with constant gap to be QMA-complete since the very same Kitaev Hamiltonian is constructed in the QMA-hardness proof for the Local Hamiltonian problem.


title: A model for periodic blazars

abstract: We describe a scenario to explain blazar periodicities with timescales of $\sim$ few years. The scenario is based on a binary super-massive black hole (SMBH) system in which one of the two SMBH carries a jet. We discuss the various mechanisms that can cause the jet to precess and produce corkscrew patterns through space with a scale of $\sim$ few pc. It turns out that the dominant mechanism responsible for the precession is simply the imprint of the jet-carrying SMBH orbital speed on the jet. Gravitational deflection and Lense-Thirring precession (due to the gravitational field of the other SMBH) are second order effects. We complement the scenario with a kinematical jet model which is inspired to the spine-sheath structure observed in M87. One of the main advantages of such a structure is that it allows the peak of the synchrotron emission to scale with frequency according to $\nu F\propto \nu^{\xi}$ as the viewing angle is changed, where $\xi$ is not necessarily 3 or 4 as in the case of jets with uniform velocity, but can be $\xi \sim 1$. Finally, we apply the model to the source PG1553+113, which has been recently claimed to show a $T_{\rm obs}=(2.18\pm 0.08)\text{ yr}$ periodicity. We are able to reproduce the optical and gamma-ray light curves and multiple synchrotron spectra simultaneously. We also give estimates of the source mass and size.


title: Criteria of irreducibility of the Koopman representations for the group ${\rm GL}_0(2\infty,{\mathbb R})$

abstract: Our aim is to find the irreducibility criteria for the Koopman representation, when the group acts on some space with a measure (Conjecture 1.5). Some general necessary conditions of the irreducibility of this representation are established. In the particular case of the group ${\rm GL}_0(2\infty,{\mathbb R})$ $= \varinjlim_{n}{\rm GL}(2n-1,{\mathbb R})$, the inductive limit of the general linear groups we prove that these conditions are also the necessary ones. The corresponding measure is infinite tensor products of one-dimensional arbitrary Gaussian non-centered measures. The corresponding $G$-space $X_m$ is a subspace of the space ${\rm Mat}(2\infty,{\mathbb R})$ of infinite in both directions real matrices. In fact, $X_m$ is a collection of $m$ infinite in both directions rows. This result was announced in [20]. We give the proof only for $m\leq 2$. The general case will be studied later.


title: A $2\ell k$ Kernel for $\ell$-Component Order Connectivity

abstract: In the $\ell$-Component Order Connectivity problem ($\ell \in \mathbb{N}$), we are given a graph $G$ on $n$ vertices, $m$ edges and a non-negative integer $k$ and asks whether there exists a set of vertices $S\subseteq V(G)$ such that $|S|\leq k$ and the size of the largest connected component in $G-S$ is at most $\ell$. In this paper, we give a kernel for $\ell$-Component Order Connectivity with at most $2\ell k$ vertices that takes $n^{\mathcal{O}(\ell)}$ time for every constant $\ell$. On the way to obtaining our kernel, we prove a generalization of the $q$-Expansion Lemma to weighted graphs. This generalization may be of independent interest.


title: A Near-Linear Pseudopolynomial Time Algorithm for Subset Sum

abstract: Given a set $Z$ of $n$ positive integers and a target value $t$, the Subset Sum problem asks whether any subset of $Z$ sums to $t$. A textbook pseudopolynomial time algorithm by Bellman from 1957 solves Subset Sum in time $O(nt)$. This has been improved to $O(n \max Z)$ by Pisinger [J. Algorithms'99] and recently to $\tilde O(\sqrt{n} t)$ by Koiliaris and Xu [SODA'17]. Here we present a simple randomized algorithm running in time $\tilde O(n+t)$. This improves upon a classic algorithm and is likely to be near-optimal, since it matches conditional lower bounds from Set Cover and k-Clique. We then use our new algorithm and additional tricks to improve the best known polynomial space solution from time $\tilde O(n^3 t)$ and space $\tilde O(n^2)$ to time $\tilde O(nt)$ and space $\tilde O(n \log t)$, assuming the Extended Riemann Hypothesis. Unconditionally, we obtain time $\tilde O(n t^{1+\varepsilon})$ and space $\tilde O(n t^\varepsilon)$ for any constant $\varepsilon > 0$.


title: SP-rings with zero-divisors

abstract: We characterize the commutative rings whose ideals (resp. regular ideals) are products of radical ideals.


title: A New Perspective on Randomized Gossip Algorithms

abstract: In this short note we propose a new approach for the design and analysis of randomized gossip algorithms which can be used to solve the average consensus problem. We show how that Randomized Block Kaczmarz (RBK) method - a method for solving linear systems - works as gossip algorithm when applied to a special system encoding the underlying network. The famous pairwise gossip algorithm arises as a special case. Subsequently, we reveal a hidden duality of randomized gossip algorithms, with the dual iterative process maintaining a set of numbers attached to the edges as opposed to nodes of the network. We prove that RBK obtains a superlinear speedup in the size of the block, and demonstrate this effect through experiments.


title: Spin excitations in an all-organic double quantum dot molecule

abstract: We realize a strongly coupled double quantum dot in a single all-organic molecule by introducing a non-conjugated bridge in between two identical conjugated moieties. Spin-1/2 Kondo and Kondo enhanced low-energy excitations for respectively the odd and even electron occupation are observed in off-resonant transport. The ground state in the even occupation can be the singlet or the triplet state varying between samples. This observation suggests that both anti-ferromagnetic and ferromagnetic interactions between spins are of the same order of magnitude.


title: Bloch Oscillations in Optical and Zeeman Lattices in the Presence of Spin-Orbit Coupling

abstract: We address Bloch oscillations of a spin-orbit coupled atom in periodic potentials of two types: Optical and Zeeman lattices. We show that in optical lattices the spin-orbit coupling allows controlling the direction of atomic motion and may lead to complete suppression of the oscillations at specific values of the coupling strength. In Zeeman lattices the energy bands are found to cross each other at the boundaries of the Brillouin zone, resulting in period-doubling of the oscillations. In all cases, the oscillations are accompanied by rotation of the pseudo-spin, with a dynamics that is determined by the strength of the spin-orbit coupling. The predicted effects are discussed also in terms of a Wannier-Stark ladder, which in optical lattices consist of two mutually-shifted equidistant sub-ladders.


title: Study of Lambda polarization in relativistic nuclear collisions at sqrt{s}_NN = 7.7-200 GeV

abstract: We present a calculation of the global polarization of Lambda hyperons in relativistic Au-Au collisions at RHIC Beam Energy Scan range sqrt{s}_NN = 7.7 - 200 GeV with a 3+1 dimensional cascade + viscous hydro + cascade model, UrQMD+vHLLE. Within this model, the mean polarization of Lambda in the out-of-plane direction is predicted to decrease rapidly with collision energy from a top value of about 2% at the lowest energy examined. We explore the connection between the polarization signal and thermal vorticity and estimate the feed-down contribution to Lambda polarization due to the decay of higher mass hyperons.


title: Generalization of metric classification algorithms for sequences classification and labelling

abstract: The article deals with the issue of modification of metric classification algorithms. In particular, it studies the algorithm k-Nearest Neighbours for its application to sequential data. A method of generalization of metric classification algorithms is proposed. As a part of it, there has been developed an algorithm for solving the problem of classification and labelling of sequential data. The advantages of the developed algorithm of classification in comparison with the existing one are also discussed in the article. There is a comparison of the effectiveness of the proposed algorithm with the algorithm of CRF in the task of chunking in the open data set CoNLL2000.


title: Ideal Relaxation of the Hopf Fibration

abstract: Ideal MHD relaxation is the topology-conserving reconfiguration of a magnetic field into a lower energy state where the net force is zero. This is achieved by modeling the plasma as perfectly conducting viscous fluid. It is an important tool for investigating plasma equilibria and is often used to study the magnetic configurations in fusion devices and astrophysical plasmas. We study the equilibrium reached by a localized magnetic field through the topology conserving relaxation of a magnetic field based on the Hopf fibration in which magnetic field lines are closed circles that are all linked with one another. Magnetic fields with this topology have recently been shown to occur in non-ideal numerical simulations. Our results show that any localized field can only attain equilibrium if there is a finite external pressure, and that for such a field a Taylor state is unattainable. We find an equilibrium plasma configuration that is characterized by a lowered pressure in a toroidal region, with field lines lying on surfaces of constant pressure. Therefore, the field is in a Grad-Shafranov equilibrium. Localized helical magnetic fields are found when plasma is ejected from astrophysical bodies and subsequently relaxes against the background plasma, as well as on earth in plasmoids generated by e.g.\ a Marshall gun. This work shows under which conditions an equilibrium can be reached and identifies a toroidal depression as the characteristic feature of such a configuration.


title: Vibration-induced coherence enhancement of the performance of a biological quantum heat engine

abstract: Photosynthesis has been a long-standing research interest due to its fundamental importance. Recently, studies on photosynthesis processes also have inspired attention from a thermodynamical aspect when considering photosynthetic apparatuses as biological quantum heat engines. Quantum coherence is shown to play a crucial role in enhancing the performance of these quantum heat engines. Based on the experimentally reported structure, we propose a quantum heat engine model with a non-Markovian vibrational mode. We show that one can obtain a performance enhancement easily for a wide range of parameters in the presence of the vibrational mode. Our results provide insights into the photosynthetic processes and a design principle mimicking natural organisms.


title: Acoustically tuneable optical transmission through a subwavelength hole with a bubble

abstract: Efficient manipulation of light with sound in subwavelength-sized volumes is important for applications in photonics, phononics and biophysics, but remains elusive. We theoretically demonstrate the control of light with MHz-range ultrasound in a subwavelength, 300 nm wide water-filled hole with a 100 nm radius air bubble. Ultrasound-driven pulsations of the bubble modulate the effective refractive index of the hole aperture, which gives rise to spectral tuning of light transmission through the hole. This control mechanism opens up novel opportunities for tuneable acousto-optic and optomechanical metamaterials, and all-optical ultrasound transduction.

categories: physics.optics cond-mat.other physics.bio-ph physics.flu-dyn


title: Control of birhythmicity through conjugate self-feedback: Theory and experiment

abstract: Birhythmicity arises in several physical, biological and chemical systems. Although, many control schemes are proposed for various forms of multistability, only a few exist for controlling birhythmicity. In this paper we investigate the control of birhythmic oscillation by introducing a self-feedback mechanism that incorporates the variable to be controlled and its canonical conjugate. Using a detailed analytical treatment, bifurcation analysis and experimental demonstrations we establish that the proposed technique is capable of eliminating birhythmicity and generates monorhythmic oscillation. Further, the detailed parameter space study reveals that, apart from monorhythmicity, the system shows transition between birhythmicity and other dynamical forms of bistability. This study may have practical applications in controlling birhythmic behavior of several systems, in particular in biochemical and mechanical processes.


title: Assembly and speed control in ion exchange based modular phoretic micro-swimmers

abstract: We report an experimental study on ion-exchange based modular micro-swimmers in low-salt water. Cationic ion-exchange particles and passive cargo particles assemble into self-propelling complexes, showing self-propulsion at speeds of several microns per second over extended distances and times. We quantify the assembly and speed of the complexes for different combinations of ion exchange particles and cargo particles, substrate types, salt types and concentrations, and cell geometries. Irrespective of experimental boundary conditions, we observe a regular development of the assembly shape with increasing number of cargo. Moreover, the swimming speed increases stepwise upon increasing the number of cargo and then saturates at a maximum speed, indicating an active role of cargo in modular swimming. We propose a geometric model of self-assembly to describe the experimental observations in a qualitative way. Our study also provides some constraints for future theoretical modelling and simulation.


title: Concentration-compactness at the mountain pass level for nonlocal Schr\"{o}dinger equations

abstract: The aim of this paper is to study a concentration-compactness principle for inhomogeneous fractional Sobolev space $H^s (\mathbb{R}^N)$ for $0<s\leq N/2.$ As an application we establish Palais-Smale compactness for the Lagrangian associated to the fractional Schr\"{o}dinger equation $(-\Delta)^{s} u + a(x)u= f(x,u)$ for $0<s<1.$ Moreover, we prove the existence of nontrivial nonnegative solutions to this class of elliptic equations for a wide class of possible singular potentials $a(x)$; not necessarily bounded away from zero. We consider possible oscillatory nonlinearities and that may not satisfy the Ambrosetti-Rabinowitz condition and for both cases; subcritical and critical growth range which are superlinear at origin.


title: Incremental One-Class Models for Data Classification

abstract: In this paper we outline a PhD research plan. This research contributes to the field of one-class incremental learning and classification in case of non-stationary environments. The goal of this PhD is to define a new classification framework able to deal with very small learning dataset at the beginning of the process and with abilities to adjust itself according to the variability of the incoming data which create large scale datasets. As a preliminary work, incremental Covariance-guided One-Class Support Vector Machine is proposed to deal with sequentially obtained data. It is inspired from COSVM which put more emphasis on the low variance directions while keeping the basic formulation of incremental One-Class Support Vector Machine untouched. The incremental procedure is introduced by controlling the possible changes of support vectors after the addition of new data points, thanks to the Karush-Kuhn-Tucker conditions, that have to be maintained on all previously acquired data. Comparative experimental results with contemporary incremental and non-incremental one-class classifiers on numerous artificial and real data sets show that our method results in significantly better classification performance.


title: Microscopic folding model analysis of the radiative $(n,\gamma)$ reactions near the $Z=28$ shell-closure and the weak s-process

abstract: The radiative thermal neutron capture cross sections over the range of thermal energies from 1 keV to 1 MeV are studied in statistical Hauser-Feshbach formalism. The optical model potential is constructed by folding the density dependent M3Y nucleon-nucleon interaction with radial matter densities of target nuclei obtained from relativistic-mean-field (RMF) theory. The standard nuclear reaction code TALYS1.8 is used for calculation of cross sections. The nuclei studied in the present work reside near the $Z=28$ proton shell closure and are of astrophysical interests taking part in p-, s-, and r-process of nucleosynthesis. The Maxwellian-averaged cross-section (MACS) values for energies important for astrophysical applications are presented.


title: Gapless spin-liquid ground state in the $S = 1/2$ kagome antiferromagnet

abstract: Frustrated quantum magnetism has moved to the forefront of physics research, posing fundamental questions concerning quantum disordered states, entanglement, topology and the nature of the quantum wavefunction. The defining problem in the field is one of the simplest, the ground state of the nearest-neighbour $S = 1/2$ antiferromagnetic Heisenberg model on the kagome lattice, but has defied all theoretical and numerical methods employed to date. We apply the formalism of tensor-network states (TNS), specifically the method of projected entangled simplex states (PESS), whose combination of a correct accounting for multipartite entanglement and infinite system size provides qualitatively new insight. By studying the ground-state energy, the staggered magnetization we find at all finite tensor bond dimensions and the effects of a second-neighbour coupling, we demonstrate that the ground state is a gapless spin liquid. We discuss the comparison with other numerical studies and the physical interpretation of the gapless ground state.


title: Shadows and quantum invariants

abstract: This is a PhD thesis about low dimensional topology, in particular knot thory in 3-manifolds also different from the 3-sphere, topological applications of quantum invariants, and Turaev's shadows. There is an introduction and a survey for these topics. The thesis uses skein theory and focues on the connected sum of copies of S^1xS^2 and on the 3-tours as ambient manifolds. The skein space of the 3-torus is computed. The Tait conjecture and Eisermann's theorem have been extended to the connected sums of copies of S^1xS^2. The knots and links in S^1xS^2 whose crossing number is at most 3 are tabulated.


title: A Radiative Neutrino Model with $SU(2)_L$ Triplet Fields

abstract: We propose a loop induced neutrino mass model, in which we introduce several exotic fermions and bosons with $SU(2)_L$ multiplet, and discuss various phenomenologies such as lepton flavor violations, muon anomalous magnetic moment, nonstandard interacting neutrinoless double beta decay, relic density of dark matter, and the possibility of the spin independent direct detection searches, imposing the constraints of oblique parameters. And we show a benchmark point to satisfy all the constraints and discuss our predictions.


title: Inferring Person-to-person Proximity Using WiFi Signals

abstract: Today's societies are enveloped in an ever-growing telecommunication infrastructure. This infrastructure offers important opportunities for sensing and recording a multitude of human behaviors. Human mobility patterns are a prominent example of such a behavior which has been studied based on cell phone towers, Bluetooth beacons, and WiFi networks as proxies for location. However, while mobility is an important aspect of human behavior, understanding complex social systems requires studying not only the movement of individuals, but also their interactions. Sensing social interactions on a large scale is a technical challenge and many commonly used approaches---including RFID badges or Bluetooth scanning---offer only limited scalability. Here we show that it is possible, in a scalable and robust way, to accurately infer person-to-person physical proximity from the lists of WiFi access points measured by smartphones carried by the two individuals. Based on a longitudinal dataset of approximately 800 participants with ground-truth interactions collected over a year, we show that our model performs better than the current state-of-the-art. Our results demonstrate the value of WiFi signals in social sensing as well as potential threats to privacy that they imply.


title: Proton-proton correlations in distinguishing the two-proton emission mechanism of $^{23}$Al and $^{22}$Mg

abstract: The proton-proton momentum correlation functions ($C_{pp}(q)$) for kinematically complete decay channels of $^{23}$Al $\rightarrow$ p + p + $^{21}$Na and $^{22}$Mg $\rightarrow$ p + p + $^{20}$Ne have been measured at the RIKEN RI Beam Factory. From the very different correlation strength of $C_{pp}(q)$ for $^{23}$Al and $^{22}$Mg, the source size and emission time information were extracted from the $C_{pp}(q)$ data by assuming a Gaussian source profile in the correlation function calculation code (CRAB). The results indicated that the mechanism of two-proton emission from $^{23}$Al was mainly sequential emission, while that of $^{22}$Mg was mainly three-body simultaneous emission. By combining our earlier results of the two-proton relative momentum and the opening angle, it is pointed out that the mechanism of two-proton emission could be distinguished clearly.


title: Anomalous Cooling of Coronal Loops with Turbulent Suppression of Thermal Conduction

abstract: We investigate the impact of turbulent suppression of parallel heat conduction on the cooling of post-flare coronal loops. Depending on the value of the mean free path $\lambda_T$ associated with the turbulent scattering process, we identify four main cooling scenarios. The overall temperature evolution, from an initial temperature in excess of $10^7$~K, is modeled in each case, highlighting the evolution of the dominant cooling mechanism throughout the cooling process. Comparison with observed cooling times allows the value of $\lambda_T$ to be constrained, and interestingly this range corresponds to situations where collision-dominated conduction plays a very limited role, or even no role at all, in the cooling of post-flare coronal loops.


title: Initial partonic eccentricity fluctuations in a multi-phase transport model

abstract: Initial partonic eccentricities in Au+Au collisions at center-of-mass energy $\sqrt{s_{NN}}$ = 200 GeV are investigated using a multi-phase transport model with string melting scenario. The initial eccentricities in different order of harmonics are studied using participant and cumulant definitions. Eccentricity in terms of second-, fourth- and sixth order cumulants as a function of number of participant nucleons are compared systematically with the traditional participant definition. The ratio of the cumulant eccentricities $\varepsilon\left\{4\right\}/\varepsilon\left\{2\right\}$ and $\varepsilon\left\{6\right\}/\varepsilon\left\{4\right\}$ are studied in comparison with the ratio of the corresponding flow harmonics. The conversion coefficients ($v_n/\varepsilon_n$) are explored up to fourth order harmonic based on cumulant method. Furthermore, studies on transverse momentum ($p_T$) and pseudo-rapidity ($\eta$) dependencies of eccentricities and their fluctuations are presented. As in ideal hydrodynamics initial eccentricities are expected to be closely related to the final flow harmonics in relativistic heavy-ion collisions, studies of the fluctuating initial condition in the AMPT model will shed light on the tomography properties of the initial source geometry.


title: Geometric Model of Topological Insulators from the Maxwell Algebra

abstract: We propose a novel geometric model of three-dimensional topological insulators in presence of an external electromagnetic field. The gapped boundary of these systems supports relativistic quantum Hall states and is described by a Chern-Simons theory with a gauge connection that takes values in the Maxwell algebra. This represents a non-central extension of the Poincar\'e algebra and takes into account both the Lorentz and magnetic-translation symmetries of the surface states. In this way, we derive a relativistic version of the Wen-Zee term, and we show that the non-minimal coupling between the background geometry and the electromagnetic field in the model is in agreement with the main properties of the relativistic quantum Hall states in the flat space.


title: Critical behaviour in the optimal generation of multipartite entanglement

abstract: Two systems whose correlations cannot be classically accounted for display the simplest instance of quantum entanglement. Although this two-party association has caused a revolution in the foundations and uses of quantum mechanics, genuine entanglement among several parties unveils a whole new class of phenomena and applications. In this work we suggest a way to prepare Dicke states from a tunable source of bipartite entanglement to investigate foundational issues. The scheme has the following distinctive features: (i) it relies on controlled information loss and {\it unentangled} measurements; (ii) irrespective of the source entanglement, whenever a Dicke state is produced, it is ideal; (iii) the optimal entanglement of the bipartite source undergoes a second-order-like transition depending on the parameters of the Dicke state to be produced. These properties lead to asymptotic results on the entanglement between any qubit belonging to a Dicke state and the remaining qubits.


title: Thermal and transport properties in central heavy-ion reactions around a few hundred MeV/nucleon

abstract: Thermalization process of nuclear matter in central fireball region of heavy-ion collisions is investigated by employing an extension model of Boltzmann-Uehling-Uhlenbeck, namely the Van der Waals Boltzmann-Uehling-Uhlenbeck (VdWBUU) model. Temperature ($T$) is extracted by the quantum Fermion fluctuation approach and other thermodynamic quantities, such as density ($\rho$), entropy density ($s$), shear viscosity ($\eta$), isospin diffusivity ($D_{I}$) and heat conductivity ($\kappa$), are also deduced. The liquid-like and gas-like phase signs are discussed through the behavior of shear viscosity during heavy-ion collisions process with the VdWBUU model.


title: Coherent spin dynamics in gadolinium-doped CaWO4 crystal

abstract: We report the first observation of Rabi oscillations in the spin-7/2 ensemble of trivalent gadolinium ions hosted in CaWO$_4$ single crystal. A number of transitions within the lowest electronic multiplet $^8S_{7/2}$ of Gd$^{3+}$ ion are studied using a combination of continuous-wave and pulsed electron paramagnetic resonance spectroscopy. The corresponding Rabi damping curves and the spin coherence times are detected at varying strengths of the microwave field. These data are well reproduced by a theoretical model which accounts for the intrinsic inhomogeneity of the microwave field within the microwave resonator and the magnetic dipole interactions in the diluted spin ensemble. The results indicate that the studied 8-level ground manifold of Gd$^{3+}$ ion can represent an effective three qubit quantum system.


title: A $p$-Laplacian Neumann problem with a possibly supercritical nonlinearity

abstract: We look for nonconstant, positive, radially nondecreasing solutions of the quasilinear equation $-\Delta_p u+u^{p-1}=f(u)$ with $p>2$, in the unit ball $B$ of $\mathbb R^N$, subject to homogeneous Neumann boundary conditions. The assumptions on the nonlinearity $f$ are very mild and allow it to be possibly supercritical in the sense of Sobolev embeddings. The main tools used are the truncation method and a mountain pass-type argument. In the pure power case, i.e., $f(u)=u^{q-1}$, we detect the limit profile of the solutions of the problems as $q\to\infty$.


title: Thermodynamic bounds on equilibrium fluctuations of a global or local order parameter

abstract: We analyze thermodynamic bounds on equilibrium fluctuations of an order parameter, which are analogous to relations, which have been derived recently in the context of non-equilibrium fluctuations of currents. We discuss the case of {\it global} fluctuations when the order parameter is measured in the full system of interest, and {\it local} fluctuations, when the order parameter is evaluated only in a sub-part of the system. Using isometric fluctuation theorems, we derive thermodynamic bounds on the fluctuations of the order parameter in both cases. These bounds could be used to infer the value of symmetry breaking field or the relative size of the observed sub-system to the full system from {\it local} fluctuations.


title: A Scalar Curvature Formula For the Noncommutative 3-Torus

abstract: We compute the scalar curvature of a curved noncommutative 3-torus. To perturb the flat metric, the standard volume form on the noncommutative 3-torus is conformally perturbed and the corresponding perturbed Laplacian is analyzed. Using Connes' pseudodifferential calculus for the noncommutative 3-torus, we explicitly compute the first three terms of the small time heat kernel expansion for the perturbed Laplacian. The third term of the expansion gives a local formula for the scalar curvature. Finally, we show that in the classical limit when the deformation parameters vanish, our formula coincides with the formula for the commutative case.


title: Drawing graphs using a small number of obstacles

abstract: An obstacle representation of a graph $G$ is a set of points in the plane representing the vertices of $G$, together with a set of polygonal obstacles such that two vertices of $G$ are connected by an edge in $G$ if and only if the line segment between the corresponding points avoids all the obstacles. The obstacle number ${\rm obs}(G)$ of $G$ is the minimum number of obstacles in an obstacle representation of $G$. We provide the first non-trivial general upper bound on the obstacle number of graphs by showing that every $n$-vertex graph $G$ satisfies ${\rm obs}(G) \leq n\lceil\log{n}\rceil-n+1$. This refutes a conjecture of Mukkamala, Pach, and P\'alv\"olgyi. For $n$-vertex graphs with bounded chromatic number, we improve this bound to $O(n)$. Both bounds apply even when the obstacles are required to be convex. We also prove a lower bound $2^{\Omega(hn)}$ on the number of $n$-vertex graphs with obstacle number at most $h$ for $h<n$ and a lower bound $\Omega(n^{4/3}M^{2/3})$ for the complexity of a collection of $M \geq \Omega(n\log^{3/2}{n})$ faces in an arrangement of line segments with $n$ endpoints. The latter bound is tight up to a multiplicative constant.


title: Effective-one-body waveforms for binary neutron stars using surrogate models

abstract: Gravitational-wave observations of binary neutron star systems can provide information about the masses, spins, and structure of neutron stars. However, this requires accurate and computationally efficient waveform models that take <1s to evaluate for use in Bayesian parameter estimation codes that perform 10^7 - 10^8 waveform evaluations. We present a surrogate model of a nonspinning effective-one-body waveform model with l = 2, 3, and 4 tidal multipole moments that reproduces waveforms of binary neutron star numerical simulations up to merger. The surrogate is built from compact sets of effective-one-body waveform amplitude and phase data that each form a reduced basis. We find that 12 amplitude and 7 phase basis elements are sufficient to reconstruct any binary neutron star waveform with a starting frequency of 10Hz. The surrogate has maximum errors of 3.8% in amplitude (0.04% excluding the last 100M before merger) and 0.043 radians in phase. The version implemented in the LIGO Algorithm Library takes ~0.07s to evaluate for a starting frequency of 30Hz and ~0.8s for a starting frequency of 10Hz, resulting in a speed-up factor of ~10^3 - 10^4 relative to the original Matlab code. This allows parameter estimation codes to run in days to weeks rather than years, and we demonstrate this with a Nested Sampling run that recovers the masses and tidal parameters of a simulated binary neutron star system.


title: Star Streams and the Assembly History of the Galaxy

abstract: Thin halo star streams originate from the evaporation of globular clusters and therefore provide information about the early epoch globular cluster population. The observed tidal tails from halo globular clusters in the Milky Way are much more shorter than expected from a star cluster orbiting for 10~Gyr. The discrepancy is likely the result of the assumption that the clusters have been orbiting in a non-evolving galactic halo for a Hubble time. As a first step towards more realistic stream histories, a toy model that combines an idealized merger model with a simplified model of the internal collisional relaxation of individual star clusters is developed. On the average, the velocity dispersion increases with distance causing the density of the stream to decline with distance. Consequently, the streams visible in current data will normally be some fraction of the entire stream. Nevertheless, the high surface density segment of the stellar streams created from the evaporation of the more massive globular clusters should all be visible in low obscuration parts of the sky if closer than about 30~kpc. The Pan-STARRS1 halo volume is used to compare the numbers of halo streams and globular clusters.


title: Hook formulas for skew shapes II. Combinatorial proofs and enumerative applications

abstract: The Naruse hook-length formula is a recent general formula for the number of standard Young tableaux of skew shapes, given as a positive sum over excited diagrams of products of hook-lengths. In 2015 we gave two different $q$-analogues of Naruse's formula: for the skew Schur functions, and for counting reverse plane partitions of skew shapes. In this paper we give an elementary proof of Naruse's formula based on the case of border strips. For special border strips, we obtain curious new formulas for the Euler and $q$-Euler numbers in terms of certain Dyck path summations.


title: Heat, temperature and Clausius inequality in a model for active brownian particles

abstract: Methods of stochastic thermodynamics and hydrodynamics are applied to the a recently introduced model of active particles. The model consists of an overdamped particle subject to Gaussian coloured noise. Inspired by stochastic thermodynamics, we derive from the system's Fokker-Planck equation the average exchanges of heat and work with the active bath and the associated entropy production. We show that a Clausius inequality holds, with the local (non-uniform) temperature of the active bath replacing the uniform temperature usually encountered in equilibrium systems. Furthermore, by restricting the dynamical space to the first velocity moments of the local distribution function we derive a hydrodynamic description where local pressure, kinetic temperature and internal heat fluxes appear and are consistent with the previous thermodynamic analysis. The procedure also shows under which conditions one obtains the unified coloured noise approximation (UCNA): such an approximation neglects the fast relaxation to the active bath and therefore yields detailed balance and zero entropy production. In the last part, by using multiple time-scale analysis, we provide a constructive method (alternative to UCNA) to determine the solution of the Kramers equation and go beyond the detailed balance condition determining negative entropy production.


title: Measurements of the Inclusive Neutrino and Antineutrino Charged Current Cross Sections in MINERvA Using the Low-$\nu$ Flux Method

abstract: The total cross sections are important ingredients for the current and future neutrino oscillation experiments. We present measurements of the total charged-current neutrino and antineutrino cross sections on scintillator (CH) in the NuMI low-energy beamline using an {\em in situ} prediction of the shape of the flux as a function of neutrino energy from 2--50 GeV. This flux prediction takes advantage of the fact that neutrino and antineutrino interactions with low nuclear recoil energy ($\nu$) have a nearly constant cross section as a function of incident neutrino energy. This measurement is the lowest energy application of the low-$\nu$ flux technique, the first time it has been used in the NuMI antineutrino beam configuration, and demonstrates that the technique is applicable to future neutrino beams operating at multi-GeV energies. The cross section measurements presented are the most precise measurements to date below 5 GeV.


title: Comment on "Critical point scaling of Ising spin glasses in a magnetic field" by J. Yeo and M.A. Moore

abstract: In a section of a recent publication, [J. Yeo and M.A. Moore, Phys. Rev. B 91, 104432 (2015)], the authors discuss some of the arguments in the paper by Parisi and Temesv\'ari [Nuclear Physics B 858, 293 (2012)]. In this comment, it is shown how these arguments are misinterpreted, and the existence of the Almeida-Thouless transition in the upper critical dimension 6 reasserted.


title: On thermal production of self-interacting dark matter

abstract: We consider thermal production mechanisms of self-interacting dark matter in models with gauged $Z_3$ symmetry. A complex scalar dark matter is stabilized by the $Z_3$, that is the remnant of a local dark $U(1)_d$. Light dark matter with large self-interaction can be produced from thermal freeze-out in the presence of SM-annihilation, SIMP and/or forbidden channels. We show that dark photon and/or dark Higgs should be relatively light for unitarity and then assist the thermal freeze-out. We identify the constraints on the parameter space of dark matter self-interaction and mass in cases that one or some of the channels are important in determining the relic density.


title: On the frequency distribution of neutral particles from low-energy strong interactions

abstract: The rejection of the contamination, or background, from low-energy strong interactions at hadron collider experiments is a topic that has received significant attention in the field of particle physics. This article builds on a particle-level view of collision events, in line with recently-proposed subtraction methods. While conventional techniques in the field usually concentrate on probability distributions, our study is, to our knowledge, the first attempt at estimating the frequency distribution of background particles across the kinematic space inside individual collision events. In fact, while the probability distribution can generally be estimated given a model of low-energy strong interactions, the corresponding frequency distribution inside a single event typically deviates from the average and cannot be predicted a priori. We present preliminary results in this direction, and establish a connection between our technique and the particle weighting methods that have been the subject of recent investigation at the Large Hadron Collider.


title: On generalized trigonometric functions and series of rational functions

abstract: Here we introduce a way to construct generalized trigonometric functions associated with any complex polynomials, and the well known trigonometric functions can be seen to associate with polynomial $x^2-1$. We will show that those generalized trigonometric functions have algebraic identities which generalizes the well known $\sin^2(x)+\cos^2(x)=1$. One application of the generalized trigonometric functions is evaluating infinite series of rational functions.


title: Unsupervised clustering under the Union of Polyhedral Cones (UOPC) model

abstract: In this paper, we consider clustering data that is assumed to come from one of finitely many pointed convex polyhedral cones. This model is referred to as the Union of Polyhedral Cones (UOPC) model. Similar to the Union of Subspaces (UOS) model where each data from each subspace is generated from a (unknown) basis, in the UOPC model each data from each cone is assumed to be generated from a finite number of (unknown) \emph{extreme rays}.To cluster data under this model, we consider several algorithms - (a) Sparse Subspace Clustering by Non-negative constraints Lasso (NCL), (b) Least squares approximation (LSA), and (c) K-nearest neighbor (KNN) algorithm to arrive at affinity between data points. Spectral Clustering (SC) is then applied on the resulting affinity matrix to cluster data into different polyhedral cones. We show that on an average KNN outperforms both NCL and LSA and for this algorithm we provide the deterministic conditions for correct clustering. For an affinity measure between the cones it is shown that as long as the cones are not very coherent and as long as the density of data within each cone exceeds a threshold, KNN leads to accurate clustering. Finally, simulation results on real datasets (MNIST and YaleFace datasets) depict that the proposed algorithm works well on real data indicating the utility of the UOPC model and the proposed algorithm.


title: Preserving the value of large scale data analytics over time through selective re-computation

abstract: A pervasive problem in Data Science is that the knowledge generated by possibly expensive analytics processes is subject to decay over time, as the data used to compute it drifts, the algorithms used in the processes are improved, and the external knowledge embodied by reference datasets used in the computation evolves. Deciding when such knowledge outcomes should be refreshed, following a sequence of data change events, requires problem-specific functions to quantify their value and its decay over time, as well as models for estimating the cost of their re-computation. What makes this problem challenging is the ambition to develop a decision support system for informing data analytics re-computation decisions over time, that is both generic and customisable. With the help of a case study from genomics, in this vision paper we offer an initial formalisation of this problem, highlight research challenges, and outline a possible approach based on the collection and analysis of metadata from a history of past computations.


title: Black hole feeding and feedback: the physics inside the "subgrid"

abstract: Black holes (BHs) are believed to be a key ingredient of galaxy formation. However, the galaxy-BH interplay is challenging to study due to the large dynamical range and complex physics involved. As a consequence, hydrodynamical cosmological simulations normally adopt sub-grid models to track the unresolved physical processes, in particular BH accretion; usually the spatial scale where the BH dominates the hydrodynamical processes (the Bondi radius) is unresolved, and an approximate Bondi-Hoyle accretion rate is used to estimate the growth of the BH. By comparing hydrodynamical simulations at different resolutions (300, 30, 3 pc) using a Bondi-Hoyle approximation to sub-parsec runs with non-parameterized accretion, our aim is to probe how well an approximated Bondi accretion is able to capture the BH accretion physics and the subsequent feedback on the galaxy. We analyse an isolated galaxy simulation that includes cooling, star formation, Type Ia and Type II supernovae, BH accretion and AGN feedback (radiation pressure, Compton heating/cooling) where mass, momentum, and energy are deposited in the interstellar medium through conical winds. We find that on average the approximated Bondi formalism can lead to both over- and under-estimations of the BH growth, depending on resolution and on how the variables entering into the Bondi-Hoyle formalism are calculated.


title: Three-dimensional contact of TITH cartilage layers, a closed-form solution

abstract: Inhomogeneity and anisotropy play a crucial role in attributing articular cartilage its properties. The frictionless contact model constructed here consists in two thin biphasic transversely isotropic transversely homogeneous (TITH) cartilage layers firmly attached onto rigid substrates and shaped as elliptic paraboloids of different radii. Using asymptotic techniques, a solution to the deformation problem of such material has been recently obtained extending previous ones referred to homogeneous materials. The layer itself is thin in comparison with the size of the contact area and the observed time is shorter than the hydrogel characteristic time. The emerging three-dimensional contact problem is solved in closed-form and numerical benchmarks for constant and oscillating loads are given. The results are shown in terms of contact pressure and approach of the bones. The latter is derived to be directly proportional to the contact area. Existing experimental data are reinterpreted in view of the current model formulation. Comparisons are made with existing solutions for homogeneous biphasic materials in order to underline the functional importance of inhomogeneity in spreading the contact pressure distribution across the contact area. Particular attention is paid to the applicability of the retrieved formulas for interpreting measurements of in vivo experiments. Future directions are also prospected.


title: An extremality of the translation-invariant Gibbs measure for the HC-model on a Cayley tree

abstract: In this paper we study the extremality of translation-invariant Gibbs measure for the HC-model on a Cayley tree. It is known that for this model the translation-invariant measure is unique. We give a new proof of this statement and found regions of the extremality of this measure on the Cayley tree of order k.


title: Stability, instability, and blowup for time fractional and other non-local in time semilinear subdiffusion equations

abstract: We consider non-local in time semilinear subdiffusion equations on a bounded domain, where the kernel in the integro-differential operator belongs to a large class, which covers many relevant cases from physics applications, in particular the important case of fractional dynamics. The elliptic operator in the equation is given in divergence form with bounded measurable coefficients. We prove a well-posedness result in the setting of bounded weak solutions and study the stability and instability of the zero function in the special case where the nonlinearity vanishes at 0. We also establish a blowup result for positive convex and superlinear nonlinearities.


title: Magnetic liquids under high electric fields as optical diodes

abstract: We show and give examples of how unidirectional propagation of light rays in the limit of geometric optics could arise in some magnetic fluids due to the magnetoelectric effect under weak DC magnetic fields and strong DC electric fields around half of their dielectric breakdown. For such liquids as kerosene and transformer oils, one-way propagation of light may occur for 30 nm diameter magnetic nanoparticles (e.g. cobalt) and concentrations of 2% or larger.


title: Sensing Emotions in Text Messages: An Application and Deployment Study of EmotionPush

abstract: Instant messaging and push notifications play important roles in modern digital life. To enable robust sense-making and rich context awareness in computer mediated communications, we introduce EmotionPush, a system that automatically conveys the emotion of received text with a colored push notification on mobile devices. EmotionPush is powered by state-of-the-art emotion classifiers and is deployed for Facebook Messenger clients on Android. The study showed that the system is able to help users prioritize interactions.


title: Localized attack on clustering networks

abstract: Clustering network is one of which complex network attracting plenty of scholars to discuss and study the structures and cascading process. We primarily analyzed the effect of clustering coefficient to other various of the single clustering network under localized attack. These network models including double clustering network and star-like NON with clustering and random regular (RR) NON of ER networks with clustering are made up of at least two networks among which exist interdependent relation among whose degree of dependence is measured by coupling strength. We show both analytically and numerically, how the coupling strength and clustering coefficient effect the percolation threshold, size of giant component, critical coupling point where the behavior of phase transition changes from second order to first order with the increase of coupling strength between the networks. Last, we study the two types of clustering network: one type is same with double clustering network in which each subnetwork satisfies identical degree distribution and the other is that their subnetwork satisfies different degree distribution. The former type is treated both analytically and numerically while the latter is treated only numerically. In each section, we compared two results obtained from localized attack and random attack according to Shao et al:[22].


title: Uncertainty Estimates in the Heston Model via Fisher Information

abstract: We address the information content of European option prices about volatility in terms of the Fisher information matrix. We assume that observed option prices are centred on the theoretical price provided by Heston's model disturbed by additive Gaussian noise. We fit the likelihood function on the components of the VIX, i.e., near- and next-term put and call options on the S&P 500 with more than 23 days and less than 37 days to expiration and non-vanishing bid, and compute their Fisher information matrices from the Greeks in the Heston model. We find that option prices allow reliable estimates of volatility with negligible uncertainty as long as volatility is large enough. Interestingly, if volatility drops below a critical value, inferences from option prices become impossible because Vega, the derivative of a European option w.r.t. volatility, nearly vanishes.