title: Laser-induced ultrafast demagnetization time and spin moment in ferromagnets: First-principles calculation

abstract: When a laser pulse excites a ferromagnet, its spin undergoes a dramatic change. The initial demagnetization process is very fast. Experimentally, it is found that the demagnetization time is related to the spin moment in the sample. In this study, we employ the first-principles method to directly simulate such a process. We use the fixed spin moment method to change the spin moment in ferromagnetic nickel, and then we employ the Liouville equation to couple the laser pulse to the system. We find that in general the dependence of demagnetization time on the spin moment is nonlinear: It decreases with the spin moment up to a point, after which an increase with the spin moment is observed, followed by a second decrease. To understand this, we employ an extended Heisenberg model, which includes both the exchange interaction and spin-orbit coupling. The model directly links the demagnetization rate to the spin moment itself and demonstrates analytically that the spin relaxes more slowly with a small spin moment. A future experimental test of our predictions is needed.


title: Adaptive Controls of FWER and FDR Under Block Dependence

abstract: Often in multiple testing, the hypotheses appear in non-overlapping blocks with the associated $p$-values exhibiting dependence within but not between blocks. We consider adapting the Benjamini-Hochberg method for controlling the false discovery rate (FDR) and the Bonferroni method for controlling the familywise error rate (FWER) to such dependence structure without losing their ultimate controls over the FDR and FWER, respectively, in a non-asymptotic setting. We present variants of conventional adaptive Benjamini-Hochberg and Bonferroni methods with proofs of their respective controls over the FDR and FWER. Numerical evidence is presented to show that these new adaptive methods can capture the present dependence structure more effectively than the corresponding conventional adaptive methods. This paper offers a solution to the open problem of constructing adaptive FDR and FWER controlling methods under dependence in a non-asymptotic setting and providing real improvements over the corresponding non-adaptive ones.


title: Ultrafast reduction of exchange splitting in ferromagnetic nickel

abstract: A decade ago Rhie \et (Phys. Rev. Lett. {\bf 90}, 247201 (2003)) reported that when ferromagnetic nickel is subject to an intense ultrashort laser pulse, its exchange splitting is reduced quickly. But to simulate such reduction remains a big challenge. The popular rigid band approximation (RBA), where both the band structure and the exchange splitting are held fixed before and after laser excitation, is unsuitable for this purpose, while the time-dependent density functional theory could be time-consuming. To overcome these difficulties, we propose a time-dependent Liouville and density functional theory (TDLDFT) that integrates the time-dependent Liouville equation into the density functional theory. As a result, the excited charge density is reiterated back into the Kohn-Sham equation, and the band structure is allowed to change dynamically. Even with the ground-state density functional, a larger demagnetization than RBA is found; after we expand Ortenzi's spin scaling method into an excited-state (laser) density functional, we find that the exchange splitting is indeed strongly reduced, as seen in the experiment. Both the majority and minority bands are shifted toward the Fermi level, but the majority shifts a lot more. The ultrafast reduction in exchange splitting occurs concomitantly with demagnetization. While our current theory is still unable to yield the same percentage loss in the spin moment as observed in the experiment, it predicts a correct trend that agrees with the experiments. With a better functional, we believe that our results can be further improved.


title: Analysis of the nonleptonic charmonium modes $B_s^0 \to J/\psi f_{2}^{\prime}(1525)$ and $B_s^0 \to J/\psi K^+K^-$

abstract: In this work, we present an analysis of the nonleptonic charmonium modes $B_s^0 \to J / \psi f_2^\prime(1525)$ and $B_s^0 \to J/\psi K^+K^-$. Within the framework of the factorization approach and using the perturbative QCD for the evaluation of the relevant form factors, we find a branching fraction for the two-body channel of BR$(B_s^{0} \to J/\psi f_2^{\prime}(1525)) = (1.6 {}^{+0.9}_{-0.7})\times 10^{-4}$ which is in agreement with the experimental values reported by the LHCb and Belle Collaborations. The associated polarization fractions to this vector-tensor mode are also presented. On the other hand, non-resonant and resonant contributions to the three-body decay $B_s^0 \to J/\psi K^+K^-$ are carefully investigated. The dominant contributions of the resonances $ \phi(1020)$ and $f_2^{\prime}(1525)$ are properly taken into account. A detailed analysis of the $K^+ K^-$ invariant mass distributions and Dalitz plot are also performed. The overall result BR$(B_s^{0} \to J/\psi K^+ K^-) = (9.3^{+1.3}_{-1.1})\times 10^{-4}$ is also in satisfactory agreement with the experimental information reported by LHCb and Belle.


title: Using Neural Networks to Compute Approximate and Guaranteed Feasible Hamilton-Jacobi-Bellman PDE Solutions

abstract: To sidestep the curse of dimensionality when computing solutions to Hamilton-Jacobi-Bellman partial differential equations (HJB PDE), we propose an algorithm that leverages a neural network to approximate the value function. We show that our final approximation of the value function generates near optimal controls which are guaranteed to successfully drive the system to a target state. Our framework is not dependent on state space discretization, leading to a significant reduction in computation time and space complexity in comparison with dynamic programming-based approaches. Using this grid-free approach also enables us to plan over longer time horizons with relatively little additional computation overhead. Unlike many previous neural network HJB PDE approximating formulations, our approximation is strictly conservative and hence any trajectories we generate will be strictly feasible. For demonstration, we specialize our new general framework to the Dubins car model and discuss how the framework can be applied to other models with higher-dimensional state spaces.


title: Scalable Compression of a Weighted Graph

abstract: Graph is a useful data structure to model various real life aspects like email communications, co-authorship among researchers, interactions among chemical compounds, and so on. Supporting such real life interactions produce a knowledge rich massive repository of data. However, efficiently understanding underlying trends and patterns is hard due to large size of the graph. Therefore, this paper presents a scalable compression solution to compute summary of a weighted graph. All the aforementioned interactions from various domains are represented as edge weights in a graph. Therefore, creating a summary graph while considering this vital aspect is necessary to learn insights of different communication patterns. By experimenting the proposed method on two real world and publically available datasets against a state of the art technique, we obtain order of magnitude performance gain and better summarization accuracy.


title: Discovery of a fourth arc in Abell 2626 at 610 MHz with the GMRT: Spectral properties and possibilities for the origin

abstract: We report the discovery of a fourth eastern arc (Arc E) towards the cool-core cluster Abell 2626 using 610 MHz Giant Metrewave Radio Telescope observations. Three arcs towards north, west and south were known from earlier works at 1400 MHz and proposed to have originated in precessing radio jets of the central active galactic nucleus. The 610 - 1400 MHz integrated spectral indices of the arcs are in the range 3.2 - 3.6 and the spectral index map shows uniform distribution along the lengths of the arcs. If associated with A2626, the arcs have linear extents in the range 79 - 152 kpc. The detection of Arc E favours the scenario in which a pair of bipolar precessing jets were active and halted to produce the arc system. Based on the morphological symmetry and spectral similarity, we indicate a possible role of gravitational lensing. Further high resolution low frequency observations and measurements of the mass of the system are needed to disentangle the mystery of this source.


title: HoTCoffeeh: Hanbury Brown-Twiss correlation functions and radii from event-by-event hydrodynamics

abstract: $\textbf{HoTCoffeeh}$ ($\textbf{H}$anbury Br$\textbf{o}$wn-$\textbf{T}$wiss $\textbf{Co}$rrelation $\textbf{f}$unctions and radii $\textbf{f}$rom $\textbf{e}$vent-by-$\textbf{e}$vent $\textbf{h}$ydrodynamics) is a new computational tool which determines Hanbury Brown-Twiss (HBT) charged pion ($\pi^+$) correlation functions and radii for event-by-event (EBE) hydrodynamics with fluctuating initial conditions in terms of Cooper-Frye integrals, including resonance decay contributions. In this paper, we review the basic formalism for computing the HBT correlation functions and radii with resonance decay contributions included, and discuss our implementation of this formalism in the form of HoTCoffeeh. This tool may be easily integrated with other numerical packages for the purpose of simulating the evolution of heavy-ion collisions and thereby extracting predictions for heavy-ion observables.


title: Autonomous learning by simple dynamical systems with a discrete-time formulation

abstract: We present a discrete-time formulation for the autonomous learning conjecture. The main feature of this formulation is the possibility to apply the autonomous learning scheme to systems in which the errors with respect to target functions are not well-defined for all times. This restriction for the evaluation of functionality is a typical feature in systems that need a finite time interval to process a unit piece of information. We illustrate its application on an artificial neural network with feed-forward architecture for classification and a phase oscillator system with synchronization properties. The main characteristics of the discrete-time formulation are shown by constructing these systems with predefined functions.


title: Mutual information of excited states and relative entropy of two disjoint subsystems in CFT

abstract: In this paper, we first study mutual information of excited states in the small subsystem size limit in generic conformal field theory. We then discuss relative entropy of two disjoint subsystems in the same limit.


title: Massless black holes and charged wormholes in string theory

abstract: We discuss the zero mass pointlike solutions and charged Einstein-Rosen bridges (wormholes) that arise from the dyonic black hole solution of the Einstein-Maxwell-dilaton theory. These massless black holes exist individually in spacetime, different from the known massless solutions, which come in pairs with opposite signs for their masses. In order to construct a massless object, we choose the integration constants of the solution to have specific values. The massless solutions present some problems: in one case the dilaton field is complex (or the gauge field has negative kinetic energy), and in the other case the solution has negative entropy and temperature or it is naked singularity in the extremal limit. For the first case, the observables computed are real quantities. This massless solution also allow the bridge construction, and we obtain an analytical and static charged wormhole solution, which satisfies the null energy condition.


title: Typical length scales in conducting disorderless networks

abstract: We take advantage of a recently established equivalence, between the intermittent dynamics of a deterministic nonlinear map and the scattering matrix properties of a disorderless double Cayley tree lattice of connectivity $K$, to obtain general electronic transport expressions and expand our knowledge of the scattering properties at the mobility edge. From this we provide a physical interpretation of the generalized localization length.


title: Non-crossing chords of a polygon with forbidden positions

abstract: In this paper, we systematically study non-crossing chords of simple polygons in the plane. We first introduce the reduced Euler characteristic of a family of line-segments, and subsequently investigate the structure of the diagonals and epigonals of a polygon. Interestingly enough, the reduced Euler characteristic of a subfamily of diagonals and epigonals characterizes the geometric convexity of polygons. In particular, an alternative and complete answer is given for a problem proposed by G. C. Shephard. Meanwhile, we extend such research to non-crossing diagonals and epigonals with forbidden positions in some appropriate sense. We prove that the reduced Euler characteristic of diagonals with forbidden positions only depends on the information involving convex partitions by those forbidden diagonals, and it determines the shapes of polygons in a surprising way. Incidentally, some kinds of generalized Catalan's numbers naturally arise.


title: Weighted iteration complexity of the sPADMM on the KKT residuals for convex composite optimization

abstract: In this paper we establish an $\mathcal{O}({1}/{k})$ weighted iteration complexity on the KKT residuals yielded by the sPADMM (semi-proximal alternating direction method of multiplier) for the convex composite optimization problem. This result, which is derived with the help of a novel generalized HPE (hybrid proximal extra-gradient) iteration formula, first fills the gap on the ergodic iteration complexity of the classic ADMM with a large step-size and its many proximal variants.


title: Tactics and Tallies: Inferring Voter Preferences in the 2016 U.S. Presidential Primaries Using Sparse Learning

abstract: In this paper, we propose a web-centered framework to infer voter preferences for the 2016 U.S. presidential primaries. Using Twitter data collected from Sept. 2015 to March 2016, we first uncover the tweeting tactics of the candidates and then exploit the variations in the number of 'likes' to infer voters' preference. With sparse learning, we are able to reveal neutral topics as well as positive and negative ones. Methodologically, we are able to achieve a higher predictive power with sparse learning. Substantively, we show that for Hillary Clinton the (only) positive issue area is women's rights. We demonstrate that Hillary Clinton's tactic of linking herself to President Obama resonates well with her supporters but the same is not true for Bernie Sanders. In addition, we show that Donald Trump is a major topic for all the other candidates, and that the women's rights issue is equally emphasized in Sanders' campaign as in Clinton's.


title: On $\mathbb{Z}_{2}\mathbb{Z}_{2}[u]$-$(1+u)$-additive constacyclic

abstract: In this paper, we study $\mathbb{Z}_{2}\mathbb{Z}_{2}[u]$-$(1+u)$-additive constacyclic code of arbitrary length. Firstly, we study the algebraic structure of this family of codes and a set of generator polynomials for this family as a $(\mathbb{Z}_{2}+u\mathbb{Z}_{2})[x]$-submodule of the ring $R_{\alpha,\beta}$. Secondly, we give the minimal generating sets of this family codes, and we determine the relationship of generators between the $\mathbb{Z}_{2}\mathbb{Z}_{2}[u]$-$(1+u)$-additive constacyclic codes and its dual and give the parameters in terms of the degrees of the generator polynomials of the code. Lastly, we also study $\mathbb{Z}_{2}\mathbb{Z}_{2}[u]$-$(1+u)$-additive constacyclic code in terms of the Gray images.


title: Tycho-Gaia Astrometric Solution Parallaxes and Proper Motions for Five Galactic Globular Clusters

abstract: We present a pilot study of Galactic globular cluster (GC) proper motion (PM) determinations using Gaia data. We search for GC stars in the Tycho-Gaia Astrometric Solution (TGAS) catalogue from Gaia Data Release 1 (DR1), and identify five members of NGC104 (47 Tucanae), one member of NGC5272 (M3), five members of NGC6121 (M4), seven members of NGC6397, and two members of NGC6656 (M22). By taking a weighted average of member stars, fully accounting for the correlations between parameters, we estimate the parallax (and, hence, distance) and PM of the GCs. This provides a homogeneous PM study of multiple GCs based on an astrometric catalog with small and well-controlled systematic errors, and yields random PM errors similar to existing measurements. Detailed comparison to the available Hubble Space Telescope (HST) measurements generally shows excellent agreement, validating the astrometric quality of both TGAS and HST. By contrast, comparison to ground-based measurements shows that some of those must have systematic errors exceeding the random errors. Our parallax estimates have uncertainties an order of magnitude larger than previous studies, but nevertheless imply distances consistent with previous estimates. By combining our PM measurements with literature positions, distances, and radial velocities, we measure Galactocentric space motions for the clusters and find that these also agree well with previous analyses. Our analysis provides a framework for determining more accurate distances and PMs of Galactic GCs using future Gaia data releases. This will provide crucial constraints on the near end of the cosmic distance ladder and provide accurate GC orbital histories.


title: Faint Object Detection in Multi-Epoch Observations via Catalog Data Fusion

abstract: Observational astronomy in the time-domain era faces several new challenges. One of them is the efficient use of observations obtained at multiple epochs. The work presented here addresses faint object detection with multi-epoch data, and describes an incremental strategy for separating real objects from artifacts in ongoing surveys, in situations where the single-epoch data are summaries of the full image data, such as single-epoch catalogs of flux and direction estimates for candidate sources. The basic idea is to produce low-threshold single-epoch catalogs, and use a probabilistic approach to accumulate catalog information across epochs; this is in contrast to more conventional strategies based on co-added or stacked image data across all epochs. We adopt a Bayesian approach, addressing object detection by calculating the marginal likelihoods for hypotheses asserting there is no object, or one object, in a small image patch containing at most one cataloged source at each epoch. The object-present hypothesis interprets the sources in a patch at different epochs as arising from a genuine object; the no-object (noise) hypothesis interprets candidate sources as spurious, arising from noise peaks. We study the detection probability for constant-flux objects in a simplified Gaussian noise setting, comparing results based on single exposures and stacked exposures to results based on a series of single-epoch catalog summaries. Computing the detection probability based on catalog data amounts to generalized cross-matching: it is the product of a factor accounting for matching of the estimated fluxes of candidate sources, and a factor accounting for matching of their estimated directions. We find that probabilistic fusion of multi-epoch catalog information can detect sources with only modest sacrifice in sensitivity and selectivity compared to stacking.


title: 2-Selmer groups of hyperelliptic curves with two marked points

abstract: We consider the family of hyperelliptic curves over $\Q$ of fixed genus along with a marked rational Weierstrass point and a marked rational non-Weierstrass point. When these curves are ordered by height, we prove that the average Mordell-Weil rank of their Jacobians is bounded above by 5/2. We prove this by showing that the average rank of the 2-Selmer groups is bounded above by 6. We also prove that the average size of the $\phi$-Selmer groups of a family of isogenies associated to this family is exactly 2.


title: The Response of an Elastic Splitter Plate Attached to a Cylinder to Laminar Pulsatile Flow

abstract: The flow-induced deformation of a thin, elastic splitter plate attached to the rear of a circular cylinder and subjected to laminar pulsatile inflow is investigated. The cylinder and elastic splitter plate are contained within a narrow channel and the Reynolds number is mostly restricted to Re = 100, primarily covering the two-dimensional flow regime. An in-house fluid-structure interaction code is employed for simulations, which couples a sharp-interface immersed boundary method for the fluid dynamics with a finite-element method to treat the structural dynamics. The structural solver is implicitly (two-way) coupled with the flow solver using a partitioned approach. This implicit coupling ensures numerical stability at low structure-fluid density ratios. A power spectrum analysis of the time-varying plate displacement shows that the plate oscillates at more than a single frequency for pulsatile inflow, compared to a single frequency observed for steady inflow. The multiple frequencies obtained for the former case can be explained by beating between the applied and plate oscillatory signals. The plate attains a self-sustained time-periodic oscillation with a plateau amplitude in the case of steady flow, while the superimposition of pulsatile inflow with induced plate oscillation affects the plateau amplitude. Lock-in of the plate oscillation with the pulsatile inflow occurs at a forcing frequency that is twice of the plate natural frequency in a particular mode and this mode depends on the plate length. The plate displacement as well as pressure drag increases at the lock-in condition. The percentage change in the maximum plate displacement, and skin-friction and pressure drag coefficients on the plate, due to pulsatile inflow is quantified. The non-linear dynamics of the plate and its coupling with the pulsatile flow are briefly discussed.


title: Spectral and pseudospectral functions of various dimensions for symmetric systems

abstract: The main object of the paper is a symmetric system $J y'-B(t)y=\l\D(t) y$ defined on an interval $\cI=[a,b) $ with the regular endpoint $a$. Let $\f(\cd,\l)$ be a matrix solution of this system of an arbitrary dimension and let $(Vf)(s)=\int\limits_\cI \f^*(t,s)\D(t)f(t)\,dt$ be the Fourier transform of the function $f(\cd)\in L_\D^2(\cI)$. We define a pseudospectral function of the system as a matrix-valued distribution function $\s(\cd)$ of the dimension $n_\s$ such that $V$ is a partial isometry from $L_\D^2(\cI)$ to $L^2(\s;\bC^{n_\s})$ with the minimally possible kernel. Moreover, we find the minimally possible value of $n_\s$ and parameterize all spectral and pseudospectral functions of every possible dimensions $n_\s$ by means of a Nevanlinna boundary parameter. The obtained results develop the results by Arov and Dym; A.~Sakhnovich, L.~Sakhnovich and Roitberg; Langer and Textorius.


title: Unified Music Theories for General Equal-Temperament Systems

abstract: Why are white and black piano keys in an octave arranged as they are today? This article examines the relations between abstract algebra and key signature, scales, degrees, and keyboard configurations in general equal-temperament systems. Without confining the study to the twelve-tone equal-temperament (12-TET) system, we propose a set of basic axioms based on musical observations. The axioms may lead to scales that are reasonable both mathematically and musically in any equal-temperament system. We reexamine the mathematical understandings and interpretations of ideas in classical music theory, such as the circle of fifths, enharmonic equivalent, degrees such as the dominant and the subdominant, and the leading tone, and endow them with meaning outside of the 12-TET system. In the process of deriving scales, we create various kinds of sequences to describe facts in music theory, and we name these sequences systematically and unambiguously with the aim to facilitate future research.


title: 5G Multimedia Massive MIMO Communications Systems

abstract: In the Fifth generation (5G) wireless communication systems, a majority of the traffic demands is contributed by various multimedia applications. To support the future 5G multimedia communication systems, the massive multiple-input multiple-output (MIMO) technique is recognized as a key enabler due to its high spectral efficiency. The massive antennas and radio frequency (RF) chains not only improve the implementation cost of 5G wireless communication systems but also result in an intense mutual coupling effect among antennas because of the limited space for deploying antennas. To reduce the cost, an optimal equivalent precoding matrix with the minimum number of RF chains is proposed for 5G multimedia massive MIMO communication systems considering the mutual coupling effect. Moreover, an upper bound of the effective capacity is derived for 5G multimedia massive MIMO communication systems. Two antenna receive diversity gain models are built and analyzed. The impacts of the antenna spacing, the number of antennas, the quality of service (QoS) statistical exponent, and the number of independent incident directions on the effective capacity of 5G multimedia massive MIMO communication systems are analyzed. Comparing with the conventional zero-forcing precoding matrix, simulation results demonstrate that the proposed optimal equivalent precoding matrix can achieve a higher achievable rate for 5G multimedia massive MIMO communication systems.


title: A Note on Random Walks with Absorbing barriers and Sequential Monte Carlo Methods

abstract: In this article we consider importance sampling (IS) and sequential Monte Carlo (SMC) methods in the context of 1-dimensional random walks with absorbing barriers. In particular, we develop a very precise variance analysis for several IS and SMC procedures. We take advantage of some explicit spectral formulae available for these models to derive sharp and explicit estimates; this provides stability properties of the associated normalized Feynman-Kac semigroups. Our analysis allows one to compare the variance of SMC and IS techniques for these models. The work in this article, is one of the few to consider an in-depth analysis of an SMC method for a particular model-type as well as variance comparison of SMC algorithms.


title: Noise reduction combining microphone and piezoelectric device

abstract: It is often required to extract the sound of an objective instrument played in concert with other instruments. Microphone array is one of the effective ways to enhance a sound from a specific direction. However it is not effective in an echoic room such as concert hall. The pickup microphone attached on the specific musical instrument is often employed to obtain the sound exclusively from other instrumental sounds. The obtained timbre differ from the one we hear at the usual listening position. The purpose of this paper is to propose a new method of sound separation that utilizes the piezoelectric device attached on the body of the instrument. The signal from the attached device has a different spectrum from the sound heard by the audience but has the same frequency components as the instrumental sound. Our idea is to use the device signal as a modifier of the sound focusing filter applied to the microphone sound at the listening position. The proposed method firstly estimates the frequency components of the signal from the piezoelectric device. The frequency characteristics for filtering the microphone sound are changed so that it pass the estimated frequency components. Thus we can extractthe target sound without distortion. The proposed method is a sort of dynamic sparseness approach. It was found that SNR is improved by 8.7dB through the experiments.


title: Extended period domains, algebraic groups, and higher Albanese manifolds

abstract: For a linear algebraic group G over the field of rational numbers, we consider the period domains D classifying G-mixed Hodge structures, and construct the extended period domains as toroidal partial compactifications. We give an interpretation of higher Albanese manifolds by Hain and Zucker by using the above D for some G, and extend them via the toroidal partial compactifications.


title: Design and performance of SiPM-based readout of PbF2 crystals for high-rate, precision timing applications

abstract: We have developed a custom amplifier board coupled to a large-format 16-channel Hamamatsu silicon photomultiplier device for use as the light sensor for the electromagnetic calorimeters in the Muon g-2 experiment at Fermilab. The calorimeter absorber is an array of lead-fluoride crystals, which produces short-duration Cherenkov light. The detector sits in the high magnetic field of the muon storage ring. The SiPMs selected, and their accompanying custom electronics, must preserve the short pulse shape, have high quantum efficiency, be non-magnetic, exhibit gain stability under varying rate conditions, and cover a fairly large fraction of the crystal exit surface area. We describe an optimized design that employs the new-generation of thru-silicon via devices. The performance is documented in a series of bench and beam tests.


title: On the algorithmic complexity of adjacent vertex closed distinguishing colorings number of graphs

abstract: An assignment of numbers to the vertices of graph G is closed distinguishing if for any two adjacent vertices v and u the sum of labels of the vertices in the closed neighborhood of the vertex v differs from the sum of labels of the vertices in the closed neighborhood of the vertex u unless they have the same closed neighborhood (i.e. N[u]=N[v]). The closed distinguishing number of G, denoted by dis[G], is the smallest integer k such that there is a closed distinguishing labeling for G using integers from the set[k].Also, for each vertex $v \in V(G)$, let L(v) denote a list of natural numbers available at v. A list closed distinguishing labeling is a closed distinguishing labeling f such that $f(v)\in L(v)$ for each $v \in V(G)$.A graph G is said to be closed distinguishing k-choosable if every k-list assignment of natural numbers to the vertices of G permits a list closed distinguishing labeling of G. The closed distinguishing choice number of G, $dis_{\ell}[G]$, is the minimum number k such that G is closed distinguishing k-choosable. We show that for each integer t there is a bipartite graph G such that $dis[G] > t$.It was shown that for every graph G with $\Delta\geq 2$, $dis[G]\leq dis_{\ell}[G]\leq \Delta^2-\Delta+1$ and there are infinitely values of $\Delta$ for which G might be chosen so that $dis[G] =\Delta^2-\Delta+1$. We show that the difference between $dis[G]$ and $dis_{\ell}[G]$ can be arbitrary large and for every positive integer t there is a graph G such that $dis_{\ell}[G]-dis[G]\geq t$. We improve the current upper bound and give some number of upper bounds for the closed distinguishing choice number by using the Combinatorial Nullstellensatz. We show that it is $\mathbf{NP}$-complete to decide for a given planar subcubic graph G, whether dis[G]=2. Also, we prove that for every $k\geq 3$, it is {\bf NP}-complete to decide whether $dis[G]=k$ for a given graph G


title: Search for GeV and X-ray flares associated with the IceCube track-like neutrinos

abstract: Dozens of high-energy neutrinos have been detected by the IceCube neutrino telescope, but no clear association with any classes of astrophysical sources has been identified so far. Recently, Kadler et al. (2016) report that a PeV cascade-like neutrino event occurs in positional and temporal coincidence with a giant gamma-ray flare of the blazar PKS B1424-418. Since IceCube track-like events have much better angular resolution, we here search for possible short-term gamma-ray flares that are associated with the IceCube track-like events with Fermi Large Area Telescope (LAT) observations. Among them, three track-like neutrino events occur within the field of view of Fermi-LAT at the time of the detection, so search for the {\em prompt} gamma-ray emission associated with neutrinos are possible. Assuming a point source origin and a single power law spectrum for the possible gamma-ray sources associated with neutrinos, a likelihood analysis of 0.2-100 GeV photons observed by Fermi-LAT on the timescales of $\sim 12$ hours and one year are performed, and for the three special neutrinos, the analysis are also performed on the timescales of thousand of seconds before and after the neutrino detection. No significant GeV excesses over the background are found and the upper limit fluxes at 95\% confidence level are obtained for different timescales. We also search for possible hard X-ray transient sources associated with the IceCube track-like neutrino events, but the search also yields null results. We discuss the implication of the non-detection of gamma-ray flares for the constraints on the neutrino source density.


title: Massive Machine Type Communication with Data Aggregation and Resource Scheduling

abstract: To enable massive machine type communication (mMTC), data aggregation is a promising approach to reduce the congestion caused by a massive number of machine type devices (MTDs). In this work, we consider a two-phase cellular-based mMTC network where MTDs transmit to aggregators (i.e., aggregation phase) and the aggregated data is then relayed to base stations (i.e., relaying phase). Due to the limited resources, the aggregators not only aggregate data, but also schedule resources among MTDs. We consider two scheduling schemes: random resource scheduling (RRS) and channel-aware resource scheduling (CRS). By leveraging the stochastic geometry, we present a tractable analytical framework to investigate the signal-to-interference ratio (SIR) for each phase, thereby computing the MTD success probability, the average number of successful MTDs and probability of successful channel utilization, which are the key metrics characterizing the overall mMTC performance. Our numerical results show that, although the CRS outperforms the RRS in terms of SIR at the aggregation phase, the simpler RRS has almost the same performance as the CRS for most cases with regards to the overall mMTC performance. Furthermore, the provision of more resources at the aggregation phase is not always beneficial to the mMTC performance.


title: Searching for Dark Matter Annihilation in Recently Discovered Milky Way Satellites with Fermi-LAT

abstract: We search for excess gamma-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using 6 years of data from the Fermi Large Area Telescope (LAT). Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted gamma-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are dark-matter-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each ~$2 \sigma$ local) for gamma-ray emission in excess of the background. However, the ensemble of derived gamma-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance $<1 \sigma$). Our analysis has increased sensitivity compared to the analysis of 15 confirmed dSphs by Ackermann et al. 2015. The observed constraints on the dark matter annihilation cross section are statistically consistent with the background expectation, improving by a factor of ~2 for large dark matter masses ($m_{{\rm DM},b \bar b} \gtrsim 1$ TeV and $m_{{\rm DM},\tau^{+}\tau^{-}} \gtrsim 70$ GeV) and weakening by a factor of ~1.5 at lower masses relative to previously observed limits.


title: A quantum phase transition detected through one dimensional ballistic conductance

abstract: A quantum phase transition is an unequivocal signature of strongly correlated many-body physics. Signatures of such phenomena are yet to be observed in ballistic transport through quantum wires. Recent developments in quantum wires have made it possible to enhance the interaction between the electrons. Here we show that hitherto unexplained anticrossing between conduction energy sub-bands, observed in such experiments, can be explained through a simple yet effective discretised model which undergoes a second order quantum phase transition within the Ising universality class. Accordingly, we observe how the charge distribution, transverse to the direction of the wire will vary across the phase transition. We show that data coming from three different samples with differing electron densities and gate voltages show a remarkable universal scaling behavior, determined by the relevant critical exponent, which is only possible near a quantum phase transition.


title: SoK: Applying Machine Learning in Security - A Survey

abstract: The idea of applying machine learning(ML) to solve problems in security domains is almost 3 decades old. As information and communications grow more ubiquitous and more data become available, many security risks arise as well as appetite to manage and mitigate such risks. Consequently, research on applying and designing ML algorithms and systems for security has grown fast, ranging from intrusion detection systems(IDS) and malware classification to security policy management(SPM) and information leak checking. In this paper, we systematically study the methods, algorithms, and system designs in academic publications from 2008-2015 that applied ML in security domains. 98 percent of the surveyed papers appeared in the 6 highest-ranked academic security conferences and 1 conference known for pioneering ML applications in security. We examine the generalized system designs, underlying assumptions, measurements, and use cases in active research. Our examinations lead to 1) a taxonomy on ML paradigms and security domains for future exploration and exploitation, and 2) an agenda detailing open and upcoming challenges. Based on our survey, we also suggest a point of view that treats security as a game theory problem instead of a batch-trained ML problem.


title: Universal Hinge Patterns for Folding Strips Efficiently into Any Grid Polyhedron

abstract: We present two universal hinge patterns that enable a strip of material to fold into any connected surface made up of unit squares on the 3D cube grid--for example, the surface of any polycube. The folding is efficient: for target surfaces topologically equivalent to a sphere, the strip needs to have only twice the target surface area, and the folding stacks at most two layers of material anywhere. These geometric results offer a new way to build programmable matter that is substantially more efficient than what is possible with a square $N \times N$ sheet of material, which can fold into all polycubes only of surface area $O(N)$ and may stack $\Theta(N^2)$ layers at one point. We also show how our strip foldings can be executed by a rigid motion without collisions (albeit assuming zero thickness), which is not possible in general with 2D sheet folding. To achieve these results, we develop new approximation algorithms for milling the surface of a grid polyhedron, which simultaneously give a 2-approximation in tour length and an 8/3-approximation in the number of turns. Both length and turns consume area when folding a strip, so we build on past approximation algorithms for these two objectives from 2D milling.


title: Emergence of small numbers in complex systems and the origin of the electroweak scale

abstract: In sufficiently complex models with many parameters that are unknown or undetermined from first principles, a small coupling or mass can naturally arise even if it is not protected by a symmetry or a result of some dynamics. For the naturalness criterion, we advocate specifying all the model parameters with one significant figure. This automatically avoids outcomes of a given observable that would require special highly-tuned choices for parameters. The obtained range of outcomes is an attribute of a given model resulting from its specific characteristics and complexity. Using these criteria in the minimal supersymmetric model, we demonstrate that the electroweak scale up to 3 orders of magnitude below superpartner masses naturally occurs.


title: Sparse Hierarchical Solvers with Guaranteed Convergence

abstract: Solving sparse linear systems from discretized PDEs is challenging. Direct solvers have in many cases quadratic complexity (depending on geometry), while iterative solvers require problem dependent preconditioners to be robust and efficient. Approximate factorization preconditioners, such as incomplete LU factorization, provide cheap approximations to the system matrix. However, even a highly accurate preconditioner may have deteriorating performance when the condition number of the system matrix increases. By increasing the accuracy on low-frequency errors, we propose a novel hierarchical solver with improved robustness with respect to the condition number of the linear system. This solver retains the linear computational cost and memory footprint of the original algorithm.


title: On estimating the force-freeness based on observed magnetograms

abstract: It is a common practice in the solar physics community to test whether or not measured photospheric or chromospheric vector magnetograms are force-free, using the Maxwell stress as a measure. Some previous studies have suggested that magnetic fields of active regions in the solar chromosphere are close to be force-free whereas there is no consistency among previous studies on whether magnetic fields of active regions in the solar photosphere are force-free or not. Here we use three kinds of representative magnetic fields (analytical force-free solutions, modeled solar-like force-free fields and observed non-force-free fields) to discuss on how the measurement issues such as limited field of view, instrument sensitivity and measurement error could affect the estimation of force-freeness based on observed magnetograms. Unlike previous studies that focus on discussing the effect of limited field of view or instrument sensitivity, our calculation shows that just measurement error alone can significantly influence the results of force-freeness estimate, due to the fact that measurement errors in horizontal magnetic fields are usually ten times larger than that of the vertical fields. This property of measurement errors, interacting with the particular form of force-freeness estimate formula, would result in wrong judgments of the force-freeness: a truly force-free field may be mistakenly estimated as being non-force-free and a true non-force-free field may be estimated as being force-free. Our analysis calls for caution when interpreting the force-freeness estimates based on measured magnetograms, and also suggests that the true photospheric magnetic field may be further away from being force-free than they currently appear to be.


title: Exploring shared genetic bases and causal relationships of schizophrenia and bipolar disorder with 28 cardiovascular and metabolic traits

abstract: Cardiovascular diseases (CVD) represent a major health issue in patients with schizophrneia (SCZ) and bipolar disorder (BD), but the exact nature of cardiometabolic (CM) abnormalities involved and the underlying mechanisms remain unclear. Using polygenic risk scores (PRS) and LD score regression, we investigated the shared genetic bases of SCZ and BD with a panel of 28 cardiometabolic traits. We performed Mendelian randomization (MR) to elucidate casual relationships between the two groups of disorders. The analysis was based on large-scale meta-analyses of genome-wide association studies (GWAS). We also identified the potential shared genetic variants by a statistical approach based on local true discovery rates, and inferred the pathways involved. We found polygenic associations of SCZ with glucose metabolism abnormalities, adverse adipokine profiles, increased wait-hip ratio and raised visceral adiposity. However, BMI showed inverse genetic correlation and polygenic link with SCZ. On the other hand, we observed polygenic associations with an overall favorable CM profile in BD. MR analysis showed that SCZ may be causally linked to raised triglyceride and that lower fasting glucose may be linked to BD; otherwise MR did not reveal other significant causal relationships in general. We also identified numerous SNPs and pathways shared between SCZ/BD with cardiometabolic traits, some of which are related to inflammation or the immune system. In conclusion, SCZ patients may be genetically associated with several CM abnormalities independent of medication side-effects, and proper surveillance and management of CV risk factors may be required from the onset of the disease. On the other hand, CM abnormalities in BD are more likely to be secondary.


title: Transitional steady states of exchange dynamics between finite quantum systems

abstract: We examine energy and particle exchange between finite-sized quantum systems and find a new form of nonequilibrium states. The exchange rate undergoes stepwise evolution in time, and its magnitude and sign dramatically change according to system size differences. The origin lies in interference effects contributed by multiply scattered waves at system boundaries. Although such characteristics are utterly different from those of true steady state for infinite systems, Onsager's reciprocal relation remains universally valid.


title: Mahalanobis Distance for Class Averaging of Cryo-EM Images

abstract: Single particle reconstruction (SPR) from cryo-electron microscopy (EM) is a technique in which the 3D structure of a molecule needs to be determined from its contrast transfer function (CTF) affected, noisy 2D projection images taken at unknown viewing directions. One of the main challenges in cryo-EM is the typically low signal to noise ratio (SNR) of the acquired images. 2D classification of images, followed by class averaging, improves the SNR of the resulting averages, and is used for selecting particles from micrographs and for inspecting the particle images. We introduce a new affinity measure, akin to the Mahalanobis distance, to compare cryo-EM images belonging to different defocus groups. The new similarity measure is employed to detect similar images, thereby leading to an improved algorithm for class averaging. We evaluate the performance of the proposed class averaging procedure on synthetic datasets, obtaining state of the art classification.


title: Full counting statistics of phonon-assisted Andreev tunneling through a quantum dot coupled to normal and superconducting leads

abstract: We present a theoretical investigation for the full counting statistics of the Andreev tunneling through a quantum dot (QD) embedded between superconducting (SC) and normal leads in the presence of a strong on-site electron-phonon interaction using nonequilibrium Green function method. For this purpose, we generalize the dressed tunneling approximation (DTA) recently developed in dealing with inelastic tunneling in a normal QD system to the Andreev transport issue, which takes account of vibrational effect in evaluation of electronic tunneling self energy in comparison with other simple approaches and meanwhile allows us to derive an explicit analytical formula for the cumulant generating function at the subgap region. We then analyze the interplay of polaronic and SC proximity effects on the Andreev reflection spectrum, current-voltage characteristics, and current fluctuations of the hybrid system. Our main findings include: (1) no phonon side peaks in the linear Andreev conductance; (2) a negative differential conductance stemming from the suppressed Andreev reflection spectrum; (3) a novel inelastic resonant peak in the differential conductance due to phonon assisted Andreev reflection; (4) enhancement or suppression of shot noise for the symmetric or asymmetric tunnel-coupling system respectively.


title: Complete Affine Connection in the Causal Boundary: Static, Spherically Symmetric Spacetimes

abstract: The boundary at $\Cal I^+$, future null infinity, for a standard static, spherically symmetric spactime is examined for possible linear connections. Two independent methods are employed, one for treating $\Cal I^+$ as the future causal boundary, and one for treating it as a conformal boundary (the latter is subsumed in the former, which is of greater generality). Both methods provide the same result: a constellation of various possible connections, depending on an arbitrary choice of a certain function, a sort of gauge freedom in obtaining a natural connection on $\Cal I^+$; choosing that function to be constant (for instance) results in a complete connection. Treating $\Cal I^+$ as part of the future causal boundary, the method is to impute affine connections on null hypersurfaces going out to $\Cal I^+$, in terms of a transverse vector field on each null hypersurface (there is much gauge freedom on choice of the transverse vector fields). Treating $\Cal I^+$ as part of a conformal boundary, the method is to make a choice of conformal factor that makes the boundary totally geodesic in the enveloping manifold (there is much gauge freedom in choice of that conformal factor). Similar examination is made of other boundaries, such as timelike infinity and timelike and spacelike singularities. These are much simpler, as they admit a unique connection from a similar limiting process (i.e., no gauge freedom); and that connection is complete.


title: Fair representation by independent sets

abstract: For a hypergraph $H$ let $\beta(H)$ denote the minimal number of edges from $H$ covering $V(H)$. An edge $S$ of $H$ is said to represent {\em fairly} (resp. {\em almost fairly}) a partition $(V_1,V_2, \ldots, V_m)$ of $V(H)$ if $|S\cap V_i|\ge \lfloor\frac{|V_i|}{\beta(H)}\rfloor$ (resp. $|S\cap V_i|\ge \lfloor\frac{|V_i|}{\beta(H)}\rfloor-1$) for all $i \le m$. In matroids any partition of $V(H)$ can be represented fairly by some independent set. We look for classes of hypergraphs $H$ in which any partition of $V(H)$ can be represented almost fairly by some edge. We show that this is true when $H$ is the set of independent sets in a path, and conjecture that it is true when $H$ is the set of matchings in $K_{n,n}$. We prove that partitions of $E(K_{n,n})$ into three sets can be represented almost fairly. The methods of proofs are topological.


title: Multiple Service Chain Placement and Routing in a Network-enabled Cloud

abstract: Network Function Virtualization (NFV) aims to abstract the functionality of traditional proprietary hardware into software as Virtual Network Functions (VNFs), which can run on commercial off the shelf (COTS) servers. Besides reducing dependency on proprietary support, NFV helps network operators to deploy multiple services in a agile fashion. Service deployment involves placement and in sequence routing through VNFs comprising a Service Chain (SC). Our study is the first to focus on the computationally complex problem of multiple VNF SC placement and routing while considering VNF service chaining explicitly. We propose a novel column generation model for placing multiple VNF SCs and routing, which reduces the computational complexity of the problem significantly. Our aim here is to determine the ideal NFV Infrastructure (NFVI) for minimizing network resource consumption. Our results indicate that a Network enabled Cloud (NeC) results in lower networkresource consumption than a centralized NFVI (e.g., Data Center) while avoiding infeasibility with a distributed NFVI.


title: Observation of resonant scattering between ultracold heteronuclear Feshbach molecules

abstract: We report the observation of a dimer-dimer inelastic collision resonance for ultracold Feshbach molecules made of bosonic sodium and rubidium atoms. This resonance, which we attribute to the crossing of the dimer-dimer threshold with a heteronuclear tetramer state, manifests itself as a pronounced inelastic loss peak of dimers when the interspecies scattering length between the constituent atoms is tuned. Near this resonance, a strong modification of the temperature dependence of the dimer-dimer scattering is observed. Our result provides insight into the heteronuclear four-body system consisting of heavy and light bosons and offers the possibility of investigating ultracold molecules with tunable interactions.


title: Low Data Drug Discovery with One-shot Learning

abstract: Recent advances in machine learning have made significant contributions to drug discovery. Deep neural networks in particular have been demonstrated to provide significant boosts in predictive power when inferring the properties and activities of small-molecule compounds. However, the applicability of these techniques has been limited by the requirement for large amounts of training data. In this work, we demonstrate how one-shot learning can be used to significantly lower the amounts of data required to make meaningful predictions in drug discovery applications. We introduce a new architecture, the residual LSTM embedding, that, when combined with graph convolutional neural networks, significantly improves the ability to learn meaningful distance metrics over small-molecules. We open source all models introduced in this work as part of DeepChem, an open-source framework for deep-learning in drug discovery.


title: Finite-size corrections and scaling for the dimer model on the checkerboard lattice

abstract: Lattice models are useful for understanding behaviors of interacting complex many-body systems. The lattice dimer model has been proposed to study the adsorption of diatomic molecules on a substrate. Here we analyze the partition function of the dimer model on an $2 M \times 2 N$ checkerboard lattice wrapped on a torus and derive the exact asymptotic expansion of the logarithm of the partition function. We find that the internal energy at the critical point is equal to zero. We also derive the exact finite-size corrections for the free energy, the internal energy, and the specific heat. Using the exact partition function and finite-size corrections for the dimer model on finite checkerboard lattice we obtain finite-size scaling functions for the free energy, the internal energy, and the specific heat of the dimer model. We investigate the properties of the specific heat near the critical point and find that specific-heat pseudocritical point coincides with the critical point of the thermodynamic limit, which means that the specific-heat shift exponent $\lambda$ is equal to $\infty$. We have also considered the limit $N \to \infty$ for which we obtain the expansion of the free energy for the dimer model on the infinitely long cylinder. From a finite-size analysis we have found that two conformal field theories with the central charges $c = 1$ for the height function description and $c = -2$ for the construction using a mapping of spanning trees can be used to describe the dimer model on the checkerboard lattice.


title: Twist-3 Distribution Amplitudes of Pion in the Light-Front Quark Model

abstract: We analyzed two twist-3 distribution amplitudes of pion, i.e. pseudoscalar $\phi^P_{3;\pi}(x)$ and pseudotensor $\phi^\sigma_{3;\pi}(x)$, within the LFQM. Our LFQM descriptions both for twist-3 $\phi^P_{3;\pi}$ and $\phi^\sigma_{3;\pi}$ not only satisfy the fundamental constraint required from the isospin symmetry, but also reproduce exactly the asymptotic forms anticipated from QCD's conformal limit.


title: Fast Adaptation of Activity Sensing Policies in Mobile Devices

abstract: With the proliferation of sensors, such as accelerometers, in mobile devices, activity and motion tracking has become a viable technology to understand and create an engaging user experience. This paper proposes a fast adaptation and learning scheme of activity tracking policies when user statistics are unknown a priori, varying with time, and inconsistent for different users. In our stochastic optimization, user activities are required to be synchronized with a backend under a cellular data limit to avoid overcharges from cellular operators. The mobile device is charged intermittently using wireless or wired charging for receiving the required energy for transmission and sensing operations. Firstly, we propose an activity tracking policy by formulating a stochastic optimization as a constrained Markov decision process (CMDP). Secondly, we prove that the optimal policy of the CMDP has a threshold structure using a Lagrangian relaxation approach and the submodularity concept. We accordingly present a fast Q-learning algorithm by considering the policy structure to improve the convergence speed over that of conventional Q-learning. Finally, simulation examples are presented to support the theoretical findings of this paper.


title: High orbital angular momentum harmonic generation

abstract: We identify and explore a high orbital angular momentum (OAM) harmonics generation and amplification mechanism that manipulates the OAM independently of any other laser property, by preserving the initial laser wavelength, through stimulated Raman backscattering in a plasma. The high OAM harmonics spectra can extend at least up to the limiting value imposed by the paraxial approximation. We show with theory and particle-in-cell simulations that the orders of the OAM harmonics can be tuned according to a selection rule that depends on the initial OAM of the interacting waves. We illustrate the high OAM harmonics generation in a plasma using several examples including the generation of prime OAM harmonics. The process can also be realised in any nonlinear optical Kerr media supporting three-wave interactions.


title: Top-k Spatial-keyword Publish/Subscribe Over Sliding Window

abstract: With the prevalence of social media and GPS-enabled devices, a massive amount of geo-textual data has been generated in a stream fashion, leading to a variety of applications such as location-based recommendation and information dissemination. In this paper, we investigate a novel real-time top-k monitoring problem over sliding window of streaming data; that is, we continuously maintain the top-k most relevant geo-textual messages (e.g., geo-tagged tweets) for a large number of spatial-keyword subscriptions (e.g., registered users interested in local events) simultaneously. To provide the most recent information under controllable memory cost, sliding window model is employed on the streaming geo-textual data. To the best of our knowledge, this is the first work to study top-k spatial-keyword publish/subscribe over sliding window. A novel centralized system, called Skype (Topk Spatial-keyword Publish/Subscribe), is proposed in this paper. In Skype, to continuously maintain top-k results for massive subscriptions, we devise a novel indexing structure upon subscriptions such that each incoming message can be immediately delivered on its arrival. To reduce the expensive top-k re-evaluation cost triggered by message expiration, we develop a novel cost-based k-skyband technique to reduce the number of re-evaluations in a cost-effective way. Extensive experiments verify the great efficiency and effectiveness of our proposed techniques. Furthermore, to support better scalability and higher throughput, we propose a distributed version of Skype, namely, DSkype, on top of Storm, which is a popular distributed stream processing system. With the help of fine-tuned subscription/message distribution mechanisms, DSkype can achieve orders of magnitude speed-up than its centralized version.


title: Generalized thermalization for integrable system under quantum quench

abstract: We investigate equilibration and generalized thermalization of the quantum Harmonic chain under local quantum quench. The quench action we consider is connecting two disjoint harmonic chains of different sizes and the system jumps between two integrable settings. We verify the validity of the Generalized Gibbs Ensemble description for this infinite dimensional Hilbert space system and also identify equilibration between the subsystems as in classical systems. Using Bogoliubov transformations, we show that the eigenstates of the system prior to the quench evolve towards the Gibbs Generalized Ensemble description. Eigenstates that are more delocalized (in the sense of inverse participation ratio) prior to the quench, tend to equilibrate more rapidly. Further, through the phase space properties of a Generalized Gibbs Ensemble and the strength of stimulated emission, we identify the necessary criterion on the initial states for such relaxation at late times and also find out the states which would potentially not be described by the Gibbs Generalized Ensemble description.

categories: quant-ph cond-mat.stat-mech gr-qc hep-th


title: Probabilistic Energy Management for Building Climate Comfort in Smart Thermal Grids with Seasonal Storage Systems

abstract: This paper presents an energy management framework for building climate comfort (BCC) systems interconnected in a grid via aquifer thermal energy storage (ATES) systems in the presence of two types of uncertainty (private and common). ATES can be used either as a heat source (hot well) or sink (cold well) depending on the season. We consider the uncertain thermal energy demand of individual buildings as a private uncertainty source and the uncertain common resource pool (ATES) between neighbors as a common uncertainty source. We develop a large-scale stochastic hybrid dynamical model to predict the thermal energy imbalance in a network of interconnected BCC systems together with mutual interactions between their local ATES. We formulate a finite-horizon mixed-integer quadratic optimization problem with multiple chance constraints at each sampling time, which is in general a non-convex problem and difficult to solve. We then provide a computationally tractable framework by extending the so-called robust randomized approach and offering a less conservative solution for a problem with multiple chance constraints. A simulation study is provided to compare completely decoupled, centralized and move-blocking centralized solutions. We also present a numerical study using a geohydrological simulation environment (MODFLOW) to illustrate the advantages of our proposed framework.


title: Ultrafast dynamics of non-equilibrium electrons and strain generation under femtosecond laser irradiation of Nickel

abstract: We present a theoretical study of the ultrafast electron dynamics in transition metals using ultrashort pulsed laser beams. The significant influence of the contribution of the dynamics of produced nonthermal electrons to electron thermalisation and electron-phonon interaction is thoroughly investigated within a range of values of the pulse duration (i.e. from 10fs to 2.3ps). The theoretical model elaborates firstly on possible transient changes of optical parameters during irradiation due to the presence of the non-equilibrium electrons and the modification of the electron distribution function in a material characterized by a variable electron density of states around the Fermi energy. The model correlates the role of nonthermal electrons, relaxation processes and induced stress-strain fields. Simulations are presented by choosing Nickel (Ni) as a test material to compute electron-phonon relaxation time due to its large electron-phonon coupling constant that expects to reveal more conclusively the influence of the non-equilibrium electrons. We demonstrate the consideration of the above factors leads to significant changes compared to the results the traditional Two Temperature Model (TTM) provides. The proposed model predicts a substantially (~33%) smaller damage threshold and a large increase of the stress (~20%, at early times) which firstly underlines the role of the nonthermal electron interactions and secondly enhances its importance with respect to the precise determination of laser specifications in material micromachining techniques.


title: SU(N) affine Toda solitons and breathers from transparent Dirac potentials

abstract: Transparent scalar and pseudoscalar potentials in the one-dimensional Dirac equation play an important role as self-consistent mean fields in 1+1 dimensional four-fermion theories (Gross-Neveu, Nambu-Jona Lasinio models) and quasi-one dimensional superconductors (Bogoliubov-De Gennes equation). Here, we show that they also serve as seed to generate a large class of classical multi-soliton and multi-breather solutions of su(N) affine Toda field theories, including the Lax representation and the corresponding vector. This generalizes previous findings about the relationship between real kinks in the Gross-Neveu model and classical solitons of the sinh-Gordon equation to complex twisted kinks.


title: Vacuum Measurements and Quantum State Reconstruction of Phonons

abstract: A quantum state is fully characterized by its density matrix or equivalently by its quasiprobabilities in phase space. A scheme to identify the quasiprobabilities of a quantum state is an important tool in the recent development of quantum technologies. Based on our highly efficient vacuum measurement scheme, we measure the quasiprobability $Q$-function of the vibrational motion for a \Yb ion {\it resonantly} interacting with its internal energy states. This interaction model is known as the Jaynes-Cummings model which is one of the fundamental models in quantum electrodynamics. We apply the capability of the vacuum measurement to study the Jaynes-Cummings dynamics, where the Gaussian peak of the initial coherent state is known to bifurcate and rotate around the origin of phase space. They merge at the so-called revival time at the other side of phase space. The measured $Q$-function agrees with the theoretical prediction. Moreover, we reconstruct the Wigner function by deconvoluting the $Q$-function and observe the quantum interference in the Wigner function at half of the revival time, where the vibrational state becomes nearly disentangled from the internal energy states and forms a superposition of two composite states. The scheme can be applied to other physical setups including cavity or circuit-QED and optomechanical systems.


title: Benchmark of the bootstrap current simulation in helical plasmas

abstract: The importance of the parallel momentum conservation on the bootstrap current evaluation in nonaxisymmetric systems is demonstrated by the benchmarks among the local drift-kinetic equation solvers, i.e., the Zero-Orbit-width(ZOW), DKES, and PENTA codes. The ZOW model is extended to include the ion parallel mean flow effect on the electron-ion parallel friction. Compared to the DKES model in which only the pitch-angle-scattering term is included in the collision operator, the PENTA model employs the Sugama-Nishimura method to correct the momentum balance. The ZOW and PENTA models agree each other well on the calculations of the bootstrap current. The DKES results without the parallel momentum conservation deviates significantly from those from the ZOW and PENTA models. This work verifies the reliability of the bootstrap current calculation with the ZOW and PENTA models for the helical plasmas.


title: Test beam studies of the light yield, time and coordinate resolutions of scintillator strips with WLS fibers and SiPM readout

abstract: Prototype scintilator+WLS strips with SiPM readout for large muon detection systems were tested in the muon beam of the Fermilab Test Beam Facility. Light yield of up to 137 photoelectrons per muon per strip has been observed, as well as time resolution of 330 ps and position resolution along the strip of 5.4 cm.


title: Boxy/Peanut/X-shape bulges: steep inner rotation curve leads to barlens face-on morphology

abstract: We use stellar dynamical bulge/disk/halo simulations to study whether barlenses (lens-like structures embedded in the narrow bar component) are just the face-on counterparts of Boxy/Peanut/X-shapes (B/P/X) seen in edge-on bars, or if some additional physical parameter affects that morphology. A range of bulge-to-disk mass and size ratios are explored: our nominal parameters ($B/D=0.08$, $r_{\rm eff}/h_r=0.07$, disk comprising 2/3 of total force at $2.2h_r$) correspond to typical MW mass galaxies. In all models a bar with pronounced B/P/X forms in a few Gyrs, visible in edge-on view. However, the pure barlens morphology forms only in models with sufficiently steep inner rotation curves, $dV_{cir}/dr\gtrsim5V_{max}/h_r$, achieved when including a small classical bulge with $B/D\gtrsim0.02$ and $r_{\rm eff}/h_r\lesssim0.1$. For shallower slopes the central structure still resembles a barlens, but shows a clear X-signature even in low inclinations. Similar result holds for bulgeless simulations, where the central slope is modified by changing the halo concentration. The predicted sensitivity on inner rotation curve is consistent with the slopes estimated from gravitational potentials calculated from the 3.6$\mu$m images, for the observed barlens and X-shape galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G). For inclinations $<60^\circ$ the galaxies with barlenses have on average twice steeper inner rotation curves than galaxies with X-shapes: the limiting slope is $\sim250$km/s/kpc. Among barred galaxies, those with barlenses have both the strongest bars and the largest relative excess of inner surface density, both in barlens region ($\lesssim0.5h_r$) and near the center ($\lesssim0.1h_r$); this provides evidence for bar-driven secular evolution in galaxies.


title: Length Matters: Clustering System Log Messages using Length of Words

abstract: The analysis techniques of system log messages (syslog messages) have a long history from when the syslog mechanism was invented. Typically, the analysis consists of two parts, one is a message template generation, and the other is finding something interesting using the messages classified by the inferred templates. It is important to generate better templates to achieve better, precise, or convincible analysis results. In this paper, we propose a classification methodology using the length of words of each message. Our method is suitable for online template generation because it does not require two-pass analysis to generate template messages, that is an important factor considering increasing amount of log messages produced by a large number of system components such as cloud infrastructure.


title: Ultimate tensorization: compressing convolutional and FC layers alike

abstract: Convolutional neural networks excel in image recognition tasks, but this comes at the cost of high computational and memory complexity. To tackle this problem, [1] developed a tensor factorization framework to compress fully-connected layers. In this paper, we focus on compressing convolutional layers. We show that while the direct application of the tensor framework [1] to the 4-dimensional kernel of convolution does compress the layer, we can do better. We reshape the convolutional kernel into a tensor of higher order and factorize it. We combine the proposed approach with the previous work to compress both convolutional and fully-connected layers of a network and achieve 80x network compression rate with 1.1% accuracy drop on the CIFAR-10 dataset.


title: CMS software and computing for LHC Run 2

abstract: The CMS offline software and computing system has successfully met the challenge of LHC Run 2. In this presentation, we will discuss how the entire system was improved in anticipation of increased trigger output rate, increased rate of pileup interactions and the evolution of computing technology. The primary goals behind these changes was to increase the flexibility of computing facilities where ever possible, as to increase our operational efficiency, and to decrease the computing resources needed to accomplish the primary offline computing workflows. These changes have resulted in a new approach to distributed computing in CMS for Run 2 and for the future as the LHC luminosity should continue to increase. We will discuss changes and plans to our data federation, which was one of the key changes towards a more flexible computing model for Run 2. Our software framework and algorithms also underwent significant changes. We will summarize the our experience with a new multi-threaded framework as deployed on our prompt reconstruction farm for 2015 and across the CMS WLCG Tier-1 facilities. We will discuss our experience with a analysis data format which is ten times smaller than our primary Run 1 format. This "miniAOD" format has proven to be easier to analyze while be extremely flexible for analysts. Finally, we describe improvements to our workflow management system that have resulted in increased automation and reliability for all facets of CMS production and user analysis operations.


title: Application of Tomographic Reconstruction Techniques for Density Analysis of Green Bodies

abstract: Progress in the manufacturing of ceramics, but also of sintered metals, strongly relies on the evaluation of the density distribution in green bodies. This evaluation is crucial from many points of view, including the calibration of constitutive models for \textit{in-silico} simulation of densification processes. To this end, X-ray tomography and other techniques are possible but can be unmanageable for some institutions. Therefore, a destructive method is introduced in the present article to measure the density field of a green body sample using a CNC mill, an analytical balance, and analysis techniques from the field of computational tomography. A virtual experiment is presented where the method is used to reconstruct a simulated green body density field and is found to satisfactorily correspond to the original solution. The green body density field of a truncated cylinder made of alumina powder is evaluated using this method and the reconstructed field is presented.


title: Real Time Video Analysis using Smart Phone Camera for Stroboscopic Image

abstract: Motion capturing and there by segmentation of the motion of any moving object from a sequence of continuous images or a video is not an exceptional task in computer vision area. Smart-phone camera application is an added integration for the development of such tasks and it also provides for a smooth testing. A new approach has been proposed for segmenting out the foreground moving object from the background and then masking the sequential motion with the static background which is commonly known as stroboscopic image. In this paper the whole process of the stroboscopic image construction technique has been clearly described along with some necessary constraints which is due to the traditional problem of estimating and modeling dynamic background changes. The background subtraction technique has been properly estimated here and number of sequential motion have also been calculated with the correlation between the motion of the object and its time of occurrence. This can be a very effective application that can replace the traditional stroboscopic system using high end SLR cameras, tripod stand, shutter speed control and position etc.


title: Learning to Play Guess Who? and Inventing a Grounded Language as a Consequence

abstract: Acquiring your first language is an incredible feat and not easily duplicated. Learning to communicate using nothing but a few pictureless books, a corpus, would likely be impossible even for humans. Nevertheless, this is the dominating approach in most natural language processing today. As an alternative, we propose the use of situated interactions between agents as a driving force for communication, and the framework of Deep Recurrent Q-Networks for evolving a shared language grounded in the provided environment. We task the agents with interactive image search in the form of the game Guess Who?. The images from the game provide a non trivial environment for the agents to discuss and a natural grounding for the concepts they decide to encode in their communication. Our experiments show that the agents learn not only to encode physical concepts in their words, i.e. grounding, but also that the agents learn to hold a multi-step dialogue remembering the state of the dialogue from step to step.


title: Statistical regionalization for estimation of extreme river discharges

abstract: Regionalization methods have long been used to estimate high return levels of river discharges at ungauged locations on a river network. In these methods, the recorded discharge measurements of a group of similar, gauged, stations is used to estimate high quantiles at the target catchment that has no observations. This group is called the region of influence and its similarity to the ungauged location is measured in terms of physical and meteorological catchment attributes. We develop a statistical method for estimation of high return levels based on regionalizing the parameters of a generalized extreme value distribution. The region of influence is chosen in an optimal way, ensuring similarity and in-group homogeneity. Our method is applied to discharge data from the Rhine basin in Switzerland, and its performance at ungauged locations is compared to that of classical regionalization methods. For gauged locations we show how our approach improves the estimation uncertainty for long return periods by combining local measurements with those from the region of influence.


title: Faster Kernel Ridge Regression Using Sketching and Preconditioning

abstract: Kernel Ridge Regression (KRR) is a simple yet powerful technique for non-parametric regression whose computation amounts to solving a linear system. This system is usually dense and highly ill-conditioned. In addition, the dimensions of the matrix are the same as the number of data points, so direct methods are unrealistic for large-scale datasets. In this paper, we propose a preconditioning technique for accelerating the solution of the aforementioned linear system. The preconditioner is based on random feature maps, such as random Fourier features, which have recently emerged as a powerful technique for speeding up and scaling the training of kernel-based methods, such as kernel ridge regression, by resorting to approximations. However, random feature maps only provide crude approximations to the kernel function, so delivering state-of-the-art results by directly solving the approximated system requires the number of random features to be very large. We show that random feature maps can be much more effective in forming preconditioners, since under certain conditions a not-too-large number of random features is sufficient to yield an effective preconditioner. We empirically evaluate our method and show it is highly effective for datasets of up to one million training examples.


title: Indecomposable $1$-factorizations of the complete multigraph $\lambda K_{2n}$ for every $\lambda\leq 2n$

abstract: A $1$-factorization of the complete multigraph $\lambda K_{2n}$ is said to be indecomposable if it cannot be represented as the union of $1$-factorizations of $\lambda_0 K_{2n}$ and $(\lambda-\lambda_0) K_{2n}$, where $\lambda_0<\lambda$. It is said to be simple if no $1$-factor is repeated. For every $n\geq 9$ and for every $(n-2)/3\leq\lambda\leq 2n$, we construct an indecomposable $1$-factorization of $\lambda K_{2n}$ which is not simple. These $1$-factorizations provide simple and indecomposable $1$-factorizations of $\lambda K_{2s}$ for every $s\geq 18$ and $2\leq\lambda\leq 2\lfloor s/2\rfloor-1$. We also give a generalization of a result by Colbourn et al. which provides a simple and indecomposable $1$-factorization of $\lambda K_{2n}$, where $2n=p^m+1$, $\lambda=(p^m-1)/2$, $p$ prime.


title: A note on the Monge-Amp\`ere type equations with general source terms

abstract: In this paper we consider the generalised solutions to the Monge-Amp{\`{e}}re type equations with general source terms. We firstly prove the so-called comparison principle and then give some important propositions for the border of generalised solutions. Furthermore, we design well-posed finite element methods for the generalised solutions with the classical and weak Dirichlet boundary conditions respectively.


title: Chiral Separation by Flows: The Role of Flow Symmetry and Dimensionality

abstract: Separation of enantiomers by flows is a promising chiral resolution method using cost-effective microfluidics. Notwithstanding a number of experimental and numerical studies, a fundamental understanding still remains elusive, and an important question as to whether it is possible to specify common physical properties of flows that induce separation has not been addressed. Here, we study the separation of rigid chiral objects of an arbitrary shape induced by a linear flow field at low Reynolds numbers. Based on a symmetry property under parity inversion, we show that the rate-of-strain field is essential to drift the objects in opposite directions according to chirality. From eigenmode analysis, we also derive an analytic expression for the separation conditions which shows that the flow field should be quasi-two-dimensional for the precise and efficient resolutions of microscopic enantiomers. We demonstrate this prediction by Langevin dynamics simulations with hydrodynamic interactions fully implemented. Finally, we discuss the practical feasibility of the linear flow analysis, considering separations by a vortex flow or an extensional flow under a confining potential.


title: Rough I-statistical convergence of sequences

abstract: The concept of I-statistical convergence of sequence was first defined by Das et.al [2]. In this paper we introduce and study the notion of rough I-statistical convergence of sequence in normed linear Spaces. We also define the set of rough I-statistical limits of a sequence and discuss some topological properties of this set.


title: Sharper Bounds for Regularized Data Fitting

abstract: We study matrix sketching methods for regularized variants of linear regression, low rank approximation, and canonical correlation analysis. Our main focus is on sketching techniques which preserve the objective function value for regularized problems, which is an area that has remained largely unexplored. We study regularization both in a fairly broad setting, and in the specific context of the popular and widely used technique of ridge regularization; for the latter, as applied to each of these problems, we show algorithmic resource bounds in which the {\em statistical dimension} appears in places where in previous bounds the rank would appear. The statistical dimension is always smaller than the rank, and decreases as the amount of regularization increases. In particular, for the ridge low-rank approximation problem $\min_{Y,X} \lVert YX - A \rVert_F^2 + \lambda \lVert Y\rVert_F^2 + \lambda\lVert X \rVert_F^2$, where $Y\in\mathbb{R}^{n\times k}$ and $X\in\mathbb{R}^{k\times d}$, we give an approximation algorithm needing \[ O(\mathtt{nnz}(A)) + \tilde{O}((n+d)\varepsilon^{-1}k \min\{k, \varepsilon^{-1}\mathtt{sd}_\lambda(Y^*)\})+ \mathtt{poly}(\mathtt{sd}_\lambda(Y^*) \varepsilon^{-1}) \] time, where $s_{\lambda}(Y^*)\le k$ is the statistical dimension of $Y^*$, $Y^*$ is an optimal $Y$, $\varepsilon$ is an error parameter, and $\mathtt{nnz}(A)$ is the number of nonzero entries of $A$.This is faster than prior work, even when $\lambda=0$. We also study regularization in a much more general setting. For example, we obtain sketching-based algorithms for the low-rank approximation problem $\min_{X,Y} \lVert YX - A \rVert_F^2 + f(Y,X)$ where $f(\cdot,\cdot)$ is a regularizing function satisfying some very general conditions (chiefly, invariance under orthogonal transformations).


title: Executing Dynamic Data Rate Actor Networks on OpenCL Platforms

abstract: Heterogeneous computing platforms consisting of general purpose processors (GPPs) and graphics processing units (GPUs) have become commonplace in personal mobile devices and embedded systems. For years, programming of these platforms was very tedious and simultaneous use of all available GPP and GPU resources required low-level programming to ensure efficient synchronization and data transfer between processors. However, in the last few years several high-level programming frameworks have emerged, which enable programmers to describe applications by means of abstractions such as dataflow or Kahn process networks and leave parallel execution, data transfer and synchronization to be handled by the framework. Unfortunately, even the most advanced high-level programming frameworks have had shortcomings that limit their applicability to certain classes of applications. This paper presents a new, dataflow-flavored programming framework targeting heterogeneous platforms, and differs from previous approaches by allowing GPU-mapped actors to have data dependent consumption of inputs / production of outputs. Such flexibility is essential for configurable and adaptive applications that are becoming increasingly common in signal processing. In our experiments it is shown that this feature allows up to 5x increase in application throughput. The proposed framework is validated by application examples from the video processing and wireless communications domains. In the experiments the framework is compared to a well-known reference framework and it is shown that the proposed framework enables both a higher degree of flexibility and better throughput.


title: Feature Selection with the R Package MXM: Discovering Statistically-Equivalent Feature Subsets

abstract: The statistically equivalent signature (SES) algorithm is a method for feature selection inspired by the principles of constrained-based learning of Bayesian Networks. Most of the currently available feature-selection methods return only a single subset of features, supposedly the one with the highest predictive power. We argue that in several domains multiple subsets can achieve close to maximal predictive accuracy, and that arbitrarily providing only one has several drawbacks. The SES method attempts to identify multiple, predictive feature subsets whose performances are statistically equivalent. Under that respect SES subsumes and extends previous feature selection algorithms, like the max-min parent children algorithm. SES is implemented in an homonym function included in the R package MXM, standing for mens ex machina, meaning 'mind from the machine' in Latin. The MXM implementation of SES handles several data-analysis tasks, namely classification, regression and survival analysis. In this paper we present the SES algorithm, its implementation, and provide examples of use of the SES function in R. Furthermore, we analyze three publicly available data sets to illustrate the equivalence of the signatures retrieved by SES and to contrast SES against the state-of-the-art feature selection method LASSO. Our results provide initial evidence that the two methods perform comparably well in terms of predictive accuracy and that multiple, equally predictive signatures are actually present in real world data.


title: Phase diagram of the two-dimensional O(3) model from dual lattice simulations

abstract: We have simulated the asymptotically free two-dimensional O(3) model at nonzero chemical potential using the model's dual representation. We first demonstrate how the latter solves the sign (complex action) problem. The system displays a crossover at nonzero temperature, while at zero temperature it undergoes a quantum phase transition when mu reaches the particle mass (generated dynamically similar to QCD). The density follows a square root behavior universal for repulsive bosons in one spatial dimension. We have also measured the spin stiffness, known to be sensitive to the spatial correlation length, using different scaling trajectories to zero temperature and infinite size. It points to a dynamical critical exponent z=2. Comparisons to thermodynamic Bethe ansaetze are shown as well.


title: Higher-Order Multilevel Framework for ADER Scheme in Computational Aeroacoustics

abstract: The versatile Arbitrary-DERivative (ADER) scheme is cast in a multilevel framework (ML-ADER) for fast solution of system of linear hyperbolic partial differential equations. The solution is cycled through spatial operators of varying accuracy while retaining highest-order accuracy by the use of a forcing function. Accuracy analysis of the multilevel framework including in the ML-ADER form is carried out in time-domain as well as frequency-domain. Results are obtained for benchmark problems in computational aeroacoustics at a much reduced computational cost.


title: A new migration mechanism for oversized solutes in cubic lattices: correlation effects

abstract: Recent ab-initio calculations showed that oversized solute atoms (OSA) in BCC and FCC iron exhibit a very strong attraction with a nearby vacancy (V) at first neighbour distance (1NN). The attraction it is so large that the 1NN pair OSA+V is no longer stable and relaxes spontaneously towards a new configuration where the OSA sits in the middle of the bond, the two ends of which are decorated with two half-vacancies (V/2). Taking into account the formation of this complex V/2+OSA+V/2 in the migration process has never been done before: it requires a new formulation of correlation effects. This is the aim of the present contribution which establishes the expression of the diffusion coefficient of the OSA and evaluates exactly the correlation factor in BCC and FCC lattices, thanks to the introduction of macrojumps. The formalism is applied to the case of Y in BCC iron and its results are compared to those obtained by previous authors. It is shown that the smallness of the correlation factor is not due to the transport mechanism of the OSA but only to the jump frequencies of the vacancy in the vicinity of the OSA. This smallness of the correlation factor is however overcompensated by a large macrojump frequency, and, as a result, the Y atom diffuses more rapidly than iron at all temperatures by orders of magnitude in the two structures.


title: Policy Search with High-Dimensional Context Variables

abstract: Direct contextual policy search methods learn to improve policy parameters and simultaneously generalize these parameters to different context or task variables. However, learning from high-dimensional context variables, such as camera images, is still a prominent problem in many real-world tasks. A naive application of unsupervised dimensionality reduction methods to the context variables, such as principal component analysis, is insufficient as task-relevant input may be ignored. In this paper, we propose a contextual policy search method in the model-based relative entropy stochastic search framework with integrated dimensionality reduction. We learn a model of the reward that is locally quadratic in both the policy parameters and the context variables. Furthermore, we perform supervised linear dimensionality reduction on the context variables by nuclear norm regularization. The experimental results show that the proposed method outperforms naive dimensionality reduction via principal component analysis and a state-of-the-art contextual policy search method.


title: Central loops in random planar graphs

abstract: Random planar graphs appear in a variety of context and it is important for many different applications to be able to characterize their structure. Local quantities fail to give interesting information and it seems that path-related measures are able to convey relevant information about the organization of these structures. In particular, nodes with a large betweenness centrality (BC) display non-trivial patterns, such as central loops. We first discuss empirical results for different random planar graphs and we then propose a toy model which allows us to discuss the condition for the emergence of non-trivial patterns such as central loops. This toy model is made of a star network with $N_b$ branches of size $n$ and links of weight $1$, superimposed to a loop at distance $\ell$ from the center and with links of weight $w$. We estimate for this model the BC at the center and on the loop and we show that the loop can be more central than the origin if $w<w_c$ where the threshold of this transition scales as $w_c\sim n/N_b$. In this regime, there is an optimal position of the loop that scales as $\ell_{opt}\sim N_b w/4$. This simple model sheds some light on the organization of these random structures and allows us to discuss the effect of randomness on the centrality of loops. In particular, it suggests that the number and the spatial extension of radial branches are the crucial ingredients that control the existence of central loops.


title: Large-scale JPEG steganalysis using hybrid deep-learning framework

abstract: Adoption of deep learning in image steganalysis is still in its initial stage. In this paper we propose a generic hybrid deep-learning framework for JPEG steganalysis incorporating the domain knowledge behind rich steganalytic models. Our proposed framework involves two main stages. The first stage is hand-crafted, corresponding to the convolution phase and the quantization & truncation phase of the rich models. The second stage is a compound deep neural network containing multiple deep subnets in which the model parameters are learned in the training procedure. We provided experimental evidences and theoretical reflections to argue that the introduction of threshold quantizers, though disable the gradient-descent-based learning of the bottom convolution phase, is indeed cost-effective. We have conducted extensive experiments on a large-scale dataset extracted from ImageNet. The primary dataset used in our experiments contains 500,000 cover images, while our largest dataset contains five million cover images. Our experiments show that the integration of quantization and truncation into deep-learning steganalyzers do boost the detection performance by a clear margin. Furthermore, we demonstrate that our framework is insensitive to JPEG blocking artifact alterations, and the learned model can be easily transferred to a different attacking target and even a different dataset. These properties are of critical importance in practical applications.


title: High-pressure study of the ground- and superconducting-state properties of CeAu$_2$Si$_2$

abstract: The pressure-temperature-phase diagram of the new heavy-fermion superconductor CeAu$_2$Si$_2$ is markedly different from those studied previously. Indeed, superconductivity emerges, not on the verge, but deep inside the magnetic phase. In this context, we have carried out ac-calorimetry, resistivity and thermoelectric power measurements on a CeAu$_2$Si$_2$-single crystal under high-pressure. The principal novelties of this experiment are the observation of a new transition line, presumably a magnetic structure rearrangement at $T_{\mathrm{M}}^{\mathrm{mod}}$ inside the antiferromagnetic phase, and the occurrence of quantum critical behavior in resistivity linked to superconductivity. Strong non-Fermi-liquid behavior is observed around the maximum of superconductivity and enhanced scattering rates are observed close to both the emergence and the maximum of superconductivity. Intriguingly, $T_{\mathrm{M}}^{\mathrm{mod}}$ almost coincides with the onset of the superconducting transition over a broad pressure range, where mysteriously Tc increases with the strengthening of magnetism. A comparison of the features in CeAu$_2$Si$_2$ and its parent compounds CeCu$_2$Si$_2$ and CeCu$_2$Ge$_2$ plotted as function of the unit-cell volume leads us to promote that critical fluctuations of a valence crossover play a crucial role in the superconducting paring mechanism. Our study reveals new intriguing features of magnetism and illustrates a complex interplay between magnetism, quantum criticality and superconductivity.


title: On the Light Curves of AM CVn

abstract: Light curves of AM CVn are analyzed by decomposing them into their Fourier components. The amplitudes of the fundamental mode and overtones of the three components: the superhumps, the negative superhumps and the orbital variations, are found to be variable. This implies that variations in the shape of the observed light curve of AM CVn are not only due to the interference between those components, but also due to the variability of their parameters.


title: The number of satisfying assignments of random regular k-SAT formulas

abstract: Let $\Phi$ be a random $k$-SAT formula in which every variable occurs precisely $d$ times positively and $d$ times negatively. Assuming that $k$ is sufficiently large and that $d$ is slightly below the critical degree where the formula becomes unsatisfiable with high probability, we determine the limiting distribution of the logarithm of the number of satisfying assignments.


title: Competition in periodic media: II -- Segregative limit of pulsating fronts and "Unity is not Strength"-type result

abstract: This paper is concerned with the limit, as the interspecific competition rate goes to infinity, of pulsating front solutions in space-periodic media for a bistable two-species competition--diffusion Lotka--Volterra system. We distinguish two important cases: null asymptotic speed and non-null as-ymptotic speed. In the former case, we show the existence of a segregated stationary equilibrium. In the latter case, we are able to uniquely characterize the segregated pulsating front, and thus full convergence is proved. The segregated pulsating front solves an interesting free boundary problem. We also investigate the sign of the speed as a function of the parameters of the competitive system. We are able to determine it in full generality, with explicit conditions depending on the various parameters of the problem. In particular, if one species is sufficiently more motile or competitive than the other, then it is the invader. This is an extension of our previous work in space-homogeneous media.