title: Supervised detection of exoplanets in high-contrast imaging sequences

abstract: Post-processing algorithms play a key role in pushing the detection limits of high-contrast imaging (HCI) instruments. State-of-the-art image processing approaches for HCI enable the production of science-ready images relying on unsupervised learning techniques, such as low-rank approximations, for generating a model PSF and subtracting the residual starlight and speckle noise. In order to maximize the detection rate of HCI instruments and survey campaigns, advanced algorithms with higher sensitivities to faint companions are needed, especially for the speckle-dominated innermost region of the images. We propose a reformulation of the exoplanet detection task (for ADI sequences) that builds on well-established machine learning techniques to take HCI post-processing from an unsupervised to a supervised learning context. In this new framework, we present algorithmic solutions using two different discriminative models: SODIRF (random forests) and SODINN (neural networks). We test these algorithms on real ADI datasets from VLT/NACO and VLT/SPHERE HCI instruments. We then assess their performances by injecting fake companions and using receiver operating characteristic analysis. This is done in comparison with state-of-the-art ADI algorithms, such as ADI-PCA. This study shows the improved sensitivity vs specificity trade-off of the proposed approach. At the diffraction limit, SODINN improves the true positive rate by a factor ranging from ~2 to ~10 (depending on the dataset and angular separation) with respect to ADI-PCA when working at the same false positive level. The proposed supervised detection framework outperforms state-of-the-art techniques in the task of discriminating planet signal from speckles. In addition, it offers the possibility of re-processing existing HCI databases to maximize their scientific return and potentially improve the demographics of directly imaged exoplanets.


title: Simulation of Scalar Field Theories with Complex Actions

abstract: We develop a method for the simulation of scalar field theories with complex actions which is local, simple to implement and can be used in any number of space-time dimensions. For model systems satisfying the $\mathcal{PT}$ symmetry condition $L^{*}(\phi)=L(-\phi)$, the complex weight problem is reduced to a sign problem. The sign problem is eliminated completely for a large subclass of these models; this class includes models within the $i\phi^{3}$ universality class, and also models with nonzero chemical potential. Simulations of models from this subclass show a rich set of behaviors. Propagators may exhibit damped oscillations, indicating a clear violation of spectral positivity. Modulated phases occur in some models, exhibiting striping and other pattern-forming behaviors. These field theory models are connected to complex systems where pattern formation occurs because of competition between interactions at two different length scales.


title: Theoretical study of asymmetric superrotors: alignment and orientation

abstract: We report a theoretical study of the optical centrifuge acceleration of an asymmetric top molecule interacting with an electric static field by solving the time-dependent Schr\"odinger equation in the rigid rotor approximation. A detailed analysis of the mixing of the angular momentum in both the molecular and the laboratory fixed frames allow us to deepen the understanding of the main features of the acceleration process, for instance, the effective angular frequency of the molecule at the end of the pulse. In addition, we prove numerically that the asymmetric superrotors rotate around one internal axis and that their dynamics is confined to the plane defined by the polarization axis of the laser, in agreement with experimental findings. Furthermore, we consider the orientation patterns induced by the dc field, showing the characteristics of their structure as a function of the strength of the static field and the initial configuration of the fields.


title: Local nets of von Neumann algebras in the Sine-Gordon model

abstract: The Haag-Kastler net of local von Neumann algebras is constructed in the ultraviolet finite regime of the sine-Gordon model, and its equivalence with the massive Thirring model is proved. In contrast to other authors, we do not add an auxiliary mass term, and we work completely in Lorentzian signature. The construction is based on the functional formalism for perturbative Algebraic Quantum Field Theory together with estimates originally derived within Constructive Quantum Field Theory and adapted to Lorentzian signature. The paper extends previous work by two of us.


title: The Shrinkage Variance Hotelling $T^2$ Test for Genomic Profiling Studies

abstract: Designed gene expression micro-array experiments, consisting of several treatment levels with a number of replicates per level, are analyzed by applying simple tests for group differences at the per gene level. The gene level statistics are sorted and a criterion for selecting important genes which takes into account multiplicity is applied. A caveat arises in that true signals (genes truly over or under expressed) are "competing" with fairly large type I error signals. False positives near the top of a sorted list can occur when genes having very small fold-change are compensated by small enough variance to yield a large test statistic. One of the first attempts around this caveat as the development of "significance analysis of micro-arrays (SAM)", which used a modified t-type statistic thresholded against its permutation distribution. The key innovation of the modified t-statistic was the addition of a constant to the per gene standard errors in order to stabilize the coefficient of variation of the resulting test statistic. Since then, several authors have proposed the use of shrinkage variance estimators in conjunction with t-type, and more generally, ANOVA type tests at the gene level. Our new approach proposes the use of a shrinkage variance Hotelling T-squared statistic in which the per gene sample covariance matrix is replaced by a shrinkage estimate borrowing strength from across all genes. It is demonstrated that the new statistic retains the F-distribution under the null, with added degrees of freedom in the denominator. Advantages of this class of tests are (i) flexibility in that a whole family of hypothesis tests is possible (ii) the gains of the above-mentioned earlier innovations are enjoyed more fully. This paper summarizes our results and presents a simulation study benchmarking the new statistic against another recently proposed statistic.


title: Resonance phenomena in the $\varphi^8$ kinks scattering

abstract: We study the scattering of the $\varphi^8$ kinks with power-law asymptotics. We found two critical values of the initial velocity, $v_{cr}^{(1)}$ and $v_{cr}^{(2)}$, which separate different regimes of the kink-antikink collision. At the initial velocities $v_{in}< v_{cr}^{(1)}$ kinks can not collide due to repulsive force between them. At $v_{in}>v_{cr}^{(2)}$ the kinks escape to spatial infinities after one collision. In the range $v_{cr}^{(1)}\le v_{in}\le v_{cr}^{(2)}$ we observed kinks capture and formation of their bound state. Besides that, at these initial velocities we found resonance phenomena -- escape windows.

categories: hep-th cond-mat.mtrl-sci math-ph math.MP nlin.PS


title: Stability of Control Systems with Feedback from Energy Harvesting Sensors

abstract: In this paper, we study the problem of certifying the stability of a closed loop system which receives feedback from an energy harvesting sensor. This is important, as energy harvesting sensors are recharged stochastically, and may only be able to provide feedback intermittently. Thus, stabilizing plants with feedback provided by energy harvesting sensors is challenging in that the feedback signal is only available stochastically, complicating the analysis of the closed-loop system. As the main contribution of the paper, we show that for a broad class of energy harvesting processes and transmission policies, the system can be modeled as a Markov jump linear system (MJLS), which thereby enables a rigorous stability analysis. We discuss the types of transmission policies and energy harvesting processes which can be accommodated in detail, demonstrating the generality of the results.


title: Strong convergence of quantum random walks via semigroup decomposition

abstract: We give a simple and direct treatment of the strong convergence of quantum random walks to quantum stochastic operator cocycles, via the semigroup decomposition of such cocycles. Our approach also delivers convergence of the pointwise product of quantum random walks to the quantum stochastic Trotter product of the respective limit cocycles, thereby revealing the algebraic structure of the limiting procedure. The repeated quantum interactions model is shown to fit nicely into the convergence scheme described.


title: Sketched Clustering via Hybrid Approximate Message Passing

abstract: In sketched clustering, a dataset of $T$ samples is first sketched down to a vector of modest size, from which the centroids are subsequently extracted. Advantages include i) reduced storage complexity and ii) centroid extraction complexity independent of $T$. For the sketching methodology recently proposed by Keriven, et al., which can be interpreted as a random sampling of the empirical characteristic function, we propose a sketched clustering algorithm based on approximate message passing. Numerical experiments suggest that our approach is more efficient than the state-of-the-art sketched clustering algorithm "CL-OMPR" (in both computational and sample complexity) and more efficient than k-means++ when $T$ is large.


title: t-Private Information Retrieval Schemes Using Transitive Codes

abstract: This paper presents private information retrieval (PIR) schemes for coded storage with colluding servers, which are not restricted to maximum distance separable (MDS) codes. PIR schemes for general linear codes are constructed and the resulting PIR rate is calculated explicitly. It is shown that codes with transitive automorphism groups yield the highest possible rates obtainable with the proposed scheme. This rate coincides with the known asymptotic PIR capacity for MDS-coded storage systems without collusion. While many PIR schemes in the literature require field sizes that grow with the number of servers and files in the system, we focus especially on the case of a binary base field, for which Reed- Muller codes serve as an important and explicit class of examples.


title: Assessing the capacity and flow of ecosystem services in multifunctional landscapes: evidence of a rural-urban gradient in a Mediterranean small island state

abstract: Distinguishing between the capacity of ecosystems to generate ecosystem services (ES) and the actual use of these service (ES flow) in ES assessment and mapping is important to develop an understanding of the sustainability of ES use. This study assesses the spatial variation in ES capacity and flow in the Mediterranean small island state of Malta. The services included in this study were crop provisioning, beekeeping and honey production, fodder and livestock production, crop pollination, air quality regulation, and aesthetic ES. This assessment develops different spatial models, which make use of available datasets, causal relationships between datasets, including a generated land use land cover (LULC) map, and statistical models and indicators based on direct measurements. Individual ES indicators were mapped to visualise and compare their spatial patterns across the case study area. Subsequently, an analysis of ES associations and bundles was carried out using Pearson parametric correlation test, for both ES capacity and flow indicators generated from this study, and through Principal Component Analysis. Results demonstrate several significant synergistic interactions between ES capacity and flow in rural landscapes characterised with agricultural and semi-natural LULC categories, indicating high landscape multifunctionality. In contrast, predominantly urban areas tend to be characterised with a low ecosystem capacity and ES flow, suggesting that ES delivery in the landscapes of the study area is determined by land use intensity. These findings support the notion that multifunctional rural landscapes provide multiple ES, making an important contribution to human well-being, and that land use planning that develops green infrastructure in urban areas can significantly contribute to support biodiversity and ES delivery.


title: Exponential Stability and Supporting Spectral Analysis of a Linearized Compressible Flow-Structure PDE Model

abstract: In this work, a result of exponential stability is obtained for solutions of a compressible flow-structure partial differential equation (PDE) model which has recently appeared in the literature. In particular, a compressible flow PDE and its associated state equation for the associated pressure variable, each evolving within a three dimensional domain $\mathcal{O}$, are coupled to a fourth order plate equation which holds on a flat portion $% \Omega $ of the boundary $\partial \mathcal{O}$. Moreover, since this coupled PDE model is the result of a linearization of the compressible Navier-Stokes equations about an arbitrary state, the flow PDE component contains a generally nonzero ambient flow profile $\mathbf{U}$. By way of obtaining the aforesaid exponential stability, a \textquotedblleft frequency domain\textquotedblright\ approach is adopted here, an approach which is predicated on obtaining a uniform estimate on the resolvent of the associated flow-structure semigroup generator.


title: Orphan spins and bound in-gap states in the $S={5\over2}$ antiferromagnet CaFe$_{2}$O$_{4}$

abstract: CaFe$_{2}$O$_{4}$ is an anisotropic $S={5\over 2}$ antiferromagnet with two competing $A$ ($\uparrow \uparrow \downarrow \downarrow$) and $B$ ($\uparrow \downarrow \uparrow \downarrow$) magnetic order parameters separated by static antiphase boundaries at low temperatures. Neutron diffraction and bulk susceptibility measurements, show that the spins near these boundaries are weakly correlated and a carry an uncompensated ferromagnetic moment that can be tuned with a magnetic field. Spectroscopic measurements find these spins are bound with excitation energies less than the bulk magnetic spin-waves and resemble the spectra from isolated spin-clusters. Localized bound orphaned spins separate the two competing magnetic order parameters in CaFe$_{2}$O$_{4}$.


title: Stochastic reconstruction of an oolitic limestone by generative adversarial networks

abstract: Stochastic image reconstruction is a key part of modern digital rock physics and materials analysis that aims to create numerous representative samples of material micro-structures for upscaling, numerical computation of effective properties and uncertainty quantification. We present a method of three-dimensional stochastic image reconstruction based on generative adversarial neural networks (GANs). GANs represent a framework of unsupervised learning methods that require no a priori inference of the probability distribution associated with the training data. Using a fully convolutional neural network allows fast sampling of large volumetric images.We apply a GAN based workflow of network training and image generation to an oolitic Ketton limestone micro-CT dataset. Minkowski functionals, effective permeability as well as velocity distributions of simulated flow within the acquired images are compared with the synthetic reconstructions generated by the deep neural network. While our results show that GANs allow a fast and accurate reconstruction of the evaluated image dataset, we address a number of open questions and challenges involved in the evaluation of generative network-based methods.


title: Stability of the Couette flow at high Reynolds number in 2D and 3D

abstract: We review works on the asymptotic stability of the Couette flow. The majority of the paper is aimed towards a wide range of applied mathematicians with an additional section aimed towards experts in the mathematical analysis of PDEs.


title: Effective Neural Solution for Multi-Criteria Word Segmentation

abstract: We present a simple yet elegant solution to train a single joint model on multi-criteria corpora for Chinese Word Segmentation (CWS). Our novel design requires no private layers in model architecture, instead, introduces two artificial tokens at the beginning and ending of input sentence to specify the required target criteria. The rest of the model including Long Short-Term Memory (LSTM) layer and Conditional Random Fields (CRFs) layer remains unchanged and is shared across all datasets, keeping the size of parameter collection minimal and constant. On Bakeoff 2005 and Bakeoff 2008 datasets, our innovative design has surpassed both single-criterion and multi-criteria state-of-the-art learning results. To the best knowledge, our design is the first one that has achieved the latest high performance on such large scale datasets. Source codes and corpora of this paper are available on GitHub.


title: On 10D SYM superamplitudes

abstract: Recently the spinor helicity and (two types of) superamplitude formalisms for 11D supergravity and 10D supersymmetric Yang-Mills theories were proposed in [1,2,3]. In this contribution we describe briefly the basic properties of these superamplitudes for the simpler case of 10D SYM.


title: Computing the Photometric and Polarimetric Variability of Be Stars

abstract: We investigate variations in the linear polarization as well as in the V-band and B-band colour-magnitudes for classical Be star disks. We present two models: disks with enhanced disk density and disks that are tilted or warped from the stellar equatorial plane. In both cases, we predict variation in observable properties of the system as the disk rotates. We use a non-LTE radiative transfer code BEDISK (Sigut & Jones) in combination with a Monte Carlo routine that includes multiple scattering (Halonen et al.) to model classical Be star systems. We find that a disk with an enhanced density region that is one order of magnitude denser than the disk's base density shows as much as $\sim$0.2% variability in the polarization while the polarization position angle varies by $\sim$8\deg. The $\Delta$V magnitude for the same system shows variations of up to $\sim$0.4 magnitude while the $\Delta$(B-V) colour varies by at most $\sim$0.01 magnitude. We find that disks tilted from the equatorial plane at small angles of $\sim$30\deg more strongly reflect the values of polarization and colour-magnitudes reported in the literature than disks tilted at larger angles. For this model, the linear polarization varies by $\sim$0.3%, the polarization position angle varies by $\sim$60\deg, the $\Delta$V magnitude varies up to 0.35 magnitude, and the $\Delta$(B-V) colour varies up to 0.1 magnitude. We find that the enhanced disk density models show ranges of polarization and colour-magnitudes that are commensurate with what is reported in the literature for all sizes of the density enhanced regions. From this, we cannot determine any preference for small or large density enhanced regions.


title: Self-supervised Multi-level Face Model Learning for Monocular Reconstruction at over 250 Hz

abstract: The reconstruction of dense 3D models of face geometry and appearance from a single image is highly challenging and ill-posed. To constrain the problem, many approaches rely on strong priors, such as parametric face models learned from limited 3D scan data. However, prior models restrict generalization of the true diversity in facial geometry, skin reflectance and illumination. To alleviate this problem, we present the first approach that jointly learns 1) a regressor for face shape, expression, reflectance and illumination on the basis of 2) a concurrently learned parametric face model. Our multi-level face model combines the advantage of 3D Morphable Models for regularization with the out-of-space generalization of a learned corrective space. We train end-to-end on in-the-wild images without dense annotations by fusing a convolutional encoder with a differentiable expert-designed renderer and a self-supervised training loss, both defined at multiple detail levels. Our approach compares favorably to the state-of-the-art in terms of reconstruction quality, better generalizes to real world faces, and runs at over 250 Hz.


title: Remarks on Bayesian Control Charts

abstract: There is a considerable amount of ongoing research on the use of Bayesian control charts for detecting a shift from a good quality distribution to a bad quality distribution in univariate and multivariate processes. It is widely claimed that Bayesian control charts are economically optimal; see, for example, Calabrese (1995) [Bayesian process control for attributes. Management Science, DOI: 10.1287/mnsc.41.4.637] and Makis (2008) [Multivariate Bayesian control chart. Operations Research, DOI: 10.1287/opre.1070.0495]. Some researchers also generalize the optimality of controls defined based on posterior probabilities to the class of partially observable Markov decision processes. This note points out that the existing Bayesian control charts cannot generally be optimal because many years ago an analytical counterexample was provided by Taylor (1965) [Markovian sequential replacement processes. The Annals of Mathematical Statistics, DOI: 10.1214/aoms/1177699796].


title: Per-Pixel Feedback for improving Semantic Segmentation

abstract: Semantic segmentation is the task of assigning a label to each pixel in the image.In recent years, deep convolutional neural networks have been driving advances in multiple tasks related to cognition. Although, DCNNs have resulted in unprecedented visual recognition performances, they offer little transparency. To understand how DCNN based models work at the task of semantic segmentation, we try to analyze the DCNN models in semantic segmentation. We try to find the importance of global image information for labeling pixels. Based on the experiments on discriminative regions, and modeling of fixations, we propose a set of new training loss functions for fine-tuning DCNN based models. The proposed training regime has shown improvement in performance of DeepLab Large FOV(VGG-16) Segmentation model for PASCAL VOC 2012 dataset. However, further test remains to conclusively evaluate the benefits due to the proposed loss functions across models, and data-sets. Submitted in part fulfillment of the requirements for the degree of Integrated Masters of Science in Applied Mathematics. Update: Further Experiment showed minimal benefits. Code Available [here](https://github.com/BardOfCodes/Seg-Unravel).


title: MoDL: Model Based Deep Learning Architecture for Inverse Problems

abstract: We introduce a model-based image reconstruction framework with a convolution neural network (CNN) based regularization prior. The proposed formulation provides a systematic approach for deriving deep architectures for inverse problems with the arbitrary structure. Since the forward model is explicitly accounted for, a smaller network with fewer parameters is sufficient to capture the image information compared to black-box deep learning approaches, thus reducing the demand for training data and training time. Since we rely on end-to-end training, the CNN weights are customized to the forward model, thus offering improved performance over approaches that rely on pre-trained denoisers. The main difference of the framework from existing end-to-end training strategies is the sharing of the network weights across iterations and channels. Our experiments show that the decoupling of the number of iterations from the network complexity offered by this approach provides benefits including lower demand for training data, reduced risk of overfitting, and implementations with significantly reduced memory footprint. We propose to enforce data-consistency by using numerical optimization blocks such as conjugate gradients algorithm within the network; this approach offers faster convergence per iteration, compared to methods that rely on proximal gradients steps to enforce data consistency. Our experiments show that the faster convergence translates to improved performance, especially when the available GPU memory restricts the number of iterations.


title: Double-negative electromagnetic metamaterials due to chirality

abstract: The aim of this paper is to provide a mathematical theory for understanding the mechanism behind the double-negative refractive index phenomenon in chiral materials. The design of double-negative metamaterials generally requires the use of two different kinds of subwavelength resonators, which may limit the applicability of double-negative metamaterials. Herein, we rely on media that consist of only a single type of dielectric resonant element, and show how the chirality of the background medium induces double-negative refractive index metamaterial, which refracts waves negatively, hence acting as a superlens. Using plasmonic dielectric particles, it is proved that both the effective electric permittivity and the magnetic permeability can be negative near some resonant frequencies. A justification of the approximation of a plasmonic particle in a chiral medium by the sum of a resonant electric dipole and a resonant magnetic dipole, is provided. Moreover, the set of resonant frequencies is characterized. For an appropriate volume fraction of plasmonic particles with certain conditions on their configuration, a double-negative effective medium can be obtained when the frequency is near one of the resonant frequencies.


title: Learned Perceptual Image Enhancement

abstract: Learning a typical image enhancement pipeline involves minimization of a loss function between enhanced and reference images. While L1 and L2 losses are perhaps the most widely used functions for this purpose, they do not necessarily lead to perceptually compelling results. In this paper, we show that adding a learned no-reference image quality metric to the loss can significantly improve enhancement operators. This metric is implemented using a CNN (convolutional neural network) trained on a large-scale dataset labelled with aesthetic preferences of human raters. This loss allows us to conveniently perform back-propagation in our learning framework to simultaneously optimize for similarity to a given ground truth reference and perceptual quality. This perceptual loss is only used to train parameters of image processing operators, and does not impose any extra complexity at inference time. Our experiments demonstrate that this loss can be effective for tuning a variety of operators such as local tone mapping and dehazing.


title: Large isotropic negative thermal expansion above a structural quantum phase transition

abstract: Perovskite structured materials contain myriad tunable ordered phases of electronic and magnetic origin with proven technological importance and strong promise for a variety of energy solutions. An always-contributing influence beneath these cooperative and competing interactions is the lattice, whose physics may be obscured in complex perovskites by the many coupled degrees of freedom which makes these systems interesting. Here we report signatures of an approach to a quantum phase transition very near the ground state of the nonmagnetic, ionic insulating, simple cubic perovskite material ScF3 and show that its physical properties are strongly effected as much as 100 K above the putative transition. Spatial and temporal correlations in the high-symmetry cubic phase determined using energy- and momentum-resolved inelastic X-ray scattering as well as X-ray diffraction reveal that soft mode, central peak and thermal expansion phenomena are all strongly influenced by the transition.


title: Quantum Calculations on the Kv1.2 Channel Voltage Sensing Domain Show H+ Transfer Provides the Gating Current

abstract: Quantum calculations on the voltage sensing domain (VSD) of the Kv1.2 potassium channel (pdb: 3Lut)have been carried out on a 904 atoms subset of the VSD, plus 24 water molecules. Side chains pointing away from the center of the VSD were truncated; S1,S2,S3 end atoms were were fixed (all calculations); S4 end atoms could be fixed or free. Open conformations (membrane potentials >= 0) closely match the known X-ray structure of the open state with salt bridges in the in the VSD not ionized (H+ on the acid) whether S4 end atoms were fixed or free (slightly closer fixed than free).The S4 segment backbone, free or not, moves less than 2.5 A for positive to negative membrane potential switches, not entirely in the expected direction, leaving H+ motion as the principal component of the gating current. Groups of 3 - 5 side chains are important for proton transport, based on the calculations. A proton transfers from tyrosine (Y266), through arginine (R300), to glutamate (E183), accounting for approximately 20 - 25% of the gating charge. Clusters of amino acids that can transfer protons (acids, bases, tyrosine, histidine) are the main paths for proton transport. A group of five amino acids, bounded by the conserved aromatic F233, appears to exchange a proton. Dipole rotations may also contribute. A proton path (calculations still in progress) is proposed for the remainder of the VSD, suggesting a hypothesis for a complete gating mechanism.


title: Alternative Model of Space-Charge Limited Thermionic Current Flow Through a Plasma

abstract: It is widely assumed that thermionic current flow through a plasma is limited by a "space-charge limited" (SCL) cathode sheath that consumes the hot cathode's negative bias and accelerates upstream ions into the cathode. Here, we formulate a fundamentally different current-limited mode. In the "inverse" mode, the potentials of both electrodes are above the plasma potential, so that the plasma ions are confined. The bias is consumed by the anode sheath. There is no potential gradient in the neutral plasma region from resistivity or presheath. The inverse cathode sheath pulls some thermoelectrons back to the cathode, thereby limiting the circuit current. Thermoelectrons entering the zero-field plasma region that undergo collisions may also be sent back to the cathode, further attenuating the circuit current. In planar geometry, the plasma density is shown to vary linearly across the electrode gap. A continuum kinetic planar plasma diode simulation model is set up to compare the properties of current modes with classical, conventional SCL, and inverse cathode sheaths. SCL modes can exist only if charge-exchange collisions are turned off in the potential well of the virtual cathode to prevent ion trapping. With the collisions, the current-limited equilibrium must be inverse. Inverse operating modes should therefore be present or possible in many plasma devices that rely on hot cathodes. Evidence from past experiments is discussed. The inverse mode may offer opportunities to minimize sputtering and power consumption that were not previously explored due to the common assumption of SCL sheaths.


title: Ultracold molecules for measuring the electron's electric dipole moment

abstract: We demonstrate one-dimensional sub-Doppler laser cooling of a beam of YbF molecules to 100 $\mu$K. This is a key step towards a measurement of the electron's electric dipole moment using ultracold molecules. We compare the effectiveness of magnetically-assisted and polarization-gradient sub-Doppler cooling mechanisms. We model the experiment and find good agreement with our data.


title: Perspectival Knowledge in PSOA RuleML: Representation, Model Theory, and Translation

abstract: In Positional-Slotted Object-Applicative (PSOA) RuleML, a predicate application (atom) can have an Object IDentifier (OID) and descriptors that may be positional arguments (tuples) or attribute-value pairs (slots). PSOA RuleML explicitly specifies for each descriptor whether it is to be interpreted under the perspective of the predicate in whose scope it occurs. This predicate-dependency dimension refines the space between oidless, positional atoms (relationships) and oidful, slotted atoms (framepoints): While relationships use only a predicate-scope-sensitive (predicate-dependent) tuple and framepoints use only predicate-scope-insensitive (predicate-independent) slots, PSOA uses a systematics of orthogonal constructs also permitting atoms with (predicate-)independent tuples and atoms with (predicate-)dependent slots. This supports data and knowledge representation where a slot attribute can have different values depending on the predicate. PSOA thus extends object-oriented multi-membership and multiple inheritance. Based on objectification, PSOA laws are given: Besides unscoping and centralization, the semantic restriction and transformation of describution permits rescoping of one atom's independent descriptors to another atom with the same OID but a different predicate. For inheritance, default descriptors are realized by rules. On top of a metamodel and a Grailog visualization, PSOA's atom systematics for facts, queries, and rules is explained. The presentation and (XML-)serialization syntaxes of PSOA RuleML are introduced. Its model-theoretic semantics is formalized by extending the interpretation functions for dependent descriptors. The open-source PSOATransRun system realizes PSOA RuleML by a translator to runtime predicates, including for dependent tuples (prdtupterm) and slots (prdsloterm). Our tests show efficiency advantages of dependent and tupled modeling.


title: The Complexity of Maximum $k$-Order Bounded Component Set Problem

abstract: Given a graph $G=(V, E)$ and a positive integer $k$, in Maximum $k$-Order Bounded Component Set (Max-$k$-OBCS), it is required to find a vertex set $S \subseteq V$ of maximum size such that each component in the induced graph $G[S]$ has at most $k$ vertices. We prove that for constant $k$, Max-$k$-OBCS is hard to approximate within a factor of $n^{1 -\epsilon}$, for any $\epsilon > 0$, unless $\mathsf{P} = \mathsf{NP}$. This is an improvement on the previous lower bound of $\sqrt{n}$ for Max-2-OBCS due to Orlovich et al. We provide lower bounds on the approximability when $k$ is not a constant as well. Max-$k$-OBCS can be seen as a generalization of Maximum Independent Set (Max-IS). We generalize Tur\'an's greedy algorithm for Max-IS and prove that it approximates Max-$k$-OBCS within a factor of $(2k - 1)\overline{d} + k$, where $\overline{d}$ is the average degree of the input graph $G$. This approximation factor is a generalization of Tur\'an's approximation factor for Max-IS.


title: Approximate public-signal correlated equilibria for nonzero-sum differential games

abstract: We construct an approximate public-signal correlated equilibrium for a nonzero-sum differential game in the class of stochastic strategies with memory. The construction is based on a solution of an auxiliary nonzero-sum continuous-time stochastic game. This class of games includes stochastic differential games and continuous-time Markov games. Moreover, we study the limit of approximate equilibrium outcomes in the case when the auxiliary stochastic games tend to the original deterministic one. We show that it lies in the convex hull of the set of equilibrium values provided by deterministic punishment strategies.


title: Dynamic Fault Trees Analysis using an Integration of Theorem Proving and Model Checking

abstract: Dynamic fault trees (DFTs) have emerged as an important tool for capturing the dynamic behavior of system failure. These DFTs are then analyzed qualitatively and quantitatively using stochastic or algebraic methods to judge the failure characteristics of the given system in terms of the failures of its sub-components. Model checking has been recently proposed to conduct the failure analysis of systems using DFTs with the motivation to provide a rigorous failure analysis of safety-critical systems. However, model checking has not been used for the DFT qualitative analysis and the reduction algorithms used in model checking are usually not formally verified. Moreover, the analysis time grows exponentially with the increase of the number of states. These issues limit the usefulness of model checking for analyzing complex systems used in safety-critical domains, where the accuracy and completeness of analysis matters the most. To overcome these limitations, we propose a comprehensive methodology to perform the qualitative and quantitative analysis of DFTs using an integration of theorem proving and model checking based approaches. For this purpose, we formalized all the basic dynamic fault tree gates using higher-order logic based on the algebraic approach and formally verified some of the simplification properties. This formalization allows us to formally verify the equivalence between the original and reduced DFTs using a theorem prover, and conduct the qualitative analysis. We then use model checking to perform the quantitative analysis of the formally verified reduced DFT. We applied our methodology to five benchmarks and the results show that the formally verified reduced DFT was analyzed using model checking with up to six times less states and up to 133000 times faster.


title: Shear-banding and superdiffusivity in entangled polymer solutions

abstract: Using high-resolution confocal rheometry, we study the shear profiles of well-entangled DNA solutions under large amplitude oscillatory shear (LAOS) in a rectilinear planar shear cell. With increasing Weissenberg number (Wi), we observe successive transitions from normal Newtonian linear shear profiles to wall-slip dominant shear profiles and finally to shear-banding profiles at high Wi. To investigate the microscopic origin of the observed shear banding, we study the dynamics of micron-sized tracers embedded in DNA solutions. Surprisingly, tracer particles in the shear frame exhibit transient super-diffusivity and strong dynamic heterogeneity. The probability distribution functions of particle displacements follow a power-law scaling at large displacements, indicating a Levy-walk-type motion, reminiscent of tracer dynamics in entangled wormlike micelle solutions and sheared colloidal glasses. We further characterize the length and time scales associated with the abnormal dynamics of tracer particles. We hypothesize that the unusual particle dynamics arise from localized shear-induced chain disentanglement.


title: Multi-Scale Video Frame-Synthesis Network with Transitive Consistency Loss

abstract: Traditional approaches to interpolate/extrapolate frames in a video sequence require accurate pixel correspondences between images, e.g., using optical flow. Their results stem on the accuracy of optical flow estimation, and could generate heavy artifacts when flow estimation failed. Recently methods using auto-encoder has shown impressive progress, however they are usually trained for specific interpolation/extrapolation settings and lack of flexibility and In order to reduce these limitations, we propose a unified network to parameterize the interest frame position and therefore infer interpolate/extrapolate frames within the same framework. To achieve this, we introduce a transitive consistency loss to better regularize the network. We adopt a multi-scale structure for the network so that the parameters can be shared across multi-layers. Our approach avoids expensive global optimization of optical flow methods, and is efficient and flexible for video interpolation/extrapolation applications. Experimental results have shown that our method performs favorably against state-of-the-art methods.


title: One More Way to Encrypt a Message

abstract: This work describes an example of an application of a novel method for symmetric cryptography. Its purpose is to show how a regular message can be encrypted and then decrypted in an easy, yet secure way. The encrypting method introduced in this work is different from others because it involves decimals as well as integers, encrypting the same initial message differently every time, and inserting misleading digits into every encrypted message, thus making the task of breaking the code even harder. A C++ program was written to support each chapter.


title: Optimal Scale Invariant Wigner Spectrum Estimation of Gaussian Locally Self-Similar Processes Using Hermite Functions

abstract: This paper investigates the mean square error optimal estimation of scale invariant Wigner spectrum for the class of Gaussian locally self-similar processes, by the multitaper method. In this method, the spectrum is estimated as a weighted sum of scale invariant windowed spectrograms. Moreover, it is shown that the optimal multitapers are approximated by the quasi Lamperti transformation of Hermite functions, which is computationally more efficient. Finally, the performance and accuracy of the estimation is studied via simulation.


title: An End to End Deep Neural Network for Iris Segmentation in Unconstraint Scenarios

abstract: With the increasing imaging and processing capabilities of today's mobile devices, user authentication using iris biometrics has become feasible. However, as the acquisition conditions become more unconstrained and as image quality is typically lower than dedicated iris acquisition systems, the accurate segmentation of iris regions is crucial for these devices. In this work, an end to end Fully Convolutional Deep Neural Network (FCDNN) design is proposed to perform the iris segmentation task for lower-quality iris images. The network design process is explained in detail, and the resulting network is trained and tuned using several large public iris datasets. A set of methods to generate and augment suitable lower quality iris images from the high-quality public databases are provided. The network is trained on Near InfraRed (NIR) images initially and later tuned on additional datasets derived from visible images. Comprehensive inter-database comparisons are provided together with results from a selection of experiments detailing the effects of different tunings of the network. Finally, the proposed model is compared with SegNet-basic, and a near-optimal tuning of the network is compared to a selection of other state-of-art iris segmentation algorithms. The results show very promising performance from the optimized Deep Neural Networks design when compared with state-of-art techniques applied to the same lower quality datasets.


title: Median inverse problem and approximating the number of $k$-median inverses of a permutation

abstract: We introduce the "Median Inverse Problem" for metric spaces. In particular, having a permutation $\pi$ in the symmetric group $S_n$ (endowed with the breakpoint distance), we study the set of all $k$-subsets $\{x_1,...,x_k\}\subset S_n$ for which $\pi$ is a breakpoint median. The set of all $k$-tuples $(x_1,...,x_k)$ with this property is called the $k$-median inverse of $\pi$. Finding an upper bound for the cardinality of this set, we provide an asymptotic upper bound for the probability that $\pi$ is a breakpoint median of $k$ permutations $\xi_1^{(n)},...,\xi_k^{(n)}$ chosen uniformly and independently at random from $S_n$.


title: Microscopic optical potentials derived from ab initio translationally invariant nonlocal one-body densities

abstract: We derive a microscopic optical potential for intermediate energies using ab initio translationally invariant nonlocal one-body nuclear densities computed within the no-core shell model (NCSM) approach utilizing two- and three-nucleon chiral interactions as the only input. The optical potential is derived at first-order within the spectator expansion of the non-relativistic multiple scattering theory by adopting the impulse approximation and using the same chiral nucleon-nucleon interaction as that used to compute densities. The ground state local and nonlocal densities of 4,6,8He, 12C, and 16O are calculated and applied to optical potential construction. The differential cross sections and the analyzing powers for the elastic proton scattering off of these nuclei are then calculated for different values of the incident proton energy. The impact of nonlocality and the COM removal is discussed. The use of nonlocal densities has a substantial impact on the differential cross sections and improves agreement with experiment in comparison to results generated with the local densities especially for light nuclei.


title: Universal classes near $\aleph_1$

abstract: Shelah has provided sufficient conditions for an $L_{\omega_1, \omega}$-sentence $\psi$ to have arbitrarily large models and for a Morley-like theorem to hold of $\psi$. These conditions involve structural and set-theoretic assumptions on all the $\aleph_n$'s. Using tools of Boney, Shelah, and the second author, we give assumptions on $\aleph_0$ and $\aleph_1$ which suffice when $\psi$ is restricted to be universal: $\mathbf{Theorem}$ Assume $2^{\aleph_{0}} < 2 ^{\aleph_{1}}$. Let $\psi$ be a universal $L_{\omega_{1}, \omega}$-sentence. - If $\psi$ is categorical in $\aleph_{0}$ and $1 \leq I(\psi, \aleph_{1}) < 2 ^{\aleph_{1}}$, then $\psi$ has arbitrarily large models and categoricity of $\psi$ in some uncountable cardinal implies categoricity of $\psi$ in all uncountable cardinals. - If $\psi$ is categorical in $\aleph_1$, then $\psi$ is categorical in all uncountable cardinals. The theorem generalizes to the framework of $L_{\omega_1, \omega}$-definable tame abstract elementary classes with primes.


title: A Pilot Study on Using an Intelligent Life-like Robot as a Companion for Elderly Individuals with Dementia and Depression

abstract: This paper presents the design, development, methodology, and the results of a pilot study on using an intelligent, emotive and perceptive social robot (aka Companionbot) for improving the quality of life of elderly people with dementia and/or depression. Ryan Companionbot prototyped in this project, is a rear-projected life-like conversational robot. Ryan is equipped with features that can (1) interpret and respond to users' emotions through facial expressions and spoken language, (2) proactively engage in conversations with users, and (3) remind them about their daily life schedules (e.g. taking their medicine on time). Ryan engages users in cognitive games and reminiscence activities. We conducted a pilot study with six elderly individuals with moderate dementia and/or depression living in a senior living facility in Denver. Each individual had 24/7 access to a Ryan in his/her room for a period of 4-6 weeks. Our observations of these individuals, interviews with them and their caregivers, and analyses of their interactions during this period revealed that they established rapport with the robot and greatly valued and enjoyed having a Companionbot in their room.


title: Columnar Database Techniques for Creating AI Features

abstract: Recent advances with in-memory columnar database techniques have increased the performance of analytical queries on very large databases and data warehouses. At the same time, advances in artificial intelligence (AI) algorithms have increased the ability to analyze data. We use the term AI to encompass both Deep Learning (DL or neural network) and Machine Learning (ML aka Big Data analytics). Our exploration of the AI full stack has led us to a cross-stack columnar database innovation that efficiently creates features for AI analytics. The innovation is to create Augmented Dictionary Values (ADVs) to add to existing columnar database dictionaries in order to increase the efficiency of featurization by minimizing data movement and data duplication. We show how various forms of featurization (feature selection, feature extraction, and feature creation) can be efficiently calculated in a columnar database. The full stack AI investigation has also led us to propose an integrated columnar database and AI architecture. This architecture has information flows and feedback loops to improve the whole analytics cycle during multiple iterations of extracting data from the data sources, featurization, and analysis.


title: Perpendicular relativistic shocks in magnetized pair plasma

abstract: Perpendicular relativistic ($\gamma_0=10$) shocks in magnetized pair plasmas are investigated using two dimensional Particle-in-Cell simulations. A systematic survey, from unmagnetized to strongly magnetized shocks, is presented accurately capturing the transition from Weibel-mediated to magnetic-reflection-shaped shocks. This transition is found to occur for upstream flow magnetizations $10^{-3}<\sigma<10^{-2}$ at which a strong perpendicular net current is observed in the precursor, driving the so-called current-filamentation instability. The global structure of the shock and shock formation time are discussed. The MHD shock jump conditions are found in good agreement with the numerical results, except for $10^{-4} < \sigma < 10^{-2}$ where a deviation up to 10\% is observed. The particle precursor length converges toward the Larmor radius of particles injected in the upstream magnetic field at intermediate magnetizations. For $\sigma>10^{-2}$, it leaves place to a purely electromagnetic precursor following from the strong emission of electromagnetic waves at the shock front. Particle acceleration is found to be efficient in weakly magnetized perpendicular shocks in agreement with previous works, and is fully suppressed for $\sigma > 10^{-2}$. Diffusive Shock Acceleration is observed only in weakly magnetized shocks, while a dominant contribution of Shock Drift Acceleration is evidenced at intermediate magnetizations. The spatial diffusion coefficients are extracted from the simulations allowing for a deeper insight into the self-consistent particle kinematics and scale with the square of the particle energy in weakly magnetized shocks. These results have implications for particle acceleration in the internal shocks of AGN jets and in the termination shocks of Pulsar Wind Nebulae.


title: Parton correlations in double parton scattering

abstract: Double parton scattering events are directly sensitive to the correlations between two partons inside a proton and can answer fundamental questions on the connections between the proton constituents. In this chapter, the different types of possible correlations, our present knowledge of them, and the processes where they are likely to be important, are introduced and explained. The increasing integrated luminosity at the LHC and the refinements of the theory of double parton scattering, lead to interesting prospects for measuring, or severely constraining, two-parton correlations in the near future.


title: The Large Synoptic Survey Telescope and Milky Way Science

abstract: The Large Synoptic Survey Telescope (LSST) surveys have initially been optimized to omit the inner part of the Milky Way disk/bar from deep and cadence observations. However it is now clear that the LSST will be powerful for Galactic astronomy and may play a crucial role in continuing to extend the Gaia astrometric catalog until a future satellite, either optical or IR, carries on. LSST will provide metallicities and kinematics for the bulge, and will map halo structures to as distant as 450 kpc, nearly half the distance to the Andromeda galaxy. Thanks to the unprecedented calibration effort for its photometric system, and surprisingly good astrometry (transverse velocity measurements of 0.2 mas/yr at r=21; 1 mas/yr at r=24). LSST will provide photometric abundances and distance constraints for a billion or more Milky Way stars to distances of 450 kpc, and kinematics from proper motions to $\sim100$ kpc. Single observation depths reach $\sim 24$ in the $ugrizy$ bands, while depths at end of mission reach $\sim 27$. Although halo structures such as streams and dwarf galaxies are initially identified by the RR Lyrae and giants, their structure will be fleshed out by the 100$\times$ more abundant dwarfs that will be detected to 100 kpc (single observation) and $\sim 300$ kpc by end of mission. More complete mapping of stream structures may constrain the mass distribution of dark matter and perhaps confirm the interaction of dark matter halos and streams. I also describe the Blanco DECam Bulge Survey, a 200 deg$^2$ LSST pathfinder survey of the bulge in $ugrizy$ using the Dark Energy Camera on the Blanco 4m telescope. The purpose of this article is to encourage active workers on the Milky Way and Local Volume to participate in the LSST project, in particular to urge that the Galactic Plane receive the same cadence and depth coverage as the rest of the extragalactic sky.


title: Cosmological constraints from galaxy clustering in the presence of massive neutrinos

abstract: The clustering ratio is defined as the ratio between the correlation function and the variance of the smoothed overdensity field. In LCDM cosmologies not accounting for massive neutrinos, it has already been proved to be independent from bias and redshift space distortions on a range of linear scales. It therefore allows for a direct comparison of measurements (from galaxies in redshift space) to predictions (for matter in real space). In this paper we first extend the applicability of such properties of the clustering ratio to cosmologies that include massive neutrinos, by performing tests against simulated data. We then investigate the constraining power of the clustering ratio when cosmological parameters such as the total neutrino mass and the equation of state of dark energy are left free. We analyse the joint posterior distribution of the parameters that must satisfy, at the same time, the measurements of the galaxy clustering ratio in the SDSS DR12, and the angular power spectrum of temperature and polarization anisotropies of the CMB measured by the Planck satellite. We find the clustering ratio to be very sensitive to the CDM density parameter, but not very much so to the total neutrino mass. Lastly, we forecast the constraining power the clustering ratio will achieve with forthcoming surveys, predicting the amplitude of its errors in a Euclid-like galaxy survey. In this case, we find it is expected to improve the constraint at 95% level on the CDM density by 40% and on the total neutrino mass by 14%.


title: Enhanced signal-to-noise ratio in Hanbury Brown Twiss interferometry by parametric amplification

abstract: The Hanbury Brown Twiss (HBT) interferometer was proposed to observe intensity correlations of starlight to measure a star's angular diameter. As the intensity of light that reaches the detector from a star is very weak, one cannot usually get a workable signal-to-noise ratio. We propose an improved HBT interferometric scheme introducing optical parametric amplifiers into the system, to amplify the correlation signal, which is used to calculate the angular diameter. With the use of optical parametric amplifiers, the signal-to-noise ratio can be increased up to 400 percent.


title: Caustic Formation upon Shift Symmetry Breaking

abstract: We examine how the breaking of shift symmetry affects the formation of caustics for the standard canonical kinetic theory as well as for the DBI theory. We show in this case, that the standard canonical kinetic theory is caustic free but the same does not always apply for the DBI model. We make similar arguments about the conformal Galileons. Finally, it is shown that the simple wave condition is not invariant under field redefinitions and the meaning of simple waves after symmetry breaking is discussed.


title: Whole-Body Nonlinear Model Predictive Control Through Contacts for Quadrupeds

abstract: In this work we present a whole-body Nonlinear Model Predictive Control approach for Rigid Body Systems subject to contacts. We use a full dynamic system model which also includes explicit contact dynamics. Therefore, contact locations, sequences and timings are not prespecified but optimized by the solver. Yet, thorough numerical and software engineering allows for running the nonlinear Optimal Control solver at rates up to 190 Hz on a quadruped for a time horizon of half a second. This outperforms the state of the art by at least one order of magnitude. Hardware experiments in form of periodic and non-periodic tasks are applied to two quadrupeds with different actuation systems. The obtained results underline the performance, transferability and robustness of the approach.


title: Network Analysis for Explanation

abstract: Safety critical systems strongly require the quality aspects of artificial intelligence including explainability. In this paper, we analyzed a trained network to extract features which mainly contribute the inference. Based on the analysis, we developed a simple solution to generate explanations of the inference processes.


title: Circuit Preserving Edge Maps II

abstract: The results obtained in this paper grew from an attempt to generalize the main theorem of [1]. There it was shown that any circuit injection (a 1-1 onto edge map f such that if C is a circuit then f(C) is a circuit) from a 3-connected, not necessarily finite graph G onto a graph H is induced by a vertex isomorphism, where H is assumed to not have any isolated vertices. In the present article we examine the situation when the 1-1 condition is dropped (Chapter 1). An interesting result then is that the theorem remains true for finite (3-connected) graphs G but not for infinite G. In Chapter 2 we retain the 1-1 condition but allow the image of f to be first an arbitrary matroid and second a binary matroid. An interesting result then is the following. Let G be a graph of even order. Then the statement "no nontrivial map f:=>M exists, where M is a binary matroid" is equivalent to "G is Hamiltonian". If G is a graph of odd order, then the statement "no nontrivial map f:G=>M exists, where M is a binary matroid" is equivalent to "G is almost Hamiltonian", where we define a graph G of order n to be almost Hamiltonian if every subset of vertices of order n-1 is contained in some circuit of G. [1] J.H. Sanders and D. Sanders, Circuit preserving edge maps, J. Combin. Theory Ser. B 22 (1977),91-96.


title: Gouy Phase Radial Mode Sorter for Light: Concepts and Experiments

abstract: We present an in principle lossless sorter for radial modes of light, using accumulated Gouy phases. The experimental setups have been found by a computer algorithm, and can be intuitively understood in a geometric way. Together with the ability to sort angular-momentum modes, we now have access to the complete 2-dimensional transverse plane of light. The device can readily be used in multiplexing classical information. On a quantum level, it is an analog of the Stern-Gerlach experiment -- significant for the discussion of fundamental concepts in quantum physics. As such, it can be applied in high-dimensional and multi-photonic quantum experiments.


title: Deep Texture and Structure Aware Filtering Network for Image Smoothing

abstract: Image smoothing is a fundamental task in computer vision, that aims to retain salient structures and remove insignificant textures. In this paper, we aim to address the fundamental shortcomings of existing image smoothing methods, which cannot properly distinguish textures and structures with similar low-level appearance. While deep learning approaches have started to explore the preservation of structure through image smoothing, existing work does not yet properly address textures. To this end, we generate a large dataset by blending natural textures with clean structure-only images, and then build a texture prediction network (TPN) that predicts the location and magnitude of textures. We then combine the TPN with a semantic structure prediction network (SPN) so that the final texture and structure aware filtering network (TSAFN) is able to identify the textures to remove ("texture-awareness") and the structures to preserve ("structure-awareness"). The proposed model is easy to understand and implement, and shows excellent performance on real images in the wild as well as our generated dataset.


title: "Space is blue and birds fly through it"

abstract: Quantum mechanics is not about 'quantum states': it is about values of physical variables. I give a short fresh presentation and update on the $relational$ perspective on the theory, and a comment on its philosophical implications.


title: Families of Quantum Fingerprinting Protocols

abstract: We introduce several families of quantum fingerprinting protocols to evaluate the equality function on two $n$-bit strings in the simultaneous message passing model. The original quantum fingerprinting protocol uses a tensor product of a small number of $\mathcal{O}(\log n)$-qubit high dimensional signals [Buhrman et al. 2001], whereas a recently-proposed optical protocol uses a tensor product of $\mathcal{O}(n)$ single-qubit signals, while maintaining the $\mathcal{O}(\log n)$ information leakage of the original protocol [Arrazola and L\"utkenhaus 2014]. We find a family of protocols which interpolate between the original and optical protocols while maintaining the $\mathcal{O}(\log n)$ information leakage, thus demonstrating a trade-off between the number of signals sent and the dimension of each signal. There has been interest in experimental realization of the recently-proposed optical protocol using coherent states [Xu et al. 2015, Guan et al. 2016], but as the required number of laser pulses grows linearly with the input size $n$, eventual challenges for the long-time stability of experimental set-ups arise. We find a coherent state protocol which reduces the number of signals by a factor $1/2$ while also reducing the information leakage. Our reduction makes use of a simple modulation scheme in optical phase space, and we find that more complex modulation schemes are not advantageous. Using a similar technique, we improve a recently-proposed coherent state protocol for evaluating the Euclidean distance between two real unit vectors [Kumar et al. 2017] by reducing the number of signals by a factor $1/2$ and also reducing the information leakage.


title: Audio-Visual Sentiment Analysis for Learning Emotional Arcs in Movies

abstract: Stories can have tremendous power -- not only useful for entertainment, they can activate our interests and mobilize our actions. The degree to which a story resonates with its audience may be in part reflected in the emotional journey it takes the audience upon. In this paper, we use machine learning methods to construct emotional arcs in movies, calculate families of arcs, and demonstrate the ability for certain arcs to predict audience engagement. The system is applied to Hollywood films and high quality shorts found on the web. We begin by using deep convolutional neural networks for audio and visual sentiment analysis. These models are trained on both new and existing large-scale datasets, after which they can be used to compute separate audio and visual emotional arcs. We then crowdsource annotations for 30-second video clips extracted from highs and lows in the arcs in order to assess the micro-level precision of the system, with precision measured in terms of agreement in polarity between the system's predictions and annotators' ratings. These annotations are also used to combine the audio and visual predictions. Next, we look at macro-level characterizations of movies by investigating whether there exist `universal shapes' of emotional arcs. In particular, we develop a clustering approach to discover distinct classes of emotional arcs. Finally, we show on a sample corpus of short web videos that certain emotional arcs are statistically significant predictors of the number of comments a video receives. These results suggest that the emotional arcs learned by our approach successfully represent macroscopic aspects of a video story that drive audience engagement. Such machine understanding could be used to predict audience reactions to video stories, ultimately improving our ability as storytellers to communicate with each other.


title: Bounding periods of subvarieties of (P^1)^n

abstract: Using methods of p-adic analysis, along with the powerful result of Medvedev-Scanlon (Annals of Mathematics, 2014) for the classification of periodic subvarieties of (P^1)^n, we bound the length of the orbit of a periodic subvariety Y of (P^1)^n under the action of a dominant endomorphism.


title: Representations of Sound in Deep Learning of Audio Features from Music

abstract: The work of a single musician, group or composer can vary widely in terms of musical style. Indeed, different stylistic elements, from performance medium and rhythm to harmony and texture, are typically exploited and developed across an artist's lifetime. Yet, there is often a discernable character to the work of, for instance, individual composers at the perceptual level - an experienced listener can often pick up on subtle clues in the music to identify the composer or performer. Here we suggest that a convolutional network may learn these subtle clues or features given an appropriate representation of the music. In this paper, we apply a deep convolutional neural network to a large audio dataset and empirically evaluate its performance on audio classification tasks. Our trained network demonstrates accurate performance on such classification tasks when presented with 5 s examples of music obtained by simple transformations of the raw audio waveform. A particularly interesting example is the spectral representation of music obtained by application of a logarithmically spaced filter bank, mirroring the early stages of auditory signal transduction in mammals. The most successful representation of music to facilitate discrimination was obtained via a random matrix transform (RMT). Networks based on logarithmic filter banks and RMT were able to correctly guess the one composer out of 31 possibilities in 68 and 84 percent of cases respectively.


title: A Taxonomy and Future Directions for Sustainable Cloud Computing: 360 Degree View

abstract: The cloud computing paradigm offers on-demand services over the Internet and supports a wide variety of applications. With the recent growth of Internet of Things (IoT) based applications the usage of cloud services is increasing exponentially. The next generation of cloud computing must be energy-efficient and sustainable to fulfil the end-user requirements which are changing dynamically. Presently, cloud providers are facing challenges to ensure the energy efficiency and sustainability of their services. The usage of large number of cloud datacenters increases cost as well as carbon footprints, which further effects the sustainability of cloud services. In this paper, we propose a comprehensive taxonomy of sustainable cloud computing. The taxonomy is used to investigate the existing techniques for sustainability that need careful attention and investigation as proposed by several academic and industry groups. Further, the current research on sustainable cloud computing is organized into several categories: application design, sustainability metrics, capacity planning, energy management, virtualization, thermal-aware scheduling, cooling management, renewable energy and waste heat utilization. The existing techniques have been compared and categorized based on the common characteristics and properties. A conceptual model for sustainable cloud computing has been proposed along with discussion on future research directions.


title: Addressing the gas kinetics Boltzmann equation with branching-path statistics

abstract: This article proposes a new statistical numerical method to address gas kinetics problems obeying the Boltzmann equation. This method is inspired from some Monte-Carlo algorithms used in linear transport physics, where virtual particles are followed backwards in time along their paths. The non-linear character of gas kinetics translates, in the numerical simulations presented here, in branchings of the virtual particle paths. The obtained algorithms have displayed in the few tests presented here two noticeable qualities: (1) They involve no mesh. (2) They allow to easily compute the gas density at rarefied places of the phase space, for example at high kinetic energy.


title: Disintegration of an Eruptive Filament via Interactions with Quasi-Separatrix Layers

abstract: The disintegration of solar filaments via mass drainage is a frequently observed phenomenon during a variety of filament activities. It is generally considered that the draining of dense filament material is directed by both gravity and magnetic field, yet the detailed process remains elusive. Here we report on a partial filament eruption during which filament material drains downward to the surface not only along the filament's legs, but to a remote flare ribbon through a fan-out curtain-like structure. It is found that the magnetic configuration is characterized by two conjoining dome-like quasi-sepratrix layers (QSLs). The filament is located underneath one QSL dome, whose footprint apparently bounds the major flare ribbons resulting from the filament eruption, whereas the remote flare ribbon matches well with the other QSL dome's far-side footprint. We suggest that the interaction of the filament with the overlying QSLs results in the splitting and disintegration of the filament.


title: Multiple Adaptive Bayesian Linear Regression for Scalable Bayesian Optimization with Warm Start

abstract: Bayesian optimization (BO) is a model-based approach for gradient-free black-box function optimization. Typically, BO is powered by a Gaussian process (GP), whose algorithmic complexity is cubic in the number of evaluations. Hence, GP-based BO cannot leverage large amounts of past or related function evaluations, for example, to warm start the BO procedure. We develop a multiple adaptive Bayesian linear regression model as a scalable alternative whose complexity is linear in the number of observations. The multiple Bayesian linear regression models are coupled through a shared feedforward neural network, which learns a joint representation and transfers knowledge across machine learning problems.


title: Blind Multiclass Ensemble Classification

abstract: The rising interest in pattern recognition and data analytics has spurred the development of innovative machine learning algorithms and tools. However, as each algorithm has its strengths and limitations, one is motivated to judiciously fuse multiple algorithms in order to find the "best" performing one, for a given dataset. Ensemble learning aims at such high-performance meta-algorithm, by combining the outputs from multiple algorithms. The present work introduces a blind scheme for learning from ensembles of classifiers, using a moment matching method that leverages joint tensor and matrix factorization. Blind refers to the combiner who has no knowledge of the ground-truth labels that each classifier has been trained on. A rigorous performance analysis is derived and the proposed scheme is evaluated on synthetic and real datasets.


title: Experimental demonstration of an isotope-sensitive warhead verification technique using nuclear resonance fluorescence

abstract: Future nuclear arms reduction efforts will require technologies to verify that warheads slated for dismantlement are authentic without revealing any sensitive weapons design information to international inspectors. Despite several decades of research, no technology has met these requirements simultaneously. Recent work by Kemp et al. [Kemp RS, Danagoulian A, Macdonald RR, Vavrek JR (2016) Proc Natl Acad Sci USA 113:8618--8623] has produced a novel physical cryptographic verification protocol that approaches this treaty verification problem by exploiting the isotope-specific nature of nuclear resonance fluorescence (NRF) measurements to verify the authenticity of a warhead. To protect sensitive information, the NRF signal from the warhead is convolved with that of an encryption foil that contains key warhead isotopes in amounts unknown to the inspector. The convolved spectrum from a candidate warhead is statistically compared against that from an authenticated template warhead to determine whether the candidate itself is authentic. Here we report on recent proof-of-concept warhead verification experiments conducted at the Massachusetts Institute of Technology. Using high-purity germanium (HPGe) detectors, we measured NRF spectra from the interrogation of proxy 'genuine' and 'hoax' objects by a 2.52 MeV endpoint bremsstrahlung beam. The observed differences in NRF intensities near 2.2 MeV indicate that the physical cryptographic protocol can distinguish between proxy genuine and hoax objects with high confidence in realistic measurement times.


title: Injection locking of optomechanical oscillators via acoustic waves

abstract: Injection locking is a powerful technique for synchronization of oscillator networks and controlling the phase and frequency of individual oscillators using similar or other types of oscillators. Here, we present the first demonstration of injection locking of a radiation-pressure driven optomechanical oscillator (OMO) via acoustic waves. As opposed to previously reported techniques (based on pump modulation or direct application of a modulated electrostatic force), injection locking of OMO via acoustic waves does not require optical power modulation or physical contact with the OMO and it can easily be implemented on various platforms. Using this approach we have locked the phase and frequency of two distinct modes of a microtoroidal silica OMO to a piezoelectric transducer (PZT). We have characterized the behavior of the injection locked OMO with three acoustic excitation configurations and showed that even without proper acoustic impedance matching the OMO can be locked to the PZT and tuned over 17 kHz with only -30 dBm of RF power fed to the PZT. The high efficiency, simplicity and scalability of the proposed approach paves the road toward a new class of photonic systems that rely on synchronization of several OMOs to a single or multiple RF oscillators with applications in optical communication, metrology and sensing. Beyond its practical applications, injection locking via acoustic waves can be used in fundamental studies in quantum optomechanics where thermal and optical isolation of the OMO are critical.


title: Class numbers and $p$-ranks in ${\mathbb Z}_p^d$-towers

abstract: To extend Iwasawa's classical theorem from ${\mathbb Z}_p$-towers to ${\mathbb Z}_p^d$-towers, Greenberg conjectured that the exponent of $p$ in the $n$-th class number in a ${\mathbb Z}_p^d$-tower of a global field $K$ ramified at finitely many primes is given by a polynomial in $p^n$ and $n$ of total degree at most $d$ for sufficiently large $n$. This conjecture remains open for $d\geq 2$. In this paper, we prove that this conjecture is true in the function field case. Further, we propose a series of general conjectures on $p$-adic stability of zeta functions in a $p$-adic Lie tower of function fields.


title: Equidistribution of dilated curves on nilmanifolds

abstract: Generalizing classic results for a family of measures in the torus, for a family $(\mu_t)_{t\geq 0}$ of measures defined on a nilmanifold $X$, we study conditions under which the family equidistributes, meaning conditions under which the measures $\mu_t$ converge as $t\to\infty$ in the weak$^\ast$ topology to the Haar measure on $X$. We give general conditions on a family of measures defined by a dilation process, showing necessary and sufficient conditions for equidistribution as the family dilates, along with conditions such that this holds for all dilates outside some set of density zero. Furthermore, we show that these two types of equidistribution are different.


title: Two-dimensional modeling of density and thermal structure of dense circumstellar outflowing disks

abstract: Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. Aims. We study the hydrodynamic and thermal structure of optically thick, dense parts of outflowing circumstellar disks that may be formed around various types of critically rotating massive stars, for example, Be stars, B[e] supergiant (sgB[e]) stars or Pop III stars. Methods. We specify the optical depth of the disk along the line-of-sight from stellar poles. Within the optically thick dense region we calculate the vertical disk thermal structure using the diffusion approximation while for the optically thin outer layers we assume a local thermodynamic equilibrium with the impinging stellar irradiation. We use two of our own types of hydrodynamic codes: two-dimensional operator-split numerical code and unsplit code based on the Roe's method. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk. In the models of dense viscous disks the viscosity increases the central temperature up to several tens of thousands of Kelvins. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen.


title: Sharing Storage in a Smart Grid: A Coalitional Game Approach

abstract: Sharing economy is a transformative socio-economic phenomenon built around the idea of sharing underused resources and services, e.g. transportation and housing, thereby reducing costs and extracting value. Anticipating continued reduction in the cost of electricity storage, we look into the potential opportunity in electrical power system where consumers share storage with each other. We consider two different scenarios. In the first scenario, consumers are assumed to already have individual storage devices and they explore cooperation to minimize the realized electricity consumption cost. In the second scenario, a group of consumers is interested to invest in joint storage capacity and operate it cooperatively. The resulting system problems are modeled using cooperative game theory. In both cases, the cooperative games are shown to have non-empty cores and we develop efficient cost allocations in the core with analytical expressions. Thus, sharing of storage in cooperative manner is shown to be very effective for the electric power system.


title: Physarum-inspired Network Optimization: A Review

abstract: The popular Physarum-inspired Algorithms (PAs) have the potential to solve challenging network optimization problems. However, the existing researches on PAs are still immature and far from being fully recognized. A major reason is that these researches have not been well organized so far. In this paper, we aim to address this issue. First, we introduce Physarum and its intelligence from the biological perspective. Then, we summarize and group four types of Physarum-inspired networking models. After that, we analyze the network optimization problems and applications that have been challenged by PAs based on these models. Ultimately, we discuss the existing researches on PAs and identify two fundamental questions: 1) What are the characteristics of Physarum networks? 2) Why can Physarum solve some network optimization problems? Answering these two questions is essential to the future development of Physarum-inspired network optimization.


title: Linked systems of symmetric designs

abstract: A linked system of symmetric designs (LSSD) is a $w$-partite graph ($w\geq 2$) where the incidence between any two parts corresponds to a symmetric design and the designs arising from three parts are related. The original construction for LSSDs by Goethals used Kerdock sets, in which $v$ is a power of two. Some four decades later, new examples were given by Davis et.\ al.\ and Jedwab et.\ al.\ using difference sets, again with $v$ a power of two. In this paper we develop a connection between LSSDs and "linked simplices", full-dimensional regular simplices with two possible inner products between vertices of distinct simplices. We then use this geometric connection to construct sets of equiangular lines and to find an equivalence between regular unbiased Hadamard matrices and certain LSSDs with Menon parameters. We then construct examples of non-trivial LSSDs in which $w$ can be made arbitarily large for fixed even part of $v$. Finally we survey the known infinite families of symmetric designs and show, using basic number theoretic conditions, that $w=2$ in most cases.


title: Exploiting Modern Hardware for High-Dimensional Nearest Neighbor Search

abstract: Many multimedia information retrieval or machine learning problems require efficient high-dimensional nearest neighbor search techniques. For instance, multimedia objects (images, music or videos) can be represented by high-dimensional feature vectors. Finding two similar multimedia objects then comes down to finding two objects that have similar feature vectors. In the current context of mass use of social networks, large scale multimedia databases or large scale machine learning applications are more and more common, calling for efficient nearest neighbor search approaches. This thesis builds on product quantization, an efficient nearest neighbor search technique that compresses high-dimensional vectors into short codes. This makes it possible to store very large databases entirely in RAM, enabling low response times. We propose several contributions that exploit the capabilities of modern CPUs, especially SIMD and the cache hierarchy, to further decrease response times offered by product quantization.


title: Decomposing arrangements of hyperplanes: VC-dimension, combinatorial dimension, and point location

abstract: $\renewcommand{\Re}{\mathbb{R}}$ We re-examine parameters for the two main space decomposition techniques---bottom-vertex triangulation, and vertical decomposition, including their explicit dependence on the dimension $d$, and discover several unexpected phenomena, which show that, in both techniques, there are large gaps between the VC-dimension (and primal shatter dimension), and the combinatorial dimension. For vertical decomposition, the combinatorial dimension is only $2d$, the primal shatter dimension is at most $d(d+1)$, and the VC-dimension is at least $1 + d(d+1)/2$ and at most $O(d^3)$. For bottom-vertex triangulation, both the primal shatter dimension and the combinatorial dimension are $\Theta(d^2)$, but there seems to be a significant gap between them, as the combinatorial dimension is $\frac12d(d+3)$, whereas the primal shatter dimension is at most $d(d+1)$, and the VC-dimension is between $d(d+1)$ and $5d^2 \log{d}$ (for $d\ge 9$). Our main application is to point location in an arrangement of $n$ hyperplanes is $\Re^d$, in which we show that the query cost in Meiser's algorithm can be improved if one uses vertical decomposition instead of bottom-vertex triangulation, at the cost of some increase in the preprocessing cost and storage. The best query time that we can obtain is $O(d^3\log n)$, instead of $O(d^4\log d\log n)$ in Meiser's algorithm. For these bounds to hold, the preprocessing and storage are rather large (super-exponential in $d$). We discuss the tradeoff between query cost and storage (in both approaches, the one using bottom-vertex trinagulation and the one using vertical decomposition).


title: Evidence for nodal superconductivity in a layered compound Ta$_4$Pd$_3$Te$_{16}$

abstract: We report an investigation of the London penetration depth $\Delta\lambda(T)$ on single crystals of the layered superconductor Ta$_4$Pd$_3$Te$_{16}$, where the crystal structure has quasi-one-dimensional characteristics. A linear temperature dependence of $\Delta\lambda(T)$ is observed for $T\ll T_c$, in contrast to the exponential decay of fully gapped superconductors. This indicates the existence of line nodes in the superconducting energy gap. A detailed analysis shows that the normalized superfluid density $\rho_s(T)$, which is converted from $\Delta\lambda(T)$, can be well described by a multigap scenario, with nodes in one of the superconducting gaps, providing clear evidence for nodal superconductivity in Ta$_4$Pd$_3$Te$_{16}$.


title: A Quasi-one-dimensional Analytic Model of Rotating Detonation Combustors

abstract: A quasi-one-dimensional analytic model is proposed for the internal fluid of rotating detonation combustors (RDCs). This model uses the shock-tube model that constrains the flow to have only a longitudinal component, while admitting the propagation of the detonation wave in the azimuthal direction. The proposed model is able to compute the thruster performance and two-dimensional distributions of gas properties. The calculation process of the model excludes iterative calculation or space discretization. The case calculations of the hydrogen-air RDC and the ethylene-oxygen RDC are conducted, and the results calculated by the analytic model are compared with those simulated by computational fluid dynamics (CFD). Good agreement has been observed between the results obtained with the proposed model and CFD, in terms of both of the qualitative and quantitative comparisons. The proposed model is simple and fast, and also maintains the fundamental characteristics of RDCs.


title: Enhanced cooperativity for quantum-nondemolition-measurement--induced spin squeezing of atoms coupled to a nanophotonic waveguide

abstract: We study the enhancement of cooperativity in the atom-light interface near a nanophotonic waveguide for application to quantum nondemolition (QND) measurement of atomic spins. Here the cooperativity per atom is determined by the ratio between the measurement strength and the decoherence rate. Counterintuitively, we find that by placing the atoms at an azimuthal position where the guided probe mode has the lowest intensity, we increase the cooperativity. This arises because the QND measurement strength depends on the interference between the probe and scattered light guided into an orthogonal polarization mode, while the decoherence rate depends on the local intensity of the probe. Thus, by proper choice of geometry, the ratio of good to bad scattering can be strongly enhanced for highly anisotropic modes. We apply this to study spin squeezing resulting from QND measurement of spin projection noise via the Faraday effect in two nanophotonic geometries, a cylindrical nanofiber and a square waveguide. We find that, with about 2500 atoms and using realistic experimental parameters, $ \sim 6.3 $ and $ \sim 13 $ dB of squeezing can be achieved on the nanofiber and square waveguide, respectively.


title: Using Intermittent Synchronization to Compensate for Rhythmic Body Motion During Autonomous Surgical Cutting and Debridement

abstract: Anatomical structures are rarely static during a surgical procedure due to breathing, heartbeats, and peristaltic movements. Inspired by observing an expert surgeon, we propose an intermittent synchronization with the extrema of the rhythmic motion (i.e., the lowest velocity windows). We performed 2 experiments: (1) pattern cutting, and (2) debridement. In (1), we found that the intermittent synchronization approach, while 1.8x slower than tracking motion, was significantly more robust to noise and control latency, and it reduced the max cutting error by 2.6x In (2), a baseline approach with no synchronization achieves 62% success rate for each removal, while intermittent synchronization achieves 80%.


title: On-chip self-sensing enhancing actuation for MEMS/NEMS with large efficiency

abstract: The implementation of on-chip MEMS/NEMS transducers for arbitrary resonators is difficult due to a number of difficulties such as material choice, large dissipation, restriction in high frequency, low sensitivity, poor reliability, and poor integrability. We show a universal on-chip transduction scheme, which can be adapted to any MEMS/NEMS resonator. We achieve all electrical, on-chip MEMS/NEMS for any resonator.


title: On Topology Optimization and Canonical Duality Method

abstract: Topology optimization for general materials is correctly formulated as a bi-level knapsack problem, which is considered to be NP-hard in global optimization and computer science. By using canonical duality theory (CDT) developed by the author, the linear knapsack problem can be solved analytically to obtain global optimal solution at each design iteration. Both uniqueness, existence, and NP-hardness are discussed. The novel CDT method for general topology optimization is refined and tested by both 2-D and 3-D benchmark problems. Numerical results show that without using filter and any other artificial technique, the CDT method can produce exactly 0-1 optimal density distribution with almost no checkerboard pattern. Its performance and novelty are compared with the popular SIMP and BESO approaches. Additionally, some mathematical and conceptual mistakes in literature are explicitly pointed out. A brief review on the canonical duality theory for solving a unified problem in multi-scale nonconvex/discrete systems is given in Appendix.


title: A Further Search for Galactic Stars with Double Radio Lobes

abstract: Over a thousand stars in our Galaxy have been detected as radio emitters, but no normal stars are known to possess radio-emitting lobes similar to radio galaxies. Several recent attempts by us and other authors to find such objects remained inconclusive. Here we present a further search for double-lobed radio stars in two large samples of spectroscopic stars: over 20,000 white dwarves from the Sloan Digital Sky Survey (SDSS) DR12, and 2.5 million stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). These were cross-matched with sources from the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey at 1.4 GHz to look for source pairs straddling the stars with moderate symmetry about the stars. We found only four promising candidates for double-lobed radio stars, confirming they must be extremely rare. By comparison with SDSS, we inferred that about 16 per cent of LAMOST spectra may have erroneous classifications. We also rediscovered the giant radio galaxy J0927+3510 and propose a different, more distant host, suggesting a much larger radio size of 2.7 Mpc.


title: Stability of fractional-order nonlinear systems by Lyapunov direct method

abstract: In this paper, by using a characterization of functions having fractional derivative, we propose a rigorous fractional Lyapunov function candidate method to analyze stability of fractional-order nonlinear systems. First, we prove an inequality concerning the fractional derivatives of convex Lyapunov functions without the assumption on the existence of derivative of pseudo-states. Second, we establish fractional Lyapunov functions to fractional-order systems without the assumption on the global existence of solutions. Our theorems fill the gaps and strengthen results in some existing papers.


title: Real Zero Polynomials and A. Horn's Problem

abstract: A. Horn's problem concerns find two self adjoint matrices A, B, so that A, B, and A+B have prescribed spectrum. In this paper, we show how it connects to an interpolation problem for two variable real zero polynomials and a tracial moment problem. In addition, we outline an algorithm to construct a pair (A,B).


title: SPRK: A Low-Cost Stewart Platform For Motion Study In Surgical Robotics

abstract: To simulate body organ motion due to breathing, heart beats, or peristaltic movements, we designed a low-cost, miniaturized SPRK (Stewart Platform Research Kit) to translate and rotate phantom tissue. This platform is 20cm x 20cm x 10cm to fit in the workspace of a da Vinci Research Kit (DVRK) surgical robot and costs $250, two orders of magnitude less than a commercial Stewart platform. The platform has a range of motion of +/- 1.27 cm in translation along x, y, and z directions and has motion modes for sinusoidal motion and breathing-inspired motion. Modular platform mounts were also designed for pattern cutting and debridement experiments. The platform's positional controller has a time-constant of 0.2 seconds and the root-mean-square error is 1.22 mm, 1.07 mm, and 0.20 mm in x, y, and z directions respectively. All the details, CAD models, and control software for the platform is available at github.com/BerkeleyAutomation/sprk.


title: Superior valley polarization and coherence of 2s excitons in monolayer WSe2

abstract: We report experimental observation of 2s exciton radiative emission from monolayer tungsten diselenide, enabled by hexagonal boron nitride protected high-quality samples. The 2s luminescence is highly robust and persists up to 150K, offering a new quantum entity for manipulating the valley degree of freedom. Remarkably, the 2s exciton displays superior valley polarization and coherence than 1s under similar experimental conditions. This observation provides evidence that the Coulomb-exchange-interaction-driven valley-depolarization process, the Maialle-Silva-Sham mechanism, plays an important role in valley excitons of monolayer transition metal dichalcogenides.


title: Nonperturbative QCD thermodynamics in the external magnetic field

abstract: The thermodynamics of quarks and gluons strongly depends on the vacuum colormagnetic field, which grows with the temperature $T$, so that spatial string tension $\sigma_s ={\rm const}~ g^4 (T) T^2$. We investigate below what happens when one imposes in addition constant magnetic field and discover remarkable structure of the resulting thermodynamic potential.


title: Online Red Packets: A Large-scale Empirical Study of Gift Giving on WeChat

abstract: Gift giving is a ubiquitous social phenomenon, and red packets have been used as monetary gifts in Asian countries for thousands of years. In recent years, online red packets have become widespread in China through the WeChat platform. Exploiting a unique dataset consisting of 61 million group red packets and seven million users, we conduct a large-scale, data-driven study to understand the spread of red packets and the effect of red packets on group activity. We find that the cash flows between provinces are largely consistent with provincial GDP rankings, e.g., red packets are sent from users in the south to those in the north. By distinguishing spontaneous from reciprocal red packets, we reveal the behavioral patterns in sending red packets: males, seniors, and people with more in-group friends are more inclined to spontaneously send red packets, while red packets from females, youths, and people with less in-group friends are more reciprocal. Furthermore, we use propensity score matching to study the external effects of red packets on group dynamics. We show that red packets increase group participation and strengthen in-group relationships, which partly explain the benefits and motivations for sending red packets.


title: Microscopic Description of Electric and Magnetic Toroidal Multipoles in Hybrid Orbitals

abstract: We present a general formalism of multipole descriptions under the space-time inversion group. We elucidate that two types of atomic toroidal multipoles, i.e., electric and magnetic, are fundamental pieces to express electronic order parameters in addition to ordinary electric and magnetic multipoles. By deriving quantum-mechanical operators for both toroidal multipoles, we show that electric (magnetic) toroidal multipole higher than dipole (monopole) can become a primary order parameter in a hybridized-orbital system. We also demonstrate emergent cross-correlated couplings between electric, magnetic, and elastic degrees of freedom, such as magneto-electric and magneto(electro)-elastic couplings, under toroidal multipole orders.


title: A hyperbolic-equation system approach for magnetized electron fluids in quasi-neutral plasmas

abstract: A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional approaches using an elliptic equation (EE). A test calculation reveals that the HES approach can robustly solve problems of strong magnetic confinement by using an upwind method. The computation time of the HES approach is compared with that of the EE approach in terms of the size of the problem and the strength of magnetic confinement. The results indicate that the HES approach can be used to solve problems in a simple structured mesh without increasing computational time compared to the EE approach and that it features fast convergence in conditions of strong magnetic confinement.


title: Generation and amplification of high-order sideband induced by two-level atoms in a hybrid optomechanical system

abstract: It is quite important to enhance and control the optomechanically induced high-order sideband generation to achieve low-power optical comb and high-sensitivity sensing with an integratable structure. Here we present and analyze a proposal for enhancement and manipulation of optical nonlinearity and high-order sideband generation in a hybrid atom-cavity optomechanical system that is coherently driven by a bichromatic input field consisting of a control field and a probe field and works beyond the perturbative regime. Our numerical analysis with experimentally achievable parameters confirms that robust high-order sideband generation and typical spectral structures with non-perturbative features can be created even under weak driven fields. The dependence of the high-order sideband generation on the atomic parameters are also discussed in detail, including the decay rate of the atoms and the coupling parameter between the atoms and the cavity field. We show that the cutoff order as well as the amplitude of the higher order sidebands can be well tuned by the atomic coupling strength and the atomic decay rate. The proposed mechanism of enhancing optical nonlinearity is quite general and can be adopted to optomechanical systems with different types of cavity.