title: Graph Partitioning using Quantum Annealing on the D-Wave System

abstract: In this work, we explore graph partitioning (GP) using quantum annealing on the D-Wave 2X machine. Motivated by a recently proposed graph-based electronic structure theory applied to quantum molecular dynamics (QMD) simulations, graph partitioning is used for reducing the calculation of the density matrix into smaller subsystems rendering the calculation more computationally efficient. Unconstrained graph partitioning as community clustering based on the modularity metric can be naturally mapped into the Hamiltonian of the quantum annealer. On the other hand, when constraints are imposed for partitioning into equal parts and minimizing the number of cut edges between parts, a quadratic unconstrained binary optimization (QUBO) reformulation is required. This reformulation may employ the graph complement to fit the problem in the Chimera graph of the quantum annealer. Partitioning into 2 parts, 2^N parts recursively, and k parts concurrently are demonstrated with benchmark graphs, random graphs, and small material system density matrix based graphs. Results for graph partitioning using quantum and hybrid classical-quantum approaches are shown to equal or out-perform current "state of the art" methods.


title: Logarithmic Hennings invariants for restricted quantum sl(2)

abstract: We construct a Hennings type logarithmic invariant for restricted quantum $\mathfrak{sl}(2)$ at a $2\mathsf{p}$-th root of unity. This quantum group $U$ is not braided, but factorizable. The invariant is defined for a pair: a 3-manifold $M$ and a colored link $L$ inside $M$. The link $L$ is split into two parts colored by central elements and by trace classes, or elements in the $0^{\text{th}}$ Hochschild homology of $U$, respectively. The two main ingredients of our construction are the universal invariant of a string link with values in tensor powers of $U$, and the modified trace introduced by the third author with his collaborators and computed on tensor powers of the regular representation. Our invariant is a colored extension of the logarithmic invariant constructed by Jun Murakami.


title: Calculations of Shear, Bulk viscosities and Electrical conductivity in Polyakov-Quark-Meson model

abstract: We have evaluated the transport coefficients of quark and hadronic matter in the frame work of Polyakov-Quark-Meson model. The thermal widths of quarks and mesons, which inversely control the strength of these transport coefficients, are obtained from the imaginary part of their respective self-energies at finite temperature. Due to the threshold conditions of their self energies, some limited temperature regions of quark and hadronic phase become relevant for our numerical predictions on transport coefficients, which are grossly in agreement with earlier results.


title: On Lipschitz rigidity of complex analytic sets

abstract: We prove that any complex analytic set in $\mathbb{C}^n$ which is Lipschitz normally embedded at infinity and has tangent cone at infinity that is a linear subspace of $\mathbb{C}^n$ must be an affine linear subspace of $\mathbb{C}^n$ itself. No restrictions on the singular set, dimension nor codimension are required. In particular, a complex algebraic set in $\mathbb{C}^n$ which is Lipschitz regular at infinity is an affine linear subspace.


title: Integral Field Spectroscopy of Balmer-Dominated Shocks in the Magellanic Cloud Supernova Remnant N103B

abstract: We present results of integral field spectroscopy of Balmer-dominated shocks in the LMC supernova remnant (SNR) N103B, carried out using the Wide Field Integral Spectrograph (WiFeS) on the 2.3 m telescope at Siding Spring Observatory in Australia. Existing X-ray studies of N103B have indicated a Type Ia supernova origin. Radiative shock emission from clumpy material surrounding the SNR may result from interaction of the forward shock with relic stellar wind material, possibly implicating a thermonuclear explosion in a single-degenerate binary system. The recently discovered Balmer-dominated shocks mark the impact of the forward shock with low density, partially neutral CSM gas, and form a partial shell encircling clumps of material exhibiting radiative shocks. The WiFeS spectra of N103B cover the range 3500 - 7000 Angstroms and reveal broad Halpha profiles up to 2360 km/s in width along the northern rim, and both Halpha and Hbeta broad profiles approximately 1300 km/s in width along the southern rim. Our fits to the Halpha line profiles indicate that in addition to from the usual broad and narrow emission components, a third component of intermediate width exists in the Balmer-dominated shocks, ranging from around 125 km/s up to 225 km/s in width. Recent Balmer-dominated shock models have predicted such an intermediate-width component will be generated in a fast neutral precursor. We derive a Sedov age of approximately 685 years for N103B from the Balmer-dominated spectra, consistent with the young age of 380 - 860 years estimated from light echo studies.


title: The effect of hardware-computed travel-time on localization accuracy in the inversion of experimental (acoustic) waveform data

abstract: This study aims to advance hardware-level computations for travel-time tomography applications in which the wavelength is close to the diameter of the information that has to be recovered. Such can be the case, for example, in the imaging applications of (1) biomedical physics, (2) astro-geophysics and (3) civil engineering. Our aim is to shed light on the effect of that preprocessing the digital waveform signal has on the inversion results and to find computational solutions that guarantee robust inversion when there are incomplete and/or noisy measurements. We describe a hardware-level implementation for integrated and thresholded travel-time computation (ITT and TTT). We compare the ITT and TTT approaches in inversion analysis with experimental acoustic travel-time data recorded using a ring geometry for the transmission and measurement points. The results obtained suggest that ITT is essential for maintaining the robustness of the inversion with imperfect signal digitization and sparsity. In order to ensure the relevance of the results, the specifications of the test setup were related to those of applications (1)-(3).


title: Trimming the Hill estimator: robustness, optimality and adaptivity

abstract: We introduce a trimmed version of the Hill estimator for the index of a heavy-tailed distribution, which is robust to perturbations in the extreme order statistics. In the ideal Pareto setting, the estimator is essentially finite-sample efficient among all unbiased estimators with a given strict upper break-down point. For general heavy-tailed models, we establish the asymptotic normality of the estimator under second order conditions and discuss its minimax optimal rate in the Hall class. We introduce the so-called trimmed Hill plot, which can be used to select the number of top order statistics to trim. We also develop an automatic, data-driven procedure for the choice of trimming. This results in a new type of robust estimator that can {\em adapt} to the unknown level of contamination in the extremes. As a by-product we also obtain a methodology for identifying extreme outliers in heavy tailed data. The competitive performance of the trimmed Hill and adaptive trimmed Hill estimators is illustrated with simulations.


title: Agatha: disentangling periodic signals from correlated noise in a periodogram framework

abstract: Periodograms are used as a key significance assessment and visualisation tool to display the significant periodicities in unevenly sampled time series. We introduce a framework of periodograms, called "Agatha", to disentangle periodic signals from correlated noise and to solve the 2-dimensional model selection problem: signal dimension and noise model dimension. These periodograms are calculated by applying likelihood maximization and marginalization and combined in a self-consistent way. We compare Agatha with other periodograms for the detection of Keplerian signals in synthetic radial velocity data produced for the Radial Velocity Challenge as well as in radial velocity datasets of several Sun-like stars. In our tests we find Agatha is able to recover signals to the adopted detection limit of the radial velocity challenge. Applied to real radial velocity, we use Agatha to confirm previous analysis of CoRoT-7 and to find two new planet candidates with minimum masses of 15.1 $M_\oplus$ and 7.08 $M_\oplus$ orbiting HD177565 and HD41248, with periods of 44.5 d and 13.4 d, respectively. We find that Agatha outperforms other periodograms in terms of removing correlated noise and assessing the significances of signals with more robust metrics. Moreover, it can be used to select the optimal noise model and to test the consistency of signals in time. Agatha is intended to be flexible enough to be applied to time series analyses in other astronomical and scientific disciplines. Agatha is available at http://www.agatha.herts.ac.uk.


title: Stable splitting of mapping spaces via nonabelian Poincar\'e duality

abstract: We use nonabelian Poincar\'e duality to recover the stable splitting of compactly supported mapping spaces, $\rm{Map_c}$$(M,\Sigma^nX)$, where $M$ is a parallelizable $n$-manifold. Our method for deriving this splitting is new, and naturally extends to give a more general stable splitting of the space of compactly supported sections of a certain bundle on $M$ with fibers $\Sigma^nX$, twisted by the tangent bundle of $M$. This generalization incorporates possible $O(n)$-actions on $X$ as well as accommodating non-parallelizable manifolds.


title: Monopole-antimonopole Interaction Potential

abstract: We numerically study the interactions of twisted monopole-antimonopole pairs in the 't Hooft-Polyakov model for a range of values of the scalar to vector mass ratio. We also recover the sphaleron solution at maximum twist discovered by Taubes, and map out its energy and size as functions of parameters.


title: Nonlinear diffusion & thermo-electric coupling in a two-variable model of cardiac action potential

abstract: This work reports the results of the theoretical investigation of nonlinear dynamics and spiral wave breakup in a generalized two-variable model of cardiac action potential accounting for thermo-electric coupling and diffusion nonlinearities. As customary in excitable media, the common Q10 and Moore factors are used to describe thermo-electric feedback in a 10-degrees range. Motivated by the porous nature of the cardiac tissue, in this study we also propose a nonlinear Fickian flux formulated by Taylor expanding the voltage dependent diffusion coefficient up to quadratic terms. A fine tuning of the diffusive parameters is performed a priori to match the conduction velocity of the equivalent cable model. The resulting combined effects are then studied by numerically simulating different stimulation protocols on a one-dimensional cable. Model features are compared in terms of action potential morphology, restitution curves, frequency spectra and spatio-temporal phase differences. Two-dimensional long-run simulations are finally performed to characterize spiral breakup during sustained fibrillation at different thermal states. Temperature and nonlinear diffusion effects are found to impact the repolarization phase of the action potential wave with non-monotone patterns and to increase the propensity of arrhythmogenesis.


title: Stress theory for classical fields

abstract: Classical field theories together with the Lagrangian and Eulerian approaches to continuum mechanics are embraced under a geometric setting of a fiber bundle. The base manifold can be either the body manifold of continuum mechanics, space manifold, or space-time. Differentiable sections of the fiber bundle represent configurations of the system and the configuration space containing them is given the structure of an infinite dimensional manifold. Elements of the cotangent bundle of the configuration space are interpreted as generalized forces and a representation theorem implies that there exist a stress object representing forces, non-uniquely. The properties of stresses are studies as well as the role of constitutive relations in the present general setting.


title: DeepMetabolism: A Deep Learning System to Predict Phenotype from Genome Sequencing

abstract: Life science is entering a new era of petabyte-level sequencing data. Converting such big data to biological insights represents a huge challenge for computational analysis. To this end, we developed DeepMetabolism, a biology-guided deep learning system to predict cell phenotypes from transcriptomics data. By integrating unsupervised pre-training with supervised training, DeepMetabolism is able to predict phenotypes with high accuracy (PCC>0.92), high speed (<30 min for >100 GB data using a single GPU), and high robustness (tolerate up to 75% noise). We envision DeepMetabolism to bridge the gap between genotype and phenotype and to serve as a springboard for applications in synthetic biology and precision medicine.


title: Requirements Model for Cyber-Physical System

abstract: The development of cyber-physical system (CPS) is a big challenge because of its complexity and its complex requirements. Especially in Requirements Engineering (RE), there exist many redundant and conflict requirements. Eliminating conflict requirements and merged redundant/common requirements lead a challenging task at the elicitation phase in the requirements engineering process for CPS. Collecting and optimizing requirements through appropriate process reduce both development time and cost as every functional requirement gets refined and optimized at very first stage (requirements elicitation phase) of the whole development process. Existing researches have focused on requirements those have already been collected. However, none of the researches have worked on how the requirements are collected and refined. This paper provides a requirements model for CPS that gives a direction about the requirements be gathered, refined and cluster in order to developing the CPS independently. The paper also shows a case study about the application of the proposed model to transport system.


title: Partial Domination in Graphs

abstract: A set $S\subseteq V$ is a dominating set of $G$ if every vertex in $V - S$ is adjacent to at least one vertex in $S$. The domination number $\gamma(G)$ of $G$ equals the minimum cardinality of a dominating set $S$ in $G$; we say that such a set $S$ is a $\gamma$-set. The single greatest focus of research in domination theory is the determination of the value of $\gamma(G)$. By definition, all vertices must be dominated by a $\gamma$-set. In this paper we propose relaxing this requirement, by seeking sets of vertices that dominate a prescribed fraction of the vertices of a graph. We focus particular attention on $1/2$ domination, that is, sets of vertices that dominate at least half of the vertices of a graph $G$. Keywords: partial domination, dominating set, partial domination number, domination number


title: On the pointwise iteration-complexity of a dynamic regularized ADMM with over-relaxation stepsize

abstract: In this paper, we extend the improved pointwise iteration-complexity result of a dynamic regularized alternating direction method of multipliers (ADMM) for a new stepsize domain. In this complexity analysis, the stepsize parameter can even be chosen in the interval $(0,2)$ instead of interval $(0,(1+\sqrt{5})/2)$. As usual, our analysis is established by interpreting this ADMM variant as an instance of a hybrid proximal extragradient framework applied to a specific monotone inclusion problem.


title: A simple yet effective baseline for 3d human pose estimation

abstract: Following the success of deep convolutional networks, state-of-the-art methods for 3d human pose estimation have focused on deep end-to-end systems that predict 3d joint locations given raw image pixels. Despite their excellent performance, it is often not easy to understand whether their remaining error stems from a limited 2d pose (visual) understanding, or from a failure to map 2d poses into 3-dimensional positions. With the goal of understanding these sources of error, we set out to build a system that given 2d joint locations predicts 3d positions. Much to our surprise, we have found that, with current technology, "lifting" ground truth 2d joint locations to 3d space is a task that can be solved with a remarkably low error rate: a relatively simple deep feed-forward network outperforms the best reported result by about 30\% on Human3.6M, the largest publicly available 3d pose estimation benchmark. Furthermore, training our system on the output of an off-the-shelf state-of-the-art 2d detector (\ie, using images as input) yields state of the art results -- this includes an array of systems that have been trained end-to-end specifically for this task. Our results indicate that a large portion of the error of modern deep 3d pose estimation systems stems from their visual analysis, and suggests directions to further advance the state of the art in 3d human pose estimation.


title: Lower Bound on the Localization Error in Infinite Networks with Random Sensor Locations

abstract: We present novel lower bounds on the mean square error (MSE) of the location estimation of an emitting source via a network where the sensors are deployed randomly. The sensor locations are modeled as a homogenous Poisson point process. In contrast to previous bounds which are a function of the specific locations of all the sensors, we present CRB-type bounds on the expected mean square error; that is, we first derive the CRB on the MSE as a function of the sensors' location, and then take expectation with respect to the distribution of the sensors' location. Thus, these bounds are not a function of a particular sensor configuration, but rather of the sensor statistics. Hence, these novel bounds can be evaluated prior to sensor deployment and provide insights into design issues such as the necessary sensor density, the effect of the channel model, the effect of the signal power, and others. The derived bounds are simple to evaluate and provide a good prediction of the actual network performance.


title: Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited

abstract: The problem of two van der Pol oscillators coupled by velocity delay terms was studied by Wirkus and Rand in 2002. The small-epsilon analysis resulted in a slow flow which contained delay terms. To simplify the analysis, Wirkus and Rand followed a common procedure of replacing the delay terms by non-delayed terms, a step said to be valid for small epsilon, resulting in a slow flow which was an ODE rather than a DDE (delay-differential equation). In the present paper we consider the same problem but leave the delay terms in the slow flow, thereby offering an evaluation of the approximate simplification made previously.


title: Approximate scattering state solutions of DKPE and SSE with Hellmann Potential

abstract: We study the approximate scattering state solutions of the Duffin-Kemmer-Petiau equation (DKPE) and the spinless Salpeter equation (SSE) with the Hellmann potential. The eigensolutions, scattering phase shifts, partial-waves transitions and the total cross section for all the partial waves are obtained and discussed. The dependence of partial-waves transitions on total angular momentum number, angular momentum number, mass combination and potential parameters were presented in the figures.


title: Scalable System Scheduling for HPC and Big Data

abstract: In the rapidly expanding field of parallel processing, job schedulers are the "operating systems" of modern big data architectures and supercomputing systems. Job schedulers allocate computing resources and control the execution of processes on those resources. Historically, job schedulers were the domain of supercomputers, and job schedulers were designed to run massive, long-running computations over days and weeks. More recently, big data workloads have created a need for a new class of computations consisting of many short computations taking seconds or minutes that process enormous quantities of data. For both supercomputers and big data systems, the efficiency of the job scheduler represents a fundamental limit on the efficiency of the system. Detailed measurement and modeling of the performance of schedulers are critical for maximizing the performance of a large-scale computing system. This paper presents a detailed feature analysis of 15 supercomputing and big data schedulers. For big data workloads, the scheduler latency is the most important performance characteristic of the scheduler. A theoretical model of the latency of these schedulers is developed and used to design experiments targeted at measuring scheduler latency. Detailed benchmarking of four of the most popular schedulers (Slurm, Son of Grid Engine, Mesos, and Hadoop YARN) are conducted. The theoretical model is compared with data and demonstrates that scheduler performance can be characterized by two key parameters: the marginal latency of the scheduler $t_s$ and a nonlinear exponent $\alpha_s$. For all four schedulers, the utilization of the computing system decreases to < 10\% for computations lasting only a few seconds. Multilevel schedulers that transparently aggregate short computations can improve utilization for these short computations to > 90\% for all four of the schedulers that were tested.


title: Dispersion-minimizing quadrature rules for $C^1$ quadratic isogeometric analysis

abstract: We develop quadrature rules for the isogeometric analysis of wave propagation and structural vibrations that minimize the discrete dispersion error of the approximation. The rules are optimal in the sense that they only require two quadrature points per element to minimize the dispersion error [1], and they are equivalent to the optimized blending rules we recently described. Our approach further simplifies the numerical integration: instead of blending two three-point standard quadrature rules, we construct directly a single two-point quadrature rule that reduces the dispersion error to the same order for uniform meshes with periodic boundary conditions. Also, we present a 2.5-point rule for both uniform and non-uniform meshes with arbitrary boundary conditions. Consequently, we reduce the computational cost by using the proposed quadrature rules. Various numerical examples demonstrate the performance of these quadrature rules.


title: Sharp phase transition for the random-cluster and Potts models via decision trees

abstract: We prove an inequality on decision trees on monotonic measures which generalizes the OSSS inequality on product spaces. As an application, we use this inequality to prove a number of new results on lattice spin models and their random-cluster representations. More precisely, we prove that 1. For the Potts model on transitive graphs, correlations decay exponentially fast for $\beta<\beta_c$. 2. For the random-cluster model with cluster weight $q\geq1$ on transitive graphs, correlations decay exponentially fast in the subcritical regime and the cluster-density satisfies the mean-field lower bound in the supercritical regime. 3. For the random-cluster models with cluster weight $q\geq1$ on planar quasi-transitive graphs $\mathbb{G}$, $$\frac{p_c(\mathbb{G})p_c(\mathbb{G}^*)}{(1-p_c(\mathbb{G}))(1-p_c(\mathbb{G}^*))}~=~q.$$ As a special case, we obtain the value of the critical point for the square, triangular and hexagonal lattices (this provides a short proof of the result of Beffara and Duminil-Copin [Probability Theory and Related Fields, 153(3-4):511--542, 2012]). These results have many applications for the understanding of the subcritical (respectively disordered) phase of all these models. The techniques developed in this paper have potential to be extended to a wide class of models including the Ashkin-Teller model, continuum percolation models such as Voronoi percolation and Boolean percolation, super-level sets of massive Gaussian Free Field, and random-cluster and Potts model with infinite range interactions.


title: A Nekhoroshev type theorem for the nonlinear Klein-Gordon equation with potential

abstract: We study the one-dimensional nonlinear Klein-Gordon (NLKG) equation with a convolution potential, and we prove that solutions with small analytic norm remain small for exponentially long times. The result is uniform with respect to $c \geq 1$, which however has to belong to a set of large measure.


title: The Weak Lensing Mass of Cosmic Web Filaments and Modified Gravity (MOG)

abstract: The weak lensing of cosmic web filaments is investigated in modified gravity (MOG) and it is demonstrated that the detected galaxies and baryonic gas of order $5-15 \%$ in filaments can with the enhanced value of the strength of gravity agree with the lensing data for filaments without dark matter.


title: Exploring many body localization and thermalization using semiclassical method

abstract: The Discrete Truncated Wigner Approximation (DTWA) is a semi-classical phase space method useful for the exploration of Many-body quantum dynamics. In this work we investigate Many-Body Localization (MBL) and thermalization using DTWA and compare its performance to exact numerical solutions. By taking as a benchmark case a 1D random field Heisenberg spin chain with short range interactions, and by comparing to numerically exact techniques, we show that DTWA is able to reproduce dynamical signatures that characterize both the thermal and the MBL phases. It exhibits the best quantitative agreement at short times deep in each of the phases and larger mismatches close to the phase transition. The DTWA captures the logarithmic growth of entanglement in the MBL phase, even though a pure classical mean-field analysis would lead to no dynamics at all. Our results suggest the DTWA can become a useful method to investigate MBL and thermalization in experimentally relevant settings intractable with exact numerical techniques, such as systems with long range interactions and/or systems in higher dimensions.


title: Wirtinger systems of generators of knot groups

abstract: We define the {\it Wirtinger number} of a link, an invariant closely related to the meridional rank. The Wirtinger number is the minimum number of generators of the fundamental group of the link complement over all meridional presentations in which every relation is an iterated Wirtinger relation arising in a diagram. We prove that the Wirtinger number of a link equals its bridge number. This equality can be viewed as establishing a weak version of Cappell and Shaneson's Meridional Rank Conjecture, and suggests a new approach to this conjecture. Our result also leads to a combinatorial technique for obtaining strong upper bounds on bridge numbers. This technique has so far allowed us to add the bridge numbers of approximately 50,000 prime knots of up to 14 crossings to the knot table. As another application, we use the Wirtinger number to show there exists a universal constant $C$ with the property that the hyperbolic volume of a prime alternating link $L$ is bounded below by $C$ times the bridge number of $L$.


title: Distributed Control for Spatial Self-Organization of Multi-Agent Swarms

abstract: In this work, we design distributed control laws for spatial self-organization of multi-agent swarms in 1D and 2D spatial domains. The objective is to achieve a desired density distribution over a simply-connected spatial domain. Since individual agents in a swarm are not themselves of interest and we are concerned only with the macroscopic objective, we view the network of agents in the swarm as a discrete approximation of a continuous medium and design control laws to shape the density distribution of the continuous medium. The key feature of this work is that the agents in the swarm do not have access to position information. Each individual agent is capable of measuring the current local density of agents and can communicate with its spatial neighbors. The network of agents implement a Laplacian-based distributed algorithm, which we call pseudo-localization, to localize themselves in a new coordinate frame, and a distributed control law to converge to the desired spatial density distribution. We start by studying self-organization in one-dimension, which is then followed by the two-dimensional case.


title: Completely Automated Equivalence Proofs

abstract: Verifying partial (i.e., termination-insensitive) equivalence of programs has significant practical applications in software development and education. Conventional equivalence verifiers typically rely on a combination of given relational summaries and suggested synchronization points; such information can be extremely difficult for programmers without a background in formal methods to provide for pairs of programs with dissimilar logic. In this work, we propose a completely automated verifier for determining partial equivalence, named Pequod. Pequod automatically synthesizes expressive proofs of equivalence conventionally only achievable via careful, manual constructions of product programs To do so, Pequod syntheses relational proofs for selected pairs of program paths and combines the per-path relational proofs to synthesize relational program invariants. To evaluate Pequod, we implemented it as a tool that targets Java Virtual Machine bytecode and applied it to verify the equivalence of hundreds of pairs of solutions submitted by students for problems hosted on popular online coding platforms, most of which could not be verified by existing techniques.


title: CAD Priors for Accurate and Flexible Instance Reconstruction

abstract: We present an efficient and automatic approach for accurate reconstruction of instances of big 3D objects from multiple, unorganized and unstructured point clouds, in presence of dynamic clutter and occlusions. In contrast to conventional scanning, where the background is assumed to be rather static, we aim at handling dynamic clutter where background drastically changes during the object scanning. Currently, it is tedious to solve this with available methods unless the object of interest is first segmented out from the rest of the scene. We address the problem by assuming the availability of a prior CAD model, roughly resembling the object to be reconstructed. This assumption almost always holds in applications such as industrial inspection or reverse engineering. With aid of this prior acting as a proxy, we propose a fully enhanced pipeline, capable of automatically detecting and segmenting the object of interest from scenes and creating a pose graph, online, with linear complexity. This allows initial scan alignment to the CAD model space, which is then refined without the CAD constraint to fully recover a high fidelity 3D reconstruction, accurate up to the sensor noise level. We also contribute a novel object detection method, local implicit shape models (LISM) and give a fast verification scheme. We evaluate our method on multiple datasets, demonstrating the ability to accurately reconstruct objects from small sizes up to $125m^3$.


title: Multi-objective integer programming: Synergistic parallel approaches

abstract: Exactly solving multi-objective integer programming (MOIP) problems is often a very time consuming process, especially for large and complex problems. Parallel computing has the potential to significantly reduce the time taken to solve such problems, but only if suitable algorithms are used. The first of our new algorithms follows a simple technique that demonstrates impressive performance for its design. We then go on to introduce new theory for developing more efficient parallel algorithms. The theory utilises elements of the symmetric group to apply a permutation to the objective functions to assign different workloads, and applies to algorithms that order the objective functions lexicographically. As a result, information and updated bounds can be shared in real time, creating a synergy between threads. We design and implement two algorithms that take advantage of such theory. To properly analyse the running time of our three algorithms, we compare them against two existing algorithms from the literature, and against using multiple threads within our chosen IP solver, CPLEX. This survey of six different parallel algorithms, the first of its kind, demonstrates the advantages of parallel computing. Across all problem types tested, our new algorithms are on par with existing algorithms on smaller cases and massively outperform the competition on larger cases. These new algorithms, and freely available implementations, allows the investigation of complex MOIP problems with four or more objectives.


title: K\"ulshammer ideals of graded categories and Hochschild cohomology

abstract: We generalize the notion of K\"ulshammer ideals to the setting of a graded category. This allows us to define and study some properties of K\"ulshammer type ideals in the graded center of a triangulated category and in the Hochschild cohomology of an algebra, providing new derived invariants. Further properties of K\"ulshammer ideals are studied in the case where the category is $d$-Calabi-Yau.


title: Turbulent stagnation in a z-pinch plasma

abstract: The ion kinetic energy in a stagnating plasma was previously determined by Kroupp et al. (Kroupp et al., 2011), from Doppler-dominated lineshapes augmented by measurements of plasma properties and assuming a uniform-plasma model. Notably, the energy was found to be dominantly stored in hydrodynamic flow. Here, we advance a new description of this stagnation as supersonically turbulent. Such turbulence implies a non-uniform density distribution. We demonstrate how to re-analyze the spectroscopic data consistent with the turbulent picture, and show this leads to better concordance of the overconstrained spectroscopic measurements, while also substantially lowering the inferred mean density.


title: Classification of Fixed Point Network Dynamics From Multiple Node Timeseries Data

abstract: Fixed point networks are dynamic networks that encode stimuli via distinct output patterns. Although such networks are omnipresent in neural systems, their structures are typically unknown or poorly characterized. It is therefore valuable to use a supervised approach for resolving how a network encodes distinct inputs of interest, and the superposition of those inputs from sampled multiple node time series. In this paper we show that accomplishing such a task involves finding a low-dimensional state space from supervised recordings. We demonstrate that standard methods for dimension reduction are unable to provide the desired functionality of optimal separation of the fixed points and transient trajectories to them. However, the combination of dimension reduction with selection and optimization can successfully provide such functionality. Specifically, we propose two methods: Exclusive Threshold Reduction (ETR) and Optimal Exclusive Threshold Reduction (OETR) for finding a basis for the classification state space. We show that the classification space constructed upon combination of dimension reduction optimal separation can directly facilitate recognition of stimuli, and classify complex inputs (mixtures) into similarity classes. We test our methodology and compare it to standard state-of-the-art methods on a benchmark dataset - an experimental neuronal network (the olfactory system) that we recorded from to test these methods. We show that our methods are capable of providing a basis for the classification space in such network, and to perform recognition at a significantly better rate than previously proposed approaches.


title: Spectral performance of Square Kilometre Array Antennas II: Calibration performance

abstract: We test the bandpass smoothness performance of two prototype Square Kilometre Array (SKA) SKA1-Low log-periodic dipole antennas, the SKALA2 and SKALA3 (`SKA Log-periodic Antenna'), and the current dipole from the Murchison Widefield Array (MWA) precursor telescope. Throughout this paper, we refer to the output complex-valued voltage response of an antenna when connected to a low noise amplifier (LNA), as the dipole bandpass. In Paper I (de Lera Acedo et al. 2017), the bandpass spectral response of the log-periodic antenna being developed for the SKA1-Low was estimated using numerical electromagnetic simulations and analyzed using low-order polynomial fittings and it was compared with the HERA antenna against the delay spectrum metric. In this work, realistic simulations of the SKA1-Low instrument, including frequency-dependent primary beams and array configuration, are used with a weighted least-squares polynomial estimator to assess the ability of prototype antennas to perform the SKA Epoch of Reionisation (EoR) statistical experiments. This work complements the ideal estimator tolerances computed for the proposed EoR science experiments in Trott & Wayth (2016), with the realised performance of an optimal and standard estimation (calibration) procedure. With a sufficient sky calibration model at higher frequencies, all antennas have bandpasses that are sufficiently smooth to meet the tolerances described in Trott & Wayth (2016) to perform the EoR statistical experiments, and these are primarily limited by an adequate sky calibration model, and the thermal noise level in the calibration data. At frequencies of the Cosmic Dawn (CD), which is of principal interest to SKA as one of the first next-generation telescopes capable of accessing higher redshifts, the MWA dipole and SKALA3 antenna have adequate performance, while the SKALA2 design will impede the ability to explore this era.


title: Polarized relations on horizontal SL(2)s

abstract: We introduce a relation on real conjugacy classes of SL(2)-orbits in a Mumford-Tate domain D which is compatible with natural partial orders on the sets of nilpotent orbits in the corresponding Lie algebra and boundary orbits in the compact dual. A generalization of the SL(2)-orbit theorem to such domains leads to an algorithm for computing this relation, which is worked out in several examples and special cases including period domains, Hermitian symmetric domains, and complete flag domains, and used to define a poset of equivalence classes of multivariable nilpotent orbits on D.


title: A 3D Ginibre point field

abstract: We introduce a three-dimensional random point field using the concept of the quaternion determinant. Orthogonal polynomials on the space of pure quaternions are defined, and used to construct a kernel function similar to the Ginibre kernel. We find explicit formulas for the polynomials and the kernel, and calculate their asymptotics in the bulk and at the center of coordinates.


title: The Future of Black Hole Astrophysics in the LIGO-VIRGO-LPF Era

abstract: There is a resurgence of interest in black holes sparked by the LIGO-VIRGO detection of stellar black hole mergers and recent astronomical investigations of jets and accretion disks which probe the spacetime geometry of black holes with masses ranging from a few times the mass of the sun to tens of billions of solar masses. Many of these black holes appear to be spinning rapidly. Some new approaches are described to studying how accreting black holes function as cosmic machines paying special attention to observations of AGN jets, especially with VLBI and $\gamma$-ray telescopes. It is assumed that these jets are powered by the electromagnetic extraction of the spin energy of their associated black holes, which are described by the Kerr metric, and that they become simpler and more electromagnetically dominated as the event horizon is approached. The major uncertainty in these models is in describing acceleration and transport of relativistic electrons and positrons and simple phenomenological prescriptions are proposed. The application of these ideas to M87 and 3C279 is outlined and the prospects for learning more, especially from the Event Horizon Telescope and the Cerenkov Telescope Array, are discussed. The main benefit of a better understanding of black hole astrophysics to the LISA mission should be a firmer understanding of the source demographics.


title: Multisolitons for the defocusing energy critical wave equation with potentials

abstract: We construct multisoliton solutions to the defocusing energy critical wave equation with potentials in $\mathbb{R}^{3}$ based on regular and reversed Strichartz estimates developed in \cite{GC3} for wave equations with charge transfer Hamiltonians. We also show the asymptotic stability of multisoliton solutions. The multisoliton structures with both stable and unstable solitons are covered. Since each soliton decays slowly with rate $\frac{1}{\left\langle x\right\rangle }$, the interactions among the solitons are strong. Some reversed Strichartz estimates and local decay estimates for the charge transfer model are established to handle strong interactions.


title: Angular dependence of magnetization reversal in epitaxial chromium telluride thin films with perpendicular magnetic anisotropy

abstract: We investigate magnetic anisotropy and magnetization reversal mechanism in chromium telluride thin films grown by molecular beam epitaxy. We report existence of strong perpendicular anisotropy in these thin films, along with a relatively strong second order anisotropy contribution. The angular variation of the switching field observed from the magnetoresistance measurement is explained quantitatively using a one-dimensional defect model. The model reveals the relative roles of nucleation and pinning in the magnetization reversal, depending on the applied field orientation. Micromagnetic simulations are performed to visualize the domain structure and switching process.


title: Convolutional Sequence to Sequence Learning

abstract: The prevalent approach to sequence to sequence learning maps an input sequence to a variable length output sequence via recurrent neural networks. We introduce an architecture based entirely on convolutional neural networks. Compared to recurrent models, computations over all elements can be fully parallelized during training and optimization is easier since the number of non-linearities is fixed and independent of the input length. Our use of gated linear units eases gradient propagation and we equip each decoder layer with a separate attention module. We outperform the accuracy of the deep LSTM setup of Wu et al. (2016) on both WMT'14 English-German and WMT'14 English-French translation at an order of magnitude faster speed, both on GPU and CPU.


title: A revision of the subtract-with-borrow random number generators

abstract: The most popular and widely used subtract-with-borrow generator, also known as RANLUX, is reimplemented as a linear congruential generator using large integer arithmetic with the modulus size of 576 bits. Modern computers, as well as the specific structure of the modulus inferred from RANLUX, allow for the development of a fast modular multiplication -- the core of the procedure. This was previously believed to be slow and have too high cost in terms of computing resources. Our tests show a significant gain in generation speed which is comparable with other fast, high quality random number generators. An additional feature is the fast skipping of generator states leading to a seeding scheme which guarantees the uniqueness of random number sequences.


title: A Mathematical Theory of Human Machine Teaming

abstract: We begin with a disquieting paradox: human machine teaming (HMT) often produces results worse than either the human or machine would produce alone. Critically, this failure is not a result of inferior human modeling or a suboptimal autonomy: even with perfect knowledge of human intention and perfect autonomy performance, prevailing teaming architectures still fail under trivial stressors~\cite{trautman-smc-2015}. This failure is instead a result of deficiencies at the \emph{decision fusion level}. Accordingly, \emph{efforts aimed solely at improving human prediction or improving autonomous performance will not produce acceptable HMTs: we can no longer model humans, machines and adversaries as distinct entities.} We thus argue for a strong but essential condition: HMTs should perform no worse than either member of the team alone, and this performance bound should be independent of environment complexity, human-machine interfacing, accuracy of the human model, or reliability of autonomy or human decision making.


title: Affinity Scheduling and the Applications on Data Center Scheduling with Data Locality

abstract: MapReduce framework is the de facto standard in Hadoop. Considering the data locality in data centers, the load balancing problem of map tasks is a special case of affinity scheduling problem. There is a huge body of work on affinity scheduling, proposing heuristic algorithms which try to increase data locality in data centers like Delay Scheduling and Quincy. However, not enough attention has been put on theoretical guarantees on throughput and delay optimality of such algorithms. In this work, we present and compare different algorithms and discuss their shortcoming and strengths. To the best of our knowledge, most data centers are using static load balancing algorithms which are not efficient in any ways and results in wasting the resources and causing unnecessary delays for users.


title: Analysis of Approximate Message Passing with a Class of Non-Separable Denoisers

abstract: Approximate message passing (AMP) is a class of efficient algorithms for solving high-dimensional linear regression tasks where one wishes to recover an unknown signal \beta_0 from noisy, linear measurements y = A \beta_0 + w. When applying a separable denoiser at each iteration, the performance of AMP (for example, the mean squared error of its estimates) can be accurately tracked by a simple, scalar iteration referred to as state evolution. Although separable denoisers are sufficient if the unknown signal has independent and identically distributed entries, in many real-world applications, like image or audio signal reconstruction, the unknown signal contains dependencies between entries. In these cases, a coordinate-wise independence structure is not a good approximation to the true prior of the unknown signal. In this paper we assume the unknown signal has dependent entries, and using a class of non-separable sliding-window denoisers, we prove that a new form of state evolution still accurately predicts AMP performance. This is an early step in understanding the role of non-separable denoisers within AMP, and will lead to a characterization of more general denoisers in problems including compressive image reconstruction.


title: Word and Phrase Translation with word2vec

abstract: Word and phrase tables are key inputs to machine translations, but costly to produce. New unsupervised learning methods represent words and phrases in a high-dimensional vector space, and these monolingual embeddings have been shown to encode syntactic and semantic relationships between language elements. The information captured by these embeddings can be exploited for bilingual translation by learning a transformation matrix that allows matching relative positions across two monolingual vector spaces. This method aims to identify high-quality candidates for word and phrase translation more cost-effectively from unlabeled data. This paper expands the scope of previous attempts of bilingual translation to four languages (English, German, Spanish, and French). It shows how to process the source data, train a neural network to learn the high-dimensional embeddings for individual languages and expands the framework for testing their quality beyond the English language. Furthermore, it shows how to learn bilingual transformation matrices and obtain candidates for word and phrase translation, and assess their quality.


title: Protein and hydration-water dynamics are decoupled: A new model connecting dynamics and biochemical function is required

abstract: Water plays a major role in bio-systems, greatly contributing to determine their structure, stability and even function. It is well know, for instance, that proteins require a minimum amount of water to be functionally active. Since the biological functions of proteins involve changes of conformation, and sometimes chemical reactions, it is natural to expect a connection of these functions with dynamical properties of the coupled system of proteins and their hydration water. However, despite many years of intensive research, the detailed nature of protein - hydration water interactions, and their effect on the biochemical activity of proteins through peculiar dynamical effects, is still partly unknown. In particular, models proposed so far, fail to explain the full set of experimental data. The well-accepted 'protein dynamical transition' scenario is based on perfect coupling between the dynamics of proteins and of their hydration water, which has never been confuted experimentally. We present high-energy resolution elastic neutron scattering measurements of the atomistic dynamics of the model protein, lysozyme, in water that were carried out on the IN16B spectrometer at the Institut Laue-Langevin in Grenoble, France. These show for the first time that the dynamics of proteins and of their hydration water are actually de-coupled. This important result militates against the well-accepted scenario, and requires a new model to link protein dynamics to the dynamics of its hydration water, and, in turn, to biochemical function.


title: Integrable Discrete Model for One-dimensional Soil Water Infiltration

abstract: We propose an integrable discrete model of one-dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection-diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time-dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self-adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the na\"ive discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.


title: Flexible Clustering for High-Dimensional Data via Mixtures of Joint Generalized Hyperbolic Models

abstract: A mixture of joint generalized hyperbolic distributions (MJGHD) is introduced for asymmetric clustering for high-dimensional data. The MJGHD approach takes into account the cluster-specific subspace, thereby limiting the number of parameters to estimate while also facilitating visualization of results. Identifiability is discussed, and a multi-cycle ECM algorithm is outlined for parameter estimation. The MJGHD approach is illustrated on two real data sets, where the Bayesian information criterion is used for model selection.


title: Extended Weyl fractional integrals and their inequalities

abstract: In the paper, the authors establish some interesting identities and inequalities involving the extended Weyl type fractional integrals.


title: Mechanisms of near-surface structural evolution in nanocrystalline materials during sliding contact

abstract: The wear-driven structural evolution of nanocrystalline Cu was simulated with molecular dynamics under constant normal loads, followed by a quantitative analysis. While the microstructure far away from the sliding contact remains unchanged, grain growth accompanied by partial dislocations and twin formation was observed near the contact surface, with more rapid coarsening promoted by higher applied normal loads. The structural evolution continues with increasing number of sliding cycles and eventually saturates to a stable distinct layer of coarsened grains, separated from the finer matrix by a steep gradient in grain size. The coarsening process is balanced by the rate of material removal when the normal load is high enough. The observed structural evolution leads to an increase in hardness and decrease in friction coefficient, which also saturate after a number of sliding cycles. This work provides important mechanistic understanding of nanocrystalline wear, while also introducing a methodology for atomistic simulations of cyclic wear damage under constant applied normal loads.


title: Finslerian dipolar modulation of the CMB power spectra at scales $2<l<600$

abstract: A common explanation for the CMB power asymmetry is to introduce a dipolar modulation at the stage of inflation, where the primordial power spectrum is spatially varying. If the universe in the stage of inflation is Finslerian, and if the Finsler spacetime is non-reversible under parity flip, $x\rightarrow-x$, then a three dimensional spectrum which is the function of wave vector and direction is valid. In this paper, a three dimensional primordial power spectrum with preferred direction is derived in the framework of Finsler spacetime. It is found that the amplitude of dipolar modulation is related to the Finslerian parameter, which in turn is a function of wave vector. The angular correlation coefficients are presented, and the numerical results for the anisotropic correlation coefficients over the multipole range $2<l<600$ are given.


title: Clustering Airbnb Reviews

abstract: In the last decade, online customer reviews increasingly exert influence on consumers' decision when booking accommodation online. The renewal importance to the concept of word-of mouth is reflected in the growing interests in investigating consumers' experience by analyzing their online reviews through the process of text mining and sentiment analysis. A clustering approach is developed for Boston Airbnb reviews submitted in the English language and collected from 2009 to 2016. This approach is based on a mixture of latent variable models, which provides an appealing framework for handling clustered binary data. We address here the problem of discovering meaningful segments of consumers that are coherent from both the underlying topics and the sentiment behind the reviews. A penalized mixture of latent traits approach is developed to reduce the number of parameters and identify variables that are not informative for clustering. The introduction of component-specific rate parameters avoids the over-penalization that can occur when inferring a shared rate parameter on clustered data. We divided the guests into four groups -- property driven guests, host driven guests, guests with recent overall negative stay and guests with some negative experiences.


title: A Model to Study Finite-Size and Magnetic Effects on the Phase Transition of a Fermion Interacting System

abstract: We present a model to study effects from an external magnetic field, chemical potential, and finite size, on the phase structure of a massive four- and six-fermion interacting system. These effects are introduced by a method of compactification of coordinates, a generalization of the standard Matsubara prescription. Through the compactification of the $z$ coordinate and of imaginary time, we describe a heated system with the shape of a film of thickness $L$, at temperature $\beta^{-1}$ undergoing first- or second-order phase transition. We have found a strong dependence of the temperature transition on the constants couplings $\lambda$ and $\eta$. Besides magnetic catalysis and symmetry breaking for both kinds of transition, we have found an inverse symmetry breaking phenomenon with respect to first-order phase transition.


title: Structure-preserving operators for thermal-nonequilibrium hydrodynamics

abstract: Radiation hydrodynamics simulations based on the one-fluid two-temperature model may violate the law of energy conservation because the governing equations are expressed in a nonconservative formulation. Here, we maintain the important physical requirements by employing a strategy based on the key concept that the mathematical structures associated with the conservative and nonconservative equations are preserved, even at the discrete level. To this end, we discretize the conservation laws and transform them via exact algebraic operations. The proposed scheme maintains the global conservation errors within the round-off level. In addition, a numerical experiment concerning the shock tube problem suggests that the proposed scheme well agrees with the jump conditions at the discontinuities regulated by the Rankine-Hugoniot relationship. The generalized derivation allows us to employ arbitrary central difference, artificial dissipation, and Runge-Kutta methods.


title: Discrete Games in Endogenous Networks: Equilibria and Policy

abstract: In games of friendship links and behaviors, I propose $k$-player Nash stability---a family of equilibria, indexed by a measure of robustness given by the number of permitted link changes, which is (ordinally and cardinally) ranked in a probabilistic sense. Application of the proposed framework to adolescents' tobacco smoking and friendship decisions suggests that: (a.) friendship networks respond to increases of tobacco prices and this response amplifies the intended policy effect on smoking, (b.) racially desegregating high-schools, via stimulating the social interactions of students with different intrinsic propensity to smoke, decreases the overall smoking prevalence, (c.) adolescents are averse to sharing friends so that there is a rivalry for friendships, (d.) when data on individuals' friendship network is not available, the importance of price centered policy tools is underestimated.


title: Characterization for entropy of shifts of finite type on Cayley trees

abstract: The notion of tree-shifts constitutes an intermediate class in between one-sided shift spaces and multidimensional ones. This paper proposes an algorithm for computing of the entropy of a tree-shift of finite type. Meanwhile, the entropy of a tree-shift of finite type is $\dfrac{1}{p} \ln \lambda$ for some $p \in \mathbb{N}$, where $\lambda$ is a Perron number. This extends Lind's work on one-dimensional shifts of finite type. As an application, the entropy minimality problem is investigated, and we obtain the necessary and sufficient condition for a tree-shift of finite type being entropy minimal with some additional conditions.


title: Large magnetic anisotropy in $Fe_xTaS_2$ single crystals

abstract: In intercalated transition metal dichalcogenide $Fe_xTaS_2$ (0.2 $\leq$ x $\leq$ 0.4) single crystals, large magnetic anisotropy is observed. Transport property measurements indicate that heavy Fe-doping leads to a large anisotropy of resistivity ($\rho$$_{c}$/$\rho$$_{ab}$). A sharp M-H hysteresis curve is observed with magnetic field along c-axis, while a linear magnetization appears with magnetic field applied in the ab-plane. The angular dependent magnetic susceptibility from in-plane to out-of-plane indicates that magnetic moments are strongly pinned along the c-axis in an unconventional manner and the coercive field reaches as large as 6 T at T = 5 K. First-principles calculation clearly suggests that the strong spin-orbital coupling give rise to such a large anisotropy of magnetism. The strong pinning effect of magnetic moments along c-axis makes this material a very promising candidate for the development of spin-aligner in spintronics devices.


title: Quarkonium $h$ States As Arbiters of Exoticity

abstract: The mass splitting between the quarkonium spin-singlet state $h$ ($J^{PC} = 1^{+-}$) and the spin average of the quarkonium spin-triplet states $\chi$ ($J^{PC} = 0^{++}, 1^{++}, 2^{++}$) is seen to be astonishingly small, not only in the charmonium and bottomonium cases where the relevant masses have been measured, but in positronium as well. We find, both in nonrelativistic quark models and in NRQCD, that this hyperfine splitting is so small that it can be used as a test of the pure $Q\bar Q$ content of the states. We discuss the $2P$ states of charmonium in the vicinity of $3.9$ GeV, where the putative exotics $X(3872)$ and $X(3915)$ have been seen and a new $\chi_{c0}(2P)$ candidate has been observed at Belle.


title: Lattice Boltzmann modeling and simulation of liquid jet breakup

abstract: A three-dimensional color-fluid lattice Boltzmann model for immiscible two-phase flows is developed in the framework of a three-dimensional 27-velocity (D3Q27) lattice. The collision operator comprises the D3Q27 versions of three sub-operators: a multiple-relaxation-time (MRT) collision operator, a generalized Liu--Valocchi--Kang perturbation operator, and a Latva-Kokko--Rothman recoloring operator. A D3Q27 version of an enhanced equilibrium distribution function is also incorporated into this model to improve the Galilean invariance. Three types of numerical tests, namely, a static droplet, an oscillating droplet, and the Rayleigh--Taylor instability, show a good agreement with analytical solutions and numerical simulations. Following these numerical tests, this model is applied to liquid-jet-breakup simulations. The simulation conditions are matched to the conditions of the previous experiments. In this case, numerical stability is maintained throughout the simulation, although the kinematic viscosity for the continuous phase is set as low as $1.8\times10^{-4}$, in which case the corresponding Reynolds number is $3.4\times10^{3}$; the developed lattice Boltzmann model based on the D3Q27 lattice enables us to perform the simulation with parameters directly matched to the experiments. The jet's liquid column transitions from an asymmetrical to an axisymmetrical shape, and entrainment occurs from the side of the jet. The measured time history of the jet's leading-edge position shows a good agreement with the experiments. Finally, the reproducibility of the regime map for liquid-liquid systems is assessed. The present lattice Boltzmann simulations well reproduce the characteristics of predicted regimes, including varicose breakup, sinuous breakup, and atomization.


title: Scars in Dirac fermion systems: the influence of an Aharonov--Bohm flux

abstract: Time-reversal ($\mathcal{T}$-) symmetry is fundamental to many physical processes. Typically, $\mathcal{T}$-breaking for microscopic processes requires the presence of magnetic field. However, for 2D massless Dirac billiards, $\mathcal{T}$-symmetry is broken automatically by the mass confinement, leading to chiral quantum scars. In this paper, we investigate the mechanism of $\mathcal{T}$-breaking by analyzing the local current of the scarring eigenstates and their magnetic response to an Aharonov--Bohm flux. Our results unveil the complete understanding of the subtle $\mathcal{T}$-breaking phenomena from both the semiclassical formula of chiral scars and the microscopic current and spin reflection at the boundaries, leading to a controlling scheme to change the chirality of the relativistic quantum scars. Our findings not only have significant implications on the transport behavior and spin textures of the relativistic pseudoparticles, but also add basic knowledge to relativistic quantum chaos.


title: Native surface oxide turns alloyed silicon membranes into nanophononic metamaterials with ultra-low thermal conductivity

abstract: A detailed understanding of the relation between microscopic structure and phonon propagation at the nan oscale is essential to design materials with desired phononic and thermal properties.Here we uncover a new mechanism of phonon interaction in surface oxidized membranes, i.e., native oxide layers interact with phonons in ultra-thin silicon membranes through local resonances. The local resonances reduce the low frequency phonon group velocities and shorten their mean free path. This effect opens up a new strategy for ultralow thermal conductivity design as it complements the scattering mechanism which scatters higher frequency modes effectively. The combination of native oxide layer and alloying with germanium in concentration as small as 5% reduces the thermal conductivity of silicon membranes to 100 time lower than the bulk. In addition, the resonance mechanism produced by native oxide surface layers is particularly effective for thermal condutivity reduction even at very low temperatures, at which only low frequency modes are populated.


title: A systematic study of mass spectra and strong decay of strange mesons

abstract: The mass spectrum of the kaon family is analyzed by the {modified} Godfrey-Isgur model with a color screening effect approximating the kaon as a heavy-light meson system. This analysis gives us the structure and possible assignments of the observed kaon candidates, which can be tested by comparing the theoretical results of their two-body strong decays with the experimental data. Additionally, prediction of some partial decay widths is made on the kaons still missing in experiment. This study is crucial to establishing the kaon family and searching for their higher excitations in the future.


title: Cyclic period oscillation of the eclipsing dwarf nova DV UMa

abstract: DV UMa is an eclipsing dwarf nova with an orbital period of $\sim2.06$ h, which lies just at the bottom edge of the period gap. To detect its orbital period changes we present 12 new mid-eclipse times by using our CCD photometric data and archival data. Combining with the published mid-eclipse times in quiescence, spanning $\sim30$ yr, the latest version of the $O-C$ diagram was obtained and analyzed. The best fit to those available eclipse timings shows that the orbital period of DV UMa is undergoing a cyclic oscillation with a period of $17.58(\pm0.52)$ yr and an amplitude of $71.1(\pm6.7)$ s. The periodic variation most likely arises from the light-travel-time effect via the presence of a circumbinary object because the required energy to drive the Applegate mechanism is too high in this system. The mass of the unseen companion was derived as $M_{3}\sin{i'}=0.025(\pm0.004)M_{\odot}$. If the third body is in the orbital plane (i.e. $i'=i=82.9^{\circ}$) of the eclipsing pair, it would match to a brown dwarf. This hypothetical brown dwarf is orbiting its host star at a separation of $\sim8.6$ AU in an eccentric orbit ($e=0.44$).


title: CHAM: action recognition using convolutional hierarchical attention model

abstract: Recently, the soft attention mechanism, which was originally proposed in language processing, has been applied in computer vision tasks like image captioning. This paper presents improvements to the soft attention model by combining a convolutional LSTM with a hierarchical system architecture to recognize action categories in videos. We call this model the Convolutional Hierarchical Attention Model (CHAM). The model applies a convolutional operation inside the LSTM cell and an attention map generation process to recognize actions. The hierarchical architecture of this model is able to explicitly reason on multi-granularities of action categories. The proposed architecture achieved improved results on three publicly available datasets: the UCF sports dataset, the Olympic sports dataset and the HMDB51 dataset.


title: ALMA Discovery of Solar Umbral Brightness Enhancement at {\lambda}=3 mm

abstract: We report the discovery of a brightness enhancement in the center of a large sunspot umbra at a wavelength of 3 mm using the Atacama Large Millimeter/sub-millimeter Array (ALMA). Sunspots are amongst the most prominent features on the solar surface, but many of their aspects are surprisingly poorly understood. We analyzed a {\lambda}=3 mm (100 GHz) mosaic image obtained by ALMA, which includes a large sunspot within the active region AR12470 on December 16, 2015. The 3 mm map has a field-of-view and spatial resolution, which is the highest spatial-resolution map of an entire sunspot in this frequency range. We find a gradient of 3 mm brightness from a high value in the outer penumbra to a low value in the inner penumbra/outer umbra. Within the inner umbra, there is a marked increase in 3mm brightness temperature, which we call an umbral brightness enhancement. This enhanced emission corresponds to a temperature excess of 800 K relative to the surrounding inner penumbral region and coincides with excess brightness in the 1330 and 1400 {\AA} slitjaw images of the Interface Region Imaging Spectrograph (IRIS), adjacent to a partial lightbridge. This {\lambda}=3 mm brightness enhancement may be an intrinsic feature of the sunspot umbra at chromospheric heights, such as a manifestation of umbral flashes, or it could be related to a coronal plume since the brightness enhancement was coincident with the footpoint of a coronal loop observed at 171 {\AA}.


title: Deep Spatio-temporal Manifold Network for Action Recognition

abstract: Visual data such as videos are often sampled from complex manifold. We propose leveraging the manifold structure to constrain the deep action feature learning, thereby minimizing the intra-class variations in the feature space and alleviating the over-fitting problem. Considering that manifold can be transferred, layer by layer, from the data domain to the deep features, the manifold priori is posed from the top layer into the back propagation learning procedure of convolutional neural network (CNN). The resulting algorithm --Spatio-Temporal Manifold Network-- is solved with the efficient Alternating Direction Method of Multipliers and Backward Propagation (ADMM-BP). We theoretically show that STMN recasts the problem as projection over the manifold via an embedding method. The proposed approach is evaluated on two benchmark datasets, showing significant improvements to the baselines.


title: Determination of spin relaxation times in heavy metals via 2nd harmonic spin injection magnetoresistance

abstract: In tunnel junctions between ferromagnets and heavy elements with strong spin orbit coupling the magnetoresistance is often dominated by tunneling anisotropic magnetoresistance (TAMR). This makes conventional DC spin injection techniques impractical for determining the spin relaxation time ($\tau_s$). Here, we show that this obstacle for measurements of $\tau_s$ can be overcome by 2nd harmonic spin-injection-magnetoresistance (SIMR). In the 2nd harmonic signal the SIMR is comparable in magnitude to TAMR, thus enabling Hanle-induced SIMR as a powerful tool to directly determine $\tau_s$. Using this approach we determined the spin relaxation time of Pt and Ta and their temperature dependences. The spin relaxation in Pt seems to be governed by Elliott-Yafet mechanism due to a constant resistivity $\times$spin relaxation time product over a wide temperature range.


title: Binary de Bruijn Sequences via Zech's Logarithms

abstract: The focus of this work is to show how to combine Zech's logarithms and each of the cycle joining and cross-join pairing methods to construct binary de Bruijn sequences of any order. A basic implementation is supplied as a proof-of-concept. The cycles, in the cycle joining method, are typically generated by a linear feedback shift register. We prove a crucial characterization that determining Zech's logarithms is equivalent to identifying conjugate pairs shared by any two distinct cycles. This speeds up the task of building a connected adjacency subgraph that contains all vertices of the complete adjacency graph. Distinct spanning trees in either graph correspond to cyclically inequivalent de Bruijn sequences. As the cycles are being joined, guided by the conjugate pairs, we track the changes in the feedback function. Certificates of star or almost-star spanning trees conveniently handle large order cases. The characterization of conjugate pairs via Zech's logarithms, as positional markings, is then adapted to identify cross-join pairs. A modified $m$-sequence is initially used, for ease of generation. The process can be repeated on each of the resulting de Bruijn sequences. We show how to integrate an analytic tool, attributed to Fryers, in the process. Most prior constructions in the literature measure the complexity of the corresponding bit-by-bit algorithms. Our approach is different. We aim first to build a connected adjacency subgraph that is certified to contain all of the cycles as vertices. The ingredients are computed just once and concisely stored. Simple strategies are offered to keep the complexities low as the order grows.


title: Phonetic Temporal Neural Model for Language Identification

abstract: Deep neural models, particularly the LSTM-RNN model, have shown great potential for language identification (LID). However, the use of phonetic information has been largely overlooked by most existing neural LID methods, although this information has been used very successfully in conventional phonetic LID systems. We present a phonetic temporal neural model for LID, which is an LSTM-RNN LID system that accepts phonetic features produced by a phone-discriminative DNN as the input, rather than raw acoustic features. This new model is similar to traditional phonetic LID methods, but the phonetic knowledge here is much richer: it is at the frame level and involves compacted information of all phones. Our experiments conducted on the Babel database and the AP16-OLR database demonstrate that the temporal phonetic neural approach is very effective, and significantly outperforms existing acoustic neural models. It also outperforms the conventional i-vector approach on short utterances and in noisy conditions.


title: Phone-aware Neural Language Identification

abstract: Pure acoustic neural models, particularly the LSTM-RNN model, have shown great potential in language identification (LID). However, the phonetic information has been largely overlooked by most of existing neural LID models, although this information has been used in the conventional phonetic LID systems with a great success. We present a phone-aware neural LID architecture, which is a deep LSTM-RNN LID system but accepts output from an RNN-based ASR system. By utilizing the phonetic knowledge, the LID performance can be significantly improved. Interestingly, even if the test language is not involved in the ASR training, the phonetic knowledge still presents a large contribution. Our experiments conducted on four languages within the Babel corpus demonstrated that the phone-aware approach is highly effective.


title: On GMRES for singular EP and GP systems

abstract: In this contribution, we study the numerical behavior of the Generalized Minimal Residual (GMRES) method for solving singular linear systems. It is known that GMRES determines a least squares solution without breakdown if the coefficient matrix is range-symmetric (EP), or if its range and nullspace are disjoint (GP) and the system is consistent. We show that the accuracy of GMRES iterates may deteriorate in practice due to three distinct factors: (i) the inconsistency of the linear system; (ii) the distance of the initial residual to the nullspace of the coefficient matrix; (iii) the extremal principal angles between the ranges of the coefficient matrix and its transpose. These factors lead to poor conditioning of the extended Hessenberg matrix in the Arnoldi decomposition and affect the accuracy of the computed least squares solution. We also compare GMRES with the range restricted GMRES (RR-GMRES) method. Numerical experiments show typical behaviors of GMRES for small problems with EP and GP matrices.


title: Cultural Values and Cross-cultural Video Consumption on YouTube

abstract: Video-sharing social media like YouTube provide access to diverse cultural products from all over the world, making it possible to test theories that the Web facilitates global cultural convergence. Drawing on a daily listing of YouTube's most popular videos across 58 countries, we investigate the consumption of popular videos in countries that differ in cultural values, language, gross domestic product, and Internet penetration rate. Although online social media facilitate global access to cultural products, we find this technological capability does not result in universal cultural convergence. Instead, consumption of popular videos in culturally different countries appears to be constrained by cultural values. Cross-cultural convergence is more advanced in cosmopolitan countries with cultural values that favor individualism and power inequality.


title: WZ Sge: an eclipsing cataclysmic variable evolving towards the period minimum

abstract: We present the photometric results of the eclipsing cataclysmic variable (CV) WZ Sge near the period minimum ($P_{min}$). Eight new mid-eclipse times were determined and the orbital ephemeris was updated. Our result shows that the orbital period of WZ Sge is decreasing at a rate of $\dot{P}=-2.72(\pm0.23)\times{10^{-13}}\,s s^{-1}$. This secular decrease, coupled with previous detection of its donor, suggest that WZ Sge is a pre-bounce system. Further analysis indicates that the observed period decrease rate is about $1.53$ times higher than pure gravitational radiation (GR) driving. We constructed the evolutionary track of WZ Sge, which predicts that $P_{min}$ of WZ Sge is $\sim77.98 (\pm0.90)$ min. If the orbital period decreases at the current rate, WZ Sge will evolve past its $P_{min}$ after $\sim25.3$ Myr. Based on the period evolution equation we find $\dot{M}_{2}\simeq4.04(\pm0.10)\times10^{-11}M_{\odot}yr^{-1}$, which is compatible with the current concept of CV evolution at ultrashort orbital periods.


title: Entropic repulsion and lack of the $g$-measure property for Dyson models

abstract: We consider Dyson models, Ising models with slow polynomial decay, at low temperature and show that its Gibbs measures deep in the phase transition region are not $g$-measures. The main ingredient in the proof is the occurrence of an entropic repulsion effect, which follows from the mesoscopic stability of a (single-point) interface for these long-range models in the phase transition region.

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


title: The number of eigenvalues of the matrix Schr\"odinger operator on the half line with general boundary conditions

abstract: We prove a bound, of Bargmann- Birman-Schwinger type, on the number of eigenvalues of the matrix Schr\"odinger operator on the half line, with the most general self adjoint boundary condition at the origin, and with selfadjoint matrix potentials that are integrable and have a finite first moment.


title: Manifold Relevance Determination: Learning the Latent Space of Robotics

abstract: In this article we present the basics of manifold relevance determination (MRD) as introduced in \cite{mrd}, and some applications where the technology might be of particular use. Section 1 acts as a short tutorial of the ideas developed in \cite{mrd}, while Section 2 presents possible applications in sensor fusion, multi-agent SLAM, and "human-appropriate" robot movement (e.g. legibility and predictability~\cite{dragan-hri-2013}). In particular, we show how MRD can be used to construct the underlying models in a data driven manner, rather than directly leveraging first principles theories (e.g., physics, psychology) as is commonly the case for sensor fusion, SLAM, and human robot interaction. We note that [Bekiroglu et al., 2016] leveraged MRD for correcting unstable robot grasps to stable robot grasps.


title: Contour Detection from Deep Patch-level Boundary Prediction

abstract: In this paper, we present a novel approach for contour detection with Convolutional Neural Networks. A multi-scale CNN learning framework is designed to automatically learn the most relevant features for contour patch detection. Our method uses patch-level measurements to create contour maps with overlapping patches. We show the proposed CNN is able to to detect large-scale contours in an image efficienly. We further propose a guided filtering method to refine the contour maps produced from large-scale contours. Experimental results on the major contour benchmark databases demonstrate the effectiveness of the proposed technique. We show our method can achieve good detection of both fine-scale and large-scale contours.


title: Physical parameters and orbital period variation of a newly discovered cataclysmic variable GSC 4560-02157

abstract: GSC 4560-02157 is a new eclipsing cataclysmic variable with an orbital period of $0.265359$ days. By using the published $V-$ and $R-$band data together with our observations, we discovered that the $O-C$ curve of GSC 4560-02157 may shows a cyclic variation with the period of $3.51$ years and an amplitude of $1.40$ min. If this variation is caused by a light travel-time effect via the existence of a third body, its mass can be derived as $M_{3}sini'\approx91.08M_{Jup}$, it should be a low-mass star. In addition, several physical parameters were measured. The colour of the secondary star was determined as $V-R=0.77(\pm0.03)$ which corresponds to a spectral type of K2-3. The secondary star's mass was estimated as $M_{2}=0.73(\pm0.02)M_{\odot}$ by combing the derived $V-R$ value around phase 0 with the assumption that it obeys the mass-luminosity relation of the main sequence stars. This mass is consistent with the mass$-$period relation of CV donor stars. For the white dwarf, the eclipse durations and contacts of the white dwarf yield an upper limit of the white dwarf's radius corresponding to a lower limit mass of $M_{1}\approx0.501M_{\odot}$. The overestimated radius and previously published spectral data indicate that the boundary layer may has a very high temperature.


title: Violating the Quantum Focusing Conjecture and Quantum Covariant Entropy Bound in $d\ge 5$ dimensions

abstract: We study the Quantum Focussing Conjecture (QFC) in curved spacetime. Noting that quantum corrections from integrating out massive fields generally induce a Gauss-Bonnet term, we study Einstein-Hilbert-Gauss-Bonnet gravity and show for $d\ge 5$ spacetime dimensions that weakly-curved solutions can violate the associated QFC for either sign of the Gauss-Bonnet coupling. The nature of the violation shows that -- so long as the Gauss-Bonnet coupling is non-zero -- it will continue to arise for local effective actions containing arbitrary further higher curvature terms, and when gravity is coupled to generic $d\ge 5$ theories of massive quantum fields. The argument also implies violations of a recently-conjectured form of the generalized covariant entropy bound. The possible validity of the QFC and covariant entropy bound in $d\le 4$ spacetime dimensions remains open.


title: Accelerating solutions of one-dimensional unsteady PDEs with GPU-based swept time-space decomposition

abstract: The expedient design of precision components in aerospace and other high-tech industries requires simulations of physical phenomena often described by partial differential equations (PDEs) without exact solutions. Modern design problems require simulations with a level of resolution difficult to achieve in reasonable amounts of time---even in effectively parallelized solvers. Though the scale of the problem relative to available computing power is the greatest impediment to accelerating these applications, significant performance gains can be achieved through careful attention to the details of memory communication and access. The swept time-space decomposition rule reduces communication between sub-domains by exhausting the domain of influence before communicating boundary values. Here we present a GPU implementation of the swept rule, which modifies the algorithm for improved performance on this processing architecture by prioritizing use of private (shared) memory, avoiding interblock communication, and overwriting unnecessary values. It shows significant improvement in the execution time of finite-difference solvers for one-dimensional unsteady PDEs, producing speedups of 2--9$\times$ for a range of problem sizes, respectively, compared with simple GPU versions and 7--300$\times$ compared with parallel CPU versions. However, for a more sophisticated one-dimensional system of equations discretized with a second-order finite-volume scheme, the swept rule performs 1.2--1.9$\times$ worse than a standard implementation for all problem sizes.


title: Long-term photometric behaviour of the eclipsing Z Cam-type dwarf nova AY Psc

abstract: New eclipse timings of the Z Cam-type dwarf nova AY Psc were measured and the orbital ephemeris was revised. Based on the long-term AAVSO data, moreover, the outburst behaviors were also explored. Our analysis suggests that the normal outbursts are quasi-periodic, with an amplitude of $\sim2.5(\pm0.1)$ mag and a period of $\sim18.3(\pm0.7)$ days. The amplitude vs. recurrence-time relation of AY Psc is discussed, and we concluded that this relation may represents general properties of dwarf nova (DN) outbursts. The observed standstill ends with an outburst, which is inconsistent with the general picture of Z Cam-type stars. This unusual behavior was considered to be related to the mass-transfer outbursts. Moreover, the average luminosity is brighter during standstills than during outburst cycles. The changes in brightness marks the variations in $\dot{M}_{2}$ due to the disc of AY Psc is nearly steady state. $\dot{M}_{2}$ value was limited to the range from $6.35\times10^{-9}$ to $1.18\times10^{-8}M_{\odot}yr^{-1}$. More detailed examination shows that there are a few small outbursts presence during standstills. These events with amplitudes of $\sim0.5-0.9$ mag are very similar to the stunted outbursts reported in some NLs. We discussed several possible mechanisms and suggested that the most reasonable mechanism for these stunted outbursts is a changing mass-transfer rate.


title: Long-term photometric behavior of the eclipsing cataclysmic variable V729 Sgr

abstract: We present the analysis results of an eclipsing cataclysmic variable (CV) V729 Sgr, based on our observations and AAVSO data. Some outburst parameters were determined such as outburst amplitude ($A_{n}$) and recurrence time ($T_{n}$), and then the relationship between $A_{n}$ and $T_{n}$ is discussed. A cursory examination for the long-term light curves reveals that there are small-amplitude outbursts and dips present, which is similar to the behaviors seen in some nova-like CVs (NLs). More detailed inspection suggests that the outbursts in V729 Sgr may be Type A (outside-in) with a rise time $\sim1.76$ d. Further analysis also shows that V729 Sgr is an intermediate between dwarf nova and NLs, and we constrain its mass transfer rate to $1.59\times10^{-9} < \dot{M}_{2} < 5.8\times10^{-9}M_{\odot}yr^{-1}$ by combining the theory for Z Cam type stars with observations. Moreover, the rapid oscillations in V729 Sgr were detected and analyzed for the first time. Our results indicate that the oscillation at $\sim 25.5$ s is a true DNO, being associated with the accretion events. The classification of the oscillations at $\sim 136$ and $154$ s as lpDNOs is based on the relation between $P_{lpDNOs}$ and $P_{DNOs}$. Meanwhile, the QPOs at the period of hundreds of seconds are also detected.


title: Phase diagram of the sulfur hydride: transition into high Tc state

abstract: Remarkable feature of the phase diagram of sulfur hydrides is a sharp increase in critical temperature from 120K up to200K. This increase is a signature of the structural transition. The study described below is concerned with the nature of this phase transition. One can demonstrate that the symmetry analysis along with an analysis of the impact of lattice deformations lead to the conclusion that we are dealing with the first order transition. It leads to an abrupt appearance of small pockets on the Fermi surface and to the two-gap energy spectrum.


title: Double cyclic variations in orbital period of the eclipsing cataclysmic variable EX Dra

abstract: EX Dra is a long-period eclipsing dwarf nova with $\sim2-3$ mag amplitude outbursts. This star has been monitored photometrically from November, 2009 to March, 2016 and 29 new mid-eclipse times were obtained. By using new data together with the published data, the best fit to the $O-C$ curve indicate that the orbital period of EX Dra have an upward parabolic change while undergoing double-cyclic variations with the periods of 21.4 and 3.99 years, respectively. The upward parabolic change reveals a long-term increase at a rate of $\dot{P}={+7.46}\times10^{-11}{s} {s^{-1}}$. The evolutionary theory of cataclysmic variables (CVs) predicts that, as a CV evolves, the orbital period should be decreasing rather than increasing. Secular increase can be explained as the mass transfer between the secondary and primary or may be just an observed part of a longer cyclic change. Most plausible explanation for the double-cyclic variations is a pair of light travel-time effect via the presence of two companions. Their masses are determined to be $M_{A}sini'_{A}=29.3(\pm0.6) M_{Jup}$ and $M_{B}sini'_{B}=50.8(\pm0.2) M_{Jup}$. When the two companions are coplanar to the orbital plane of the central eclipsing pair, their masses would match to brown dwarfs.


title: Solving Constrained Horn Clauses Using Dependence-Disjoint Expansions

abstract: Recursion-free Constrained Horn Clauses (CHCs) are logic-programming problems that can model safety properties of programs with bounded iteration and recursion. In addition, many CHC solvers reduce recursive systems to a series of recursion-free CHC systems that can each be solved efficiently. In this paper, we define a novel class of recursion-free systems, named Clause-Dependence Disjoint (CDD), that generalizes classes defined in previous work. The advantage of this class is that many CDD systems are smaller than systems which express the same constraints but are part of a different class. This advantage in size allows CDD systems to be solved more efficiently than their counterparts in other classes. We implemented a CHC solver named Shara. Shara solves arbitrary CHC systems by reducing the input to a series of CDD systems. Our evaluation indicates that Shara outperforms state-of-the-art implementations in many practical cases.


title: Shortcuts to adiabaticity in the infinite-range Ising model by mean-field counter-diabatic driving

abstract: The strategy of shortcuts to adiabaticity enables us to realize adiabatic dynamics in finite time. In the counter-diabatic driving approach, an auxiliary Hamiltonian which is called the counter-diabatic Hamiltonian is appended to an original Hamiltonian to cancel out diabatic transitions. The counter-diabatic Hamiltonian is constructed by using the eigenstates of the original Hamiltonian. Therefore, it is in general difficult to construct the counter-diabatic Hamiltonian for quantum many-body systems. Even if the counter-diabatic Hamiltonian for quantum many-body systems is obtained, it is generally non-local and even diverges at critical points. We construct an approximated counter-diabatic Hamiltonian for the infinite-range Ising model by making use of the mean-field approximation. An advantage of this method is that the mean-field counter-diabatic Hamiltonian is constructed by only local operators. We numerically demonstrate the effectiveness of this method through quantum annealing processes going the vicinity of the critical point. It is also confirmed that the mean-field counter-diabatic Hamiltonian is still well-defined in the limit to the critical point. The present method can take higher order contributions into account and is consistent with the variational approach for local counter-diabatic driving.


title: Quasi-nonlinear approach to the Weibel instability

abstract: Astrophysical and high-energy-density laboratory plasmas often have large-amplitude, sub-Larmor-scale electromagnetic fluctuations excited by various kinetic-streaming or anisotropy-driven instabilities. The Weibel (or the filamentation) instability is particularly important because it can rapidly generate strong magnetic fields, even in the absence of seed fields. Particles propagating in collisionless plasmas with such small-scale magnetic fields undergo stochastic deflections similar to Coulomb collisions, with the magnetic pitch-angle diffusion coefficient representing the effective "collision" frequency. We show that this effect of the plasma "quasi-collisionality" can strongly affect the growth rate and evolution of the Weibel instability in the deeply nonlinear regime. This result is especially important for understanding cosmic-ray-driven turbulence in an upstream region of a collisionless shock of a gamma-ray burst or a supernova. We demonstrate that the quasi-collisions caused by the fields generated in the upstream suppress the instability slightly but can never shut it down completely. This confirms the assumptions made in the self-similar model of the collisionless foreshock.


title: Indistinguishable Encoding for Bidirectional Quantum Key Distribution: Theory to Experiment

abstract: We present for the first time, a bidirectional Quantum Key Distribution protocol with minimal encoding operations derived from the use of two `nonorthogonal' unitary transformations selected from two mutually unbiased unitary bases; which are indistinguishable in principle for a single use. Along with its decoding procedure, it is a stark contrast to its `orthogonal encoding' predecessors. Defining a more relevant notion of security threshold for such protocols, the current protocol outperforms its predecessor in terms of security as the maximal amount of information an eavesdropper can glean is essentially limited by the indistinguishability of the transformations. We further propose adaptations for a practical scenario and report on a proof of concept experimental scheme based on polarised photons from an attenuated pulsed laser for qubits, demonstrating the feasibility of such a protocol.


title: Lazarsfeld-Mukai Reflexive Sheaves and their Stability

abstract: Consider an ample and globally generated line bundle $L$ on a smooth projective variety $X$ of dimension $N\geq 2$ over $\mathbb{C}$. Let $D$ be a smooth divisor in the complete linear system of $L$. We construct reflexive sheaves on $X$ by an elementary transformation of a trivial bundle on $X$ along certain globally generated torsion-free sheaves on $D$. The dual reflexive sheaves are called the Lazarsfeld-Mukai reflexive sheaves. We prove the $\mu_L$-(semi)stability of such reflexive sheaves under certain conditions.


title: A Survey of Location Prediction on Twitter

abstract: Locations, e.g., countries, states, cities, and point-of-interests, are central to news, emergency events, and people's daily lives. Automatic identification of locations associated with or mentioned in documents has been explored for decades. As one of the most popular online social network platforms, Twitter has attracted a large number of users who send millions of tweets on daily basis. Due to the world-wide coverage of its users and real-time freshness of tweets, location prediction on Twitter has gained significant attention in recent years. Research efforts are spent on dealing with new challenges and opportunities brought by the noisy, short, and context-rich nature of tweets. In this survey, we aim at offering an overall picture of location prediction on Twitter. Specifically, we concentrate on the prediction of user home locations, tweet locations, and mentioned locations. We first define the three tasks and review the evaluation metrics. By summarizing Twitter network, tweet content, and tweet context as potential inputs, we then structurally highlight how the problems depend on these inputs. Each dependency is illustrated by a comprehensive review of the corresponding strategies adopted in state-of-the-art approaches. In addition, we also briefly review two related problems, i.e., semantic location prediction and point-of-interest recommendation. Finally, we list future research directions.


title: Star Formation in Simulated Galaxies: Understanding the Transition to Quiescence at $3\times10^{10}$ M$_\odot$

abstract: Star formation in galaxies relies on the availability of cold, dense gas, which, in turn, relies on factors internal and external to the galaxies. In order to provide a simple model for how star formation is regulated by various physical processes in galaxies, we analyse data at redshift $z=0$ from a hydrodynamical cosmological simulation that includes prescriptions for star formation and stellar evolution, active galactic nuclei (AGN), and their associated feedback processes. This model can determine the star formation rate (SFR) as a function of galaxy stellar mass, gas mass, black hole mass, and environment. We find that gas mass is the most important quantity controlling star formation in low-mass galaxies, and star-forming galaxies in dense environments have higher SFR than their counterparts in the field. In high-mass galaxies, we find that black holes more massive than $\sim10^{7.5}$ M$_\odot$ can be triggered to quench star formation in their host; this mass scale is emergent in our simulations. Furthermore, this black hole mass corresponds to a galaxy bulge mass $\sim2\times10^{10}$ M$_\odot$, consistent with the mass at which galaxies start to become dominated by early types ($\sim3\times10^{10}$ M$_\odot$, as previously shown in observations by Kauffmann et al.). Finally, we demonstrate that our model can reproduce well the SFR measured from observations of galaxies in the GAMA and ALFALFA surveys.


title: Pyramids and 2-representations

abstract: We describe a diagrammatic procedure which lifts strict monoidal actions from additive categories to categories of complexes avoiding any use of direct sums. As an application, we prove that every simple transitive $2$-representation of the $2$-category of projective bimodules over a finite dimensional algebra is equivalent to a cell $2$-representation.


title: Emotional Metaheuristics For in-situ Foraging Using Sensor Constrained Robot Swarms

abstract: We present a new social animal inspired emotional swarm intelligence technique. This technique is used to solve a variant of the popular collective robots problem called foraging. We show with a simulation study how simple interaction rules based on sensations like hunger and loneliness can lead to globally coherent emergent behavior which allows sensor constrained robots to solve the given problem