title: Duality between Final-Seed and Initial-Seed Mutations in Cluster Algebras

abstract: We study the duality between the mutations and the initial-seed mutations in cluster algebras, where the initial-seed mutations are the transformations of rational expressions of cluster variables in terms of the initial cluster under the change of the initial cluster. In particular, we define the maximal degree matrices of the $F$-polynomials called the $F$-matrices and show that the $F$-matrices have the self-duality which is analogous to the duality between the $C$- and $G$-matrices.


title: Experimental Design Modulates Variance in BOLD Activation: The Variance Design General Linear Model

abstract: Typical fMRI studies have focused on either the mean trend in the blood-oxygen-level-dependent (BOLD) time course or functional connectivity (FC). However, other statistics of the neuroimaging data may contain important information. Despite studies showing links between the variance in the BOLD time series (BV) and age and cognitive performance, a formal framework for testing these effects has not yet been developed. We introduce the Variance Design General Linear Model (VDGLM), a novel framework that facilitates the detection of variance effects. We designed the framework for general use in any fMRI study by modeling both mean and variance in BOLD activation as a function of experimental design. The flexibility of this approach allows the VDGLM to i) simultaneously make inferences about a mean or variance effect while controlling for the other and ii) test for variance effects that could be associated with multiple conditions and/or noise regressors. We demonstrate the use of the VDGLM in a working memory application and show that engagement in a working memory task is associated with whole-brain decreases in BOLD variance.


title: A Novel Optimization Approach to Fictitious Domain Methods

abstract: A new approach to the solution of boundary value problems within the so-called fictitious domain methods philosophy is proposed which avoids well known shortcomings of other fictitious domain methods, including the need to generate extensions of the data. The salient feature of the novel method, which we refer to as SSEM (Smooth Selection Embedding Method), is that it reduces the whole boundary value problem to a linear constraint for an appropriate optimization problem formulated in a larger simpler set containing the domain on which the boundary value problem is posed and which allows for the use of straightforward discretizations. The proposed method in essence computes a (discrete) extension of the solution to the boundary value problem by selecting it as a smooth element of the complete affine family of solutions of the extended, yet unmodified, under-determined problem. The actual regularity of this extension is determined by that of the analytic solution and the choice of obejctive functional. Numerical experiments will demonstrate that it can be stably used to efficiently deal with non-constant coefficients, general geometries, and different boundary conditions in dimensions d=1,2,3 and that it produces solutions of tunable (and high) accuracy.


title: Connections Between Complex PT-Invariant Solutions and Complex Periodic Solutions of Several Nonlinear Equations

abstract: We point out novel connections between complex PT-invariant solutions of several nonlinear equations such as $\phi^4$, $\phi^6$, sine-Gordon, hyperbolic sine-Gordon, double sine-Gordon, double hyperbolic sine-Gordon, mKdV, etc. We then use these connections to obtain several new complex PT-invariant periodic solutions of some of these equations.


title: Representations of simple noncommutative Jordan superalgebras I

abstract: In this article we begin the study of representations of simple finite-dimensional noncommutative Jordan superalgebras. In the case of degree $\geq 3$ we show that any finite-dimensional representation is completely reducible and, depending on the superalgebra, quasiassociative or Jordan. Then we study representations of superalgebras $D_t(\alpha,\beta,\gamma)$ and $K_3(\alpha, \beta, \gamma)$ and prove the Kronecker factorization theorem for superalgebras $D_t(\alpha,\beta,\gamma)$. In the last section we use a new approach to study noncommutative Jordan representations of simple Jordan superalgebras.


title: Thermohaline mixing in extremely metal-poor stars

abstract: Extremely metal-poor (EMP) stars are an integral piece in the puzzle that is the early Universe, and although anomolous subclasses of EMP stars such as carbon-enhanced metal-poor (CEMP) stars are well-studied, they make up less than half of all EMP stars with [Fe/H] $\sim -3.0$. The amount of carbon depletion occurring on the red giant branch (carbon offset) is used to determine the evolutionary status of EMP stars, and this offset will differ between CEMP and normal EMP stars. The depletion mechanism employed in stellar models (from which carbon offfsets are derived) is very important, however the only widely available carbon offsets in the literature are derived from stellar models using a thermohaline mixing mechanism that cannot simultaneously match carbon and lithium abundances to observations for a single diffusion coeffcient. Our stellar evolution models utilise a modified thermohaline mixing model that can match carbon and lithium in the metal-poor globular cluster NGC 6397. We compare our models to the bulk of the EMP star sample at [Fe/H] $= -3$ and show that our modified models follow the trend of the observations and deplete less carbon compared to the standard thermohaline mixing theory. We conclude that stellar models that employ the standard thermohaline mixing formalism overestimate carbon offsets and hence CEMP star frequencies, particularly at metallicities where carbon-normal stars dominate the EMP star population.


title: The Parameterized Complexity of Finding Point Sets with Hereditary Properties

abstract: We consider problems where the input is a set of points in the plane and an integer $k$, and the task is to find a subset $S$ of the input points of size $k$ such that $S$ satisfies some property. We focus on properties that depend only on the order type of the points and are monotone under point removals. We show that not all such problems are fixed-parameter tractable parameterized by $k$, by exhibiting a property defined by three forbidden patterns for which finding a $k$-point subset with the property is $\mathrm{W}[1]$-complete and (assuming the exponential time hypothesis) cannot be solved in time $n^{o(k/\log k)}$. However, we show that problems of this type are fixed-parameter tractable for all properties that include all collinear point sets, properties that exclude at least one convex polygon, and properties defined by a single forbidden pattern.


title: Holographic $J/\psi$ production near threshold and the proton mass problem

abstract: It has been suggested that the production of a heavy quarkonium near threshold in electron-proton scattering can shed light on the origin of the proton mass via the QCD trace anomaly. We study the photoproduction of $J/\psi$ off the proton using gauge/string duality and demonstrate that the $t$-dependence of the differential cross section $d\sigma/dt$ at small-$t$ is a sensitive probe of the trace anomaly.


title: A Note on Transportation Cost Inequalities for Diffusions with Reflections

abstract: We prove that reflected Brownian motion with normal reflections in a convex domain satisfies a dimension free Talagrand type transportation cost-information inequality. The result is generalized to other reflected diffusions with suitable drifts and diffusions. We apply this to get such an inequality for interacting Brownian particles with rank-based drift and diffusion coefficients such as the infinite Atlas model. This is an improvement over earlier dimension-dependent results.


title: Scalar Gravitational Radiation from Binaries: Vainshtein Mechanism in Time-dependent Systems

abstract: We develop a full four-dimensional numerical code to study scalar gravitational radiation emitted from binary systems and probe the Vainshtein mechanism in situations that break the static and spherical symmetry, relevant for binary pulsars as well as black holes and neutron stars binaries. The present study focuses on the cubic Galileon which arises as the decoupling limit of massive theories of gravity. Limitations associated with the numerical methods prevent us from reaching a physically realistic hierarchy of scales; nevertheless, within this context we observe the same power law scaling of the radiated power as previous analytic estimates, and confirm a strong suppression of the power emitted in the monopole and dipole as compared with quadrupole radiation. Following the trend to more physically realistic parameters, we confirm the suppression of the power emitted in scalar gravitational radiation and the recovery of General Relativity with good accuracy. This paves the way for future numerical work, probing more generic, physically relevant situations and sets of interactions that may exhibit the Vainshtein mechanism.


title: Measurement of the phase between strong and electromagnetic amplitudes of $J/\psi$ decays

abstract: Using 16 energy points of $e^{+}e^{-}$ annihilation data collected in the vicinity of the $J/\psi$ resonance with the BESIII detector and with a total integrated luminosity of around 100 pb$^{-1}$, we study the relative phase between the strong and electromagnetic amplitudes of $J/\psi$ decays. The relative phase between $J\psi$ electromagnetic decay and the continuum process ($e^{+}e^{-}$ annihilation without the $J/\psi$ resonance) is confirmed to be zero by studying the cross section lineshape of $\mu^{+}\mu^{-}$ production. The relative phase between $J/\psi$ strong and electromagnetic decays is then measured to be $(84.9\pm3.6)^\circ$ or $(-84.7\pm3.1)^\circ$ for the $2(\pi^{+}\pi^{-})\pi^{0}$ final state by investigating the interference pattern between the $J/\psi$ decay and the continuum process. This is the first measurement of the relative phase between $J/\psi$ strong and electromagnetic decays into a multihadron final state using the lineshape of the production cross section. We also study the production lineshape of the multihadron final state $\eta\pi^{+}\pi^{-}$ with $\eta\to\pi^{+}\pi^{-}\pi^{0}$, which provides additional information about the phase between the $J/\psi$ electromagnetic decay amplitude and the continuum process. Additionally, the branching fraction of $J/\psi\to 2(\pi^{+}\pi^{-})\pi^{0}$ is measured to be $(4.73\pm0.44)\%$ or $(4.85\pm0.45)\%$, and the branching fraction of $J/\psi\to\eta\pi^{+}\pi^{-}$ is measured to be $(3.78\pm0.68)\times10^{-4}$. Both of them are consistent with the world average values. The quoted uncertainties include both statistical and systematic uncertainties, which are mainly caused by the low statistics.


title: Efficient Fusion of Sparse and Complementary Convolutions

abstract: We propose a new method to create compact convolutional neural networks (CNNs) by exploiting sparse convolutions. Different from previous works that learn sparsity in models, we directly employ hand-crafted kernels with regular sparse patterns, which result in the computational gain in practice without sophisticated and dedicated software or hardware. The core of our approach is an efficient network module that linearly combines sparse kernels to yield feature representations as strong as those from regular kernels. We integrate this module into various network architectures and demonstrate its effectiveness on three vision tasks, object classification, localization and detection. For object classification and localization, our approach achieves comparable or better performance than several baselines and related works while providing lower computational costs with fewer parameters (on average, a $2-4\times$ reduction of convolutional parameters and computation). For object detection, our approach leads to a VGG-16-based Faster RCNN detector that is 12.4$\times$ smaller and about 3$\times$ faster than the baseline.


title: Rotationally induced failure of irregularly shaped asteroids

abstract: Many asteroids are rubble piles with irregular shapes. While the irregular shapes of large asteroids may be attributed to collisional events, those of small asteroids may result from not only impact events but also rotationally induced failure, a long-term consequence of small torques caused by, for example, solar radiation pressure. A better understanding of shape deformation induced by such small torques will allow us to give constraints on the evolution process of an asteroid and its structure. However, no quantitative study has been reported to provide the relationship between an asteroid's shape and its failure mode due to its fast rotation. Here, we use a finite element model (FEM) technique to analyze the failure modes and conditions of 24 asteroids with diameters less than 30 - 40 km, which were observed at high resolution by ground radar or asteroid exploration missions. Assuming that the material distribution is uniform, we investigate how these asteroids fail structurally at different spin rates. Our FEM simulations describe the detailed deformation mode of each irregularly shaped asteroid at fast spin. The failed regions depend on the original shape. Spheroidal objects structurally fail from the interior, while elongated objects experience structural failure on planes perpendicular to the minimum moment of inertia axes in the middle of their structure. Contact binary objects have structural failure across their most sensitive cross sections. We further investigate if our FEM analysis is consistent with earlier works that theoretically explored a uniformly rotating triaxial ellipsoid. The results show that global shape variations may significantly change the failure condition of an asteroid. Our work suggests that it is critical to take into account the actual shapes of asteroids to explore their failure modes in detail.


title: Fast Variance Reduction Method with Stochastic Batch Size

abstract: In this paper we study a family of variance reduction methods with randomized batch size---at each step, the algorithm first randomly chooses the batch size and then selects a batch of samples to conduct a variance-reduced stochastic update. We give the linear convergence rate for this framework for composite functions, and show that the optimal strategy to achieve the optimal convergence rate per data access is to always choose batch size of 1, which is equivalent to the SAGA algorithm. However, due to the presence of cache/disk IO effect in computer architecture, the number of data access cannot reflect the running time because of 1) random memory access is much slower than sequential access, 2) when data is too big to fit into memory, disk seeking takes even longer time. After taking these into account, choosing batch size of $1$ is no longer optimal, so we propose a new algorithm called SAGA++ and show how to calculate the optimal average batch size theoretically. Our algorithm outperforms SAGA and other existing batched and stochastic solvers on real datasets. In addition, we also conduct a precise analysis to compare different update rules for variance reduction methods, showing that SAGA++ converges faster than SVRG in theory.


title: A new class of semi-implicit methods with linear complexity for nonlinear fractional differential equations

abstract: We propose a new class of semi-implicit methods for solving nonlinear fractional differential equations and study their stability. Several versions of our new schemes are proved to be unconditionally stable by choosing suitable parameters. Subsequently, we develop an efficient strategy to calculate the discrete convolution for the approximation of the fractional operator in the semi-implicit method and we derive an error bound of the fast convolution. The memory requirement and computational cost of the present semi-implicit methods with a fast convolution are about $O(N\log n_T)$ and $O(Nn_T\log n_T)$, respectively, where $N$ is a suitable positive integer and $n_T$ is the final number of time steps. Numerical simulations, including the solution of a system of two nonlinear fractional diffusion equations with different fractional orders in two-dimensions, are presented to verify the effectiveness of the semi-implicit methods.


title: Dialog-context aware end-to-end speech recognition

abstract: Existing speech recognition systems are typically built at the sentence level, although it is known that dialog context, e.g. higher-level knowledge that spans across sentences or speakers, can help the processing of long conversations. The recent progress in end-to-end speech recognition systems promises to integrate all available information (e.g. acoustic, language resources) into a single model, which is then jointly optimized. It seems natural that such dialog context information should thus also be integrated into the end-to-end models to improve further recognition accuracy. In this work, we present a dialog-context aware speech recognition model, which explicitly uses context information beyond sentence-level information, in an end-to-end fashion. Our dialog-context model captures a history of sentence-level context so that the whole system can be trained with dialog-context information in an end-to-end manner. We evaluate our proposed approach on the Switchboard conversational speech corpus and show that our system outperforms a comparable sentence-level end-to-end speech recognition system.


title: Algebraic tangent cones of reflexive sheaves

abstract: We study the notion of algebraic tangent cones at singularities of reflexive sheaves. These correspond to extensions of reflexive sheaves across a negative divisor. We show the existence of optimal extensions in a constructive manner, and we prove the uniqueness in a suitable sense. The results here are an algebro-geometric counterpart of our previous study on singularities of Hermitian-Yang-Mills connections.


title: Adapted $\theta$-Scheme and Its Error Estimates for Backward Stochastic Differential Equations

abstract: In this paper we propose a new kind of high order numerical scheme for backward stochastic differential equations(BSDEs). Unlike the traditional $\theta$-scheme, we reduce truncation errors by taking $\theta$ carefully for every subinterval according to the characteristics of integrands. We give error estimates of this nonlinear scheme and verify the order of scheme through a typical numerical experiment.


title: Test without Trust: Optimal Locally Private Distribution Testing

abstract: We study the problem of distribution testing when the samples can only be accessed using a locally differentially private mechanism and focus on two representative testing questions of identity (goodness-of-fit) and independence testing for discrete distributions. We are concerned with two settings: First, when we insist on using an already deployed, general-purpose locally differentially private mechanism such as the popular RAPPOR or the recently introduced Hadamard Response for collecting data, and must build our tests based on the data collected via this mechanism; and second, when no such restriction is imposed, and we can design a bespoke mechanism specifically for testing. For the latter purpose, we introduce the Randomized Aggregated Private Testing Optimal Response (RAPTOR) mechanism which is remarkably simple and requires only one bit of communication per sample. We propose tests based on these mechanisms and analyze their sample complexities. Each proposed test can be implemented efficiently. In each case (barring one), we complement our performance bounds for algorithms with information-theoretic lower bounds and establish sample optimality of our proposed algorithm. A peculiar feature that emerges is that our sample-optimal algorithm based on RAPTOR uses public-coins, and any test based on RAPPOR or Hadamard Response, which are both private-coin mechanisms, requires significantly more samples.


title: Solar Thermochemical Water-Splitting Reaction Enhanced by Hydrogen Permeation Membrane

abstract: The low conversion rate and efficiency always weaken the performance of thermochemical water-splitting reaction. Herein, we, for the first time, conducted the thermodynamic study of a hydrogen permeable membrane (HPM) in an isothermal thermochemical water-splitting reaction driven by solar energy, which has showed a sharply enhanced conversion rate of 87.8% at 1500 oC and 10-5 atm at permeated side (versus 1.26% with oxygen permeation membrane). According to thermodynamic analysis, the first-law thermodynamic efficiency can reach as high as 59.1%. When taking solar-to-electric efficiency and vacuum pump efficiency into account (for converting separation work into solar energy), we simulated the appreciable efficiency of 3.05% at 1500 oC. The numerical model will provide guidance for the actual production of hydrogen by high temperature solar water splitting. Such novel work manifests the great significance of constructing a HPM reactor for efficient solar thermochemical water splitting, which shows a novel approach for high-temperature solar water splitting.


title: Independent Analysis of the Distance to NGC1052-DF2

abstract: The conclusion that the diffuse galaxy NGC1052-DF2 is deficient in dark matter depends critically on the distance, as noted in the original publication and subsequent works. The results of our preliminary assessment of the SBF distance have already been referenced on the web. Here we present the details of a more careful analysis, supply some context, and comment on the remaining uncertainty.


title: Regularized Single and Double Layer Integrals in 3D Stokes Flow

abstract: We present a numerical method for computing the single layer (Stokeslet) and double layer (stresslet) integrals in Stokes flow. The method applies to smooth, closed surfaces in three dimensions, and achieves high accuracy both on and near the surface. The singular Stokeslet and stresslet kernels are regularized and, for the nearly singular case, corrections are added to reduce the regularization error. These corrections are derived analytically for both the Stokeslet and the stresslet using local asymptotic analysis. For the case of evaluating the integrals on the surface, as needed when solving integral equations, we design high order regularizations for both kernels that do not require corrections. This approach is direct in that it does not require grid refinement or special quadrature near the singularity, and therefore does not increase the computational complexity of the overall algorithm. Numerical tests demonstrate the uniform convergence rates for several surfaces in both the singular and near singular cases, as well as the importance of corrections when two surfaces are close to each other.


title: Random conductance models with stable-like jumps: heat kernel estimates and Harnack inequalities

abstract: We establish two-sided heat kernel estimates for random conductance models with non-uniformly elliptic (possibly degenerate) stable-like jumps on graphs. These are long range counterparts of well known two-sided Gaussian heat kernel estimates by M.T. Barlow for nearest neighbor (short range) random walks on the supercritical percolation cluster. Unlike the cases for nearest neighbor conductance models, the idea through parabolic Harnack inequalities does not work, since even elliptic Harnack inequalities do not hold in the present setting. As an application, we establish the local limit theorem for the models.


title: Nonpositive curvature is not coarsely universal

abstract: We prove that not every metric space embeds coarsely into an Alexandrov space of nonpositive curvature. This answers a question of Gromov (1993) and is in contrast to the fact that any metric space embeds coarsely into an Alexandrov space of nonnegative curvature, as shown by Andoni, Naor and Neiman (2015). We establish this statement by proving that a metric space which is $q$-barycentric for some $q\in [1,\infty)$ has metric cotype $q$ with sharp scaling parameter. Our proof utilizes nonlinear (metric space-valued) martingale inequalities and yields sharp bounds even for some classical Banach spaces. This allows us to evaluate the bi-Lipschitz distortion of the $\ell_\infty$ grid $[m]_\infty^n=(\{1,\ldots,m\}^n,\|\cdot\|_\infty)$ into $\ell_q$ for all $q\in (2,\infty)$, from which we deduce the following discrete converse to the fact that $\ell_\infty^n$ embeds with distortion $O(1)$ into $\ell_q$ for $q=O(\log n)$. A rigidity theorem of Ribe (1976) implies that for every $n\in \mathbb{N}$ there exists $m\in \mathbb{N}$ such that if $[m]_\infty^n$ embeds into $\ell_q$ with distortion $O(1)$, then $q$ is necessarily at least a universal constant multiple of $\log n$. Ribe's theorem does not give an explicit upper bound on this $m$, but by the work of Bourgain (1987) it suffices to take $m=n$, and this was the previously best-known estimate for $m$. We show that the above discretization statement actually holds when $m$ is a universal constant.


title: Instance-Dependent PU Learning by Bayesian Optimal Relabeling

abstract: When learning from positive and unlabelled data, it is a strong assumption that the positive observations are randomly sampled from the distribution of $X$ conditional on $Y = 1$, where X stands for the feature and Y the label. Most existing algorithms are optimally designed under the assumption. However, for many real-world applications, the observed positive examples are dependent on the conditional probability $P(Y = 1|X)$ and should be sampled biasedly. In this paper, we assume that a positive example with a higher $P(Y = 1|X)$ is more likely to be labelled and propose a probabilistic-gap based PU learning algorithms. Specifically, by treating the unlabelled data as noisy negative examples, we could automatically label a group positive and negative examples whose labels are identical to the ones assigned by a Bayesian optimal classifier with a consistency guarantee. The relabelled examples have a biased domain, which is remedied by the kernel mean matching technique. The proposed algorithm is model-free and thus do not have any parameters to tune. Experimental results demonstrate that our method works well on both generated and real-world datasets.


title: Observation of a metastable intermediate during solid-solid phase transformation in response to rapid compression

abstract: In order to probe the mechanism of solid-solid phase transformations, we have applied ultrafast shock wave compression (120 picosecond duration) and ultrashort (130 femtosecond) x-ray diffraction at the Linac Coherent Light Source (LCLS) to probe the compression-induced phase transition pathway in zirconium. Surprisingly, rather than transform from alpha-Zr to the more disordered hex-3 equilibrium omega-Zr phase, in its place we find the formation of a non-equilibrium body-centered cubic (bcc) metastable intermediate. Theoretically hypothesized for several decades, this bcc intermediate state has now been found to be dynamically stabilized under uniaxial loading at sub-nanosecond timescales. Molecular dynamics simulations of shock-wave propagation in zirconium predict this transformation via the dynamical intermediate state. In contrast with longer timescale experiments where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time-dependence in the kinetics of phase transformation at extreme conditions.


title: Optimal bail-out dividends problem with transaction cost and capital injection constraint

abstract: We consider the bail-out optimal dividend problem under fixed transaction costs for a L\'evy risk model. Furthermore, we consider the version with a constraint expected net present value of injected capital. To characterize the solution to the aforementioned models, we first solve the bail-out optimal dividend problem under transaction costs and capital injection and show the optimality of reflected (c1; c2)- policies. Next, we introduce the dual Lagrangian problem and show that the complementary slackness conditions are satisfied, characterizing the optimal Lagrange multiplier. Finally, we illustrate our findings with a series of numerical examples.


title: Synopsis of Duality

abstract: Properties and examples of the dual transformation between two planes, which is such that the coordinates of a point in the original plane give the coefficients of a line in the dual plane and the coefficients of a line in the original plane give the coordinates of a point in the dual plane, are examined.


title: Dualities and Phases of 3D N=1 SQCD

abstract: We study gauge theories with N=1 supersymmetry in 2+1 dimensions. We start by calculating the 1-loop effective superpotential for matter in an arbitrary representation. We then restrict ourselves to gauge theories with fundamental matter. Using the 1-loop superpotential, we find a universal form for the phase diagrams of many such gauge theories, which is proven to persist to all orders in perturbation theory using a symmetry argument. This allows us to conjecture new dualities for N=1 gauge theories with fundamental matter. We also show that these dualities are related to results in N=2 supersymmetric gauge theories, which provides further evidence for them.


title: Round-Table Group Optimization for Sequencing Problems

abstract: In this paper, a round-table group optimization (RTGO) algorithm is presented. RTGO is a simple metaheuristic framework using the insights of research on group creativity. In a cooperative group, the agents work in iterative sessions to search innovative ideas in a common problem landscape. Each agent has one base idea stored in its individual memory, and one social idea fed by a round-table group support mechanism in each session. The idea combination and improvement processes are respectively realized by using a recombination search (XS) strategy and a local search (LS) strategy, to build on the base and social ideas. RTGO is then implemented for solving two difficult sequencing problems, i.e., the flowshop scheduling problem and the quadratic assignment problem. The domain-specific LS strategies are adopted from existing algorithms, whereas a general XS class, called socially biased combination (SBX), is realized in a modular form. The performance of RTGO is then evaluated on commonly-used benchmark datasets. Good performance on different problems can be achieved by RTGO using appropriate SBX operators. Furthermore, RTGO is able to outperform some existing methods, including methods using the same LS strategies.


title: Some remarks on bi-f-harmonic maps and f-biharmonic maps

abstract: In this paper, we prove that the class of bi-f-harmonic maps and that of f-biharmonic maps from a conformal manifold of dimension not equal to 2 are the same (Theorem 1.1). We also give several results on nonexistence of proper bi-f-harmonic maps and f-biharmonic maps from complete Riemannian manifolds into nonpositively curved Riemannian manifolds. These include: any bi-f-harmonic map from a compact manifold into a non-positively curved manifold is f-harmonic (Theorem 1.6), and any f-biharmonic (respectively, bi-f-harmonic) map with bounded f and bounded f-bienrgy (respectively, bi-f-energy) from a complete Riemannian manifold into a manifold of strictly negative curvature has rank < 2 everywhere (Theorems 2.2 and 2.3).


title: Extreme precision photometry from the ground with beam-shaping diffusers for K2, TESS and beyond

abstract: The Transiting Exoplanet Survey Satellite (TESS, launched early 2018) is expected to find a multitude of new transiting planet candidates around the nearest and brightest stars. Timely high-precision follow-up observations from the ground are essential in confirming and further characterizing the planet candidates that TESS will find. However, achieving extreme photometric precisions from the ground is challenging, as ground-based telescopes are subject to numerous deleterious atmospheric effects. Beam-shaping diffusers are emerging as a low-cost technology to achieve hitherto unachievable differential photometric precisions from the ground. These diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. In this paper, we expand on our previous work (Stefansson et al. 2017; Stefansson et al. 2018 [submitted]), providing a further detailed discussion of key guidelines when sizing a diffuser for use on a telescope. Furthermore, we present our open source Python package iDiffuse which can calculate the expected PSF size of a diffuser in a telescope system, along with its expected on-sky diffuser-assisted photometric precision for a host star of a given magnitude. We use iDiffuse to show that most ($\sim$80\%) of the planet hosts that TESS will find will be scintillation limited in transit observations from the ground. Although iDiffuse has primarily been developed to plan challenging transit observations using the diffuser on the ARCTIC imager on the ARC 3.5m Telescope at Apache Point observatory, iDiffuse is modular and can be easily extended to calculate the expected diffuser-assisted photometric precisions on other telescopes.


title: Probing the effects of new physics in $\bar{B}^* \to P \ell \bar{\nu}_\ell$ decays

abstract: The significant divergence between the SM predictions and experimental measurements for the ratios, $R_{D^{(*)}}\equiv \mathcal{B}(\bar{B}\to D^{(*)}\tau^- \bar{\nu}_\tau)/ \mathcal{B}(\bar{B}\to D^{(*)} \ell^{\prime-} \bar{\nu}_{\ell^{\prime}})$ with $(\ell^{\prime}=e\,,\mu)$, implies possible hint of new physics in the flavour sector. In this paper, motivated by the "$R_{D^{(*)}}$ puzzle" and abundant $B^*$ data samples at high-luminosity heavy-flavor experiments in the future, we try to probe possible effects of new physics in the semileptonic $\bar{B}^*_{u,d,s} \to P \ell^- \bar{\nu}_\ell$ $(P=D\,,D_s\,,\pi\,,K)$ decays induced by $b \to (u,\,c)\ell^- \bar{\nu}_{\ell}$ transitions in the model-independent vector and scalar scenarios. Using the spaces of NP parameters obtained by fitting to the data of $R_{D}$ and $R_{D^{*}}$, the NP effects on the observables including branching fraction, ratio $R_P^{\ast}$, lepton spin asymmetry and lepton forward-backward asymmetry are studied in detail. We find that the vector type couplings have large effects on the branching fraction and ratio $R_P^{\ast}$. Meanwhile, the scalar type couplings provide significant contributions to all of the observables. The future measurements of these observables in the $\bar{B}^* \to P \ell^- \bar{\nu}_\ell$ decays at the LHCb and Belle-II could provide a way to crosscheck the various NP solutions to the "$R_{D^{(*)}}$ puzzle".


title: Robust path-following control for articulated heavy-duty vehicles

abstract: Path following and lateral stability are crucial issues for autonomous vehicles. Moreover, these problems increase in complexity when handling articulated heavy-duty vehicles due to their poor manoeuvrability, large sizes and mass variation. In addition, uncertainties on mass may have the potential to significantly decrease the performance of the system, even to the point of destabilising it. These parametric variations must be taken into account during the design of the controller. However, robust control techniques usually require offline adjustment of auxiliary tuning parameters, which is not practical, leading to sub-optimal operation. Hence, this paper presents an approach to path-following and lateral control for autonomous articulated heavy-duty vehicles subject to parametric uncertainties by using a robust recursive regulator. The main advantage of the proposed controller is that it does not depend on the offline adjustment of tuning parameters. Parametric uncertainties were assumed to be on the payload, and an $\mathcal{H}_{\infty}$ controller was used for performance comparison. The performance of both controllers is evaluated in a double lane-change manoeuvre. Simulation results showed that the proposed method had better performance in terms of robustness, lateral stability, driving smoothness and safety, which demonstrates that it is a very promising control technique for practical applications.


title: A Bayesian Downscaler Model to Estimate Daily PM2.5 levels in the Continental US

abstract: There has been growing interest in extending the coverage of ground PM2.5 monitoring networks based on satellite remote sensing data. With broad spatial and temporal coverage, satellite based monitoring network has a strong potential to complement the ground monitor system in terms of the spatial-temporal availability of the air quality data. However, most existing calibration models focused on a relatively small spatial domain and cannot be generalized to national-wise study. In this paper, we proposed a statistically reliable and interpretable national modeling framework based on Bayesian downscaling methods with the application to the calibration of the daily ground PM2.5 concentrations across the Continental U.S. using satellite-retrieved aerosol optical depth (AOD) and other ancillary predictors in 2011. Our approach flexibly models the PM2.5 versus AOD and the potential related geographical factors varying across the climate regions and yields spatial and temporal specific parameters to enhance the model interpretability. Moreover, our model accurately predicted the national PM2.5 with a R2 at 70% and generates reliable annual and seasonal PM2.5 concentration maps with its SD. Overall, this modeling framework can be applied to the national scale PM2.5 exposure assessments and also quantify the prediction errors.


title: Privacy in Social Media: Identification, Mitigation and Applications

abstract: The increasing popularity of social media has attracted a huge number of people to participate in numerous activities on a daily basis. This results in tremendous amounts of rich user-generated data. This data provides opportunities for researchers and service providers to study and better understand users' behaviors and further improve the quality of the personalized services. Publishing user-generated data risks exposing individuals' privacy. Users privacy in social media is an emerging task and has attracted increasing attention in recent years. These works study privacy issues in social media from the two different points of views: identification of vulnerabilities, and mitigation of privacy risks. Recent research has shown the vulnerability of user-generated data against the two general types of attacks, identity disclosure and attribute disclosure. These privacy issues mandate social media data publishers to protect users' privacy by sanitizing user-generated data before publishing it. Consequently, various protection techniques have been proposed to anonymize user-generated social media data. There is a vast literature on privacy of users in social media from many perspectives. In this survey, we review the key achievements of user privacy in social media. In particular, we review and compare the state-of-the-art algorithms in terms of the privacy leakage attacks and anonymization algorithms. We overview the privacy risks from different aspects of social media and categorize the relevant works into five groups 1) graph data anonymization and de-anonymization, 2) author identification, 3) profile attribute disclosure, 4) user location and privacy, and 5) recommender systems and privacy issues. We also discuss open problems and future research directions for user privacy issues in social media.


title: A Submillimeter Burst of S255IR~SMA1 - The Rise And Fall Of Its Luminosity

abstract: Temporal photometric variations at near infrared to submillimeter wavelengths have been found in low-mass young stellar objects. These phenomena are generally interpreted as accretion events of star-disk systems with varying accretion rates. There is growing evidence suggesting that similar luminosity flaring also occurs in high-mass star/cluster-forming regions. We report in this Letter the rise and fall of the 900 ${\mu}$m continuum emission and the newly found 349.1 GHz methanol maser emission in the massive star forming region S255IR~SMA1 observed with the Submillimeter Array and the Atacama Large Millimeter/submillimeter Array. The level of flux variation at a factor of $\sim$ 2 at the submillimeter band and the relatively short 2-year duration of this burst suggest that the event is probably similar to those milder and more frequent minor bursts seen in 3D numerical simulations.


title: Generalized Integrative Principal Component Analysis for Multi-Type Data with Block-Wise Missing Structure

abstract: High-dimensional multi-source data are encountered in many fields. Despite recent developments on the integrative dimension reduction of such data, most existing methods cannot easily accommodate data of multiple types (e.g., binary or count-valued). Moreover, multi-source data often have block-wise missing structure, i.e., data in one or more sources may be completely unobserved for a sample. The heterogeneous data types and presence of block-wise missing data pose significant challenges to the integration of multi-source data and further statistical analyses. In this paper, we develop a low-rank method, called Generalized Integrative Principal Component Analysis (GIPCA), for the simultaneous dimension reduction and imputation of multi-source block-wise missing data, where different sources may have different data types. We also devise an adapted BIC criterion for rank estimation. Comprehensive simulation studies demonstrate the efficacy of the proposed method in terms of rank estimation, signal recovery, and missing data imputation. We apply GIPCA to a mortality study. We achieve accurate block-wise missing data imputation and identify intriguing latent mortality rate patterns with sociological relevance.


title: Quantized Densely Connected U-Nets for Efficient Landmark Localization

abstract: In this paper, we propose quantized densely connected U-Nets for efficient visual landmark localization. The idea is that features of the same semantic meanings are globally reused across the stacked U-Nets. This dense connectivity largely improves the information flow, yielding improved localization accuracy. However, a vanilla dense design would suffer from critical efficiency issue in both training and testing. To solve this problem, we first propose order-K dense connectivity to trim off long-distance shortcuts; then, we use a memory-efficient implementation to significantly boost the training efficiency and investigate an iterative refinement that may slice the model size in half. Finally, to reduce the memory consumption and high precision operations both in training and testing, we further quantize weights, inputs, and gradients of our localization network to low bit-width numbers. We validate our approach in two tasks: human pose estimation and face alignment. The results show that our approach achieves state-of-the-art localization accuracy, but using ~70% fewer parameters, ~98% less model size and saving ~75% training memory compared with other benchmark localizers. The code is available at https://github.com/zhiqiangdon/CU-Net.


title: Fisher information matrix for single molecules with stochastic trajectories

abstract: Tracking of objects in cellular environments has become a vital tool in molecular cell biology. A particularly important example is single molecule tracking which enables the study of the motion of a molecule in cellular environments and provides quantitative information on the behavior of individual molecules in cellular environments, which were not available before through bulk studies. Here, we consider a dynamical system where the motion of an object is modeled by stochastic differential equations (SDEs), and measurements are the detected photons emitted by the moving fluorescently labeled object, which occur at discrete time points, corresponding to the arrival times of a Poisson process, in contrast to uniform time points which have been commonly used in similar dynamical systems. The measurements are distributed according to optical diffraction theory, and therefore, they would be modeled by different distributions, e.g., a Born and Wolf profile for an out-of-focus molecule. For some special circumstances, Gaussian image models have been proposed. In this paper, we introduce a stochastic framework in which we calculate the maximum likelihood estimates of the biophysical parameters of the molecular interactions, e.g., diffusion and drift coefficients. More importantly, we develop a general framework to calculate the Cram\'er-Rao lower bound (CRLB), given by the inverse of the Fisher information matrix, for the estimation of unknown parameters and use it as a benchmark in the evaluation of the standard deviation of the estimates. There exists no established method, even for Gaussian measurements, to systematically calculate the CRLB for the general motion model that we consider in this paper. We apply the developed methodology to simulated data of a molecule with linear trajectories and show that the standard deviation of the estimates matches well with the square root of the CRLB.


title: Fractal Dimension and the Persistent Homology of Random Geometric Complexes

abstract: We prove that the fractal dimension of a metric space equipped with an Ahlfors regular measure can be recovered from the persistent homology of random samples. Our main result is that if $x_1,\ldots, x_n$ are i.i.d. samples from a $d$-Ahlfors regular measure on a metric space, and $E^0_\alpha\left(x_1,\ldots,x_n\right)$ denotes the $\alpha$-weight of the minimum spanning tree on $x_1,\ldots,x_n:$ \[E_\alpha^0\left(x_1,\ldots,x_n\right)=\sum_{e\in T\left(x_1,\ldots,x_n\right)} |e|^\alpha\,,\] then there exist constants $0<C_1\leq C_2$ so that \[C_1\leq n^{-\frac{d-\alpha}{d}} E^0_\alpha\left(x_1,\ldots,x_n\right)\leq C_2\,\] with high probability as $n\rightarrow \infty.$ In particular, \[\log\big(E^0_\alpha(x_1,\ldots,x_n)\big)/\log(n)\longrightarrow (d-\alpha)/d\,.\] This is a generalization of a result of Steele (1988) from the non-singular case to the fractal setting. Our result is best possible, in the sense that there exist Ahlfors regular measures for which the limit $\lim_{n\rightarrow\infty} n^{-\frac{d-\alpha}{d}} E^0_\alpha\left(x_1,\ldots,x_n\right)$ does not exist with high probability. We also prove analogous results for weighted sums defined in terms of higher dimensional persistent homology.


title: Strong light-field effects driven by nearly single-cycle 7-fs light field in correlated organic conductors

abstract: We have demonstrated transient charge localization effects with a driving high-frequency field of 7-fs, 1.5-cycle near infrared light in correlated organic conductors. In a layered organic conductor alpha-(BEDT-TTF)2I3 (BEDT-TTF: bis[ethylenedithio]-tetrathiafulvalene), a transient short-range charge order (CO) state is induced in a metallic phase. In contrast to such drastic change in the electronic state from the metal to the transient CO in alpha-(BEDT-TTF)2I3, dynamics of a field-induced reduction of a transfer integral are captured as a red shift of the plasma-like reflectivity edge in a quasi-one-dimensional organic conductor (TMTTF)2AsF6 (TMTTF: tetramethyltetrathiafulvalene). These studies on the field-induced charge localization have been motivated by the theory of dynamical localization on the basis of tight-binding models with no electron correlation under a strong continuous field. However, the results of pump-probe transient reflectivity measurements using nearly single-cycle 7-fs, 11 MV/cm pulses and the theoretical studies which are presented in this review indicate that the pulsed field contributes to the similar phenomenon with the help of a characteristic lattice structure and Coulomb repulsion.


title: Absence of 15-dimensional subalgebras at 16-dimensional Clifford algebra

abstract: It is shown that the 16-dimensional Clifford algebra over any field has no one 15-dimensional subalgebra. Canonical bases are used throughout the determination.


title: Classification of 7-dimensional subalgebras of 8-dimensional Clifford algebra

abstract: All 7-dimensional subalgebras of the 8-dimensional Clifford algebra over the field C of complex numbers are found. Canonical bases are used throughout the determination. It is found that the 8-dimensional Clifford algebra over C has exactly eight 7-dimensional subalgebras and each of these eight is over the complex numbers. It follows that the 8-dimensional Clifford algebra over C has no 7-dimensional subalgebra over the real numbers. Another consequence is that the eight 7-dimensional subalgebras of the 8-dimentional Clifford over C are in fact 7-dimensional subalgebras of all greater dimensional Clifford algebras over C.


title: Deep Learning with Predictive Control for Human Motion Tracking

abstract: We propose to combine model predictive control with deep learning for the task of accurate human motion tracking with a robot. We design the MPC to allow switching between the learned and a conservative prediction. We also explored online learning with a DyBM model. We applied this method to human handwriting motion tracking with a UR-5 robot. The results show that the framework significantly improves tracking performance.


title: Holistic 3D Scene Parsing and Reconstruction from a Single RGB Image

abstract: We propose a computational framework to jointly parse a single RGB image and reconstruct a holistic 3D configuration composed by a set of CAD models using a stochastic grammar model. Specifically, we introduce a Holistic Scene Grammar (HSG) to represent the 3D scene structure, which characterizes a joint distribution over the functional and geometric space of indoor scenes. The proposed HSG captures three essential and often latent dimensions of the indoor scenes: i) latent human context, describing the affordance and the functionality of a room arrangement, ii) geometric constraints over the scene configurations, and iii) physical constraints that guarantee physically plausible parsing and reconstruction. We solve this joint parsing and reconstruction problem in an analysis-by-synthesis fashion, seeking to minimize the differences between the input image and the rendered images generated by our 3D representation, over the space of depth, surface normal, and object segmentation map. The optimal configuration, represented by a parse graph, is inferred using Markov chain Monte Carlo (MCMC), which efficiently traverses through the non-differentiable solution space, jointly optimizing object localization, 3D layout, and hidden human context. Experimental results demonstrate that the proposed algorithm improves the generalization ability and significantly outperforms prior methods on 3D layout estimation, 3D object detection, and holistic scene understanding.


title: On second-order sufficient optimality conditions for $C^1$ vector optimization problems

abstract: In this paper, we present some second-order sufficient conditions in terms of the Demyanov-Pevnyi's second-order directional derivatives for efficiency of $C^1$ vector optimization problems with constraints. Our results improve and generalize conditions obtained by various authors in recent papers.


title: The multiplicity problems for the unitary Ginzburg-Rallis models

abstract: We consider the local multiplicity problems of the analogy of the Ginzburg-Rallis model for the unitary group and the unitary similitude group cases. For the unitary similitude group case, by proving a local trace formula for the model, we are able to prove a multiplicity formula for all tempered representations, which implies that the summation of the multiplicities is equal to $1$ over every tempered local Vogan $L$-packet. For the unitary group case, we also prove a multiplicity formula for all tempered representations which implies that the summation of the multiplicities is equal to $2$ over every tempered local Vogan $L$-packet.


title: How initial size governs core collapse in globular clusters

abstract: Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called "core-collapsed" and "non-core-collapsed" clusters. Here, we use our H\'enon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, $r_v$) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe ($r_v \approx 0.5-5$ pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify "best-fit" models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of $\sim 50-100$ stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main sequence stars, white dwarfs, and sub-subgiants.


title: Remarks on the Casimir force for magnetodielectric media

abstract: Boyer showed that a perfect electrically conducting slab repels a perfect magnetically conducting slab, in contrast to the attractive Casimir force between two identical perfect electrically/magnetically conducting slabs. To gain insight for the difference between the Boyer force and the Casimir force, we present the derivation of the Boyer force using the stress tensor method and then using the method of variation in dielectric. The Green dyadic, in terms of electric and magnetic Green's functions, is presented for an electric medium filling half of space and another magnetic medium filling another half of space such that the two half-spaces are parallel and separated by a distance $a$. We make the observation that the spectral distribution of scattering in a Boyer cavity is that of Fermi-Dirac type, while the spectral distribution of scattering in a Casimir cavity is that of Bose-Einstein type. Based on this observation we conclude that the difference between the Boyer force and the Casimir force is governed by the statistics of the possible scattering in the respective cavities.


title: Reliable Determination of Contact Angle from the Height and Volume of Sessile Drops

abstract: Contact angle is an important parameter in characterizing the wetting properties of fluids. The most common methods for measuring the contact angle is to measure it directly from the profile curve of a sessile drop, a method with certain inherent drawbacks. Here we describe an alternative method that uses the height and volume of a sessile drop as constraints to construct its profile by numerical integration of its two governing differential equations. The integration yields, self consistently, the average value of the contact angle along the entire contact line as well as the footprint radius of the drop and its crown radius of curvature. As a test case, the new method is used to obtain the contact angle of pure water on two different substrates, Teflon and Lucite. For each substrate, four drops ranging in volume from 10 {\mu}l to 40 {\mu}l are used. The computed contact angles are consistent across the four different drop sizes for each substrate and are in agreement with typical literature values.


title: Strangeness and $\Delta$ resonance in compact stars with relativistic-mean-field models

abstract: We explore the effects of strangeness and $\Delta$ resonance in baryonic matter and compact stars within the relativistic-mean-field (RMF) models. The covariant density functional PKDD is adopted for $N$-$N$ interaction, parameters fixed based on finite hypernuclei and neutron stars are taken for the hyperon-meson couplings, and the universal baryon-meson coupling scheme is adopted for the $\Delta$-meson couplings. In light of the recent observations of GW170817 with the dimensionless combined tidal deformability $197 \leq \bar{\Lambda}\leq 720$, we find it is essential to include the $\Delta$ resonances in compact stars, and small $\Delta$-$\rho$ coupling $g_{\rho \Delta}$ is favored if the mass $2.27{}_{-0.15}^{+0.17}\ M_\odot$ of PSR J2215+5135 is confirmed.


title: Generative Adversarial Estimation of Channel Covariance in Vehicular Millimeter Wave Systems

abstract: Enabling highly-mobile millimeter wave (mmWave) systems is challenging because of the huge training overhead associated with acquiring the channel knowledge or designing the narrow beams. Current mmWave beam training and channel estimation techniques do not normally make use of the prior beam training or channel estimation observations. Intuitively, though, the channel matrices are functions of the various elements of the environment. Learning these functions can dramatically reduce the training overhead needed to obtain the channel knowledge. In this paper, a novel solution that exploits machine learning tools, namely conditional generative adversarial networks (GAN), is developed to learn these functions between the environment and the channel covariance matrices. More specifically, the proposed machine learning model treats the covariance matrices as 2D images and learns the mapping function relating the uplink received pilots, which act as RF signatures of the environment, and these images. Simulation results show that the developed strategy efficiently predicts the covariance matrices of the large-dimensional mmWave channels with negligible training overhead.


title: A Constraint on Chern Classes of Strictly Pseudoconvex CR Manifolds

abstract: This short paper gives a constraint on Chern classes of closed strictly pseudoconvex CR manifolds (or equivalently, closed holomorphically fillable contact manifolds) of dimension at least five. We also see that our result is ''optimal'' through some examples.


title: Resonant Hawking radiation as an instability

abstract: We consider a simple model for a black-hole laser: a Bose-Einstein condensate with uniform speed of sound and partially uniform flow, establishing two horizons, a black-hole and a white-hole horizon. Waves confined between the horizons are amplified similar to radiation in a laser cavity. Black-hole lasing appears as an instability with discrete sets of modes given approximately by a round-trip condition. We found that, in addition to the regular Hawking radiation, trans-Planckian radiation does tunnel out of the black-hole laser.


title: Completely Positive Binary Tensors

abstract: A symmetric tensor is completely positive (CP) if it is a sum of tensor powers of nonnegative vectors. This paper characterizes completely positive binary tensors. We show that a binary tensor is completely positive if and only if it satisfies two linear matrix inequalities. This result can be used to determine whether a binary tensor is completely positive or not. When it is, we give an algorithm for computing its cp-rank and the decomposition. When the order is odd, we show that the cp-rank decomposition is unique. When the order is even, we completely characterize when the cp-rank decomposition is unique. We also discuss how to compute the nearest cp-approximation when a binary tensor is not completely positive.


title: Contemplating Visual Emotions: Understanding and Overcoming Dataset Bias

abstract: While machine learning approaches to visual emotion recognition offer great promise, current methods consider training and testing models on small scale datasets covering limited visual emotion concepts. Our analysis identifies an important but long overlooked issue of existing visual emotion benchmarks in the form of dataset biases. We design a series of tests to show and measure how such dataset biases obstruct learning a generalizable emotion recognition model. Based on our analysis, we propose a webly supervised approach by leveraging a large quantity of stock image data. Our approach uses a simple yet effective curriculum guided training strategy for learning discriminative emotion features. We discover that the models learned using our large scale stock image dataset exhibit significantly better generalization ability than the existing datasets without the manual collection of even a single label. Moreover, visual representation learned using our approach holds a lot of promise across a variety of tasks on different image and video datasets.


title: Importance of the Mathematical Foundations of Machine Learning Methods for Scientific and Engineering Applications

abstract: There has been a lot of recent interest in adopting machine learning methods for scientific and engineering applications. This has in large part been inspired by recent successes and advances in the domains of Natural Language Processing (NLP) and Image Classification (IC). However, scientific and engineering problems have their own unique characteristics and requirements raising new challenges for effective design and deployment of machine learning approaches. There is a strong need for further mathematical developments on the foundations of machine learning methods to increase the level of rigor of employed methods and to ensure more reliable and interpretable results. Also as reported in the recent literature on state-of-the-art results and indicated by the No Free Lunch Theorems of statistical learning theory incorporating some form of inductive bias and domain knowledge is essential to success. Consequently, even for existing and widely used methods there is a strong need for further mathematical work to facilitate ways to incorporate prior scientific knowledge and related inductive biases into learning frameworks and algorithms. We briefly discuss these topics and discuss some ideas proceeding in this direction.


title: Allocations of Cold Standbys to Series and Parallel Systems with Dependent Components

abstract: In the context of industrial engineering, cold-standby redundancies allocation strategy is usually adopted to improve the reliability of coherent systems. This paper investigates optimal allocation strategies of cold standbys for series and parallel systems comprised of dependent components with left/right tail weakly stochastic arrangement increasing lifetimes. For the case of heterogeneous and independent matched cold standbys, it is proved that better redundancies should be put in the nodes having weaker [better] components for series [parallel] systems. For the case of homogeneous and independent cold standbys, it is shown that more redundancies should be put in standby with weaker [better] components to enhance the reliability of series [parallel] systems. The results developed here generalize and extend those corresponding ones in the literature to the case of series and parallel systems with dependent components. Numerical examples are also presented to provide guidance for the practical use of our theoretical findings.


title: STTM: A Tool for Short Text Topic Modeling

abstract: Along with the emergence and popularity of social communications on the Internet, topic discovery from short texts becomes fundamental to many applications that require semantic understanding of textual content. As a rising research field, short text topic modeling presents a new and complementary algorithmic methodology to supplement regular text topic modeling, especially targets to limited word co-occurrence information in short texts. This paper presents the first comprehensive open-source package, called STTM, for use in Java that integrates the state-of-the-art models of short text topic modeling algorithms, benchmark datasets, and abundant functions for model inference and evaluation. The package is designed to facilitate the expansion of new methods in this research field and make evaluations between the new approaches and existing ones accessible. STTM is open-sourced at https://github.com/qiang2100/STTM.


title: Contention resolution on a restrained channel

abstract: We examine deterministic broadcasting on multiple-access channels for a scenario when packets are injected continuously by an adversary to the buffers of the devices at rate $\rho$ packages per round. The aim is to maintain system stability, that is, bounded queues. In contrast to previous work we assume that there is a strict limit of available power, defined as the total number of stations allowed to transmit or listen to the channel at a given time, that can never be exceeded. We study how this constraint influences the quality of services with particular focus on stability. We show that in the regime of deterministic algorithms, the significance of energy restriction depends strongly on communication capabilities of broadcasting protocols. For the adaptive and full-sensing protocols, wherein stations may substantially adopt their behavior to the injection pattern, one can construct efficient algorithms using very small amounts of power without sacrificing throughput or stability of the system. In particular, we construct constant-energy adaptive and full sensing protocols stable for $\rho=1$ and any $\rho<1$, respectively, even for worst case (adversarial) injection patterns. Surprisingly, for the case of acknowledgment based algorithms that cannot adopt to the situation on the channel (i.e., their transmitting pattern is fixed in advance), limiting power leads to reducing the throughput. That is, for this class of protocols in order to preserve stability we need to reduce injection rate significantly. We support our theoretical analysis by simulation results of algorithms constructed in the paper. We depict how they work for systems of moderate, realistic sizes. We also provide a comprehensive simulation to compare our algorithms with backoff algorithms, which are common in real-world implementations, in terms of queue sizes and energy consumption.


title: Simulation of the electrical conductivity of two-dimensional films with aligned rod-like conductive fillers: effect of the filler length dispersity

abstract: Using Monte Carlo simulation, we studied the electrical conductivity of two-dimensional films. The films consisted of a poorly conductive host matrix and highly conductive rod-like fillers (rods). The rods were of various lengths fitting a log-normal distribution. They could be aligned along a direction. Special attention was paid to films having completely aligned rods. The impact of length dispersity and the extent of rod alignment on the insulator-to-conductor phase transition was studied. The greater the length dispersity the smaller the critical concentration. The anisotropy of the electrical conductivity was more pronounced in the vicinity of the phase transition. A finite size effect was found to be significant only in the vicinity of the phase transition.


title: On a generalised Blaschke-Santal\`o inequality

abstract: In this paper, we establish a generalised Blaschke-Santal\`o inequality for convex bodies in $\mathbb R^{n+1}$. This inequality gives an upper bound estimate for the product of dual quermassintegrals of convex body and its polar set. Our argument is based on induction on dimensions.


title: Frustrated Magnetism of Pharmacosiderite Comprising Tetrahedral Clusters Arranged in the Primitive Cubic Lattice

abstract: We show that pharmacosiderite is a novel cluster antiferromagnet comprising frustrated regular tetrahedra made of spin-5/2 Fe3+ ions that are arranged in the primitive cubic lattice. The connectivity of the tetrahedra and the inter-cluster interaction of 2.9 K, which is significantly large compared with the intra-cluster interaction of 10.6 K, gives a unique playground for frustration physics. An unconventional antiferromagnetic order is observed below TN ~ 6 K, which is accompanied by a weak ferromagnetic moment and a large fluctuation as evidenced by Mossbauer spectroscopy. A q = 0 magnetic order with the total S = 0 for the tetrahedral cluster is proposed based on the irreducible representation analysis, which may explain the origin of the weak ferromagnetism and fluctuation.


title: CO destruction in protoplanetary disk midplanes: inside versus outside the CO snow surface

abstract: CO has long been thought to be the best tracer to measure gas masses as it is readily detected at (sub)mm wavelengths in many disks. Inferred gas masses from CO in recent ALMA observations of large samples of disks seem inconsistent with their inferred dust masses. The derived gas-to-dust mass ratios from CO are 1-2 orders of magnitude lower than $\sim$100 even if photodissociation and freeze-out are included. Herschel measurements of HD line emission imply gas masses in line with gas-to-dust mass ratios of 100. This suggests that at least one additional mechanism is removing gaseous CO. Here we test the suggestion that the bulk of the CO is chemically processed and that the carbon is sequestered into less volatile species such as CO2, CH3OH and CH4 in the dense, shielded midplane regions of the disk. Using our gas-grain chemical code we perform a parameter exploration and follow the CO abundance evolution over a range of conditions representative of shielded disk midplanes. We find that no chemical processing of CO takes place on 1-3 Myr timescales for low cosmic-ray ionisation rates, $< 5\times 10^{-18}$ s$^{-1}$. Assuming an ionisation rate of $10^{-17}$ s$^{-1}$, more than 90% of the CO is destroyed, but only in the parts of the disk below 30 K. This order of magnitude destruction of CO is robust against the choice of grain-surface reaction rate parameters except for temperatures below 20 K, where there is a strong dependence on the H tunnelling barrier. The low temperatures needed for CO chemical processing indicate that the exact disk temperature structure is important, with warm disks around luminous Herbig stars expected to have little to no CO conversion. For cold disks around T Tauri stars, a large fraction of the emitting CO layer is affected unless the disks are young (< 1 Myr). This can lead to inferred gas masses that are up to two orders of magnitude too low.


title: Real light curves of FK Comae Berenices: Farewell flip-flop

abstract: For seven decades, the widely held view has been that the formation, the migration and the decay of short-lived starspots explain the constantly changing light curves of chromospherically active stars. Our hypothesis is that these deceptive observed light curves are interference of two real constant period light curves of long-lived starspots. The slow motion of these long-lived starspots with respect to each other causes the observed light curve changes. This hypothesis contradicts the current views of starspots. Therefore, we subject it to eight reproducible tests. Our new period finding method detects the two real light curves of FK Com. Our hypothesis is a total success: all real light curve parameters are directly connected to the long-lived starspots which are also seen in the Doppler images of FK Com.These parameters are spatially and temporally correlated just like in the Sun, including weak solar-like surface differential rotation. As for other chromospherically active stars, all eight reproducible tests also support our hypothesis. It explains many spurious phenomena: the rapid light curve changes, the short starspot life-times, the rapid rotation period changes, the active longitudes, the starspot migration, the period cycles, the amplitude cycles and the minimum epoch cycles. It also explains why the light curves and the Doppler images give contradicting surface differential rotation estimates even for the same individual star, as well as the abrupt 180 degrees shifts of activity (the flip-flop events) and the long-term mean light curves. We argue that the current views of starspots need to be revised.


title: Coherence Depletion in Quantum Algorithms

abstract: Besides the superior efficiency compared to their classical counterparts, quantum algorithms known so far are basically task-dependent, and scarcely any common features are shared between them. In this work, however, we show that the depletion of quantum coherence turns out to be a common phenomenon in these algorithms. For all the quantum algorithms that we investigated including Grover's algorithm, Deutsch-Jozsa algorithm and Shor's algorithm, quantum coherence of the system states reduces to the minimum along with the successful execution of the respective processes. Notably, a similar conclusion cannot be drawn using other quantitative measures such as quantum entanglement. Thus, we expect that coherence depletion as a common feature can be useful for devising new quantum algorithms in the future.


title: Self-testing using only marginal information

abstract: The partial states of a multipartite quantum state may carry a lot of information: in some cases, they determine the global state uniquely. This result is known for tomographic information, that is for fully characterized measurements. We extend it to the device-independent framework by exhibiting sets of two-party correlations that self-test pure three-qubit states.


title: Products of involutions of an infinite-dimensional vector space

abstract: We prove that every automorphism of an infinite-dimensional vector space over a field is the product of four involutions, a result that is optimal in the general case. We also characterize the automorphisms that are the product of three involutions. More generally, we study decompositions of automorphisms into three or four factors with prescribed split annihilating polynomials of degree $2$.


title: Neutrino Losses in Type I Thermonuclear X-ray Bursts: An Improved Nuclear Energy Generation Approximation

abstract: Type I X-ray bursts are thermonuclear explosions on the surface of accreting neutron stars. Hydrogen rich X-ray bursts burn protons far from the line of stability and can release energy in the form of neutrinos from $\beta$-decays. We have estimated, for the first time, the neutrino fluxes of Type I bursts for a range of initial conditions based on the predictions of a 1D implicit hydrodynamics code, KEPLER, which calculates the complete nuclear reaction network. We find that neutrino losses are between $6.7 \times 10^{-5}$ and $0.14$ of the total energy per nucleon, Q$_{nuc}$, depending upon the hydrogen fraction in the fuel. These values are significantly below the $35\,\%$ value for neutrino losses often adopted in recent literature for the rp-process. The discrepancy arises because it is only at $\beta$-decays that $\approx35\,\%$ of energy is lost due to neutrino emission, whereas there are no neutrino losses in $(p,\gamma)$ and $(\alpha,p)$ reactions. Using the total measured burst energies from KEPLER for a range of initial conditions, we have determined an approximation formula for the total energy per nucleon released during an X-ray burst, Q$_{nuc}$=1.31+6.95$\bar{X} - 1.92\bar{X}^2 $MeV/nucleon, where $\bar{X}$ is the average hydrogen mass fraction of the ignition column, with an RMS error of $0.052\,$Mev/nucleon. We provide a detailed analysis of the nuclear energy output of a burst and find an incomplete extraction of mass excess in the burst fuel, with $14\,\%$ of the mass excess in the fuel not being extracted.


title: Local numerical equivalences and Okounkov bodies in higher dimensions

abstract: We continue to explore the numerical nature of the Okounkov bodies focusing on the local behaviors near given points. More precisely, we show that the set of Okounkov bodies of a pseudoeffective divisor with respect to admissible flags centered at a fixed point determines the local numerical equivalence class of divisors which is defined in terms of refined divisorial Zariski decompositions. Our results extend Ro\'{e}'s work on surfaces to higher dimensional varieties although our proof is essentially different in nature.


title: Hierarchical Clustering better than Average-Linkage

abstract: Hierarchical Clustering (HC) is a widely studied problem in exploratory data analysis, usually tackled by simple agglomerative procedures like average-linkage, single-linkage or complete-linkage. In this paper we focus on two objectives, introduced recently to give insight into the performance of average-linkage clustering: a similarity based HC objective proposed by [Moseley and Wang, 2017] and a dissimilarity based HC objective proposed by [Cohen-Addad et al., 2018]. In both cases, we present tight counterexamples showing that average-linkage cannot obtain better than 1/3 and 2/3 approximations respectively (in the worst-case), settling an open question raised in [Moseley and Wang, 2017]. This matches the approximation ratio of a random solution, raising a natural question: can we beat average-linkage for these objectives? We answer this in the affirmative, giving two new algorithms based on semidefinite programming with provably better guarantees.


title: Segmental Audio Word2Vec: Representing Utterances as Sequences of Vectors with Applications in Spoken Term Detection

abstract: While Word2Vec represents words (in text) as vectors carrying semantic information, audio Word2Vec was shown to be able to represent signal segments of spoken words as vectors carrying phonetic structure information. Audio Word2Vec can be trained in an unsupervised way from an unlabeled corpus, except the word boundaries are needed. In this paper, we extend audio Word2Vec from word-level to utterance-level by proposing a new segmental audio Word2Vec, in which unsupervised spoken word boundary segmentation and audio Word2Vec are jointly learned and mutually enhanced, so an utterance can be directly represented as a sequence of vectors carrying phonetic structure information. This is achieved by a segmental sequence-to-sequence autoencoder (SSAE), in which a segmentation gate trained with reinforcement learning is inserted in the encoder. Experiments on English, Czech, French and German show very good performance in both unsupervised spoken word segmentation and spoken term detection applications (significantly better than frame-based DTW).


title: Grassmannian Learning: Embedding Geometry Awareness in Shallow and Deep Learning

abstract: Modern machine learning algorithms have been adopted in a range of signal-processing applications spanning computer vision, natural language processing, and artificial intelligence. Many relevant problems involve subspace-structured features, orthogonality constrained or low-rank constrained objective functions, or subspace distances. These mathematical characteristics are expressed naturally using the Grassmann manifold. Unfortunately, this fact is not yet explored in many traditional learning algorithms. In the last few years, there have been growing interests in studying Grassmann manifold to tackle new learning problems. Such attempts have been reassured by substantial performance improvements in both classic learning and learning using deep neural networks. We term the former as shallow and the latter deep Grassmannian learning. The aim of this paper is to introduce the emerging area of Grassmannian learning by surveying common mathematical problems and primary solution approaches, and overviewing various applications. We hope to inspire practitioners in different fields to adopt the powerful tool of Grassmannian learning in their research.

categories: cs.LG cs.CV cs.IT eess.SP math.IT stat.ML


title: Eigenvector sensitivity under general and structured perturbations of tridiagonal Toeplitz-type matrices

abstract: The sensitivity of eigenvalues of structured matrices under general or structured perturbations of the matrix entries has been thoroughly studied in the literature. Error bounds are available and the pseudospectrum can be computed to gain insight. Few investigations have focused on analyzing the sensitivity of eigenvectors under general or structured perturbations. The present paper discusses this sensitivity for tridiagonal Toeplitz and Toeplitz-type matrices.


title: Anonymity and Confidentiality in Secure Distributed Simulation

abstract: Research on data confidentiality, integrity and availability is gaining momentum in the ICT community, due to the intrinsically insecure nature of the Internet. While many distributed systems and services are now based on secure communication protocols to avoid eavesdropping and protect confidentiality, the techniques usually employed in distributed simulations do not consider these issues at all. This is probably due to the fact that many real-world simulators rely on monolithic, offline approaches and therefore the issues above do not apply. However, the complexity of the systems to be simulated, and the rise of distributed and cloud based simulation, now impose the adoption of secure simulation architectures. This paper presents a solution to ensure both anonymity and confidentiality in distributed simulations. A performance evaluation based on an anonymized distributed simulator is used for quantifying the performance penalty for being anonymous. The obtained results show that this is a viable solution.


title: Short-range and tensor correlations in 4He and 8Be studied with antisymmetrized quasi cluster model

abstract: We apply tensor version of antisymmetrized quasi cluster model (AQCM-T) to 4He and 8Be while focusing on the NN correlations in alpha clusters. We adopt the NN interactions including realistic ones containing a repulsive core for the central part in addition to the tensor part. In 4He, the pn pair in the 3D channel has been known to play a decisive role in the tensor correlation and the framework is capable of treating not only this channel but also the NN correlations in the 1S and 3S channels. In 8Be, when two alpha clusters approach, the 3D pair is suppressed because of the Pauli blocking effect, which induces the decrease of the 3S component through the 3S-3D coupling. This effect results in the reduction of the attractive effect of the central-even interaction in the middle-range region.


title: Robust Pricing with Refunds

abstract: Before purchase, a buyer of an experience good learns about the product's fit using various information sources, including some of which the seller may be unaware of. The buyer, however, can conclusively learn the fit only after purchasing and trying out the product. We show that the seller can use a simple mechanism to best take advantage of the buyer's post-purchase learning to maximize his guaranteed-profit. We show that this mechanism combines a generous refund, which performs well when the buyer is relatively informed, with non-refundable random discounts, which work well when the buyer is relatively uninformed.


title: Deep Stacked Stochastic Configuration Networks for Lifelong Learning of Non-Stationary Data Streams

abstract: The concept of SCN offers a fast framework with universal approximation guarantee for lifelong learning of non-stationary data streams. Its adaptive scope selection property enables for proper random generation of hidden unit parameters advancing conventional randomized approaches constrained with a fixed scope of random parameters. This paper proposes deep stacked stochastic configuration network (DSSCN) for continual learning of non-stationary data streams which contributes two major aspects: 1) DSSCN features a self-constructing methodology of deep stacked network structure where hidden unit and hidden layer are extracted automatically from continuously generated data streams; 2) the concept of SCN is developed to randomly assign inverse covariance matrix of multivariate Gaussian function in the hidden node addition step bypassing its computationally prohibitive tuning phase. Numerical evaluation and comparison with prominent data stream algorithms under two procedures: periodic hold-out and prequential test-then-train processes demonstrate the advantage of proposed methodology.


title: The perverse filtration for the Hitchin fibration is locally constant

abstract: We prove that the perverse Leray filtration for the Hitchin morphism is locally constant in families, thus providing some evidence towards the validity of the $P=W$ conjecture due to de Cataldo, Hausel and Migliorini in non Abelian Hodge theory.


title: Focusing NLS with inverse square potential

abstract: In this paper, we utilize the method in Dodson-Murphy [4] to establish the radial scattering result for the focusing nonlinear Schr\"odinger equation with inverse square potential $i\pa_tu-\la u=-|u|^{p-1}u$ in the energy space $H^1_a(\R^d)$ in dimensions $d\geq3$, which extends the result of [10,11] to higher dimensions cases but with radial initial data. The new ingredient is to establish the dispersive estimate for radial function and overcome the weak dispersive estimate when $a<0$.


title: DeePathology: Deep Multi-Task Learning for Inferring Molecular Pathology from Cancer Transcriptome

abstract: Despite great advances, molecular cancer pathology is often limited to the use of a small number of biomarkers rather than the whole transcriptome, partly due to computational challenges. Here, we introduce a novel architecture of Deep Neural Networks (DNNs) that is capable of simultaneous inference of various properties of biological samples, through multi-task and transfer learning. It encodes the whole transcription profile into a strikingly low-dimensional latent vector of size 8, and then recovers mRNA and miRNA expression profiles, tissue and disease type from this vector. This latent space is significantly better than the original gene expression profiles for discriminating samples based on their tissue and disease. We employed this architecture on mRNA transcription profiles of 10787 clinical samples from 34 classes (one healthy and 33 different types of cancer) from 27 tissues. Our method significantly outperforms prior works and classical machine learning approaches in predicting tissue-of-origin, normal or disease state and cancer type of each sample. For tissues with more than one type of cancer, it reaches 99.4\% accuracy in identifying the correct cancer subtype. We also show this system is very robust against noise and missing values. Collectively, our results highlight applications of artificial intelligence in molecular cancer pathology and oncological research. DeePathology is freely available at \url{https://github.com/SharifBioinf/DeePathology}.


title: Toward an algebraic theory of Welschinger invariants

abstract: Let $S$ be a smooth del Pezzo surface over a field $k$ of characteristic $\neq 2, 3$. We define an invariant in the Grothendieck-Witt ring $GW(k)$ for "counting" rational curves in a curve class $D$ of fixed positive degree (with respect to the anti-canonical bundle $-K_S$) and containing a collection of distinct closed points $\mathfrak{p}=\sum_ip_i$ of total degree $r:=-D\cdot K_S-1$ on $S$. This recovers Welschinger's invariant in case $k=\mathbb{R}$ by applying the signature map. The main result is that this quadratic invariant depends only on the $\mathbb{A}^1$-connected component containing $\mathfrak{p}$ in $Sym^r(S)^0(k)$, where $Sym^r(S)^0$ is the open subscheme of $Sym^r(S)$ parametrizing geometrically reduced 0-cycles.


title: Biodiversity and robustness of large ecosystems

abstract: We study the biodiversity problem for resource competition systems with extinctions and self-limitation effects. Our main result establishes estimates of biodiversity in terms of the fundamental parameters of the model. We also prove the global stability of solutions for systems with extinctions and large turnover rate. We show that when the extinction threshold is distinct from zero, the large time dynamics of system is fundamentally non-predictable. In the last part of the paper we obtain explicit analytical estimates of ecosystem robustness with respect to variations of resource supply which support the $R^*$ rule for a system with random parameters.


title: Speeding Up Distributed Gradient Descent by Utilizing Non-persistent Stragglers

abstract: Distributed gradient descent (DGD) is an efficient way of implementing gradient descent (GD), especially for large data sets, by dividing the computation tasks into smaller subtasks and assigning to different computing servers (CSs) to be executed in parallel. In standard parallel execution, per-iteration waiting time is limited by the execution time of the straggling servers. Coded DGD techniques have been introduced recently, which can tolerate straggling servers via assigning redundant computation tasks to the CSs. In most of the existing DGD schemes, either with coded computation or coded communication, the non-straggling CSs transmit one message per iteration once they complete all their assigned computation tasks. However, although the straggling servers cannot complete all their assigned tasks, they are often able to complete a certain portion of them. In this paper, we allow multiple transmissions from each CS at each iteration in order to make sure a maximum number of completed computations can be reported to the aggregating server (AS), including the straggling servers. We numerically show that the average completion time per iteration can be reduced significantly by slightly increasing the communication load per server.

categories: cs.IT cs.DC cs.LG eess.SP math.IT stat.ML


title: Signal propagation on $\kappa$-Minkowski spacetime and non-local two-point functions

abstract: We study the propagation of quantum fields on $\kappa$-Minkowsi spacetime. Starting from the non-commutative partition function for a free field written in momentum space we derive the Feynman propagator and analyze the non-trivial singularity structure determined by the group manifold geometry of momentum space. The additional contributions due to such singularity structure result in a deformed field propagation which can be alternatively described in terms of an ordinary field propagation determined by a source with a blurred spacetime profile. We show that the $\kappa$-deformed Feynman propagator can be written in terms of vacuum expectation values of a commutative non-local quantum field. For sub-Planckian modes the $\kappa$-deformed propagator corresponds to the vacuum expectation value of the time-ordered product of non-local field operators while for trans-Plankian modes this is replaced by the Hadamard two-point function, the vacuum expectation value of the anti-commutator of non-local field operators.


title: Optimal Subpattern Assignment Metric for Multiple Tracks (OSPAMT Metric)

abstract: In this paper, we propose a new metric which measures the distance between two finite sets of tracks (a track is a path of either a real or estimated target). This metric is based on the same principle as the Optimal Subpattern Assignment (OSPA) metric devised by Schuhmacher et al. Importantly however, the new metric measures the distance between two finite sets of tracks whereas the OSPA metric measures the distance between two finite sets of target states. By also considering the properties of false tracks, missed tracks and many tracks assigned to one track situations caused by missed detections and false alarms, the minimization of all distances between tracks across two finite sets of tracks employed by the new OSPAMT metric enables performance evaluation of multi-target tracking (MTT) algorithms in a more comprehensive and accurate manner than existing metrics such as the OSPA metric and the enhanced OSPAT metric introduced by Ristic et al which measures the distance between two finite sets of labeled target states.


title: Modularity of Erd\H{o}s-R\'enyi random graphs

abstract: For a given graph $G$, each partition of the vertices has a modularity score, with higher values indicating that the partition better captures community structure in $G$. The modularity $q^*(G)$ of the graph $G$ is defined to be the maximum over all vertex partitions of the modularity score, and satisfies $0\leq q^*(G) < 1$. Modularity is at the heart of the most popular algorithms for community detection. We investigate the behaviour of the modularity of the Erd\H{o}s-R\'enyi random graph $G_{n,p}$ with $n$ vertices and edge-probability $p$. Two key findings are that the modularity is $1+o(1)$ with high probability (whp) for $np$ up to $1+o(1)$ and no further; and when $np \geq 1$ and $p$ is bounded below 1, it has order $(np)^{-1/2}$ whp, in accord with a conjecture by Reichardt and Bornholdt in 2006. We also show that the modularity of a graph is robust to changes in a few edges, in contrast to the sensitivity of optimal vertex partitions.


title: A Generic Multi-Projection-Center Model and Calibration Method for Light Field Cameras

abstract: Light field cameras can capture both spatial and angular information of light rays, enabling 3D reconstruction by a single exposure. The geometry of 3D reconstruction is affected by intrinsic parameters of a light field camera significantly. In the paper, we propose a multi-projection-center (MPC) model with 6 intrinsic parameters to characterize light field cameras based on traditional two-parallel-plane (TPP) representation. The MPC model can generally parameterize light field in different imaging formations, including conventional and focused light field cameras. By the constraints of 4D ray and 3D geometry, a 3D projective transformation is deduced to describe the relationship between geometric structure and the MPC coordinates. Based on the MPC model and projective transformation, we propose a calibration algorithm to verify our light field camera model. Our calibration method includes a close-form solution and a non-linear optimization by minimizing re-projection errors. Experimental results on both simulated and real scene data have verified the performance of our algorithm.


title: Direction curves of tangent indicatrix of a curve

abstract: In this paper, we define some new associated curves as integral curves of a vector field generated by Frenet vectors of tangent indicatrix of a curve in Euclidean 3-space. We give some relationships between curvatures of these curves. By using these associated curves, we give some methods to construct helices and slant helices from some special spherical curves such as circles on unit sphere, spherical helices, and spherical slant helices. Finally, we give some related examples.


title: SAM-RCNN: Scale-Aware Multi-Resolution Multi-Channel Pedestrian Detection

abstract: Convolutional neural networks (CNN) have enabled significant improvements in pedestrian detection owing to the strong representation ability of the CNN features. Recently, aggregating features from multiple layers of a CNN has been considered as an effective approach, however, the same approach regarding feature representation is used for detecting pedestrians of varying scales. Consequently, it is not guaranteed that the feature representation for pedestrians of a particular scale is optimised. In this paper, we propose a Scale-Aware Multi-resolution (SAM) method for pedestrian detection which can adaptively select multi-resolution convolutional features according to pedestrian sizes. The proposed SAM method extracts the appropriate CNN features that have strong representation ability as well as sufficient feature resolution, given the size of the pedestrian candidate output from a region proposal network. Moreover, we propose an enhanced SAM method, termed as SAM+, which incorporates complementary features channels and achieves further performance improvement. Evaluations on the challenging Caltech and KITTI pedestrian benchmarks demonstrate the superiority of our proposed method.


title: Extending some results on the second neighborhood conjecture

abstract: A vertex in a directed graph is said to have a large second neighborhood if it has at least as many second out-neighbors as out-neighbors. The Second Neighborhood Conjecture, first stated by Seymour, asserts that there is a vertex having a large second neighborhood in every oriented graph (a directed graph without loops or digons). We prove that oriented graphs whose missing edges can be partitioned into a (possibly empty) matching and a (possibly empty) star satisfy this conjecture. This generalizes a result of Fidler and Yuster. An implication of our result is that every oriented graph without a sink and whose missing edges form a (possibly empty) matching has at least two vertices with large second neighborhoods. This is a strengthening of a theorem of Havet and Thomasse, who showed that the same holds for tournaments without a sink. Moreover, we also show that the conjecture is true for oriented graphs whose vertex set can be partitioned into an independent set and a 2-degenerate graph.


title: Generation of second-mode internal solitary waves during the winter by a shoaling internal tide in the northern South China Sea

abstract: Field measurements of second-mode internal solitary waves (mode-2 ISWs) during the winter on the upper continental slope of the northern South China Sea were reported in Yang et al. (2009), but their generation mechanism remains elusive. We investigated this issue with a multi-modal evolution model and theoretical analysis, which suggest that the observed mode-2 ISWs were generated by a shoaling mode-2 semidiurnal internal tide (IT) that emanated from the Luzon Strait. The results show that two groups of mode-2 ISWs usually appear within one semidiurnal tidal period, successively riding on expanded and subsequently compressed pycnoclines. The number of wave groups largely depends on the amplitudes of the ITs; that is, a larger IT produces larger mode-2 ISWs. Furthermore, intermodal coupling dominates the evolution of a mode-1 IT, highlighting the importance of considering mode scattering in the propagation of low-mode ITs.


title: Entanglement Protection via Periodic Environment Resetting in Continuous Time Quantum Dynamical Processes

abstract: The temporal evolution of entanglement between a noisy system and an ancillary system is analyzed in the context of continuous time open quantum system dynamics. Focusing on a couple of analytically solvable models for qubit systems, we study how Markovian and non-Markovian characteristics influence the problem, discussing in particular their associated entanglement-breaking regimes. These performances are compared with those one could achieve when the environment of the system is forced to return to its input configuration via periodic instantaneous resetting procedures.