title: Algorithms for Scheduling Malleable Cloud Tasks

abstract: Due to the ubiquity of batch data processing in cloud computing, the related problem of scheduling malleable batch tasks and its extensions have received significant attention recently. In this paper, we consider a fundamental model where a set of n tasks is to be processed on C identical machines and each task is specified by a value, a workload, a deadline and a parallelism bound. Within the parallelism bound, the number of machines assigned to a task can vary over time without affecting its workload. For this model, we obtain two core results: a sufficient and necessary condition such that a set of tasks can be finished by their deadlines on C machines, and an algorithm to produce such a schedule. These core results provide a conceptual tool and an optimal scheduling algorithm that enable proposing new algorithmic analysis and design and improving existing algorithms under various objectives.


title: Asymptotic bounds of depth for a reversible circuit consisting of NOT, CNOT and 2-CNOT gates

abstract: The paper discusses the asymptotic depth of a reversible circuit consisting of NOT, CNOT and 2-CNOT gates. Reversible circuit depth function $D(n, q)$ for a circuit implementing a transformation $f\colon \mathbb Z_2^n \to \mathbb Z_2^n$ is introduced as a function of $n$ and the number of additional inputs $q$. It is proved that for the case of implementing a permutation from $A(\mathbb Z_2^n)$ with a reversible circuit having no additional inputs the depth is bounded as $D(n, 0) \gtrsim 2^n / (3\log_2 n)$. It is proved that for the case of implementing a transformation $f\colon \mathbb Z_2^n \to \mathbb Z_2^n$ with a reversible circuit having $q_0 \sim 2^n$ additional inputs the depth is bounded as $D(n, q_0) \lesssim 3n$.


title: Efficient fourth order symplectic integrators for near-harmonic separable Hamiltonian systems

abstract: Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of the simple harmonic oscillator. The performance of the methods is evaluated for various Hamiltonian systems: Integration errors are compared to those of acclaimed integrators composed by S. Blanes et al. (2013), W. Kahan et al. (1999) and H. Yoshida (1990). Numerical tests indicate that the integrators obtained in this paper perform significantly better than previous integrators for common Hamiltonian systems.


title: Mathematical Language Processing: Automatic Grading and Feedback for Open Response Mathematical Questions

abstract: While computer and communication technologies have provided effective means to scale up many aspects of education, the submission and grading of assessments such as homework assignments and tests remains a weak link. In this paper, we study the problem of automatically grading the kinds of open response mathematical questions that figure prominently in STEM (science, technology, engineering, and mathematics) courses. Our data-driven framework for mathematical language processing (MLP) leverages solution data from a large number of learners to evaluate the correctness of their solutions, assign partial-credit scores, and provide feedback to each learner on the likely locations of any errors. MLP takes inspiration from the success of natural language processing for text data and comprises three main steps. First, we convert each solution to an open response mathematical question into a series of numerical features. Second, we cluster the features from several solutions to uncover the structures of correct, partially correct, and incorrect solutions. We develop two different clustering approaches, one that leverages generic clustering algorithms and one based on Bayesian nonparametrics. Third, we automatically grade the remaining (potentially large number of) solutions based on their assigned cluster and one instructor-provided grade per cluster. As a bonus, we can track the cluster assignment of each step of a multistep solution and determine when it departs from a cluster of correct solutions, which enables us to indicate the likely locations of errors to learners. We test and validate MLP on real-world MOOC data to demonstrate how it can substantially reduce the human effort required in large-scale educational platforms.


title: Quantum criticality of the two-channel pseudogap Anderson model: Universal scaling in linear and non-linear conductance

abstract: The quantum criticality of the two-lead two-channel pseudogap Anderson model is studied. Based on the non-crossing approximation, we calculate both the linear and nonlinear conductance of the model at finite temperatures with a voltage bias and a power-law vanishing conduction electron density of states, $\propto |\omega-\mu_F|^r$ ($0<r<1$) near the Fermi energy. Equilibrium and non-equilibrium quantum critical properties at the two-channel Kondo (2CK) to local moment (LM) phase transition are addressed by extracting universal scaling functions in both linear and non-linear conductances, respectively. Clear distinctions are found on the critical exponents between linear and non-linear conductance. The implications of these two distinct quantum critical properties for the non-equilibrium quantum criticality in general are discussed.


title: The cost of attack in competing networks

abstract: Real-world attacks can be interpreted as the result of competitive interactions between networks, ranging from predator-prey networks to networks of countries under economic sanctions. Although the purpose of an attack is to damage a target network, it also curtails the ability of the attacker, which must choose the duration and magnitude of an attack to avoid negative impacts on its own functioning. Nevertheless, despite the large number of studies on interconnected networks, the consequences of initiating an attack have never been studied. Here, we address this issue by introducing a model of network competition where a resilient network is willing to partially weaken its own resilience in order to more severely damage a less resilient competitor. The attacking network can take over the competitor nodes after their long inactivity. However, due to a feedback mechanism the takeovers weaken the resilience of the attacking network. We define a conservation law that relates the feedback mechanism to the resilience dynamics for two competing networks. Within this formalism, we determine the cost and optimal duration of an attack, allowing a network to evaluate the risk of initiating hostilities.


title: Quantum Logic with Interacting Bosons in 1D

abstract: We present a scheme for implementing high-fidelity quantum logic gates using the quantum walk of a few interacting bosons on a one-dimensional lattice. The gate operation is carried out by a single compact lattice described by a one-dimensional Bose-Hubbard model with only nearest-neighbor hopping and on-site interactions. We find high-fidelity deterministic logic operations for a gate set (including the CNOT gate) that is universal for quantum information processing. We discuss the applicability of this scheme in light of recent developments in controlling and monitoring cold-atoms in optical lattices, as well as an implementation with realistic nonlinear quantum photonic devices.


title: Ultra-high mechanical stretchability and controllable topological phase transitions in two-dimensional arsenic

abstract: The mechanical stretchability is the magnitude of strain which a material can suffer before it breaks. Materials with high mechanical stretchability, which can reversibly withstand extreme mechanical deformation and cover arbitrary surfaces and movable parts, are used for stretchable display devices, broadband photonic tuning and aberration-free optical imaging. Strain can be utilised to control the band structures of materials and can even be utilised to induce a topological phase transition, driving the normal insulators to topological non-trivial materials with non-zero Chern number or Z2 number. Here, we propose a new two-dimensional topological material with ultra-high mechanical stretchability - the ditch-like 2D arsenic. This new anisotropic material possesses a large Poisson's ratio 1.049, which is larger than any other reported inorganic materials and has a ultra-high stretchability 44% along the armchair direction, which is unprecedent in inorganic materials as far as we know. Its minimum bend radius of this material can be as low as 0.66 nm, which is comparable to the radius of carbon-nanotube. Such mechanical properties make this new material be a stretchable semiconductor which could be used to construct flexible display devices and stretchable sensors. Axial strain will make a conspicuous affect on the band structure of the system, and a proper strain along the zigzag direction will drive the 2D arsenic into the topological insulator in which the topological edge state can host dissipation-less spin current and spin transfer toque, which are useful in spintronics devices such as dissipation transistor, interconnect channels and spin valve devices.


title: Bayesian derivation of plasma equilibrium distribution function for tokamak scenarios and the associated Landau collision operator

abstract: A class of parametric distribution functions has been proposed in [C.DiTroia, Plasma Physics and Controlled Fusion,54,2012] as equilibrium distribution functions (EDFs) for charged particles in fusion plasmas, representing supra-thermal particles in anisotropic equilibria for Neutral Beam Injection, Ion Cyclotron Heating scenarios. Moreover, the EDFs can also represent nearly isotropic equilibria for Slowing-Down $alpha$ particles and core thermal plasma populations. These EDFs depend on constants of motion (COMs). Assuming an axisymmetric system with no equilibrium electric field, the EDF depends on the toroidal canonical momentum $P_\phi$, the kinetic energy $w$ and the magnetic moment \mu. In the present work, the EDFs are obtained from first principles and general hypothesis. The derivation is probabilistic and makes use of the Bayes' Theorem. The bayesian argument allows us to describe how far from the prior probability distribution function (pdf), e.g. Maxwellian, the plasma is, based on the information obtained from magnetic moment and GC velocity pdf. Once the general functional form of the EDF has been settled, it is shown how to associate a Landau collision operator and a Fokker-Planck equation that ensures the system relaxation towards the proposed EDF.


title: Solving the Infinite-horizon Constrained LQR Problem using Accelerated Dual Proximal Methods

abstract: This work presents an algorithmic scheme for solving the infinite-time constrained linear quadratic regulation problem. We employ an accelerated version of a popular proximal gradient scheme, commonly known as the Forward-Backward Splitting (FBS), and prove its convergence to the optimal solution in our infinite-dimensional setting. Each iteration of the algorithm requires only finite memory, is computationally cheap, and makes no use of terminal invariant sets; hence, the algorithm can be applied to systems of very large dimensions. The acceleration brings in optimal convergence rates O(1/k^2) for function values and O(1/k) for primal iterates and renders the proposed method a practical alternative to model predictive control schemes for setpoint tracking. In addition, for the case when the true system is subject to disturbances or modelling errors, we propose an efficient warm-starting procedure, which significantly reduces the number of iterations when the algorithm is applied in closed-loop. Numerical examples demonstrate the approach.


title: From charge motion in general magnetic fields to the non perturbative gyrokinetic equation

abstract: The exact analytical description of non relativistic charge motion in general magnetic fields is, apparently, a simple problem but, it has not been solved up to now apart for rare cases. The key feature of the present derivation is to adopt a non perturbative magnetic field description to find new solutions of motion. Among all solutions, two are particularly important: guiding particle and gyro-particle solutions. The guiding particle has been characterized to be minimally coupled to the magnetic field, the gyro-particle has been defined to be maximally coupled to the magnetic field and, also, to move on a closed orbit. The generic charged particle motion is shown to be expressed as the sum of such particular solutions. This non perturbative approach corresponds to the description of the particle motion in the gyro-center and/or guiding center reference frame obtained at all the orders of the modern gyro-center transformation. The Boltzmann equation is analyzed with the described exact guiding center coordinates. The obtained gyrokinetic equation is solved for the Boltzmann equation at marginal stability conditions.


title: Coinduction: an elementary approach

abstract: The main aim of this paper is to promote a certain style of doing coinductive proofs, similar to inductive proofs as commonly done by mathematicians. For this purpose, we provide a reasonably direct justification for coinductive proofs written in this style, i.e., converting a coinductive proof into a non-coinductive argument is purely a matter of routine. In this way, we provide an elementary explanation of how to interpret coinduction in set theory.


title: Warped AdS$_3$, dS$_3$ and flows from $\mathcal{N} = (0,2)$ SCFTs

abstract: We present the general form of all timelike supersymmetric solutions to 3D U(1)^3 gauged supergravity, a known consistent truncation of string theory. We uncover a rich vacuum structure, including an infinite class of new timelike-warped AdS_3 (Goedel) and timelike-warped dS_3 critical points. We outline the construction of supersymmetric flows, driven by irrelevant scalar operators in the SCFT, which interpolate between critical points. For flows from AdS_3 to Goedel, the natural candidate for the central charge decreases along the flow. Flows to timelike-warped dS_3 exhibit topology change.


title: Dynamics of Dirac strings and monopole-like excitations in chiral magnets under a current drive

abstract: Skyrmion lines in metallic chiral magnets carry an emergent magnetic field experienced by the conduction electrons. The inflow and outflow of this field across a closed surface is not necessarily equal, thus it allows for the existence of emergent monopoles. One example is a segment of skyrmion line inside a crystal, where a monopole and antimonopole pair is connected by the emergent magnetic flux line. This is a realization of Dirac string-like excitations. Here we study the dynamics of monopoles in chiral magnets under an electric current. We show that in the process of creation of skyrmion lines, skyrmion line segments are first created via the proliferation of monopoles and antimonopoles. Then these line segments join and span the whole system through the annihilation of monopoles. The skyrmion lines are destroyed via the proliferation of monopoles and antimonopoles at high currents, resulting in a chiral liquid phase. We also propose to create the monopoles in a controlled way by applying an inhomogenous current to a crystal. Remarkably, an electric field component in the magnetic field direction proportional to the current squared in the low current region is induced by the motion of distorted skyrmion lines, in addition to the Hall and longitudinal voltage. The existence of monopoles can be inferred from transport or imaging measurements.


title: A Note On Nilpotent Representations

abstract: Let $\Gamma$ be a finitely generated nilpotent group and let G be a complex reductive algebraic group. The representation variety $\mathrm{Hom}(\Gamma,G)$ and the character variety $\mathrm{Hom}(\Gamma,G)//G$ each carry a natural topology, and we describe the topology of their connected components in terms of representations factoring through quotients of $\Gamma$ by elements of its lower central series.


title: Does Learning Imply a Decrease in the Entropy of Behavior?

abstract: Shannon's information entropy measures of the uncertainty of an event's outcome. If learning about a system reflects a decrease in uncertainty, then a plausible intuition is that learning should be accompanied by a decrease in the entropy of the organism's actions and/or perceptual states. To address whether this intuition is valid, I examined an artificial organism -- a simple robot -- that learned to navigate in an arena and analyzed the entropy of the outcome variables action, state, and reward. Entropy did indeed decrease in the initial stages of learning, but two factors complicated the scenario: (1) the introduction of new options discovered during the learning process and (2) the shifting patterns of perceptual and environmental states resulting from changes to the robot's learned movement strategies. These factors lead to a subsequent increase in entropy as the agent learned. I end with a discussion of the utility of information-based characterizations of learning.


title: Realization and Extension of Abstract Operation Contracts for Program Logic

abstract: For engineering software with formal correctness proofs it is crucial that proofs can be efficiently reused in case the software or its specification is changed. Unfortunately, in reality even slight changes in the code or its specification often result in disproportionate waste of verification effort: For instance, whenever a method's specification is modified and as a consequence the proof of its correctness breaks, all other proofs based on this specification break too. Abstract method calls is a recently proposed verification rule for method calls that allows for efficient systematic reuse of proofs. In this thesis, we implement, extend and evaluate this approach within the KeY verification system.


title: Non-constant volume exponential solutions in higher-dimensional Lovelock cosmologies

abstract: In this paper we propose a scheme which allows one to find all possible exponential solutions of special class -- non-constant volume solutions -- in Lovelock gravity in arbitrary number of dimensions and with arbitrate combinations of Lovelock terms. We apply this scheme to (6+1)- and (7+1)-dimensional flat anisotropic cosmologies in Einstein-Gauss-Bonnet and third-order Lovelock gravity to demonstrate how our scheme does work. In course of this demonstration we derive all possible solutions in (6+1) and (7+1) dimensions and compare solutions and their abundance between cases with different Lovelock terms present. As a special but more "physical" case we consider spaces which allow three-dimensional isotropic subspace for they could be viewed as examples of compactification schemes. Our results suggest that the same solution with three-dimensional isotropic subspace is more "probable" to occur in the model with most possible Lovelock terms taken into account, which could be used as kind of anthropic argument for consideration of Lovelock and other higher-order gravity models in multidimensional cosmologies.


title: Consumption-Investment Problem with Transaction Costs for L\'evy-Driven Price Processes

abstract: We consider an optimal control problem for a linear stochastic integro-diffe\-rential equation with conic constraints on the phase variable and the control of singular-regular type. Our setting includes consumption-investment problems for models of financial markets in the presence of proportional transaction costs where the price of the assets are given by a geometric L\'evy process and the investor is allowed to take short positions. We prove that the Bellman function of the problem is a viscosity solution of the HJB equation. A uniqueness theorem for the solution of the latter is established. Special attention is paid to the Dynamic Programming Principle.


title: Constrained BSDEs representation of the value function in optimal control of pure jump Markov processes

abstract: We consider a classical finite horizon optimal control problem for continuous-time pure jump Markov processes described by means of a rate transition measure depending on a control parameter and controlled by a feedback law. For this class of problems the value function can often be described as the unique solution to the corresponding Hamilton-Jacobi-Bellman equation. We prove a probabilistic representation for the value function, known as nonlinear Feynman-Kac formula. It relates the value function with a backward stochastic differential equation (BSDE) driven by a random measure and with a sign constraint on its martingale part. We also prove existence and uniqueness results for this class of constrained BSDEs. The connection of the control problem with the constrained BSDE uses a control randomization method recently developed in the works of I. Kharroubi and H. Pham and their co-authors. This approach also allows to prove that the value function of the original non-dominated control problem coincides with the value function of an auxiliary dominated control problem, expressed in terms of equivalent changes of probability measures.


title: No arbitrage and local martingale deflators

abstract: A supermartingale deflator (resp., local martingale deflator) multiplicatively transforms nonnegative wealth processes into supermartingales (resp., local martingales). The supermartingale numeraire (resp., local martingale numeraire) is the wealth processes whose reciprocal is a supermartingale deflator (resp., local martingale deflator). It has been established in previous literature that absence of arbitrage of the first kind (NA1) is equivalent to existence of the supermartingale numeraire, and further equivalent to existence of a strictly positive local martingale deflator; however, under NA1, the local martingale numeraire may fail to exist. In this work, we establish that, under NA1, any total-variation neighbourhood of the original probability has an equivalent probability under which the local martingale numeraire exists. This result, available previously only for single risky-asset models, is in striking resemblance with the fact that any total-variation neighbourhood of a separating measure contains an equivalent $\sigma$-martingale measure. The presentation of our main result is relatively self-contained, including a proof of existence of the supermartingale numeraire under NA1. We further show that, if the Levy measures of the asset-price process have finite support, NA1 is equivalent to existence of the local martingale numeraire with respect to the original probability.


title: Thirty Years of Turnstiles and Transport

abstract: To characterize transport in a deterministic dynamical system is to compute exit time distributions from regions or transition time distributions between regions in phase space. This paper surveys the considerable progress on this problem over the past thirty years. Primary measures of transport for volume-preserving maps include the exiting and incoming fluxes to a region. For area-preserving maps, transport is impeded by curves formed from invariant manifolds that form partial barriers, e.g., stable and unstable manifolds bounding a resonance zone or cantori, the remnants of destroyed invariant tori. When the map is exact volume preserving, a Lagrangian differential form can be used to reduce the computation of fluxes to finding a difference between the action of certain key orbits, such as homoclinic orbits to a saddle or to a cantorus. Given a partition of phase space into regions bounded by partial barriers, a Markov tree model of transport explains key observations, such as the algebraic decay of exit and recurrence distributions.


title: General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

abstract: We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.


title: Source Compression with a Quantum Helper

abstract: We study classical source coding with quantum side-information where the quantum side-information is observed by a helper and sent to the decoder via a classical channel. We derive a single-letter characterization of the achievable rate region for this problem. The direct part of our result is proved via the measurement compression theory by Winter. Our result reveals that a helper's scheme that separately conducts a measurement and a compression is suboptimal, and the measurement compression is fundamentally needed to achieve the optimal rate region.


title: Reconstruction-free action inference from compressive imagers

abstract: Persistent surveillance from camera networks, such as at parking lots, UAVs, etc., often results in large amounts of video data, resulting in significant challenges for inference in terms of storage, communication and computation. Compressive cameras have emerged as a potential solution to deal with the data deluge issues in such applications. However, inference tasks such as action recognition require high quality features which implies reconstructing the original video data. Much work in compressive sensing (CS) theory is geared towards solving the reconstruction problem, where state-of-the-art methods are computationally intensive and provide low-quality results at high compression rates. Thus, reconstruction-free methods for inference are much desired. In this paper, we propose reconstruction-free methods for action recognition from compressive cameras at high compression ratios of 100 and above. Recognizing actions directly from CS measurements requires features which are mostly nonlinear and thus not easily applicable. This leads us to search for such properties that are preserved in compressive measurements. To this end, we propose the use of spatio-temporal smashed filters, which are compressive domain versions of pixel-domain matched filters. We conduct experiments on publicly available databases and show that one can obtain recognition rates that are comparable to the oracle method in uncompressed setup, even for high compression ratios.


title: Synergy, suppression and immorality: forward differences of the entropy function

abstract: Conditional mutual information is important in the selection and interpretation of graphical models. Its empirical version is well known as a generalised likelihood ratio test and that it may be represented as a difference in entropy. We consider the forward difference expansion of the entropy function defined on all subsets of the variables under study. The elements of this expansion are invariant to permutation of their suffices and relate higher order mutual informations to lower order ones. The third order difference is expressible as an, apparently assymmetric, difference between a marginal and a conditional mutual information. Its role in the decomposition for explained information provides a technical definition for synergy between three random variables. Positive values occur when two variables provide alternative explanations for a third; negative values, termed synergies, occur when the sum of explained information is greater than the sum of its parts. Synergies tend to be infrequent; they connect the seemingly unrelated concepts of suppressor variables in regression, on the one hand, and unshielded colliders in Bayes networks (immoralities), on the other. We give novel characterizations of these phenomena that generalise to categorical variables and to higher dimensions. We propose an algorithm for systematically computing low order differences from a given graph. Examples from small scale real-life studies indicate the potential of these techniques for empirical statistical analysis.


title: One-dimensional Fermi gas with a single impurity in a harmonic trap: Perturbative description of the upper branch

abstract: The transition from "few to many" has recently been probed experimentally in an ultra cold harmonically confined one-dimensional lithium gas, in which a single impurity atom interacts with a background gas consisting of one, two, or more identical fermions [A. N. Wenz {\em{et al.}}, Science {\bf{342}}, 457 (2013)]. For repulsive interactions between the background or majority atoms and the impurity, the interaction energy for relatively moderate system sizes was analyzed and found to converge toward the corresponding expression for an infinitely large Fermi gas. Motivated by these experimental results, we apply perturbative techniques to determine the interaction energy for weak and strong coupling strengths and derive approximate descriptions for the interaction energy for repulsive interactions with varying strength between the impurity and the majority atoms and any number of majority atoms.


title: Structure Learning in Bayesian Networks of Moderate Size by Efficient Sampling

abstract: We study the Bayesian model averaging approach to learning Bayesian network structures (DAGs) from data. We develop new algorithms including the first algorithm that is able to efficiently sample DAGs according to the exact structure posterior. The DAG samples can then be used to construct estimators for the posterior of any feature. We theoretically prove good properties of our estimators and empirically show that our estimators considerably outperform the estimators from the previous state-of-the-art methods.


title: Universal formula for the energy--momentum tensor via a flow equation in the Gross--Neveu model

abstract: For the fermion field in the two-dimensional Gross--Neveu model, we introduce a flow equation that allows a simple $1/N$ expansion. By employing the $1/N$ expansion, we examine the validity of a universal formula for the energy--momentum tensor which is based on the small flow-time expansion. We confirm that the formula reproduces a correct normalization and the conservation law of the energy--momentum tensor by computing the translation Ward--Takahashi relation in the leading non-trivial order in the $1/N$ expansion. Also, we confirm that the expectation value at finite temperature correctly reproduces thermodynamic quantities. These observations support the validity of a similar construction of the energy--momentum tensor via the gradient/Wilson flow in lattice gauge theory.


title: What can we learn from higher multipole power spectra of galaxy distribution in redshift space?

abstract: We investigate a potential of the higher multipole power spectra of the galaxy distribution in redshift space as a cosmological probe on halo scales. Based on the fact that a halo model explains well the multipole power spectra of the luminous red galaxy (LRG) sample in the Sloan Digital Sky Survey (SDSS), we focus our investigation on the random motions of the satellite LRGs that determine the higher multipole spectra at large wavenumbers. We show that our theoretical model fits the higher multipole spectra at large wave numbers from N-body numerical simulations and we apply these results for testing the gravity theory and the velocity structure of galaxies on the halo scales. In this analysis, we use the multipole spectra P_4(k) and P_6(k) on the small scales of the range of wavenumber 0.3<k/[h{Mpc}^{-1}]<0.6, which is in contrast to the usual method of testing gravity by targeting the linear growth rate on very large scales. We demonstrate that our method could be useful for testing gravity on the halo scales.


title: A topological semigroup structure on the space of actions modulo weak equivalence

abstract: We introduce a topology on the space of actions modulo weak equivalence finer than the one previously studied in the literature. We show that the product of actions is a continuous operation with respect to this topology, so that the space of actions modulo weak equivalence becomes a topological semigroup.


title: Spatial Auditory Brain-computer Interface using Head Related Impulse Response

abstract: This study reports on a head related impulse response (HRIR) application to an auditory spatial brain-computer interface (BCI) speller paradigm. Six experienced and five BCI-naive users participated in an experimental spelling set up based on five Japanese vowels. Obtained auditory evoked potentials resulted with encouragingly good and stable P300-responses in online BCI experiments. Our case study indicated that the headphone reproduced auditory (HRIR- based) spatial sound paradigm could be a viable alternative to the established multi-loudspeaker surround sound BCI-speller applications, as far as healthy pilot study users are concerned.


title: On Normalizers of $C^*$-Subalgebras in the Cuntz Algebra ${\mathcal O}_n$. II

abstract: We investigate subalgebras $A$ of the Cuntz algebra ${\mathcal O}_n$ that arise as finite direct sums of corners of the UHF-subalgebra ${\mathcal F}_n$. For such an $A$, we completely determine its normalizer group inside ${\mathcal O}_n$.


title: Optimal Power Allocation for Secure Communications in Large-Scale MIMO Relaying Systems

abstract: In this paper, we address the problem of optimal power allocation at the relay in two-hop secure communications. In order to solve the challenging issue of short-distance interception in secure communications, the benefit of large-scale MIMO (LS-MIMO) relaying techniques is exploited to improve the secrecy performance significantly, even in the case without eavesdropper channel state information (CSI). The focus of this paper is on the analysis and design of optimal power allocation for the relay, so as to maximize the secrecy outage capacity. We reveal the condition that the secrecy outage capacity is positive, prove that there is one and only one optimal power, and present an optimal power allocation scheme. Moreover, the asymptotic characteristics of the secrecy outage capacity is carried out to provide some clear insights for secrecy performance optimization. Finally, simulation results validate the effectiveness of the proposed scheme.


title: A gate defined quantum dot on the two-dimensional transition metal dichalcogenide semiconductor WSe2

abstract: Two-dimensional layered materials, such as transition metal dichalcogenides (TMDCs), are promising materials for future electronics owing to their unique electronic properties. With the presence of a band gap, atomically thin gate defined quantum dots (QDs) can be achieved on TMDCs. Here, standard semiconductor fabrication techniques are used to demonstrate quantum confined structures on WSe2 with tunnel barriers defined by electric fields, thereby eliminating the edge states induced by etching steps, which commonly appear in gapless graphene QDs. Over 40 consecutive Coulomb diamonds with a charging energy of approximately 2 meV were observed, showing the formation of a QD, which is consistent with the simulations. The size of the QD could be tuned over a factor of 2 by changing the voltages applied to the top gates. These results shed light on quantum nano-devices on TMDCs for further researches.


title: Instance Significance Guided Multiple Instance Boosting for Robust Visual Tracking

abstract: Multiple Instance Learning (MIL) recently provides an appealing way to alleviate the drifting problem in visual tracking. Following the tracking-by-detection framework, an online MILBoost approach is developed that sequentially chooses weak classifiers by maximizing the bag likelihood. In this paper, we extend this idea towards incorporating the instance significance estimation into the online MILBoost framework. First, instead of treating all instances equally, with each instance we associate a significance-coefficient that represents its contribution to the bag likelihood. The coefficients are estimated by a simple Bayesian formula that jointly considers the predictions from several standard MILBoost classifiers. Next, we follow the online boosting framework, and propose a new criterion for the selection of weak classifiers. Experiments with challenging public datasets show that the proposed method outperforms both existing MIL based and boosting based trackers.


title: Twist-3 fragmentation functions in a spectator model with gluon rescattering

abstract: We study the twist-3 fragmentation functions $H$ and $\tilde{H}$, by applying a spectator model. In the calculation we consider the effect of the gluon rescattering at one loop level. We find that in this case the hard-vertex diagram, which gives zero contribution to the Collins function, does contribute to the fragmentation function $H$. The calculation shows that the twist-3 T-odd fragmentation functions are free of light-cone divergences. The parameters of the model are fitted from the known parametrization of the unpolarized fragmentation $D_1$ and the Collins function $H_1^\perp$. We find our result for the favored fragmentation function is consistent with the recent extraction on $H$ and $\tilde{H}$ from pp data. We also check numerically the equation of motion relation for $H$, $\tilde{H}$ and find that relation holds fairly well in the spectator model.


title: The ranks of homotopy groups of Kac-Moody groups

abstract: Let $A$ be a Cartan matrix and $G(A)$ be the Kac-Moody group associated to Cartan matrix $A$. In this paper, we discuss the computation of the rank $i_k$ of homotopy group $\pi_k(G(A))$. For a large class of Kac-Moody groups, we construct Lie algebras with grade from the Poincar\'{e} series of their flag manifolds. And we interpret $i_{2k}$ as the dimension of the degree $2k$ homogeneous component of the Lie algebra we constructed. Since the computation of $i_{2k-1}$ is trivial, this gives a combinatorics interpretation of $i_k$ for all $k>0$.


title: Low scale quantum gravity in gauge-Higgs unified models

abstract: We consider the scale at which gravity becomes strong in linearized General Relativity coupled to the gauge-Higgs unified(GHU) model. We also discuss the unitarity of S-matrix in the same framework. The Kaluza-Klein(KK) gauge bosons, KK scalars and KK fermions in the GHU models can drastically change the strong gravity scale and the unitarity violation scale. In particular we consider two models GHU_SM and GHU_MSSM which have the zero modes corresponding to the particle content of the Standard Model and the Minimal Supersymmetric Standard Model, respectively. We find that the strong gravity scale could be lowered as much as 10^13 (10^14) GeV in the GHU_SM (GHU_MSSM) for one extra dimension taking 1 TeV as the compactification scale. It is also shown that these scales are proportional to the inverse of the number of extra dimensions d. In the d=10 case, they could be lowered up to 10^5 GeV for both models. We also find that the maximum compactification scales of extra dimensions quickly converge into one special scale M_O near Planck scale or equivalently into one common radius R_0 irrespectively of d as the number of zero modes increases. It may mean that all extra dimensions emerge with the same radius near Planck scale. In addition, it is shown that the supersymmetry can help to remove the discordance between the strong gravity scale and the unitarity violation scale.


title: The accretion history of dark matter halos II: The connections with the mass power spectrum and the density profile

abstract: We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^{\alpha}e^{\beta z}$, where the parameters $\alpha$ and $\beta$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $\alpha$, $\beta$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.


title: The 'Power of Then': The Uniquely Human Capacity to Imagine Beyond the Present

abstract: Humans appear to be uniquely in their ability to transcend the present and reflect on current ideas in terms of what was experienced in the past, or fantasize about the future. This paper presents an account, in layperson terms, of how this ability came about, its importance in modern life, and why it defines our human-ness.


title: On triangle-free graphs of order 10 with prescribed 1-defective chromatic number

abstract: A graph is (m, k)-colourable if its vertices can be coloured with m colours such that the maximum degree of any subgraph induced on ver- tices receiving the same colour is at most k. The k-defective chromatic number for a graph is the least positive integer m for which the graph is (m, k)-colourable. All triangle-free graphs on 8 or fewer vertices are (2, 1)-colourable. There are exactly four triangle-free graphs of order 9 which have 1-defective chromatic number 3. We show that these four graphs appear as subgraphs in almost all triangle-free graphs of order 10 with 1-defective chromatic number equal to 3. In fact there is a unique triangle-free (3, 1)-critical graph on 10 vertices and we exhibit this graph.


title: Investigation of Dark Matter in Minimal 3-3-1 Models

abstract: It is shown that the 3-3-1 model with the minimal lepton content can work as two-Higgs-triplet 3-3-1 model while leaving the other scalars as inert particles responsible for dark matter. We study two cases of dark matter corresponding to the doublet and singlet scalar candidates. We figure out the parameter spaces in the WMAP allowed region of the relic density. The indirect and direct searches for dark matter in both cases are investigated by using micrOMEGAs.


title: Scaling of nuclear modification factors for hadrons and light nuclei

abstract: The number of constituent quarks (NCQ-) scaling of hadrons and the number of constituent nucleons (NCN-) scaling of light nuclei are proposed for nuclear modification factors ($R_{cp}$) of hadrons and light nuclei, respectively, according to the experimental investigations in relativistic heavy-ion collisions. Based on coalescence mechanism the scalings are performed for pions and protons in quark level, and light nuclei $d (\bar d)$ and $^3$He for nucleonic level, respectively, formed in Au + Au and Pb + Pb collisions and nice scaling behaviour emerges. NCQ or NCN scaling law of $R_{cp}$ can be respectively taken as a probe for quark or nucleon coalescence mechanism for the formation of hadron or light nuclei in relativistic heavy-ion collisions.


title: On minimal triangle-free planar graphs with prescribed 1-defective chromatic number

abstract: A graph is (m,k)-colourable if its vertices can be coloured with m colours such that the maximum degree of the subgraph induced on the set of all vertices receiving the same colour is at most k. The k-defective chromatic number $\chi_k(G)$ is the least positive integer $m$ for which graph G is (m,k)-colourable. Let f(m,k;tfp) be the smallest order of a triangle-free planar graph such that $\chi_k(G)$=m. In this paper we show that f(3,1;tfp)=11.


title: Chromatic and flow polynomials of generalized vertex join graphs and outerplanar graphs

abstract: A generalized vertex join of a graph is obtained by joining an arbitrary multiset of its vertices to a new vertex. We present a low-order polynomial time algorithm for finding the chromatic polynomials of generalized vertex joins of trees, and by duality we find the flow polynomials of arbitrary outerplanar graphs. We also present closed formulas for the chromatic and flow polynomials of vertex joins of cliques and cycles, otherwise known as "generalized wheel" graphs.


title: Revisiting Entanglement Entropy of Lattice Gauge Theories

abstract: Casini et al raise the issue that the entanglement entropy in gauge theories is ambiguous because its definition depends on the choice of the boundary between two regions.; even a small change in the boundary could annihilate the otherwise finite topological entanglement entropy between two regions. In this article, we first show that the topological entanglement entropy in the Kitaev model which is not a true gauge theory, is free of ambiguity. Then, we give a physical interpretation, from the perspectives of what can be measured in an experiement, to the purported ambiguity of true gauge theories, where the topological entanglement arises as redundancy in counting the degrees of freedom along the boundary separating two regions. We generalize these discussions to non-Abelian gauge theories.


title: Search for {\eta}'(958)-nucleus bound states by (p,d) reaction at GSI and FAIR

abstract: The mass of the {\eta}' meson is theoretically expected to be reduced at finite density, which indicates the existence of {\eta}'-nucleus bound states. To investigate these states, we perform missing-mass spectroscopy for the (p, d) reaction near the {\eta}' production threshold. The overview of the experimental situation is given and the current status is discussed.


title: How a Generation Was Misled About Natural Selection

abstract: This article explains how natural selection works and how it has been inappropriately applied to the description of cultural change. It proposes an alternative evolutionary explanation for cultural evolution that describes it in terms of communal exchange.


title: Isolation in the construction of natural experiments

abstract: A natural experiment is a type of observational study in which treatment assignment, though not randomized by the investigator, is plausibly close to random. A process that assigns treatments in a highly nonrandom, inequitable manner may, in rare and brief moments, assign aspects of treatments at random or nearly so. Isolating those moments and aspects may extract a natural experiment from a setting in which treatment assignment is otherwise quite biased, far from random. Isolation is a tool that focuses on those rare, brief instances, extracting a small natural experiment from otherwise useless data. We discuss the theory behind isolation and illustrate its use in a reanalysis of a well-known study of the effects of fertility on workforce participation. Whether a woman becomes pregnant at a certain moment in her life and whether she brings that pregnancy to term may reflect her aspirations for family, education and career, the degree of control she exerts over her fertility, and the quality of her relationship with the father; moreover, these aspirations and relationships are unlikely to be recorded with precision in surveys and censuses, and they may confound studies of workforce participation. However, given that a women is pregnant and will bring the pregnancy to term, whether she will have twins or a single child is, to a large extent, simply luck. Given that a woman is pregnant at a certain moment, the differential comparison of two types of pregnancies on workforce participation, twins or a single child, may be close to randomized, not biased by unmeasured aspirations. In this comparison, we find in our case study that mothers of twins had more children but only slightly reduced workforce participation, approximately 5% less time at work for an additional child.


title: The multi-frequency parsec-scale structure of PKS 2254-367 (IC 1459): a luminosity-dependent break in morphology for the precursors of radio galaxies?

abstract: We present the first multi-frequency VLBI images of PKS 2254-367, a Giga-hertz Peaked Spectrum (GPS) radio source hosted by the nearby galaxy IC 1459 (D=20.5 Mpc). PKS 2254-367 and the radio source in NGC 1052 (PKS 0238-084; D=17.2 Mpc) are the two closest GPS radio sources to us, far closer than the next closest example, PKS 1718-649 (D=59 Mpc). As such, IC 1459 and NGC 1052 offer opportunities to study the details of the pc-scale radio sources as well as the environments that the radio sources inhabit, across the electromagnetic spectrum. Given that some models for the origin and evolution of GPS radio sources require a strong connection between the radio source morphology and the gaseous nuclear environment, such opportunities for detailed study are important. Our VLBI images of PKS 2254-367 show that the previously identified similarities between IC 1459 and NGC 1052 continue onto the pc-scale. Both compact radio sources appear to have symmetric jets of approximately the same luminosity, much lower than typically noted in compact double GPS sources. Similarities between PKS 2254-367 and NGC 1052, and differences with respect to other GPS galaxies, lead us to speculate that a sub-class of GPS radio sources, with low luminosity and with jet-dominated morphologies, exists and would be largely absent from radio source surveys with ~1 Jy flux density cutoffs. We suggest that this possible low-luminosity, jet-dominated population of GPS sources could be an analog of the FR-I radio galaxies, with the higher luminosity lobe-dominated GPS sources being the analog of the FR-II radio galaxies.


title: Band Splitting Permutations for Spatially Coupled LDPC Codes Enhancing Burst Erasure Immunity

abstract: It is well known that spatially coupled (SC) codes with erasure-BP decoding have powerful error correcting capability over memoryless erasure channels. However, the decoding performance of SC-codes significantly degrades when they are used over burst erasure channels. In this paper, we propose band splitting permutations (BSP) suitable for $(l,r,L)$ SC-codes. The BSP splits a diagonal band in a base matrix into multiple bands in order to enhance the span of the stopping sets in the base matrix. As theoretical performance guarantees, lower and upper bounds on the maximal burst correctable length of the permuted $(l,r,L)$ SC-codes are presented. Those bounds indicate that the maximal correctable burst ratio of the permuted SC-codes converges to 1/k where k=r/l. This implies the asymptotic optimality of the permuted SC-codes in terms of burst erasure correction.


title: Spherical Harmonic Expansion of Fisher-Bingham Distribution and 3D Spatial Fading Correlation for Multiple-Antenna Systems

abstract: This paper considers the 3D spatial fading correlation (SFC) resulting from an angle-of-arrival (AoA) distribution that can be modelled by a mixture of Fisher-Bingham distributions on the sphere. By deriving a closed-form expression for the spherical harmonic transform for the component Fisher-Bingham distributions, with arbitrary parameter values, we obtain a closed-form expression of the 3D-SFC for the mixture case. The 3D-SFC expression is general and can be used in arbitrary multi-antenna array geometries and is demonstrated for the cases of a 2D uniform circular array in the horizontal plane and a 3D regular dodecahedral array. In computational aspects, we use recursions to compute the spherical harmonic coefficients and give pragmatic guidelines on the truncation size in the series representations to yield machine precision accuracy results. The results are further corroborated through numerical experiments to demonstrate that the closed-form expressions yield the same results as significantly more computationally expensive numerical integration methods.


title: Perfect state transfer in products and covers of graphs

abstract: A continuous-time quantum walk on a graph $X$ is represented by the complex matrix $\exp (-\mathrm{i} t A)$, where $A$ is the adjacency matrix of $X$ and $t$ is a non-negative time. If the graph models a network of interacting qubits, transfer of state among such qubits throughout time can be formalized as the action of the continuous-time quantum walk operator in the characteristic vectors of the vertices. Here we are concerned with the problem of determining which graphs admit a perfect transfer of state. More specifically, we will study graphs whose adjacency matrix is a sum of tensor products of $01$-matrices, focusing on the case where a graph is the tensor product of two other graphs. As a result, we will construct many new examples of perfect state transfer.


title: Nuclear Equation of State and Neutron Star Cooling

abstract: We investigate the cooling of neutron stars with relativistic and non-relativistic models of dense nuclear matter. We focus on the effects of uncertainties originated from the nuclear models, the composition of elements in the envelope region, and the formation of superfluidity in the core and the crust of neutron stars. Discovery of $2 M_\odot$ neutron stars PSR J1614-2230 and PSR J0343+0432 has triggered the revival of stiff nuclear equation of state at high densities. In the mean time, observation of a neutron star in Cassiopeia A for more than 10 years has provided us with very accurate data for the thermal evolution of neutron stars. Both mass and temperature of neutron stars depend critically on the equation of state of nuclear matter, so we first search for nuclear models that satisfy the constraints from mass and temperature simultaneously within a reasonable range. With selected models, we explore the effects of element composition in the evenlope region, and the existence of superfluidity in the core and the crust of neutron stars. Due to uncertainty in the composition of particles in the envelope region we obtain a range of cooling curves that can cover substantial region of observation data.


title: Spectrally extremal vertices, strong cospectrality and state transfer

abstract: In order to obtain perfect state transfer between two sites in a network of interacting qubits, their corresponding vertices in the underlying graph must satisfy a combinatorial property called strong cospectrality. Here we determine the structure of graphs containing pairs of vertices which are strongly cospectral and satisfy a certain extremal property related to the spectrum of the graph. If the graph satisfies this property globally and is regular, we also show that the existence of a partition of the vertex set into pairs of vertices at maximum distance admitting perfect state transfer forces the graph to be distance-regular.


abstract: We consider Diophantine quintuples $\{a, b, c, d, e\}$, sets of distinct positive integers the product of any two elements of which is one less than a perfect square. Triples of the first kind are the subsets $\{a, b, d\}$ with $d> b^{5}$. We show that there are no Diophantine quintuples containing triples of the first kind.


title: A counterexample shows that not every locally $L^0$--convex topology is necessarily induced by a family of $L^0$--seminorms

abstract: This paper constructs a counterexample showing that not every locally $L^0$--convex topology is necessarily induced by a family of $L^0$--seminorms. Random convex analysis is the analytic foundation for $L^0$--convex conditional risk measures, this counterexample, however, shows that a locally $L^0$--convex module is not a proper framework for random convex analysis. Further, this paper also gives a necessary and sufficient condition for a locally $L^0$--convex topology to be induced by a family of $L^0$--seminorms. Finally, we give some comments showing that based on random locally convex modules, we can establish a perfect random convex analysis to meet the needs of the study of $L^0$--convex conditional risk measures.


title: Bounds on the number of Diophantine quintuples

abstract: We consider Diophantine quintuples $\{a, b, c, d, e\}$. These are sets of distinct positive integers, the product of any two elements of which is one less than a perfect square. It is conjectured that there are no Diophantine quintuples; we improve on current estimates to show that there are at most $1.9\cdot 10^{29}$ Diophantine quintuples.


title: Subgraph Domatic Problem and Writing Capacity of Memory Devises with Restricted State Transitions

abstract: A code design problem for memory devises with restricted state transitions is formulated as a combinatorial optimization problem that is called a subgraph domatic partition (subDP) problem. If any neighbor set of a given state transition graph contains all the colors, then the coloring is said to be valid. The goal of a subDP problem is to find a valid coloring with the largest number of colors for a subgraph of a given directed graph. The number of colors in an optimal valid coloring gives the writing capacity of a given state transition graph. The subDP problems are computationally hard; it is proved to be NP-complete in this paper. One of our main contributions in this paper is to show the asymptotic behavior of the writing capacity $C(G)$ for sequences of dense bidirectional graphs, that is given by C(G)=Omega(n/ln n) where n is the number of nodes. A probabilistic method called Lovasz local lemma (LLL) plays an essential role to derive the asymptotic expression.


title: Surface Wave Propagation in non--ideal plasmas

abstract: The properties of surface waves in a partially ionized, compressible magnetized plasma slab are investigated in this work. The waves are affected by the nonideal magnetohydrodynamic effects which causes finite drift of the magnetic field in the medium. When the magnetic field drift is ignored, the characteristics of the wave propagation in a partially ionized plasma fluid is similar to the fully ionized ideal MHD except now the propagation properties depend on the fractional ionization as well as on the compressibility of the medium. The phase velocity of the sausage and kink waves increases marginally (by a few percent) due to the compressibility of the medium in both ideal as well as Hall diffusion dominated regimes. However, unlike ideal regime, only waves below certain cut off frequency can propagate in the medium in Hall dominated regime. This cut off for a thin slab has a weak dependence on the plasma beta whereas for thick slab no such dependence exists. More importantly, since the cut off is introduced by the Hall diffusion, the fractional ionization of the medium is more important than the plasma compressibility in determining such a cut off. We discuss the relevance of these results in the context of solar photosphere-chromosphere.


title: Bayesian Ages for Early-Type Stars from Isochrones Including Rotation, and a Possible Old Age for the Hyades

abstract: We combine recently computed models of stellar evolution using a new treatment of rotation with a Bayesian statistical framework to constrain the ages and other properties of early-type stars. We find good agreement for early-type stars and clusters with known young ages, including beta Pictoris, the Pleiades, and the Ursa Majoris Moving Group. However, we derive a substantially older age for the Hyades open cluster (750+/-100 Myr compared to 625+/-50 Myr). This older age results from both the increase in main-sequence lifetime with stellar rotation and from the fact that rotating models near the main-sequence turnoff are more luminous, overlapping with slightly more massive (and shorter-lived) nonrotating ones. Our method uses a large grid of nonrotating models to interpolate between a much sparser rotating grid, and also includes a detailed calculation of synthetic magnitudes as a function of orientation. We provide a web interface at www.bayesianstellarparameters.info where the results of our analysis may be downloaded for individual early-type (B-V<~0.25) Hipparcos stars. The web interface accepts user-supplied parameters for a Gaussian metallicity prior and returns posterior probability distributions on mass, age, and orientation.


title: Suzaku X-ray Observations of the Fermi Bubbles: Northernmost Cap and Southeast Claw Discovered with MAXI-SSC

abstract: We report on Suzaku observations of large-scale X-ray structures possibly related with the Fermi Bubbles obtained in 2013 with a total duration of ~ 80 ks. The observed regions were the: (i) northern cap (N-cap; l ~ 0 deg, 45 deg < b < 55 deg) seen in the Mid-band (1.7-4.0 keV) map recently provided by MAXI-SSC and (ii) southeast claw (SE-claw; l ~ 10 deg, -20 deg < b < -10 deg) seen in the ROSAT all-sky map and MAXI-SSC Low-band (0.7-1.7 keV) map. In each region, we detected diffuse X-ray emissions which are represented by a three component plasma model consisting of an unabsorbed thermal component (kT ~ 0.1 keV) from the Local Bubble, absorbed kT = 0.30+/-0.05 keV emission representing the Galactic Halo, and a power-law component due to the isotropic cosmic X-ray background radiation. The emission measure of the GH component in the SE-claw shows an excess by a factor of ~ 2.5 over the surrounding emission at 2 deg away. We also found a broad excess in the 1.7-4.0 keV count rates across the N-cap after compiling other archival data from Suzaku and Swift. The spectral stacking analysis of the N-cap data indicates the presence of another thermal component with kT = 0.70 (+0.22,-0.11) keV. The temperature of kT ~ 0.3 keV of the Galactic Halo is higher than the ubiquitous value of kT ~ 0.2 keV near the Fermi Bubbles, and can be even higher (~ 0.7 keV). We discuss our findings in the context of bubble-halo interaction.


title: How Creative Ideas Take Shape

abstract: According to the honing theory of creativity, creative thought works not on individually considered, discrete, predefined representations but on a contextually-elicited amalgam of items which exist in a state of potentiality and may not be readily separable. This leads to the prediction that analogy making proceeds not by mapping correspondences from candidate sources to target, as predicted by the structure mapping theory of analogy, but by weeding out non-correspondences, thereby whittling away at potentiality. Participants were given an analogy problem, interrupted before they had time to solve it, and asked to write down what they had by way of a solution. Na\"ive judges categorized responses as significantly more supportive of the predictions of honing theory than those of structure mapping.


title: Lead Iodide Perovskite Light-Emitting Field-Effect Transistor

abstract: Despite the widespread use of solution-processable hybrid organic-inorganic perovskites in photovoltaic and light-emitting applications, determination of their intrinsic charge transport parameters has been elusive due to the variability of film preparation and history-dependent device performance. Here we show that screening effects associated to ionic transport can be effectively eliminated by lowering the operating temperature of methylammonium lead iodide perovskite (CH3NH3PbI3) field-effect transistors (FETs). Field-effect carrier mobility is found to increase by almost two orders of magnitude below 200 K, consistent with phonon scattering limited transport. Under balanced ambipolar carrier injection, gate-dependent electroluminescence is also observed from the transistor channel, with spectra revealing the tetragonal to orthorhombic phase transition. This first demonstration of CH3NH3PbI3 light-emitting FETs provides intrinsic transport parameters to guide materials and solar cell optimization, and will drive the development of new electro-optic device concepts, such as gated light emitting diodes and lasers operating at room temperature.


title: A patient-specific scatter artifacts correction method

abstract: This paper provides a fast and patient-specific scatter artifact correction method for cone-beam computed tomography (CBCT) used in image-guided interventional procedures. Due to increased irradiated volume of interest in CBCT imaging, scatter radiation has increased dramatically compared to 2D imaging, leading to a degradation of image quality. In this study, we propose a scatter artifact correction strategy using an analytical convolution-based model whose free parameters are estimated using a rough estimation of scatter profiles from the acquired cone-beam projections. It was evaluated using Monte Carlo simulations with both monochromatic and polychromatic X-ray sources. The results demonstrated that the proposed method significantly reduced the scatter-induced shading artifacts and recovered CT numbers.


title: Fast Scatter Artifacts Correction for Cone-Beam CT without System Modification and Repeat Scan

abstract: We provide a fast and accurate scatter artifacts correction algorithm for cone beam CT (CBCT) imaging. The method starts with an estimation of coarse scatter profile for a set of CBCT images. A total-variation denoising algorithm designed specifically for Poisson signal is then applied to derive the final scatter distribution. Qualitatively and quantitatively evaluations using Monte Carlo (MC) simulations, experimental CBCT phantom data, and \emph{in vivo} human data acquired for a clinical image guided radiation therapy were performed. Results show that the proposed algorithm can significantly reduce scatter artifacts and recover the correct HU within either projection domain or image domain. Further test shows the method is robust with respect to segmentation procedure.


title: Special birational transformations of type (2,1)

abstract: A birational transformation f: P^n --> Z, where Z is a nonsingular variety of Picard number 1, is called a special birational transformation of type (a, b) if f is given by a linear system of degree a, its inverse is given by a linear system of degree b and the base locus S \subset P^n of f is irreducible and nonsingular. In this paper, we classify special birational transformations of type (2,1). In addition to previous works Alzati-Sierra and Russo on this topic, our proof employs natural C^*-actions on Z in a crucial way. These C^*-actions also relate our result to the problem studied in our previous work on smooth projective varieties with nonzero prolongations.


title: Active coupling control in densely packed subwavelength waveguides via dark mode

abstract: The ever growing need for energy-efficient and fast communications is driving the development of highly integrated photonic circuits where controlling light at the nanoscale becomes the most critical aspect of information transfer . Directional couplers, two interacting optical waveguides placed in close proximity, are important building blocks in these integrated photonics circuits and have been employed as optical modulators and switches for high speed communication, data processing and integrated quantum operations. However, active control over the coupling between closely packed waveguides is highly desirable and yet remains a critical barrier towards ultra small footprint devices. A general approach to achieve active control in waveguide systems is to exploit optical nonlinearities enabled by a strong control pulse. However these devices suffer from the nonlinear absorption induced by the intense control pulse as the signal and its control propagate in the same waveguide. Here we experimentally demonstrate a unique scheme based on adiabatic elimination (AE) concept that effectively manipulates the coupling between densely packed waveguides. We demonstrate active coupling control between two closely packed waveguides by tuning the mode index of an in-between decoupled waveguide. This is achieved via a dark mode and thus leaves the signal unaffected by the induced losses. Such a scheme is a promising candidate for ultra-dense integrated nano-photonics such as on-chip ultrafast modulators and tunable filters for optical communication and quantum computing.


title: Power Allocation Games on Interference Channels with Complete and Partial Information

abstract: We consider a wireless channel shared by multiple transmitter-receiver pairs. Their transmissions interfere with each other. Each transmitter-receiver pair aims to maximize its long-term average transmission rate subject to an average power constraint. This scenario is modeled as a stochastic game under different assumptions. We first assume that each transmitter and receiver has knowledge of all direct and cross link channel gains. We later relax the assumption to the knowledge of incident channel gains and then further relax to the knowledge of the direct link channel gains only. In all the cases, we formulate the problem of finding the Nash equilibrium as a variational inequality (VI) problem and present an algorithm to solve the VI.


title: Statistical-mechanical analysis of pre-training and fine tuning in deep learning

abstract: In this paper, we present a statistical-mechanical analysis of deep learning. We elucidate some of the essential components of deep learning---pre-training by unsupervised learning and fine tuning by supervised learning. We formulate the extraction of features from the training data as a margin criterion in a high-dimensional feature-vector space. The self-organized classifier is then supplied with small amounts of labelled data, as in deep learning. Although we employ a simple single-layer perceptron model, rather than directly analyzing a multi-layer neural network, we find a nontrivial phase transition that is dependent on the number of unlabelled data in the generalization error of the resultant classifier. In this sense, we evaluate the efficacy of the unsupervised learning component of deep learning. The analysis is performed by the replica method, which is a sophisticated tool in statistical mechanics. We validate our result in the manner of deep learning, using a simple iterative algorithm to learn the weight vector on the basis of belief propagation.

categories: stat.ML cond-mat.dis-nn cond-mat.stat-mech cs.AI cs.LG


title: Hierarchical Electrochemical Modeling and Simulation of Bio-Hybrid Interfaces

abstract: In this article we propose and investigate a hierarchy of mathematical models based on partial differential equations (PDE) and ordinary differential equations (ODE) for the simulation of the biophysical phenomena occurring in the electrolyte fluid that connects a biological component (a single cell or a system of cells) and a solid-state device (a single silicon transistor or an array of transistors). The three members of the hierarchy, ordered by decreasing complexity, are: (i) a 3D Poisson-Nernst-Planck (PNP) PDE system for ion concentrations and electric potential; (ii) a 2D reduced PNP system for the same dependent variables as in (i); (iii) a 2D area-contact PDE system for electric potential coupled with a system of ODEs for ion concentrations. The backward Euler method is adopted for temporal semi-discretization and a fixed-point iteration based on Gummel's map is used to decouple system equations. Spatial discretization is performed using piecewise linear triangular finite elements stabilized via edge-based exponential fitting. Extensively conducted simulation results are in excellent agreement with existing analytical solutions of the PNP problem in radial coordinates and experimental and simulated data using simplified lumped parameter models.


title: Imputation of truncated p-values for meta-analysis methods and its genomic application

abstract: Microarray analysis to monitor expression activities in thousands of genes simultaneously has become routine in biomedical research during the past decade. A tremendous amount of expression profiles are generated and stored in the public domain and information integration by meta-analysis to detect differentially expressed (DE) genes has become popular to obtain increased statistical power and validated findings. Methods that aggregate transformed $p$-value evidence have been widely used in genomic settings, among which Fisher's and Stouffer's methods are the most popular ones. In practice, raw data and $p$-values of DE evidence are often not available in genomic studies that are to be combined. Instead, only the detected DE gene lists under a certain $p$-value threshold (e.g., DE genes with $p$-value${}<0.001$) are reported in journal publications. The truncated $p$-value information makes the aforementioned meta-analysis methods inapplicable and researchers are forced to apply a less efficient vote counting method or na\"{i}vely drop the studies with incomplete information. The purpose of this paper is to develop effective meta-analysis methods for such situations with partially censored $p$-values. We developed and compared three imputation methods - mean imputation, single random imputation and multiple imputation - for a general class of evidence aggregation methods of which Fisher's and Stouffer's methods are special examples. The null distribution of each method was analytically derived and subsequent inference and genomic analysis frameworks were established. Simulations were performed to investigate the type I error, power and the control of false discovery rate (FDR) for (correlated) gene expression data. The proposed methods were applied to several genomic applications in colorectal cancer, pain and liquid association analysis of major depressive disorder (MDD). The results showed that imputation methods outperformed existing na\"{i}ve approaches. Mean imputation and multiple imputation methods performed the best and are recommended for future applications.


title: Messier 35 (NGC2168) DANCe I. Membership, proper motions and multi-wavelength photometry

abstract: We aim at identifying the cluster's members by deriving membership probabilities for the sources within 1 degree of the cluster's center, going further away than equivalent previous studies. We measure accurate proper motions and multi-wavelength (optical and near-infrared) photometry using ground based archival images of the cluster. We use these measurements to compute membership probabilities. The list of candidate members from Barrado+2001 is used as training set to identify the cluster's locus in a multi-dimensional space made of proper motions, luminosities and colors. The final catalog includes 338892 sources with multi-wavelength photometry. Approximately half (194452) were detected at more than two epochs and we measured their proper motion and used it to derive membership probability. A total of 4349 candidate members with membership probabilities greater than 50% are found in this sample in the luminosity range between 10 and 22mag. The slow proper motion of the cluster and the overlap of its sequence with the field and background sequences in almost all color-magnitude and color-color diagrams complicate the analysis and the contamination level is expected to be significant. Our study nevertheless provides a coherent and quantitative membership analysis of Messier 35 based on a large fraction of the best ground-based data sets obtained over the past 18 years. As such, it represents a valuable input for follow-up studies using in particular the Kepler K2 photometric time series.


title: Continuous Multi-line Queues and TASEP

abstract: In this paper, we study a distribution of labeled particles on a continuous ring. It arises in three different ways, all related to the multi-type TASEP on a ring. We prove formulas for the probability density function for some permutations and give conjectures for a larger class. We give a complete conjecture for the probability of two particles i, j being next to each other on the cycle, for which we prove some cases. We also find that two natural events associated to the process have exactly the same probability expressed as a Vandermonde determinant. It is unclear whether this is just a coincidence or a consequence of a deeper connection.


title: Template-space metric for searches for gravitational waves from the inspiral, merger and ringdown of binary black holes

abstract: Searches for gravitational waves (GWs) from binary black holes using interferometric GW detectors require the construction of template banks for performing matched filtering while analyzing the data. Placement of templates over the parameter space of binaries, as well as coincidence tests of GW triggers from multiple detectors make use of the definition of a metric over the space of gravitational waveforms. Although recent searches have employed waveform templates coherently describing the inspiral, merger and ringdown (IMR) of the coalescence, the metric used in the template banks and coincidence tests was derived from post-Newtonian inspiral waveforms. In this paper, we compute the template-space metric of the IMR waveform family IMRPhenomB over the parameter space of masses and the effective spin parameter. We also propose a coordinate system, which is a modified version of post-Newtonian chirp time coordinates, in which the metric is slowly varying over the parameter space. The match function analytically computed using the metric has excellent agreement with the "exact" match function computed numerically. We show that the metric is able to provide a reasonable approximation to the match function of other IMR waveform families, such that the effective-one-body model calibrated to numerical relativity (EOBNRv2). The availability of this metric can contribute to improving the sensitivity of searches for GWs from binary black holes in the advanced detector era.


title: Fully Bayesian binary Markov random field models: Prior specification and posterior simulation

abstract: We propose a flexible prior model for the parameters of binary Markov random fields (MRF) defined on rectangular lattices and with maximal cliques defined from a template maximal clique. The prior model allows higher-order interactions to be included. We also define a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to sample from the associated posterior distribution. The number of possible parameters for an MRF with for instance k x l maximal cliques becomes high even for small values of k and l. To get a flexible model which may adapt to the structure of a particular observed image we do not put any absolute restrictions on the parametrisation. Instead we define a parametric form for the MRF where the parameters have interpretation as potentials for the various clique configurations, and limit the effective number of parameters by assigning apriori discrete probabilities for events where groups of parameter values are equal. To run our RJMCMC algorithm we have to cope with the computationally intractable normalising constant of MRFs. For this we adopt a previously defined approximation for binary MRFs, but we also briefly discuss other alternatives. We demonstrate the flexibility of our prior formulation with simulated and real data examples.


title: Longitudinal high-dimensional principal components analysis with application to diffusion tensor imaging of multiple sclerosis

abstract: We develop a flexible framework for modeling high-dimensional imaging data observed longitudinally. The approach decomposes the observed variability of repeatedly measured high-dimensional observations into three additive components: a subject-specific imaging random intercept that quantifies the cross-sectional variability, a subject-specific imaging slope that quantifies the dynamic irreversible deformation over multiple realizations, and a subject-visit-specific imaging deviation that quantifies exchangeable effects between visits. The proposed method is very fast, scalable to studies including ultrahigh-dimensional data, and can easily be adapted to and executed on modest computing infrastructures. The method is applied to the longitudinal analysis of diffusion tensor imaging (DTI) data of the corpus callosum of multiple sclerosis (MS) subjects. The study includes $176$ subjects observed at $466$ visits. For each subject and visit the study contains a registered DTI scan of the corpus callosum at roughly 30,000 voxels.


title: Codimension one attracting sets in $\mathbb{P}^k(\mathbb{C})$

abstract: We are interested in attracting sets of $\mathbb{P}^k(\mathbb{C})$ which are of small topological degree and of codimension $1.$ We first show that there exists a large family of examples. Then we study their ergodic and pluripotential theoretic properties.


title: Limit laws on extremes of non-homogeneous Gaussian random fields

abstract: In this paper, by using the exact tail asymptotics derived by Debicki, Hashorva and Ji (Ann. Probab. 2014), we proved the Gumbel limit theorem for the maximum of a class of non-homogeneous Gaussian random fields. By using the obtained results, we also derived the Gumbel laws for Shepp statistics of fractional Brownian motion and Gaussian integrated process as well as the Gumbel law for Storage process with fractional Brownian motion as input.


title: Photochemical Escape of Oxygen from Early Mars

abstract: Photochemical escape is an important process for oxygen escape from present Mars. In this work, a 1-D Monte-Carlo Model is developed to calculate escape rates of energetic oxygen atoms produced from O2+ dissociative recombination reactions (DR) under 1, 3, 10, and 20 times present solar XUV fluxes. We found that although the overall DR rates increase with solar XUV flux almost linearly, oxygen escape rate increases from 1 to 10 times present solar XUV conditions but decreases when increasing solar XUV flux further. Analysis shows that atomic species in the upper thermosphere of early Mars increases more rapidly than O2+ when increasing XUV fluxes. While the latter is the source of energetic O atoms, the former increases the collision probability and thus decreases the escape probability of energetic O. Our results suggest that photochemical escape be a less important escape mechanism than previously thought for the loss of water and/or CO2 from early Mars.


title: Search for the U-boson in the process $e^ +e^- \to \mu^+ \mu^- \gamma, \,\, U \to \mu^+\mu^-$ with the KLOE detector

abstract: We present a search for a new light vector boson, carrier of a "dark force" between WIMPs, with the KLOE detector at \DA$\Phi$NE. We analyzed $e^+ e^- \to \mu^+ \mu^- \gamma$ ISR events corresponding to an integrated luminosity of $239$ pb$^{-1}$ to find evidence for the $e^+ e^- \to U\gamma ,\,\, U\to\mu^+\mu^-$ process. We found no $U$ vector boson signal and set a 90% CL upper limit on the ratio of the U boson and photon coupling constants between 1.6$\times10^{-5}$ to 8.6$\times10^{-7}$ in the mass region $520<M_{\rm U}<980$ MeV. A projection of the KLOE sensitivity for the $\mu\mu\gamma$ and $\pi\pi\gamma$ channels at full statistics and extended muon acceptance is also presented.


title: 21cm Forest with the SKA

abstract: An alternative to both the tomography technique and the power spectrum approach is to search for the 21cm forest, that is the 21cm absorption features against high-z radio loud sources caused by the intervening cold neutral intergalactic medium (IGM) and collapsed structures. Although the existence of high-z radio loud sources has not been confirmed yet, SKA-low would be the instrument of choice to find such sources as they are expected to have spectra steeper than their lower-z counterparts. Since the strongest absorption features arise from small scale structures (few tens of physical kpc, or even lower), the 21cm forest can probe the HI density power spectrum on small scales not amenable to measurements by any other means. Also, it can be a unique probe of the heating process and the thermal history of the early universe, as the signal is strongly dependent on the IGM temperature. Here we show what SKA1-low could do in terms of detecting the 21cm forest in the redshift range z = 7.5-15.


title: Comment on `Update of 40K and 226Ra and 232Th series $\gamma$-to-dose conversion factors for soil'

abstract: A letter to the editor of the Journal of Environmental Radioactivity on the article: E. Gasser, A. Nachab, A. Nourreddine, Ch. Roy, and A. Sellam, `Update of 40K and 226Ra and 232Th series $\gamma$-to-dose conversion factors for soil', J. Environ. Radioactiv. 138, 68-71 (2014), DOI: 10.1016/j.jenvrad.2014.08.002.


title: Quantum dynamics of propagating photons with strong interactions: a generalized input-output formalism

abstract: There has been rapid development of systems that yield strong interactions between freely propagating photons in one dimension via controlled coupling to quantum emitters. This raises interesting possibilities such as quantum information processing with photons or quantum many-body states of light, but treating such systems generally remains a difficult task theoretically. Here, we describe a novel technique in which the dynamics and correlations of a few photons can be exactly calculated, based upon knowledge of the initial photonic state and the solution of the reduced effective dynamics of the quantum emitters alone. We show that this generalized "input-output" formalism allows for a straightforward numerical implementation regardless of system details, such as emitter positions, external driving, and level structure. As a specific example, we apply our technique to show how atomic systems with infinite-range interactions and under conditions of electromagnetically induced transparency enable the selective transmission of correlated multi-photon states.


title: Predicting binding free energies in solution

abstract: Recent predictions of absolute binding free energies of host-guest complexes in aqueous solution using electronic structure theory have been encouraging for some systems, while other systems remain problematic for others. In paper I summarize some of the many factors that could easily contribute 1-3 kcal/mol errors at 298 K: three-body dispersion effects, molecular symmetry, anharmonicity, spurious imaginary frequencies, insufficient conformational sampling, wrong or changing ionization states, errors in the solvation free energy of ions, and explicit solvent (and ion) effects that are not well-represented by continuum models. While I focus on binding free energies in aqueous solution the approach also applies (with minor adjustments) to any free energy difference such as conformational or reaction free energy differences or activation free energies in any solvent.


title: Cosmic Dawn and Epoch of Reionization Foreground Removal with the SKA

abstract: The exceptional sensitivity of the SKA will allow observations of the Cosmic Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both spectrally and spatially. This wealth of information is buried under Galactic and extragalactic foregrounds, which must be removed accurately and precisely in order to reveal the cosmological signal. This problem has been addressed already for the previous generation of radio telescopes, but the application to SKA is different in many aspects. In this chapter we summarise the contributions to the field of foreground removal in the context of high redshift and high sensitivity 21-cm measurements. We use a state-of-the-art simulation of the SKA Phase 1 observations complete with cosmological signal, foregrounds and frequency-dependent instrumental effects to test both parametric and non-parametric foreground removal methods. We compare the recovered cosmological signal using several different statistics and explore one of the most exciting possibilities with the SKA --- imaging of the ionized bubbles. We find that with current methods it is possible to remove the foregrounds with great accuracy and to get impressive power spectra and images of the cosmological signal. The frequency-dependent PSF of the instrument complicates this recovery, so we resort to splitting the observation bandwidth into smaller segments, each of a common resolution. If the foregrounds are allowed a random variation from the smooth power law along the line of sight, methods exploiting the smoothness of foregrounds or a parametrization of their behaviour are challenged much more than non-parametric ones. However, we show that correction techniques can be implemented to restore the performances of parametric approaches, as long as the first-order approximation of a power law stands.


title: Direct measurements of the magnetocaloric effect in pulsed magnetic fields: The example of the Heusler alloy Ni$_{50}$Mn$_{35}$In$_{15}$

abstract: We have studied the magnetocaloric effect (MCE) in the shape-memory Heusler alloy Ni$_{50}$Mn$_{35}$In$_{15}$ by direct measurements in pulsed magnetic fields up to 6 and 20 T. The results in 6 T are compared with data obtained from heat-capacity experiments. We find a saturation of the inverse MCE, related to the first-order martensitic transition, with a maximum adiabatic temperature change of $\Delta T_{ad} = -7$ K at 250 K and a conventional field-dependent MCE near the second-order ferromagnetic transition in the austenitic phase. The pulsed magnetic field data allow for an analysis of the temperature response of the sample to the magnetic field on a time scale of $\sim 10$ to 100 ms which is on the order of typical operation frequencies (10 to 100 Hz) of magnetocaloric cooling devices. Our results disclose that in shape-memory alloys the different contributions to the MCE and hysteresis effects around the martensitic transition have to be carefully considered for future cooling applications.


title: Field-normalized citation impact indicators and the choice of an appropriate counting method

abstract: Bibliometric studies often rely on field-normalized citation impact indicators in order to make comparisons between scientific fields. We discuss the connection between field normalization and the choice of a counting method for handling publications with multiple co-authors. Our focus is on the choice between full counting and fractional counting. Based on an extensive theoretical and empirical analysis, we argue that properly field-normalized results cannot be obtained when full counting is used. Fractional counting does provide results that are properly field normalized. We therefore recommend the use of fractional counting in bibliometric studies that require field normalization, especially in studies at the level of countries and research organizations. We also compare different variants of fractional counting. In general, it seems best to use either the author-level or the address-level variant of fractional counting.


title: Quantum group-twisted tensor products of C*-algebras II

abstract: For a quasitriangular C*-quantum group, we enrich the twisted tensor product constructed in the first part of this series to a monoidal structure on the category of its continuous coactions on C*-algebras. We define braided C*-quantum groups, where the comultiplication takes values in a twisted tensor product. We show that compact braided C*-quantum groups yield compact quantum groups by a semidirect product construction.


title: A note on weighted badly approximable linear forms

abstract: We prove a result in the area of twisted Diophantine approximation related to the theory of Schmidt games. In particular, under certain restrictions we give a affirmative answer to the analogue in this setting of a famous conjecture of Schmidt from Diophantine approximation.


title: "`They brought in the horrible key ring thing!" Analysing the Usability of Two-Factor Authentication in UK Online Banking

abstract: To prevent password breaches and guessing attacks, banks increasingly turn to two-factor authentication (2FA), requiring users to present at least one more factor, such as a one-time password generated by a hardware token or received via SMS, besides a password. We can expect some solutions -- especially those adding a token -- to create extra work for users, but little research has investigated usability, user acceptance, and perceived security of deployed 2FA. This paper presents an in-depth study of 2FA usability with 21 UK online banking customers, 16 of whom had accounts with more than one bank. We collected a rich set of qualitative and quantitative data through two rounds of semi-structured interviews, and an authentication diary over an average of 11 days. Our participants reported a wide range of usability issues, especially with the use of hardware tokens, showing that the mental and physical workload involved shapes how they use online banking. Key targets for improvements are (i) the reduction in the number of authentication steps, and (ii) removing features that do not add any security but negatively affect the user experience.


title: Swift follow-up of IceCube triggers, and implications for the Advanced-LIGO era

abstract: Between 2011 March and 2014 August Swift responded to 20 triggers from the IceCube neutrino observatory, observing the IceCube 50% confidence error circle in X-rays, typically within 5 hours of the trigger. No confirmed counterpart has been detected. We describe the Swift follow up strategy and data analysis and present the results of the campaign. We discuss the challenges of distinguishing the X-ray counterpart to a neutrino trigger from serendipitous uncatalogued X-ray sources in the error circle, and consider the implications of our results for future strategies for multi-messenger astronomy, with particular reference to the follow up of gravitational wave triggers from the advanced-era detectors.


title: Numerical study of the depinning transition of a ferromagnetic magnetic domain wall in films

abstract: We report first principle numerical study of domain wall (DW) depinning in two-dimensional magnetic film, which is modeled by 2D random-field Ising system with the dipole-dipole interaction. We observe nonconventional activation-type motion of DW and reveal its fractal structure of DW near the depinning transition. We determine scaling functions describing critical dynamics near the transition and obtain universal exponents establishing connection between thermal softening of pinning potential and critical dynamics. We observe that tuning the strength of the dipole-dipole interaction switches DW dynamics between two different universality classes corresponding to two distinct dynamic regimes, motion in the random potential and that in the random force.


title: A Wasserstein gradient flow approach to Poisson-Nernst-Planck equations

abstract: The Poisson-Nernst-Planck system of equations used to model ionic transport is interpreted as a gradient flow for the Wasserstein distance and a free energy in the space of probability measures with finite second moment. A variational scheme is then set up and is the starting point of the construction of global weak solutions in a unified framework for the cases of both linear and nonlinear diffusion. The proof of the main results relies on the derivation of extra estimates based on the flow interchange technique developed by Matthes, McCann, and Savar\'e.


title: Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation

abstract: We introduce an accurate and efficient method for a class of nonlocal potential evaluations with free boundary condition, including the 3D/2D Coulomb, 2D Poisson and 3D dipolar potentials. Our method is based on a Gaussian-sum approximation of the singular convolution kernel and Taylor expansion of the density. Starting from the convolution formulation, for smooth and fast decaying densities, we make a full use of the Fourier pseudospectral (plane wave) approximation of the density and a separable Gaussian-sum approximation of the kernel in an interval where the singularity (the origin) is excluded. Hence, the potential is separated into a regular integral and a near-field singular correction integral, where the first integral is computed with the Fourier pseudospectral method and the latter singular one can be well resolved utilizing a low-order Taylor expansion of the density. Both evaluations can be accelerated by fast Fourier transforms (FFT). The new method is accurate (14-16 digits), efficient ($O(N \log N)$ complexity), low in storage, easily adaptable to other different kernels, applicable for anisotropic densities and highly parallelable.


title: Finite field calculations of static polarizabilities and hyperpolarizabilities of In$^{+}$ and Sr

abstract: The finite field calculations are performed for two heavy frequency-standard candidates In$^+$ and Sr. The progressive hierarchy of electron correlations is implemented by the relativistic coupled-cluster and configuration interaction methods combined with basis set of increasing size. The dipole polarizabilities, dipole hyperpolarizabilities, quadrupole moments, and quadrupole polarizabilities are recommended for the ground state 5s$^2$ $^1S_0$ and low-lying states 5s5p $^3P^{\rm o}_{0,1,2}$ of In$^+$ and Sr. Comparative study of the fully and scalar relativistic electron correlation calculations reveals the effect of the spin-orbit interaction on the dipole polarizabilities of In$^{+}$ and Sr. Finally, the blackbody-radiation shifts due to the dipole polarizability, dipole hyperpolarizability, and quadrupole polarizability are evaluated for the clock transition 5s$^2$ $^1S_0$ - 5s5p $^3P^{\rm o}_0$ of In$^+$ and Sr.