title: Topology by Design in Magnetic nano-Materials: Artificial Spin Ice

abstract: Artificial Spin Ices are two dimensional arrays of magnetic, interacting nano-structures whose geometry can be chosen at will, and whose elementary degrees of freedom can be characterized directly. They were introduced at first to study frustration in a controllable setting, to mimic the behavior of spin ice rare earth pyrochlores, but at more useful temperature and field ranges and with direct characterization, and to provide practical implementation to celebrated, exactly solvable models of statistical mechanics previously devised to gain an understanding of degenerate ensembles with residual entropy. With the evolution of nano--fabrication and of experimental protocols it is now possible to characterize the material in real-time, real-space, and to realize virtually any geometry, for direct control over the collective dynamics. This has recently opened a path toward the deliberate design of novel, exotic states, not found in natural materials, and often characterized by topological properties. Without any pretense of exhaustiveness, we will provide an introduction to the material, the early works, and then, by reporting on more recent results, we will proceed to describe the new direction, which includes the design of desired topological states and their implications to kinetics.


title: Binary Bouncy Particle Sampler

abstract: The Bouncy Particle Sampler is a novel rejection-free non-reversible sampler for differentiable probability distributions over continuous variables. We generalize the algorithm to piecewise differentiable distributions and apply it to generic binary distributions using a piecewise differentiable augmentation. We illustrate the new algorithm in a binary Markov Random Field example, and compare it to binary Hamiltonian Monte Carlo. Our results suggest that binary BPS samplers are better for easy to mix distributions.


title: Fractional differential equations with mixed boundary conditions

abstract: In this paper, we discuss the existence and uniqueness of solutions of a boundary value problem for a fractional differential equation of order $\alpha\in(2,3)$, involving a general form of fractional derivative. First, we prove an equivalence between the Cauchy problem and the Volterra equation. Then, two results on the existence of solutions are proven, and we end with some illustrative examples.


title: Approximation properties of q-Bernoulli polynomials

abstract: We study the q-analogue of Euler-Maclaurin formula and by introducing a new q-operator we drive to this form. Moreover, approximation properties of q-Bernoulli polynomials is discussed. We estimate the suitable functions as a combination of truncated series of q-Bernoulli polynomials and the error is calculated. This paper can be helpful in two different branches, first solving differential equations by estimating functions and second we may apply these techniques for operator theory.


title: The scaling properties and the multiple derivative of Legendre polynomials

abstract: In this paper, we study the scaling properties of Legendre polynomials Pn(x). We show that Pn(ax), where a is a constant, can be expanded as a sum of either Legendre polynomials Pn(x) or their multiple derivatives dkPn(x)/dxk, and we derive a general expression for the expansion coefficients. In addition, we demonstrate that the multiple derivative dkPn(x)/dxk can also be expressed as a sum of Legendre polynomials and we obtain a recurrence relation for the coefficients.


title: On choosing the start time of binary black hole ringdown

abstract: The final stage of a binary black hole merger is ringdown, in which the system is described by a Kerr black hole with quasinormal mode perturbations. It is far from straightforward to identify the time at which the ringdown begins. Yet determining this time is important for precision tests of the general theory of relativity that compare an observed signal with quasinormal mode descriptions of the ringdown, such as tests of the no-hair theorem. We present an algorithmic method to analyze the choice of ringdown start time in the observed waveform. This method is based on determining how close the strong field is to a Kerr black hole (Kerrness). Using numerical relativity simulations, we characterize the Kerrness of the strong-field region close to the black hole using a set of local, gauge-invariant geometric and algebraic conditions that measure local isometry to Kerr. We produce a map that associates each time in the gravitational waveform with a value of each of these Kerrness measures; this map is produced by following outgoing null characteristics from the strong and near-field regions to the wave zone. We perform this analysis on a numerical relativity simulation with parameters consistent with GW150914- the first gravitational wave detection. We find that the choice of ringdown start time of $3\,\mathrm{ms}$ after merger used in the GW150914 study to test general relativity corresponds to a high dimensionless perturbation amplitude of $ \sim 7.5 \times 10^{-3}$ in the strong-field region. This suggests that in higher signal-to-noise detections, one would need to start analyzing the signal at a later time for studies that depend on the validity of black hole perturbation theory.


title: Audio Set classification with attention model: A probabilistic perspective

abstract: This paper investigates the classification of the Audio Set dataset. Audio Set is a large scale weakly labelled dataset of sound clips. Previous work used multiple instance learning (MIL) to classify weakly labelled data. In MIL, a bag consists of several instances, and a bag is labelled positive if at least one instances in the audio clip is positive. A bag is labelled negative if all the instances in the bag are negative. We propose an attention model to tackle the MIL problem and explain this attention model from a novel probabilistic perspective. We define a probability space on each bag, where each instance in the bag has a trainable probability measure for each class. Then the classification of a bag is the expectation of the classification output of the instances in the bag with respect to the learned probability measure. Experimental results show that our proposed attention model modeled by fully connected deep neural network obtains mAP of 0.327 on Audio Set dataset, outperforming the Google's baseline of 0.314 and recurrent neural network of 0.325.


title: Physiological Tremor Increases when Skeletal Muscle is Shortened: Implications for Fusimotor Control

abstract: The involuntary force fluctuations associated with physiological (as distinct from pathological) tremor are an unavoidable component of human motor control. While the origins of the physiological tremor are known to depend on muscle afferentation, it is possible that the mechanical properties of muscle-tendon systems also affect its generation, amplification and maintenance. In this paper, we investigated the dependence of physiological tremor during tonic, isometric plantarflextion torque at 30% of maximum at three ankle angles. The amplitude of physiological tremor increased as calf muscles shortened in contrast to the stretch reflex whose amplitude decreases as muscle shortens. We used a closed-loop simulation model of afferented muscle to explore the mechanisms responsible for this behavior. We demonstrate that changing muscle lengths does not suffice to explain our experimental findings. Rather, the model consistently required the modulation of gamme-static fusimotor drive to produce increases in physiological tremor with muscle shortening--while successfully replicating the concomitant reduction in stretch reflex amplitude. This need to control gamma-static fusimotor drive explicitly as a function of muscle length has important implication. First, it permits the amplitudes of physiological tremor and stretch reflex to be decoupled. Second, it postulates neuromechanical interaction that require length-dependent gamma drive modulation to be independent from alpha drive to the parent muscle. Lastly, it suggests that physiological tremor can be used as a simple, non-invasive measure of the afferent mechanisms underlying healthy motor function, and their disruption in neurological conditions.


title: Positive projectively flat manifolds are locally conformally flat-K\"ahler Hopf manifolds

abstract: We define a partition of the space of projectively flat metrics in three classes according to the sign of the Chern scalar curvature; we prove that the class of negative projectively flat metrics is empty, and that the class of positive projectively flat metrics consists precisely of locally conformally flat-K\"ahler metrics on Hopf manifolds, explicitly characterized by Vaisman. Finally, we review the properties of zero projectively flat metrics. As applications, we refine a list of possible projectively flat metrics by Li, Yau, and Zheng; moreover we prove that projectively flat astheno-K\"ahler metrics are in fact K\"ahler and globally conformally flat.


title: A new Measurement of the Intergalactic Temperature at $z\sim 2.55-2.95$

abstract: We present two measurements of the temperature-density relationship (TDR) of the intergalactic medium (IGM) in the redshift range $2.55 < z < 2.95$ using a sample of 13 high-quality quasar spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{HI}$ and the Doppler parameter $b$ of the absorption lines in the \mlya\ forest. The first measurement is obtained using a novel Bayesian scheme which takes into account the statistical correlations between the parameters characterising the lower cut-off of the $b-N_{HI}$ distribution and the power-law parameters $T_0$ and $\gamma$ describing the TDR. This approach yields $T_0/ 10^3\, {\rm K}=15.6 \pm 4.4 $ and $\gamma=1.45 \pm 0.17$ independent of the assumed pressure smoothing of the small scale density field. In order to explore the information contained in the overall $b-N_{HI}$ distribution rather than only the lower cut-off, we obtain a second measurement based on a similar Bayesian analysis of the median Doppler parameter for separate column-density ranges of the absorbers. In this case we obtain $T_0/ 10^3\, {\rm K}=14.6 \pm 3.7$ and $\gamma=1.37 \pm 0.17$ in good agreement with the first measurement. Our Bayesian analysis reveals strong anti-correlations between the inferred $T_0$ and $\gamma$ for both methods as well as an anti-correlation of the inferred $T_0$ and the pressure smoothing length for the second method, suggesting that the measurement accuracy can in the latter case be substantially increased if independent constraints on the smoothing are obtained. Our results are in good agreement with other recent measurements of the thermal state of the IGM probing similar (over-)density ranges.


title: A Denotational Semantics for SPARC TSO

abstract: The SPARC TSO weak memory model is defined axiomatically, with a non-compositional formulation that makes modular reasoning about programs difficult. Our denotational approach uses pomsets to provide a compositional semantics capturing exactly the behaviours permitted by SPARC TSO. It uses buffered states and an inductive definition of execution to assign an input-output meaning to pomsets. We show that our denotational account is sound and complete relative to the axiomatic account, that is, that it captures exactly the behaviours permitted by the axiomatic account. Our compositional approach facilitates the study of SPARC TSO and supports modular analysis of program behaviour.


title: Interior C2 regularity of convex solutions to prescribing scalar curvature equations

abstract: We establish interior $C^2$ estimates for convex solutions of scalar curvature equation and $\sigma_2$-Hessian equation. We also prove interior curvature estimate for isometrically immersed hypersurfaces $(M^n,g)\subset \mathbb R^{n+1}$ with positive scalar curvature. These estimates are consequences of an interior estimate for these equations obtained under a weakened condition.


title: BPS Sphalerons in the $F_2$ Non-Linear Sigma Model

abstract: We construct static and also time-dependent solutions in a non-linear sigma model with target space being the flag manifold $F_2=SU(3)/U(1)^2$ on the four dimensional Minkowski space-time by analytically solving the second order Euler-Lagrange equation. We show the static solutions saturate an energy lower bound and can be derived from coupled first order equations though they are saddle point solutions. We also discuss basic properties of the time-dependent solutions.


title: A manual proof of the Four-colour Theorem

abstract: In this paper, we provide an easy proof of the Four-colour Theorem in a special case indeed.


title: A Note on Robust Biarc Computation

abstract: A new robust algorithm for the numerical computation of biarcs, i.e. $G^1$ curves composed of two arcs of circle, is presented. Many algorithms exist but are based on geometric constructions, which must consider many geometrical configurations. The proposed algorithm uses an algebraic construction which is reduced to the solution of a single $2$ by $2$ linear system. Singular angles configurations are treated smoothly by using the pseudoinverse matrix when solving the linear system. The proposed algorithm is compared with the Matlab's routine \texttt{rscvn} that solves geometrically the same problem. Numerical experiments show that Matlab's routine sometimes fails near singular configurations and does not select the correct solution for large angles, whereas the proposed algorithm always returns the correct solution. The proposed solution smoothly depends on the geometrical parameters so that it can be easily included in more complex algorithms like splines of biarcs or least squares data fitting.


title: NMR evidence for static local nematicity and its cooperative interplay with low-energy magnetic fluctuations in FeSe under pressure

abstract: We present $^{77}$Se-NMR measurements on single-crystalline FeSe under pressures up to 2 GPa. Based on the observation of the splitting and broadening of the NMR spectrum due to structural twin domains, we discovered that static, local nematic ordering exists well above the bulk nematic ordering temperature, $T_{\rm s}$. The static, local nematic order and the low-energy stripe-type antiferromagnetic spin fluctuations, as revealed by NMR spin-lattice relaxation rate measurements, are both insensitive to pressure application. These NMR results provide clear evidence for the microscopic cooperation between magnetism and local nematicity in FeSe.


title: Neural Discrete Representation Learning

abstract: Learning useful representations without supervision remains a key challenge in machine learning. In this paper, we propose a simple yet powerful generative model that learns such discrete representations. Our model, the Vector Quantised-Variational AutoEncoder (VQ-VAE), differs from VAEs in two key ways: the encoder network outputs discrete, rather than continuous, codes; and the prior is learnt rather than static. In order to learn a discrete latent representation, we incorporate ideas from vector quantisation (VQ). Using the VQ method allows the model to circumvent issues of "posterior collapse" -- where the latents are ignored when they are paired with a powerful autoregressive decoder -- typically observed in the VAE framework. Pairing these representations with an autoregressive prior, the model can generate high quality images, videos, and speech as well as doing high quality speaker conversion and unsupervised learning of phonemes, providing further evidence of the utility of the learnt representations.


title: A Comparison of Feature-Based and Neural Scansion of Poetry

abstract: Automatic analysis of poetic rhythm is a challenging task that involves linguistics, literature, and computer science. When the language to be analyzed is known, rule-based systems or data-driven methods can be used. In this paper, we analyze poetic rhythm in English and Spanish. We show that the representations of data learned from character-based neural models are more informative than the ones from hand-crafted features, and that a Bi-LSTM+CRF-model produces state-of-the art accuracy on scansion of poetry in two languages. Results also show that the information about whole word structure, and not just independent syllables, is highly informative for performing scansion.


title: Deep Air Learning: Interpolation, Prediction, and Feature Analysis of Fine-grained Air Quality

abstract: The interpolation, prediction, and feature analysis of fine-gained air quality are three important topics in the area of urban air computing. The solutions to these topics can provide extremely useful information to support air pollution control, and consequently generate great societal and technical impacts. Most of the existing work solves the three problems separately by different models. In this paper, we propose a general and effective approach to solve the three problems in one model called the Deep Air Learning (DAL). The main idea of DAL lies in embedding feature selection and semi-supervised learning in different layers of the deep learning network. The proposed approach utilizes the information pertaining to the unlabeled spatio-temporal data to improve the performance of the interpolation and the prediction, and performs feature selection and association analysis to reveal the main relevant features to the variation of the air quality. We evaluate our approach with extensive experiments based on real data sources obtained in Beijing, China. Experiments show that DAL is superior to the peer models from the recent literature when solving the topics of interpolation, prediction, and feature analysis of fine-gained air quality.


title: Monotone bargaining is Nash-solvable

abstract: Given two finite ordered sets $A = \{a_1, \ldots, a_m\}$ and $B = \{b_1, \ldots, b_n\}$, introduce the set of $m n$ outcomes of the game $O = \{(a, b) \mid a \in A, b \in B\} = \{(a_i, b_j) \mid i \in I = \{1, \ldots, m\}, j \in J = \{1, \ldots, n\}$. Two players, Alice and Bob, have the sets of strategies $X$ and $Y$ that consist of all monotone non-decreasing mappings $x: A \rightarrow B$ and $y: B \rightarrow A$, respectively. It is easily seen that each pair $(x,y) \in X \times Y$ produces at least one {\em deal}, that is, an outcome $(a,b) \in O$ such that $x(a) = b$ and $y(b) = a$. Denote by $G(x,y) \subseteq O$ the set of all such deals related to $(x,y)$. The obtained mapping $G = G_{m,n}: X \times Y \rightarrow 2^O$ is a game correspondence. Choose an arbitrary deal $g(x,y) \in G(x,y)$ to obtained a mapping $g : X \times Y \rightarrow O$, which is a game form. We will show that each such game form is tight and, hence, Nash-solvable, that is, for any pair $u = (u_A, u_B)$ of utility functions $u_A : O \rightarrow \mathbb R$ of Alice and $u_B: O \rightarrow \mathbb R$ of Bob, the obtained monotone bargaining game $(g, u)$ has at least one Nash equilibrium in pure strategies. Moreover, the same equilibrium can be chosen for all selections $g(x,y) \in G(x,y)$. We also obtain an efficient algorithm that determines such an equilibrium in time linear in $m n$, although the numbers of strategies $|X| = \binom{m+n-1}{m}$ and $|Y| = \binom{m+n-1}{n}$ are exponential in $m n$. Our results show that, somewhat surprising, the players have no need to hide or randomize their bargaining strategies, even in the zero-sum case.


title: Deep Active Learning over the Long Tail

abstract: This paper is concerned with pool-based active learning for deep neural networks. Motivated by coreset dataset compression ideas, we present a novel active learning algorithm that queries consecutive points from the pool using farthest-first traversals in the space of neural activation over a representation layer. We show consistent and overwhelming improvement in sample complexity over passive learning (random sampling) for three datasets: MNIST, CIFAR-10, and CIFAR-100. In addition, our algorithm outperforms the traditional uncertainty sampling technique (obtained using softmax activations), and we identify cases where uncertainty sampling is only slightly better than random sampling.


title: Compressing Large-Scale Wave Propagation Models via Phase-Preconditioned Rational Krylov Subspaces

abstract: Rational Krylov subspace (RKS) techniques are well-established and powerful tools for projection-based model reduction of time-invariant dynamic systems. For hyperbolic wavefield problems, such techniques perform well in configurations where only a few modes contribute to the field. RKS methods, however, are fundamentally limited by the Nyquist-Shannon sampling rate, making them unsuitable for the approximation of wavefields in configuration characterized by large travel times and propagation distances, since wavefield responses in such configurations are highly oscillatory in the frequency-domain. To overcome this limitation, we propose to precondition the RKSs by factoring out the rapidly varying frequency-domain field oscillations. The remaining amplitude functions are generally slowly varying functions of source position and spatial coordinate and allow for a significant compression of the approximation subspace. Our one-dimensional analysis together with numerical experiments for large scale 2D acoustic models show superior approximation properties of preconditioned RKS compared with the standard RKS model-order reduction. The preconditioned RKS results in a reduction of the frequency sampling well below the Nyquist-Shannon rate, a weak dependence of the RKS size on the number of inputs and outputs for multiple-input/multiple-output (MIMO) problems, and, most importantly, in a significant coarsening of the finite-difference grid used to generate the RKS. A prototype implementation indicates that the preconditioned RKS algorithm is competitive in the modern high performance computing environment.


title: Performance of 20:1 multiplexer for large area charge readouts in directional dark matter TPC detectors

abstract: More target mass is required in current TPC based directional dark matter detectors for improved detector sensitivity. This can be achieved by scaling up the detector volumes, but this results in the need for more analogue signal channels. A possible solution to reducing the overall cost of the charge readout electronics is to multiplex the signal readout channels. Here, we present a multiplexer system in expanded mode based on LMH6574 chips produced by Texas Instruments, originally designed for video processing. The setup has a capability of reducing the number of readouts in such TPC detectors by a factor of 20. Results indicate that the important charge distribution asymmetry along an ionization track is retained after multiplexed signals are demultiplexed.

categories: physics.ins-det astro-ph.IM hep-ex physics.app-ph


title: Controlled Growth of a Large-Size 2D Selenium Nanosheet and Its Electronic and Optoelectronic Applications

abstract: Selenium has attracted intensive attention as a promising material candidate for future optoelectronic applications. However, selenium has a strong tendency to grow into nanowire forms due to its anisotropic atomic structure, which has largely hindered the exploration of its potential applications. In this work, using a physical vapor deposition method, we have demonstrated the synthesis of large-size, high-quality 2D selenium nanosheets, the minimum thickness of which could be as thin as 5 nm. The Se nanosheet exhibits a strong in-plane anisotropic property, which is determined by angle-resolved Raman spectroscopy. Back-gating field-effect transistors based on a Se nanosheet exhibit p-type transport behaviors with on-state current density around 20 mA/mm at Vds = 3 V. Four-terminal field effect devices are also fabricated to evaluate the intrinsic hole mobility of the selenium nanosheet, and the value is determined to be 0.26 cm2 Vs at 300 K. The selenium nanosheet phototransistors show an excellent photoresponsivity of up to 263 A/W, with a rise time of 0.1 s and fall time of 0.12 s. These results suggest that crystal selenium as a 2D form of a 1D van der Waals solid opens up the possibility to explore device applications.


title: Partition function zeros of adsorbing Dyck paths

abstract: The zeros of the size-$n$ partition functions for a statistical mechanical model can be used to help understand the critical behaviour of the model as $n\to\infty$. Here we use weighted Dyck paths as a simple model of two-dimensional polymer adsorption, and study the behaviour of the partition function zeros, particularly in the thermodynamic limit. The exact solvability of the model allows for a precise calculation of the locus of the zeros and the way in which an edge-singularity on the positive real axis is formed. We also show that in the limit $n\to\infty$ the zeros converge on a lima\c{c}on in the complex plane.


title: Learning Linear Dynamical Systems via Spectral Filtering

abstract: We present an efficient and practical algorithm for the online prediction of discrete-time linear dynamical systems with a symmetric transition matrix. We circumvent the non-convex optimization problem using improper learning: carefully overparameterize the class of LDSs by a polylogarithmic factor, in exchange for convexity of the loss functions. From this arises a polynomial-time algorithm with a near-optimal regret guarantee, with an analogous sample complexity bound for agnostic learning. Our algorithm is based on a novel filtering technique, which may be of independent interest: we convolve the time series with the eigenvectors of a certain Hankel matrix.


title: Phase-Resolved Spectroscopy of the Low-Mass X-ray Binary V801 Ara

abstract: We present phase-resolved optical spectra of the low mass X-ray binary system V801 Ara. The spectra, obtained in 2014 with IMACS on the Magellan/Baade telescope at Las Campanas Observatory, cover the full binary orbit of 3.8 hours. They contain strong emission features allowing us to map the emission of HAlpha, HBeta, He II 4686, and the Bowen blend at 4640. The radial velocity curves of the Bowen blend shows significantly stronger modulation at the orbital period than HAlpha as expected for the former originating on the secondary with the latter consistent with emission dominated by the disk. Our tomograms of HAlpha and HBeta are the most detailed studies of these lines for V801 to date and they clearly detect the accretion disk. The HBeta emission extends to higher velocities than HAlpha, suggesting emission from closer to the neutron star and differentiating temperature variance in the accretion disk for the first time. The center of the accretion disk appears offset from the center-of-mass of the neutron star as has been seen in several other X-ray binaries. This is often interpreted to imply disk eccentricity. Our tomograms do not show strong evidence for a hot spot at the point where the accretion stream hits the disk. This could imply a reduced accretion rate or could be due to the spot being drowned out by bright accretion flow around it. There is enhanced emission further along the disk, however, which implies gas stream interaction downstream of the hot spot.


title: Interior Hessian estimates for Sigma-2 equations in dimension three

abstract: We prove a priori interior C2 estimate for \sigma_2 = f in R3, which generalizes Warren-Yuan's result.


title: Selective inference for the problem of regions via multiscale bootstrap

abstract: A general approach to selective inference is considered for hypothesis testing of the null hypothesis represented as an arbitrary shaped region in the parameter space of multivariate normal model. This approach is useful for hierarchical clustering where confidence levels of clusters are calculated only for those appeared in the dendrogram, thus subject to heavy selection bias. Our computation is based on a raw confidence measure, called bootstrap probability, which is easily obtained by counting how many times the same cluster appears in bootstrap replicates of the dendrogram. We adjust the bias of the bootstrap probability by utilizing the scaling-law in terms of geometric quantities of the region in the abstract parameter space, namely, signed distance and mean curvature. Although this idea has been used for non-selective inference of hierarchical clustering, its selective inference version has not been discussed in the literature. Our bias-corrected $p$-values are asymptotically second-order accurate in the large sample theory of smooth boundary surfaces of regions, and they are also justified for nonsmooth surfaces such as polyhedral cones. The $p$-values are asymptotically equivalent to those of the iterated bootstrap but with less computation.


title: Beyond normality: Learning sparse probabilistic graphical models in the non-Gaussian setting

abstract: We present an algorithm to identify sparse dependence structure in continuous and non-Gaussian probability distributions, given a corresponding set of data. The conditional independence structure of an arbitrary distribution can be represented as an undirected graph (or Markov random field), but most algorithms for learning this structure are restricted to the discrete or Gaussian cases. Our new approach allows for more realistic and accurate descriptions of the distribution in question, and in turn better estimates of its sparse Markov structure. Sparsity in the graph is of interest as it can accelerate inference, improve sampling methods, and reveal important dependencies between variables. The algorithm relies on exploiting the connection between the sparsity of the graph and the sparsity of transport maps, which deterministically couple one probability measure to another.


title: Superfluid Onset and Compressibility of $^4$He Films Adsorbed on Carbon Nanotubes

abstract: Third sound measurements of superfluid $^4$He thin films adsorbed on 10 nm diameter multiwall carbon nanotubes are used to probe the superfluid onset temperature as a function of the film thickness, and to study the temperature dependence of the film compressibility. The nanotubes provide a highly ordered carbon surface, with layer-by-layer growth of the adsorbed film as shown by oscillation peaks in the third sound velocity at the completion of the third, fourth, and fifth atomic layers, arising from oscillations in the compressibility. In temperature sweeps the third sound velocity at very low temperatures is found to be linear with temperature, but oscillating between positive and negative slope depending on the film thickness. Analysis shows that this can be attributed to a linearly decreasing compressibility of the film with temperature that appears to hold even near zero temperature. The superfluid onset temperature is found to be linear in the film thickness, as predicted by the Kosterlitz-Thouless theory, but the slope is anomalous, a factor of three smaller than the predicted universal value.


title: Radial terrace solutions and propagation profile of multistable reaction-diffusion equations over $\mathbb R^N$

abstract: We study the propagation profile of the solution $u(x,t)$ to the nonlinear diffusion problem $u_t-\Delta u=f(u)\; (x\in \mathbb R^N,\;t>0)$, $u(x,0)=u_0(x) \; (x\in\mathbb R^N)$, where $f(u)$ is of multistable type: $f(0)=f(p)=0$, $f'(0)<0$, $f'(p)<0$, where $p$ is a positive constant, and $f$ may have finitely many nondegenerate zeros in the interval $(0, p)$. The class of initial functions $u_0$ includes in particular those which are nonnegative and decay to 0 at infinity. We show that, if $u(\cdot, t)$ converges to $p$ as $t\to\infty$ in $L^\infty_{loc}(\mathbb R^N)$, then the long-time dynamical behavior of $u$ is determined by the one dimensional propagating terraces introduced by Ducrot, Giletti and Matano [DGM]. For example, we will show that in such a case, in any given direction $\nu\in\mathbb{S}^{N-1}$, $u(x\cdot \nu, t)$ converges to a pair of one dimensional propagating terraces, one moving in the direction of $x\cdot \nu>0$, and the other is its reflection moving in the opposite direction $x\cdot\nu<0$. Our approach relies on the introduction of the notion "radial terrace solution", by which we mean a special solution $V(|x|, t)$ of $V_t-\Delta V=f(V)$ such that, as $t\to\infty$, $V(r,t)$ converges to the corresponding one dimensional propagating terrace of [DGM]. We show that such radial terrace solutions exist in our setting, and the general solution $u(x,t)$ can be well approximated by a suitablly shifted radial terrace solution $V(|x|, t)$. These will enable us to obtain better convergence result for $u(x,t)$. We stress that $u(x,t)$ is a high dimensional solution without any symmetry. Our results indicate that the one dimensional propagating terrace is a rather fundamental concept; it provides the basic structure and ingredients for the long-time profile of solutions in all space dimensions.


title: Automatic Query Image Disambiguation for Content-Based Image Retrieval

abstract: Query images presented to content-based image retrieval systems often have various different interpretations, making it difficult to identify the search objective pursued by the user. We propose a technique for overcoming this ambiguity, while keeping the amount of required user interaction at a minimum. To achieve this, the neighborhood of the query image is divided into coherent clusters from which the user may choose the relevant ones. A novel feedback integration technique is then employed to re-rank the entire database with regard to both the user feedback and the original query. We evaluate our approach on the publicly available MIRFLICKR-25K dataset, where it leads to a relative improvement of average precision by 23% over the baseline retrieval, which does not distinguish between different image senses.


title: Efficient construction of tensor ring representations from sampling

abstract: In this paper we propose an efficient method to compress a high dimensional function into a tensor ring format, based on alternating least-squares (ALS). Since the function has size exponential in $d$ where $d$ is the number of dimensions, we propose efficient sampling scheme to obtain $O(d)$ important samples in order to learn the tensor ring. Furthermore, we devise an initialization method for ALS that allows fast convergence in practice. Numerical examples show that to approximate a function with similar accuracy, the tensor ring format provided by the proposed method has less parameters than tensor-train format and also better respects the structure of the original function.


title: Random Walks on Lattices. Influence of Competing Reaction Centers on Diffusion-Controlled Processes

abstract: We study diffusion-reaction processes on periodic square planar lattices and simple cubic (sc) lattices. Considered first is a single diffusing reactant undergoing an irreversible reaction upon first encounter with a stationary co-reactant ["one-walker (1W) problem"]. We then generalize this scenario to allow for a competing reaction, i.e., instantaneous trapping of the diffusing reactant with probability $s$ at any vacant site before interacting with a (stationary) co-reactant at a target site. We determine the mean walklength of the diffusing reactant until irreversible reaction occurs. We use generating functions and the theory of finite Markov processes, as well as MC simulations. To investigate the dependence of walklength on lattice size we compute the first, finite size corrections to the Green function of the sc lattice, and provide a Pad\'e approximation for this quantity. Finally, we consider the case where both reactant and co-reactant undergo synchronous nearest-neighbor displacements ["two-walker (2W) problem"]. In this case, reactant and co-reactant can individually be trapped with probability $s$ at any vacant lattice site, or can undergo an irreversible reaction on first encounter at any site. When $s=0$ we find that, both for the 1W and the 2W problem, for lattices with (approximately) the same number of sites, the mean walklength is smaller (and hence the reaction efficiency greater) in $d=3$ than in $d=2$. Increasing $s$ tends to reduce differences in system dimensionality, and distinctions between the 1W problem and the 2W problem. Our model provides a good starting point to develop studies on the efficiency of apparently diverse diffusion-reaction processes, such as diffusion on a partially poisoned catalytic substrate or photosynthetic trapping of excitations.


title: Running Time Analysis of the (1+1)-EA for OneMax and LeadingOnes under Bit-wise Noise

abstract: In many real-world optimization problems, the objective function evaluation is subject to noise, and we cannot obtain the exact objective value. Evolutionary algorithms (EAs), a type of general-purpose randomized optimization algorithm, have been shown to be able to solve noisy optimization problems well. However, previous theoretical analyses of EAs mainly focused on noise-free optimization, which makes the theoretical understanding largely insufficient for the noisy case. Meanwhile, the few existing theoretical studies under noise often considered the one-bit noise model, which flips a randomly chosen bit of a solution before evaluation; while in many realistic applications, several bits of a solution can be changed simultaneously. In this paper, we study a natural extension of one-bit noise, the bit-wise noise model, which independently flips each bit of a solution with some probability. We analyze the running time of the (1+1)-EA solving OneMax and LeadingOnes under bit-wise noise for the first time, and derive the ranges of the noise level for polynomial and super-polynomial running time bounds. The analysis on LeadingOnes under bit-wise noise can be easily transferred to one-bit noise, and improves the previously known results. Since our analysis discloses that the (1+1)-EA can be efficient only under low noise levels, we also study whether the sampling strategy can bring robustness to noise. We prove that using sampling can significantly increase the largest noise level allowing a polynomial running time, that is, sampling is robust to noise.


title: Thermally-Induced Structural Evolution of Silicon- and Oxygen-Containing Hydrogenated Amorphous Carbon

abstract: The thermally-induced structural evolution of silicon- and oxygen-containing hydrogenated amorphous carbon (a-C:H:Si:O) was investigated by X-ray photoelectron and absorption spectroscopy, as well as molecular dynamics (MD) simulations. The spectroscopic results indicate that the introduction of Si and O in hydrogenated amorphous carbon (a-C:H) increases the activation energy for the conversion of sp3- to sp2-bonded C. MD simulations indicate that the higher thermal stability of a-C:H:Si:O compared to a-C:H derives from the lower fraction of strained C-C sp3 bonds in a-C:H:Si:O.


title: Measurement of a Magnonic Crystal at Millikelvin Temperatures

abstract: Hybrid systems combining magnons and superconducting quantum circuits have attracted increasing interest in recent years. Magnonic crystals (MCs) are one of the building blocks of room-temperature magnonics and can be used to create devices with an engineered band structure. These devices, exhibiting tunable frequency selectivity and the ability to store travelling excitations in the microwave regime, may form the basis of a set of new tools to be used in the context of quantum information processing. In order to ascertain the feasibility of such plans, MCs must be demonstrated to work at the low temperatures required for microwave-frequency quantum experiments. We report the first measurements of the transmission of microwave signals through an MC at 20 mK and observe a magnonic bandgap in both continuous-wave and pulsed excitation experiments. The spin-wave damping at low temperatures in our yttrium iron garnet MC is higher than expected, indicating that further work is necessary before the full potential of quantum experiments using magnonic crystals can be realised.


title: Transition from Plan Driven to SAFe : Periodic Team Self-Assessment

abstract: Context: How to adopt, scale and tailor agile methods depends on several factors such as the size of the organization, business goals, operative model, and needs. The Scaled Agile Framework (SAFe) was developed to support organizations to scale agile practices across the enterprise. Problem: Early adopters of SAFe tend to be large multi-national enterprises who report that the adoption of SAFe has led to significant productivity and quality gains. However, little is known about whether these benefits translate to small to medium sized enterprises (SMEs). Method: As part of a longitudinal study of an SME transitioning to SAFe we ask, to what extent are SAFe practices adopted at the team level? We targeted all team members and administrated a mixed method survey in February, 2017 and in July, 2017 to identify and evaluate the adoption rate of SAFe practices. Results: Initially in Quarter 1, teams were struggling with PI/Release health and Technical health throughout the organization as most of the teams were transitioning from plan-driven to SAFe . But, during the transition period in Quarter 3, we observed discernible improvements in different areas of SAFe practice adoption. Conclusion: The observed improvement might be due to teams merely becoming more familiar with the practices over-time. However, management had also made some structural changes to the teams that may account for the change.


title: Isospin Dependence of the EMC effect and Short range Correlations

abstract: Recent studies have shown that the per-nucleon number of Short-Range Correlated (SRC) pairs in nuclei and the strength of the EMC effect are linearly correlated, increasing from light nuclei up to iron and then saturating. This paper shows that the per-proton number of SRC pairs and the strength of the EMC effect are linearly correlated and increase from light to heavy asymmetric neutron-rich nuclei without saturation. These quantities, calculated per neutron, are also linearly related, but saturate (have constant value) remarkably early, starting with 12C. We propose that the observed phenomenological relationships indicate an isospin dependence of the EMC effect which is associated with the dominance of SRCs by neutron-proton pairs from 3He to heavy asymmetric neutron-rich nuclei.


title: A note on a boundary sine-Gordon model at the free-Fermion point

abstract: We investigate the free-Fermion point of a boundary sine-Gordon model with nondiagonal boundary interactions for the ground state using auxiliary functions, obtained from T-Q equations of a corresponding inhomogeneous open spin-1/2 XXZ chain with nondiagonal boundary terms. In particular, we obtain the Casimir energy. Our result for the Casimir energy is shown to agree with the result from the TBA approach. The analytical result for the effective central charge in the ultraviolet (UV) limit is also verified from the plots of effective central charge for intermediate values of volume.


title: Energy-Delay Efficient Power Control in Wireless Networks

abstract: This work aims at developing a power control framework to jointly optimize energy efficiency (measured in bit/Joule) and delay in wireless networks. A multi-objective approach is taken to deal with both performance metrics, while ensuring a minimum quality-of-service to each user in the network. Each user in the network is modeled as a rational agent that engages in a generalized non-cooperative game. Feasibility conditions are derived for the existence of each player's best response, and used to show that if these conditions are met, the game best response dynamics will converge to a unique Nash equilibrium. Based on these results, a convergent power control algorithm is derived, which can be implemented in a fully decentralized fashion. Next, a centralized power control algorithm is proposed, which also serves as a benchmark for the proposed decentralized solution. Due to the non-convexity of the centralized problem, the tool of maximum block improvement is used, to trade-off complexity with optimality.


title: The Complexity of Finding Small Separators in Temporal Graphs

abstract: Temporal graphs are graphs with time-stamped edges. We study the problem of finding a small vertex set (the separator) with respect to two designated terminal vertices such that the removal of the set eliminates all temporal paths connecting one terminal to the other. Herein, we consider two models of temporal paths: paths that pass through arbitrarily many edges per time step (non-strict) and paths that pass through at most one edge per time step (strict). Regarding the number of time steps of a temporal graph, we show a complexity dichotomy (NP-hardness versus polynomial-time solvability) for both problem variants. Moreover we prove both problem variants to be NP-complete even on temporal graphs whose underlying graph is planar. We further show that, on temporal graphs with planar underlying graph, if additionally the number of time steps is constant, then the problem variant for strict paths is solvable in quasi-linear time. Finally, we introduce and motivate the notion of a temporal core (vertices whose incident edges change over time). We prove that the non-strict variant is fixed-parameter tractable when parameterized by the size of the temporal core, while the strict variant remains NP-complete, even for constant-size temporal cores.


title: Ergodicity breaking dynamics of arch collapse

abstract: Gravity driven flows such as in hoppers and silos are susceptible to clogging due to the formation of arches at the exit whose failure is the key to re-initiation of flow. In vibrated hoppers, clog durations exhibit a broad distribution, which poses a challenge for devising efficient unclogging protocols. Using numerical simulations, we demonstrate that the dynamics of arch shapes preceding failure can be modeled as a continuous time random walk (CTRW) with a broad distribution of waiting times, which breaks ergodicity. Treating arch failure as a first passage process of this random walk, we argue that the distribution of unclogging times is determined by this waiting time distribution. We hypothesize that this is a generic feature of unclogging, and that specific characteristics, such as hopper geometry, and mechanical properties of the grains modify the waiting time distribution.


title: A free boundary problem with facets

abstract: We study a free boundary problem on the lattice whose scaling limit is a harmonic free boundary problem with a discontinuous Hamiltonian. We find an explicit formula for the Hamiltonian, prove the solutions are unique, and prove that the limiting free boundary has a facets in every rational direction. Our choice of problem presents difficulties that require the development of a new uniqueness proof for certain free boundary problems. The problem is motivated by physical experiments involving liquid drops on patterned solid surfaces.


title: A note on maximal plurifinely plurisubharmonic functions

abstract: In this note we study the plurifinely locally maximal plurifinely plurisubharmonic functions and improve some known results on these functions. We prove in particular that any locally bounded plurifinely locally maximal plurifinely plurisubharmonic function is maximal, and that any finite plurisubharmonic function on an open subset of $\mathbb{C}^n$ is maximal if and only if it is plurifinely locally maximal.


title: Dynamic Influence Networks for Rule-based Models

abstract: We introduce the Dynamic Influence Network (DIN), a novel visual analytics technique for representing and analyzing rule-based models of protein-protein interaction networks. Rule-based modeling has proved instrumental in developing biological models that are concise, comprehensible, easily extensible, and that mitigate the combinatorial complexity of multi-state and multi-component biological molecules. Our technique visualizes the dynamics of these rules as they evolve over time. Using the data produced by KaSim, an open source stochastic simulator of rule-based models written in the Kappa language, DINs provide a node-link diagram that represents the influence that each rule has on the other rules. That is, rather than representing individual biological components or types, we instead represent the rules about them (as nodes) and the current influence of these rules (as links). Using our interactive DIN-Viz software tool, researchers are able to query this dynamic network to find meaningful patterns about biological processes, and to identify salient aspects of complex rule-based models. To evaluate the effectiveness of our approach, we investigate a simulation of a circadian clock model that illustrates the oscillatory behavior of the KaiC protein phosphorylation cycle.


title: Deep Reinforcement Learning for Resource Allocation in V2V Communications

abstract: In this article, we develop a decentralized resource allocation mechanism for vehicle-to-vehicle (V2V) communication systems based on deep reinforcement learning. Each V2V link is considered as an agent, making its own decisions to find optimal sub-band and power level for transmission. Since the proposed method is decentralized, the global information is not required for each agent to make its decisions, hence the transmission overhead is small. From the simulation results, each agent can learn how to satisfy the V2V constraints while minimizing the interference to vehicle-to-infrastructure (V2I) communications.


title: Type Ia supernovae have two physical width-luminosity relations and they favor sub-Chandrasekhar and direct collision models. I. Bolometric

abstract: While the width-luminosity relation (WLR) among type Ia supernovae (slower is brighter) has been extensively studied, its physical basis has not been convincingly identified. In particular, the 'width' has not been quantitatively linked yet to a physical time scale. We demonstrate that there are two robust fundamental time scales that 1. can be calculated based on integral quantities of the ejecta, with little dependence on radiation transfer modeling and 2. can be inferred from observations. The first is the gamma-ray escape time $t_0$, which determines the long-term evolution of the bolometric light curve and is studied in this Paper I. The second is the recombination time of $^{56}$Fe and $^{56}$Co, which sets the long-term color evolution of the emitted light and is studied in Paper II. Here we show that the gamma-ray escape time $t_0$ can be derived with $\sim 15\%$ accuracy from bolometric observations based on first principles. When applied to a sample of supernovae, the observed values of $t_0$ span a narrow range of $30-45$ days for the wide range of observed $^{56}\rm Ni$ masses $0.1M_{\odot}\lesssim M_{^{56}\rm Ni}\lesssim 1 M_{\odot}$. %This (trivial) bolometric WLR This narrow range of the gamma-ray escape time across the range of luminosities (a trivial WLR) is consistent with central detonations and direct collisions of sub-Chandrasekhar mass white dwarfs (WDs) but not with delayed detonation models for explosions of Chandrasekhar mass WDs, which are therefore disfavored as the primary channel for the population of type Ia supernovae. Computer codes for extracting $t_0$ from observations and models and for calculating gamma-ray transfer in 1D-3D are provided.


title: A Classification-Based Study of Covariate Shift in GAN Distributions

abstract: A basic, and still largely unanswered, question in the context of Generative Adversarial Networks (GANs) is whether they are truly able to capture all the fundamental characteristics of the distributions they are trained on. In particular, evaluating the diversity of GAN distributions is challenging and existing methods provide only a partial understanding of this issue. In this paper, we develop quantitative and scalable tools for assessing the diversity of GAN distributions. Specifically, we take a classification-based perspective and view loss of diversity as a form of covariate shift introduced by GANs. We examine two specific forms of such shift: mode collapse and boundary distortion. In contrast to prior work, our methods need only minimal human supervision and can be readily applied to state-of-the-art GANs on large, canonical datasets. Examining popular GANs using our tools indicates that these GANs have significant problems in reproducing the more distributional properties of their training dataset.


title: Laser electron acceleration on curved surfaces

abstract: Electron acceleration by relativistically intense laser beam propagating along a curved surface allows to split softly the accelerated electron bunch and the laser beam. The presence of a curved surface allows to switch an adiabatic invariant of electrons in the wave instantly leaving the gained energy to the particles. The efficient acceleration is provided by the presence of strong transient quasistationary fields in the interaction region and a long efficient acceleration length. The curvature of the surface allows to select the accelerated particles and provides their narrow angular distribution. The mechanism at work is explicitly demonstrated in theoretical models and experiments.


title: The Achievement of Higher Flexibility in Multiple Choice-based Tests Using Image Classification Techniques

abstract: In spite of the high accuracy of the existing optical mark reading (OMR) systems and devices, a few restrictions remain existent. In this work, we aim to reduce the restrictions of multiple choice questions (MCQ) within tests. We use an image registration technique to extract the answer boxes from answer sheets. Unlike other systems that rely on simple image processing steps to recognize the extracted answer boxes, we address the problem from another perspective by training a machine learning classifier to recognize the class of each answer box (i.e., confirmed, crossed out, or blank answer). This gives us the ability to deal with a variety of shading and mark patterns, and distinguish between chosen (i.e., confirmed) and canceled answers (i.e., crossed out). All existing machine learning techniques require a large number of examples in order to train a model for classification, therefore we present a dataset including six real MCQ assessments with different answer sheet templates. We evaluate two strategies of classification: a straight-forward approach and a two-stage classifier approach. We test two handcrafted feature methods and a convolutional neural network. In the end, we present an easy-to-use graphical user interface of the proposed system. Compared with existing OMR systems, the proposed system has the least constraints and achieves a high accuracy. We believe that the presented work will further direct the development of OMR systems towards reducing the restrictions of the MCQ tests.


title: Energy-Aware Virtual Machine Management in Inter-datacenter Networks over Elastic Optical Infrastructure

abstract: Datacenters (DCs), deployed in a large scale to support the ever increasing demand for data processing applications, consume tremendous energy. Powering DCs with renewable energy can effectively reduce the brown energy consumption. Owing to geographically distributed deployment of DCs, the renewable energy generation and the data processing demands usually vary in different DCs. Migrating virtual machines (VMs) among DCs according to the availability of renewable energy helps match the energy demands and the renewable energy generation in DCs, and thus maximizes the utilization of renewable energy. We first elicit the renewable energy-aware inter-datacenter (inter-DC) VM migration problem in an inter-DC network over the elastic optical infrastructure, present it as a many-manycast communications problem, and then formulate it as an integer linear programming problem. The objective is to minimize the total cost of the brown energy consumption of DCs in such inter-DC network via VM migration. We use CVX and Gurobi to solve this problem for small network configurations, and we propose a few heuristic algorithms that approximate the optimal solution for large network configurations. Through extensive simulations, we show that the proposed algorithms, by migrating VM among DCs, can reduce up to 19.7% cost of the brown energy consumption.


title: An inverse free electron laser acceleration-driven Compton scattering X-ray source

abstract: The generation of X-rays and {\gamma}-rays based on synchrotron radiation from free electrons, emitted in magnet arrays such as undulators, forms the basis of much of modern X-ray science. This approach has the drawback of requiring very high energy, up to the multi-GeV-scale, electron beams, to obtain the required photon energy. Due to the limit in accelerating gradients in conventional particle accelerators, reaching high energy typically demands use of instruments exceeding 100's of meters in length. Compact, less costly, monochromatic X-ray sources based on very high field acceleration and very short period undulators, however, may revolutionize diverse advanced X-ray applications ranging from novel X-ray therapy techniques to active interrogation of sensitive materials, by making them accessible in cost and size. Such compactness may be obtained by an all-optical approach, which employs a laser-driven high gradient accelerator based on inverse free electron laser (IFEL), followed by a collision point for inverse Compton scattering (ICS), a scheme where a laser is used to provide undulator fields. We present an experimental proof-of-principle of this approach, where a TW-class CO2 laser pulse is split in two, with half used to accelerate a high quality electron beam up to 84 MeV through the IFEL interaction, and the other half acts as an electromagnetic undulator to generate up to 13 keV X-rays via ICS. These results demonstrate the feasibility of this scheme, which can be joined with other techniques such as laser recirculation to yield very compact, high brilliance photon sources, extending from the keV to MeV scale. Furthermore, use of the IFEL acceleration with the ICS interaction produces a train of very high intensity X-ray pulses, thus also permitting a unique tool that can be phase-locked to a laser pulse in frontier pump-probe experimental scenarios.


title: Scalable Streaming Tools for Analyzing $N$-body Simulations: Finding Halos and Investigating Excursion Sets in One Pass

abstract: Cosmological $N$-body simulations play a vital role in studying models for the evolution of the Universe. To compare to observations and make a scientific inference, statistic analysis on large simulation datasets, e.g., finding halos, obtaining multi-point correlation functions, is crucial. However, traditional in-memory methods for these tasks do not scale to the datasets that are forbiddingly large in modern simulations. Our prior paper proposes memory-efficient streaming algorithms that can find the largest halos in a simulation with up to $10^9$ particles on a small server or desktop. However, this approach fails when directly scaling to larger datasets. This paper presents a robust streaming tool that leverages state-of-the-art techniques on GPU boosting, sampling, and parallel I/O, to significantly improve performance and scalability. Our rigorous analysis of the sketch parameters improves the previous results from finding the centers of the $10^3$ largest halos to $\sim 10^4-10^5$, and reveals the trade-offs between memory, running time and number of halos. Our experiments show that our tool can scale to datasets with up to $\sim 10^{12}$ particles while using less than an hour of running time on a single GPU Nvidia GTX 1080.


title: The global existence of solutions and asymptotic stability of a reaction-diffusion system

abstract: This paper studies the solutions of a reaction--diffusion system with nonlinearities that generalise the Lengyel--Epstein and FitzHugh--Nagumo nonlinearities. Sufficient conditions are derived for the global asymptotic stability of the system's solutions and confirmed through numerical Examples.


title: Global regolith thermophysical properties of the Moon from the Diviner Lunar Radiometer Experiment

abstract: We used infrared data from the Lunar Reconnaissance Orbiter (LRO) Diviner Lunar Radiometer Experiment to globally map thermophysical properties of the Moon's regolith fines layer. Thermal conductivity varies from 7.4$\times$10$^{-4}$ W m$^{-1}$ K$^{-1}$ at the surface, to 3.4$\times$10$^{-3}$ W m$^{-1}$ K$^{-1}$ at depths of $\sim$1 m, given density values of 1100 kg m$^{-3}$ at the surface, to 1800 kg m$^{-3}$ at 1-m depth. On average, the scale height of these profiles is $\sim$7 cm, corresponding to a thermal inertia of 55 $\pm$2 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$ at 273 K, relevant to the diurnally active near-surface layer, $\sim$4-7 cm. The temperature-dependence of thermal conductivity and heat capacity leads to a $\sim$2$\times$ diurnal variation in thermal inertia at the equator. On global scales, the regolith fines are remarkably uniform, implying rapid homogenization by impact gardening of this layer on timescales $<$ 1 Gyr. Regional and local scale variations show prominent impact features $<$ 1 Gyr old, including higher thermal inertia ($>$ 100 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) in the interiors and ejecta of Copernican-aged impact craters, and lower thermal inertia ($<$ 50 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) within the lunar cold spots identified by Bandfield et al. (2014). Observed trends in ejecta thermal inertia provide a potential tool for age-dating craters of previously unknown age, complementary to the approach suggested by Ghent et al. (2014). Several anomalous regions are identified in the global 128 pixels-per-degree maps presented here, including a high-thermal inertia deposit near the antipode of Tycho crater.


title: Notes on nonlocal dispersal equations in a periodic habitat

abstract: In this paper, we prove that the solution maps of a large class of nonlocal dispersal equations are $\alpha$-contractions, where $\alpha$ is the Kuratowski measure of noncompactness. Then we give some remarks on the spreading speeds and traveling waves for such evolution equations in a periodic habitat.


title: Using reweighting and free energy surface interpolation to predict solid-solid phase diagrams

abstract: Many physical properties of small organic molecules are dependent on the current crystal packing, or polymorph, of the material, including bioavailability of pharmaceuticals, optical properties of dyes, and charge transport properties of semiconductors. Predicting the most stable crystalline form requires determining the crystalline form with the lowest relative Gibbs free energy. Effective computational prediction of the most stable polymorph could save significant time and effort in the design of novel molecular crystalline solids or predict their behavior under new conditions. In this study, we introduce a new approach using multistate reweighting to address the problem of determining solid-solid phase diagrams, and apply this approach to the phase diagram of solid benzene. For this approach, we perform sampling at a selection of temperature and pressure states in the region of interest. We use multistate reweighting methods to determine the reduced free energy differences between $T$ and $P$ states within a given polymorph. The relative stability of the polymorphs at the sampled states can be successively interpolated from these points to create the phase diagram by combining these reduced free energy differences with a reference Gibbs free energy difference between polymorphs. The method also allows for straightforward estimation of uncertainties in the phase boundary. We also find that when properly implemented, multistate reweighting for phase diagram determination scales better with size of system than previously estimated.


title: A representation theorem for the $p^n$ torsion of the Brauer group in characteristic $p$

abstract: If $K$ is a field of characteristic $p$ then the $p$-torsion of the Brauer group, ${}_p{\rm Br\,}(K)$, is represented by a quotient of the group of $1$-forms, $\Omega^1(K)$. Namely, we have a group isomorphism $$\alpha_p:\Omega^1(K)/\langle{\rm d}a,\, (a^p-a){\rm dlog}b\, :\, a,b\in K,\, b\neq 0\rangle\to{}_p{\rm Br\,}(K),$$ given by $a{\rm d}b\mapsto [ab,b)_p$ $\forall a,b\in K$, $b\neq 0$. Here $[\cdot,\cdot )_p:K/\wp (K)\times K^\times/K^{\times p}\to{}_p{\rm Br\,}(K)$ denotes the Artin-Schreier symbol. In this paper we generalize this result. Namely, we prove that for every $n\geq 1$ we have a representation of ${}_{p^n}{\rm Br\,}(K)$ by a quotient of $\Omega^1(W_n(K))$, where $W_n(K)$ is the truncation of length $n$ of the ring of $p$-typical Witt vectors, i.e. $W_{\{1,p,\ldots,p^{n-1}\}}(K)$. Explicitly, we have a group isomorphism $$\alpha_{p^n}:\Omega^1(W_{p^n}(K))/\langle Fa{\rm d}b-a{\rm d}Vb\, :\, a,b\in W_n(K),\, ([a^p]-[a]){\rm dlog}[b]\, :\, a,b\in K,\, b\neq 0\rangle\to{}_{p^n}{\rm Br\,}(K).$$ Here $F$ is the Frobenius isomorphism, $V$ is the Verschiebung map and $[a]$ is the Teichm\"uller representative of $a\in K$, $[a]=(a,0,0,\ldots )$.


title: The Impacts of Three Flamelet Burning Regimes in Nonlinear Combustion Dynamics

abstract: Axisymmetric simulations of a liquid rocket engine are performed using a delayed detached-eddy-simulation (DDES) turbulence model with the Compressible Flamelet Progress Variable (CFPV) combustion model. Three different pressure instability domains are simulated: completely unstable, semi-stable, and fully stable. The different instability domains are found by varying the combustion chamber and oxidizer post length. Laminar flamelet solutions with a detailed chemical mechanism are examined. The $\beta$ Probability Density Function (PDF) for the mixture fraction and Dirac $\delta$ PDF for both the pressure and the progress variable are used. A coupling mechanism between the Heat Release Rate (HRR) and the pressure in an unstable cycle is demonstrated. Local extinction and reignition is investigated for all the instability domains using the full S-curve approach. A monotonic decrease in the amount of local extinctions and reignitions occurs when pressure oscillation amplitude becomes smaller. The flame index is used to distinguish between the premixed and non-premixed burning mode in different stability domains. An additional simulation of the unstable pressure oscillation case using only the stable flamelet burning branch of the S-curve is performed. Better agreement with experiments in terms of pressure oscillation amplitude is found when the full S-curve is used.


title: From which world is your graph?

abstract: Discovering statistical structure from links is a fundamental problem in the analysis of social networks. Choosing a misspecified model, or equivalently, an incorrect inference algorithm will result in an invalid analysis or even falsely uncover patterns that are in fact artifacts of the model. This work focuses on unifying two of the most widely used link-formation models: the stochastic blockmodel (SBM) and the small world (or latent space) model (SWM). Integrating techniques from kernel learning, spectral graph theory, and nonlinear dimensionality reduction, we develop the first statistically sound polynomial-time algorithm to discover latent patterns in sparse graphs for both models. When the network comes from an SBM, the algorithm outputs a block structure. When it is from an SWM, the algorithm outputs estimates of each node's latent position.


title: Compressible Ferrimagnetism in the depleted Periodic Anderson Model

abstract: Tight-binding Hamiltonians with single and multiple orbitals exhibit an intriguing array of magnetic phase transitions. In most cases the spin ordered phases are insulating, while the disordered phases may be either metallic or insulating. In this paper we report a Determinant Quantum Monte Carlo study of interacting electrons in a geometry which can be regarded as a two-dimensional Periodic Anderson Model with depleted interacting ($f$) orbitals. For a single depletion, we observe an enhancement of antiferromagnetic correlations and formation of localized states. For half of the $f$-orbitals regularly depleted, the system exhibits a ferrimagnetic ground state. We obtain a quantitative determination of the nature of magnetic order, which we discuss in the context of Tsunetsugu's theorem, and show that, although the dc conductivity indicates insulating behavior at half-filling, the compressibility remains finite.


title: Fast integration of DPG matrices based on tensorization

abstract: Numerical integration of the stiffness matrix in higher order finite element (FE) methods is recognized as one of the heaviest computational tasks in a FE solver. The problem becomes even more relevant when computing the Gram matrix in the algorithm of the Discontinuous Petrov Galerkin (DPG) FE methodology. Making use of 3D tensor-product shape functions, and the concept of sum factorization, known from standard high order FE and spectral methods, here we take advantage of this idea for the entire exact sequence of FE spaces defined on the hexahedron. The key piece to the presented algorithms is the exact sequence for the one-dimensional element, and use of hierarchical shape functions. Consistent with existing results, the presented algorithms for the integration of $H^1$, $H(\text{curl})$, $H(\text{div})$, and $L^2$ inner products, have the $O(p^7)$ computational complexity. Additionally, a modified version of the algorithms is proposed when the element map can be simplified, resulting in the reduced $O(p^6)$ complexity. Use of Legendre polynomials for shape functions is critical in this implementation. Computational experiments performed with $H^1$, $H(\text{div})$ and $H(\text{curl})$ test shape functions show good correspondence with the expected rates.


title: Modular Virasoro Vertex Algebras and Affine Vertex Algebras

abstract: In this paper, we study Virasoro vertex algebras and affine vertex algebras over a general field of characteristic $p>2$. More specifically, we study certain quotients of the universal Virasoro and affine vertex algebras by ideals related to the $p$-centers of the Virasoro algebra and affine Lie algebras. Among the main results, we classify their irreducible $\mathbb{N}$-graded modules by explicitly determining their Zhu algebras and show that these vertex algebras have only finitely many irreducible $\mathbb{N}$-graded modules and they are $C_2$-cofinite.


title: Symmetric Invariant Bilinear Forms on Modular Vertex Algebras

abstract: In this paper, we study contragredient duals and invariant bilinear forms for modular vertex algebras (in characteristic $p$). We first introduce a bialgebra $\mathcal{H}$ and we then introduce a notion of $\mathcal{H}$-module vertex algebra and a notion of $(V,\mathcal{H})$-module for an $\mathcal{H}$-module vertex algebra $V$. Then we give a modular version of Frenkel-Huang-Lepowsky's theory and study invariant bilinear forms on an $\mathcal{H}$-module vertex algebra. As the main results, we obtain an explicit description of the space of invariant bilinear forms on a general $\mathcal{H}$-module vertex algebra, and we apply our results to affine vertex algebras and Virasoro vertex algebras.


title: Analysis of Biased Stochastic Gradient Descent Using Sequential Semidefinite Programs

abstract: We present a convergence rate analysis for biased stochastic gradient descent (SGD), where individual gradient updates are corrupted by computation errors. We develop stochastic quadratic constraints to formulate a small linear matrix inequality (LMI) whose feasible points lead to convergence bounds of biased SGD. Based on this LMI condition, we develop a sequential minimization approach to analyze the intricate trade-offs that couple stepsize selection, convergence rate, optimization accuracy, and robustness to gradient inaccuracy. We also provide feasible points for this LMI and obtain theoretical formulas that quantify the convergence properties of biased SGD under various assumptions on the loss functions.


title: Electrical-driven Plasmon Source on Silicon based on Quantum Tunneling

abstract: An efficient silicon based light source presents an unreached goal in the field of photonics, due to Silicons indirect electronic band structure preventing direct carrier recombination and subsequent photon emission. Here we utilize inelastically tunneling electrons to demonstrate an electrically driven light emitting silicon based tunnel junction operating at room temperature. We show that such a junction is a source for plasmons driven by the electrical tunnel current. We find that the emission spectrum is not given by the quantum condition where the emission frequency would be proportional to the applied voltage, but the spectrum is determined by the spectral overlap between the energy dependent tunnel current and the modal dispersion of the plasmon. Experimentally we find the highest light outcoupling efficiency corresponding to the skin depth of the metallic contact of this metal insulator semiconductor junction. Distinct from LEDs, the temporal response of this tunnel source is not governed by nanosecond carrier lifetimes known to semiconductors, but rather by the tunnel event itself and Heisenbergs uncertainty principle.


title: The cardiac Na$^+$/K$^+$ ATPase: An updated, thermodynamically consistent model

abstract: The Na$^+$/K$^+$ ATPase is an essential component of cardiac electrophysiology, maintaining physiological Na$^+$ and K$^+$ concentrations over successive heart beats. Terkildsen et al. (2007) developed a model of the ventricular myocyte Na$^+$/K$^+$ ATPase to study extracellular potassium accumulation during ischaemia, demonstrating the ability to recapitulate a wide range of experimental data, but unfortunately there was no archived code associated with the original manuscript. Here we detail an updated version of the model and provide CellML and MATLAB code to ensure reproducibility and reusability. We note some errors within the original formulation which have been corrected to ensure that the model is thermodynamically consistent, and although this required some reparameterisation, the resulting model still provides a good fit to experimental measurements that demonstrate the dependence of Na$^+$/K$^+$ ATPase pumping rate upon membrane voltage and metabolite concentrations. To demonstrate thermodynamic consistency we also developed a bond graph version of the model. We hope that these models will be useful for community efforts to assemble a whole-cell cardiomyocyte model which facilitates the investigation of cellular energetics.


title: Cogenesis of LIGO Primordial Black Holes and Dark Matter

abstract: In this letter, we propose a novel scenario which simultaneously explains $\mathcal{O}(10)M_\odot$ primordial black holes (PBHs) and dark matter in the minimally supersymmetric standard model. Gravitational waves (GWs) events detected by LIGO-Virgo collaboration suggest an existence of black holes as heavy as $\sim 30M_\odot$. In our scenario, as seeds of the PBHs, we make use of the baryon number perturbations which are induced by the special type of Affleck-Dine mechanism. Furthermore, the scenario does not suffer from the stringent constraints from CMB $\mu$-distortion due to the Silk damping and pulsar timing. We find the scenario can explain not only the current GWs events consistently, but also dark matter abundance by the non-topological solitons formed after Affleck-Dine mechanism, called Q-balls.


title: The neighborhood lattice for encoding partial correlations in a Hilbert space

abstract: Neighborhood regression has been a successful approach in graphical and structural equation modeling, with applications to learning undirected and directed graphical models. We extend these ideas by defining and studying an algebraic structure called the neighborhood lattice based on a generalized notion of neighborhood regression. We show that this algebraic structure has the potential to provide an economic encoding of all conditional independence statements in a Gaussian distribution (or conditional uncorrelatedness in general), even in the cases where no graphical model exists that could "perfectly" encode all such statements. We study the computational complexity of computing these structures and show that under a sparsity assumption, they can be computed in polynomial time, even in the absence of the assumption of perfectness to a graph. On the other hand, assuming perfectness, we show how these neighborhood lattices may be "graphically" computed using the separation properties of the so-called partial correlation graph. We also draw connections with directed acyclic graphical models and Bayesian networks. We derive these results using an abstract generalization of partial uncorrelatedness, called partial orthogonality, which allows us to use algebraic properties of projection operators on Hilbert spaces to significantly simplify and extend existing ideas and arguments. Consequently, our results apply to a wide range of random objects and data structures, such as random vectors, data matrices, and functions.


title: Shelstad's character identity from the point of view of index theory

abstract: Shelstad's character identity is an equality between sums of characters of tempered representations in corresponding $L$-packets of two real, semisimple, linear, algebraic groups that are inner forms to each other. We reconstruct this character identity in the case of the discrete series, using index theory of elliptic operators in the framework of $K$-theory. Our geometric proof of the character identity is evidence that index theory can play a role in the classification of group representations via the Langlands program.


title: Tau neutrinos from ultracompact dark matter minihalos and constraints on the primordial curvature perturbations

abstract: The observations and research on the neutrinos provide a kind of indirect way of revealing the properties of dark matter particles. For the detection of muon neutrinos, the main issue is the large atmospheric background, which is caused by the interactions between the cosmic rays and atoms within the atmosphere. Compared with muon neutrinos, tau neutrinos have a smaller atmospheric background especially for the downward-going direction. Except for the classical neutrino sources, dark matter particles can also annihilate into the neutrinos and are the potential high energy astrophysical sources. The annihilation rate of dark matter particles is proportional to the square of number density; therefore, the annihilation rate is large near the center of dark matter halos especially for the new kind of dark matter structures named ultracompact dark matter minihalos (UCMHs). In previous works, we have investigated the potential muon neutrino flux from UCMHs due to dark matter annihilation. Moreover, since the formation of UCMHs is related to the primordial density perturbations of small scales, we get the constraints on the amplitude of the primordial curvature perturbations of small scales, $1 \lesssim k \lesssim 10^{7} ~\rm Mpc^{-1}$. In this work, we focus on the downward-going tau neutrinos from UCMHs due to dark matter annihilation. Compared with the background of tau neutrino flux we get the constraints on the mass fraction of UCMHs. Then using the limits on the mass fraction of UCMHs we got the constraints on the amplitude of the primordial curvature perturbations which are extended to the scale $k \sim 10^{8} ~ \rm Mpc^{-1}$ compared with previous results.


title: Remark on upper bound for lifespan of solutions to semilinear evolution equations in a two-dimensional exterior domain

abstract: In this paper we consider the initial-boundary value problem for the heat, damped wave, complex-Ginzburg-Landau and Schr"odinger equations with the power type nonlinearity $|u|^p$ with $p in (1,2]$ in a two-dimensional exterior domain. Remark that $2=1+2/N$ is well-known as the Fujita exponent. If $p>2$, then there exists a small global-in-time solution of the damped wave equation for some initial data small enough (see Ikehata'05), and if $p<2$, then global-in-time solutions cannot exist for any positive initial data (see Ogawa-Takeda'09 and Lai-Yin'17). The result is that for given initial data $(f,tau g)in H^1_0(Omega)times L^2(Omega)$ satisfying $(f+tau g)log |x|in L^1(Omega)$ with some requirement, the solution blows up at finite time, and moreover, the upper bound for lifespan of solutions to the problem is given as the following {it double exponential type} when $p=2$: [ lifespan(u) leq exp[exp(Cep^{-1})] . ] The crucial idea is to use test functions which approximates the harmonic function $log |x|$ satisfying Dirichlet boundary condition and the technique for derivation of lifespan estimate in Ikeda-Sobajima(arXiv:1710.06780).


title: Blowing Polar Skyrmion Bubbles in Oxide Superlattices

abstract: Particle-like topological structures such as skyrmions and vortices have garnered ever-increasing interests due to the rich physical insights and potential broad applications. Here we discover the reversible switching between polar skyrmion bubbles and ordered vortex arrays in ferroelectric superlattices under an electric field, reminiscent of the Plateau-Raleigh instability in fluid mechanics. Electric field phase diagram is constructed, showing wide stability window for the observed polar skyrmions. This study is a demonstration for the computational design of ferroelectric topological structures and field-induced topological phase transitions.


title: Black hole and cosmos with multiple horizons and multiple singularities in vector-tensor theories

abstract: A stationary and spherically symmetric black hole (For example, Reissner-Nordstrom black hole or Kerr-Newman black hole) has at most one singularity and two horizons. One horizon is the outer event horizon and the other is the inner Cauchy horizon. Can we construct static and spherically symmetric black hole solutions with N horizons and M singularities? De Sitter cosmos has only one apparent horizon. Can we construct cosmos solutions with N horizons? In this article, we present the static and spherically symmetric black hole and cosmos solutions with N horizons and M singularities in the vector-tensor theories. Following these motivations, we also construct the black hole solutions with a firewall. The deviation of these black hole solutions from the usual ones can be potentially tested by future measurements of gravitational waves.


title: A Software-Based Lock-In Measurement for Student Laboratories

abstract: A student laboratory experiment is presented that introduces the concept of a lock-in measurement through the exploration of the relationship between the power detected from a modulated light source and the distance between the source and the detector. The lock-in measurement is done in software using time streams for both the reference and the detector signals. The achievable experimental sensitivity is shown to increase several orders of magnitude over the simple case without using modulation. In addition, the student becomes familiar with computational techniques pertinent to large data sets.


title: Gr\"unbaum's inequality for sections

abstract: We show \begin{align*} \frac{ \int_{E \cap \theta^+} f(x) dx }{ \int_E f(x) dx } \geq \left(\frac{k \gamma+1}{(n+1) \gamma+1}\right)^{\frac{k \gamma+1}{\gamma}} \end{align*} for all $k$-dimensional subspaces $E\subset\mathbb{R}^n$, $\theta\in E\cap S^{n-1}$, and all $\gamma$-concave functions $f:\mathbb{R}^n\rightarrow [0,\infty)$ with $\gamma >0$, $0< \int_{\mathbb{R}^n} f(x)\, dx <\infty$, and $\int_{\mathbb{R}^n} x f(x)\, dx$ at the origin $o\in\mathbb{R}^n$. Here, $\theta^+ := \lbrace x\, : \, \langle x,\theta\rangle \geq 0 \rbrace$. As a consequence of this result, we get the following generalization of Gr\"unbaum's inequality: \begin{align*} \frac{ \mbox{vol}_k(K\cap E\cap\theta^+) }{ \mbox{vol}_k(K\cap E) } \geq \left( \frac{k}{n+1} \right)^k \end{align*} for all convex bodies $K\subset\mathbb{R}^n$ with centroid at the origin, $k$-dimensional subspaces $E\subset\mathbb{R}^n$, and $\theta\in E\cap S^{n-1}$. The lower bounds in both of our inequalities are the best possible, and we discuss the equality conditions.


title: Thick Bi2Sr2CaCu2O8+{\delta} films grown by liquid-phase epitaxy for Josephson THz applications

abstract: Theoretical and experimental studies of intrinsic Josephson junctions that naturally occur in high-Tc superconducting Bi2Sr2CaCu2O8+{\delta} (Bi-2212) have demonstrated their potential for novel types of compact devices for the generation and sensing of electromagnetic radiation in the THz range. Here, we show that the THz-on-a-chip concept may be realized in liquid phase epitaxial-grown (LPE) thick Bi-2212 films. We have grown {\mu}m-thick Bi-2212 LPE films on MgO substrates. These films display excellent c-axis alignment and single crystal grains of about 650x150 {\mu}m2 in size. A branched current-voltage characteristic was clearly observed in c-axis transport, which is a clear signature of underdamped intrinsic Josephson junctions, and a prerequisite for THz-generation. We discuss LPE growth conditions allowing improvement of the structural quality and superconducting properties of Bi-2212 films for THz applications.


title: Switchable polar spirals in tricolor oxide superlattices

abstract: There are increasing evidences that ferroelectric states at the nanoscale can exhibit fascinating topological structures including polar vortices and skyrmions,akin to those observed in the ferromagnetic systems. Here we report the discovery of a new type of polar topological structure,ordered array of nanoscale spirals in the tricolor ferroelectric superlattice system via phase field simulations. This polar spiral structure is composed of fine ordered semivortex arrays with vortex cores forming a wavy distribution. It is demonstrated that this tricolor system has an ultrahigh Curie temperature of 1000 K and a temperature of 650 K for the phase transformation from spiral structure to inplane orthorhombic domain structure, showing a greatly enhanced thermal stability than the recently discovered polar vortex lattices in the superlattice system. Moreover, the spiral structure has a net inplane polarization that could be switched by an experimentally feasible irrotational inplane field. The switching process involves a metastable vortex state, and is fully reversible.


title: Thermodynamics of Spin-1/2 Kagom\'e Heisenberg Antiferromagnet: Algebraic Paramagnetic Liquid and Finite-Temperature Phase Diagram

abstract: Quantum fluctuations from frustration can trigger quantum spin liquids (QSLs) at zero temperature. However, it is unclear how thermal fluctuations affect a QSL. We employ state-of-the-art tensor network-based methods to explore the ground state and thermodynamic properties of the spin-1/2 kagome Heisenberg antiferromagnet (KHA). Its ground state is shown to be consistent with a gapless QSL by observing the absence of zero-magnetization plateau as well as the algebraic behaviors of susceptibility and specific heat at low temperatures, respectively. We show that there exists an \textit{algebraic paramagnetic liquid} (APL) that possesses both the paramagnetic properties and the algebraic behaviors inherited from the QSL. The APL is induced under the interplay between quantum fluctuations from geometrical frustration and thermal fluctuations. By studying the temperature-dependent behaviors of specific heat and magnetic susceptibility, a finite-temperature phase diagram in a magnetic field is suggested, where various phases are identified. This present study gains useful insight into the thermodynamic properties of the spin-1/2 KHA with or without a magnetic field and is helpful for relevant experimental studies.


title: Infrared Emission from Supernova Remnants: Formation and Destruction of Dust

abstract: We review the observations of dust emission in supernova rem- nants (SNRs) and supernovae (SNe). Theoretical calculations suggest that SNe, particularly core-collapse, should make significant quantities of dust, perhaps as much as a solar mass. Observations of extragalactic SNe have yet to find anywhere near this amount, but this may be the result of observa- tional limitations. SN 1987A, in the process of transitioning from a SN to an SNR, does show signs of a significant amount of dust forming in its ejecta, but whether this dust will survive the passage of the reverse shock to be injected into the ISM is unknown. IR observations of SNRs have not turned up significant quantities of dust, and the dust that is observed is generally swept-up by the forward shock, rather than created in the ejecta. Because the shock waves also destroy dust in the ISM, we explore the question of whether SNe might be net destroyers, rather than net creators of dust in the universe.


abstract: A minor of a matrix is quasi-principal if it is a principal or an almost-principal minor. The quasi principal rank characteristic sequence (qpr-sequence) of an $n\times n$ symmetric matrix is introduced, which is defined as $q_1 q_2 \cdots q_n$, where $q_k$ is $\tt A$, $\tt S$, or $\tt N$, according as all, some but not all, or none of its quasi-principal minors of order $k$ are nonzero. This sequence extends the principal rank characteristic sequences in the literature, which only depend on the principal minors of the matrix. A necessary condition for the attainability of a qpr-sequence is established. Using probabilistic techniques, a complete characterization of the qpr-sequences that are attainable by symmetric matrices over fields of characteristic $0$ is given.


title: AxonDeepSeg: automatic axon and myelin segmentation from microscopy data using convolutional neural networks

abstract: Segmentation of axon and myelin from microscopy images of the nervous system provides useful quantitative information about the tissue microstructure, such as axon density and myelin thickness. This could be used for instance to document cell morphometry across species, or to validate novel non-invasive quantitative magnetic resonance imaging techniques. Most currently-available segmentation algorithms are based on standard image processing and usually require multiple processing steps and/or parameter tuning by the user to adapt to different modalities. Moreover, only few methods are publicly available. We introduce AxonDeepSeg, an open-source software that performs axon and myelin segmentation of microscopic images using deep learning. AxonDeepSeg features: (i) a convolutional neural network architecture; (ii) an easy training procedure to generate new models based on manually-labelled data and (iii) two ready-to-use models trained from scanning electron microscopy (SEM) and transmission electron microscopy (TEM). Results show high pixel-wise accuracy across various species: 85% on rat SEM, 81% on human SEM, 95% on mice TEM and 84% on macaque TEM. Segmentation of a full rat spinal cord slice is computed and morphological metrics are extracted and compared against the literature. AxonDeepSeg is freely available at https://github.com/neuropoly/axondeepseg


title: In-Bed Pose Estimation: Deep Learning with Shallow Dataset

abstract: Although human pose estimation for various computer vision (CV) applications has been studied extensively in the last few decades, yet in-bed pose estimation using camera-based vision methods has been ignored by the CV community because it is assumed to be identical to the general purpose pose estimation methods. However, in-bed pose estimation has its own specialized aspects and comes with specific challenges including the notable differences in lighting conditions throughout a day and also having different pose distribution from the common human surveillance viewpoint. In this paper, we demonstrate that these challenges significantly lessen the effectiveness of existing general purpose pose estimation models. In order to address the lighting variation challenge, infrared selective (IRS) image acquisition technique is proposed to provide uniform quality data under various lighting conditions. In addition, to deal with unconventional pose perspective, a 2-end histogram of oriented gradient (HOG) rectification method is presented. In this work, we explored the idea of employing a pre-trained convolutional neural network (CNN) model trained on large public datasets of general human poses and fine-tuning the model using our own shallow in-bed IRS dataset. We developed an IRS imaging system and collected IRS image data from several realistic life-size mannequins in a simulated hospital room environment. A pre-trained CNN called convolutional pose machine (CPM) was repurposed for in-bed pose estimation by fine-tuning its specific intermediate layers. Using the HOG rectification method, the pose estimation performance of CPM significantly improved by 26.4% in PCK0.1 criteria compared to the model without such rectification.


title: Towards Neural Machine Translation with Partially Aligned Corpora

abstract: While neural machine translation (NMT) has become the new paradigm, the parameter optimization requires large-scale parallel data which is scarce in many domains and language pairs. In this paper, we address a new translation scenario in which there only exists monolingual corpora and phrase pairs. We propose a new method towards translation with partially aligned sentence pairs which are derived from the phrase pairs and monolingual corpora. To make full use of the partially aligned corpora, we adapt the conventional NMT training method in two aspects. On one hand, different generation strategies are designed for aligned and unaligned target words. On the other hand, a different objective function is designed to model the partially aligned parts. The experiments demonstrate that our method can achieve a relatively good result in such a translation scenario, and tiny bitexts can boost translation quality to a large extent.


title: Wireless Network Simplification: The Performance of Routing

abstract: Consider a wireless Gaussian network where a source wishes to communicate with a destination with the help of N full-duplex relay nodes. Most practical systems today route information from the source to the destination using the best path that connects them. In this paper, we show that routing can in the worst case result in an unbounded gap from the network capacity - or reversely, physical layer cooperation can offer unbounded gains over routing. More specifically, we show that for $N$-relay Gaussian networks with an arbitrary topology, routing can in the worst case guarantee an approximate fraction $\frac{1}{\left\lfloor N/2 \right\rfloor + 1}$ of the capacity of the full network, independently of the SNR regime. We prove that this guarantee is fundamental, i.e., it is the highest worst-case guarantee that we can provide for routing in relay networks. Next, we consider how these guarantees are refined for Gaussian layered relay networks with $L$ layers and $N_L$ relays per layer. We prove that for arbitrary $L$ and $N_L$, there always exists a route in the network that approximately achieves at least $\frac{2}{(L-1)N_L + 4}$ $\left(\mbox{resp.}\frac{2}{LN_L+2}\right)$ of the network capacity for odd $L$ (resp. even $L$), and there exist networks where the best routes exactly achieve these fractions. These results are formulated within the network simplification framework, that asks what fraction of the capacity we can achieve by using a subnetwork (in our case, a single path). A fundamental step in our proof is a simplification result for MIMO antenna selection that may also be of independent interest. To the best of our knowledge, this is the first result that characterizes, for general wireless network topologies, what is the performance of routing with respect to physical layer cooperation techniques that approximately achieve the network capacity.


title: Cost-Efficient and Robust On-Demand Video Transcoding Using Heterogeneous Cloud Services

abstract: Video streams usually have to be transcoded to match the characteristics of viewers' devices. Streaming providers have to store numerous transcoded versions of a given video to serve various display devices. Given the fact that viewers' access pattern to video streams follows a long tail distribution, for the video streams with low access rate, we propose to transcode them in an on-demand manner using cloud computing services. The challenge in utilizing cloud services for on-demand video transcoding is to maintain a robust QoS for viewers and cost-efficiency for streaming service providers. To address this challenge, we present the Cloud-based Video Streaming Services (CVS2) architecture. It includes a QoS-aware scheduling that maps transcoding tasks to the VMs by considering the affinity of the transcoding tasks with the allocated heterogeneous VMs. To maintain robustness in the presence of varying streaming requests, the architecture includes a cost-efficient VM Provisioner. This component provides a self- configurable cluster of heterogeneous VMs. The cluster is reconfigured dynamically to maintain the maximum affinity with the arriving workload. Results obtained under diverse workload conditions demonstrate that CVS2 architecture can maintain a robust QoS for viewers while reducing the incurred cost of the streaming service provider up to 85%


title: Isogeometric B\'ezier dual mortaring: Refineable higher-order spline dual bases and weakly continuous geometry

abstract: In this paper we develop the isogeometric B\'ezier dual mortar method. It is based on B\'ezier extraction and projection and is applicable to any spline space which can be represented in B\'ezier form (i.e., NURBS, T-splines, LR-splines, etc.). The approach weakly enforces the continuity of the solution at patch interfaces and the error can be adaptively controlled by leveraging the refineability of the underlying dual spline basis without introducing any additional degrees of freedom. We also develop weakly continuous geometry as a particular application of isogeometric B\'ezier dual mortaring. Weakly continuous geometry is a geometry description where the weak continuity constraints are built into properly modified B\'ezier extraction operators. As a result, multi-patch models can be processed in a solver directly without having to employ a mortaring solution strategy. We demonstrate the utility of the approach on several challenging benchmark problems. Keywords: Mortar methods, Isogeometric analysis, B\'ezier extraction, B\'ezier projection