title: A scheme for determining fundamental interactions and the universality principle

abstract: The Standard Model of particle physics was established based on the equivalence principle and gauge invariance. The Lagrangians were built upon experimental data demonstrating the violation of discrete symmetries together with ideas of spontaneous symmetry breaking. The features of vector and axial-vector interactions and chirality principle are manually added on in order to fit in the observation data, rather than explained by the model as one may hope. Here we develop a theory of interactions based on the Dirac algebra together with the exterior operations the wedge product and the star operator on the space-time, and develop a scheme for determining all possible interaction formulations. In our scheme, the chirality transformation and the left and right handedness come up naturally. A simple postulate about the universality of physical interactions based on the space-time causal structure yields the attracting features mentioned above as consequences.


title: Entanglement and quantum state geometry of spin system with all-range Ising-type interaction

abstract: The evolution of $N$ spin-$1/2$ system with all-range Ising-type interaction is considered. For this system we study the entanglement of one spin with the rest spins. It is shown that the entanglement depends on the amount of spins and the initial state. Also the geometry of manifold which contains entangled states is obtained. Finally we find the dependence of entanglement on the scalar curvature of manifold and examine it for different number of spins in the system.


title: Survey on Semi-Explicit Time Integration of Eddy Current Problems

abstract: The spatial discretization of the magnetic vector potential formulation of magnetoquasistatic field problems results in an infinitely stiff differential-algebraic equation system. It is transformed into a finitely stiff ordinary differential equation system by applying a generalized Schur complement. Applying the explicit Euler time integration scheme to this system results in a small maximum stable time step size. Fast computations are required in every time step to yield an acceptable overall simulation time. Several acceleration methods are presented.


title: Relativistic disc line: a tool to constrain neutron star equation of state models

abstract: Relativistic iron K$\alpha$ spectral emission line from the inner disc of a neutron star low-mass X-ray binary (LMXB) was first detected in 2007. This discovery opened up new ways to probe strong gravity and dense matter. The past decade has seen detections of such a line from many neutron star LMXBs, and confirmation of this line from the same source with several X-ray satellites. These have firmly established the new field of relativistic disc line from neutron star systems in only ten years. Fitting the shape of such a line with an appropriate general relativistic model provides the accretion disc inner edge radius to the stellar mass ratio. In this review, we briefly discuss how an accurate measurement of this ratio with a future larger area X-ray instrument can be used to constrain neutron star equation of state models.


title: Accretion of Planetary Material onto Host Stars

abstract: Accretion of planetary material onto host stars may occur throughout a star's life. Especially prone to accretion, extrasolar planets in short-period orbits, while relatively rare, constitute a significant fraction of the known population, and these planets are subject to dynamical and atmospheric influences that can drive significant mass loss. Theoretical models frame expectations regarding the rates and extent of this planetary accretion. For instance, tidal interactions between planets and stars may drive complete orbital decay during the main sequence. Many planets that survive their stars' main sequence lifetime will still be engulfed when the host stars become red giant stars. There is some observational evidence supporting these predictions, such as a dearth of close-in planets around fast stellar rotators, which is consistent with tidal spin-up and planet accretion. There remains no clear chemical evidence for pollution of the atmospheres of main sequence or red giant stars by planetary materials, but a wealth of evidence points to active accretion by white dwarfs. In this article, we review the current understanding of accretion of planetary material, from the pre- to the post-main sequence and beyond. The review begins with the astrophysical framework for that process and then considers accretion during various phases of a host star's life, during which the details of accretion vary, and the observational evidence for accretion during these phases.


title: A Global Enhanced Vibrational Kinetic Model for Radio-Frequency Hydrogen Discharges and Application to the Simulation of a High Current Negative Hydrogen Ion Source

abstract: A Global Enhanced Vibrational Kinetic (GEVKM) model is presented for a new High Current Negative Hydrogen Ion Source (HCNHIS) developed by Busek Co. Inc. and Worcester Polytechnic Institute. The HCNHIS consists of a high-pressure radio-frequency discharge (RFD) chamber where high-lying vibrational states of the hydrogen molecules are produced, a bypass system, a nozzle, and a low-pressure negative hydrogen ion production region where $\text{H}^-$ are generated by the dissociative attachment of low energy electrons to rovibrationally excited hydrogen molecules. The GEVKM is developed from moment equations for multi-temperature chemically reacting plasmas derived from the Wang Chang-Uhlenbeck equations for a cylindrical geometry of an inductively coupled RFD chamber. The species included into the model are $\text{H}(n), n=0-3, \text{H}_2(v), v=0-14, \text{H}^+, \text{H}_2^+, \text{H}_3^+, \text{H}^-$, and $e$. The volume-averaged steady-state continuity equations coupled with the electron energy equation, the total energy equation and heat transfer to the chamber walls are solved simultaneously to obtain the volume-averaged number densities of the plasma components, the electron, heavy-particle and wall temperatures. The GEVKM is supplemented by a comprehensive set of surface and volumetric chemical processes governing vibrational and ionization kinetics of hydrogen plasmas. Analytic boundary conditions developed for the flow in the bypass tubes and in the nozzle take into account rarefaction effects and are validated by comparing pressures predictions in the RFD chamber with the pressure measurements. The GEVKM is verified and validated in the wide range of pressures and absorbed powers of different plasma reactors and is used in a parametric investigation of the RFD chamber of the HCNHIS with hydrogen inlet flow rates 5-1000 sccm and absorbed powers 200-1000 W.


title: The double obstacle problem on non divergence form

abstract: We study the regularity of the solution of the double obstacle problem form for fully non linear parabolic and elliptic operators. We show that when the obstacles are sufficiently regular the solution is $C^{1,\alpha}$ in the interior for both the parabolic and the elliptic cases.


title: A Unified Approach to the Global Exactness of Penalty and Augmented Lagrangian Functions I: Parametric Exactness

abstract: In this two-part study we develop a unified approach to the analysis of the global exactness of various penalty and augmented Lagrangian functions for finite-dimensional constrained optimization problems. This approach allows one to verify in a simple and straightforward manner whether a given penalty/augmented Lagrangian function is exact, i.e. whether the problem of unconstrained minimization of this function is equivalent (in some sense) to the original constrained problem, provided the penalty parameter is sufficiently large. Our approach is based on the so-called localization principle that reduces the study of global exactness to a local analysis of a chosen merit function near globally optimal solutions. In turn, such local analysis can usually be performed with the use of sufficient optimality conditions and constraint qualifications. In the first paper we introduce the concept of global parametric exactness and derive the localization principle in the parametric form. With the use of this version of the localization principle we recover existing simple necessary and sufficient conditions for the global exactness of linear penalty functions, and for the existence of augmented Lagrange multipliers of Rockafellar-Wets' augmented Lagrangian. Also, we obtain completely new necessary and sufficient conditions for the global exactness of general nonlinear penalty functions, and for the global exactness of a continuously differentiable penalty function for nonlinear second-order cone programming problems. We briefly discuss how one can construct a continuously differentiable exact penalty function for nonlinear semidefinite programming problems, as well.


title: Etude et mise en oeuvre des techniques Temps-Distance en Heliosismologie

abstract: In the same way that seismologists study the interior of the earth from the waves generated by earthquakes, heliosismology is the science that is interested in the study of oscillatory waves inside the Sun. Temperature, chemical composition, rotational velocities at different depths are all factors that influence the oscillation frequencies of the waves trapped inside the Sun. The observation and measurement of the frequencies (or velocities) of oscillations of the principal modes on the surface of the Sun allow us to study the internal properties of the star. In our study, we use a local helioseismic method, called "Time-Distance"; the latter, after appropriate treatment (remapping, tracking, and filtering) of the data, which carried out from the GONG (Global Oscillation Network Group) observation network, allows us to deduce, by correlation of the observed signals, the relationship between the time of the wave path and its crossed distance (hence the name of the method) between different points on the surface as a function of the different solar sub-surface parameters, via the dispersion relation. Once the travel times have been established by approximation of the correlation function by a Gaussian wave packet and based on the Fermat principle dealing with acoustic waves, we find the relationship between these times and the internal parameters of the crossed medium. Two models are considered in this study. The first does not include the effects of the magnetic field and takes into account the flow velocity and the sound speed. The second, on the other hand, includes the magnetic effects, by means of the Alfven velocity (magnetic field perturbation velocity), in addition to the parameters mentioned above. The last part of this paper was devoted to the inversion of the obtained times, which allows us to go back to the solar internal parameters.


title: A QM/MM equation-of-motion coupled-cluster approach for predicting semiconductor color-center structure and emission frequencies

abstract: Valence excitation spectra are computed for all deep-center silicon-vacancy defect types in 3C, 4H, and 6H silicon carbide (SiC) and comparisons are made with literature photoluminescence measurements. Nuclear geometries surrounding the defect centers are optimized within a Gaussian basis-set framework using many-body perturbation theory or density functional theory (DFT) methods, with computational expenses minimized by a QM/MM technique called SIMOMM. Vertical excitation energies are subsequently obtained by applying excitation-energy, electron-attached, and ionized equation-of-motion coupled-cluster (EOMCC) methods, where appropriate, as well as time-dependent (TD) DFT, to small models including only a few atoms adjacent to the defect center. We consider the relative quality of various EOMCC and TD-DFT methods for (i) energy-ordering potential ground states differing incrementally in charge and multiplicity, (ii) accurately reproducing experimentally measured photoluminescence peaks, and (iii) energy-ordering defects of different types occurring within a given polytype. The extensibility of this approach to transition-metal defects is also tested by applying it to silicon-substitutional chromium defects in SiC and comparing with measurements. It is demonstrated that, when used in conjunction with SIMOMM-optimized geometries, EOMCC-based methods can provide a reliable prediction of the ground-state charge and multiplicity, while also giving a quantitative description of the photoluminescence spectra, accurate to within 0.1 eV of measurement in all cases considered.


title: Rayleigh scattering in dense fluid helium

abstract: Iglesias et al. (2002) showed that the Rayleigh scattering from helium atoms decreases by collective effects in the atmospheres of cool white dwarf stars. Their study is here extended to consider an accurate evaluation of the atomic polarizability and the density effects involved in the Rayleigh cross section over a wide density-temperature region. The dynamic dipole polarizability of helium atoms in the ground state is determinated with the oscillator-strength distribution approach. The spectral density of oscillator strength considered includes most significant single and doubly excited transitions to discrete and continuum energies. Static and dynamic polarizability results are confronted with experiments and other theoretical evaluations shown a very good agreement. In addition, the refractive index of helium is evaluated with the Lorentz-Lorenz equation and shows a satisfactory agreement with the most recent experiments. The effect of spatial correlation of atoms on the Rayleigh scattering is calculated with Monte Carlo simulations and effective energy potentials that represent the particle interactions, covering fluid densities between 0.005 and a few g/cm$^3$ and temperatures between $1000$ K and $15000$ K. We provide analytical fits from which the Rayleigh cross section of fluid helium can be easily calculated at wavelength $\lambda>505.35$ \AA. Collision-induced light scattering was estimated to be the dominant scattering process at densities greater than 1-2 g/cm$^3$ depending on the temperature.


title: Estimated Depth Map Helps Image Classification

abstract: We consider image classification with estimated depth. This problem falls into the domain of transfer learning, since we are using a model trained on a set of depth images to generate depth maps (additional features) for use in another classification problem using another disjoint set of images. It's challenging as no direct depth information is provided. Though depth estimation has been well studied, none have attempted to aid image classification with estimated depth. Therefore, we present a way of transferring domain knowledge on depth estimation to a separate image classification task over a disjoint set of train, and test data. We build a RGBD dataset based on RGB dataset and do image classification on it. Then evaluation the performance of neural networks on the RGBD dataset compared to the RGB dataset. From our experiments, the benefit is significant with shallow and deep networks. It improves ResNet-20 by 0.55% and ResNet-56 by 0.53%. Our code and dataset are available publicly.


title: Max-Min Fair Resource Allocation in Millimetre-Wave Backhauls

abstract: 5G mobile networks are expected to provide pervasive high speed wireless connectivity, to support increasingly resource intensive user applications. Network hyper-densification therefore becomes necessary, though connecting to the Internet tens of thousands of base stations is non-trivial, especially in urban scenarios where optical fibre is difficult and costly to deploy. The millimetre wave (mm-wave) spectrum is a promising candidate for inexpensive multi-Gbps wireless backhauling, but exploiting this band for effective multi-hop data communications is challenging. In particular, resource allocation and scheduling of very narrow transmission/ reception beams requires to overcome terminal deafness and link blockage problems, while managing fairness issues that arise when flows encounter dissimilar competition and traverse different numbers of links with heterogeneous quality. In this paper, we propose WiHaul, an airtime allocation and scheduling mechanism that overcomes these challenges specific to multi-hop mm-wave networks, guarantees max-min fairness among traffic flows, and ensures the overall available backhaul resources are fully utilised. We evaluate the proposed WiHaul scheme over a broad range of practical network conditions, and demonstrate up to 5 times individual throughput gains and a fivefold improvement in terms of measurable fairness, over recent mm-wave scheduling solutions.


title: Search for Nonstandard Neutrino Interactions with IceCube DeepCore

abstract: As atmospheric neutrinos propagate through the Earth, vacuum-like oscillations are modified by Standard-Model neutral- and charged-current interactions with electrons. Theories beyond the Standard Model introduce heavy, TeV-scale bosons that can produce nonstandard neutrino interactions. These additional interactions may modify the Standard Model matter effect producing a measurable deviation from the prediction for atmospheric neutrino oscillations. The result described in this paper constrains nonstandard interaction parameters, building upon a previous analysis of atmospheric muon-neutrino disappearance with three years of IceCube-DeepCore data. The best fit for the muon to tau flavor changing term is $\epsilon_{\mu \tau}=-0.0005$, with a 90\% C.L. allowed range of $-0.0067 <\epsilon_{\mu \tau}< 0.0081$. This result is more restrictive than recent limits from other experiments for $\epsilon_{\mu \tau}$. Furthermore, our result is complementary to a recent constraint on $\epsilon_{\mu \tau}$ using another publicly available IceCube high-energy event selection. Together, they constitute the world's best limits on nonstandard interactions in the $\mu-\tau$ sector.


title: A Deep-Reinforcement Learning Approach for Software-Defined Networking Routing Optimization

abstract: In this paper we design and evaluate a Deep-Reinforcement Learning agent that optimizes routing. Our agent adapts automatically to current traffic conditions and proposes tailored configurations that attempt to minimize the network delay. Experiments show very promising performance. Moreover, this approach provides important operational advantages with respect to traditional optimization algorithms.


title: Performance of the diamond active target prototype for the PADME experiment at the DA$\Phi$NE BTF

abstract: The PADME experiment at the DA$\Phi$NE Beam-Test Facility (BTF) is designed to search for the gauge boson of a new $\rm U(1)$ interaction in the process e$^+$e$^-\rightarrow\gamma$+$\rm A'$, using the intense positron beam hitting a light target. The $\rm A'$, usually referred as dark photon, is assumed to decay into invisible particles of a secluded sector and it can be observed by searching for an anomalous peak in the spectrum of the missing mass measured in events with a single photon in the final state. The measurement requires the determination of the 4-momentum of the recoil photon, performed by a homogeneous, highly segmented BGO crystals calorimeter. A significant improvement of the missing mass resolution is possible using an active target capable to determine the average position of the positron bunch with a resolution of less than 1 mm. This report presents the performance of a real size $\rm (2x2 cm^2)$ PADME active target made of a thin (50 $\mu$m) diamond sensor, with graphitic strips produced via laser irradiation on both sides. The measurements are based on data collected in a beam test at the BTF in November 2015.


title: Machine learning of molecular properties: locality and active learning

abstract: In recent years the machine learning techniques have shown a great potential in various problems from a multitude of disciplines, including materials design and drug discovery. The high computational speed on the one hand and the accuracy comparable to that of DFT on another hand make machine learning algorithms efficient for high-throughput screening through chemical and configurational space. However, the machine learning algorithms available in the literature require large training datasets to reach the chemical accuracy and also show large errors for the so-called outliers - the out-of-sample molecules, not well-represented in the training set. In the present paper we propose a new machine learning algorithm for predicting molecular properties that addresses these two issues: it is based on a local model of interatomic interactions providing high accuracy when trained on relatively small training sets and an active learning algorithm of optimally choosing the training set that significantly reduces the errors for the outliers. We compare our model to the other state-of-the-art algorithms from the literature on the widely used benchmark tests.


title: On recognizing shapes of polytopes from their shadows

abstract: Let $P$ and $Q$ be two convex polytopes both contained in the interior of an Euclidean ball $r\textbf{B}^{d}$. We prove that $P=Q$ provided that their sight cones from any point on the sphere $rS^{d-1}$ are congruent. We also prove an analogous result for spherical projections.


title: Solving the flatness problem with an anisotropic instanton in Ho\v{r}ava-Lifshitz gravity

abstract: In Ho\v{r}ava-Lifshitz gravity a scaling isotropic in space but anisotropic in spacetime, often called anisotropic scaling with the dynamical critical exponent z=3, lies at the base of its renormalizability. This scaling also leads to a novel mechanism of generating scale-invariant cosmological perturbations, solving the horizon problem without inflation. In this paper we propose a possible solution to the flatness problem, in which we assume that the initial condition of the Universe is set by a small instanton respecting the same scaling. We argue that the mechanism may be more general than the concrete model presented here, and rely simply on the deformed dispersion relations of the theory, and on equipartition of the various forms of energy at the starting point.


title: A Flocking-based Approach for Distributed Stochastic Optimization

abstract: In recent years, the paradigm of cloud computing has emerged as an architecture for computing that makes use of distributed (networked) computing resources. In this paper, we consider a distributed computing algorithmic scheme for stochastic optimization which relies on modest communication requirements amongst processors and most importantly, does not require synchronization. Specifically, we analyze a scheme with $N>1$ independent threads implementing each a stochastic gradient algorithm. The threads are coupled via a perturbation of the gradient (with attractive and repulsive forces) in a similar manner to mathematical models of flocking, swarming and other group formations found in nature with mild communication requirements. When the objective function is convex, we show that a flocking-like approach for distributed stochastic optimization provides a noise reduction effect similar to that of a centralized stochastic gradient algorithm based upon the average of $N$ gradient samples at each step. The distributed nature of flocking makes it an appealing computational alternative. We show that when the overhead related to the time needed to gather $N$ samples and synchronization is not negligible, the flocking implementation outperforms a centralized stochastic gradient algorithm based upon the average of $N$ gradient samples at each step. When the objective function is not convex, the flocking-based approach seems better suited to escape locally optimal solutions due to the repulsive force which enforces a certain level of diversity in the set of candidate solutions. Here again, we show that the noise reduction effect is similar to that associated to the centralized stochastic gradient algorithm based upon the average of $N$ gradient samples at each step.


title: (m,n)-Quasitilted and (m,n)-Almost Hereditary Algebras

abstract: Motivated by the study of (m,n)-quasitilted algebras, which are the piecewise hereditary algebras obtained from quasitilted algebras of global dimension two by a sequence of (co)tiltings involving n-1 tilting modules and m-1 cotilting modules, we introduce (m,n)-almost hereditary algebras. These are the algebras with global dimension m+n and such that any indecomposable module has projective dimension at most m, or else injective dimension at most n. We relate these two classes of algebras, among which (m,1)-almost hereditary ones play a special role. For these, we prove that any indecomposable module lies in the right part of the module category, or else in an m-analog of the left part. This is based on the more general study of algebras the module categories of which admit a torsion-free subcategory such that any indecomposable module lies in that subcategory, or else has injective dimension at most n.


title: Equivalence of cosmological observables in conformally related scalar tensor theories

abstract: Scalar tensor theories can be expressed in different frames, such as the commonly-used Einstein and Jordan frames, and it is generally accepted that cosmological observables are the same in these frames. We revisit this by making a detailed side-by-side comparison of the quantities and equations in two conformally related frames, from the actions and fully covariant field equations to the linearised equations in both real and Fourier spaces. This confirms that the field and conservation equations are equivalent in the two frames, in the sense that we can always re-express equations in one frame using relevant transformations of variables to derive the corresponding equations in the other. We show, with both analytical derivation and a numerical example, that the line-of-sight integration to calculate CMB temperature anisotropies can be done using either Einstein frame or Jordan frame quantities, and the results are identical, provided the correct redshift is used in the Einstein frame ($1+z\neq1/a$).


title: Two-dimensional quantum droplets in dipolar Bose gases

abstract: We calculate analytically the quantum and thermal fluctuations corrections of a dilute quasi-two-dimensional Bose-condensed dipolar gas. We show that these fluctuations may change their character from repulsion to attraction in the density-temperature plane owing to the striking momentum dependence of the dipole-dipole interactions. The dipolar instability is halted by such unconventional beyond mean field corrections leading to the formation of a droplet phase. The equilibrium features and coherence properties exhibited by such droplets are deeply discussed. At finite temperature, we find that the equilibrium density crucially depends on the temperature and on the confinement strength and thus, a stable droplet can exist only at ultralow temperature due to the strong thermal fluctuations.


title: On the Design of LQR Kernels for Efficient Controller Learning

abstract: Finding optimal feedback controllers for nonlinear dynamic systems from data is hard. Recently, Bayesian optimization (BO) has been proposed as a powerful framework for direct controller tuning from experimental trials. For selecting the next query point and finding the global optimum, BO relies on a probabilistic description of the latent objective function, typically a Gaussian process (GP). As is shown herein, GPs with a common kernel choice can, however, lead to poor learning outcomes on standard quadratic control problems. For a first-order system, we construct two kernels that specifically leverage the structure of the well-known Linear Quadratic Regulator (LQR), yet retain the flexibility of Bayesian nonparametric learning. Simulations of uncertain linear and nonlinear systems demonstrate that the LQR kernels yield superior learning performance.


title: Photoassisted physical vapor epitaxial growth of semiconductors: a review of light-induced modifications to growth processes

abstract: Herein we review the remarkable range of modifications to materials properties associated with photoexcitation of the growth surface during physical vapor epitaxy of semiconductors. We concentrate on mechanisms producing measureable, utilizable changes in crystalline perfection, phase, composition, doping, and defect distribution. We outline the relevant physics of different mechanisms, concentrating on those yielding effects orthogonal to the primary growth variables of temperature and atomic or molecular fluxes and document the phenomenological effects reported. Based on experimental observations from a range of semiconductor systems and growth conditions, the primary effects include enhanced anion desorption, molecular dissociation, increased doping efficiency, modification to defect populations and improvements to the crystalline quality of epilayers grown at low temperatures. Future research directions and technological applications are also discussed.


title: Axion as a non-WIMP dark matter candidate

abstract: The axion arises in well-motivated extensions of the Standard Model of particle physics and is regarded as an alternative to the weakly interacting massive particle paradigm to explain the nature of dark matter. In this contribution, we review theoretical aspects of dark matter axions, particularly focusing on recent developments in the estimation of their relic abundance. A closer look at their non-thermal production mechanisms in the early universe reveals the possibility of explaining the observed dark matter abundance in various mass ranges. The mass ranges predicted in various cosmological scenarios are briefly summarized.


title: On Compiling DNNFs without Determinism

abstract: State-of-the-art knowledge compilers generate deterministic subsets of DNNF, which have been recently shown to be exponentially less succinct than DNNF. In this paper, we propose a new method to compile DNNFs without enforcing determinism necessarily. Our approach is based on compiling deterministic DNNFs with the addition of auxiliary variables to the input formula. These variables are then existentially quantified from the deterministic structure in linear time, which would lead to a DNNF that is equivalent to the input formula and not necessarily deterministic. On the theoretical side, we show that the new method could generate exponentially smaller DNNFs than deterministic ones, even by adding a single auxiliary variable. Further, we show that various existing techniques that introduce auxiliary variables to the input formulas can be employed in our framework. On the practical side, we empirically demonstrate that our new method can significantly advance DNNF compilation on certain benchmarks.


title: Near Optimal Sketching of Low-Rank Tensor Regression

abstract: We study the least squares regression problem \begin{align*} \min_{\Theta \in \mathcal{S}_{\odot D,R}} \|A\Theta-b\|_2, \end{align*} where $\mathcal{S}_{\odot D,R}$ is the set of $\Theta$ for which $\Theta = \sum_{r=1}^{R} \theta_1^{(r)} \circ \cdots \circ \theta_D^{(r)}$ for vectors $\theta_d^{(r)} \in \mathbb{R}^{p_d}$ for all $r \in [R]$ and $d \in [D]$, and $\circ$ denotes the outer product of vectors. That is, $\Theta$ is a low-dimensional, low-rank tensor. This is motivated by the fact that the number of parameters in $\Theta$ is only $R \cdot \sum_{d=1}^D p_d$, which is significantly smaller than the $\prod_{d=1}^{D} p_d$ number of parameters in ordinary least squares regression. We consider the above CP decomposition model of tensors $\Theta$, as well as the Tucker decomposition. For both models we show how to apply data dimensionality reduction techniques based on {\it sparse} random projections $\Phi \in \mathbb{R}^{m \times n}$, with $m \ll n$, to reduce the problem to a much smaller problem $\min_{\Theta} \|\Phi A \Theta - \Phi b\|_2$, for which if $\Theta'$ is a near-optimum to the smaller problem, then it is also a near optimum to the original problem. We obtain significantly smaller dimension and sparsity in $\Phi$ than is possible for ordinary least squares regression, and we also provide a number of numerical simulations supporting our theory.


title: Decomposing GR(1) Games with Singleton Liveness Guarantees for Efficient Synthesis

abstract: Temporal logic based synthesis approaches are often used to find trajectories that are correct-by-construction for tasks in systems with complex behavior. Some examples of such tasks include synchronization for multi-agent hybrid systems, reactive motion planning for robots. However, the scalability of such approaches is of concern and at times a bottleneck when transitioning from theory to practice. In this paper, we identify a class of problems in the GR(1) fragment of linear-time temporal logic (LTL) where the synthesis problem allows for a decomposition that enables easy parallelization. This decomposition also reduces the alternation depth, resulting in more efficient synthesis. A multi-agent robot gridworld example with coordination tasks is presented to demonstrate the application of the developed ideas and also to perform empirical analysis for benchmarking the decomposition-based synthesis approach.


title: Practical Machine Learning for Cloud Intrusion Detection: Challenges and the Way Forward

abstract: Operationalizing machine learning based security detections is extremely challenging, especially in a continuously evolving cloud environment. Conventional anomaly detection does not produce satisfactory results for analysts that are investigating security incidents in the cloud. Model evaluation alone presents its own set of problems due to a lack of benchmark datasets. When deploying these detections, we must deal with model compliance, localization, and data silo issues, among many others. We pose the problem of "attack disruption" as a way forward in the security data science space. In this paper, we describe the framework, challenges, and open questions surrounding the successful operationalization of machine learning based security detections in a cloud environment and provide some insights on how we have addressed them.


title: Covert Wireless Communication with Artificial Noise Generation

abstract: Covert communication conceals the transmission of the message from an attentive adversary. Recent work on the limits of covert communication in additive white Gaussian noise (AWGN) channels has demonstrated that a covert transmitter (Alice) can reliably transmit a maximum of $\mathcal{O}\left(\sqrt{n}\right)$ bits to a covert receiver (Bob) without being detected by an adversary (Warden Willie) in $n$ channel uses. This paper focuses on the scenario where other friendly nodes distributed according to a two-dimensional Poisson point process with density $m$ are present in the environment. We propose a strategy where the friendly node closest to the adversary, without close coordination with Alice, produces artificial noise. We show that this method allows Alice to reliably and covertly send $\mathcal{O}(\min\{{n,m^{\gamma/2}\sqrt{n}}\})$ bits to Bob in $n$ channel uses, where $\gamma$ is the path-loss exponent. Moreover, we also consider a setting where there are $N_{\mathrm{w}}$ collaborating adversaries uniformly and randomly located in the environment and show that in $n$ channel uses, Alice can reliably and covertly send $\mathcal{O}\left(\min\left\{n,\frac{m^{\gamma/2} \sqrt{n}}{N_{\mathrm{w}}^{\gamma}}\right\}\right)$ bits to Bob when $\gamma >2$, and $\mathcal{O}\left(\min\left\{n,\frac{m \sqrt{n}}{N_{\mathrm{w}}^{2}\log^2 {N_{\mathrm{w}}}}\right\}\right)$ when $\gamma = 2$. Conversely, we demonstrate that no higher covert throughput is possible for $\gamma>2$.


title: Persistence Flamelets: multiscale Persistent Homology for kernel density exploration

abstract: In recent years there has been noticeable interest in the study of the "shape of data". Among the many ways a "shape" could be defined, topology is the most general one, as it describes an object in terms of its connectivity structure: connected components (topological features of dimension 0), cycles (features of dimension 1) and so on. There is a growing number of techniques, generally denoted as Topological Data Analysis, aimed at estimating topological invariants of a fixed object; when we allow this object to change, however, little has been done to investigate the evolution in its topology. In this work we define the Persistence Flamelets, a multiscale version of one of the most popular tool in TDA, the Persistence Landscape. We examine its theoretical properties and we show how it could be used to gain insights on KDEs bandwidth parameter.


title: Talagrand Concentration Inequalities for Stochastic Partial Differential Equations

abstract: One way to define the concentration of measure phenomenon is via Talagrand inequalities, also called transportation-information inequalities. That is, a comparison of the Wasserstein distance from the given measure to any other absolutely continuous measure with finite relative entropy. Such transportation-information inequalities were recently established for some stochastic differential equations. Here, we develop a similar theory for some stochastic partial differential equations.


title: Emulating galaxy clustering and galaxy-galaxy lensing into the deeply nonlinear regime: methodology, information, and forecasts

abstract: The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $\xi_{\text{gg}}(r)$ and galaxy-matter cross-correlation $\xi_{\text{gm}}(r)$ as a function of $\sigma_8$, $\Omega_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $\sigma_8\Omega_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $\gamma_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $\sigma_8\Omega_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.


title: Supervised Learning with Indefinite Topological Kernels

abstract: Topological Data Analysis (TDA) is a recent and growing branch of statistics devoted to the study of the shape of the data. In this work we investigate the predictive power of TDA in the context of supervised learning. Since topological summaries, most noticeably the Persistence Diagram, are typically defined in complex spaces, we adopt a kernel approach to translate them into more familiar vector spaces. We define a topological exponential kernel, we characterize it, and we show that, despite not being positive semi-definite, it can be successfully used in regression and classification tasks.


title: FairFuzz: Targeting Rare Branches to Rapidly Increase Greybox Fuzz Testing Coverage

abstract: In recent years, fuzz testing has proven itself to be one of the most effective techniques for finding correctness bugs and security vulnerabilities in practice. One particular fuzz testing tool, American Fuzzy Lop or AFL, has become popular thanks to its ease-of-use and bug-finding power. However, AFL remains limited in the depth of program coverage it achieves, in particular because it does not consider which parts of program inputs should not be mutated in order to maintain deep program coverage. We propose an approach, FairFuzz, that helps alleviate this limitation in two key steps. First, FairFuzz automatically prioritizes inputs exercising rare parts of the program under test. Second, it automatically adjusts the mutation of inputs so that the mutated inputs are more likely to exercise these same rare parts of the program. We conduct evaluation on real-world programs against state-of-the-art versions of AFL, thoroughly repeating experiments to get good measures of variability. We find that on certain benchmarks FairFuzz shows significant coverage increases after 24 hours compared to state-of-the-art versions of AFL, while on others it achieves high program coverage at a significantly faster rate.


title: Automatic Detection of Malware-Generated Domains with Recurrent Neural Models

abstract: Modern malware families often rely on domain-generation algorithms (DGAs) to determine rendezvous points to their command-and-control server. Traditional defence strategies (such as blacklisting domains or IP addresses) are inadequate against such techniques due to the large and continuously changing list of domains produced by these algorithms. This paper demonstrates that a machine learning approach based on recurrent neural networks is able to detect domain names generated by DGAs with high precision. The neural models are estimated on a large training set of domains generated by various malwares. Experimental results show that this data-driven approach can detect malware-generated domain names with a F_1 score of 0.971. To put it differently, the model can automatically detect 93 % of malware-generated domain names for a false positive rate of 1:100.


title: Automated Adaptive Optics

abstract: Large area surveys will dominate the forthcoming decades of astronomy and their success requires characterizing thousands of discoveries through additional observations at higher spatial or spectral resolution, and at complementary cadences or periods. Only the full automation of adaptive optics systems will enable high-acuity, high-sensitivity follow-up observations of several tens of thousands of these objects per year, maximizing on-sky time. Automation will also enable rapid response to target-of-opportunity events within minutes, minimizing the time between discovery and characterization. In June 2012, we demonstrated the first fully automated operation of an astronomical adaptive optics system by observing 125 objects in succession with the Robo-AO system. Efficiency has increased ever since, with a typical night comprising 200-250 automated observations at the visible diffraction limit. By observing tens of thousands of targets in the largest-ever adaptive-optics surveys, Robo-AO has demonstrated the ability to address the follow-up needs of current and future large astronomical surveys.


title: On the Use of Machine Translation-Based Approaches for Vietnamese Diacritic Restoration

abstract: This paper presents an empirical study of two machine translation-based approaches for Vietnamese diacritic restoration problem, including phrase-based and neural-based machine translation models. This is the first work that applies neural-based machine translation method to this problem and gives a thorough comparison to the phrase-based machine translation method which is the current state-of-the-art method for this problem. On a large dataset, the phrase-based approach has an accuracy of 97.32% while that of the neural-based approach is 96.15%. While the neural-based method has a slightly lower accuracy, it is about twice faster than the phrase-based method in terms of inference speed. Moreover, neural-based machine translation method has much room for future improvement such as incorporating pre-trained word embeddings and collecting more training data.


title: High Coherence Mid-Infrared Dual Comb Spectroscopy Spanning 2.6 to 5.2 microns

abstract: Mid-infrared dual-comb spectroscopy has the potential to supplant conventional high-resolution Fourier transform spectroscopy in applications that require high resolution, accuracy, signal-to-noise ratio, and speed. Until now, dual-comb spectroscopy in the mid-infrared has been limited to narrow optical bandwidths or to low signal-to-noise ratios. Using a combination of digital signal processing and broadband frequency conversion in waveguides, we demonstrate a mid-infrared dual-comb spectrometer that can measure comb-tooth resolved spectra across an octave of bandwidth in the mid-infrared from 2.6-5.2 $\mu$m with sub-MHz frequency precision and accuracy and with a spectral signal-to-noise ratio as high as 6500. As a demonstration, we measure the highly structured, broadband cross-section of propane (C3H8) in the 2860-3020 cm-1 region, the complex phase/amplitude spectrum of carbonyl sulfide (COS) in the 2000 to 2100 cm-1 region, and the complex spectra of methane, acetylene, and ethane in the 2860-3400 cm-1 region.


title: Magnetic Fields of Neutron Stars

abstract: This article briefly reviews our current understanding of the evolution of magnetic fields in neutron stars, which basically defines the evolutionary pathways between different observational classes of neutron stars. The emphasis here is on the evolution in binary systems and the newly emergent classes of millisecond pulsars.


title: Statistical methods for estimating ecological breakpoints and prediction intervals

abstract: The relationships among ecological variables are usually obtained by fitting statistical models that go through the conditional means of the dependent variables. For example, the nonparametric loess and the parametric piecewise linear regression models, which pass through the conditional mean of the response variable given the predictor, are used to analyze simple to complex relationships among variables. We used loess and bootstrapped confidence interval to subjectively identify the number and positions of potential ecological breakpoints in a bivariate relationship, and a piecewise linear regression model (PLRM) to quantitatively estimate the location of breakpoints and the associated precision. We also estimated breakpoint location and precision using a piecewise linear quantile regression model (PQRM), which is fitted to the quantiles of the conditional distribution of the response variable given the predictor and provides much richer information in terms of estimating relationships and breakpoints. We compared the precision of breakpoints estimated by PQRM relative to PLRM. We compared the precision of the methods using two examples from the ecological literature suspected to exhibit multiple breakpoints: relating a Fish Index of Biotic Integrity (an index of wetlands' fish community 'health') to the amount of human activity in wetlands' adjacent watersheds; and relating the biomass of cyanobacteria to the total phosphorus concentration in Canadian lakes. Statistically significant breakpoints were detected for both datasets, demarcating the boundaries of three line segments with markedly different slopes. We recommend the piecewise linear quantile regression as an effective means of characterizing bivariate environmental relationships where the scatter of points represents natural environmental variation rather than measurement error.


title: Spectral signatures of many-body localization with interacting photons

abstract: Statistical mechanics is founded on the assumption that a system can reach thermal equilibrium, regardless of the starting state. Interactions between particles facilitate thermalization, but, can interacting systems always equilibrate regardless of parameter values\,? The energy spectrum of a system can answer this question and reveal the nature of the underlying phases. However, most experimental techniques only indirectly probe the many-body energy spectrum. Using a chain of nine superconducting qubits, we implement a novel technique for directly resolving the energy levels of interacting photons. We benchmark this method by capturing the intricate energy spectrum predicted for 2D electrons in a magnetic field, the Hofstadter butterfly. By increasing disorder, the spatial extent of energy eigenstates at the edge of the energy band shrink, suggesting the formation of a mobility edge. At strong disorder, the energy levels cease to repel one another and their statistics approaches a Poisson distribution - the hallmark of transition from the thermalized to the many-body localized phase. Our work introduces a new many-body spectroscopy technique to study quantum phases of matter.


title: Deconvolutional Latent-Variable Model for Text Sequence Matching

abstract: A latent-variable model is introduced for text matching, inferring sentence representations by jointly optimizing generative and discriminative objectives. To alleviate typical optimization challenges in latent-variable models for text, we employ deconvolutional networks as the sequence decoder (generator), providing learned latent codes with more semantic information and better generalization. Our model, trained in an unsupervised manner, yields stronger empirical predictive performance than a decoder based on Long Short-Term Memory (LSTM), with less parameters and considerably faster training. Further, we apply it to text sequence-matching problems. The proposed model significantly outperforms several strong sentence-encoding baselines, especially in the semi-supervised setting.


title: LHC limits on axion-like particles from heavy-ion collisions

abstract: In these proceedings we use recent LHC heavy-ion data to set a limit on axion-like particles coupling to electromagnetism with mass in the range 10-100 GeV. We recast ATLAS data as per the strategy proposed in 1607.06083, and find results in-line with the projections given there.


title: Fidelity Susceptibility in the Quantum Rabi Model

abstract: Quantum criticality usually occurs in many-body systems. Recently it was shown that the quantum Rabi model, which describes a two-level atom coupled to a single model cavity field, presents quantum phase transitions from a normal phase to a superradiate phase when the ratio between the frequency of the two level atom and the frequency of the cavity field extends to infinity. In this work, we study quantum phase transitions in the quantum Rabi model from the fidelity susceptibility perspective. We found that the fidelity susceptibility and the generalized adiabatic susceptibility present universal finite size scaling behaviors near the quantum critical point of the Rabi model if the ratio between frequency of the two level atom and frequency of the cavity field is finite. From the finite size scaling analysis of the fidelity susceptibility, we found that the adiabatic dimension of the fidelity susceptibility and the generalized adiabatic susceptibility of fourth order in the Rabi model are $4/3$ and $2$, respectively. Meanwhile the correlation length critical exponent and the dynamical critical exponent in the quantum critical point of the Rabi model are found to be $3/2$ and $1/3$ respectively. Since the fidelity susceptibility and the generalized adiabatic susceptibility are the moments of the quantum noise spectrum which is directly measurable by experiments in linear response regime, the scaling behavior of the fidelity susceptibility in the Rabi model could be tested experimentally. The simple structure of the quantum Rabi model paves the way for experimentally observing the universal scaling behavior of the fidelity susceptibility at a quantum phase transition.


title: Cooperative Adaptive Control for Cloud-Based Robotics

abstract: This paper studies collaboration through the cloud in the context of cooperative adaptive control for robot manipulators. We first consider the case of multiple robots manipulating a common object through synchronous centralized update laws to identify unknown inertial parameters. Through this development, we introduce a notion of Collective Sufficient Richness, wherein parameter convergence can be enabled through teamwork in the group. The introduction of this property and the analysis of stable adaptive controllers that benefit from it constitute the main new contributions of this work. Building on this original example, we then consider decentralized update laws, time-varying network topologies, and the influence of communication delays on this process. Perhaps surprisingly, these nonidealized networked conditions inherit the same benefits of convergence being determined through collective effects for the group. Simple simulations of a planar manipulator identifying an unknown load are provided to illustrate the central idea and benefits of Collective Sufficient Richness.


title: A Detached Protostellar Disk around a $\sim$0.2$M_{\odot}$ protostar in a Possible Site of a Multiple Star Formation in a Dynamical Environment in Taurus

abstract: We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.


title: Cost Adaptation for Robust Decentralized Swarm Behaviour

abstract: Decentralized receding horizon control (D-RHC) provides a mechanism for coordination in multi-agent settings without a centralized command center. However, combining a set of different goals, costs, and constraints to form an efficient optimization objective for D-RHC can be difficult. To allay this problem, we use a meta-learning process -- cost adaptation -- which generates the optimization objective for D-RHC to solve based on a set of human-generated priors (cost and constraint functions) and an auxiliary heuristic. We use this adaptive D-RHC method for control of mesh-networked swarm agents. This formulation allows a wide range of tasks to be encoded and can account for network delays, heterogeneous capabilities, and increasingly large swarms through the adaptation mechanism. We leverage the Unity3D game engine to build a simulator capable of introducing artificial networking failures and delays in the swarm. Using the simulator we validate our method on an example coordinated exploration task. We demonstrate that cost adaptation allows for more efficient and safer task completion under varying environment conditions and increasingly large swarm sizes. We release our simulator and code to the community for future work.


title: Steady double vortex patches with opposite signs in a planar ideal fluid

abstract: In this paper we consider steady vortex flows for the incompressible Euler equations in a planar bounded domain. By solving a variational problem for the vorticity, we construct steady double vortex patches with opposite signs concentrating at a strict local minimum point of the Kirchhoff-Routh function with $k=2$. Moreover, we show that such steady solutions are in fact local maximizers of the kinetic energy among isovortical patches, which correlates stability to uniqueness.


title: Variational Memory Addressing in Generative Models

abstract: Aiming to augment generative models with external memory, we interpret the output of a memory module with stochastic addressing as a conditional mixture distribution, where a read operation corresponds to sampling a discrete memory address and retrieving the corresponding content from memory. This perspective allows us to apply variational inference to memory addressing, which enables effective training of the memory module by using the target information to guide memory lookups. Stochastic addressing is particularly well-suited for generative models as it naturally encourages multimodality which is a prominent aspect of most high-dimensional datasets. Treating the chosen address as a latent variable also allows us to quantify the amount of information gained with a memory lookup and measure the contribution of the memory module to the generative process. To illustrate the advantages of this approach we incorporate it into a variational autoencoder and apply the resulting model to the task of generative few-shot learning. The intuition behind this architecture is that the memory module can pick a relevant template from memory and the continuous part of the model can concentrate on modeling remaining variations. We demonstrate empirically that our model is able to identify and access the relevant memory contents even with hundreds of unseen Omniglot characters in memory


title: A stabilized trace finite element method for partial differential equations on evolving surfaces

abstract: In this paper, we study a numerical method for the solution of partial differential equations on evolving surfaces. The numerical method is built on the stabilized trace finite element method (TraceFEM) for the spatial discretization and finite differences for the time discretization. The TraceFEM uses a stationary background mesh, which can be chosen independent of time and the position of the surface. The stabilization ensures well-conditioning of the algebraic systems and defines a regular extension of the solution from the surface to its volumetric neighborhood. Having such an extension is essential for the numerical method to be well-defined. The paper proves numerical stability and optimal order error estimates for the case of simplicial background meshes and finite element spaces of order $m\ge1$. For the algebraic condition numbers of the resulting systems we prove estimates, which are independent of the position of the interface. The method allows that the surface and its evolution are given implicitly with the help of an indicator function. Results of numerical experiments for a set of 2D evolving surfaces are provided.


title: A Mathematical Aspect of Hohenberg-Kohn Theorem

abstract: The Hohenberg-Kohn theorem plays a fundamental role in density functional theory, which has become a basic tool for the study of electronic structure of matter. In this article, we study the Hohenberg-Kohn theorem for a class of external potentials based on a unique continuation principle.

categories: physics.chem-ph math-ph math.MP quant-ph


title: A Practical Quantum Algorithm for the Schur Transform

abstract: We describe an efficient quantum algorithm for the quantum Schur transform. The Schur transform is an operation on a quantum computer that maps the standard computational basis to a basis composed of irreducible representations of the unitary and symmetric groups. We simplify and extend the algorithm of Bacon, Chuang, and Harrow, and provide a new practical construction as well as sharp theoretical and practical analyses. Our algorithm decomposes the Schur transform on $n$ qubits into $O\left(n^4\log\left(\frac{n}{\epsilon}\right)\right)$ operators in the Clifford+T fault-tolerant gate set and uses exactly $2\lfloor\log_2(n)\rfloor-1$ ancillary qubits. We extend our qubit algorithm to decompose the Schur transform on $n$ qudits of dimension $d$ into $O\left(d^{1+p}n^{3d}\log^p\left(\frac{d n}{\epsilon}\right)\right)$ primitive operators from any universal gate set, for $p\approx3.97$.


title: Irreversibility of mechanical and hydrodynamic instabilities

abstract: The literature on dynamical systems has, for the most part, considered self-oscillators (i.e., systems capable of generating and maintaining a periodic motion at the expense of an external energy source with no corresponding periodicity) either as applications of the concepts of limit cycle and Hopf bifurcation in the theory of differential equations, or else as instability problems in feedback control systems. Here we outline a complementary approach, based on physical considerations of work extraction and thermodynamic irreversibility. We illustrate the power of this method with two concrete examples: the mechanical instability of rotors that spin at super-critical speeds, and the hydrodynamic Kelvin-Helmholtz instability of the interface between fluid layers with different tangential velocities. Our treatment clarifies the necessary role of frictional or viscous dissipation (and therefore of irreversibility), while revealing an underlying unity to the physics of many irreversible processes that generate mechanical work and an autonomous temporal structure (periodic, quasi-periodic, or chaotic) in the presence of an out-of-equilibrium background.


title: Discrete-Time Polar Opinion Dynamics with Susceptibility

abstract: This paper considers a discrete-time opinion dynamics model in which each individual's susceptibility to being influenced by others is dependent on her current opinion. We assume that the social network has time-varying topology and that the opinions are scalars on a continuous interval. We first propose a general opinion dynamics model based on the DeGroot model, with a general function to describe the functional dependence of each individual's susceptibility on her own opinion, and show that this general model is analogous to the Friedkin-Johnsen model, which assumes a constant susceptibility for each individual. We then consider two specific functions in which the individual's susceptibility depends on the \emph{polarity} of her opinion, and provide motivating social examples. First, we consider stubborn positives, who have reduced susceptibility if their opinions are at one end of the interval and increased susceptibility if their opinions are at the opposite end. A court jury is used as a motivating example. Second, we consider stubborn neutrals, who have reduced susceptibility when their opinions are in the middle of the spectrum, and our motivating examples are social networks discussing established social norms or institutionalized behavior. For each specific susceptibility model, we establish the initial and graph topology conditions in which consensus is reached, and develop necessary and sufficient conditions on the initial conditions for the final consensus value to be at either extreme of the opinion interval. Simulations are provided to show the effects of the susceptibility function when compared to the DeGroot model.


title: Accelerating PageRank using Partition-Centric Processing

abstract: PageRank is a fundamental link analysis algorithm that also functions as a key representative of the performance of Sparse Matrix-Vector (SpMV) multiplication. The traditional PageRank implementation generates fine granularity random memory accesses resulting in large amount of wasteful DRAM traffic and poor bandwidth utilization. In this paper, we present a novel Partition-Centric Processing Methodology (PCPM) to compute PageRank, that drastically reduces the amount of DRAM communication while achieving high sustained memory bandwidth. PCPM uses a Partition-centric abstraction coupled with the Gather-Apply-Scatter (GAS) programming model. By carefully examining how a PCPM based implementation impacts communication characteristics of the algorithm, we propose several system optimizations that improve the execution time substantially. More specifically, we develop (1) a new data layout that significantly reduces communication and random DRAM accesses, and (2) branch avoidance mechanisms to get rid of unpredictable data-dependent branches. We perform detailed analytical and experimental evaluation of our approach using 6 large graphs and demonstrate an average 2.7x speedup in execution time and 1.7x reduction in communication volume, compared to the state-of-the-art. We also show that unlike other GAS based implementations, PCPM is able to further reduce main memory traffic by taking advantage of intelligent node labeling that enhances locality. Although we use PageRank as the target application in this paper, our approach can be applied to generic SpMV computation.


title: Height, trunk and representativity of knots

abstract: In this paper, we investigate three geometrical invariants of knots, the height, the trunk and the representativity. First, we give a conterexample for the conjecture which states that the height is additive under connected sum of knots. We also define the minimal height of a knot and give a potential example which has a gap between the height and the minimal height. Next, we show that the representativity is bounded above by a half of the trunk. We also define the trunk of a tangle and show that if a knot has an essential tangle decomposition, then the representativity is bounded above by half of the trunk of either of the two tangles. Finally, we remark on the difference among Gabai's thin position, ordered thin position and minimal critical position. We also give an example of a knot which bounds an essential non-orientable spanning surface, but has arbitrarily large representativity.


title: Deep Recurrent NMF for Speech Separation by Unfolding Iterative Thresholding

abstract: In this paper, we propose a novel recurrent neural network architecture for speech separation. This architecture is constructed by unfolding the iterations of a sequential iterative soft-thresholding algorithm (ISTA) that solves the optimization problem for sparse nonnegative matrix factorization (NMF) of spectrograms. We name this network architecture deep recurrent NMF (DR-NMF). The proposed DR-NMF network has three distinct advantages. First, DR-NMF provides better interpretability than other deep architectures, since the weights correspond to NMF model parameters, even after training. This interpretability also provides principled initializations that enable faster training and convergence to better solutions compared to conventional random initialization. Second, like many deep networks, DR-NMF is an order of magnitude faster at test time than NMF, since computation of the network output only requires evaluating a few layers at each time step. Third, when a limited amount of training data is available, DR-NMF exhibits stronger generalization and separation performance compared to sparse NMF and state-of-the-art long-short term memory (LSTM) networks. When a large amount of training data is available, DR-NMF achieves lower yet competitive separation performance compared to LSTM networks.


title: Simulating the effect of boron doping in superconducting carbon

abstract: We examine the effect of boron doping in superconducting forms of amorphous carbon. By judiciously optimizing boron substitutional sites in simulated amorphous carbon we predict a superconducting transition temperature near 37 K at 14 % boron concentration. Our findings have direct implications for understanding the recently discovered high T_$c$ superconductivity in Q-carbon.


title: Measurement of Longitudinal Single-Spin Asymmetry for W$^\pm$ Production in Polarized Proton+Proton Collisions at STAR

abstract: The sea quark contribution to the nucleon spin is an important piece for a complete understanding of the nucleon spin structure. The production of W bosons in longitudinally polarized $p + p$ collisions at RHIC provides an unique probe for the sea quark polarization, through the parity-violating single-spin asymmetry, $A_L$. At the STAR experiment, W bosons through the leptonic decay channel $W \to e \nu $ can be effectively determined with the Electromagnetic Calorimeters and Time Projection Chamber at mid-rapidity. The previous STAR measurements of $A_L$ for W boson production from datasets taken in 2011 and 2012, have provided significant constraints on the helicity distribution functions of $\bar u$ and $\bar d$ quarks. In 2013 the STAR experiment collected $p+p$ data with an integrated luminosity of about 300 pb$^{-1}$ at $\sqrt s$ = 510 GeV with an average beam polarization of about $56\%$, which is more than three times larger than the total integrated luminosity of previous years. The preliminary results of W-boson $A_L$ from 2013 data sample will be presented.


title: Current-induced magnetization switching using electrically-insulating spin-torque generator

abstract: Current-induced magnetization switching through spin-orbit torques (SOTs) is the fundamental building block of spin-orbitronics. The SOTs generally arise from the spin-orbit coupling of heavy metals. However, even in a heterostructure where a metallic magnet is sandwiched by two different insulators, a nonzero current-induced SOT is expected because of the broken inversion symmetry; an electrical insulator can be a spin-torque generator. Here, we demonstrate current-induced magnetization switching using an insulator. We show that oxygen incorporation into the most widely used spintronic material, Pt, turns the heavy metal into an electrically-insulating generator of the SOTs, enabling the electrical switching of perpendicular magnetization in a ferrimagnet sandwiched by electrically-insulating oxides. We further found that the SOTs generated from the Pt oxide can be controlled electrically through voltage-driven oxygen migration. These findings open a route towards energy-efficient, voltage-programmable spin-orbit devices based on solid-state switching of heavy metal oxidation.


title: The toric Frobenius morphism and a conjecture of Orlov

abstract: We combine the Bondal-Uehara method for producing exceptional collections on toric varieties with a result of the first author and Favero to expand the set of varieties satisfying Orlov's Conjecture on derived dimension.


title: Hypergraph Theory: Applications in 5G Heterogeneous Ultra-Dense Networks

abstract: Heterogeneous ultra-dense network (HUDN) can significantly increase the spectral efficiency of cellular networks and cater for the explosive growth of data traffic in the fifth-generation (5G) communications. Due to the dense deployment of small cells (SCs), interference among neighboring cells becomes severe. As a result, the effective resource allocation and user association algorithms are essential to minimize inter-cell interference and optimize network performance. However, optimizing network resources in HUDN is extremely complicated as resource allocation and user association are coupled. Therefore, HUDN requires low-complexity but effective resource allocation schemes to address these issues. Hypergraph theory has been recognized as a useful mathematical tool to model the complex relations among multiple entities. In this article, we show how the hypergraph models can be used to effectively tackle resource allocation problems in HUDN. We also discuss several potential research issues in this field.


title: Modeling and Quantifying the Forces Driving Online Video Popularity Evolution

abstract: Video popularity is an essential reference for optimizing resource allocation and video recommendation in online video services. However, there is still no convincing model that can accurately depict a video's popularity evolution. In this paper, we propose a dynamic popularity model by modeling the video information diffusion process driven by various forms of recommendation. Through fitting the model with real traces collected from a practical system, we can quantify the strengths of the recommendation forces. Such quantification can lead to characterizing video popularity patterns, user behaviors and recommendation strategies, which is illustrated by a case study of TV episodes.


title: Fast and accurate computation of normalized Bargmann transform

abstract: The linear canonical transform (LCT) was extended to complex-valued parameters, called complex LCT, to describe the complex amplitude propagation through lossy or lossless optical systems. Bargmann transform is a special case of the complex LCT. In this paper, we normalize the Bargmann transform such that it can be bounded near infinity. We derive the relationships of the normalized Bargmann transform to Gabor transform, Hermite-Gaussian functions, gyrator transform, and 2D nonseparable LCT. Several kinds of fast and accurate computational methods of the normalized Bargmann transform and its inverse are proposed based on these relationships.


title: Engineering topological phases with a three-dimensional nodal-loop semimetal

abstract: A three-dimensional (3D) nodal-loop semimetal phase is exploited to engineer a number of intriguing phases featuring different peculiar topological surface states. In particular, by introducing various two-dimensional gap terms to a 3D tight-binding model of a nodal-loop semimetal, we obtain a rich variety of topological phases of great interest to ongoing theoretical and experimental studies, including chiral insulator, degenerate-surface-loop insulator, second-order topological insulator, as well as Weyl semimetal with tunable Fermi arc profiles. The unique concept underlying our approach is to engineer topological surface states that inherit their dispersion relations from a gap term. The results provide one rather unified principle for the creation of novel topological phases and can guide the search for new topological materials. Two-terminal transport studies are also carried out to distinguish the engineered topological phases.


title: Correlated In-Situ Low-Frequency Noise and Impedance Spectroscopy Reveal Recombination Dynamics in Organic Solar Cells using Fullerene and Non-Fullerene Acceptors

abstract: Non-fullerene acceptors based on perylenediimides (PDIs) have garnered significant interest as an alternative to fullerene acceptors in organic photovoltaics (OPVs), but their charge transport phenomena are not well understood, especially in bulk heterojunctions (BHJs). Here, we investigate charge transport and current fluctuations by performing correlated low-frequency noise and impedance spectroscopy measurements on two BHJ OPV systems, one employing a fullerene acceptor and the other employing a dimeric PDI acceptor. In the dark, these measurements reveal that PDI-based OPVs have a greater degree of recombination in comparison to fullerene-based OPVs. Furthermore, for the first time in organic solar cells, 1/f noise data are fit to the Kleinpenning model to reveal underlying current fluctuations in different transport regimes. Under illumination, 1/f noise increases by approximately four orders of magnitude for the fullerene-based OPVs and three orders of magnitude for the PDI-based OPVs. An inverse correlation is also observed between noise spectral density and power conversion efficiency. Overall, these results show that low-frequency noise spectroscopy is an effective in-situ diagnostic tool to assess charge transport in emerging photovoltaic materials, thereby providing quantitative guidance for the design of next-generation solar cell materials and technologies.


title: On the Schr\"odinger equations with time-dependent potentials growing polynomially in the spatial direction

abstract: The Cauchy problem for the Schr\"odinger equations is studied with time-dependent potentials growing polynomially in the spatial direction. First the existence and the uniqueness of solutions are shown in the weighted Sobolev spaces. In addition, we suppose that our potentials are depending on a parameter. Secondly it is shown that if potentials depend continuously and differentiably on the parameter, the solutions to the Schr\"odinger equations respectively become continuous and differentiable with respect to its parameter.


title: Maximal Moments and Uniform Modulus of Continuity for Stable Random Fields

abstract: In this work, we solve an open problem mentioned in Xiao (2010) and provide sharp bounds on the rate of growth of maximal moments for stationary symmetric stable random fields using structure theorem of finitely generated abelian groups and ergodic theory of quasi-invariant group actions. We also investigate the relationship between this rate of growth and the path regularity properties of self-similar stable random fields with stationary increments, and establish uniform modulus of continuity of such fields. In the process, a new notion of weak effective dimension is introduced for stable random fields and is connected to maximal moments and path properties. Our results establish a boundary between shorter and longer memory in relation to Holder continuity of S alpha S random fields confirming a conjecture of Samorodnitsky (2004).


title: Impacts of Surface Depletion on the Plasmonic Properties of Doped Semiconductor Nanocrystals

abstract: Degenerately doped semiconductor nanocrystals (NCs) exhibit a localized surface plasmon resonance (LSPR) in the infrared range of the electromagnetic spectrum. Unlike metals, semiconductor NCs offer tunable LSPR characteristics enabled by doping, or via electrochemical or photochemical charging. Tuning plasmonic properties through carrier density modulation suggests potential applications in smart optoelectronics, catalysis, and sensing. Here, we elucidate fundamental aspects of LSPR modulation through dynamic carrier density tuning in Sn-doped Indium Oxide NCs. Monodisperse Sn-doped Indium Oxide NCs with various doping level and sizes were synthesized and assembled in uniform films. NC films were then charged in an in situ electrochemical cell and the LSPR modulation spectra were monitored. Based on spectral shifts and intensity modulation of the LSPR, combined with optical modeling, it was found that often-neglected semiconductor properties, specifically band structure modification due to doping and surface states, strongly affect LSPR modulation. Fermi level pinning by surface defect states creates a surface depletion layer that alters the LSPR properties; it determines the extent of LSPR frequency modulation, diminishes the expected near field enhancement, and strongly reduces sensitivity of the LSPR to the surroundings.

categories: physics.optics cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph


title: $L^0$--convex compactness and its applications to random convex optimization and random variational inequalities

abstract: In 2010, Gordan \v{Z}itkovi\'{c} introduced the notion of convex compactness for a convex subset of a linear topological space and gave some important applications to both nonlinear analysis and mathematical economics in [ Gordan \v{Z}itkovi\'{c}, Convex compactness and its applications, Math. Finance Econom. 3(1) (2010) 1--12 ]. Motivated by Gordan \v{Z}itkovi\'{c}'s idea, in this paper we introduce the notion of $L^0$--convex compactness for an $L^0$--convex subset of a topological module over the topological algebra $L^0(\mathcal{F},K)$, where $L^0(\mathcal{F},K)$ is the algebra of equivalence classes of random variables from a probability space $(\Omega,\mathcal{F},P)$ to the scalar field $K$ of real numbers or complex numbers, endowed with the topology of convergence in probability. This paper continues to develop the theory of $L^0$--convex compactness by establishing various kinds of characterization theorems for $L^0$--convex subsets of a class of important topological modules--complete random normed modules, in particular, we make use of the theory of random conjugate spaces to give a characterization theorem of James type for a closed $L^0$--convex subset of a complete random normed module. As applications, we successfully generalize some basic theorems of classical convex optimization and variational inequalities from a convex function on a reflexive Banach space to an $L^0$--convex function on a random reflexive random normed module. Since the usual weak compactness method fails in the random setting of this paper and in particular, since the difficulties caused by the partial order structure of the range of an $L^0$--valued function also frequently occurs in the study of problems involved in this paper, we are forced to discover a series of new skills to meet the needs of this paper.


title: Inverse Bremsstrahlung current drive

abstract: The generation of the plasma current resulting from Bremsstrahlung absorption is considered. It is shown that the electric current is higher than the naive estimates assuming that electrons absorb only the photon momentum and using the Spitzer conductivity would suggest. The current enhancement is in part because electrons get the recoil momentum from the Coulomb field of ions during the absorption and in part because the electromagnetic power is absorbed asymmetrically within the electron velocity distribution space.


title: Learning to Prove Safety over Parameterised Concurrent Systems (Full Version)

abstract: We revisit the classic problem of proving safety over parameterised concurrent systems, i.e., an infinite family of finite-state concurrent systems that are represented by some finite (symbolic) means. An example of such an infinite family is a dining philosopher protocol with any number n of processes (n being the parameter that defines the infinite family). Regular model checking is a well-known generic framework for modelling parameterised concurrent systems, where an infinite set of configurations (resp. transitions) is represented by a regular set (resp. regular transducer). Although verifying safety properties in the regular model checking framework is undecidable in general, many sophisticated semi-algorithms have been developed in the past fifteen years that can successfully prove safety in many practical instances. In this paper, we propose a simple solution to synthesise regular inductive invariants that makes use of Angluin's classic L* algorithm (and its variants). We provide a termination guarantee when the set of configurations reachable from a given set of initial configurations is regular. We have tested L* algorithm on standard (as well as new) examples in regular model checking including the dining philosopher protocol, the dining cryptographer protocol, and several mutual exclusion protocols (e.g. Bakery, Burns, Szymanski, and German). Our experiments show that, despite the simplicity of our solution, it can perform at least as well as existing semi-algorithms.


title: Classification of rational 1-forms on the Riemann sphere up to PSL(2,C)

abstract: We study the family $\Omega^1(-1^s)$ of rational 1--forms on the Riemann sphere, having exactly $-s \leq -2$ simple poles. Three equivalent $(2s-1)$--dimensional complex atlases on $\Omega^1(-1^s)$, using coefficients, zeros--poles and residues--poles of the 1--forms, are recognized. A rational 1--form is isochronous when all their residues are purely imaginary. We prove that the subfamily $\mathcal{RI}\Omega^1(-1^s)$ of isochronous 1--forms is a $(3s-1)$--dimensional real analytic submanifold in the complex manifold $\Omega^1(-1^s)$. The complex Lie group $PSL(2,\mathbb{C})$ acts holomorphically on $\Omega^1(-1^s)$. For $s \geq 3$, the $PSL(2,\mathbb{C})$--action is proper on $\Omega^1(-1^s)$ and $\mathcal{RI}\Omega^1(-1^s)$. Therefore, the quotients $\Omega^1(-1^s)/PSL(2,\mathbb{C})$ and $\mathcal{RI}\Omega^1(-1^s)/PSL(2,\mathbb{C})$ admit a stratification by orbit types. Using an explicit set of $PSL(2,\mathbb{C})$--invariant functions, we give realizations for the quotients $\Omega^1(-1^s)/PSL(2,\mathbb{C})$ and $\mathcal{RI}\Omega^1(-1^s)/PSL(2,\mathbb{C})$.


title: The k-tacnode process

abstract: The tacnode process is a universal behavior arising in nonintersecting particle systems and tiling problems. For Dyson Brownian bridges, the tacnode process describes the grazing collision of two packets of walkers. We consider such a Dyson sea on the unit circle with drift. For any integer k, we show that an appropriate double scaling of the drift and return time leads to a generalization of the tacnode process in which k particles are expected to wrap around the circle. We derive winding number probabilities and an expression for the correlation kernel in terms of functions related to the generalized Hastings-McLeod solutions to the inhomogeneous Painleve-II equation. The method of proof is asymptotic analysis of discrete orthogonal polynomials with a complex weight.


title: The $B \to \pi K$ Puzzle Revisited

abstract: For a number of years, there has been a certain inconsistency among the measurements of the branching ratios and CP asymmetries of the four $B \to \pi K$ decays ($B^+ \to \pi^+ K^0$, $B^+ \to \pi^0 K^+$, $B^0 \to \pi^- K^+$, $B^0 \to \pi^0 K^0$). In this paper, we re-examine this $B \to \pi K$ puzzle. We find that the key unknown parameter is $|C'/T'|$, the ratio of color-suppressed and color-allowed tree amplitudes. If this ratio is large, $|C'/T'| = 0.5$, the SM can explain the data. But if it is small, $|C'/T'| = 0.2$, the SM cannot explain the $B \to \pi K$ puzzle -- new physics (NP) is needed. The two types of NP that can contribute to $B \to \pi K$ at tree level are $Z'$ bosons and diquarks. $Z'$ models can explain the puzzle if the $Z'$ couples to right-handed $u{\bar u}$ and/or $d{\bar d}$, with $g_R^{dd} \ne g_R^{uu}$. Interestingly, half of the many $Z'$ models proposed to explain the present anomalies in $b \to s \mu^+ \mu^-$ decays have the required $Z'$ couplings to $u{\bar u}$ and/or $d{\bar d}$. Such models could potentially explain both the $b \to s \mu^+ \mu^-$ anomalies and the $B \to \pi K$ puzzle. The addition of a color sextet diquark that couples to $ud$ can also explain the puzzle.


title: Fractional iterated Ornstein-Uhlenbeck Processes

abstract: In this work we present a Gaussian process that arise from the iteration of p fractional Ornstein-Uhlenbeck processes generated by the same fractional Brownian motion. This iteration results, when the values of lambdas are pairwise differents, in a particular linear combination of those processes. Although for $H>1/2$ each term of the linear combination is a long memory processes, we prove that it results in a short memory processes. We include applications to real data that show improvement in predictive performance compared with different ARMA models.


title: Radiative transfer dynamo effect

abstract: Magnetic fields in rotating and radiating astrophysical plasma can be produced due to a radiative interaction between plasma layers moving relative to each other. The efficiency of current drive, and with it the associated dynamo effect, is considered in a number of limits. It is shown here, however, that predictions for these generated magnetic fields can be significantly higher when kinetic effects, previously neglected, are taken into account.


title: Hyperbolic carbon nanoforest for phase matching of ordinary and backward electromagnetic waves: second harmonic generation

abstract: We show that deliberately engineered dispersive metamaterial slab can enable the co-existence and phase matching of contra-propagating ordinary fundamental and extraordinary backward second harmonic surface electromagnetic modes. Energy flux and phase velocity are contra-directed in the backward waves which is the phenomenon that gives rise to unique nonlinear optical propagation processes. We show that frequencies, phase, and group velocities, as well as nanowaveguide losses inherent to the electromagnetic modes supported by such metamaterial, can be tailored to maximize conversion of frequencies and to reverse propagation direction of the generated wave. Such a possibility, which is of paramount importance for nonlinear photonics, is proved with a numerical model of the hyperbolic metamaterial made of carbon nanotubes standing on the metal surface. Extraordinary properties of the backward-wave second harmonic generation in the reflection direction and of the corresponding frequency doubling metareflector in the THz are investigated with a focus on the pulsed regime.


title: Redefining climate index with sea surface temperature of key area can better express atmospheric physical activity

abstract: Since the 1920s, when the three oscillations ENSO, NAO and NPO were found and defined, the theory of ocean-atmosphere coupling is particularly important in explaining the phenomena of atmospheric physics, then defined more than ten climate indexes, such as SOI, AO, AAO and Nino monitoring area and so on, the indexes usually adopt standardized sea level pressure (SLP) difference or standardized sea surface temperature anomaly (SSTA) difference definitions. The survey is carried out with seventy stations runoff of major rivers in the world shows that the correlation coefficient between these indexes and runoff is very low, most of which is between 0.1-0.2, the impact of the index on the hydrologic climate cannot be truly reflected, which has led to doubts about the clear physical mechanism and even questions of teleological authenticity. This is due to the fact that standardization is a nonlinear process, which changes the structure of the original data, causing the correlation coefficient to be low. Using sea surface temperature difference (SSTD) to define these climate indexes, the correlation coefficient is significantly improved, the Hankou station of Yangtze river, Darband station of Ind, Hope station of Fraser, Bangui station of Oubangui generally exceeds 0.8 strong degree, which prove that these relevant elements have a high resonance effect, there must be a common driving physical quantity. We find that the sun-earth distance has a very strong resonance effect, it is the earth's rotational orbit that determines their common activity intensity and cycle characteristics. The sun-earth distance can accurately calculate the tens and hundreds years of SST and SSTD, if using the SST or SSTD as the climate indicators, which express physical process of various ocean-atmosphere coupling, we can thoroughly solve long-term ultra-long-term forecast of international difficulty problems.


title: Two-dimensional Dirac particles in a P\"oschl-Teller waveguide

abstract: We obtain exact solutions to the two-dimensional (2D) Dirac equation for the one-dimensional P\"oschl-Teller potential which contains an asymmetry term. The eigenfunctions are expressed in terms of Heun confluent functions, while the eigenvalues are determined via the solutions of a simple transcendental equation. For the symmetric case, the eigenfunctions of the supercritical states are expressed as spheroidal wave functions, and approximate analytical expressions are obtained for the corresponding eigenvalues. A universal condition for any square integrable symmetric potential is obtained for the minimum strength of the potential required to hold a bound state of zero energy. Applications for smooth electron waveguides in 2D Dirac-Weyl systems are discussed.

categories: cond-mat.mes-hall math-ph math.MP quant-ph


title: Lens elliptic gamma function solution of the Yang-Baxter equation at roots of unity

abstract: We study the root of unity limit of the lens elliptic gamma function solution of the star-triangle relation, for an integrable model with continuous and discrete spin variables. This limit involves taking an elliptic nome to a primitive $rN$-th root of unity, where $r$ is an existing integer parameter of the lens elliptic gamma function, and $N$ is an additional integer parameter. This is a singular limit of the star-triangle relation, and at subleading order of an asymptotic expansion, another star-triangle relation is obtained for a model with discrete spin variables in $\mathbb{Z}_{rN}$. Some special choices of solutions of equation of motion are shown to result in well-known discrete spin solutions of the star-triangle relation. The saddle point equations themselves are identified with three-leg forms of "3D-consistent" classical discrete integrable equations, known as $Q4$ and $Q3_{(\delta=0)}$. We also comment on the implications for supersymmetric gauge theories, and in particular comment on a close parallel with the works of Nekrasov and Shatashvili.

categories: math-ph cond-mat.stat-mech hep-th math.MP nlin.SI


title: Learning RBM with a DC programming Approach

abstract: By exploiting the property that the RBM log-likelihood function is the difference of convex functions, we formulate a stochastic variant of the difference of convex functions (DC) programming to minimize the negative log-likelihood. Interestingly, the traditional contrastive divergence algorithm is a special case of the above formulation and the hyperparameters of the two algorithms can be chosen such that the amount of computation per mini-batch is identical. We show that for a given computational budget the proposed algorithm almost always reaches a higher log-likelihood more rapidly, compared to the standard contrastive divergence algorithm. Further, we modify this algorithm to use the centered gradients and show that it is more efficient and effective compared to the standard centered gradient algorithm on benchmark datasets.


title: Feature Engineering for Predictive Modeling using Reinforcement Learning

abstract: Feature engineering is a crucial step in the process of predictive modeling. It involves the transformation of given feature space, typically using mathematical functions, with the objective of reducing the modeling error for a given target. However, there is no well-defined basis for performing effective feature engineering. It involves domain knowledge, intuition, and most of all, a lengthy process of trial and error. The human attention involved in overseeing this process significantly influences the cost of model generation. We present a new framework to automate feature engineering. It is based on performance driven exploration of a transformation graph, which systematically and compactly enumerates the space of given options. A highly efficient exploration strategy is derived through reinforcement learning on past examples.


title: Effect of Prandtl number on heat transport enhancement in Rayleigh-B\'enard convection under geometrical confinement

abstract: We study, using direct numerical simulations, the effect of geometrical confinement on heat transport and flow structure in Rayleigh-B\'enard convection in fluids with different Prandtl numbers. Our simulations span over two decades of Prandtl number $Pr$, $0.1 \leq Pr \leq 40$, with the Rayleigh number $Ra$ fixed at $10^8$. The width-to-height aspect ratio $\Gamma$ spans between $0.025$ and $0.25$ while the length-to-height aspect ratio is fixed at one. We first find that for $Pr \geq 0.5$, geometrical confinement can lead to a significant enhancement in heat transport as characterized by the Nusselt number $Nu$. For those cases, $Nu$ is maximal at a certain $\Gamma = \Gamma_{opt}$. It is found that $\Gamma_{opt}$ exhibits a power-law relation with $Pr$ as $\Gamma_{opt}=0.11Pr^{-0.06}$, and the maximal relative enhancement generally increases with $Pr$ over the explored parameter range. As opposed to the situation of $Pr \geq 0.5$, confinement-induced enhancement in $Nu$ is not realized for smaller values of $Pr$, such as $0.1$ and $0.2$. The $Pr$ dependence of the heat transport enhancement can be understood in its relation to the coverage area of the thermal plumes over the thermal boundary layer (BL) where larger coverage is observed for larger $Pr$ due to a smaller thermal diffusivity. We further show that $\Gamma_{opt}$ is closely related to the crossing of thermal and momentum BLs, and find that $Nu$ declines sharply when the thickness ratio of the thermal and momentum BLs exceeds a certain value of about one. In addition, through examining the temporally averaged flow fields and 2D mode decomposition, it is found that for smaller $Pr$ the large-scale circulation is robust against the geometrical confinement of the convection cell.


title: Observing Compact Stars with AstroSat

abstract: This article presents a brief description of India's AstroSat mission which is a powerful space based observatory for compact star research. An account is given of observational constraints and spectral and timing capabilities as realised post-launch. Some preliminary results of observations of the Crab pulsar and an X-ray binary system GX~301-2 are presented to illustrate some of the capabilities of the mission.


title: Large Vocabulary Automatic Chord Estimation Using Deep Neural Nets: Design Framework, System Variations and Limitations

abstract: In this paper, we propose a new system design framework for large vocabulary automatic chord estimation. Our approach is based on an integration of traditional sequence segmentation processes and deep learning chord classification techniques. We systematically explore the design space of the proposed framework for a range of parameters, namely deep neural nets, network configurations, input feature representations, segment tiling schemes, and training data sizes. Experimental results show that among the three proposed deep neural nets and a baseline model, the recurrent neural network based system has the best average chord quality accuracy that significantly outperforms the other considered models. Furthermore, our bias-variance analysis has identified a glass ceiling as a potential hindrance to future improvements of large vocabulary automatic chord estimation systems.


title: Gravitation as a source of decoherence

abstract: It is believed that classical behavior emerges in a quantum system due to decoherence. It has also been proposed that gravity can be a source of this decoherence. We examine this in detail by studying a number of quantum systems, including ultra-relativistic and non-relativistic particles, at low and high temperatures in an expanding Universe, and show that this proposal is valid for a large class of quantum systems.