title: On Extreme Value Index Estimation under Random Censoring

abstract: Extreme value analysis in the presence of censoring is receiving much attention as it has applications in many disciplines, including survival and reliability studies. Estimation of extreme value index (EVI) is of primary importance as it is a critical parameter needed in estimating extreme events such as quantiles and exceedance probabilities. In this paper, we review several estimators of the extreme value index when data is subject to random censoring. In addition, four estimators are proposed, one based on the exponential regression approximation of log spacings, one based on a Zipf estimator and two based on variants of the moment estimator. The proposed estimators and the existing ones are compared under the same simulation conditions. The performance measures for the estimators include confidence interval length and coverage probability. The simulation results show that no estimator is universally the best as the estimators depend on the size of the EVI parameter, percentage of censoring in the right tail and the underlying distribution. However, certain estimators such as the proposed reduced-bias estimator and the adapted moment estimator are found to perform well across most scenarios. Moreover, we present a bootstrap algorithm for obtaining samples for extreme value analysis in the context of censoring. Some of the estimators that performed well in the simulation study are illustrated using a practical dataset from medical research


title: Enhancement of Coulomb blockade in epsilon near zero and hyperbolic metamaterials

abstract: Coulomb interaction of charges in a metamaterial may be expressed via its effective dielectric response function. Coulomb interaction appears considerably enhanced in artificially engineered epsilon near zero (ENZ) and hyperbolic metamaterials. Indeed, very recently it was observed that negative ENZ conditions enhance superconducting properties of a composite metamaterial via enhancement of Cooper pairing. Here we demonstrate that positive ENZ conditions enhance Coulomb repulsion of charges, leading to considerable enhancement of the Coulomb blockade effect in such metamaterials. Similar to high Tc cuprates, in hyperbolic metamaterials superconductivity and the Coulomb blockade-induced charge ordering is predicted to compete with each other.


title: Signatures of thermal hysteresis in Tamm-wave propagation

abstract: We numerically solved the boundary-value problem for Tamm waves (which may also be classified as Uller-Zenneck waves here) guided by the planar interface of a homogeneous isotropic dissipative dielectric (HIDD) material and a periodically multilayered isotropic dielectric material. The HIDD material was chosen to be VO${}_2$ which, at optical wavelengths, has a temperature-dependent refractive index with a hysteresis feature, i.e., the temperature-dependence of its refractive index varies depending upon whether the temperature is increasing or decreasing. A numerical code was implemented to extract solutions of the dispersion equation at a fixed wavelength for both $p$- and $s$-polarization states over the temperature range [50,80] degrees. A multitude of Tamm waves of both linear polarization states were found, demonstrating a clear demarcation of the heating and cooling phases in terms of wavenumbers and propagation distances. Thereby, the signatures of thermal hysteresis in Tamm-wave propagation were revealed.


title: A Simulation Comparison of Estimators of Conditional Extreme Value Index under Right Random Censoring

abstract: In extreme value analysis, the extreme value index plays a vital role as it determines the tail heaviness of the underlying distribution and is the primary parameter required for the estimation of other extreme events. In this paper, we review the estimation of the extreme value index when observations are subject to right random censoring and the presence of covariate information. In addition, we propose some estimators of the extreme value index, including a maximum likelihood estimator from a perturbed Pareto distribution. The existing estimators and the proposed ones are compared through a simulation study under identical conditions. The results show that the performance of the estimators depend on the percentage of censoring, the underlying distribution, the size of extreme value index and the number of top order statistics. Overall, we found the proposed estimator from the perturbed Pareto distribution to be robust to censoring, size of the extreme value index and the number of top order statistics.


title: Mappings of finite distortion: compactness of the branch set

abstract: We show that an entire branched cover of finite distortion cannot have a compact branch set if its distortion satisfies a certain asymptotic growth condition. We furthermore show that this bound is strict by constructing an entire, continuous, open and discrete mapping of finite distortion which is piecewise smooth, has a branch set homeomorphic to $(n-2)$-dimensional torus and distortion arbitrarily close to the asymptotic bound.


title: A Machine Learning Framework to Forecast Wave Conditions

abstract: A~machine learning framework is developed to estimate ocean-wave conditions. By supervised training of machine learning models on many thousands of iterations of a physics-based wave model, accurate representations of significant wave heights and period can be used to predict ocean conditions. A model of Monterey Bay was used as the example test site; it was forced by measured wave conditions, ocean-current nowcasts, and reported winds. These input data along with model outputs of spatially variable wave heights and characteristic period were aggregated into supervised learning training and test data sets, which were supplied to machine learning models. These machine learning models replicated wave heights with a root-mean-squared error of 9cm and correctly identify over 90% of the characteristic periods for the test-data sets. Impressively, transforming model inputs to outputs through matrix operations requires only a fraction (<1/1,000) of the computation time compared to forecasting with the physics-based model.


title: Standard Temperature and Pressure Superconductivity

abstract: Superconductivity at standard temperature and pressure is far from the extreme conditions where new fundamental laws of physics are expected to arise. Yet it is widely believed that the Landau-Ginzburg-Wilson-Fisher paradigm of broken symmetry and renormalization does not give a satisfactory account of the phenomenon. Almost a decade ago, we used the Bardeen-Cooper-Shrieffer wavefunction to show that superconductors have topological order. Here we report progress using topological order to look for standard temperature and pressure superconductivity.


title: Wolf-Rayet stars in the Small Magellanic Cloud as testbed for massive star evolution

abstract: The majority of Wolf-Rayet (WR) stars represent the stripped cores of evolved massive stars who lost most of their hydrogen envelope. In low metallicity environments, such as the Small Magellanic Cloud (SMC), stellar winds are weaker and binary interaction is expected to dominate WR-star formation. However, the WR binary fraction appears to be ~40% at any metallicity. We use the recently determined physical properties of the twelve known SMC WR stars to explore their possible formation channels through comparisons with grids of SMC models, simulated with the detailed stellar evolution code MESA. These include models of rapidly rotating single stars, which experience (partial) chemically homogeneous evolution (CHE). We find that CHE is not able to account for the majority of the SMC WR stars. However, the apparently single WN star SMC AB12 and the double WR system SMC AB5 (HD 5980) appear consistent with this channel. We also analyze core helium burning stellar models assuming constant hydrogen gradients in their envelopes. We find a dichotomy in the envelope hydrogen gradients required to explain the observed temperatures of the SMC WR stars. Shallow gradients are found for the WR stars with O star companions, consistent with binary models where mass transfer occurs early, which is in agreement with their binary properties. On the other hand, much steeper hydrogen gradients are inferred for the group of hot apparently single WR stars. Since the hydrogen profiles in post main sequence models of massive stars become steeper with time, we conclude that these stars have likely been stripped by a companion during a phase of common envelope evolution. The companions, either main sequence stars or compact objects, are expected to still be present. A corresponding search might identify the first immediate double black hole binary progenitor with masses as high as those detected in GW150914.


title: Stochastic Nonconvex Optimization with Large Minibatches

abstract: We study stochastic optimization of nonconvex loss functions, which are typical objectives for training neural networks. We propose stochastic approximation algorithms which optimize a series of regularized, nonlinearized losses on large minibatches of samples, using only first-order gradient information. Our algorithms provably converge to an approximate critical point of the expected objective with faster rates than minibatch stochastic gradient descent, and facilitate better parallelization by allowing larger minibatches.


title: Fermi--Pasta--Ulam--Tsingou problems: Passage from Boltzmann to $q$-statistics

abstract: The Fermi-Pasta-Ulam (FPU) one-dimensional Hamiltonian includes a quartic term which guarantees ergodicity of the system in the thermodynamic limit. Consistently, the Boltzmann factor $P(\epsilon) \sim e^{-\beta \epsilon}$ describes its equilibrium distribution of one-body energies, and its velocity distribution is Maxwellian, i.e., $P(v) \sim e^{- \beta v^2/2}$. We consider here a generalized system where the quartic coupling constant between sites decays as $1/d_{ij}^{\alpha}$ $(\alpha \ge 0; d_{ij} = 1,2,\dots)$. Through {\it first-principle} molecular dynamics we demonstrate that, for large $\alpha$ (above $\alpha \simeq 1$), i.e., short-range interactions, Boltzmann statistics (based on the {\it additive} entropic functional $S_B[P(z)]=-k \int dz P(z) \ln P(z)$) is verified. However, for small values of $\alpha$ (below $\alpha \simeq 1$), i.e., long-range interactions, Boltzmann statistics dramatically fails and is replaced by q-statistics (based on the {\it nonadditive} entropic functional $S_q[P(z)]=k (1-\int dz [P(z)]^q)/(q-1)$, with $S_1 = S_B$). Indeed, the one-body energy distribution is q-exponential, $P(\epsilon) \sim e_{q_{\epsilon}}^{-\beta_{\epsilon} \epsilon} \equiv [1+(q_{\epsilon} - 1) \beta_{\epsilon}{\epsilon}]^{-1/(q_{\epsilon}-1)}$ with $q_{\epsilon} > 1$, and its velocity distribution is given by $P(v) \sim e_{q_v}^{ - \beta_v v^2/2}$ with $q_v > 1$. Moreover, within small error bars, we verify $q_{\epsilon} = q_v = q$, which decreases from an extrapolated value q $\simeq$ 5/3 to q=1 when $\alpha$ increases from zero to $\alpha \simeq 1$, and remains q = 1 thereafter.


title: Understanding a Version of Multivariate Symmetric Uncertainty to assist in Feature Selection

abstract: In this paper, we analyze the behavior of the multivariate symmetric uncertainty (MSU) measure through the use of statistical simulation techniques under various mixes of informative and non-informative randomly generated features. Experiments show how the number of attributes, their cardinalities, and the sample size affect the MSU. We discovered a condition that preserves good quality in the MSU under different combinations of these three factors, providing a new useful criterion to help drive the process of dimension reduction.


title: Computing properties of stable configurations of thermodynamic binding networks

abstract: The promise of chemical computation lies in controlling systems incompatible with traditional electronic micro-controllers, with applications in synthetic biology and nano-scale manufacturing. Computation is typically embedded in kinetics---the specific time evolution of a chemical system. However, if the desired output is not thermodynamically stable, basic physical chemistry dictates that thermodynamic forces will drive the system toward error throughout the computation. The thermodynamic binding network (TBN) model was introduced to formally study how the thermodynamic equilibrium can be made consistent with the desired computation, and it idealizes tradeoffs between configurational entropy and binding. Here we prove the computational hardness of natural questions about TBNs and develop a practical algorithm for verifying the correctness of constructions by translating the problem into propositional logic and solving the resulting formula. The TBN model together with automated verification tools will help inform strategies for error reduction in molecular computing, including the extensively studied models of strand displacement cascades and algorithmic tile assembly.


title: Non-Abelian cosmic string in the Starobinsky model of gravity

abstract: In this paper, we analyze numerically the behaviour of the solutions corresponding to a non-Abelian cosmic string in the framework of the Starobinsky model, i.e. where $f(R)=R + \zeta R^2$. We perform the calculations for both an asymptotically flat and asymptotically (anti)de Sitter spacetimes. We found that the angular deficit generated by the string decreases as the parameter $\zeta$ increases, in the case of a null cosmological constant. For a positive cosmological constant, we found that the cosmic horizon is affected in a non-trivial way by the parameter $\zeta$.


title: In search of cool flow accretion onto galaxies $-$ where does the disk gas end?

abstract: The processes taking place in the outermost reaches of spiral disks (the 'proto-disk') are intimately connected to the build-up of mass and angular momentum in galaxies. The thinness of spiral disks suggests that the activity is mostly quiescent and presumably this region is fed by cool flows coming into the halo from the intergalactic medium. While there is abundant evidence for the presence of a circumgalactic medium (CGM) around disk galaxies as traced by quasar absorption lines, it has been very difficult to connect this material to the outer gas disk. This has been a very difficult transition region to explore because baryon tracers are hard to observe. In particular, HI disks have been argued to truncate at a critical column density N(H) $\approx 3\times 10^{19}$ cm$^{-2}$ at 30 kpc for an L* galaxy where the gas is vulnerable to the external ionizing background. But new deep observations of nearby L* spirals (e.g. Milky Way, NGC 2997) suggest that HI disks may extend much further than recognised to date, up to 60 kpc at N(H) $\approx 10^{18}$ cm$^{-2}$. Motivated by these observations, here we show that a clumpy outer disk of dense clouds or cloudlets is potentially detectable to much larger radii and lower HI column densities than previously discussed. This extended proto-disk component is likely to explain some of the MgII forest seen in quasar spectra as judged from absorption-line column densities and kinematics. We fully anticipate that the armada of new radio facilities and planned HI surveys coming online will detect this extreme outer disk (scree) material. We also propose a variant on the successful 'Dragonfly' technique to go after the very weak H$\alpha$ signals expected in the proto-disk region.


title: Spaces which invert weak homotopy equivalences

abstract: It is well known that if $X$ is a CW-complex, then for every weak homotopy equivalence $f:A\to B$, the map $f_*:[X,A]\to [X,B]$ induced in homotopy classes is a bijection. For which spaces $X$ is $f^*:[B,X]\to [A,X]$ a bijection for every weak equivalence $f$? This question was considered by J. Strom and T. Goodwillie. In this note we prove that a non-empty space inverts weak equivalences if and only if it is contractible.


title: On $^{146}$Nd, $^{144}$Sm and other unexplored 2$\beta$ decay isotopes

abstract: $^{146}$Nd is one of only four (over a total of 35 existing) 2$\beta^-$ decay isotope candidates whose half-lives currently lack experimental limits. The $\alpha$ activity of the $^{146}$Sm daughter nuclide allows placement of limits on the 2$\beta$ decay using the $^{146}$Nd/$^{142}$Nd abundance ratio ($T^{\beta\beta}_{1/2} \gsim 3 \times 10^9$ yr) or direct search for $^{146}$Sm with accelerator mass spectrometers ($T^{\beta\beta}_{1/2} \gsim 4.5 \times 10^{19}$ yr). With a similar approach, a modest ($\sim$ Gyr) first limit on half-lives for the other unexplored $2\beta$ unstable isotopes and competitive limits (few $10^{13}$yr) for $^{98}$Mo and $^{122}$Sn are also given. Finally, it is shown how the limit $T^{0\nu\epsilon\beta^+}_{1/2} \gsim 10^{15}$ yr for the unexplored $^{144}$Sm $0\nu\epsilon \beta^+$ decay may be obtained from the data of a GSO crystal scintillator.


title: Decorated Dyck paths, the Delta conjecture, and a new q,t-square

abstract: We discuss the combinatorics of the decorated Dyck paths appearing in the Delta conjecture framework of Haglund, Remmel and Wilson (2015) and Zabrocki (2016), by introducing two new statistics, bounce and bounce'. We then provide plethystic formulae for their q,t-enumerators, by proving in this way a decorated version of Haglund's q,t-Schroeder theorem, answering a question in Haglund, Remmel and Wilson (2015). In particular we provide both an algebraic and a combinatorial explanation of a symmetry conjectured in Haglund, Remmel and Wilson (2015) and Zabrocki (2016). Finally, we show a link between these results and the q,t-square theorem of Can and Loehr (2006) (conjectured by Loehr and Warrington (2007)), getting a new q,t-square, whose specialization at t=1/q is identical to the classical one, and hinting to a broader phenomenon of this type.


title: A meson-baryon molecular interpretation for some $\Omega_c$ excited baryons

abstract: We explore the possibility that some of the five narrow $\Omega_c$ resonances recently observed at LHCb could correspond to pentaquark states, structured as meson-baryon bound states or molecules. The interaction of the low-lying pseudoscalar mesons with the ground-state baryons in the charm $+1$, strangeness $-2$ and isospin 0 sector is built from t-channel vector meson exchange, using effective Lagrangians. The resulting s-wave coupled-channel unitarized amplitudes show the presence of two structures with similar masses and widths to those of the observed $\Omega_c(3050)^0$ and $\Omega_c(3090)^0$. The identification of these resonances with the meson-baryon bound states found in this work would also imply assigning the values $1/2^-$ for their spin-parity. An experimental determination of the spin-parity of the $\Omega_c(3090)^0$ would help in disentangling its structure, as the quark-based models predict its spin-parity to be either $3/2^-$ or $5/2^-$.


title: A note on the uniqueness of the dRGT massive gravity theory

abstract: We revisit the problem of gravity coupled to a background metric $\eta_{\mu\nu}$, looking for ghost free interactions. It is known that elimination of the Boulware-Desewr ghost is equivalent to a certain Hessian condition on the interacting potential. We elaborate on this equation and provide a proof that dRGT potential is the unique interaction between the spacetime metric and a background field $\eta_{\mu\nu}$ which is both Lorentz invariant and ghost free. Our approach is fully non-perturbative using only the ADM formulation of general relativity.


title: Camera-Aware Multi-Resolution Analysis (CAMRA) for Raw Sensor Data Compression

abstract: We propose a novel lossless and lossy compression scheme for color filter array~(CFA) sampled images based on the wavelet transform of them. Our analysis suggests that the wavelet coefficients of HL and LH subbands are highly correlated. Hence, we decorrelate Mallat wavelet packet decomposition to further sparsify the coefficients. In addition, we develop a camera processing pipeline for compressing CFA sampled images aimed at maximizing the quality of the color images constructed from the compressed CFA sampled images. We validated our theoretical analysis and the performance of the proposed compression scheme using images of natural scenes captured in a raw format. The experimental results verify that our proposed method improves coding efficiency relative to the standard and the state-of-the-art compression schemes CFA sampled images.


title: On the error of a priori sampling: zero forcing sets and propagation time

abstract: Zero forcing is an iterative process on a graph used to bound the maximum nullity. The process begins with select vertices as colored, and the remaining vertices can become colored under a specific color change rule. The goal is to find a minimum set of vertices such that after iteratively applying the rule, all of the vertices become colored (i.e., a minimum zero forcing set). Of particular interest is the propagation time of a chosen set which is the number of steps the rule must be applied in order to color all the vertices of a graph. We give a purely linear algebraic interpretation of zero forcing: Find a set of vertices $S$ such that for any weighted adjacency matrix $\mathbf{A}$, whenever $\mathbf{Ax} = \mathbf{0}$, the entirety of of $\mathbf{x}$ can be recovered using only $\mathbf{x}_S$, the entries corresponding to $S$. The key here is that $S$ must be chosen before $\mathbf{A}$. In this light, we are able to give a linear algebraic interpretation of the propagation time: Any error in $\mathbf{x}_S$ effects the error of $\mathbf{x}$ exponentially in the propagation time. This error can be quantitatively measured using newly defined zero forcing-related parameters, the error polynomial vector and the variance polynomial vector. In this sense, the quality of two zero forcing sets can objectively be compared even if the sets are the same size and their propagation time is the same. Examples and constructions are given.


title: Light-induced switching in pDTE - FICO 1D photonic structures

abstract: We propose the design of 1D photonic crystals and microcavities in which fluorine-indium codoped cadmium oxide (FICO) nanocrystal based layers and layers of diarylethene-based polyester (pDTE) are alternated or embedded in a microcavity. The irradiation with UV light results in two different behaviours: i) it dopes the FICO nanocrystals inducing a blue shift of their plasmonic resonances; ii) it changes the real part of the refractive index of the photochromic pDTE polymer. These two behaviours are combined in the proposed photonic structures and can be useful for switchable filters and cavities for light emission.


title: POLAMI: Polarimetric Monitoring of Active Galactic Nuclei at Millimetre Wavelengths. I. The program, calibration, and calibrators data products

abstract: We describe the POLAMI program for the monitoring of all four Stokes parameters of a sample of bright radio-loud active galactic nuclei with the IRAM 30m telescope at 3.5 and 1.3mm. The program started in October 2006 and accumulated, until August 2014, 2300 observations at 3.5mm, achieving a median time sampling interval of 22 days for the sample of 37 sources. This first paper explains the source selection, mostly blazars, the observing strategy and data calibration, and gives the details of the instrumental polarisation corrections. The sensitivity (1sigma) reached at 3.5mm is 0.5% (linear polarisation degree), 4.7 deg. (polarisation angle), and 0.23% (circular polarisation), while the corresponding values at 1.3mm are 1.7%, 9.9 deg., and 0.72%, respectively. The data quality is demonstrated by the time sequences of our calibrators Mars and Uranus. For the quasar 3C286, widely used as a linear polarisation calibrator, we give improved estimates of its linear polarisation, and show for the first time occasional detections of its weak circular polarisation, which suggests a small level of variability of the source at millimeter wavelengths.


title: POLAMI: Polarimetric Monitoring of Active Galactic Nuclei at Millimetre Wavelengths. II. Widespread circular polarisation

abstract: We analyse the circular polarisation data accumulated in the first 7 years of the POLAMI project introduced in an accompanying paper (Agudo et al.). In the 3mm wavelength band, we acquired more than 2600 observations, and all but one of our 37 sample sources were detected, most of them several times. For most sources, the observed distribution of the degree of circular polarisation is broader than that of unpolarised calibrators, indicating that weak (<0.5%) circular polarisation is present most of the time. Our detection rate and the maximum degree of polarisation found, 2.0%, are comparable to previous surveys, all made at much longer wavelengths. We argue that the process generating circular polarisation must not be strongly wavelength dependent, and we propose that the widespread presence of circular polarisation in our short wavelength sample dominated by blazars is mostly due to Faraday conversion of the linearly polarised synchrotron radiation in the helical magnetic field of the jet. Circular polarisation is variable, most notably on time scales comparable to or shorter than our median sampling interval <1 month. Longer time scales of about one year are occasionally detected, but severely limited by the weakness of the signal. At variance with some longer wavelength investigations we find that the sign of circular polarisation changes in most sources, while only 7 sources, including 3 already known, have a strong preference for one sign. The degrees of circular and linear polarisation do not show any systematic correlation. We do find however one particular event where the two polarisation degrees vary in synchronism during a time span of 0.9 years. The paper also describes a novel method for calibrating the sign of circular polarisation observations.


title: POLAMI: Polarimetric Monitoring of Active Galactic Nuclei at Millimetre Wavelengths. III. Characterisation of total flux density and polarisation variability of relativistic jets

abstract: We report on the first results of the POLAMI program, a simultaneous 3.5 and 1.3mm full-Stokes-polarisation monitoring of a sample of 36 of the brightest active galactic nuclei in the northern sky with the IRAM 30m Telescope. Through a systematic statistical study of data taken from October 2006 (from December 2009 for the case of the 1.3mm observations) to August 2014, we characterise the variability of the total flux density and linear polarisation. We find that all sources in the sample are highly variable in total flux density at both 3.5 and 1.3mm, as well as in spectral index, that is found to be optically thin in general. The total flux-density variability at 1.3mm is found, in general, to be faster, and to have larger amplitude and flatter PSD slopes than 3.5mm. The polarisation degree is on average larger at 1.3mm than at 3.5mm, by a factor of 2.6. The variability of linear polarisation degree is faster and has higher fractional amplitude than for total flux density, with the typical time scales during prominent polarisation peaks being significantly faster at 1.3mm than at 3.5mm. The polarisation angle at both 3.5 and 1.3mm is highly variable. Most of the sources show one or two excursions of >180 deg. on time scales from a few weeks to about a year during the course of our observations. The 3.5 and 1.3mm polarisation angle evolution follow rather well each other, although the 1.3mm data show a clear preference to more prominent variability on the short time scales, i.e. weeks. The data are compatible with multi-zone models of conical jets involving smaller emission regions for the shortest-wavelength emitting sites. Such smaller emitting regions should also be more efficient in energising particle populations. The data also favours the integrated emission at 1.3mm to have better ordered magnetic fields than the one at 3.5mm.


title: The Replacement Property of PSL$(2,p)$ and PSL$(2,p^2)$

abstract: In 2014, Benjamin Nachman showed that when $p\equiv$1 mod 8, the 2-dimensional projective linear group over the field of $p$ elements fails the replacement property if the maximal length $m$ of an irredundant generating sequence for the group is 3. In addition, if $m=4$, the group satisfies the property for any prime $p$. In this paper, we will extend such classification for PSL$(2,q)$ where $q$ is $p$ and $p^2$ with conditions of modulo 8 and 10 on $p$.


title: DIeSEL: DIstributed SElf-Localization of a network of underwater vehicles

abstract: How can teams of artificial agents localize and position themselves in GPS-denied environments? How can each agent determine its position from pairwise ranges, own velocity, and limited interaction with neighbors? This paper addresses this problem from an optimization point of view: we directly optimize the nonconvex maximum-likelihood estimator in the presence of range measurements contaminated with Gaussian noise, and we obtain a provably convergent, accurate and distributed positioning algorithm that outperforms the extended Kalman filter, a standard centralized solution for this problem.


title: Simultaneous Measurements of Microwave Photoresistance and Cyclotron Reflection in the Multi-Photon Regime

abstract: We simultaneously measure photoresistance with electrical transport and plasmon-cyclotron resonance (PCR) using microwave reflection spectroscopy in high mobility GaAs/AlGaAs quantum wells under a perpendicular magnetic field. Multi-photon transitions are revealed as sharp peaks in the resistance and the cyclotron reflection on samples with various carrier densities. Our main finding is that plasmon coupling is relevant in the cyclotron reflection spectrum but has not been observed in the electrical conductivity signal. We discuss possible mechanisms relevant to reflection or dc conductivity signal to explain this discrepancy. We further confirm a trend that higher order multi-photon features can be observed using higher carrier density samples.


title: On Memory System Design for Stochastic Computing

abstract: Growing uncertainty in design parameters (and therefore, in design functionality) renders stochastic computing particularly promising, which represents and processes data as quantized probabilities. However, due to the difference in data representation, integrating conventional memory (designed and optimized for non-stochastic computing) in stochastic computing systems inevitably incurs a significant data conversion overhead. Barely any stochastic computing proposal to-date covers the memory impact. In this paper, as the first study of its kind to the best of our knowledge, we rethink the memory system design for stochastic computing. The result is a seamless stochastic system, StochMem, which features analog memory to trade the energy and area overhead of data conversion for computation accuracy. In this manner StochMem can reduce the energy (area) overhead by up-to 52.8% (93.7%) at the cost of at most 0.7% loss in computation accuracy.


title: The Shape of Velocity Dispersion Profiles and the Dynamical State of Galaxy Clusters

abstract: Motivated by the existence of the relationship between the dynamical state of clusters and the shape of the velocity dispersion profiles (VDP), we study the VDPs for Gaussian (G) and Non-Gaussian (NG) systems for a subsample of clusters from the Yang catalog. The groups cover a redshift interval of $0.03\leq z\leq0.1$ with halo mass $\geq 10^{14}$M$_{\odot}$. We use a robust statistical method, Hellinger Distance, to classify the dynamical state of the systems according to their velocity distribution. The stacked VDP of each class, G and NG, is then determined using either Bright or Faint galaxies. The stacked VDP for G groups displays a central peak followed by a monotonically decreasing trend which indicates a predominance of radial orbits, with the Bright stacked VDP showing lower velocity dispersions in all radii. The distinct features we find in NG systems are manifested not only by the characteristic shape of VDP, with a depression in the central region, but also by a possible higher infall rate associated with galaxies in the Faint stacked VDP.


title: Bobtail: A Proof-of-Work Target that Minimizes Blockchain Mining Variance (Draft)

abstract: Blockchain systems are designed to produce blocks at a constant average rate. The most popular systems currently employ a Proof of Work (PoW) algorithm as a means of creating these blocks. Bitcoin produces, on average, one block every 10 minutes. An unfortunate limitation of all deployed PoW blockchain systems is that the time between blocks has high variance. For example, 5% of the time, Bitcoin's inter-block time is at least 40 minutes. This variance impedes the consistent flow of validated transactions through the system. We propose an alternative process for PoW-based block discovery that results in an inter-block time with significantly lower variance. Our algorithm, called Bobtail, generalizes the current algorithm by comparing the mean of the k lowest order statistics to a target. We show that the variance of inter-block times decreases as k increases. If our approach were applied to Bitcoin, about 80% of blocks would be found within 7 to 12 minutes, and nearly every block would be found within 5 to 18 minutes; the average inter-block time would remain at 10 minutes. Further, we show that low-variance mining significantly thwarts doublespend and selfish mining attacks. For Bitcoin and Ethereum currently (k=1), an attacker with 40% of the mining power will succeed with 30% probability when the merchant sets up an embargo of 8 blocks; however, when k>=20, the probability of success falls to less than 1%. Similarly, for Bitcoin and Ethereum currently, a selfish miner with 40% of the mining power will claim about 66% of blocks; however, when k>=5, the same miner will find that selfish mining is less successful than honest mining. The cost of our approach is a larger block header.


title: Index divisibility in the orbit of 0 for integral polynomials

abstract: Let $f(x) \in \bbz[x]$ and consider the index divisibility set $D = \{n \in \bbn : n \mid f^n(0)\}$. We present a number of properties of $D$ in the case that $(f^n(0))_{n=1}^\infty$ is a rigid divisibility sequence, generalizing a number of results of Chen, Stange, and the first author. We then study the polynomial $x^d + x^e + c \in \bbz[x]$, where $d > e \ge 2$ and determine all cases where this map has a finite index divisibility set.


title: Doping Graphene via Organic Solid-Solid Wetting Deposition

abstract: Organic Solid-Solid Wetting Deposition (OSWD) enables the fabrication of supramolecular architectures without the need for solubility or vacuum conditions. The technique is based on a process which directly generates two-dimensional monolayers from three-dimensional solid organic powders. Consequently, insoluble organic pigments and semiconductors can be made to induce monolayer self-assembly on substrate surfaces, such as graphene and carbon nanotubes, under ambient conditions. The above factuality hence opens up the potential of the OSWD for bandgap engineering applications within the context of carbon based nanoelectronics. However, the doping of graphene via OSWD has not yet been verified, primarily owing to the fact that the classical OSWD preparation procedures do not allow for the analysis via Raman spectroscopy, one of the main techniques to determine graphene doping. Hence, here we describe a novel approach to induce OSWD on graphene leading to samples suitable for Raman spectroscopy. The analysis reveals peak shifts within the Raman spectrum of graphene, which are characteristics for p-type doping. Additional evidence for chemical doping is found via Scanning Tunneling Spectroscopy. The results open up a very easily applicable, low-cost, and eco-friendly way for doping graphene via commercially available organic pigments.


title: Automated Behavioral Analysis of Malware A Case Study of WannaCry Ransomware

abstract: Ransomware, a class of self-propagating malware that uses encryption to hold the victims' data ransom, has emerged in recent years as one of the most dangerous cyber threats, with widespread damage; e.g., zero-day ransomware WannaCry has caused world-wide catastrophe, from knocking U.K. National Health Service hospitals offline to shutting down a Honda Motor Company in Japan[1]. Our close collaboration with security operations of large enterprises reveals that defense against ransomware relies on tedious analysis from high-volume systems logs of the first few infections. Sandbox analysis of freshly captured malware is also commonplace in operation. We introduce a method to identify and rank the most discriminating ransomware features from a set of ambient (non-attack) system logs and at least one log stream containing both ambient and ransomware behavior. These ranked features reveal a set of malware actions that are produced automatically from system logs, and can help automate tedious manual analysis. We test our approach using WannaCry and two polymorphic samples by producing logs with Cuckoo Sandbox during both ambient, and ambient plus ransomware executions. Our goal is to extract the features of the malware from the logs with only knowledge that malware was present. We compare outputs with a detailed analysis of WannaCry allowing validation of the algorithm's feature extraction and provide analysis of the method's robustness to variations of input data\textemdash changing quality/quantity of ambient data and testing polymorphic ransomware. Most notably, our patterns are accurate and unwavering when generated from polymorphic WannaCry copies, on which 63 (of 63 tested) anti-virus (AV) products fail.


title: Theoretical analysis of electrostatic energy harvester configured as Bennet's doubler based on Q-V cycles

abstract: This paper presents theoretical analysis of a MEMS electrostatic energy harvester configured as the Bennet's doubler. Steady-state operation of the doubler circuit can be approximated by a right-angled trapezoid Q-V cycle. A similarity between voltage doubler and resistive-based charge-pump circuit is highlighted. By taking electromechanical coupling into account, the analytical solution of the saturation voltage is the first time derived, providing a greater comprehension of the system performance and multi-parameter effects. The theoretical approach is verified by results of circuit simulation for two cases of mathematically idealized diode and of Schottky diode. Development of the doubler/multiplier circuits that can further increase the saturation voltage is investigated.


title: Analytic approach to variance optimization under an $\ell_1$ constraint

abstract: The optimization of the variance supplemented by a budget constraint and an asymmetric $\ell_1$ regularizer is carried out analytically by the replica method borrowed from the theory of disordered systems. The asymmetric regularizer allows us to penalize short and long positions differently, so the present treatment includes the no-short-constrained portfolio optimization problem as a special case. Results are presented for the out-of-sample and the in-sample estimator of the regularized variance, the relative estimation error, the density of the assets eliminated from the portfolio by the regularizer, and the distribution of the optimal portfolio weights. We have studied the dependence of these quantities on the ratio $r$ of the portfolio's dimension $N$ to the sample size $T$, and on the strength of the regularizer. We have checked the analytic results by numerical simulations, and found general agreement. Regularization extends the interval where the optimization can be carried out, and suppresses the large sample fluctuations, but the performance of $\ell_1$ regularization is rather disappointing: if the sample size is large relative to the dimension, i.e. $r$ is small, the regularizer does not play any role, while for $r$'s where the regularizer starts to be felt the estimation error is already so large as to make the whole optimization exercise pointless. We find that the $\ell_1$ regularization can eliminate at most half the assets from the portfolio, corresponding to this there is a critical ratio $r=2$ beyond which the $\ell_1$ regularized variance cannot be optimized: the regularized variance becomes constant over the simplex. These facts do not seem to have been noticed in the literature.


title: Monotonicity and local uniqueness for the Helmholtz equation

abstract: This work extends monotonicity-based methods in inverse problems to the case of the Helmholtz (or stationary Schr\"odinger) equation $(\Delta + k^2 q) u = 0$ in a bounded domain for fixed non-resonance frequency $k>0$ and real-valued scattering coefficient function $q$. We show a monotonicity relation between the scattering coefficient $q$ and the local Neumann-Dirichlet operator that holds up to finitely many eigenvalues. Combining this with the method of localized potentials, or Runge approximation, adapted to the case where finitely many constraints are present, we derive a constructive monotonicity-based characterization of scatterers from partial boundary data. We also obtain the local uniqueness result that two coefficient functions $q_1$ and $q_2$ can be distinguished by partial boundary data if there is a neighborhood of the boundary where $q_1\geq q_2$ and $q_1\not\equiv q_2$.


title: Flat families of point schemes for connected graded algebras

abstract: We study truncated point schemes of connected graded algebras as families over the parameter space of varying relations for the algebras, proving that the families are flat over the open dense locus where the point schemes achieve the expected (i.e. minimal) dimension. When the truncated point scheme is zero-dimensional we obtain its number of points counted with multiplicity via a Chow ring computation. This latter application in particular confirms a conjecture of Brazfield to the effect that a generic two-generator, two-relator 4-dimensional Artin-Schelter regular algebra has seventeen truncated point modules of length six.


title: An application of kissing numbers in sum-product estimates

abstract: The boundedness of the kissing numbers of convex bodies has been known to Hadwiger for long. We present an application of it to the sum-product estimate $$\max(|\mathcal{A}+\mathcal{A}|,|\mathcal{A}\mathcal{A}|)\gg\dfrac{|\mathcal{A}|^{4/3}}{\lceil\log|\mathcal{A}|\rceil^{1/3}}$$ for finite sets $\mathcal{A}$ of quaternions and of a certain family of well-conditioned matrices.


title: Analysis of MEMS electrostatic energy harvesters electrically configured as voltage multipliers

abstract: This paper presents the analysis of an efficient alternative interface circuit for MEMS electrostatic energy harvesters. It is entirely composed by diodes and capacitors. Based on modeling and simulation, the anti-phase gap-closing structure is investigated. We find that when configured as a voltage multiplier, it can operate at very low acceleration amplitudes. In addition, the allowed maximum voltage between electrodes is barely limited by the pull-in effect. The parasitic capacitance of the harvester and non-ideal characteristics of electronic components are taken into account. A lumped-model of the harvesting system has been implemented in a circuit simulator. Simulation results show that an output voltage of 22 V is obtained with 0.15 g input acceleration. The minimum necessary ratio between the maximum and minimum capacitances of the generators which allows operation of the circuit, can be lower than 2. This overcomes a crucial obstacle in low-power energy harvesting devices. A comparison between the voltage multiplier against other current topologies is highlighted. An advantage of the former over the latter is to generate much higher saturation voltage, while the minimum required initial bias and the minimum capacitance ratio in both cases are in the similar levels.


title: On the structure of cyclotomic nilHecke algebras

abstract: In this paper we study the structure of the cyclotomic nilHecke algebras $\HH_{\ell,n}^{(0)}$, where $\ell,n\in\N$. We construct a monomial basis for $\HH_{\ell,n}^{(0)}$ which verifies a conjecture of Mathas. We show that the graded basic algebra of $\HH_{\ell,n}^{(0)}$ is commutative and hence isomorphic to the center $Z$ of $\HH_{\ell,n}^{(0)}$. We further prove that $\HH_{\ell,n}^{(0)}$ is isomorphic to the full matrix algebra over $Z$ and construct an explicit basis for the center $Z$. We also construct a complete set of pairwise orthogonal primitive idempotents of $\HH_{\ell,n}^{(0)}$. Finally, we present a new homogeneous symmetrizing form $\Tr$ on $\HH_{\ell,n}^{(0)}$ by explicitly specifying its values on a given homogeneous basis of $\HH_{\ell,n}^{(0)}$ and show that it coincides with Shan--Varagnolo--Vasserot's symmetrizing form $\Tr^{\text{SVV}}$ on $\HH_{\ell,n}^{(0)}$.


title: Image similarity using Deep CNN and Curriculum Learning

abstract: Image similarity involves fetching similar looking images given a reference image. Our solution called SimNet, is a deep siamese network which is trained on pairs of positive and negative images using a novel online pair mining strategy inspired by Curriculum learning. We also created a multi-scale CNN, where the final image embedding is a joint representation of top as well as lower layer embedding's. We go on to show that this multi-scale siamese network is better at capturing fine grained image similarities than traditional CNN's.


title: Chandra and NuSTAR follow-up observations of Swift-BAT selected AGN

abstract: Based on current models of the Cosmic X-ray Background (CXB), heavily obscured Active Galactic Nuclei (AGN) are expected to make up ~10% of the peak emission of the CXB and ~20% of the total population of AGN, yet few of these sources have been recorded and characterized in current surveys. Here we present the Chandra follow-up observation of 14 AGN detected by Swift-BAT. For five sources in the sample, NuSTAR observations in the 3-80 keV band are also available. The X-ray spectral fitting over the 0.3-150 keV energy range allows us to determine the main X-ray spectral parameters, such as the photon index and the intrinsic absorption, of these objects, and to make hypotheses on the physical structures responsible for the observed spectra. We find that 13 of the 14 objects are absorbed AGN, and one is a candidate Compton thick AGN, having intrinsic absorption NH>1E24 cm^{-2}. Finally, we verified that the use of NuSTAR observations is strategic to strongly constrain the properties of obscured AGN, since the best-fit values we obtained for parameters such as the power-law photon index and the intrinsic absorption NH changed sometimes significantly fitting the spectra with and without the use of NuSTAR data.


title: Encoding Bitrate Optimization Using Playback Statistics for HTTP-based Adaptive Video Streaming

abstract: HTTP video streaming is in wide use to deliver video over the Internet. With HTTP adaptive steaming, a video playback dynamically selects a video stream from a pre-encoded representation based on available bandwidth and viewport (screen) size. The viewer's video quality is therefore influenced by the encoded bitrates. We minimize the average delivered bitrate subject to a quality lower bound on a per-chunk basis by modeling the probability that a player selects a particular encoding. Through simulation and real-world experiments, the proposed method saves 9.6% of bandwidth while average delivered video quality comparing with state of the art while keeping average delivered video quality.


title: The importance of scale in spatially varying coefficient modeling

abstract: While spatially varying coefficient (SVC) models have attracted considerable attention in applied science, they have been criticized as being unstable. The objective of this study is to show that capturing the "spatial scale" of each data relationship is crucially important to make SVC modeling more stable, and in doing so, adds flexibility. Here, the analytical properties of six SVC models are summarized in terms of their characterization of scale. Models are examined through a series of Monte Carlo simulation experiments to assess the extent to which spatial scale influences model stability and the accuracy of their SVC estimates. The following models are studied: (i) geographically weighted regression (GWR) with a fixed distance or (ii) an adaptive distance bandwidth (GWRa), (iii) flexible bandwidth GWR (FB-GWR) with fixed distance or (iv) adaptive distance bandwidths (FB-GWRa), (v) eigenvector spatial filtering (ESF), and (vi) random effects ESF (RE-ESF). Results reveal that the SVC models designed to capture scale dependencies in local relationships (FB-GWR, FB-GWRa and RE-ESF) most accurately estimate the simulated SVCs, where RE-ESF is the most computationally efficient. Conversely GWR and ESF, where SVC estimates are naively assumed to operate at the same spatial scale for each relationship, perform poorly. Results also confirm that the adaptive bandwidth GWR models (GWRa and FB-GWRa) are superior to their fixed bandwidth counterparts (GWR and FB-GWR).


title: Network Topology and Communication-Computation Tradeoffs in Decentralized Optimization

abstract: In decentralized optimization, nodes cooperate to minimize an overall objective function that is the sum (or average) of per-node private objective functions. Algorithms interleave local computations with communication among all or a subset of the nodes. Motivated by a variety of applications---distributed estimation in sensor networks, fitting models to massive data sets, and distributed control of multi-robot systems, to name a few---significant advances have been made towards the development of robust, practical algorithms with theoretical performance guarantees. This paper presents an overview of recent work in this area. In general, rates of convergence depend not only on the number of nodes involved and the desired level of accuracy, but also on the structure and nature of the network over which nodes communicate (e.g., whether links are directed or undirected, static or time-varying). We survey the state-of-the-art algorithms and their analyses tailored to these different scenarios, highlighting the role of the network topology.


title: Can math beat gamers in Quantum Moves?

abstract: In a recent work on quantum state preparation, S{\o}rensen and colleagues explore the possibility of using video games to help design quantum control protocols. The authors present a game called "Quantum Moves" in which gamers have to move an atom from A to B by means of optical tweezers. They report that, players succeed where purely numerical optimization fails [1]. Moreover, by harnessing the player strategies they can outperform the most prominent established numerical methods [1]. The aim of this manuscript is to analyze the problem in detail and show that those claims are untenable. In fact a simple stochastic local optimization method can easily find very good solutions to this problem in a few 1000 trials rather than the astronomical $7.4\times 10^{8}$ trials of the most successful optimization method reported in [1]. Next, counter-diabatic driving is used to generate protocols without resorting to numeric optimization; the protocols are shown to outperform virtually all players. The analysis moreover results in an accurate analytic estimate of the quantum speed limit which, apart from zero-point motion, is shown to be entirely classical in nature. The latter might explain why gamers are remarkably good at the game.


title: BOSS-LDG: A Novel Computational Framework that Brings Together Blue Waters, Open Science Grid, Shifter and the LIGO Data Grid to Accelerate Gravitational Wave Discovery

abstract: We present a novel computational framework that connects Blue Waters, the NSF-supported, leadership-class supercomputer operated by NCSA, to the Laser Interferometer Gravitational-Wave Observatory (LIGO) Data Grid via Open Science Grid technology. To enable this computational infrastructure, we configured, for the first time, a LIGO Data Grid Tier-1 Center that can submit heterogeneous LIGO workflows using Open Science Grid facilities. In order to enable a seamless connection between the LIGO Data Grid and Blue Waters via Open Science Grid, we utilize Shifter to containerize LIGO's workflow software. This work represents the first time Open Science Grid, Shifter, and Blue Waters are unified to tackle a scientific problem and, in particular, it is the first time a framework of this nature is used in the context of large scale gravitational wave data analysis. This new framework has been used in the last several weeks of LIGO's second discovery campaign to run the most computationally demanding gravitational wave search workflows on Blue Waters, and accelerate discovery in the emergent field of gravitational wave astrophysics. We discuss the implications of this novel framework for a wider ecosystem of Higher Performance Computing users.


title: A recent origin for Saturn's rings from the collisional disruption of an icy moon

abstract: The disruption of an icy moon in a collision with an interloping comet a few hundred million years ago is a simple way to create Saturn's rings. A ring parent moon with a mass comparable to Mimas could be trapped in mean motion resonance with Enceladus and Dione in an orbit near the current outer edge of the rings just beyond the Roche zone. I present collisional N-body simulations of cometary impacts that lead to the partial disruption of a differentiated moon with a rocky core and icy mantle. The core can survive largely intact while the debris from the mantle settles into a ring of predominantly ice particles straddling the orbital radius of the parent moon. The nascent ring spreads radially due to collisional viscosity while mass re-accretes onto the remnant rocky core to form a new moon that can be identified as Mimas. The icy debris that migrates into the Roche zone evolves into Saturn's ring system. Torques from tidal interaction with Saturn and resonant interactions with the rings push the recently formed Mimas outward to its current position on the same timescale of a few hundred million years. This scenario accounts for the high ice fraction observed in Saturn's rings and explains why the ring mass is comparable to the mass of Mimas. The prior existence of a ring parent moon in mean motion resonance results in a tidal heating rate for Enceladus in the recent past that is significantly larger than the current rate.


title: Green rings of Drinfeld Doubles of Taft algebras

abstract: In this article, we investigate the representation ring (or Green ring) of the Drinfeld double $D(H_n(q))$ of the Taft algebra $H_n(q)$, where $n$ is an integer with $n>2$ and $q$ is a root of unity of order $n$. It is shown that the Green ring $r(D(H_n(q)))$ is a commutative ring generated by infinitely many elements subject to certain relations.


title: On the Model Shrinkage Effect of Gamma Process Edge Partition Models

abstract: The edge partition model (EPM) is a fundamental Bayesian nonparametric model for extracting an overlapping structure from binary matrix. The EPM adopts a gamma process ($\Gamma$P) prior to automatically shrink the number of active atoms. However, we empirically found that the model shrinkage of the EPM does not typically work appropriately and leads to an overfitted solution. An analysis of the expectation of the EPM's intensity function suggested that the gamma priors for the EPM hyperparameters disturb the model shrinkage effect of the internal $\Gamma$P. In order to ensure that the model shrinkage effect of the EPM works in an appropriate manner, we proposed two novel generative constructions of the EPM: CEPM incorporating constrained gamma priors, and DEPM incorporating Dirichlet priors instead of the gamma priors. Furthermore, all DEPM's model parameters including the infinite atoms of the $\Gamma$P prior could be marginalized out, and thus it was possible to derive a truly infinite DEPM (IDEPM) that can be efficiently inferred using a collapsed Gibbs sampler. We experimentally confirmed that the model shrinkage of the proposed models works well and that the IDEPM indicated state-of-the-art performance in generalization ability, link prediction accuracy, mixing efficiency, and convergence speed.


title: On the lateral instability analysis of MEMS comb-drive electrostatic transducers

abstract: This paper investigates the lateral pull-in effect of an in-plane overlap-varying transducer. The instability is induced by the translational and rotational displacements. Based on the principle of virtual work, the equilibrium conditions of force and moment in lateral directions are derived. The analytical solutions of the critical voltage, at which the pull-in phenomenon occurs, are developed when considering only the translational stiffness or only the rotational stiffness of the mechanical spring. The critical voltage in general case is numerically determined by using nonlinear optimization techniques, taking into account the combined effect of translation and rotation. The effects of possible translational offsets and angular deviations to the critical voltage are modeled and numerically analyzed. The investigation is then the first time expanded to anti-phase operation mode and Bennet's doubler configuration of the two transducers.


title: Dynamic Reconfiguration of Mission Parameters in Underwater Human-Robot Collaboration

abstract: This paper presents a real-time programming and parameter reconfiguration method for autonomous underwater robots in human-robot collaborative tasks. Using a set of intuitive and meaningful hand gestures, we develop a syntactically simple framework that is computationally more efficient than a complex, grammar-based approach. In the proposed framework, a convolutional neural network is trained to provide accurate hand gesture recognition; subsequently, a finite-state machine-based deterministic model performs efficient gesture-to-instruction mapping, and further improves robustness of the interaction scheme. The key aspect of this framework is that it can be easily adopted by divers for communicating simple instructions to underwater robots without using artificial tags such as fiducial markers, or requiring them to memorize a potentially complex set of language rules. Extensive experiments are performed both on field-trial data and through simulation, which demonstrate the robustness, efficiency, and portability of this framework in a number of different scenarios. Finally, a user interaction study is presented that illustrates the gain in usability of our proposed interaction framework compared to the existing methods for underwater domains.


title: Tensor network subspace identification of polynomial state space models

abstract: This article introduces a tensor network subspace algorithm for the identification of specific polynomial state space models. The polynomial nonlinearity in the state space model is completely written in terms of a tensor network, thus avoiding the curse of dimensionality. We also prove how the block Hankel data matrices in the subspace method can be exactly represented by low rank tensor networks, reducing the computational and storage complexity significantly. The performance and accuracy of our subspace identification algorithm are illustrated by numerical experiments, showing that our tensor network implementation is around 20 times faster than the standard matrix implementation before the latter fails due to insufficient memory, is robust with respect to noise and can model real-world systems.


title: Exploring the Design Space of Immersive Urban Analytics

abstract: Recent years have witnessed the rapid development and wide adoption of immersive head-mounted devices, such as HTC VIVE, Oculus Rift, and Microsoft HoloLens. These immersive devices have the potential to significantly extend the methodology of urban visual analytics by providing critical 3D context information and creating a sense of presence. In this paper, we propose an theoretical model to characterize the visualizations in immersive urban analytics. Further more, based on our comprehensive and concise model, we contribute a typology of combination methods of 2D and 3D visualizations that distinguish between linked views, embedded views, and mixed views. We also propose a supporting guideline to assist users in selecting a proper view under certain circumstances by considering visual geometry and spatial distribution of the 2D and 3D visualizations. Finally, based on existing works, possible future research opportunities are explored and discussed.


title: Estimation of the Hurst Exponent Using Trimean Estimators on Nondecimated Wavelet Coefficients

abstract: Hurst exponent is an important feature summarizing the noisy high-frequency data when the inherent scaling pattern cannot be described by standard statistical models. In this paper, we study the robust estimation of Hurst exponent based on non-decimated wavelet transforms (NDWT). The robustness is achieved by applying a general trimean estimator on non-decimated wavelet coefficients of the transformed data. The general trimean estimator is derived as a weighted average of the distribution's median and quantiles, combining the median's emphasis on central values with the quantiles' attention to the extremes. The properties of the proposed Hurst exponent estimators are studied both theoretically and numerically. Compared with other standard wavelet-based methods (Veitch $\&$ Abry (VA) method, Soltani, Simard, $\&$ Boichu (SSB) method, median based estimators MEDL and MEDLA), our methods reduce the variance of the estimators and increase the prediction precision in most cases. The proposed methods are applied to a data set in high frequency pupillary response behavior (PRB) with the goal to classify individuals according to a degree of their visual impairment.


title: Neglecting Model Structural Uncertainty Underestimates Upper Tails of Flood Hazard

abstract: Coastal flooding drives considerable risks to many communities, but projections of future flood risks are deeply uncertain. The paucity of observations of extreme events often motivates the use of statistical approaches to model the distribution of extreme storm surge events. A key deep uncertainty that is often overlooked is model structural uncertainty. There is currently no strong consensus among experts regarding which class of statistical model to use as a best practice. Robust management of coastal flooding risks requires coastal managers to consider the distinct possibility of non-stationarity in storm surges. This increases the complexity of the potential models to use, which tends to increase the data required to constrain the model. Here, we use a Bayesian model averaging approach to analyze the balance between model complexity sufficient to capture decision-relevant risks and data availability to constrain complex model structures. We characterize deep model structural uncertainty through a set of calibration experiments. Specifically, we calibrate a set of models ranging in complexity using long-term tide gauge observations from the Netherlands and the United States. We find that in both cases, roughly half the model weight is associated with non-stationary models. Our approach provides a formal framework to integrate information across model structures, in light of the potentially sizable modeling uncertainties. By combining information from multiple models, our inference sharpens for the projected storm surge 100-year return levels, and estimated return levels increase by several centimeters. We assess the impacts of data availability through a set of experiments with temporal subsets and model comparison metrics. Our analysis suggests about 70 years of data are required to stabilize estimates of the 100-year return level, for the locations and methods considered here.


title: Magnetic field and radial velocities of the star Chi Draconis A

abstract: We present high-resolution spectropolarimetric observations of the spectroscopic binary Chi Dra. Spectral lines in the spectrum of the main component Chi Dra A show variable Zeeman displacement, which confrms earlier suggestions about the presence of a weak magnetic field on the surface of this star. Within about 2 years of time base of our observations, the longitudinal component BL of the magnetic field exhibits variation from -11.5 +/- 2.5 G to +11.1 +/- 2.1 G with a period of about 23 days. Considering the rotational velocity of Chi DraA in the literature and that newly measured in this work, this variability may be explained by the stellar rotation under the assumption that the magnetic field is globally stable. Our new measurements of the radial velocities (RV) in high-resolution I-spectra of Chi Dra A refined the orbital parameters and reveal persistent deviations of RVs from the orbital curve. We suspect that these deviations may be due to the in uence of local magnetically generated spots, pulsations, or a Jupiter-size planet orbiting the system.


title: SiO Maser Survey towards off-plane O-rich AGBs around the orbital plane of the Sagittarius Stellar Stream

abstract: We conducted an SiO maser survey towards 221 O-rich AGB stars with the aim of identifying maser emission associated with the Sagittarius stellar stream. In this survey, maser emission was detected in 44 targets, of which 35 were new detections. All of these masers are within 5 kpc of the Sun. We also compiled a Galactic SiO maser catalogue including ~2300 SiO masers from the literature. The distribution of these SiO masers give a scale height of 0.40 kpc, while 42 sources deviate from the Galactic plane by more than 1.2 kpc, half of which were found in this survey. Regarding SiO masers in the disc, we found both the rotational speeds and the velocity dispersions vary with the Galactic plane distance. Assuming Galactic rotational speed $\Theta$0 = 240 km/s , we derived the velocity lags are 15 km/s and 55 km/s for disc and off-plane SiO masers respectively. Moreover, we identified three groups with significant peculiar motions (with 70% confidence). The most significant group is in the thick disc that might trace stream/peculiar motion of the Perseus arm. The other two groups are mainly made up of off-plane sources. The northern and southern off-plane sources were found to be moving at ~33 km/s and 54 km/s away from the Galactic plane, respectively. Causes of these peculiar motions are still unclear. For the two off-plane groups, we suspect they are thick disc stars whose kinematics affected by the Sgr stellar stream or very old Sgr stream debris.


title: Two-Dimensional Molecular Gas and ongoing star formation around H II region Sh2-104

abstract: We performed a multiwavelength study towards HII region Sh2-104. New maps of 12CO J=1-0 and 13CO J=1-0 were obtained from the Purple Mountain Observatory (PMO) 13.7 m radio telescope. Sh2-104 displays a double-ring structure. The outer ring with a radius of 4.4 pc is dominated by 12 um, 500 um, 12CO J=1-0, and 13CO J=1-0 emission, while the inner ring with a radius of 2.9 pc is dominated by 22 um and 21 cm emission. We did not detect CO emission inside the outer ring. The north-east portion of the outer ring is blueshifted, while the south-west portion is redshifted. The present observations have provided evidence that the collected outer ring around Sh2-104 is a two-dimensional structure. From the column density map constructed by the Hi-GAL survey data, we extract 21 clumps. About 90\% of all the clumps will form low-mass stars. A power-law fit to the clumps yields M=281Msun(r/pc)^1.31. The selected YSOs are associated with the collected material on the edge of Sh2-104. The derived dynamical age of Sh2-104 is 1.6*10^6 yr. Compared the Sh2-104 dynamical age with the YSOs timescale and the fragmentation time of the molecular ring, we further confirm that collect-and-collapse process operates in this region, indicating a positive feedback from a massive star for surrounding gas.


title: Convex Relaxations for Global Optimization Under Uncertainty Described by Continuous Random Variables

abstract: This article considers nonconvex global optimization problems subject to uncertainties described by continuous random variables. Such problems arise in chemical process design, renewable energy systems, stochastic model predictive control, etc. Here, we restrict our attention to problems with expected-value objectives and no recourse decisions. In principle, such problems can be solved globally using spatial branch-and-bound (B&B). However, B&B requires the ability to bound the optimal objective value on subintervals of the search space, and existing techniques are not generally applicable because expected-value objectives often cannot be written in closed-form. To address this, this article presents a new method for computing convex and concave relaxations of nonconvex expected-value functions, which can be used to obtain rigorous bounds for use in B&B. Furthermore, these relaxations obey a second-order pointwise convergence property, which is sufficient for finite termination of B&B under standard assumptions. Empirical results are shown for three simple examples.


title: A new mechanism for generating particle number asymmetry through interactions

abstract: A new mechanism for generating particle number asymmetry (PNA) has been developed. This mechanism is realized with a Lagrangian including a complex scalar field and a neutral scalar field. The complex scalar carries U(1) charge which is associated with the PNA. It is written in terms of the condensation and Green's function, which is obtained with two-particle irreducible (2PI) closed time path (CTP) effective action (EA). In the spatially flat universe with a time-dependent scale factor, the time evolution of the PNA is computed. We start with an initial condition where only the condensation of the neutral scalar is non-zero. The initial condition for the fields is specified by a density operator parameterized by the temperature of the universe. With the above initial conditions, the PNA vanishes at the initial time and later it is generated through the interaction between the complex scalar and the condensation of the neutral scalar. We investigate the case that both the interaction and the expansion rate of the universe are small and include their effects up to the first order of the perturbation. The expanding universe causes the effects of the dilution of the PNA, freezing interaction and the redshift of the particle energy. As for the time dependence of the PNA, we found that PNA oscillates at the early time and it begins to dump at the later time. The period and the amplitude of the oscillation depend on the mass spectrum of the model, the temperature and the expansion rate of the universe.


title: The Projective Class Rings of a family of pointed Hopf algebras of Rank two

abstract: In this paper, we compute the projective class rings of the tensor product $\mathcal{H}_n(q)=A_n(q)\otimes A_n(q^{-1})$ of Taft algebras $A_n(q)$ and $A_n(q^{-1})$, and its cocycle deformations $H_n(0,q)$ and $H_n(1,q)$, where $n>2$ is a positive integer and $q$ is a primitive $n$-th root of unity. It is shown that the projective class rings $r_p(\mathcal{H}_n(q))$, $r_p(H_n(0,q))$ and $r_p(H_n(1,q))$ are commutative rings generated by three elements, three elements and two elements subject to some relations, respectively. It turns out that even $\mathcal{H}_n(q)$, $H_n(0,q)$ and $H_n(1,q)$ are cocycle twist-equivalent to each other, they are of different representation types: wild, wild and tame, respectively.


title: Green Heterogeneous Cloud Radio Access Networks: Potential Techniques, Performance Tradeoffs, and Challenges

abstract: As a flexible and scalable architecture, heterogeneous cloud radio access networks (H-CRANs) inject strong vigor into the green evolution of current wireless networks. But the brutal truth is that energy efficiency (EE) improves at the cost of other indexes such as spectral efficiency (SE), fairness, and delay. It is thus important to investigate performance tradeoffs for striking flexible balances between energy-efficient transmission and excellent quality-of-service (QoS) guarantees under this new architecture. In this article, we first propose some potential techniques to energy-efficiently operate H-CRANs by exploiting their features. We then elaborate the initial ideas of modeling three fundamental tradeoffs, namely EE-SE, EE-fairness, and EE-delay tradeoffs, when applying these green techniques, and present open issues and challenges for future investigations. These related results are expected to shed light on green operation of H-CRANs from adaptive resource allocation, intelligent network control, and scalable network planning.


title: Renyi Entropy of Chaotic Eigenstates

abstract: Using arguments built on ergodicity, we derive an analytical expression for the Renyi entanglement entropies corresponding to the finite-energy density eigenstates of chaotic many-body Hamiltonians. The expression is a universal function of the density of states and is valid even when the subsystem is a finite fraction of the total system - a regime in which the reduced density matrix is not thermal. We find that in the thermodynamic limit, only the von Neumann entropy density is independent of the subsystem to the total system ratio $V_A/V$, while the Renyi entropy densities depend non-linearly on $V_A/V$. Surprisingly, Renyi entropies $S_n$ for $n > 1$ are convex functions of the subsystem size, with a volume law coefficient that depends on $V_A/V$, and exceeds that of a thermal mixed state at the same energy density. We provide two different arguments to support our results: the first one relies on a many-body version of Berry's formula for chaotic quantum mechanical systems, and is closely related to eigenstate thermalization hypothesis. The second argument relies on the assumption that for a fixed energy in a subsystem, all states in its complement allowed by the energy conservation are equally likely. We perform Exact Diagonalization study on quantum spin-chain Hamiltonians to test our analytical predictions, and find good agreement.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Mixed Precision Solver Scalable to 16000 MPI Processes for Lattice Quantum Chromodynamics Simulations on the Oakforest-PACS System

abstract: Lattice Quantum Chromodynamics (Lattice QCD) is a quantum field theory on a finite discretized space-time box so as to numerically compute the dynamics of quarks and gluons to explore the nature of subatomic world. Solving the equation of motion of quarks (quark solver) is the most compute-intensive part of the lattice QCD simulations and is one of the legacy HPC applications. We have developed a mixed-precision quark solver for a large Intel Xeon Phi (KNL) system named "Oakforest-PACS", employing the $O(a)$-improved Wilson quarks as the discretized equation of motion. The nested-BiCGSTab algorithm for the solver was implemented and optimized using mixed-precision, communication-computation overlapping with MPI-offloading, SIMD vectorization, and thread stealing techniques. The solver achieved 2.6 PFLOPS in the single-precision part on a $400^3\times 800$ lattice using 16000 MPI processes on 8000 nodes on the system.


title: Photo-induced change of refractive index and transparency in Bi2Te3 films

abstract: We report on an 800 nm femtosecond laser beam induced giant refractive index modulation and enhancement of near-infrared transparency in topological insulator material Bi2Te3 thin films. An ultrahigh refractive index of up to 5.9 was observed in the Bi2Te3 thin film in near-infrared frequency. The refractive index dramatically decreases by a factor of ~ 3 by an exposure to the 800 nm femtosecond laser beam. Simultaneously, the transmittance of the Bi2Te3 thin films markedly increases to ~ 96% in the near-infrared frequency. The Raman spectra provides strong evidences that the observed both refractive index modulation and transparency enhancement result from laser beam induced photooxidation effects in the Bi2Te3 thin films. The Bi2Te3 compound transfers into Bi2O3 and TeO2 under the laser beam illumination. These experimental results pave the way towards the design of various optical devices, such as near-infrared flat lenses, waveguide and holograms, based on topological insulator materials.


title: Liouville action and Holography on quasi-Fuchsian deformation spaces

abstract: We study the Liouville action for quasi-Fuchsian groups with parabolic and elliptic elements. In particular, when the group is Fuchsian, the contribution of elliptic elements to the classical Liouville action is derived in terms of the Bloch-Wigner functions. We prove the first and second variation formulas for the classical Liouville action on the quasi-Fuchsian deformation space. We prove an equality expressing the holography principle, which relates the Liouville action and the renormalized volume for quasi-Fuchsian groups with parabolic and elliptic elements. We also construct the potential functions of the K\"ahler forms corresponding to the Takhtajan-Zograf metrics associated to the elliptic elements in the quasi-Fuchsian groups.


title: Nanoscale Bandgap Tuning across an Inhomogeneous Ferroelectric Interface

abstract: We report nanoscale bandgap engineering via a local strain across the inhomogeneous ferroelectric interface, which is controlled by the visible-light-excited probe voltage. Switchable photovolatic effects and the spectral response of the photocurrent were explore to illustrate the reversible bandgap variation (~0.3eV). This local-strain-engineered bandgap has been further revealed by in situ probe-voltage-assisted valence electron energy-loss spectroscopy (EELS). Phase-field simulations and first-principle calculations were also employed for illustration of the large local strain and the bandgap variation in ferroelectric perovskite oxides. This reversible bandgap tuning in complex oxides demonstrates a framework for the understanding of the opticallyrelated behaviors (photovoltaic, photoemission, and photocatalyst effects) affected by order parameters such as charge, orbital, and lattice parameters.


title: Molecular Beam Epitaxy Growth of Tetragonal FeS Film on SrTiO3(001) Substrate

abstract: We report the successful growth of tetragonal FeS film with one or two unit-cell (UC) thickness on SrTiO3(001) substrate by molecular beam epitaxy. Large lattice constant mismatch with the substrate leads to high density of defects in single UC FeS, while it has been significantly reduced in double UC thick film due to the lattice relaxation. The scanning tunneling spectra on the surface of FeS thin film reveal the electronic doping effect of single UC FeS from the substrate. In addition, at the Fermi level, the energy gaps of approximate 1.5 meV are observed in films of both thicknesses at 4.6 K and below. The absence of coherence peaks of gap spectra may be related to the preformed Cooper-pairs without phase coherence.


title: Light-Enhanced Spin Fluctuations and d-Wave Superconductivity at a Phase Boundary

abstract: Time-domain techniques have shown the potential of photo-manipulating existing orders and inducing new states of matter in strongly correlated materials. Using time-resolved exact diagonalization, we perform numerical studies of pump dynamics in a Mott-Peierls system with competing charge and spin density waves. A light-enhanced $d$-wave superconductivity is observed when the system resides near a quantum phase boundary. By examining the evolution of spin, charge and superconducting susceptibilities, we show that a sub-dominant state in equilibrium can be stabilized by photomanipulating charge order to allow superconductivity to appear and dominate. This work provides an interpretation of light-induced superconductivity from the perspective of order competition, and offers a promising approach for designing novel emergent states out of equilibrium.


title: Exploiting Dual Connectivity in Heterogeneous Cellular Networks

abstract: We consider network utility maximization problems over heterogeneous cellular networks (HetNets) that permit dual connectivity. Dual connectivity (DC) is a feature that targets emerging practical HetNet deployments that will comprise of non-ideal (higher latency) connections between transmission nodes, and has been recently introduced to the LTE-Advanced standard. DC allows for a user to be simultaneously served by a macro node as well as one other (typically micro or pico) node and requires relatively coarser level coordination among serving nodes. For such a DC enabled HetNet we comprehensively analyze the problem of determining an optimal user association, where in any feasible association each user can be associated with (i.e., configured to receive data from) any one macro node (in a given set of macro nodes) and any one pico node that lies in the chosen macro node's coverage area. We consider the weighted sum rate system utility subject to per-user maximum and minimum rate constraints, as well as the proportional fairness (PF) system utility. For both utility choices we construct approximation algorithms and establish their respective approximation guarantees. We then validate the performance of our algorithms via numerical results.


title: Closure of resource-bounded randomness notions under polynomial time permutations

abstract: An infinite bit sequence is called recursively random if no computable strategy betting along the sequence has unbounded capital. It is well-known that the property of recursive randomness is closed under computable permutations. We investigate analogous statements for randomness notions defined by betting strategies that are computable within resource bounds. Suppose that S is a polynomial time computable permutation of the set of strings over the unary alphabet (identified with N). If the inverse of S is not polynomially bounded, it is easy to build a polynomial time random bit sequence Z such that Z o S is not polynomial time random. So one should only consider permutations S satisfying the extra condition that the inverse is polynomially bounded. Now the closure depends on additional assumptions in complexity theory. Our first result shows that if BPP contains a superpolynomial deterministic time class then polynomial time randomness is not preserved by some permutation S such that in fact both S and its inverse are in P. Our second result shows that polynomial space randomness is preserved by polynomial time permutations with polynomially boun\-ded inverse, so if P=PSPACE then polynomial time randomness is preserved.


title: Proton pairing in neutron stars from chiral effective field theory

abstract: We study the ${}^{1}S_0$ proton pairing gap in beta-equilibrated neutron star matter within the framework of chiral effective field theory. We focus on the role of three-body forces, which strongly modify the effective proton-proton spin-singlet interaction in dense matter. We find that three-body forces generically reduce both the size of the pairing gap and the maximum density at which proton pairing may occur. The pairing gap is computed within BCS theory, and model uncertainties are estimated by varying the nuclear potential and the choice of single-particle spectrum in the gap equation. We find that a second-order perturbative treatment of the single-particle spectrum suppresses the proton ${}^{1}S_0$ pairing gap relative to the use of a free spectrum. We estimate the critical temperature for the onset of proton superconductivity to be $T_c = (3.7 - 6.0)\times 10^{9} $ K, which is consistent with previous theoretical results in the literature and marginally within the range deduced from a recent Bayesian analysis of neutron star cooling observations.


title: Observation of "broad" d-wave Feshbach resonances with a triplet structure

abstract: High partial-wave ($l\ge2$) Feshbach resonance (FR) in an ultracold mixture of $^{85}$Rb-$^{87}$Rb atoms is investigated experimentally aided by a partial-wave insensitive analytic multichannel quantum-defect theory (MQDT). Two "broad" resonances from coupling between d-waves in both the open and closed channels are observed and characterized. One of them shows a fully resolved triplet structure with splitting ratio well explained by the perturbation to the closed channel due to interatomic spin-spin interaction. These tunable "broad" d-wave resonances, especially the one in the lowest-energy open channel, could find important applications in simulating d-wave coupling dominated many-body systems. In addition, we find that there is generally a time and temperature requirement, associated with tunneling through the angular momentum barrier, to establish and observe resonant coupling in nonzero partial waves.


title: On Stein's Identity and Near-Optimal Estimation in High-dimensional Index Models

abstract: We consider estimating the parametric components of semi-parametric multiple index models in a high-dimensional and non-Gaussian setting. Such models form a rich class of non-linear models with applications to signal processing, machine learning and statistics. Our estimators leverage the score function based first and second-order Stein's identities and do not require the covariates to satisfy Gaussian or elliptical symmetry assumptions common in the literature. Moreover, to handle score functions and responses that are heavy-tailed, our estimators are constructed via carefully thresholding their empirical counterparts. We show that our estimator achieves near-optimal statistical rate of convergence in several settings. We supplement our theoretical results via simulation experiments that confirm the theory.


title: A practical theoretical formalism for atomic multielectron processes: direct multiple ionization by a single Auger decay or by impact of a single electron or photon

abstract: Multiple electron processes occur widely in atoms, molecules, clusters, and condensed matters when they are interacting with energetic particles or intense laser fields. In this work, a unified and accurate theoretical formalism is proposed on the direct multielectron processes of atoms including the multiple Auger decay and multiple ionization by an impact of an energetic electron or a photon based on the atomic collision theory described by a correlated many-body Green's function. Such a practical treatment is made possible due to different coherence features of the particles (matter waves) in the initial and final states. We first explain how the coherence characteristics of the ejected continuum electrons is largely destructed, by taking the electron impact direct double ionization process as an example. This process is completely different from the single ionization where the complete interference can be maintained. The detailed expressions are obtained for the energy correlations among the continuum electrons and energy resolved differential and integral cross sections according to the separation of knock-out and shake-off mechanisms for the electron impact direct double ionization, direct double and triple Auger decay, and double and triple photoionization processes. The approach is applied to investigate the electron impact double ionization processes of C$^+$, N$^+$, and O$^+$, the direct double and triple Auger decay of the K-shell excited states of C$^{+}$ $1s2s^22p^2$ $^2D$ and $^2P$, and the double and triple photoionization of lithium. Comparisons with available experimental and theoretical results show that our proposed theoretical formalism is accurate and effective in treating the atomic multielectron processes.


title: Ultra-Dense HetNets Meet Big Data: Green Frameworks, Techniques, and Approaches

abstract: Ultra-dense heterogeneous networks (Ud-HetNets) have been put forward to improve the network capacity for next-generation wireless networks. However, counter to the 5G vision, ultra-dense deployment of networks would significantly increase energy consumption and thus decrease network energy efficiency suffering from the conventional worst-case network design philosophy. This problem becomes particularly severe when Ud-HetNets meet big data because of the traditional reactive request-transmit service mode. In view of these, this article first develops a big-data-aware artificial intelligent based framework for energy-efficient operations of Ud-HetNets. Based on the framework, we then identify four promising techniques, namely big data analysis, adaptive base station operation, proactive caching, and interference-aware resource allocation, to reduce energy cost on both large and small scales. We further develop a load-aware stochastic optimization approach to show the potential of our proposed framework and techniques in energy conservation. In a nutshell, we devote to constructing green Ud-HetNets of big data with the abilities of learning and inferring by improving the flexibility of control from worst-case to adaptive design and shifting the manner of services from reactive to proactive modes.


title: Pulsations of the High-Amplitude $\delta$ Scuti Star YZ Bootis

abstract: We present a study on the pulsations of the high-amplitude $\delta$ Scuti star YZ Boo based on photometric observations in Johnson $V$ and $R$ bands with both the Nanshan 1-m telescope of Xinjiang Astronomical Observatory (XAO) and the Xinglong 85-cm telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC). The Fourier analysis of the light curves reveals the fundamental radial mode and its five harmonics, with the fourth and the fifth newly detected. 39 new times of light maximum are determined from the light curves, and combined with those in the literature, we construct the $O - C$ diagram and derive a new ephemeris and the determination of a new value of the updated period 0.104091579(2). Besides, the $O - C$ diagram reveals an increasing period change rate of YZ Boo. Theoretical models are calculated and constrained with the observationally determined parameters of YZ Boo. The mass and age of YZ Boo are hence derived as $M$ = 1.61 $\pm$ 0.05 $M_\odot$ and $age$ = (1.44 $\pm$ 0.14) $\times$ 10$^9$ yr, respectively. With both the frequency of the fundamental radial mode and the period change rate, YZ Boo is located at the post-main-sequence stage.


title: Collective Resonances in Nanoparticle Oligomers

abstract: The study of nanostructured artificial media for optics has expanded rapidly over the last few decades, coupled with improvements of fabrication technology that have enabled investigation of previously unrealisable optical scattering systems. Such development is complemented by renewed impetus to understand the physics of optical scattering from complex subwavelength geometry and nanoparticle systems. Here I investigate specifically the optical properties of closely packed arrangements of nanoparticles, known as nanoparticle oligomers, which provide an intuitive platform for analytical and numerical study on the formation and interplay of collective resonances. I consider both plasmonic nanoparticles, and also high-refractive-index dielectric nanoparticles that support Mie-type electric and magnetic dipole resonances. Specific outcomes of this study are listed as follows. (i) A new model is presented for optical Fano resonances, which is based on interference between nonorthogonal eigenmodes of the associated scattering object. This is demonstrated to correctly describe Fano resonances in both plasmonic and high-refractive-index dielectric nanoparticle oligomers. (ii) Polarisation-independent scattering and absorption losses are shown to be enforced by n-fold discrete rotational symmetry, Cn (n \geq 3), and the reciprocal degeneracy of eigenmodes. (iii) A form of circular dichroism that emerges from interaction of nonorthogonal resonances is presented, and specifically impacts the ratio of radiative scattering loss to dissipative absorption loss experienced by reciprocal plane waves. The sequence of general theoretical conclusions (i)-(iii) serve to build the understanding of optical scattering from nanoparticle systems while removing existing ambiguities.


title: TuringMobile: A Turing Machine of Oblivious Mobile Robots with Limited Visibility and its Applications

abstract: In this paper we investigate the computational power of a set of mobile robots with limited visibility. At each iteration, a robot takes a snapshot of its surroundings, uses the snapshot to compute a destination point, and it moves toward its destination. Each robot is punctiform and memoryless, it operates in $\mathbb{R}^m$, it has a local reference system independent of the other robots' ones, and is activated asynchronously by an adversarial scheduler. Moreover, robots are non-rigid, in that they may be stopped by the scheduler at each move before reaching their destination (but are guaranteed to travel at least a fixed unknown distance before being stopped). We show that despite these strong limitations, it is possible to arrange $3m+3k$ of these weak entities in $\mathbb{R}^m$ to simulate the behavior of a stronger robot that is rigid (i.e., it always reaches its destination) and is endowed with $k$ registers of persistent memory, each of which can store a real number. We call this arrangement a TuringMobile. In its simplest form, a TuringMobile consisting of only three robots can travel in the plane and store and update a single real number. We also prove that this task is impossible with fewer than three robots. Among the applications of the TuringMobile, we focused on Near-Gathering (all robots have to gather in a small-enough disk) and Pattern Formation (of which Gathering is a special case) with limited visibility. Interestingly, our investigation implies that both problems are solvable in Euclidean spaces of any dimension, even if the visibility graph of the robots is initially disconnected, provided that a small amount of these robots are arranged to form a TuringMobile. In the special case of the plane, a basic TuringMobile of only three robots is sufficient.


title: Valued fields, Metastable groups

abstract: We introduce a class of theories called metastable, including the theory of algebraically closed valued fields (ACVF) as a motivating example. The key local notion is that of definable types dominated by their stable part. A theory is metastable (over a sort $\Gamma$) if every type over a sufficiently rich base structure can be viewed as part of a $\Gamma$-parametrized family of stably dominated types. We initiate a study of definable groups in metastable theories of finite rank. Groups with a stably dominated generic type are shown to have a canonical stable quotient. Abelian groups are shown to be decomposable into a part coming from $\Gamma$, and a definable direct limit system of groups with stably dominated generic. In the case of ACVF, among definable subgroups of affine algebraic groups, we characterize the groups with stably dominated generics in terms of group schemes over the valuation ring. Finally, we classify all fields definable in ACVF.


title: A Deep Learning Model for Traffic Flow State Classification Based on Smart Phone Sensor Data

abstract: This study proposes a Deep Belief Network model to classify traffic flow states. The model is capable of processing massive, high-density, and noise-contaminated data sets generated from smartphone sensors. The statistical features of Vehicle acceleration, angular acceleration, and GPS speed data, recorded by smartphone software, are analyzed, and then used as input for traffic flow state classification. Data from a five-day experiment is used to train and test the proposed model. A total of 747,856 sets of data are generated and used for both traffic flow states classification and sensitivity analysis of input variables. The result shows that the proposed Deep Belief Network model is superior to traditional machine learning methods in both classification performance and computational efficiency.


title: Three dimensional free-surface flow over arbitrary bottom topography

abstract: We consider steady nonlinear free surface flow past an arbitrary bottom topography in three dimensions, concentrating on the shape of the wave pattern that forms on the surface of the fluid. Assuming ideal fluid flow, the problem is formulated using a boundary integral method and discretised to produce a nonlinear system of algebraic equations. The Jacobian of this system is dense due to integrals being evaluated over the entire free surface. To overcome the computational difficulty and large memory requirements, a Jacobian-free Newton Krylov (JFNK) method is utilised. Using a block-banded approximation of the Jacobian from the linearised system as a preconditioner for the JFNK scheme, we find significant reductions in computational time and memory required for generating numerical solutions. These improvements also allow for a larger number of mesh points over the free surface and the bottom topography. We present a range of numerical solutions for both subcritical and supercritical regimes, and for a variety of bottom configurations. We discuss nonlinear features of the wave patterns as well as their relationship to ship wakes.


title: Probing the dark matter-electron interactions via hydrogen-atmosphere pulsating white dwarfs

abstract: In this work, we propose a novel scenario to probe the interactions between dark matter (DM) particles and electrons, via hydrogen-atmosphere pulsating white dwarfs (DAVs) in globular clusters. In this special configuration, the DM particles, which are predominantly captured by multiple scattering with the electrons in a DAV, would annihilate by pairs and provide extra energy source to the DAV. This mechanism slows down the natural cooling evolution which can be presented by the period variation rates of pulsation modes. The differences between the secular rates predicted by the precise asteroseismology and the secular rates obtained from observation can reveal the DM-electron interactions. An important observable has been proposed and corresponding estimations have been made. According to the estimation, if this scenario could be implemented in the near future, the potential sensitivity on $m_{\chi}$ (DM particle's mass) and $\sigma_{\chi,e}$ (elastic scattering cross section between DM and electron) could be hopefully extended to a region $5 \mathrm{GeV} \lesssim m_{\chi} \lesssim 10^{4} \mathrm{GeV}$ and $\sigma_{\chi,e} \gtrsim 10^{-40} \mathrm{cm}^{2}$. Combining with indirect DM detection results, this could give us a cross check on the existence of such leptonphilic DM particles to some extent.