title: Direct photon production at low transverse momentum in proton-proton collisions at $\sqrt{s}$ = 2.76 and 8 TeV

abstract: Measurements of inclusive and direct photon production at mid-rapidity in pp collisions at $\sqrt{s}=2.76$ and 8 TeV are presented by the ALICE experiment at the LHC. The results are reported in transverse momentum ranges of $0.4<p_{T}<10$ GeV/$c$ and $0.3<p_{T}<16$ GeV/$c$, respectively. Photons are detected with the electromagnetic calorimeter~(EMCal) and via reconstruction of e$^+$e$^-$ pairs from conversions in the ALICE detector material using the central tracking system. For the final measurement of the inclusive photon spectra the results are combined in the overlapping $p_{T}$ interval of both methods. Direct photon spectra, or their upper limits at 90% C.L. are extracted using the direct photon excess ratio $R_{\gamma}$, which quantifies the ratio of inclusive photons over decay photons generated with a decay-photon simulation. An additional hybrid method, combining photons reconstructed from conversions with those identified in the EMCal, is used for the combination of the direct photon excess ratio $R_{\gamma}$, as well as the extraction of direct photon spectra or their upper limits. While no significant signal of direct photons is seen over the full $p_{T}$ range, $R_{\gamma}$ for $p_{T}>7$ GeV/$c$ is at least one $\sigma$ above unity and consistent with expectations from next-to-leading order pQCD calculations.


title: Unconditional convergence of a fast two-level linearized algorithm for semilinear subdiffusion equations

abstract: A fast two-level linearized scheme with unequal time-steps is constructed and analyzed for an initial-boundary-value problem of semilinear subdiffusion equations. The two-level fast L1 formula of the Caputo derivative is derived based on the sum-of-exponentials technique. The resulting fast algorithm is computationally efficient in long-time simulations because it significantly reduces the computational cost $O(MN^2)$ and storage $O(MN)$ for the standard L1 formula to $O(MN\log N)$ and $O(M\log N)$, respectively, for $M$ grid points in space and $N$ levels in time. The nonuniform time mesh would be graded to handle the typical singularity of the solution near the time $t=0$, and Newton linearization is used to approximate the nonlinearity term. Our analysis relies on three tools: a new discrete fractional Gr\"{o}nwall inequality, a global consistency analysis and a discrete $H^2$ energy method. A sharp error estimate reflecting the regularity of solution is established without any restriction on the relative diameters of the temporal and spatial mesh sizes. Numerical examples are provided to demonstrate the effectiveness of our approach and the sharpness of error analysis.


title: WebSeg: Learning Semantic Segmentation from Web Searches

abstract: In this paper, we improve semantic segmentation by automatically learning from Flickr images associated with a particular keyword, without relying on any explicit user annotations, thus substantially alleviating the dependence on accurate annotations when compared to previous weakly supervised methods. To solve such a challenging problem, we leverage several low-level cues (such as saliency, edges, etc.) to help generate a proxy ground truth. Due to the diversity of web-crawled images, we anticipate a large amount of 'label noise' in which other objects might be present. We design an online noise filtering scheme which is able to deal with this label noise, especially in cluttered images. We use this filtering strategy as an auxiliary module to help assist the segmentation network in learning cleaner proxy annotations. Extensive experiments on the popular PASCAL VOC 2012 semantic segmentation benchmark show surprising good results in both our WebSeg (mIoU = 57.0%) and weakly supervised (mIoU = 63.3%) settings.


title: Three Birds One Stone: A General Architecture for Salient Object Segmentation, Edge Detection and Skeleton Extraction

abstract: In this paper, we aim at solving pixel-wise binary problems, including salient object segmentation, skeleton extraction, and edge detection, by introducing a unified architecture. Previous works have proposed tailored methods for solving each of the three tasks independently. Here, we show that these tasks share some similarities that can be exploited for developing a unified framework. In particular, we introduce a horizontal cascade, each component of which is densely connected to the outputs of previous component. Stringing these components together allows us to effectively exploit features across different levels hierarchically to effectively address the multiple pixel-wise binary regression tasks. To assess the performance of our proposed network on these tasks, we carry out exhaustive evaluations on multiple representative datasets. Although these tasks are inherently very different, we show that our unified approach performs very well on all of them and works far better than current single-purpose state-of-the-art methods. All the code in this paper will be publicly available.


title: A Classification Model for Sensing Human Trust in Machines Using EEG and GSR

abstract: Today, intelligent machines \emph{interact and collaborate} with humans in a way that demands a greater level of trust between human and machine. A first step towards building intelligent machines that are capable of building and maintaining trust with humans is the design of a sensor that will enable machines to estimate human trust level in real-time. In this paper, two approaches for developing classifier-based empirical trust sensor models are presented that specifically use electroencephalography (EEG) and galvanic skin response (GSR) measurements. Human subject data collected from 45 participants is used for feature extraction, feature selection, classifier training, and model validation. The first approach considers a general set of psychophysiological features across all participants as the input variables and trains a classifier-based model for each participant, resulting in a trust sensor model based on the general feature set (i.e., a "general trust sensor model"). The second approach considers a customized feature set for each individual and trains a classifier-based model using that feature set, resulting in improved mean accuracy but at the expense of an increase in training time. This work represents the first use of real-time psychophysiological measurements for the development of a human trust sensor. Implications of the work, in the context of trust management algorithm design for intelligent machines, are also discussed.


title: A Decision Tree Approach to Predicting Recidivism in Domestic Violence

abstract: Domestic violence (DV) is a global social and public health issue that is highly gendered. Being able to accurately predict DV recidivism, i.e., re-offending of a previously convicted offender, can speed up and improve risk assessment procedures for police and front-line agencies, better protect victims of DV, and potentially prevent future re-occurrences of DV. Previous work in DV recidivism has employed different classification techniques, including decision tree (DT) induction and logistic regression, where the main focus was on achieving high prediction accuracy. As a result, even the diagrams of trained DTs were often too difficult to interpret due to their size and complexity, making decision-making challenging. Given there is often a trade-off between model accuracy and interpretability, in this work our aim is to employ DT induction to obtain both interpretable trees as well as high prediction accuracy. Specifically, we implement and evaluate different approaches to deal with class imbalance as well as feature selection. Compared to previous work in DV recidivism prediction that employed logistic regression, our approach can achieve comparable area under the ROC curve results by using only 3 of 11 available features and generating understandable decision trees that contain only 4 leaf nodes.


title: The hard potential relativistic Boltzmann equation in the whole space

abstract: We study the Cauchy problem for the relativistic Boltzmann equation near relativistic Maxwellians in the whole space. The purpose of this article is to handle hard potentials, and for initial data with finite $L^\infty$ norm, to construct global in time mild solutions. We also prove the optimal time decay rates of the solutions.


title: Constraints on the Density and Internal Strength of 1I/'Oumuamua

abstract: 1I/'Oumuamua was discovered by the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS 1) on 19 October 2017. Unlike all previously discovered minor planets this object was determined to have eccentricity $e > 1.0$, suggesting an interstellar origin. Since this discovery and within the limited window of opportunity, several photometric and spectroscopic studies of the object have been made. Using the measured light curve amplitudes and rotation periods we find that, under the assumption of a triaxial ellipsoid, a density range $1500 < \rho < 2800$ kg m$^{-3}$ matches the observations and no significant cohesive strength is required. We also determine that an aspect ratio of $6\pm 1:1$ is most likely after accounting for phase-angle effects and considering the potential effect of surface properties. This elongation is still remarkable but less than some other estimates.


title: Epitaxial growth of highly strained antimonene on Ag (111)

abstract: The synthesis of antimonene, which is a promising group-V 2D material for both fundamental studies and technological applications, remains highly challenging. Thus far, it has been synthesized only by exfoliation or growth on a few substrates. In this study, we show that thin layers of antimonene can be grown on Ag (111) by molecular beam epitaxy. High-resolution scanning tunneling microscopy combined with theoretical calculations revealed that the submonolayer Sb deposited on a Ag (111) surface forms a layer of AgSb2 surface alloy upon annealing. Further deposition of Sb on the AgSb2 surface alloy causes an epitaxial layer of Sb to form, which is identified as antimonene with a buckled honeycomb structure. More interestingly, the lattice constant of the epitaxial antimonene (5 {\AA}) is much larger than that of freestanding antimonene, indicating a high tensile strain of more than 20%. This kind of large strain is expected to make the antimonene a highly promising candidate for room-temperature quantum spin Hall material.


title: Accelerating Empowerment Computation with UCT Tree Search

abstract: Models of intrinsic motivation present an important means to produce sensible behaviour in the absence of extrinsic rewards. Applications in video games are varied, and range from intrinsically motivated general game-playing agents to non-player characters such as companions and enemies. The information-theoretic quantity of Empowerment is a particularly promising candidate motivation to produce believable, generic and robust behaviour. However, while it can be used in the absence of external reward functions that would need to be crafted and learned, empowerment is computationally expensive. In this paper, we propose a modified UCT tree search method to mitigate empowerment's computational complexity in discrete and deterministic scenarios. We demonstrate how to modify a Monte-Carlo Search Tree with UCT to realise empowerment maximisation, and discuss three additional modifications that facilitate better sampling. We evaluate the approach both quantitatively, by analysing how close our approach gets to the baseline of exhaustive empowerment computation with varying amounts of computational resources, and qualitatively, by analysing the resulting behaviour in a Minecraft-like scenario.


title: Towards Human-Machine Cooperation: Self-supervised Sample Mining for Object Detection

abstract: Though quite challenging, leveraging large-scale unlabeled or partially labeled images in a cost-effective way has increasingly attracted interests for its great importance to computer vision. To tackle this problem, many Active Learning (AL) methods have been developed. However, these methods mainly define their sample selection criteria within a single image context, leading to the suboptimal robustness and impractical solution for large-scale object detection. In this paper, aiming to remedy the drawbacks of existing AL methods, we present a principled Self-supervised Sample Mining (SSM) process accounting for the real challenges in object detection. Specifically, our SSM process concentrates on automatically discovering and pseudo-labeling reliable region proposals for enhancing the object detector via the introduced cross image validation, i.e., pasting these proposals into different labeled images to comprehensively measure their values under different image contexts. By resorting to the SSM process, we propose a new AL framework for gradually incorporating unlabeled or partially labeled data into the model learning while minimizing the annotating effort of users. Extensive experiments on two public benchmarks clearly demonstrate our proposed framework can achieve the comparable performance to the state-of-the-art methods with significantly fewer annotations.


title: Bypassing Feature Squeezing by Increasing Adversary Strength

abstract: Feature Squeezing is a recently proposed defense method which reduces the search space available to an adversary by coalescing samples that correspond to many different feature vectors in the original space into a single sample. It has been shown that feature squeezing defenses can be combined in a joint detection framework to achieve high detection rates against state-of-the-art attacks. However, we demonstrate on the MNIST and CIFAR-10 datasets that by increasing the adversary strength of said state-of-the-art attacks, one can bypass the detection framework with adversarial examples of minimal visual distortion. These results suggest for proposed defenses to validate against stronger attack configurations.


title: Extensions of Bernstein's Lethargy Theorem

abstract: In this paper, we examine the aptly-named "Lethargy Theorem" of Bernstein and survey its recent extensions. We show that one of these extensions shrinks the interval for best approximation by half while the other gives a surprising connection to the space of bounded linear operators between two Banach spaces.


title: Homology cobordism and triangulations

abstract: The study of triangulations on manifolds is closely related to understanding the three-dimensional homology cobordism group. We review here what is known about this group, with an emphasis on the local equivalence methods coming from Pin(2)- equivariant Seiberg-Witten Floer spectra and involutive Heegaard Floer homology.


title: Rate-distortion functions of non-stationary Markoff chains and their block-independent approximations

abstract: It is proved that the limit of the normalized rate-distortion functions of block independent approximations of an irreducible, aperiodic Markoff chain is independent of the initial distribution of the Markoff chain and thus, is also equal to the rate-distortion function of the Markoff chain.


title: General Relativity as a Hybrid theory: The Genesis of Einstein's work on the problem of motion

abstract: In this paper I describe the genesis of Einstein's early work on the problem of motion in general relativity (GR): the question of whether the motion of matter subject to gravity can be derived directly from the Einstein field equations. In addressing this question, Einstein himself always preferred the vacuum approach to the problem: the attempt to derive geodesic motion of matter from the vacuum Einstein equations. The paper first investigates why Einstein was so skeptical of the energy-momentum tensor and its role in GR. Drawing on hitherto unknown correspondence between Einstein and George Yuri Rainich, I then show step by step how his work on the vacuum approach came about, and how his quest for a unified field theory informed his interpretation of GR. I show that Einstein saw GR as a hybrid theory from very early on: fundamental and correct as far as gravity was concerned but phenomenological and effective in how it accounted for matter. As a result, Einstein saw energy-momentum tensors and singularities in GR as placeholders for a theory of matter not yet delivered. The reason he preferred singularities was that he hoped that their mathematical treatment would give a hint as to the sought after theory of matter, a theory that would do justice to quantum features of matter.


title: A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem

abstract: Stability and convergence of a time-weighted discrete scheme with nonuniform time steps are established for linear reaction-subdiffusion equations. The Caupto derivative is approximated at an offset point by using linear and quadratic polynomial interpolation. Our analysis relies on two tools: a discrete fractional Gr\"{o}nwall inequality and the global consistency analysis. The new consistency analysis makes use of an interpolation error formula for quadratic polynomials, which leads to a convolution-type bound for the local truncation error. To exploit these two tools, some theoretical properties of the discrete kernels in the numerical Caputo formula are crucial and we investigate them intensively in the nonuniform setting. Taking the initial singularity of the solution into account, we obtain a sharp error estimate on nonuniform time meshes. The fully discrete scheme generates a second-order accurate solution on the graded mesh provided a proper grading parameter is employed. An example is presented to show the sharpness of our analysis.


title: Bernstein Lethargy Theorem and Reflexivity

abstract: In this paper, we prove the equivalence of reflexive Banach spaces and those Banach spaces which satisfy the following form of Bernstein's Lethargy Theorem. Let $X$ be an arbitrary infinite-dimensional Banach space, and let the real-valued sequence $\{d_n\}_{n\ge1}$ decrease to $0$. Suppose that $\{Y_n\}_{n\ge1}$ is a system of strictly nested subspaces of $X$ such that $\overline Y_n \subset Y_{n+1}$ for all $n\ge1$ and for each $n\ge1$, there exists $y_n\in Y_{n+1}\backslash Y_n$ such that the distance $\rho(y_n,Y_n)$ from $y_n$ to the subspace $Y_n$ satisfies $$ \rho(y_n,Y_n)=\|y_n\|. $$ Then, there exists an element $x\in X$ such that $\rho(x,Y_n)=d_n$ for all $n\ge1$.


title: A Web Scraping Methodology for Bypassing Twitter API Restrictions

abstract: Retrieving information from social networks is the first and primordial step many data analysis fields such as Natural Language Processing, Sentiment Analysis and Machine Learning. Important data science tasks relay on historical data gathering for further predictive results. Most of the recent works use Twitter API, a public platform for collecting public streams of information, which allows querying chronological tweets for no more than three weeks old. In this paper, we present a new methodology for collecting historical tweets within any date range using web scraping techniques bypassing for Twitter API restrictions.


title: JSweep: A Patch-centric Data-driven Approach for Parallel Sweeps on Large-scale Meshes

abstract: In mesh-based numerical simulations, sweep is an important computation pattern. During sweeping a mesh, computations on cells are strictly ordered by data dependencies in given directions. Due to such a serial order, parallelizing sweep is challenging, especially for unstructured and deforming structured meshes. Meanwhile, recent high-fidelity multi-physics simulations of particle transport, including nuclear reactor and inertial confinement fusion, require {\em sweeps} on large scale meshes with billions of cells and hundreds of directions. In this paper, we present JSweep, a parallel data-driven computational framework integrated in the JAxMIN infrastructure. The essential of JSweep is a general patch-centric data-driven abstraction, coupled with a high performance runtime system leveraging hybrid parallelism of MPI+threads and achieving dynamic communication on contemporary multi-core clusters. Built on JSweep, we implement a representative data-driven algorithm, Sn transport, featuring optimizations of vertex clustering, multi-level priority strategy and patch-angle parallelism. Experimental evaluation with two real-world applications on structured and unstructured meshes respectively, demonstrates that JSweep can scale to tens of thousands of processor cores with reasonable parallel efficiency.


title: DRACO: Byzantine-resilient Distributed Training via Redundant Gradients

abstract: Distributed model training is vulnerable to byzantine system failures and adversarial compute nodes, i.e., nodes that use malicious updates to corrupt the global model stored at a parameter server (PS). To guarantee some form of robustness, recent work suggests using variants of the geometric median as an aggregation rule, in place of gradient averaging. Unfortunately, median-based rules can incur a prohibitive computational overhead in large-scale settings, and their convergence guarantees often require strong assumptions. In this work, we present DRACO, a scalable framework for robust distributed training that uses ideas from coding theory. In DRACO, each compute node evaluates redundant gradients that are used by the parameter server to eliminate the effects of adversarial updates. DRACO comes with problem-independent robustness guarantees, and the model that it trains is identical to the one trained in the adversary-free setup. We provide extensive experiments on real datasets and distributed setups across a variety of large-scale models, where we show that DRACO is several times, to orders of magnitude faster than median-based approaches.


title: Neck Detection for Two-Convex Hypersurfaces Embedded in Euclidean Space undergoing Brendle-Huisken G-Flow

abstract: Recently Brendle-Huisken introduced a fully nonlinear flow $G$. Their aim was to extend the surgery algorithm of Huisken-Sinestrari, into the Riemannian setting. The aim of this paper is to go through the details on how to perform neck detection for a closed, embedded hypersurface $M_0$ in $\mathbb{R}^{n+1}$ undergoing this $G$-flow. In order to do this we make some adjustments to Brendle and Huiskens gradient estimate, after we have done this we can go on to argue as in Huisken-Sinestrari's paper Mean curvature flow with surgeries of two-convex hyperusrfaces, in order to classify two-convex surfaces undergoing $G$-flow.


title: A discrete Gr\"{o}nwall inequality with application to numerical schemes for subdiffusion problems

abstract: We consider a class of numerical approximations to the Caputo fractional derivative. Our assumptions permit the use of nonuniform time steps, such as is appropriate for accurately resolving the behavior of a solution whose derivatives are singular at~$t=0$. The main result is a type of fractional Gr\"{o}nwall inequality and we illustrate its use by outlining some stability and convergence estimates of schemes for fractional reaction-subdiffusion problems. This approach extends earlier work that used the familiar L1 approximation to the Caputo fractional derivative, and will facilitate the analysis of higher order and linearized fast schemes.


title: k-Ary spanning trees contained in tournaments

abstract: A rooted tree is called a $k$-ary tree, if all non-leaf vertices have exactly $k$ children, except possibly one non-leaf vertex has at most $k-1$ children. Denote by $h(k)$ the minimum integer such that every tournament of order at least $h(k)$ contains a $k$-ary spanning tree. It is well-known that every tournament contains a Hamiltonian path, which implies that $h(1)=1$. Lu et al. [J. Graph Theory {\bf 30}(1999) 167--176] proved the existence of $h(k)$, and showed that $h(2)=4$ and $h(3)=8$. The exact values of $h(k)$ remain unknown for $k\geq 4$. A result of Erd\H{o}s on the domination number of tournaments implies $h(k)=\Omega(k\log k)$. In this paper, we prove that $h(4)=10$ and $h(5)\geq13$.


title: Magneto-Electric Effect in a Spin-State Transition System

abstract: Magnetic, dielectric, and magnetoelectric properties in a spin-state transition system are examined, motivated by the recent discovery of a multiferroic behavior in a cobalt oxide. We construct an effective model Hamiltonian based on the two-orbital Hubbard model, in which the spin-state degrees of freedom in magnetic ions couple with ferroelectric-type lattice distortions. A phase transition occurs from the high-temperature low-spin phase to the low-temperature high-spin ferroelectric phase with accompanying an increase of the spin entropy. The calculated results are consistent with the experimental pressure-temperature phase diagram. We predict the magnetic-field induced electric polarization in the low-spin paraelectric phase near the ferroelectric phase boundary.


title: Diversity Regularized Spatiotemporal Attention for Video-based Person Re-identification

abstract: Video-based person re-identification matches video clips of people across non-overlapping cameras. Most existing methods tackle this problem by encoding each video frame in its entirety and computing an aggregate representation across all frames. In practice, people are often partially occluded, which can corrupt the extracted features. Instead, we propose a new spatiotemporal attention model that automatically discovers a diverse set of distinctive body parts. This allows useful information to be extracted from all frames without succumbing to occlusions and misalignments. The network learns multiple spatial attention models and employs a diversity regularization term to ensure multiple models do not discover the same body part. Features extracted from local image regions are organized by spatial attention model and are combined using temporal attention. As a result, the network learns latent representations of the face, torso and other body parts using the best available image patches from the entire video sequence. Extensive evaluations on three datasets show that our framework outperforms the state-of-the-art approaches by large margins on multiple metrics.


title: Extremal weight projectors II

abstract: In previous work, we have constructed diagrammatic idempotents in an affine extension of the Temperley-Lieb category, which describe extremal weight projectors for sl(2), and which categorify Chebyshev polynomials of the first kind. In this paper, we generalize the construction of extremal weight projectors to the case of gl(N) for N greater than or equal to 2, with a view towards categorifying the corresponding torus skein algebras via Khovanov-Rozansky link homology. As by-products, we obtain compatible diagrammatic presentations of the representation categories of gl(N) and its Cartan subalgebra, and a categorification of power-sum symmetric polynomials.


title: Jittering jets in cooling flow clusters

abstract: I put forth a scenario of jet-driven jittering-jets in cooling flows in galaxies and clusters of galaxies. By uplifting cool gas from the center of the cluster jets and the bubbles they inflate cause the direction of later jets to substantially change. Some uplifted cool clumps fall back to feed the active galactic nucleus, in what is termed the cold feedback mechanism, but with an average angular momentum that is largely inclined to that of the original jets. The proposed jittering jets scenario allows, over the cooling time of the intracluster medium (ICM), the heating by the mixing of hot jet-inflated bubble gas with the ICM to operate in all directions. The proposed scenario further supports the heating by mixing mechanism of the ICM. The feedback heating mechanism of the ICM chiefly refers to the energy cycle. The cold feedback mechanism introduces mass as another ingredient of the feedback cycle, as the mass accretion of future activity cycles is determined by the formation of perturbations by the present activity cycle. The jittering jets scenario that I propose here brings angular momentum to be another ingredient of the feedback cycle.


title: Lolisa: Formal syntax and semantics for a subset of the solidity programming language in Mathematical Tool Coq

abstract: This article presents the formal syntax and semantics for a large subset of the Solidity programming language developed for the Etheruem blockchain platform based on our resent work about developing a general, extensible, and reusable formal memory (GERM) framework and an extension of Curry-Howard isomorphism, denoted as execution-verification isomorphism (EVI). This subset is denoted as Lolisa, which, to our knowledge, is the first mechanized and validated formal syntax and semantics developed for Solidity. The formal syntax of Lolisa adopts a stronger static type system than Solidity for enhanced type safety. In addition, Lolisa not only includes nearly all the syntax components of Solidity, such as mapping, modifier, contract, and address types, but it also contains general-purpose programming language features, such as multiple return values, pointer arithmetic, struct, and field access. Therefore, the inherent compatibility of Lolisa allows Solidity programs to be directly translated into Lolisa with a line-by-line correspondence without rebuilding or abstracting, and, in addition, the inherent generality of Lolisa allows it to be extended to express other programming languages as well. To this end, we also present a preliminary scheme for extending Lolisa to other languages systematically.


title: The Future of CISE Distributed Research Infrastructure

abstract: Shared research infrastructure that is globally distributed and widely accessible has been a hallmark of the networking community. This paper presents an initial snapshot of a vision for a possible future of mid-scale distributed research infrastructure aimed at enabling new types of research and discoveries. The paper is written from the perspective of "lessons learned" in constructing and operating the Global Environment for Network Innovations (GENI) infrastructure and attempts to project future concepts and solutions based on these lessons. The goal of this paper is to engage the community to contribute new ideas and to inform funding agencies about future research directions to realize this vision.


title: MLE-induced Likelihood for Markov Random Fields

abstract: Due to the intractable partition function, the exact likelihood function for a Markov random field (MRF), in many situations, can only be approximated. Major approximation approaches include pseudolikelihood and Laplace approximation. In this paper, we propose a novel way of approximating the likelihood function through first approximating the marginal likelihood functions of individual parameters and then reconstructing the joint likelihood function from these marginal likelihood functions. For approximating the marginal likelihood functions, we derive a particular likelihood function from a modified scenario of coin tossing which is useful for capturing how one parameter interacts with the remaining parameters in the likelihood function. For reconstructing the joint likelihood function, we use an appropriate copula to link up these marginal likelihood functions. Numerical investigation suggests the superior performance of our approach. Especially as the size of the MRF increases, both the numerical performance and the computational cost of our approach remain consistently satisfactory, whereas Laplace approximation deteriorates and pseudolikelihood becomes computationally unbearable.


title: Three-Dimensional Wave Behaviour of Light

abstract: We describe a simple experimental apparatus which allows one to observe the wave properties of light in a new way. This apparatus also makes possible to introduce and illustrate, in a very suggestive way, some fundamental principles of quantum theory.


title: Prediction of Stable Cu-Li Binary Intermetallics From First-Principles Calculations: Stoichiometries, Crystal Structures, and Physical Properties

abstract: Towards a resolution of the longstanding controversy regarding the existence of Cu-Li intermetallic compounds, we extensively investigate the phase stability of Cu-Li intermetallics with various possible stoichiometries at zero temperature and pressure using a global structure searching method. It is found that Cu-Li intermetallics can exist stably at atmospheric pressure, and three stable intermetallics (Fmmm Cu1Li2, Fd m Cu2Li1 and P Cu7Li1,) are identified. Electronic structure analysis reveals that although the three stable phases are metallic, covalent Cu-Cu and ionic Cu-Li bonds are found in the three structures. Moreover, the 3d states of copper atoms are mostly responsible for bond formations in the Cu-Li intermetallics. For all the predicted Cu-Li intermetallics, the effect of Cu concentration on structure, mechanical and thermodynamic properties are calculated systematically.


title: POKs Based Secure and Energy-Efficient Access Control for Implantable Medical Devices

abstract: Implantable medical devices (IMDs), such as pacemakers, implanted cardiac defibrillators, and neurostimulators are medical devices implanted into patients' bodies for monitoring physiological signals and performing medical treatments. Many IMDs have built-in wireless communication modules to facilitate data collecting and device reprogramming by external programmers. The wireless communication brings significant conveniences for advanced applications such as real-time and remote monitoring but also introduces the risk of unauthorized wireless access. The absence of effective access control mechanisms exposes patients' life to cyber attacks. In this paper, we present a lightweight and universally applicable access control system for IMDs. By leveraging Physically Obfuscated Keys (POKs) as the hardware root of trust, provable security is achieved based on standard cryptographic primitives while attaining high energy efficiency. In addition, barrier-free IMD access under emergent situations is realized by utilizing the patient's biometrical information. We evaluate our proposed scheme through extensive security analysis and a prototype implementation, which demonstrates our work's superiority on security and energy efficiency.


title: Cloud-based MPC with Encrypted Data

abstract: This paper explores the privacy of cloud outsourced Model Predictive Control (MPC) for a linear system with input constraints. In our cloud-based architecture, a client sends her private states to the cloud who performs the MPC computation and returns the control inputs. In order to guarantee that the cloud can perform this computation without obtaining anything about the client's private data, we employ a partially homomorphic cryptosystem. We propose protocols for two cloud-MPC architectures motivated by the current developments in the Internet of Things: a client-server architecture and a two-server architecture. In the first case, a control input for the system is privately computed by the cloud server, with the assistance of the client. In the second case, the control input is privately computed by two independent, non-colluding servers, with no additional requirements from the client. We prove that the proposed protocols preserve the privacy of the client's data and of the resulting control input. Furthermore, we compute bounds on the errors introduced by encryption. We present numerical simulations for the two architectures and discuss the trade-off between communication, MPC performance and privacy.


title: A {\it Chandra} Survey of Milky Way Globular Clusters I: Emissivity and Abundance of Weak X-ray Sources

abstract: Based on archival {\it Chandra} data, we have carried out an X-ray survey of 69, or nearly half the known population of, Milky Way globular clusters (GCs), focusing on weak X-ray sources, mainly cataclysmic variables (CVs) and coronally active binaries (ABs). Using the cumulative X-ray luminosity per unit stellar mass (i.e., X-ray emissivity) as a proxy of the source abundance, we demonstrate a paucity (lower by $41\%\pm27\%$ on average) of weak X-ray sources in most GCs relative to the field, which is represented by the Solar neighborhood and Local Group dwarf elliptical galaxies. We also revisit the mutual correlations among the cumulative X-ray luminosity ($L_X$), cluster mass ($M$) and stellar encounter rate ($\Gamma$), finding $L_{X}\propto M^{0.74\pm 0.13}$, $L_{X}\propto\Gamma^{0.67\pm0.07}$ and $\Gamma \propto M^{1.28 \pm 0.17}$. The three quantities can further be expressed as $L_{X} \propto M^{0.64\pm0.12}\ \Gamma^{0.19 \pm0.07}$, which indicates that the dynamical formation of CVs and ABs through stellar encounters in GCs is less dominant than previously suggested, and that the primordial formation channel has a substantial contribution. Taking these aspects together, we suggest that a large fraction of primordial, soft binaries have been disrupted in binary-single or binary-binary stellar interactions before they can otherwise evolve into X-ray-emitting close binaries, whereas the same interactions also have led to the formation of new close binaries. No significant correlations between $L_{X}/L_{K}$ and cluster properties, including dynamical age, metallicity and structural parameters, are found.


title: Constraint on ion-neutral drift velocity in the Class 0 protostar B335 from ALMA observations

abstract: Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal MHD effects proposed to enable the formation of Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare gas kinematics traced by the H13CO+ (3-2) and C18O (2-1) lines in the Class 0 protostar B335 obtained with our ALMA observations. A central compact (~1"-2") component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. We constructed kinematical models to fit the observed velocity structures and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. The infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85+/-0.2 km/s at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km/s at a radius of 100 au. The Hall parameter for H13CO+ is estimated to be >>1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between ions and neutrals could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0.1") are needed to examine this possibility and ambipolar diffusion on a smaller scale.


title: Multi-Scale Structure-Aware Network for Human Pose Estimation

abstract: We develop a robust multi-scale structure-aware neural network for human pose estimation. This method improves the recent deep conv-deconv hourglass models with four key improvements: (1) multi-scale supervision to strengthen contextual feature learning in matching body keypoints by combining feature heatmaps across scales, (2) multi-scale regression network at the end to globally optimize the structural matching of the multi-scale features, (3) structure-aware loss used in the intermediate supervision and at the regression to improve the matching of keypoints and respective neighbors to infer a higher-order matching configurations, and (4) a keypoint masking training scheme that can effectively fine-tune our network to robustly localize occluded keypoints via adjacent matches. Our method can effectively improve state-of-the-art pose estimation methods that suffer from difficulties in scale varieties, occlusions, and complex multi-person scenarios. This multi-scale supervision tightly integrates with the regression network to effectively (i) localize keypoints using the ensemble of multi-scale features, and (ii) infer global pose configuration by maximizing structural consistencies across multiple keypoints and scales. The keypoint masking training enhances these advantages to focus learning on hard occlusion samples. Our method achieves the leading position in the MPII challenge leaderboard among the state-of-the-art methods.


title: The fate of atomic spin in atomic scattering off surfaces

abstract: We explore model electron dynamics of an atom scattering off a surface within the time-dependent complete active space self consistent field (TD-CASSCF) approximation. We focus especially on the scattering of a hydrogen atom and its resulting spin-dynamics starting from an initially spin-polarized state. Our results reveal competing electronic time-scales that are governed by the electronic structure of the surface as well as the character of the atom. The timescales and nonadiabaticity of the dynamics are reported on by the final spin-polarization of the scattered atom, which may be probed in future experiments.


title: Linear perturbations of low angular momentum accretion flow in the Kerr metric and the corresponding emergent gravity phenomena

abstract: For certain geometric configuration of matter falling onto a rotating black hole, we develop a novel linear perturbation analysis scheme to perform the stability analysis of stationary integral accretion solutions corresponding to the steady state low angular momentum, inviscid, barotropic, irrotational, general relativistic accretion of hydrodynamic fluid. We demonstrate that such steady states remain stable under linear perturbation, and hence the stationary solutions are reliable to probe the black hole spacetime using the accretion phenomena.We report that a relativistic acoustic geometry emerges out as the consequence of such stability analysis procedure. We study various properties of that sonic geometry in detail. We construct the causal structures to establish the one to one correspondences of the sonic points with the acoustic black hole horizons, and the shock location with an acoustic white hole horizon. The influence of the spin of the rotating black holes on the emergence of such acoustic spacetime has been discussed.


title: The effects of kinetic instabilities on the electron cyclotron emission from runaway electrons

abstract: In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron-cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway electrons is obtained using a synthetic diagnostic model based on the reciprocity method. The electron distribution function is calculated using a kinetic simulation model including the whistler wave instabilities and the quasilinear diffusion effects. Simulations based on DIII-D low-density discharge reproduces the rapid growth of the ECE signals observed in DIII-D experiments. Unlike the thermal ECE where radiation for a certain frequency is strongly localized inside the resonance region, the non-thermal ECE radiation from runaway electrons is nonlocal, and the emission-absorption ratio is higher than that of thermal electrons. The runaway electron tail is more significant for ECE with higher frequencies, and the ECE spectrum becomes flatter as RE population grows. The nonlinear behavior of the kinetic instabilities is illustrated in the osculations of the ECE waves. The good agreement with the DIII-D experimental observations after including the kinetic instabilities clearly illustrate the significance of the scattering effects from wave-particle interactions, which can also be important for runaway electrons produced in disruptions.


title: On Fairness of Systemic Risk Measures

abstract: In our previous paper, "A Unified Approach to Systemic Risk Measures via Acceptance Set" (\textit{Mathematical Finance, 2018}), we have introduced a general class of systemic risk measures that allow for random allocations to individual banks before aggregation of their risks. In the present paper, we prove the dual representation of a particular subclass of such systemic risk measures and the existence and uniqueness of the optimal allocation related to them. We also introduce an associated utility maximization problem which has the same optimal solution as the systemic risk measure. In addition, the optimizer in the dual formulation provides a \textit{risk allocation} which is fair from the point of view of the individual financial institutions. The case with exponential utilities which allows for explicit computation is treated in details.


title: On conditions for weak conservativeness of regularized explicit finite-difference schemes for 1D barotropic gas dynamics equations

abstract: We consider explicit two-level three-point in space finite-difference schemes for solving 1D barotropic gas dynamics equations. The schemes are based on special quasi-gasdynamic and quasi-hydrodynamic regularizations of the system. We linearize the schemes on a constant solution and derive the von Neumann type necessary condition and a CFL type criterion (necessary and sufficient condition) for weak conservativeness in $L^2$ for the corresponding initial-value problem on the whole line. The criterion is essentially narrower than the necessary condition and wider than a sufficient one obtained recently in a particular case; moreover, it corresponds most well to numerical results for the original gas dynamics system.


title: On Difference-of-SOS and Difference-of-Convex-SOS Decompositions for Polynomials

abstract: In this article, we are interested in developing polynomial decomposition techniques based on sums-of-squares (SOS), namely the difference-of-sums-of-squares (D-SOS) and the difference-of-convex-sums-of-squares (DC-SOS). In particular, the DC-SOS decomposition is very useful for difference-of-convex (DC) programming formulation of polynomial optimization problems. First, we introduce the cone of convex-sums-of-squares (CSOS) polynomials and discuss its relationship to the sums-of-squares (SOS) polynomials, the non-negative polynomials and the SOS-convex polynomials. Then, we propose the set of D-SOS and DC-SOS polynomials, and prove that any polynomial can be formulated as D-SOS and DC-SOS. The problem of finding D-SOS and DC-SOS decompositions can be formulated as a semi-definite program and solved for any desired precision in polynomial time using interior point methods. Some algebraic properties of CSOS, D-SOS and DC-SOS are established. Second, we focus on establishing several practical algorithms for exact D-SOS and DC-SOS polynomial decompositions without solving any SDP. The numerical performance of the proposed D-SOS and DC-SOS decomposition algorithms and their parallel versions, tested on a dataset of 1750 randomly generated polynomials, is reported.


title: Mittens: An Extension of GloVe for Learning Domain-Specialized Representations

abstract: We present a simple extension of the GloVe representation learning model that begins with general-purpose representations and updates them based on data from a specialized domain. We show that the resulting representations can lead to faster learning and better results on a variety of tasks.


title: Unitary symmetry constraints on tensorial group field theory renormalization group flow

abstract: Renormalization group methods are an essential ingredient in the study of nonperturbative problems of quantum field theory. This paper deal with the symmetry constraints on the renormalization group flow for quartic melonic tensorial group field theories. Using the unitary invariance of the interactions, we provide a set of Ward-Takahashi identities which leads to relations between correlation functions. There are numerous reasons to consider such Ward identities in the functional renormalization group. Their compatibility along the flow provides a non-trivial constraint on the reliability of the approximation schemes used in the non-perturbative regime, especially on the truncation and the choice of the regulator. We establish the so called structure equations in the melonic sector and in the symmetric phase. As an example we consider the $T^4_5$ TGFT model without gauge constraint. The Wetterich flow equation is given and the way to improve the truncation on the effective action is also scrutinized.


title: Compassionately Conservative Balanced Cuts for Image Segmentation

abstract: The Normalized Cut (NCut) objective function, widely used in data clustering and image segmentation, quantifies the cost of graph partitioning in a way that biases clusters or segments that are balanced towards having lower values than unbalanced partitionings. However, this bias is so strong that it avoids any singleton partitions, even when vertices are very weakly connected to the rest of the graph. Motivated by the B\"uhler-Hein family of balanced cut costs, we propose the family of Compassionately Conservative Balanced (CCB) Cut costs, which are indexed by a parameter that can be used to strike a compromise between the desire to avoid too many singleton partitions and the notion that all partitions should be balanced. We show that CCB-Cut minimization can be relaxed into an orthogonally constrained $\ell_{\tau}$-minimization problem that coincides with the problem of computing Piecewise Flat Embeddings (PFE) for one particular index value, and we present an algorithm for solving the relaxed problem by iteratively minimizing a sequence of reweighted Rayleigh quotients (IRRQ). Using images from the BSDS500 database, we show that image segmentation based on CCB-Cut minimization provides better accuracy with respect to ground truth and greater variability in region size than NCut-based image segmentation.


title: Two loop electroweak corrections to $\bar B\rightarrow X_s\gamma$ and $B_s^0\rightarrow \mu^+\mu^-$ in the B-LSSM

abstract: The rare decays $\bar B\rightarrow X_s\gamma$ and $B_s^0\rightarrow \mu^+\mu^-$ are important to research new physics beyond standard model. In this work, we investigate two loop electroweak corrections to $\bar B\rightarrow X_s\gamma$ and $B_s^0\rightarrow \mu^+\mu^-$ in the minimal supersymmetric extension of the SM with local $B-L$ gauge symmetry (B-LSSM), under a minimal flavor violating assumption for the soft breaking terms. In this framework, new particles and new definition of squarks can affect the theoretical predictions of these two processes, with respect to the MSSM. Considering the constraints from updated experimental data, the numerical results show that the B-LSSM can fit the experimental data for the branching ratios of $\bar B\rightarrow X_s\gamma$ and $B_s^0\rightarrow \mu^+\mu^-$. The results of the rare decays also further constrain the parameter space of the B-LSSM.


title: Conformity in Scientific Networks

abstract: Scientists are generally subject to social pressures, including pressures to conform with others in their communities, that affect achievement of their epistemic goals. Here we analyze a network epistemology model in which agents, all else being equal, prefer to take actions that conform with those of their neighbors. This preference for conformity interacts with the agents' beliefs about which of two (or more) possible actions yields the better outcome. We find a range of possible outcomes, including stable polarization in belief and action. The model results are sensitive to network structure. In general, though, conformity has a negative effect on a community's ability to reach accurate consensus about the world.


title: Large-Scale Fabrication of RF MOSFETs on Liquid-Exfoliated MoS2

abstract: For the first time, thousands of RF MOSFETs were batch-fabricated on liquid-exfoliated MoS2 below 300 {\deg}C with nearly 100% yield. The large-scale fabrication with high yield allowed the average performance instead of the best performance to be reported. The DC performance of these devices were typical, but the RF performance, enabled by buried gates and on the order of 100 MHz, was reported for the first time for liquid-exfoliated MoS2. To resolve the dilemma of thin vs. thick films, gate recess was used on 20-nm thick films to improve the gate control while keeping the contact resistance lower than that on 10-nm films. These innovations may enable thin-film transistors to operate in the microwave range.


title: Reaction Mechanism of the Selective Reduction of CO$_2$ to CO by a Tetraaza [Co$^\text{II}$N$_4$H]$^{2+}$ Complex in the Presence of Protons

abstract: The tetraaza [Co$^\text{II}$N$_4$H]$^{2+}$ complex (\textbf{1}) is remarkable for its ability to selectively reduce CO$_2$ to CO with 45\% Faradaic efficiency and a CO to H$_2$ ratio of 3:2. We employ density functional theory (DFT) to determine the reasons behind the unusual catalytic properties of \textbf{1} and the most likely mechanism for CO$_2$ reduction. The selectivity for CO$_2$ over proton reduction is explained by analyzing the catalyst's affinity for the possible ligands present under typical reaction conditions: acetonitrile, water, CO$_2$, and bicarbonate. After reduction of the catalyst by two electrons, formation of [Co$^\text{I}$N$_4$H]$^{+}$-CO$_{2}^{-}$ is strongly favored. Based on thermodynamic and kinetic data, we establish that the only likely route for producing CO from here consists of a protonation step to yield [Co$^\text{I}$N$_4$H]$^{+}$-CO$_{2}$H, followed by reaction with CO$_2$ to form [Co$^\text{II}$N$_4$H]$^{2+}$-CO and bicarbonate. This conclusion corroborates the idea of a direct role of CO$_2$ as a Lewis acid to assist in {C-O} bond dissociation, a conjecture put forward by other authors to explain recent experimental observations. The pathway to formic acid is predicted to be forbidden by high activation barriers, in accordance with the products that are known to be generated by \textbf{1}. Calculated physical observables such as standard reduction potentials and the turnover frequency for our proposed catalytic cycle are in agreement with available experimental data reported in the literature. The mechanism also makes a prediction that may be experimentally verified: that the rate of CO formation should increase linearly with the partial pressure of CO$_2$.


title: Free hyperplane arrangements over arbitrary fields

abstract: In this paper, we study the class of free hyperplane arrangements. Specifically, we investigate the relations between freeness over a field of finite characteristic and freeness over $\mathbb{Q}$.


title: A Divide-and-Conquer Approach to Compressed Sensing MRI

abstract: Compressed sensing (CS) theory assures us that we can accurately reconstruct magnetic resonance images using fewer k-space measurements than the Nyquist sampling rate requires. In traditional CS-MRI inversion methods, the fact that the energy within the Fourier measurement domain is distributed non-uniformly is often neglected during reconstruction. As a result, more densely sampled low-frequency information tends to dominate penalization schemes for reconstructing MRI at the expense of high-frequency details. In this paper, we propose a new framework for CS-MRI inversion in which we decompose the observed k-space data into "subspaces" via sets of filters in a lossless way, and reconstruct the images in these various spaces individually using off-the-shelf algorithms. We then fuse the results to obtain the final reconstruction. In this way we are able to focus reconstruction on frequency information within the entire k-space more equally, preserving both high and low frequency details. We demonstrate that the proposed framework is competitive with state-of-the-art methods in CS-MRI in terms of quantitative performance, and often improves an algorithm's results qualitatively compared with it's direct application to k-space.


title: On the the successive passage times of certain one-dimensional diffusions

abstract: We show in detail some results, outlined in a previous paper regarding the case of Brownian motion (BM), about the distribution of the $n$th-passage time of a one-dimensional diffusion obtained by a space or time transformation of BM, through a constant barrier $a.$ Some explicit examples are reported.


title: Hypernatural inflation

abstract: We constructed a model of natural inflation in the context of $\alpha$-attractor supergravity, in which both the dilaton field and the axion field are light during inflation, and the inflaton may be a combination of the two. The T-model version of this theory is defined on the Poincare disk with radius |Z| = 1. It describes a Mexican hat potential with the flat axion direction corresponding to a circle of radius |Z| < 1. The axion decay constant $f_{a}$ in this theory can be exponentially large because of the hyperbolic geometry of the Poincare disk. Depending on initial conditions, this model may describe $\alpha$-attractor inflation driven by the radial component of the inflaton field, natural inflation driven by the axion field, or a sequence of these two regimes. We also construct the E-model version of this theory, which has similar properties. In addition, we describe generalized $\alpha$-attractor models where the potential can be singular at the boundary of the moduli space, and show that they can provide a simple solution for the problem of initial conditions for the models with plateau potentials.


title: Non-LTE Calculations of the Fe I 6173 {\AA} Line in a Flaring Atmosphere

abstract: The Fe I 6173 {\AA} line is widely used in the measurements of vector magnetic fields by instruments including the Helioseismic and Magnetic Imager (HMI). We perform non-local thermodynamic equilibrium calculations of this line based on radiative hydrodynamic simulations in a flaring atmosphere. We employ both a quiet-Sun atmosphere and a penumbral atmosphere as the initial one in our simulations. We find that, in the quiet-Sun atmosphere, the line center is obviously enhanced during an intermediate flare. The enhanced emission is contributed from both radiative backwarming in the photosphere and particle beam heating in the lower chromosphere. A blue asymmetry of the line profile also appears due to an upward mass motion in the lower chromosphere. If we take a penumbral atmosphere as the initial atmosphere, the line has a more significant response to the flare heating, showing a central emission and an obvious asymmetry. The low spectral resolution of HMI would indicate some loss of information but the enhancement and line asymmetry are still kept. By calculating polarized line profiles, we find that the Stokes I and V profiles can be altered as a result of flare heating. Thus the distortion of this line has a crucial influence on the magnetic field measured from this line, and one should be cautious in interpreting the magnetic transients observed frequently in solar flares.


title: Network Science approach to Modelling Emergence and Topological Robustness of Supply Networks: A Review and Perspective

abstract: Due to the increasingly complex and interconnected nature of global supply chain networks (SCNs), a recent strand of research has applied network science methods to model SCN growth and subsequently analyse various topological features, such as robustness. This paper provides: (1) a comprehensive review of the methodologies adopted in literature for modelling the topology and robustness of SCNs; (2) a summary of topological features of the real world SCNs, as reported in various data driven studies; and (3) a discussion on the limitations of existing network growth models to realistically represent the observed topological characteristics of SCNs. Finally, a novel perspective is proposed to mimic the SCN topologies reported in empirical studies, through fitness based generative network models.


title: The Role of Intracellular Interactions in the Collective Polarization of Tissues and its Interplay with Cellular Geometry

abstract: Planar cell polarity (PCP), the coherent in-plane polarization of a tissue on multicellular length scales, provides directional information that guides a multitude of developmental processes at cellular and tissue levels. While it is manifest that cells utilize both intracellular and intercellular mechanisms, how the two produce the collective polarization remains an active area of investigation. We study the role of intracellular interactions in the large-scale spatial coherence of cell polarities, and scrutinize the role of intracellular interactions in the emergence of tissue-wide polarization. We demonstrate that nonlocal cytoplasmic interactions are necessary and sufficient for the robust long-range polarization, and are essential to the faithful detection of weak directional signals. In the presence of nonlocal interactions, signatures of geometrical information in tissue polarity become manifest. We investigate the deleterious effects of geometric disorder, and determine conditions on the cytoplasmic interactions that guarantee the stability of polarization. These conditions get progressively more stringent upon increasing the geometric disorder. Another situation where the role of geometrical information might be evident is elongated tissues. Strikingly, our model recapitulates an observed influence of tissue elongation on the orientation of polarity. Eventually, we introduce three classes of mutants: lack of membrane proteins, cytoplasmic proteins, and local geometrical irregularities. We adopt core-PCP as a model pathway, and interpret the model parameters accordingly, through comparing the in silico and in vivo phenotypes. This comparison helps us shed light on the roles of the cytoplasmic proteins in cell-cell communication, and make predictions regarding the cooperation of cytoplasmic and membrane proteins in long-range polarization.


title: Evolution of Late-type Galaxies in a Cluster Environment: Effects of High-speed Multiple Encounters with Early-type Galaxies

abstract: Late-type galaxies falling into a cluster would evolve being influenced by the interactions with both the cluster and the nearby cluster member galaxies. Most numerical studies, however, tend to focus on the effects of the former with little work done on those of the latter. We thus perform a numerical study on the evolution of a late-type galaxy interacting with neighboring early-type galaxies at high speed, using hydrodynamic simulations. Based on the information obtained from the Coma cluster, we set up the simulations for the case where a Milky Way-like late-type galaxy experiences six consecutive collisions with twice as massive early-type galaxies having hot gas in their halos at the closest approach distances of 15-65 kpc/h at the relative velocities of 1500-1600 km/s. Our simulations show that the evolution of the late-type galaxy can be significantly affected by the accumulated effects of the high-speed multiple collisions with the early-type galaxies, such as on cold gas content and star formation activity of the late-type galaxy, particularly through the hydrodynamic interactions between cold disk and hot gas halos. We find that the late-type galaxy can lose most of its cold gas after the six collisions and have more star formation activity during the collisions. By comparing our simulation results with those of galaxy-cluster interactions, we claim that the role of the galaxy-galaxy interactions on the evolution of late-type galaxies in clusters could be comparable with that of the galaxy-cluster interactions, depending on the dynamical history.


title: Hard Photons in Hadroproduction of Top Quarks with Realistic Final States

abstract: We present a complete description of top quark pair production in association with a hard photon in the dilepton channel. Our calculation is accurate to NLO in QCD. It is based on matrix elements for $e^+\nu_e \mu^-\bar{\nu}_\mu b \bar{b}\gamma$ production and includes all resonant and non-resonant diagrams, interferences, and off-shell effects of the top quarks and the $W$ gauge bosons. This calculation constitutes the first full computation for top quark pair production with a final state photon in hadronic collisions at NLO in QCD. Numerical results for total and differential cross sections are presented for the LHC at a centre-of-mass energy of $\sqrt{s}=13$ TeV. For a few observables relevant for new physics searches, beyond some kinematic bounds, we observe shape distortions of more than $100\%$. In addition, we confirm that the size of the top quark off-shell effects for the total cross section is consistent with the expected uncertainties of the narrow width approximation. Results presented here are not only relevant for beyond the Standard Model physics searches but also important for precise measurements of the top-quark fiducial cross sections and top-quark properties at the LHC.


title: Electrically-activated spin-controlled orbital angular momentum multiplexer

abstract: We present and test the integration of a static orbital angular momentum mode multiplexer with a dynamical geometric-phase optical element enabling on-demand spin-controlled angular momentum multiplexing. A diffractive optics multiplexer fabricated with 3D high-resolution electron beam lithography performs a conformal mapping for the conversion from linear to azimuthal phase gradients. The latter is functionalized by a dynamic spin-orbit add-on that consists of a self-engineered electrically-activated liquid crystal optical vortex generator having large clear-aperture and high-resolution. By combining several functionalities based on the optical angular momentum of light in a compact manner, the proposed hybrid device could find applications in next-generation high-dimensional mode switchers and routers based on orbital angular momentum.


title: Reconfigurable Antenna Multiple Access for 5G mmWave Systems

abstract: This paper aims to realize a new multiple access technique based on recently proposed millimeter-wave reconfigurable antenna architectures. To this end, first we show that integration of the existing reconfigurable antenna systems with the well-known non-orthogonal multiple access (NOMA) technique causes a significant degradation in sum rate due to the inevitable power division in reconfigurable antennas. To circumvent this fundamental limit, a new multiple access technique is proposed. The technique which is called reconfigurable antenna multiple access (RAMA) transmits only each user's intended signal at the same time/frequency/code, which makes RAMA an inter-user interference-free technique. Two different cases are considered, i.e., RAMA with partial and full channel state information (CSI). In the first case, CSI is not required and only the direction of arrival for a specific user is used. Our analytical results indicate that with partial CSI and for symmetric channels, RAMA outperforms NOMA in terms of sum rate. Further, the analytical result indicates that RAMA for asymmetric channels achieves better sum rate than NOMA when less power is assigned to users that experience better channel quality. In the second case, RAMA with full CSI allocates optimal power to each user which leads to higher achievable rates compared to NOMA for both symmetric and asymmetric channels. The numerical computations demonstrate the analytical findings.


title: InSe: a two-dimensional material with strong interlayer coupling

abstract: Atomically thin, two-dimensional (2D) indium selenide (InSe) has attracted considerable attention due to large tunability in the band gap (from 1.4 to 2.6 eV) and high carrier mobility. The intriguingly high dependence of band gap on layer thickness may lead to novel device applications, although its origin remains poorly understood, and generally attributed to quantum confinement effect. In this work, we demonstrate via first-principles calculations that strong interlayer coupling may be mainly responsible for this phenomenon, especially in the fewer-layer region, and it could also be an essential factor influencing other material properties of {\beta}-InSe and {\gamma}-InSe. Existence of strong interlayer coupling manifests itself in three aspects: (i) indirect-to-direct band gap transitions with increasing layer thickness; (ii) fan-like frequency diagrams of the shear and breathing modes of few-layer flakes; (iii) strong layer-dependent carrier mobilities. Our results indicate that multiple-layer InSe may be deserving of attention from FET-based technologies and also an ideal system to study interlayer coupling, possibly inherent in other 2D materials.


title: Subamorphous thermal conductivity of crystalline half-Heusler superlattices

abstract: The quest to improve the thermoelectric figure of merit has mainly followed the roadmap of lowering the thermal conductivity while keeping unaltered the power factor of the material. Ideally an electron-crystal phonon-glass system is desired. In this work, we report an extraordinary reduction of the cross-plane thermal conductivity in crystalline (TiNiSn):(HfNiSn) half-Heusler superlattices. We create SLs with thermal conductivities below the effective amorphous limit, which is kept in a large temperature range (120-300 K). We measured thermal conductivity at room temperature values as low as 0.75 W/(m K), the lowest thermal conductivity value reported so far for half-Heusler compounds. By changing the deposition conditions, we also demonstrate that the thermal conductivity is highly impacted by the way the single segments of the superlattice grow. These findings show a huge potential for thermoelectric generators where an extraordinary reduction of the thermal conductivity is required but without losing the crystal quality of the system.


title: On 2-Absorbing Primary Hyperideals Of Multiplicative Hyperrings

abstract: Primary hyperideals have been introduced and studied in multiplicative hyperrings. In this paper, we intend to study extensively primary hyperideals of multiplicative hyperrings with absorbing zero and prove some results regarding them. Also, we describe Cu- ideals of multiplicative hyperrings which are particular classes of hyperideals. In the last section, we introduce 2-absorbing primary hyperideals and investigate the properties of this notion in commutative multiplicative hyperrings.


title: Horizontal vector fields and Seifert fiberings

abstract: This paper gives a classification of the topology of vector fields which are nowhere tangent to the fibers of a Seifert fibering.


title: Search for flavour-changing neutral current top-quark decays $t\to qZ$ in proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

abstract: A search for flavour-changing neutral-current processes in top-quark decays is presented. Data collected with the ATLAS detector from proton-proton collisions at the Large Hadron Collider at a centre-of-mass energy of $\sqrt{s}=13$ TeV, corresponding to an integrated luminosity of 36.1 fb$^{-1}$, are analysed. The search is performed using top-quark pair events, with one top quark decaying through the $t\rightarrow qZ$ $(q = u, c)$ flavour-changing neutral-current channel, and the other through the dominant Standard Model mode $t\rightarrow bW$. Only $Z$ boson decays into charged leptons and leptonic $W$ boson decays are considered as signal. Consequently, the final-state topology is characterized by the presence of three isolated charged leptons (electrons or muons) and at least two jets, one of the jets originating from a $b$-quark, and missing transverse momentum from the undetected neutrino. The data are consistent with Standard Model background contributions, and at 95% confidence level the search sets observed (expected) upper limits of $1.7\times10^{-4}$ ($2.4\times10^{-4}$) on the $t\to uZ$ branching ratio and $2.4\times10^{-4}$ ($3.2\times10^{-4}$) on the $t\to cZ$ branching ratio, constituting the most stringent limits to date.


title: New Calder\'on Reproducing Formulae with Exponential Decay on Spaces of Homogeneous Type

abstract: Assume that $(X, d, \mu)$ is a space of homogeneous type in the sense of Coifman and Weiss. In this article, motivated by the breakthrough work of P. Auscher and T. Hyt\"onen on orthonormal bases of regular wavelets on spaces of homogeneous type, the authors introduce a new kind of approximations of the identity with exponential decay (for short, $\exp$-ATI). Via such an $\exp$-ATI, motivated by another creative idea of Y. Han et al. to merge the aforementioned orthonormal bases of regular wavelets into the frame of the existed distributional theory on spaces of homogeneous type, the authors establish the homogeneous continuous/discrete Calder\'on reproducing formulae on $(X, d, \mu)$, as well as their inhomogeneous counterparts. The novelty of this article exists in that $d$ is only assumed to be a quasi-metric and the underlying measure $\mu$ a doubling measure, not necessary to satisfy the reverse doubling condition. It is well known that Calder\'on reproducing formulae are the cornerstone to develop analysis and, especially, harmonic analysis on spaces of homogeneous type.


title: Cancellable elements of the lattice of epigroup varieties

abstract: We completely determine all commutative epigroup varieties that are cancellable elements of the lattice EPI of all epigroup varieties. In particular, we verify that a commutative epigroup variety is a cancellable element of the lattice EPI if and only if it is a modular element of this lattice.


title: Diagonalwise Refactorization: An Efficient Training Method for Depthwise Convolutions

abstract: Depthwise convolutions provide significant performance benefits owing to the reduction in both parameters and mult-adds. However, training depthwise convolution layers with GPUs is slow in current deep learning frameworks because their implementations cannot fully utilize the GPU capacity. To address this problem, in this paper we present an efficient method (called diagonalwise refactorization) for accelerating the training of depthwise convolution layers. Our key idea is to rearrange the weight vectors of a depthwise convolution into a large diagonal weight matrix so as to convert the depthwise convolution into one single standard convolution, which is well supported by the cuDNN library that is highly-optimized for GPU computations. We have implemented our training method in five popular deep learning frameworks. Evaluation results show that our proposed method gains $15.4\times$ training speedup on Darknet, $8.4\times$ on Caffe, $5.4\times$ on PyTorch, $3.5\times$ on MXNet, and $1.4\times$ on TensorFlow, compared to their original implementations of depthwise convolutions.


title: A statistical mechanics approach to de-biasing and uncertainty estimation in LASSO for random measurements

abstract: In high-dimensional statistical inference in which the number of parameters to be estimated is larger than that of the holding data, regularized linear estimation techniques are widely used. These techniques have, however, some drawbacks. First, estimators are biased in the sense that their absolute values are shrunk toward zero because of the regularization effect. Second, their statistical properties are difficult to characterize as they are given as numerical solutions to certain optimization problems. In this manuscript, we tackle such problems concerning LASSO, which is a widely used method for sparse linear estimation, when the measurement matrix is regarded as a sample from a rotationally invariant ensemble. We develop a new computationally feasible scheme to construct a de-biased estimator with a confidence interval and conduct hypothesis testing for the null hypothesis that a certain parameter vanishes. It is numerically confirmed that the proposed method successfully de-biases the LASSO estimator and constructs confidence intervals and p-values by experiments for noisy linear measurements.


title: Entropy based Independent Learning in Anonymous Multi-Agent Settings

abstract: Efficient sequential matching of supply and demand is a problem of interest in many online to offline services. For instance, Uber, Lyft, Grab for matching taxis to customers; Ubereats, Deliveroo, FoodPanda etc for matching restaurants to customers. In these online to offline service problems, individuals who are responsible for supply (e.g., taxi drivers, delivery bikes or delivery van drivers) earn more by being at the "right" place at the "right" time. We are interested in developing approaches that learn to guide individuals to be in the "right" place at the "right" time (to maximize revenue) in the presence of other similar "learning" individuals and only local aggregated observation of other agents states (e.g., only number of other taxis in same zone as current agent). A key characteristic of the domains of interest is that the interactions between individuals are anonymous, i.e., the outcome of an interaction (competing for demand) is dependent only on the number and not on the identity of the agents. We model these problems using the Anonymous MARL (AyMARL) model. The key contribution of this paper is in employing principle of maximum entropy to provide a general framework of independent learning that is both empirically effective (even with only local aggregated information of agent population distribution) and theoretically justified. Finally, our approaches provide a significant improvement with respect to joint and individual revenue on a generic simulator for online to offline services and a real world taxi problem over existing approaches. More importantly, this is achieved while having the least variance in revenues earned by the learning individuals, an indicator of fairness.


title: Neutron background measurements with a hybrid neutron detector at the Kuo-Sheng Reactor Neutrino Laboratory

abstract: We report in situ neutron background measurements at the Kuo-Sheng Reactor Neutrino Laboratory (KSNL) by a hybrid neutron detector (HND) with a data size of 33.8 days under identical shielding configurations as during the neutrino physics data taking. The HND consists of BC-501A liquid and BC-702 phosphor powder scintillation neutron detectors, which is sensitive to both fast and thermal neutrons, respectively. Neutron-induced events for the two channels are identified and differentiated by pulse shape analysis, such that background of both are simultaneously measured. The fast neutron fluxes are derived by an iterative unfolding algorithm. Neutron induced background in the germanium detector under the same fluxes, both due to cosmic-rays and ambient radioactivity, are derived and compared with the measurements. The results are valuable to background understanding of the neutrino data at the KSNL. In particular, neutron-induced background events due to ambient radioactivity as well as from reactor operation are negligible compared to intrinsic cosmogenic activity and ambient $\gamma$-activity. The detector concept and analysis procedures are applicable to neutron background characterization in similar rare-event experiments.


title: Worst-Case Optimal Join Algorithms: Techniques, Results, and Open Problems

abstract: Worst-case optimal join algorithms are the class of join algorithms whose runtime match the worst-case output size of a given join query. While the first provably worst-case optimal join algorithm was discovered relatively recently, the techniques and results surrounding these algorithms grow out of decades of research from a wide range of areas, intimately connecting graph theory, algorithms, information theory, constraint satisfaction, database theory, and geometric inequalities. These ideas are not just paperware: in addition to academic project implementations, two variations of such algorithms are the work-horse join algorithms of commercial database and data analytics engines. This paper aims to be a brief introduction to the design and analysis of worst-case optimal join algorithms. We discuss the key techniques for proving runtime and output size bounds. We particularly focus on the fascinating connection between join algorithms and information theoretic inequalities, and the idea of how one can turn a proof into an algorithm. Finally, we conclude with a representative list of fundamental open problems in this area.


title: Classical $N$-Reflection Equation and Gaudin Models

abstract: We introduce the notion of $N$-reflection equation which provides a large generalization of the usual classical reflection equation describing integrable boundary conditions. The latter is recovered as a special example of the $N=2$ case. The basic theory is established and illustrated with several examples of solutions of the $N$-reflection equation associated to the rational and trigonometric $r$-matrices. A central result is the construction of a Poisson algebra associated to a non skew-symmetric $r$-matrix whose form is specified by a solution of the $N$-reflection equation. Generating functions of quantities in involution can be identified within this Poisson algebra. As an application, we construct new classical Gaudin-type Hamiltonians, particular cases of which are Gaudin Hamiltonians of $BC_L$-type .


title: Image Semantic Transformation: Faster, Lighter and Stronger

abstract: We propose Image-Semantic-Transformation-Reconstruction-Circle(ISTRC) model, a novel and powerful method using facenet's Euclidean latent space to understand the images. As the name suggests, ISTRC construct the circle, able to perfectly reconstruct images. One powerful Euclidean latent space embedded in ISTRC is FaceNet's last layer with the power of distinguishing and understanding images. Our model will reconstruct the images and manipulate Euclidean latent vectors to achieve semantic transformations and semantic images arthimetic calculations. In this paper, we show that ISTRC performs 10 high-level semantic transformations like "Male and female","add smile","open mouth", "deduct beard or add mustache", "bigger/smaller nose", "make older and younger", "bigger lips", "bigger eyes", "bigger/smaller mouths" and "more attractive". It just takes 3 hours(GTX 1080) to train the models of 10 semantic transformations.


title: Image-based deep learning for classification of noise transients in gravitational wave detectors

abstract: The detection of gravitational waves has inaugurated the era of gravitational astronomy and opened new avenues for the multimessenger study of cosmic sources. Thanks to their sensitivity, the Advanced LIGO and Advanced Virgo interferometers will probe a much larger volume of space and expand the capability of discovering new gravitational wave emitters. The characterization of these detectors is a primary task in order to recognize the main sources of noise and optimize the sensitivity of interferometers. Glitches are transient noise events that can impact the data quality of the interferometers and their classification is an important task for detector characterization. Deep learning techniques are a promising tool for the recognition and classification of glitches. We present a classification pipeline that exploits convolutional neural networks to classify glitches starting from their time-frequency evolution represented as images. We evaluated the classification accuracy on simulated glitches, showing that the proposed algorithm can automatically classify glitches on very fast timescales and with high accuracy, thus providing a promising tool for online detector characterization.


title: Variability of Red Supergiants in M31 from the Palomar Transient Factory

abstract: Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution phase when they are known to pulsate with varying amplitudes. The RSG period-luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory (PTF) survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K~ -10 mag (log(L/L_sun)>4.8) are variable at dm_R>0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M_K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first overtone hypothesis and indicates that the variable stars in this sample have 12 M_sun<M<24 M_sun. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors.


title: A Perfect Specialization Model for Gravity Equation in Bilateral Trade based on Production Structure

abstract: Although initially originated as a totally empirical relationship to explain the volume of trade between two partners, gravity equation has been the focus of several theoretic models that try to explain it. Specialization models are of great importance in providing a solid theoretic ground for gravity equation in bilateral trade. Some research papers try to improve specialization models by adding imperfect specialization to model, but we believe it is unnecessary complication. We provide a perfect specialization model based on the phenomenon we call tradability, which overcomes the problems with simpler initial. We provide empirical evidence using estimates on panel data of bilateral trade of 40 countries over 10 years that support the theoretical model. The empirical results have implied that tradability is the only reason for deviations of data from basic perfect specialization models.


title: Limit Behavior of Mass Critical Hartree Minimization Problems with Steep Potential Wells

abstract: We consider minimizers of the following mass critical Hartree minimization problem: \[ e_\lambda(N):=\underset{\{u\in H^1(R^d),\,\|u\|^2_2=N\}}{\inf} E_\lambda(u),\,\ d\ge 3, \] where the Hartree energy functional $E_\lambda(u)$ is defined by \[ E_\lambda(u):=\int_{R ^d}|\nabla u(x)|^2dx+\lambda \int_{R ^d}g(x)u^2(x)dx-\frac{1}{2} \int_{R ^d}\int_{R ^d} \frac{u^2(x)u^2(y)}{|x-y|^2}dxdy,\,\ \lambda>0,\] and the steep potential $g(x)$ satisfies $0=g(0)=\inf _{R^d}g(x)\le g(x)\le 1$ and $1-g(x)\in L^{\frac{d}{2}}(R^d)$. We prove that there exists a constant $N^*>0$, independent of $\lambda g(x)$, such that if $N\ge N^*$, then $e_\lambda(N)$ does not admit minimizers for any $\lambda >0$; if $0<N<N^*$, then there exists a constant $\lambda ^*(N)>0$ such that $e_\lambda(N)$ admits minimizers for any $\lambda >\lambda ^*(N)$, and $e_\lambda(N)$ does not admit minimizers for $0<\lambda <\lambda ^*(N)$. For any given $0<N<N^*$, the limit behavior of positive minimizers for $e_\lambda(N)$ is also studied as $\lambda\to\infty$, where the mass concentrates at the bottom of $g(x)$.


title: Dual Attention Matching Network for Context-Aware Feature Sequence based Person Re-Identification

abstract: Typical person re-identification (ReID) methods usually describe each pedestrian with a single feature vector and match them in a task-specific metric space. However, the methods based on a single feature vector are not sufficient enough to overcome visual ambiguity, which frequently occurs in real scenario. In this paper, we propose a novel end-to-end trainable framework, called Dual ATtention Matching network (DuATM), to learn context-aware feature sequences and perform attentive sequence comparison simultaneously. The core component of our DuATM framework is a dual attention mechanism, in which both intra-sequence and inter-sequence attention strategies are used for feature refinement and feature-pair alignment, respectively. Thus, detailed visual cues contained in the intermediate feature sequences can be automatically exploited and properly compared. We train the proposed DuATM network as a siamese network via a triplet loss assisted with a de-correlation loss and a cross-entropy loss. We conduct extensive experiments on both image and video based ReID benchmark datasets. Experimental results demonstrate the significant advantages of our approach compared to the state-of-the-art methods.


title: Directional Modulation: A Secure Solution to 5G and Beyond Mobile Networks

abstract: Directional modulation (DM), as an efficient secure transmission way, offers security through its directive property and is suitable for line-of-propagation (LoP) channels such as millimeter wave (mmWave) massive multiple-input multiple-output (MIMO), satellite communication, unmanned aerial vehicle (UAV), and smart transportation. If the direction angle of the desired received is known, the desired channel gain vector is obtainable. Thus, in advance, the DM transmitter knows the values of directional angles of desired user and eavesdropper, or their direction of arrival (DOAs) because the beamforming vector of confidential messages and artificial noise (AN) projection matrix is mainly determined by directional angles of desired user and eavesdropper. For a DM transceiver, working as a receiver, the first step is to measure the DOAs of desired user and eavesdropper. Then, in the second step, using the measured DOAs, the beamforming vector of confidential messages and AN projection matrix is designed. In this paper, we describe the DOA measurement methods, power allocation, and beamforming in DM networks. A machine learning-based DOA measurement method is proposed to make a substantial SR performance gain compared to single-snapshot measurement without machine learning for a given null-space projection beamforming scheme. However, for a conventional DM network, there still exists a serious secure issue: the eavesdropper moves inside the main beam of the desired user and may intercept the confidential messages intended to the desired users because the beamforming vector of confidential messages and AN projection matrix are only angle-dependence. To address this problem, we present a new concept of secure and precise transmission, where the transmit waveform has two-dimensional even three-dimensional dependence by using DM, random frequency selection, and phase alignment at DM transmitter.


title: An Empirical Study of Fault Localization Families and Their Combinations

abstract: The performance of fault localization techniques is critical to their adoption in practice. This paper reports on an empirical study of a wide range of fault localization techniques on real-world faults. Different from previous studies, this paper (1) considers a wide range of techniques from different families, (2) combines different techniques, and (3) considers the execution time of different techniques. Our results reveal that a combined technique significantly outperforms any individual technique (200% increase in faults localized in Top 1), suggesting that combination may be a desirable way to apply fault localization techniques and that future techniques should also be evaluated in the combined setting. Our implementation is publicly available for evaluating and combining fault localization techniques.


title: Weakly Consistent Extensions of Lower Previsions

abstract: Several consistency notions are available for a lower prevision P assessed on a set D of gambles (bounded random variables), ranging from the well known coherence to convexity and to the recently introduced 2-coherence and 2-convexity. In all these instances, a procedure with remarkable features, called (coherent, convex, 2-coherent or 2-convex) natural extension, is available to extend P, preserving its consistency properties, to an arbitrary superset of gambles. We analyse the 2-coherent and 2-convex natural extensions, E_2 and E_2c respectively, showing that they may coincide with the other extensions in certain, special but rather common, cases of `full' conditional lower prevision or probability assessments. This does generally not happen if P is a(n unconditional) lower probability on the powerset of a given partition and is extended to the gambles defined on the same partition. In this framework we determine alternative formulae for E_2 and E_2c. We also show that E_2c may be nearly vacuous in some sense, while the Choquet integral extension is 2-coherent if P is, and bounds from above the 2-coherent natural extension. Relationships between the finiteness of the various natural extensions and conditions of avoiding sure loss or weaker are also pointed out.


title: Iteration complexity of first-order augmented Lagrangian methods for convex conic programming

abstract: In this paper we consider a class of convex conic programming. In particular, we first propose an inexact augmented Lagrangian (I-AL) method that resembles the classical I-AL method for solving this problem, in which the augmented Lagrangian subproblems are solved approximately by a variant of Nesterov's optimal first-order method. We show that the total number of first-order iterations of the proposed I-AL method for finding an $\epsilon$-KKT solution is at most $\mathcal{O}(\epsilon^{-7/4})$. We then propose an adaptively regularized I-AL method and show that it achieves a first-order iteration complexity $\mathcal{O}(\epsilon^{-1}\log\epsilon^{-1})$, which significantly improves existing complexity bounds achieved by first-order I-AL methods for finding an $\epsilon$-KKT solution. Our complexity analysis of the I-AL methods is based on a sharp analysis of inexact proximal point algorithm (PPA) and the connection between the I-AL methods and inexact PPA. It is vastly different from existing complexity analyses of the first-order I-AL methods in the literature, which typically regard the I-AL methods as an inexact dual gradient method.


title: Oxygen and zinc abundances in 417 Galactic bulge red giants

abstract: Oxygen and zinc in the Galactic bulge are key elements for the understanding of the bulge chemical evolution. Oxygen-to-iron abundance ratios provide a most robust indicator of the star formation rate and chemical evolution of the bulge. Zinc is enhanced in metal-poor stars, behaving as an $\alpha$-element, and its production may require nucleosynthesis in hypernovae. Most of the neutral gas at high redshift is in damped Lyman-alpha systems (DLAs), where Zn is also observed to behave as an alpha-element. The aim of this work is the derivation of the alpha-element oxygen, together with nitrogen, and the iron-peak element zinc abundances in 417 bulge giants, from moderate resolution (R~22,000) FLAMES-GIRAFFE spectra. For stars in common with a set of UVES spectra with higher resolution (R~45,000), the data are intercompared. The results are compared with literature data and chemodynamical models.


title: An efficient stochastic algorithm for the perturbative density matrix renormalization group in large active spaces

abstract: We present an efficient stochastic algorithm for the recently introduced perturbative density matrix renormalization group (p-DMRG) method for large active spaces. The stochastic implementation bypasses the computational bottleneck involved in solving the first order equation in the earlier deterministic algorithm. We demonstrate the efficiency and accuracy of the algorithm on the C$_2$ and Cr$_2$ molecular benchmark systems.


title: Tricolor Technique for Visualization of Spatial Variations of Polydisperse Dust in Gas-Dust Flows

abstract: The aim of this work is to construct an algorithm for visualizing a polydisperse phase of solid particles (dust) in an inhomogeneous flow of a two-phase gas-dust mixture that would allow us to see, within one plot, the degree of polydispersity of the dust phase and the difference in the spatial distributions of individual fractions of dust particles in the computational domain. The developed technique allows us to reproduce concentrations from one to three fractions of dust particles in each cell in the computational domain. Each of the three fractions of dust particles is mapped to one of the main channels of the RGB palette. The intensity of the color shade is set to be proportional to the relative concentration of dust particles in this fraction. The final image for a polydisperse mixture is obtained by adding images in each of the three color channels. To visualize the degree of polydispersity, I propose depicting the spatial distribution of the entropy of the dust mixture. The definition of the entropy of a mixture is generalized to take into account the states of a mixture with zero number of particles in the mixture. They correspond to dust-free sections of the computational domain (voids). The proposed method for visualizing the polydispersity of a mixture of particles is demonstrated using the example of dynamic numerical modeling of the spatial features of dust structures formed in turbulent gas-dust flows and in flows with shock waves.


title: A Model of Fast Radio Bursts: Collisions Between Episodic Magnetic Blobs

abstract: Fast radio bursts (FRBs) are bright radio pulses from the sky with millisecond durations and Jansky-level flux densities. Their origins are still largely uncertain. Here we suggest a new model for FRBs. We argue that the collision of a white dwarf with a black hole can give birth to a transient accretion disk, from which powerful episodic magnetic blobs will be launched. The collision between two consecutive magnetic blobs can result in a catastrophic magnetic reconnection, which releases a large amount of free magnetic energy and forms a forward shock. The shock propagates through the cold magnetized plasma within the blob in the collision region, radiating through synchrotron maser mechanism, which is responsible for a non-repeating FRB signal. Our calculations show that the theoretical energetics, radiation frequency, duration timescale, and event rate can be well consistent with the observational characteristic of FRBs.