title: Evaluating Reputation Management Schemes of Internet of Vehicles based on Evolutionary Game Theory

abstract: Conducting reputation management is very important for Internet of vehicles. However, most of the existing researches evaluate the effectiveness of their schemes with settled attacking behaviors in their simulation which cannot represent the scenarios in reality. In this paper, we propose to consider dynamical and diversity attacking strategies in the simulation of reputation management scheme evaluation. To that end, we apply evolutionary game theory to model the evolution process of malicious users' attacking strategies, and discuss the methodology of the evaluation simulations. We further apply our evaluation method to a reputation management scheme with multiple utility functions, and discuss the evaluation results. The results indicate that our evaluation method is able to depict the evolving process of the dynamic attacking strategies in a vehicular network, and the final state of the simulation could be used to quantify the protection effectiveness of the reputation management scheme.


title: A reconstruction theorem for varieties

abstract: We show that varieties of dimension at least 2 over infinite fields are determined as abstract schemes by their Zariski topological spaces together with the rational equivalence relation on the set of effective divisors. This gives a universal Torelli theorem in the sense of Bogomolov and Tschinkel. The proof relies heavily on a rational version of the classical Fundamental Theorem of Projective Geometry.


title: A new diagnostic tool for water, energy and entropy budgets in climate models

abstract: This work presents a novel diagnostic tool for studying the thermodynamics of the climate systems with a wide range of applications, from sensitivity studies to model tuning. It includes a number of modules for assessing the internal energy budget, the hydrological cycle, the Lorenz Energy Cycle and the material entropy production, respectively. The routine receives as inputs energy fluxes at surface and at the Top-of-Atmosphere (TOA), for the computation of energy budgets at Top-of-Atmosphere (TOA), at the surface, and in the atmosphere as a residual. Meridional enthalpy transports are also computed from the divergence of the zonal mean energy budget fluxes; location and intensity of peaks in the two hemispheres are then provided as outputs. Rainfall, snowfall and latent heat fluxes are received as inputs for computing the water mass and latent energy budgets. If a land-sea mask is provided, the required quantities are separately computed over continents and oceans. The diagnostic tool also computes the Lorenz Energy Cycle (LEC) and its storage/conversion terms as annual mean global and hemispheric values. In order to achieve this, one needs to provide as input three-dimensional daily fields of horizontal wind velocity and temperature in the troposphere. Two methods have been implemented for the computation of the material entropy production, one relying on the convergence of radiative heat fluxes in the atmosphere (indirect method), one combining the irreversible processes occurring in the climate system, particularly heat fluxes in the boundary layer, the hydrological cycle and the kinetic energy dissipation as retrieved from the residuals of the LEC. A version of the diagnostic tool is included in the Earth System Model eValuation Tool (ESMValTool) community diagnostics, in order to assess the performances of soon available CMIP6 model simulations.

categories: physics.ao-ph physics.class-ph physics.flu-dyn physics.geo-ph


title: A model of high-energy emission from jets of microquasars of Population III

abstract: Current research indicates that the radiation of the first stars alone would not have been suficient to ionize the intergalactic medium on long distances because of the high columnar density of the clouds in which they formed. It becomes necessary, then, to explore other alternatives for the reionization. The study of the contribution of the first accreting binary systems to the reionization and heating of the intergalactic medium requires the formulation of a concrete model for the first microquasars. We aim at constructing a consistent model for jets of microquasars where the donor star is of Population III. We calculate the spectral energy distribution of the radiation produced by the relativistic particles in the jets, within the framework of a lepto-hadronic model. In addition, we study the interaction between the jets and the early intergalactic medium. All calculations have been computed for a cosmological redshift value z = 10. We determine that jets of microquasars of Population III are sources of intense gamma radiation in regions of particle acceleration near the compact object and in the terminal jets.


title: A Novel Maneuvering Target Tracking Approach by Stochastic Volatility GARCH Model

abstract: In this paper, we introduce a new single model maneuvering target tracking approach using stochastic differential equation (SDE) based on GARCH volatility. The traditional input estimation (IE) techniques assume constant acceleration level which do not cover all the possible acceleration quintessence. In contrast, the multiple model (MM) algorithms that take care of some IE's shortcomings, are sensitive to the transition probability matrices. In this paper, an innovative model is proposed to overcome these drawbacks by using a new generalized dynamic modeling of acceleration and a Bayesian filter. We utilize SDE to model Markovian jump acceleration of a maneuvering target through GARCH process as the SDE volatility. In the proposed scheme, the original state and stochastic volatility (SV) are estimated simultaneously by a bootstrap particle filter (PF). We introduce the bootstrap resampling to obtain the statistical properties of a GARCH density. Due to the heavy-tailed nature of the GARCH distribution, the bootstrap PF is more effective in the presence of large errors that can occur in the state equation. We show analytically that the target tracking performance is improved by considering GARCH acceleration model. Finally, the effectiveness and capabilities of our proposed strategy (PF-AR-GARCH) are demonstrated and validated through simulation studies.


title: Free-standing 2-inch bulk GaN crystal fabrication by HVPE using a carbon buffer layer

abstract: A free-standing bulk gallium nitride layer with a thickness of 365 $\mu$m and a diameter of 50 mm was obtained by hydride vapor phase epitaxy on a sapphire substrate with a carbon buffer layer. The carbon buffer layer was deposited by thermal decomposition of methane $\textit{in situ}$ in the same process with the growth of a bulk GaN layer. The bulk GaN layer grown on the carbon buffer layer self-separated from the sapphire substrate during the cooling after the growth. The dislocation density was $8\cdot10^{6}$ cm$^{-2}$. The (0002) X-Ray rocking curve full width at half maximum was 164 arcsec.


title: Enhanced Balancing of Bias-Variance Tradeoff in Stochastic Estimation: A Minimax Perspective

abstract: Biased stochastic estimators, such as finite-differences for noisy gradient estimation, often contain parameters that need to be properly chosen to balance impacts from the bias and the variance. While the optimal order of these parameters in terms of the simulation budget can be readily established, the precise best values depend on model characteristics that are typically unknown in advance. We introduce a framework to construct new classes of estimators, based on judicious combinations of simulation runs on sequences of tuning parameter values, such that the estimators consistently outperform a given tuning parameter choice in the conventional approach, regardless of the unknown model characteristics. We argue the outperformance via what we call the asymptotic minimax risk ratio, obtained by minimizing the worst-case asymptotic ratio between the mean square errors of our estimators and the conventional one, where the worst case is over any possible values of the model unknowns. In particular, when the minimax ratio is less than 1, the calibrated estimator is guaranteed to perform better asymptotically. We identify this minimax ratio for general classes of weighted estimators, and the regimes where this ratio is less than 1. Moreover, we show that the best weighting scheme is characterized by a sum of two components with distinct decay rates. We explain how this arises from bias-variance balancing that combats the adversarial selection of the model constants, which can be analyzed via a tractable reformulation of a non-convex optimization problem.


title: Towards moderate overparameterization: global convergence guarantees for training shallow neural networks

abstract: Many modern neural network architectures are trained in an overparameterized regime where the parameters of the model exceed the size of the training dataset. Sufficiently overparameterized neural network architectures in principle have the capacity to fit any set of labels including random noise. However, given the highly nonconvex nature of the training landscape it is not clear what level and kind of overparameterization is required for first order methods to converge to a global optima that perfectly interpolate any labels. A number of recent theoretical works have shown that for very wide neural networks where the number of hidden units is polynomially large in the size of the training data gradient descent starting from a random initialization does indeed converge to a global optima. However, in practice much more moderate levels of overparameterization seems to be sufficient and in many cases overparameterized models seem to perfectly interpolate the training data as soon as the number of parameters exceed the size of the training data by a constant factor. Thus there is a huge gap between the existing theoretical literature and practical experiments. In this paper we take a step towards closing this gap. Focusing on shallow neural nets and smooth activations, we show that (stochastic) gradient descent when initialized at random converges at a geometric rate to a nearby global optima as soon as the square-root of the number of network parameters exceeds the size of the training data. Our results also benefit from a fast convergence rate and continue to hold for non-differentiable activations such as Rectified Linear Units (ReLUs).


title: Long-lived non-equilibrium superconductivity in a non-centrosymmetric Rashba semiconductor

abstract: We report non-equilibrium magnetodynamics in the Rashba-superconductor GeTe, which lacks inversion symmetry in the bulk. We find that at low temperature the system exhibits a non-equilibrium state, which decays on time scales that exceed conventional electronic scattering times by many orders of magnitude. This reveals a non-equilibrium magnetoresponse that is asymmetric under magnetic field reversal and, strikingly, induces a non-equilibrium superconducting state distinct from the equilibrium one. We develop a model of a Rashba system where non-equilibrium configurations relax on a finite timescale which captures the qualitative features of the data. We also obtain evidence for the slow dynamics in another non-superconducting Rashba system. Our work provides novel insights into the dynamics of non-centrosymmetric superconductors and Rashba systems in general.


title: Raman spectrum of 1T'-WTe2 under tensile strain: A first-principles prediction

abstract: Monolayer WTe2 attracts rapidly growing interests for its large-gap quantum spin Hall effect,which enables promising apllications in flexible logic devices. Due to one-dimensional W-W chains,1T'-WTe2 exhibits unique anisotropic structure and promising properties, which can be modified by simply applying strains. Based on the first-principles simulations, we show that phonon branch undergoes soft down to negative frequency at special q points under different critical strains, i.e., epsilon_a = 11.55 percent along a-axis (with W-W chains) direction, epsilon_b = 7.0 percent along b-axis direction and epsilon_ab = 8.44 percent along biaxial direction. Before each critical strain, the Raman-shift of A1_g, A3_g, and A4_g modes, corresponding to the main peaks in Raman spectra of 1T'-WTe2 , shows anisotropic response to uniaxial strain but most sensitive to biaxial strain. Interestingly, we find that the frequency shift of A3_g mode show parabolic characters of strained 1T'-WTe2, then we split it into two parts and it shows a Raman-shift transition at about 5 percent strains. While for the A1_g and A4_g modes,the frequencies change linearly. Through careful structure and vibration analysis, we try to explain these Raman irregularity in strained 1T'-WTe2.


title: Channel-Statistics-Based Hybrid Precoding for Millimeter-Wave MIMO Systems With Dynamic Subarrays

abstract: This paper investigates the hybrid precoding design for millimeter wave (mmWave) multiple-input multiple-output (MIMO) systems with finite-alphabet inputs. The mmWave MIMO system employs partially-connected hybrid precoding architecture with dynamic subarrays, where each radio frequency (RF) chain is connected to a dynamic subset of antennas. We consider the design of analog and digital precoders utilizing statistical and/or mixed channel state information (CSI), which involve solving an extremely difficult problem in theory: First, designing the optimal partition of antennas over RF chains is a combinatorial optimization problem, whose optimal solution requires an exhaustive search over all antenna partitioning solutions; Second, the average mutual information under mmWave MIMO channels lacks closed-form expression and involves prohibitive computational burden; Third, the hybrid precoding problem with given partition of antennas is nonconvex with respect to the analog and digital precoders. To address these issues, this study first presents a simple criterion and the corresponding low complexity algorithm to design the optimal partition of antennas using statistical CSI. Then it derives the lower bound and its approximation for the average mutual information, in which the computational complexity is greatly reduced compared to calculating the average mutual information directly. In addition, it also shows that the lower bound with a constant shift offers a very accurate approximation to the average mutual information. This paper further proposes utilizing the lower bound approximation as a low-complexity and accurate alternative for developing a manifold-based gradient ascent algorithm to find near optimal analog and digital precoders. Several numerical results are provided to show that our proposed algorithm outperforms existing hybrid precoding algorithms.


title: Chirality relaxation in low-temperature strongly Rashba-coupled systems

abstract: We study the relaxation dynamics of non-equilibrium chirality distributions of charge carriers in Rashba systems. We find that at low temperature inter-Rashba band transitions become suppressed due to the combined effect of the Rashba momentum split and the chiral spin texture of a Rashba system. Specifically, we show that momentum exchange between carriers and the phonon bath is effectively absent at temperatures where the momentum of thermal phonons is less than twice the Rashba momentum. This allows us to identify inter-carrier scattering as the dominant process by which non-equilibrium chirality distributions relax. We show that the magnitude of inter-carrier scattering is strongly influenced by the opposing spin structure of the Rashba bands. Finally, we provide an explicit result for the inter-band relaxation timescale associated with inter-carrier Coulomb scattering. We develop a general framework and assess its implications for GeTe, a bulk Rashba semiconductor with a strong Rashba momentum split.


title: High dimensionality: The latest challenge to data analysis

abstract: The advent of modern technology, permitting the measurement of thousands of characteristics simultaneously, has given rise to floods of data characterized by many large or even huge datasets. This new paradigm presents extraordinary challenges to data analysis and the question arises: how can conventional data analysis methods, devised for moderate or small datasets, cope with the complexities of modern data? The case of high dimensional data is particularly revealing of some of the drawbacks. We look at the case where the number of characteristics measured in an object is at least the number of observed objects and conclude that this configuration leads to geometrical and mathematical oddities and is an insurmountable barrier for the direct application of traditional methodologies. If scientists are going to ignore fundamental mathematical results arrived at in this paper and blindly use software to analyze data, the results of their analyses may not be trustful, and the findings of their experiments may never be validated. That is why new methods together with the wise use of traditional approaches are essential to progress safely through the present reality.


title: Quantum Engineered Kondo Lattices

abstract: Recent advances in atomic manipulation techniques have provided a novel bottom-up approach to investigating the unconventional properties and complex phases of strongly correlated electron materials. By engineering artificial condensed matter systems containing tens to thousands of atoms with tailored electronic or magnetic properties, it has become possible to explore how quantum many-body effects---whose existence lies at the heart of strongly correlated materials---emerge as the size of a system is increased from the nanoscale to the mesoscale. Here we investigate both theoretically and experimentally the quantum engineering of nanoscopic Kondo lattices -- Kondo droplets -- that exemplify nanoscopic replicas of heavy-fermion materials. We demonstrate that by changing a droplet's real space geometry, it is possible to not only create coherently coupled Kondo droplets whose properties asymptotically approach those of a quantum-coherent Kondo lattice, but also to markedly increase or decrease the droplet's Kondo temperature. Furthermore we report on the discovery of a new quantum phenomenon -- the Kondo echo -- a signature of droplets containing Kondo holes functioning as direct probes of spatially extended, quantum-coherent Kondo cloud correlations.


title: Resolving the energy levels of a nanomechanical oscillator

abstract: The coherent states that describe the classical motion of a mechanical oscillator do not have well-defined energy, but are rather quantum superpositions of equally-spaced energy eigenstates. Revealing this quantized structure is only possible with an apparatus that measures the mechanical energy with a precision greater than the energy of a single phonon, $\hbar\omega_\text{m}$. One way to achieve this sensitivity is by engineering a strong but nonresonant interaction between the oscillator and an atom. In a system with sufficient quantum coherence, this interaction allows one to distinguish different phonon number states by resolvable differences in the atom's transition frequency. Such dispersive measurements have been studied in cavity and circuit quantum electrodynamics where experiments using real and artificial atoms have resolved the photon number states of cavities. Here, we report an experiment where an artificial atom senses the motional energy of a driven nanomechanical oscillator with sufficient sensitivity to resolve the quantization of its energy. To realize this, we build a hybrid platform that integrates nanomechanical piezoelectric resonators with a microwave superconducting qubit on the same chip. We excite phonons with resonant pulses of varying amplitude and probe the resulting excitation spectrum of the qubit to observe phonon-number-dependent frequency shifts $\approx 5$ times larger than the qubit linewidth. Our result demonstrates a fully integrated platform for quantum acoustics that combines large couplings, considerable coherence times, and excellent control over the mechanical mode structure. With modest experimental improvements, we expect our approach will make quantum nondemolition measurements of phonons an experimental reality, leading the way to new quantum sensors and information processing approaches that use chip-scale nanomechanical devices.


title: Combating the Distance Problem in the Millimeter Wave and Terahertz Frequency Bands

abstract: In the millimeter wave (30-300 GHz) and Terahertz (0.1-10 THz) frequency bands, high spreading loss and molecular absorption often limit the signal transmission distance and coverage range. In this paper, four directions to tackle the crucial problem of distance limitation are investigated, namely, a physical layer distance-aware design, ultra-massive MIMO communication, reflectarrays, and intelligent surfaces. Additionally, the potential joint design of these technologies is proposed to combine the benefits and possibly further extend the communication distance. Qualitative analyses and quantitative simulations are provided to illustrate the benefits of the proposed techniques and demonstrate the feasibility of mm-wave and THz band communications up to 100 meters in both line-of-sight and non-line-of-sight areas.


title: Search for Higgs-like particles produced in association with bottom quarks in proton-antiproton collisions

abstract: We report on a search for a spin-zero non-standard-model particle in proton-antiproton collisions collected by the Collider Detector at Fermilab at a center-of-mass-energy of 1.96 TeV. This particle, the $\phi$ boson, is expected to decay into a bottom-antibottom quark pair and to be produced in association with at least one bottom quark. The data sample consists of events with three jets identified as initiated by bottom quarks and corresponds to $5.4~\text{fb}^{-1}$ of integrated luminosity. In each event, the invariant mass of the two most energetic jets is studied by looking for deviations from the multijet background, which is modeled using data. No evidence is found for such particle. Exclusion upper limits ranging from 20 to 2 pb are set for the product of production cross sections times branching fraction for hypothetical $\phi$ boson with mass between 100 and 300 GeV/$c^2$. These are the most stringent constraints to date.


title: Forensic Similarity for Digital Images

abstract: In this paper we introduce a new digital image forensics approach called forensic similarity, which determines whether two image patches contain the same forensic trace or different forensic traces. One benefit of this approach is that prior knowledge, e.g. training samples, of a forensic trace are not required to make a forensic similarity decision on it in the future. To do this, we propose a two part deep-learning system composed of a CNN-based feature extractor and a three-layer neural network, called the similarity network. This system maps pairs of image patches to a score indicating whether they contain the same or different forensic traces. We evaluated system accuracy of determining whether two image patches were 1) captured by the same or different camera model, 2) manipulated by the same or different editing operation, and 3) manipulated by the same or different manipulation parameter, given a particular editing operation. Experiments demonstrate applicability to a variety of forensic traces, and importantly show efficacy on "unknown" forensic traces that were not used to train the system. Experiments also show that the proposed system significantly improves upon prior art, reducing error rates by more than half. Furthermore, we demonstrated the utility of the forensic similarity approach in two practical applications: forgery detection and localization, and database consistency verification.


title: Detecting and determining preserved measures and integrals of birational maps

abstract: In this paper we use the method of discrete Darboux polynomials to calculate preserved measures and integrals of rational maps. The approach is based on the use of cofactors and Darboux polynomials and relies on the use of symbolic algebra tools. Given sufficient computing power, most, if not all, rational preserved integrals can be found (and even some non-rational ones). We show, in a number of examples, how it is possible to use this method to both determine and detect preserved measures and integrals of the considered rational maps. Many of the examples arise from the Kahan-Hirota-Kimura discretization of completely integrable systems of ordinary differential equations.


title: The Complexity of Making the Gradient Small in Stochastic Convex Optimization

abstract: We give nearly matching upper and lower bounds on the oracle complexity of finding $\epsilon$-stationary points ($\| \nabla F(x) \| \leq\epsilon$) in stochastic convex optimization. We jointly analyze the oracle complexity in both the local stochastic oracle model and the global oracle (or, statistical learning) model. This allows us to decompose the complexity of finding near-stationary points into optimization complexity and sample complexity, and reveals some surprising differences between the complexity of stochastic optimization versus learning. Notably, we show that in the global oracle/statistical learning model, only logarithmic dependence on smoothness is required to find a near-stationary point, whereas polynomial dependence on smoothness is necessary in the local stochastic oracle model. In other words, the separation in complexity between the two models can be exponential, and that the folklore understanding that smoothness is required to find stationary points is only weakly true for statistical learning. Our upper bounds are based on extensions of a recent "recursive regularization" technique proposed by Allen-Zhu (2018). We show how to extend the technique to achieve near-optimal rates, and in particular show how to leverage the extra information available in the global oracle model. Our algorithm for the global model can be implemented efficiently through finite sum methods, and suggests an interesting new computational-statistical tradeoff.


title: Gauge invariant renormalizability of quantum gravity

abstract: Using the Batalin-Vilkovisky technique and the background field method the proof of gauge invariant renormalizability is elaborated for a generic model of quantum gravity which is diffeomorphism invariant and has no other, potentially anomalous, symmetries. The gauge invariant renormalizability means that in all orders of loop expansion of the quantum effective action one can control deformations of the generators of gauge transformations which leave such an action invariant. In quantum gravity this means that one can maintain general covariance of the divergent part of effective action when the mean quantum field, ghosts and antifields are switched off.


title: Privacy-Utility Trade-off of Linear Regression under Random Projections and Additive Noise

abstract: Data privacy is an important concern in machine learning, and is fundamentally at odds with the task of training useful learning models, which typically require the acquisition of large amounts of private user data. One possible way of fulfilling the machine learning task while preserving user privacy is to train the model on a transformed, noisy version of the data, which does not reveal the data itself directly to the training procedure. In this work, we analyze the privacy-utility trade-off of two such schemes for the problem of linear regression: additive noise, and random projections. In contrast to previous work, we consider a recently proposed notion of differential privacy that is based on conditional mutual information (MI-DP), which is stronger than the conventional $(\epsilon, \delta)$-differential privacy, and use relative objective error as the utility metric. We find that projecting the data to a lower-dimensional subspace before adding noise attains a better trade-off in general. We also make a connection between privacy problem and (non-coherent) SIMO, which has been extensively studied in wireless communication, and use tools from there for the analysis. We present numerical results demonstrating the performance of the schemes.


title: Wideband polarized radio emission from the newly revived magnetar XTE J1810$-$197

abstract: The anomalous X-ray pulsar XTE J1810$-$197 was the first magnetar found to emit pulsed radio emission. After spending almost a decade in a quiescent, radio-silent state, the magnetar was reported to have undergone a radio outburst in December, 2018. We observed radio pulsations from XTE J1810$-$197 during this early phase of its radio revival using the Ultra-Wideband Low receiver system of the Parkes radio telescope, obtaining wideband (704 MHz to 4032 MHz) polarization pulse profiles, single pulses and flux density measurements. Dramatic changes in polarization and rapid variations of the position angle of linear polarization across the main pulse and in time have been observed. The pulse profile exhibits similar structures throughout our three observations (over a week time scale), displaying a small amount of profile evolution in terms of polarization and pulse width across the wideband. We measured a flat radio spectrum across the band with a positive spectral index, in addition to small levels of flux and spectral index variability across our observing span. The observed wideband polarization properties are significantly different compared to those taken after the 2003 outburst, and therefore provide new information about the origin of radio emission.


title: Fragmentation and inefficiencies in US equity markets: Evidence from the Dow 30

abstract: Using the most comprehensive source of commercially available data on the US National Market System, we analyze all quotes and trades associated with Dow 30 stocks in 2016 from the vantage point of a single and fixed frame of reference. We find that inefficiencies created in part by the fragmentation of the equity marketplace are relatively common and persist for longer than what physical constraints may suggest. Information feeds reported different prices for the same equity more than 120 million times, with almost 64 million dislocation segments featuring meaningfully longer duration and higher magnitude. During this period, roughly 22% of all trades occurred while the SIP and aggregated direct feeds were dislocated. The current market configuration resulted in a realized opportunity cost totaling over $160 million when compared with a single feed, single exchange alternative---a conservative estimate that does not take into account intra-day offsetting events.


title: Scaling of inefficiencies in the U.S. equity markets: Evidence from three market indices and more than 2900 securities

abstract: Using the most comprehensive, commercially-available dataset of trading activity in U.S. equity markets, we catalog and analyze quote dislocations between the SIP National Best Bid and Offer (NBBO) and a synthetic BBO constructed from direct feeds. We observe a total of over 3.1 billion dislocation segments in the Russell 3000 during trading in 2016, roughly 525 per second of trading. However, these dislocations do not occur uniformly throughout the trading day. We identify a characteristic structure that features more dislocations near the open and close. Additionally, around 23% of observed trades executed during dislocations. These trades may have been impacted by stale information, leading to estimated opportunity costs on the order of $ 2 billion USD. A subset of the constituents of the S&P 500 index experience the greatest amount of opportunity cost and appear to drive inefficiencies in other stocks. These results quantify impacts of the physical structure of the U.S. National Market System.


title: Analysis of Baseline Evolutionary Algorithms for the Packing While Travelling Problem

abstract: The performance of base-line Evolutionary Algorithms (EAs) on combinatorial problems has been studied rigorously. From the theoretical viewpoint, the literature extensively investigates the linear problems, while the theoretical analysis of the non-linear problems is still far behind. In this paper, variations of the Packing While Travelling (PWT) -- also known as the non-linear knapsack problem -- are studied as an attempt to analyse the behaviour of EAs on non-linear problems from theoretical perspective. We investigate PWT for two cities and $n$ items with correlated weights and profits, using single-objective and multi-objective algorithms. Our results show that RLS\_swap, which differs from the classical RLS by having the ability to swap two bits in one iteration, finds the optimal solution in $O(n^3)$ expected time. We also study an enhanced version of GSEMO, which a specific selection operator to deal with exponential population size, and prove that it finds the Pareto front in the same asymptotic expected time. In the case of uniform weights, (1+1)~EA is able to find the optimal solution in expected time $O(n^2\log{(\max\{n,p_{\max}\})})$, where $p_{\max}$ is the largest profit of the given items. We also perform an experimental analysis to complement our theoretical investigations and provide additional insights into the runtime behavior.


title: A positivity-preserving finite volume element method for anisotropic diffusion problems on quadrilateral meshes

abstract: In this paper, we propose a nonlinear positivity-preserving finite volume element(FVE) scheme for anisotropic diffusion problems on quadrilateral meshes. Based on an overlapping dual partition, the one-sided flux is approximated by the iso-parametric bilinear element. A positivity-preserving nonlinear scheme with vertex-centered unknowns is obtained by a new two-point flux technique, which avoids the convex decomposition of co-normals and the introduction of intermediate unknowns. The existence of a solution is proved for this nonlinear system by applying the Brouwer's theorem. Numerical results show that the proposed positivity-preserving scheme is effective on distorted quadrilateral meshes and has approximate second-order accuracy for both isotropic and anisotropic diffusion problems. Moreover, the presented scheme is applied on an equilibrium radiation diffusion problem with discontinuous coefficients. Numerical results show that the new scheme is much more efficient than the standard FVE method.


title: Adapting Quality Assurance to Adaptive Systems: The Scenario Coevolution Paradigm

abstract: From formal and practical analysis, we identify new challenges that self-adaptive systems pose to the process of quality assurance. When tackling these, the effort spent on various tasks in the process of software engineering is naturally re-distributed. We claim that all steps related to testing need to become self-adaptive to match the capabilities of the self-adaptive system-under-test. Otherwise, the adaptive system's behavior might elude traditional variants of quality assurance. We thus propose the paradigm of scenario coevolution, which describes a pool of test cases and other constraints on system behavior that evolves in parallel to the (in part autonomous) development of behavior in the system-under-test. Scenario coevolution offers a simple structure for the organization of adaptive testing that allows for both human-controlled and autonomous intervention, supporting software engineering for adaptive systems on a procedural as well as technical level.


title: Advances in the optimization of silicon-based thermoelectrics: a theory perspective

abstract: Thermoelectric devices convert temperature gradients into electrical power and vice versa, thus enabling energy scavenging from waste heat, sensing and cooling. Yet, many of these attractive applications are hindered by the limited efficiency of thermoelectric materials, especially in the low temperature regime. This review provides a summary of the recent advances in the design of new efficient silicon--based thermoelectric materials through nanostructuring, alloying and chemical optimization, emphasizing the contribution from theory and atomistic modeling.


title: Using Approximate Models in Robot Learning

abstract: Trajectory following is one of the complicated control problems when its dynamics are nonlinear, stochastic and include a large number of parameters. The problem has significant difficulties including a large number of trials required for data collection and a massive volume of computations required to find a closed-loop controller for high dimensional and stochastic domains. For solving this type of problem, if we have an appropriate reward function and dynamics model; finding an optimal control policy is possible by using model-based reinforcement learning and optimal control algorithms. However, defining an accurate dynamics model is not possible for complicated problems. Pieter Abbeel and Andrew Ng recently presented an algorithm that requires only an approximate model and only a small number of real-life trials. This algorithm has broad applicability; however, there are some problems regarding the convergence of the algorithm. In this research, required modifications are presented that provide more powerful assurance for converging to optimal control policy. Also updated algorithm is implemented to evaluate the efficiency of the new algorithm by comparing the acquired results with human expert performance. We are using differential dynamic programming (DDP) as the locally trajectory optimizer, and a 2D dynamics and kinematics simulator is used to evaluate the accuracy of the presented algorithm.


title: Rethinking Generative Mode Coverage: A Pointwise Guaranteed Approach

abstract: Many generative models have to combat $\textit{missing modes}$. The conventional wisdom to this end is by reducing through training a statistical distance (such as $f$-divergence) between the generated distribution and provided data distribution. But this is more of a heuristic than a guarantee. The statistical distance measures a $\textit{global}$, but not $\textit{local}$, similarity between two distributions. Even if it is small, it does not imply a plausible mode coverage. Rethinking this problem from a game-theoretic perspective, we show that a complete mode coverage is firmly attainable. If a generative model can approximate a data distribution moderately well under a global statistical distance measure, then we will be able to find a mixture of generators that collectively covers $\textit{every}$ data point and thus $\textit{every}$ mode, with a lower-bounded generation probability. Constructing the generator mixture has a connection to the multiplicative weights update rule, upon which we propose our algorithm. We prove that our algorithm guarantees complete mode coverage. And our experiments on real and synthetic datasets confirm better mode coverage over recent approaches, ones that also use generator mixtures but rely on global statistical distances.


title: Identity Crisis: Memorization and Generalization under Extreme Overparameterization

abstract: We study the interplay between memorization and generalization of overparameterized networks in the extreme case of a single training example and an identity-mapping task. We examine fully-connected and convolutional networks (FCN and CNN), both linear and nonlinear, initialized randomly and then trained to minimize the reconstruction error. The trained networks stereotypically take one of two forms: the constant function (memorization) and the identity function (generalization). We formally characterize generalization in single-layer FCNs and CNNs. We show empirically that different architectures exhibit strikingly different inductive biases. For example, CNNs of up to 10 layers are able to generalize from a single example, whereas FCNs cannot learn the identity function reliably from 60k examples. Deeper CNNs often fail, but nonetheless do astonishing work to memorize the training output: because CNN biases are location invariant, the model must progressively grow an output pattern from the image boundaries via the coordination of many layers. Our work helps to quantify and visualize the sensitivity of inductive biases to architectural choices such as depth, kernel width, and number of channels.


title: Differential Description Length for Hyperparameter Selection in Machine Learning

abstract: This paper introduces a new method for model selection and more generally hyperparameter selection in machine learning. Minimum description length (MDL) is an established method for model selection, which is however not directly aimed at minimizing generalization error, which is often the primary goal in machine learning. The paper demonstrates a relationship between generalization error and a difference of description lengths of the training data; we call this difference differential description length (DDL). This allows prediction of generalization error from the training data alone by performing encoding of the training data. DDL can then be used for model selection by choosing the model with the smallest predicted generalization error. We show how this method can be used for linear regression and neural networks and deep learning. Experimental results show that DDL leads to smaller generalization error than cross-validation and traditional MDL and Bayes methods.


title: Impact of Metal ns2 Lone Pair on Luminescence Efficiency in Low-Dimensional Halide Perovskites

abstract: Based on first-principles calculations, we show that chemically active metal ns2 lone pairs play an important role in exciton relaxation and dissociation in low-dimensional halide perovskites. We studied excited-state properties of several recently discovered luminescent all-inorganic and hybrid organic-inorganic zero-dimensional (0D) Sn and Pb halides. The results show that, despite the similarity in ground-state electronic structure between Sn and Pb halide perovskites, the chemically more active Sn2+ lone pair leads to stronger excited-state structural distortion and larger Stokes shift in Sn halides. The enhanced Stokes shift hinders excitation energy transport, which reduces energy loss to defects and increases the photoluminescence quantum efficiency (PLQE). The presence of the ns2 metal cations in the 0D halide perovskites also promotes the exciton dissociation into electron and hole polarons especially in all-inorganic compounds, in which the coupling between metal-halide clusters is significant.


title: Efficient Bayesian shape-restricted function estimation with constrained Gaussian process priors

abstract: This article revisits the problem of Bayesian shape-restricted inference in the light of a recently developed approximate Gaussian process that admits an equivalent formulation of the shape constraints in terms of the basis coefficients. We propose a strategy to efficiently sample from the resulting constrained posterior by absorbing a smooth relaxation of the constraint in the likelihood and using circulant embedding techniques to sample from the unconstrained modified prior. We additionally pay careful attention to mitigate the computational complexity arising from updating hyperparameters within the covariance kernel of the Gaussian process. The developed algorithm is shown to be accurate and highly efficient in simulated and real data examples.


title: The substructures in the local stellar halo from Gaia and LAMOST

abstract: Based on the second Gaia data release (Gaia DR2) and spectroscopy from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data, we identified 20,089 halo stars kinematically and chemically. The halo streams in the solar neighborhood could be detected in the space of energy and angular momentum. We reshuffle the velocities of these stars to determine the significance of substructure. Finally, we identify 4 statistically significant substructures that are labeled GL-1 through 4. Among these substructures, GL-1 is previously known stream ("N2" stream) and the rest 3 substructures are new. These substructures may be the debris of dwarf galaxies accretion event, their dynamical and chemical information can help to understand the history of the Milky Way.


title: Assessing Solution Quality of 3SAT on a Quantum Annealing Platform

abstract: When solving propositional logic satisfiability (specifically 3SAT) using quantum annealing, we analyze the effect the difficulty of different instances of the problem has on the quality of the answer returned by the quantum annealer. A high-quality response from the annealer in this case is defined by a high percentage of correct solutions among the returned answers. We show that the phase transition regarding the computational complexity of the problem, which is well-known to occur for 3SAT on classical machines (where it causes a detrimental increase in runtime), persists in some form (but possibly to a lesser extent) for quantum annealing.


title: Neural network models and deep learning - a primer for biologists

abstract: Originally inspired by neurobiology, deep neural network models have become a powerful tool of machine learning and artificial intelligence, where they are used to approximate functions and dynamics by learning from examples. Here we give a brief introduction to neural network models and deep learning for biologists. We introduce feedforward and recurrent networks and explain the expressive power of this modeling framework and the backpropagation algorithm for setting the parameters. Finally, we consider how deep neural networks might help us understand the brain's computations.


title: Self-adaptive Single and Multi-illuminant Estimation Framework based on Deep Learning

abstract: Illuminant estimation plays a key role in digital camera pipeline system, it aims at reducing color casting effect due to the influence of non-white illuminant. Recent researches handle this task by using Convolution Neural Network (CNN) as a mapping function from input image to a single illumination vector. However, global mapping approaches are difficult to deal with scenes under multi-light-sources. In this paper, we proposed a self-adaptive single and multi-illuminant estimation framework, which includes the following novelties: (1) Learning local self-adaptive kernels from the entire image for illuminant estimation with encoder-decoder CNN structure; (2) Providing confidence measurement for the prediction; (3) Clustering-based iterative fitting for computing single and multi-illumination vectors. The proposed global-to-local aggregation is able to predict multi-illuminant regionally by utilizing global information instead of training in patches, as well as brings significant improvement for single illuminant estimation. We outperform the state-of-the-art methods on standard benchmarks with the largest relative improvement of 16%. In addition, we collect a dataset contains over 13k images for illuminant estimation and evaluation. The code and dataset is available on https://github.com/LiamLYJ/KPF_WB


title: Simultaneously Learning Vision and Feature-based Control Policies for Real-world Ball-in-a-Cup

abstract: We present a method for fast training of vision based control policies on real robots. The key idea behind our method is to perform multi-task Reinforcement Learning with auxiliary tasks that differ not only in the reward to be optimized but also in the state-space in which they operate. In particular, we allow auxiliary task policies to utilize task features that are available only at training-time. This allows for fast learning of auxiliary policies, which subsequently generate good data for training the main, vision-based control policies. This method can be seen as an extension of the Scheduled Auxiliary Control (SAC-X) framework. We demonstrate the efficacy of our method by using both a simulated and real-world Ball-in-a-Cup game controlled by a robot arm. In simulation, our approach leads to significant learning speed-ups when compared to standard SAC-X. On the real robot we show that the task can be learned from-scratch, i.e., with no transfer from simulation and no imitation learning. Videos of our learned policies running on the real robot can be found at https://sites.google.com/view/rss-2019-sawyer-bic/.


title: Sampling networks by nodal attributes

abstract: In a social network individuals or nodes connect to other nodes by choosing one of the channels of communication at a time to re-establish the existing social links. Since available data sets are usually restricted to a limited number of channels or layers, these autonomous decision making processes by the nodes constitute the sampling of a multiplex network leading to just one (though very important) example of sampling bias caused by the behavior of the nodes. We develop a general setting to get insight and understand the class of network sampling models, where the probability of sampling a link in the original network depends on the attributes $h$ of its adjacent nodes. Assuming that the nodal attributes are independently drawn from an arbitrary distribution $\rho(h)$ and that the sampling probability $r(h_i , h_j)$ for a link $ij$ of nodal attributes $h_i$ and $h_j$ is also arbitrary, we derive exact analytic expressions of the sampled network for such network characteristics as the degree distribution, degree correlation, and clustering spectrum. The properties of the sampled network turn out to be sums of quantities for the original network topology weighted by the factors stemming from the sampling. Based on our analysis, we find that the sampled network may have sampling-induced network properties that are absent in the original network, which implies the potential risk of a naive generalization of the results of the sample to the entire original network. We also consider the case, when neighboring nodes have correlated attributes to show how to generalize our formalism for such sampling bias and we get good agreement between the analytic results and the numerical simulations.


title: Discorrelation between primes in short intervals and polynomial phases

abstract: Let $H = N^{\theta}, \theta > 2/3$ and $k \geq 1$. We obtain estimates for the following exponential sum over primes in short intervals: \[ \sum_{N < n \leq N+H} \Lambda(n) e(g(n)), \] where $g$ is a polynomial of degree $k$. As a consequence of this in the special case $g(n) = \alpha n^k$, we deduce a short interval version of the Waring-Goldbach problem.


title: Experimental and Theoretical Study of Thermodynamic Effects in a Quantum Annealer

abstract: Quantum devices are affected by intrinsic and environmental noises. An in-depth characterization of noise effects is essential for exploiting noisy quantum computing. To this end, we studied the energy dissipative behavior of a quantum annealer via experiments and numerical simulations. Our investigation adopts a recently proposed technique that interpolates between pure quantum dynamics and pure thermodynamics. Experiments were conducted on a quantum annealer with an anneal pause function, which inserts a thermal relaxation period into the annealing schedule by pausing the transverse field, which is a source of quantum fluctuation. After investigating the special Hamiltonian that characterizes the quantum thermodynamics of the system, we then observed enhancement of thermodynamic signature depending on the anneal pause parameter. The time development of the state vector, observed in the open quantum simulation, provides rich information for investigating phenomena beyond energy-gap analysis. We identified a special eigenstate bridges ground states far-separated in Hilbert space and the transfer probabilities from one ground state to another. This finding can improve the sampling uniformity by reducing the sampling bias in finding the classical ground states in the quantum annealer. Our study does not only characterize the open quantum phenomenon of the specific Hamiltonian but also demonstrates the usefulness of the method in investigating noisy quantum devices.


title: Tuning Dissipation and Excitations in Superfluid Fermi Gases with a Moving Impurity

abstract: We develop a method to extract the dissipation for a heavy moving impurity immersed in superfluid Fermi gases. The drag force is derived analytically. As a reward, we are able to extract the dynamical structure factor, from which density excitations of the system is carefully examined. We show that dissipations through drag force is associated with two types of excitations, one being single-particle and the other being collective. We map out the critical velocity for dissipation across the BEC-BCS crossover, consistent with existing experiments. For a magnetic impurity, we show that the dissipation is immune to collective excitations. Our study clearly manifests that dissipation and associated excitations can be controlled by coupling superfluid Fermi gases with a moving impurity, and paves the way for further exploring intriguing realm of nonequilibrium phenomena and dissipation dynamics.


title: Towards a Better Indicator for Cache Timing Channels

abstract: Recent studies highlighting the vulnerability of computer architecture to information leakage attacks have been a cause of significant concern. Among the various classes of microarchitectural attacks, cache timing channels are especially worrisome since they have the potential to compromise users' private data at high bit rates. Prior works have demonstrated the use of cache miss patterns to detect these attacks. We find that cache miss traces can be easily spoofed and thus they may not be able to identify smarter adversaries. In this work, we show that \emph{cache occupancy}, which records the number of cache blocks owned by a specific process, can be leveraged as a stronger indicator for the presence of cache timing channels. We observe that the modulation of cache access latency in timing channels can be recognized through analyzing pairwise cache occupancy patterns. Our experimental results show that cache occupancy patterns cannot be easily obfuscated even by advanced adversaries that successfully evade cache miss-based detection.


title: Theoretical predictions on polarization asymmetry for Drell-Yan process with spin-one deuteron and tensor-polarized structure function $b_1$

abstract: We report recent theoretical progress on a polarization asymmetry in the proton-deuteron Drell-Yan process with a polarized-deuteron target and the tensor-polarized structure function $b_1$. Experimental measurements are possible at JLab for $b_1$ and at Fermilab for the Drell-Yan process. First, we show a theoretical estimate for the proton-deuteron Drell-Yan asymmetry in the Fermilab-E1039 experiment. We evolved tensor-polarized parton distribution functions, which explain existing HERMES $b_1$ data, at $Q^2=2.5$ GeV$^2$ to the $Q^2$ range of the Fermilab Drell-Yan measurements. Then, we predicted that the asymmetry is of the order of a few percent. The Drell-Yan experiment has an advantage to probe the tensor-polarized antiquark distributions, which were suggested by the HERMES experiment as a finite sum for $b_1$ ($\int dx b_1 (x) \ne 0$). Second, we predicted $b_1$ for the JLab experiment by the standard convolution model of the deuteron. Our theoretical $b_1$ structure function seems to be much different from the HERMES data. Furthermore, a significant distribution exists at very large $x$ ($>1$) beyond the kinematical limit $x_{max}=1$ for the proton. Because the standard deuteron-model estimate is much different from the HERMES data, there could be an interesting development as a new hadron-physics field if future JLab data will be much different from our conventional prediction.


title: Automated Segmentation of the Optic Disk and Cup using Dual-Stage Fully Convolutional Networks

abstract: Automated segmentation of the optic cup and disk on retinal fundus images is fundamental for the automated detection / analysis of glaucoma. Traditional segmentation approaches depend heavily upon hand-crafted features and a priori knowledge of the user. As such, these methods are difficult to be adapt to the clinical environment. Recently, deep learning methods based on fully convolutional networks (FCNs) have been successful in resolving segmentation problems. However, the reliance on large annotated training data is problematic when dealing with medical images. If a sufficient amount of annotated training data to cover all possible variations is not available, FCNs do not provide accurate segmentation. In addition, FCNs have a large receptive field in the convolutional layers, and hence produce coarse outputs of boundaries. Hence, we propose a new fully automated method that we refer to as a dual-stage fully convolutional networks (DSFCN). Our approach leverages deep residual architectures and FCNs and learns and infers the location of the optic cup and disk in a step-wise manner with fine-grained details. During training, our approach learns from the training data and the estimated results derived from the previous iteration. The ability to learn from the previous iteration optimizes the learning of the optic cup and the disk boundaries. During testing (prediction), DSFCN uses test (input) images and the estimated probability map derived from previous iterations to gradually improve the segmentation accuracy. Our method achieved an average Dice co-efficient of 0.8488 and 0.9441 for optic cup and disk segmentation and an area under curve (AUC) of 0.9513 for glaucoma detection.


title: Beyond the Chinese Restaurant and Pitman-Yor processes: Statistical Models with Double Power-law Behavior

abstract: Bayesian nonparametric approaches, in particular the Pitman-Yor process and the associated two-parameter Chinese Restaurant process, have been successfully used in applications where the data exhibit a power-law behavior. Examples include natural language processing, natural images or networks. There is also growing empirical evidence that some datasets exhibit a two-regime power-law behavior: one regime for small frequencies, and a second regime, with a different exponent, for high frequencies. In this paper, we introduce a class of completely random measures which are doubly regularly-varying. Contrary to the Pitman-Yor process, we show that when completely random measures in this class are normalized to obtain random probability measures and associated random partitions, such partitions exhibit a double power-law behavior. We discuss in particular three models within this class: the beta prime process (Broderick et al. (2015, 2018), a novel process called generalized BFRY process, and a mixture construction. We derive efficient Markov chain Monte Carlo algorithms to estimate the parameters of these models. Finally, we show that the proposed models provide a better fit than the Pitman-Yor process on various datasets.


title: Using discrete Darboux polynomials to detect and determine preserved measures and integrals of rational maps

abstract: In this Letter we propose a systematic approach for detecting and calculating preserved measures and integrals of a rational map. The approach is based on the use of cofactors and Discrete Darboux Polynomials and relies on the use of symbolic algebra tools. Given sufficient computing power, all rational preserved integrals can be found. We show, in two examples, how to use this method to detect and determine preserved measures and integrals of the considered rational maps.


title: Coherent Transition Radiation from Relativistic Beam-Foil Interaction in the Terahertz and Optical Range

abstract: Coherent transition radiation (CTR) from relativistic electron beam interaction with an overdense plasma foil is investigated by making use of two-dimensional particle-in-cell simulations. Well-defined single electron beam either of uniform profile or having substructures is considered for various beam-plasma parameters. The main purpose is to mimic the complicated beam-plasma conditions that is often found, for example, in intense laser plasma interactions. Key properties of the CTR concerning their temporal, angular and spectral profiles are identified. Several saturation effects due to the beam energy, size and foil density are found for the CTR energy, and the dependences vary for different spectral components such as in the Terahertz (THz) and optical range. The detailed substructure of the beam also affects greatly the radiation generation, leading to distinctive high harmonic components. Electrons with kinetic energy from sub MeV to tens of GeV are explored. For few MeV electron beams, the effects of the foil plasma on the beam dynamics and associated CTR generation, resembles closely the CTR from hot electrons produced in intense laser-plasma interactions. These results may find important applications in beam diagnostics either in laser-plasma based acceleration or conventional accelerators. They may also be employed to design novel THz radiation sources using tunable electron beams.


title: Anisotropic magnetic-field response of quantum critical fluctuations in Ni-doped CeCoIn5

abstract: This paper demonstrates the anisotropic response of quantum critical fluctuations with respect to the direction of the magnetic field $B$ in Ni-doped CeCoIn$_5$ by measuring the magnetization $M$ and specific heat $C$. The results show that $M/B$ at $B=0.1\ {\rm T}$ for both the tetragonal $c$ and $a$ directions exhibits $T^{-\eta}$ dependencies, and that $C/T$ at $B=0$ follows a $-\ln T$ function, which are the characteristics of non-Fermi-liquid (NFL) behaviors. For $B\,||\,c$, both the $M/B\propto T^{-\eta}$ and $C/T \propto -\ln T$ dependencies change into nearly temperature-constant behaviors by increasing $B$, indicating a crossover from the NFL state to the Fermi-liquid state. For $B\,||\,a$, however, the NFL behavior in $C/T$ persists up to $B=7\ {\rm T}$, whereas $M/B$ exhibits temperature-independent behavior for $B\ge 1\ {\rm T}$. These contrasting characteristics in $M/B$ and $C/T$ reflect the anisotropic nature of quantum critical fluctuations; the $c$-axis spin component significantly contributes to the quantum critical fluctuations. We compare this anisotropic behavior of the spin fluctuations to superconducting properties in pure CeCoIn$_5$, especially to the anisotropy in the upper critical field and the Ising-like characteristics in the spin resonance excitation, and suggest a close relationship between them.


title: Bell-Instability-Mediated Spectral Modulation of Hadronic Gamma Rays from a Supernova Remnant Interacting with a Molecular Cloud

abstract: Supernova remnants (SNRs) are believed to be the site of galactic cosmic-ray acceleration. However, the details of the cosmic-ray acceleration are still not well understood. Gamma ray observation is a promising method to study cosmic-ray acceleration in the SNRs, because a hadronic gamma ray can trace high energy cosmic-rays above ~GeV energy. Conventional theory predicts that the hadronic gamma ray shows a flat vFv spectrum from the pion-creation threshold energy to the maximum energy of diffusive shock acceleration. In this paper, by employing numerical simulations that solve a hybrid system of the magnetohydrodynamics of a molecular cloud and diffusive propagation of cosmic-rays, we demonstrate that the hadronic gamma ray spectrum can be harder than the conventional one and that the modulated spectrum becomes consistent with observations. The modification mechanism is explained as follows: The cosmic-rays accelerated at the supernova blast wave shock propagate into a clump of a molecular cloud. The cosmic-ray streaming in the cloud induces the so-called Bell instability that induces Alfven waves in the cloud. The induced magnetic field fluctuations prevent further cosmic-ray incursion by diminishing the diffusion coefficient for the cosmic-rays below ~1 TeV energy. This mechanism reinforces recent claims of a similar spectral modification by magnetic field amplification around a molecular cloud by Inoue et al.(2012) and Celli et al.(2018).


title: Magnon transport in quasi-two-dimensional van der Waals antiferromagnets

abstract: The recent emergence of 2D van der Waals magnets down to atomic layer thickness provides an exciting platform for exploring quantum magnetism and spintronics applications. The van der Waals nature stabilizes the long-range ferromagnetic order as a result of magnetic anisotropy. Furthermore, giant tunneling magnetoresistance and electrical control of magnetism have been reported. However, the potential of 2D van der Waals magnets for magnonics, magnon-based spintronics, has not been explored yet. Here, we report the experimental observation of long-distance magnon transport in quasi-twodimensional van der Waals antiferromagnet MnPS3, which demonstrates the 2D magnets as promising material candidates for magnonics. As the 2D MnPS3 thickness decreases, a shorter magnon diffusion length is observed, which could be attributed to the surface-impurity-induced magnon scattering. Our results could pave the way for exploring quantum magnonics phenomena and designing future magnonics devices based on 2D van der Waals magnets.


title: A Limited Habitable Zone for Complex Life

abstract: The habitable zone (HZ) is commonly defined as the range of distances from a host star within which liquid water, a key requirement for life, may exist on a planet's surface. Substantially more CO2 than present in Earth's modern atmosphere is required to maintain clement temperatures for most of the HZ, with several bars required at the outer edge. However, most complex aerobic life on Earth is limited by CO2 concentrations of just fractions of a bar. At the same time, most exoplanets in the traditional HZ reside in proximity to M dwarfs, which are more numerous than Sun-like G dwarfs but are predicted to promote greater abundances of gases that can be toxic in the atmospheres of orbiting planets, such as carbon monoxide (CO). Here we show that the HZ for complex aerobic life is likely limited relative to that for microbial life. We use a 1D radiative-convective climate and photochemical models to circumscribe a Habitable Zone for Complex Life (HZCL) based on known toxicity limits for a range of organisms as a proof of concept. We find that for CO2 tolerances of 0.01, 0.1, and 1 bar, the HZCL is only 21%, 32%, and 50% as wide as the conventional HZ for a Sun-like star, and that CO concentrations may limit some complex life throughout the entire HZ of the coolest M dwarfs. These results cast new light on the likely distribution of complex life in the universe and have important ramifications for the search for exoplanet biosignatures and technosignatures.


title: Energy Efficiency of RSMA and NOMA in Cellular-Connected mmWave UAV Networks

abstract: Cellular-connected unmanned aerial vehicles (UAVs) are recently getting significant attention due to various practical use cases, e.g., surveillance, data gathering, purchase delivery, among other applications. Since UAVs are low power nodes, energy and spectral efficient communication is of paramount importance. To that end, multiple access (MA) schemes can play an important role in achieving high energy efficiency and spectral efficiency. In this work, we introduce rate-splitting MA (RSMA) and non-orthogonal MA (NOMA) schemes in a cellular-connected UAV network. In particular, we investigate the energy efficiency of the RSMA and NOMA schemes in a millimeter wave (mmWave) downlink transmission scenario. Furthermore, we optimize precoding vectors of both the schemes by explicitly taking into account the 3GPP antenna propagation patterns. The numerical results for this realistic transmission scheme indicate that RSMA is superior to NOMA in terms of overall energy efficiency.


title: Technical Report: All Principal Congruence Link Groups

abstract: This is a technical report accompanying the paper "All Principal Congruence Link Groups" (arXiv:1802.01275) classifying all principal congruence link complements in S^3 by the same authors. It provides a complete overview of all cases (d,I) that had to be considered, as well as describes the necessary computations and computer programs written for the classification result.


title: The matroid structure of vectors of the Mordell-Weil lattice and the topology of plane quartics and bitangent lines

abstract: In this paper, we introduce the terminology of matroids into the study of Zariski-pairs related to rational elliptic surfaces, aiming to simplify the presentation and arguments involved. As an application, we provide new examples of Zariski $N$-ples of relatively low degree. Namely we show that a Zariski 102-ple of degree 18 exists.


title: A characterisation of S-box fitness landscapes in cryptography

abstract: Substitution Boxes (S-boxes) are nonlinear objects often used in the design of cryptographic algorithms. The design of high quality S-boxes is an interesting problem that attracts a lot of attention. Many attempts have been made in recent years to use heuristics to design S-boxes, but the results were often far from the previously known best obtained ones. Unfortunately, most of the effort went into exploring different algorithms and fitness functions while little attention has been given to the understanding why this problem is so difficult for heuristics. In this paper, we conduct a fitness landscape analysis to better understand why this problem can be difficult. Among other, we find that almost each initial starting point has its own local optimum, even though the networks are highly interconnected.


title: Displacement of transport processes on networked topologies

abstract: Consider a particle whose position evolves along the edges of a network. One definition for the displacement of a particle is the length of the shortest path on the network between the current and initial positions of the particle. Such a definition fails to incorporate information of the actual path the particle traversed. In this work we consider another definition for the displacement of a particle on networked topologies. Using this definition, which we term the winding distance, we demonstrate that for Brownian particles, confinement to a network can induce a transition in the mean squared displacement from diffusive to ballistic behaviour, $\langle x^2(t) \rangle \propto t^2$ for long times. A multiple scales approach is used to derive a macroscopic evolution equation for the displacement of a particle and uncover a topological condition for whether this transition in the mean squared displacement will occur. Furthermore, for networks satisfying this topological condition, we identify a prediction of the timescale upon which the displacement transitions to long-time behaviour. Finally, we extend the investigation of displacement on networks to a class of anomalously diffusive transport processes, where we find that the mean squared displacement at long times is affected by both network topology and the character of the transport process.


title: Correlated emission of X-ray and sound from water film irradiated by femtosecond laser pulses

abstract: Simultaneous measurements of hard X-ray by a Geiger counter and audible sound (10 Hz-20kHz) by a microphone from a thin water film in air were carried out under intense single and double pulse irradiations of femtosecond laser (35 fs, 800 nm, 1 kHz). Emission profiles of X-ray and sound under the single pulse irradiation by changing the water film position along the laser incident direction (Z-axis) show the same peak positions with a broader emission in sound (403{\mu}m at FWHM) than in X-ray (37{\mu}m). Under the double pulse irradiation condition with the time delay at 0 ps and 4.6 ns, it was clearly observed that the acoustic signal intensity is enhanced in associated with X-ray intensity enhancements. The enhancements can be assigned to laser ablation dynamics such as pre-plasma formation and transient surface roughness formation induced by the pre-pulse irradiation. For the acoustic signal under the double-pulse irradiation with the time delay, there was a weak dependence observed on the pre-pulse irradiation position at the laser focus. It is consistent with a long breakdown filament formation which makes the microphone-detection less position-sensitive.


title: A simple statistical approach to prediction in open high dimensional chaotic systems

abstract: Two recent papers on prediction of chaotic systems, one on multi-view embedding1 , and the second on prediction in projection2 provide empirical evidence to support particular prediction methods for chaotic systems. Multi-view embedding1 is a method of using several multivariate time series to come up with an improved embedding based predictor of a chaotic time series. Prediction in projection2 discusses how much smaller embeddings can provide useful prediction even though they may not be able to resolve the dynamics of the system. Both papers invoke a nearest neighbor3, or Lorenz method of Analogue (LMA)4 approach to estimation. However with open high dimensional chaotic systems there may be no very close nearest neighbor trajectories in a history, so in this paper we add in the thread of linear response theory5, although our approach is quite simple assuming that linear regressions6 on multiple embeddings5 will be sufficient. The approach is thus to create linear response models5 of multiple1 low dimensional embeddings2 to provide practical methods of predicting in high dimensional, open chaotic systems. Some theory is developed, mostly around two unproven conjectures to suggest methods of prediction, and they are applied to prediction in the earth's climate system (predicting regional rainfall in a small region) multiple seasons ahead and to power a simulated automated trading program applied to three stock indexes, the Dow Jones industrial average, the Dow Jones transportation average, and the Dow Jones utility average, based on data from 1929 through 20077


title: Learning Ising Models with Independent Failures

abstract: We give the first efficient algorithm for learning the structure of an Ising model that tolerates independent failures; that is, each entry of the observed sample is missing with some unknown probability p. Our algorithm matches the essentially optimal runtime and sample complexity bounds of recent work for learning Ising models due to Klivans and Meka (2017). We devise a novel unbiased estimator for the gradient of the Interaction Screening Objective (ISO) due to Vuffray et al. (2016) and apply a stochastic multiplicative gradient descent algorithm to minimize this objective. Solutions to this minimization recover the neighborhood information of the underlying Ising model on a node by node basis.


title: Accurate 3D Cell Segmentation using Deep Feature and CRF Refinement

abstract: We consider the problem of accurately identifying cell boundaries and labeling individual cells in confocal microscopy images, specifically, 3D image stacks of cells with tagged cell membranes. Precise identification of cell boundaries, their shapes, and quantifying inter-cellular space leads to a better understanding of cell morphogenesis. Towards this, we outline a cell segmentation method that uses a deep neural network architecture to extract a confidence map of cell boundaries, followed by a 3D watershed algorithm and a final refinement using a conditional random field. In addition to improving the accuracy of segmentation compared to other state-of-the-art methods, the proposed approach also generalizes well to different datasets without the need to retrain the network for each dataset. Detailed experimental results are provided, and the source code is available on GitHub.


title: Log-decay $F$-isocrystals on higher dimensional varieties

abstract: Let $k$ be a perfect field of positive characteristic and let $X$ be a smooth irreducible quasi-compact scheme over $k$. The Drinfeld-Kedlaya theorem states that for an irreducible $F$-isocrystal on $X$, the gap between consecutive generic slopes is bounded by one. In this note we provide a new proof of this theorem. Our proof utilizes the theory of $F$-isocrystals with $r$-log decay. We first show that a rank one $F$-isocrystal with $r$-log decay is overconvergent if $r<1$. Next, we establish a connection between slope gaps and the rate of log-decay of the slope filtration. The Drinfeld-Kedlaya theorem then follows from a simple patching argument.


title: Jointly Low-Rank and Bisparse Recovery: Questions and Partial Answers

abstract: We investigate the problem of recovering jointly $r$-rank and $s$-bisparse matrices from as few linear measurements as possible, considering arbitrary measurements as well as rank-one measurements. In both cases, we show that $m \asymp r s \ln(en/s)$ measurements make the recovery possible in theory, meaning via a nonpractical algorithm. In case of arbitrary measurements, we investigate the possibility of achieving practical recovery via an iterative-hard-thresholding algorithm when $m \asymp r s^\gamma \ln(en/s)$ for some exponent $\gamma > 0$. We show that this is feasible for $\gamma = 2$, and that the proposed analysis cannot cover the case $\gamma \leq 1$. The precise value of the optimal exponent $\gamma \in [1,2]$ is the object of a question, raised but unresolved in this paper, about head projections for the jointly low-rank and bisparse structure. Some related questions are partially answered in passing. For rank-one measurements, we suggest on arcane grounds an iterative-hard-thresholding algorithm modified to exploit the nonstandard restricted isometry property obeyed by this type of measurements.


title: Statistical Failure Mechanism Analysis of Earthquakes Revealing Time Relationships

abstract: If we assume that earthquakes are chaotic, and influenced locally then chaos theory suggests that there should be a temporal association between earthquakes in a local region that should be revealed with statistical examination. To date no strong relationship has been shown (refs not prediction). However, earthquakes are basically failures of structured material systems, and when multiple failure mechanisms are present, prediction of failure is strongly inhibited without first separating the mechanisms. Here we show that by separating earthquakes statistically, based on their central tensor moment structure, along lines first suggested by a separation into mechanisms according to depth of the earthquake, a strong indication of temporal association appears. We show this in earthquakes above 200 Km along the pacific ring of fire, with a positive association in time between earthquakes of the same statistical type and a negative association in time between earthquakes of different types. Whether this can reveal either useful mechanistic information to seismologists, or can result in useful forecasts remains to be seen.


title: Learning partial differential equations for biological transport models from noisy spatiotemporal data

abstract: We investigate methods for learning partial differential equation (PDE) models from spatiotemporal data under biologically realistic levels and forms of noise. Recent progress in learning PDEs from data have used sparse regression to select candidate terms from a denoised set of data, including approximated partial derivatives. We analyze the performance in utilizing previous methods to denoise data for the task of discovering the governing system of partial differential equations (PDEs). We also develop a novel methodology that uses artificial neural networks (ANNs) to denoise data and approximate partial derivatives. We test the methodology on three PDE models for biological transport, i.e., the advection-diffusion, classical Fisher-KPP, and nonlinear Fisher-KPP equations. We show that the ANN methodology outperforms previous denoising methods, including finite differences and polynomial regression splines, in the ability to accurately approximate partial derivatives and learn the correct PDE model.


title: Circuit quantization in the presence of time-dependent external flux

abstract: Circuit quantization links a physical circuit to its corresponding quantum Hamiltonian. The standard quantization procedure generally assumes any external magnetic flux to be static. Time dependence naturally arises, however, when flux is modulated or when flux noise is considered. In this case, application of the existing quantization procedure can lead to inconsistencies. To resolve these, we generalize circuit quantization to incorporate time-dependent external flux.


title: Computing the Yolk in Spatial Voting Games without Computing Median Lines

abstract: The yolk is an important concept in spatial voting games as it generalises the equilibrium and provides bounds on the uncovered set. We present near-linear time algorithms for computing the yolk in the spatial voting model in the plane. To the best of our knowledge our algorithm is the first algorithm that does not require precomputing the median lines and hence able to break the existing $O(n^{4/3})$ bound which equals the known upper bound on the number of median lines. We avoid this requirement by using Megiddo's parametric search, which is a powerful framework that could lead to faster algorithms for many other spatial voting problems.


title: Sufficient condition for existence of special type of primitive normal elements over finite fields

abstract: Let $\mathbb{F}_{q^n}$ be the extension of the field $\mathbb{F}_q$ of degree n, where $q$ is power of prime $p$, i.e $q=p^k$, where k is a positive integer. In this paper, we provide sufficient condition for the existence of a primitive normal element $\alpha\in\mathbb{F}_{q^n} $ such that $\alpha^2+\alpha+1$ is also primitive normal element over $\mathbb{F}_{q^n}$.


title: Minimalist large scale dynamo from shear-driven inhomogeneity

abstract: We show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $\alpha$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism underlying large scale field generation in astrophysical, laboratory, and simulated shear flows, and highlights the essentiality of correctly averaging dynamo coefficients and fields.

categories: physics.flu-dyn astro-ph.HE astro-ph.SR physics.plasm-ph


title: Mesh: Compacting Memory Management for C/C++ Applications

abstract: Programs written in C/C++ can suffer from serious memory fragmentation, leading to low utilization of memory, degraded performance, and application failure due to memory exhaustion. This paper introduces Mesh, a plug-in replacement for malloc that, for the first time, eliminates fragmentation in unmodified C/C++ applications. Mesh combines novel randomized algorithms with widely-supported virtual memory operations to provably reduce fragmentation, breaking the classical Robson bounds with high probability. Mesh generally matches the runtime performance of state-of-the-art memory allocators while reducing memory consumption; in particular, it reduces the memory of consumption of Firefox by 16% and Redis by 39%.


title: Global existence and blowup for Choquard equations with an inverse-square potential

abstract: In this paper, the Choquard equation with an inverse-square potential and both focusing and defocusing nonlinearities in the energy-subcritical regime is investigated. For all the cases, the local well-posedness result in $H^1(\mathbb{R}^N)$ is established. Moreover, the global existence result for arbitrary initial values is proved in the defocusing case while a global existence/blowup dichotomy below the ground state is established in the focusing case.


title: CSPs with Global Modular Constraints: Algorithms and Hardness via Polynomial Representations

abstract: We study the complexity of Boolean constraint satisfaction problems (CSPs) when the assignment must have Hamming weight in some congruence class modulo M, for various choices of the modulus M. Due to the known classification of tractable Boolean CSPs, this mainly reduces to the study of three cases: 2-SAT, HORN-SAT, and LIN-2 (linear equations mod 2). We classify the moduli M for which these respective problems are polynomial time solvable, and when they are not (assuming the ETH). Our study reveals that this modular constraint lends a surprising richness to these classic, well-studied problems, with interesting broader connections to complexity theory and coding theory. The HORN-SAT case is connected to the covering complexity of polynomials representing the NAND function mod M. The LIN-2 case is tied to the sparsity of polynomials representing the OR function mod M, which in turn has connections to modular weight distribution properties of linear codes and locally decodable codes. In both cases, the analysis of our algorithm as well as the hardness reduction rely on these polynomial representations, highlighting an interesting algebraic common ground between hard cases for our algorithms and the gadgets which show hardness. These new complexity measures of polynomial representations merit further study. The inspiration for our study comes from a recent work by N\"agele, Sudakov, and Zenklusen on submodular minimization with a global congruence constraint. Our algorithm for HORN-SAT has strong similarities to their algorithm, and in particular identical kind of set systems arise in both cases. Our connection to polynomial representations leads to a simpler analysis of such set systems, and also sheds light on (but does not resolve) the complexity of submodular minimization with a congruency requirement modulo a composite M.


title: Learning to Screen

abstract: Imagine a large firm with multiple departments that plans a large recruitment. Candidates arrive one-by-one, and for each candidate the firm decides, based on her data (CV, skills, experience, etc), whether to summon her for an interview. The firm wants to recruit the best candidates while minimizing the number of interviews. We model such scenarios as an assignment problem between items (candidates) and categories (departments): the items arrive one-by-one in an online manner, and upon processing each item the algorithm decides, based on its value and the categories it can be matched with, whether to retain or discard it (this decision is irrevocable). The goal is to retain as few items as possible while guaranteeing that the set of retained items contains an optimal matching. We consider two variants of this problem: (i) in the first variant it is assumed that the $n$ items are drawn independently from an unknown distribution $D$. (ii) In the second variant it is assumed that before the process starts, the algorithm has an access to a training set of $n$ items drawn independently from the same unknown distribution (e.g.\ data of candidates from previous recruitment seasons). We give tight bounds on the minimum possible number of retained items in each of these variants. These results demonstrate that one can retain exponentially less items in the second variant (with the training set).


title: Uniform convergence may be unable to explain generalization in deep learning

abstract: Aimed at explaining the surprisingly good generalization behavior of overparameterized deep networks, recent works have developed a variety of generalization bounds for deep learning, all based on the fundamental learning-theoretic technique of uniform convergence. While it is well-known that many of these existing bounds are numerically large, through numerous experiments, we bring to light a more concerning aspect of these bounds: in practice, these bounds can {\em increase} with the training dataset size. Guided by our observations, we then present examples of overparameterized linear classifiers and neural networks trained by gradient descent (GD) where uniform convergence provably cannot "explain generalization" -- even if we take into account the implicit bias of GD {\em to the fullest extent possible}. More precisely, even if we consider only the set of classifiers output by GD, which have test errors less than some small $\epsilon$ in our settings, we show that applying (two-sided) uniform convergence on this set of classifiers will yield only a vacuous generalization guarantee larger than $1-\epsilon$. Through these findings, we cast doubt on the power of uniform convergence-based generalization bounds to provide a complete picture of why overparameterized deep networks generalize well.


title: Session-based Sequential Skip Prediction via Recurrent Neural Networks

abstract: The focus of WSDM cup 2019 is session-based sequential skip prediction, i.e. predicting whether users will skip tracks, given their immediately preceding interactions in their listening session. This paper provides the solution of our team \textbf{ekffar} to this challenge. We focus on recurrent-neural-network-based deep learning approaches which have previously been shown to perform well on session-based recommendation problems. We show that by choosing an appropriate recurrent architecture that properly accounts for the given information such as user interaction features and song metadata, a single neural network could achieve a Mean Average Accuracy (AA) score of 0.648 on the withheld test data. Meanwhile, by ensembling several variants of the core model, the overall recommendation accuracy can be improved even further. By using the proposed approach, our team was able to attain the 1st place in the competition. We have open-sourced our implementation at GitHub.


title: Proving Expected Sensitivity of Probabilistic Programs with Randomized Variable-Dependent Termination Time

abstract: The notion of program sensitivity (aka Lipschitz continuity) specifies that changes in the program input result in proportional changes to the program output. For probabilistic programs the notion is naturally extended to expected sensitivity. A previous approach develops a relational program logic framework for proving expected sensitivity of probabilistic while loops, where the number of iterations is fixed and bounded. In this work, we consider probabilistic while loops where the number of iterations is not fixed, but randomized and depends on the initial input values. We present a sound approach for proving expected sensitivity of such programs. Our sound approach is martingale-based and can be automated through existing martingale-synthesis algorithms. Furthermore, our approach is compositional for sequential composition of while loops under a mild side condition. We demonstrate the effectiveness of our approach on several classical examples from Gambler's Ruin, stochastic hybrid systems and stochastic gradient descent. We also present experimental results showing that our automated approach can handle various probabilistic programs in the literature.


title: A Strongly Miniaturized and Inherently Matched Folded Dipole Antenna for Narrowband Applications

abstract: A novel miniaturized printed folded dipole antenna has been designed for narrowband sensing applications. It is shown that the antenna may be systematically engineered to achieve matching to a practical source impedance at frequencies well below its half-wavelength resonance using series-LC loading and asymmetrical arm widths. This approach is used to design an inherently well-matched antenna that does not require an external matching network and demonstrates approximately 75% miniaturization, excellent co- to cross-polarization separation, over 24 dB higher realized gain, and strongly improved radiation efficiency relative to an unloaded folded dipole of the same electrically small size. Experimental validation of a fabricated prototype demonstrates excellent agreement with simulations.


title: An analytic approach to the Riemann hypothesis

abstract: In this work we consider an equation for the Riemann zeta-function in the critical half-strip. With the help of this equation we prove that finding non-trivial zeros of the Riemann zeta-function outside the critical line would be equivalent to the existence of complex numbers for which equation (5.1) in the paper holds. Such a condition is studied, and the attempt of proving the Riemann hypothesis is found to involve also the functional equation (6.26), where t is a real variable bigger than or equal to 1 and n is any natural number. The limiting behavior of the solutions as t approaches 1 is then studied in detail.


title: Unidirectional light transmission by two-layer nanostructures interacting via optical near-fields

abstract: We experimentally demonstrate unidirectional light transmission through two-layer nanostructured materials considering that the horizontal-to-vertical-polarization conversion efficiency in the forward direction and the vertical-to-horizontal efficiency in the backward direction, which are usually equivalent due to optical reciprocity, are different. The different ways of transferring light momentum in the forward and backward directions via optical near-fields between the layers are responsible for the unidirectionality of light, which was theoretically investigated in our recent work [J. Opt. Soc. Am. B 31, 2404-2413]. With two-layer metal nanostructures experimentally fabricated via a repeated lift-off technique, consistent optical characteristics are observed, verifying the utilization of the large momentum of optical near-fields.


title: Amphoteric behavior of hydrogen (H+1and H-1) in complex hydrides from van der Waals interaction included ab-initio calculation

abstract: In order to identify potential hydrogen storage materials, worldwide attention has been focused on hydrides with high gravimetric and volumetric capacity. Hydrogen is a unique element that possess positive, negative or neutral oxidation state in solids depending upon the chemical environment. If one can find hydrogen storage materials where hydrogen is present in both negative and positive oxidation state within the same structural framework then one can accommodate hydrogen with high volume density. So, it is fundamentally as well as technologically important to identify compounds in which hydrogen is in amphoteric nature and understand the necessary criteria for its origin. The experimental structural analysis of Cyclotriborazane and Diammonium dodeca hydro-closo-dodecaborate insinuate the presence of hydrogen with anionic and cationic behavior within the same structure. In order to understand the role of van der Waals (vdW) interactions on structural parameters, we have considered 12 different vdW corrected functionals and found that the optPBE-vdW functional predicts the equilibrium structural parameters reliably with less than 0.009 % accuracy. So, the optPBE-vdW functional is used to calculate charge density, electron localization function, total as well as the partial density of state, Bader and Born effective charge etc. We have observed that hydrogen is exhibiting amphoteric behavior with H closer to B in the negatively charged state and that neighboring to N is in a positively charged state.


title: Dependence Properties of B-Spline Copulas

abstract: We construct by using B-spline functions a class of copulas that includes the Bernstein copulas arising in Baker's distributions. The range of correlation of the B-spline copulas is examined, and the Frechet--Hoeffding upper bound is proved to be attained when the number of B-spline functions goes to infinity. As the B-spline functions are well-known to be an order-complete weak Tchebycheff system from which the property of total positivity of any order follows for the maximum correlation case, the results given here extend classical results for the Bernstein copulas. In addition, we derive in terms of the Stirling numbers of the second kind an explicit formula for the moments of the related B-spline functions for nonnegative real numbers.


title: I-mode investigation on the Experimental Advanced Superconducting Tokamak

abstract: By analyzing large quantities of discharges in the unfavorable ion $ \vec B\times \nabla B $ drift direction, the I-mode operation has been confirmed in EAST tokamak. During the L-mode to I-mode transition, the energy confinement has a prominent improvement by the formation of a high-temperature edge pedestal, while the particle confinement remains almost identical to that in the L-mode. Similar with the I-mode observation on other devices, the $ E_r $ profiles obtained by the eight-channel Doppler backscattering system (DBS8)\cite{J.Q.Hu} show a deeper edge $ E_r $ well in the I-mode than that in the L-mode. And a weak coherent mode (WCM) with the frequency range of 40-150 kHz is observed at the edge plasma with the radial extend of about 2-3 cm. WCM could be observed in both density fluctuation and radial electric field fluctuation, and the bicoherence analyses showed significant couplings between WCM and high frequency turbulence, implying that the $ E_r $ fluctuation and the caused flow shear from WCM should play an important role during I-mode. In addition, a low-frequency oscillation with a frequency range of 5-10 kHz is always accompanied with WCM, where GAM intensity is decreased or disappeared. Many evidences show that the a low-frequency oscillation may be a novel kind of limited cycle oscillation but further investigations are needed to explain the new properties such as the harmonics and obvious magnetical perturbations.


title: Random positive operator valued measures

abstract: We introduce several notions of random positive operator valued measures (POVMs), and we prove that some of them are equivalent. We then study statistical properties of the effect operators for the canonical examples, obtaining limiting eigenvalue distributions with the help of free probability theory. Similarly, we obtain the large system limit for several quantities of interest in quantum information theory, such as the sharpness, the noise content, and the probability range. Finally, we study different compatibility criteria, and we compare them for generic POVMs.