title: 3DOF Pedestrian Trajectory Prediction Learned from Long-Term Autonomous Mobile Robot Deployment Data

abstract: This paper presents a novel 3DOF pedestrian trajectory prediction approach for autonomous mobile service robots. While most previously reported methods are based on learning of 2D positions in monocular camera images, our approach uses range-finder sensors to learn and predict 3DOF pose trajectories (i.e. 2D position plus 1D rotation within the world coordinate system). Our approach, T-Pose-LSTM (Temporal 3DOF-Pose Long-Short-Term Memory), is trained using long-term data from real-world robot deployments and aims to learn context-dependent (environment- and time-specific) human activities. Our approach incorporates long-term temporal information (i.e. date and time) with short-term pose observations as input. A sequence-to-sequence LSTM encoder-decoder is trained, which encodes observations into LSTM and then decodes as predictions. For deployment, it can perform on-the-fly prediction in real-time. Instead of using manually annotated data, we rely on a robust human detection, tracking and SLAM system, providing us with examples in a global coordinate system. We validate the approach using more than 15K pedestrian trajectories recorded in a care home environment over a period of three months. The experiment shows that the proposed T-Pose-LSTM model advances the state-of-the-art 2D-based method for human trajectory prediction in long-term mobile robot deployments.


title: A Challenge to Identify an Optical Counterpart of the Gravitational Wave Event GW151226 with Hyper Suprime-Cam

abstract: We present the results of the detailed analysis of an optical imaging survey conducted using the Subaru / Hyper Suprime-Cam (HSC), which aims to identify an optical counterpart to the gravitational wave event GW151226. In half a night, the $i$- and $z$-band imaging survey by HSC covers 63.5deg$^2$ of the error region, which contains about 7\% of the LIGO localization probability, and the same field is observed in three different epochs. The detectable magnitude of the candidates in a differenced image is evaluated as $i \sim 23.2$ mag for the requirement of at least two 5$\sigma$ detections, and 1744 candidates are discovered. Assuming a kilonova as an optical counterpart, we compared the optical properties of the candidates with model predictions. A red and rapidly declining light curve condition enables the discrimination of a kilonova from other transients, and a small number of candidates satisfy this condition. The presence of stellar-like counterparts in the reference frame suggests that the surviving candidates are likely to be flare stars. The fact that most of those candidates are in galactic plane, $|b|<5^{\circ}$, supports this interpretation. We also checked whether the candidates are associated with the nearby GLADE galaxies, which reduces the number of contaminants even with a looser color cut. When a better probability map (with localization accuracy of $\sim50{\rm deg}^2$) is available, kilonova searches of up to approximately $200$ Mpc will become feasible by conducting immediate follow-up observations with an interval of 3--6 days.


title: Delay Embedding of Periodic Orbits Using a Fixed Observation Function

abstract: Delay coordinates are a widely used technique to pass from observations of a dynamical system to a representation of the dynamical system as an embedding in Euclidean space. Current proofs show that delay coordinates of a given dynamical system result in embeddings generically over a space of observations (Sauer, Yorke, Casdagli, J. Stat. Phys., vol. 65 (1991), p. 579-616). Motivated by applications of the embedding theory, we consider the situation where the observation function is fixed. For example, the observation function may simply be some fixed coordinate of the state vector. For a fixed observation function (any nonzero linear combination of coordinates) and for the special case of periodic solutions, we prove that delay coordinates result in an embedding generically over the space of flows in the $C^{r}$ topology with $r\geq2$.


title: Reverberation Mapping of PG 0934+013 with the Southern African Large Telescope

abstract: We present the variability and time lag measurements of PG 0934+013 based on a photometric and spectroscopic monitoring campaign over a two year period. We obtained 46 epochs of data from the spectroscopic campaign, which was carried out using the Southern African Large Telescope with $\sim$1 week cadence over two sets of 4 month-long observing period, while we obtained 80 epochs of \textit{B}-band imaging data using a few 1-m class telescopes. Due to the seven month gap between the two observing periods, we separately measured the time lags of broad emission lines including H$\beta$, by comparing the emission line light curve with the \textit{B}-band continuum light curve using the cross-correlation function techniques. We determined the H$\beta$ lag, $\tau_{\rm cent} = 8.46^{+2.08}_{-2.14}$ days in the observed-frame based on Year 2 data, while the time lag from Year 1 data was not reliably determined. Using the rms spectrum of Year 2 data, we measured the \Hb\ line dispersion \sigmaline = 668 $\pm$ 44 \kms\ after correcting for the spectral resolution. Adopting a virial factor f = 4.47 from Woo et al. 2015, we determined the black hole mass M$_{BH}$ = $3.13 ^{+0.91} _{-0.93} \times 10^{6}$ \msun based on the \Hb\ time lag and velocity.


title: A Cheeger-type exponential bound for the number of triangulated manifolds

abstract: In terms of the number of triangles, it is known that there are more than exponentially many triangulations of surfaces, but only exponentially many triangulations of surfaces with bounded genus. In this paper we provide a first geometric extension of this result to higher dimensions. We show that in terms of the number of facets, there are only exponentially many geometric triangulations of space forms with bounded geometry in the sense of Cheeger (curvature and volume bounded below, and diameter bounded above). This establishes a combinatorial version of Cheeger's finiteness theorem. Further consequences of our work are: (1) There are exponentially many geometric triangulations of $S^d$. (2) There are exponentially many convex triangulations of the d-ball.


title: Electrical Control of Linear Dichroism in Black Phosphorus from the Visible to Mid-Infrared

abstract: The incorporation of electrically tunable materials into photonic structures such as waveguides and metasurfaces enables dynamic control of light propagation by an applied potential. While many materials have been shown to exhibit electrically tunable permittivity and dispersion, including transparent conducting oxides (TCOs) and III-V semiconductors and quantum wells, these materials are all optically isotropic in the propagation plane. In this work, we report the first known example of electrically tunable linear dichroism, observed here in few-layer black phosphorus (BP), which is a promising candidate for multi-functional, broadband, tunable photonic elements. We measure active modulation of the linear dichroism from the mid-infrared to visible frequency range, which is driven by anisotropic quantum-confined Stark and Burstein-Moss effects, and field-induced forbidden-to-allowed optical transitions. Moreover, we observe high BP absorption modulation strengths, approaching unity for certain thicknesses and photon energies.


title: Dense RGB-D semantic mapping with Pixel-Voxel neural network

abstract: For intelligent robotics applications, extending 3D mapping to 3D semantic mapping enables robots to, not only localize themselves with respect to the scene's geometrical features but also simultaneously understand the higher level meaning of the scene contexts. Most previous methods focus on geometric 3D reconstruction and scene understanding independently notwithstanding the fact that joint estimation can boost the accuracy of the semantic mapping. In this paper, a dense RGB-D semantic mapping system with a Pixel-Voxel network is proposed, which can perform dense 3D mapping while simultaneously recognizing and semantically labelling each point in the 3D map. The proposed Pixel-Voxel network obtains global context information by using PixelNet to exploit the RGB image and meanwhile, preserves accurate local shape information by using VoxelNet to exploit the corresponding 3D point cloud. Unlike the existing architecture that fuses score maps from different models with equal weights, we proposed a Softmax weighted fusion stack that adaptively learns the varying contributions of PixelNet and VoxelNet, and fuses the score maps of the two models according to their respective confidence levels. The proposed Pixel-Voxel network achieves the state-of-the-art semantic segmentation performance on the SUN RGB-D benchmark dataset. The runtime of the proposed system can be boosted to 11-12Hz, enabling near to real-time performance using an i7 8-cores PC with Titan X GPU.


title: Gravitational catalysis of merons in Einstein-Yang-Mills theory

abstract: We construct regular configurations of the Einstein-Yang-Mills theory in various dimensions. The gauge field is of meron-type: it is proportional to a pure gauge (with a suitable parameter $\lambda$ determined by the field equations). The corresponding smooth gauge transformation cannot be deformed continuously to the identity. In the three-dimensional case we consider the inclusion of a Chern-Simons term into the analysis, allowing $\lambda$ to be different from its usual value of $1/2$. In four dimensions, the gravitating meron is a smooth Euclidean wormhole interpolating between different vacua of the theory. In five and higher dimensions smooth meron-like configurations can also be constructed by considering warped products of the three-sphere and lower-dimensional Einstein manifolds. In all cases merons (which on flat spaces would be singular) become regular due to the coupling with general relativity. This effect is named "gravitational catalysis of merons".


title: Alternative methods to measure global polarization of $\Lambda$ hyperons

abstract: We propose alternative methods to measure the global polarization of $\Lambda$ hyperons. These methods involve event averages of proton's and $\Lambda$'s momenta in the lab frame. We carry out simulations using these methods and show that all of them work equivalently well in obtaining the global polarization of $\Lambda$ hyperons.


title: Improved Hardy and Rellich inequalities on nonreversible Finsler manifolds

abstract: In this paper, we study the sharp constants of quantitative Hardy and Rellich inequalities on nonreversible Finsler manifolds equipped with arbitrary measures. In particular, these inequalities can be globally refined by adding remainder terms like the Brezis-V\'azquez improvement, if Finsler manifolds are of strictly negative flag curvature, vanishing S-curvature and finite uniformity constant. Furthermore, these results remain valid when Finsler metrics are reversible.


title: Hot Spaghetti: Viscous Gravitational Collapse

abstract: We explore the fate of matter falling into a macroscopic Schwarzschild black hole for the simplified case of a radially collapsing thin spherical shell for which the back reaction of the geometry can be neglected. We treat the internal dynamics of the infalling matter in the framework of viscous relativistic hydrodynamics and calculate how the internal temperature of the collapsing matter evolves as it falls toward the Schwarzschild singularity. We find that viscous hydrodynamics fails when either, the dissipative radial pressure exceeds the thermal pressure and the total radial pressure becomes negative, or the time scale of variation of the tidal forces acting on the collapsing matter becomes shorter than the characteristic hydrodynamic response time.


title: Generic Newton polygon for exponential sums in $n$ variables with parallelotope base

abstract: Let $p$ be a prime number. Every $n$-variable polynomial $f(\underline x)$ over a finite field of characteristic $p$ defines an Artin--Schreier--Witt tower of varieties whose Galois group is isomorphic to $\mathbb{Z}_p$. Our goal of this paper is to study the Newton polygon of the $L$-function associated to a finite character of $\mathbb{Z}_p$ and a generic polynomial whose convex hull is an $n$-dimensional paralleltope $\Delta$. We denote this polygon by $\mathrm{GNP}(\Delta)$. We prove a lower bound of $\mathrm{GNP}(\Delta)$, which is called the improved Hodge polygon $\mathrm{IHP}(\Delta)$. We show that $\mathrm{IHP}(\Delta)$ lies above the usual Hodge polygon $\mathrm{HP}(\Delta)$ at certain infinitely many points, and when $p$ is larger than a fixed number determined by $\Delta$, it coincides with $\mathrm{GNP}(\Delta)$ at these points. As a corollary, we roughly determine the distribution of the slopes of $\mathrm{GNP}(\Delta)$.


title: Magnetic response from constant backgrounds to Coulomb sources

abstract: Magnetically uncharged, magnetic linear response of the vacuum filled with arbitrarily combined constant electric and magnetic fields to an imposed static electric charge is found within general nonlinear electrodynamics. When the electric charge is point-like and external fields are parallel, the response found may be interpreted as a field of two point-like magnetic charges of opposite polarity in one point. Coefficients characterizing the magnetic response and induced currents are specialized to Quantum Electrodynamics, where the nonlinearity is taken as that determined by the Heisenberg-Euler effective Lagrangian.


title: Effect of three-spin interaction on thermal entanglement in Heisenberg XXZ model

abstract: The effect of three-spin interaction k on thermal entanglement between alternate qubits is studied using pairwise concurrence C and energy level diagram. It is found that k breaks the symmetry about the effect of magnetic field h on C. It shifts a dip structure and gradually effaces a boot structure when |k|<|J| (J is spin exchange coupling). A peak with C maintains 1 appears and expands, elbowing the dip backwards when |k|>|J|. A sudden change of the concurrence occurs around |k|=|J|, h=-k. Similar conclusions about nearest-neighbor qubits are directly given.


title: Fermionic Matrix Product States and One-Dimensional Short-Range Entangled Phases with Anti-Unitary Symmetries

abstract: We extend the formalism of Matrix Product States (MPS) to describe one-dimensional gapped systems of fermions with both unitary and anti-unitary symmetries. Additionally, systems with orientation-reversing spatial symmetries are considered. The short-ranged entangled phases of such systems are classified by three invariants, which characterize the projective action of the symmetry on edge states. We give interpretations of these invariants as properties of states on the closed chain. The relationship between fermionic MPS systems at an RG fixed point and equivariant algebras is exploited to derive a group law for the stacking of fermionic phases protected by general fermionic symmetry groups.


title: Black hole hyperaccretion inflow-outflow model. I. long and ultra-long gamma-ray bursts

abstract: Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.


title: Outflows from black hole hyperaccretion systems: short and long-short gamma-ray bursts and "quasi-supernovae"

abstract: The detections of some long gamma-ray bursts (LGRBs) relevant to mergers of neutron star (NS)-NS or black hole (BH)-NS, as well as some short gamma-ray bursts (SGRBs) probably produced by collapsars, muddle the boundary of two categories of gamma-ray bursts (GRBs). In both cases, a plausible candidate of central engine is a BH surrounded by a hyperaccretion disc with strong outflows, launching relativistic jets driven by Blandford-Znajek mechanism. In the framework of compact binary mergers, we test the applicability of the BH hyperaccretion inflow-outflow model on powering observed GRBs. We find that, for a low outflow ratio, $\sim 50\%$, postmerger hyperaccretion processes could power not only all SGRBs but also most of LGRBs. Some LGRBs might do originate from merger events in the BH hyperaccretion scenario, at least on the energy requirement. Moreover, kilonovae might be produced by neutron-rich outflows, and their luminosities and timescales significantly depend on the outflow strengths. GRBs and their associated kilonovae are competitive with each other on the disc mass and total energy budgets. The stronger the outflow, the more similar the characteristics of kilonovae to supernovae (SNe). This kind of `nova' might be called `quasi-SN'.


title: Sub classes of Bi-Univalent Functions Defined by Salagean type $q-$ Difference Operator

abstract: In this paper, we introduce and investigate a new subclass of the function class $\Sigma$ of bi-univalent functions defined in the open unit disk, which are associated with the S\u{a}l\u{a}gean type $q-$ difference operator and satisfy some subordination conditions. Furthermore, we find estimates on the Taylor-Maclaurin coefficients $|a_2|$ and $|a_3|$ for functions in the new subclass introduced here. Several (known or new) consequences of the results are also pointed out. Further we obtain Fekete-Szeg$\ddot{o}$ inequality for the new function class.


title: On the Long-Range Exciton Transport in H-Aggregated Heterotriangulene Chains

abstract: Self-assembled aggregates of pigment molecules are potential building blocks for excitonic circuits that find their application in energy conversion and optical signal processing. Recent experimental studies of one-dimensional heterotriangulene supramolecular aggregates suggested that singlet excitons in these structures can propagate on several micron distances. We explore this possibility theoretically by combining electronic structure calculations with microscopic models for exciton transport. Detailed characterization of the structural disorder and exciton decoherence is provided. We argue that conventional exciton transport models give about an order of magnitude shorter estimates for the exciton propagation length which suggest that there are other possible explanations of the experimental results.


title: A Game-Theoretic Framework for Multi-Period-Multi-Company Demand Response Management in the Smart Grid

abstract: By utilizing tools from game theory, we develop a novel multi-period-multi-company demand response framework considering the interactions between companies (sellers of energy) and their consumers (buyers of energy). We model the interactions in terms of a Stackelberg game, where companies set their prices and consumers respond by choosing their demands. We show that the underlying game has a unique equilibrium at which the companies maximize their revenues while the consumers maximize their utilities subject to their local constraints. Closed-form expressions are provided for the optimal strategies of all players. Based on these solutions, a power allocation game has been formulated, which is shown to admit a unique pure-strategy Nash equilibrium, for which closed-form expressions are also provided. This equilibrium is found under the assumption that companies can freely allocate their power across the time horizon, but we also demonstrate that it is possible to relax this assumption. We further provide a fast distributed algorithm for the computation of all optimal strategies using only local information. We also study the effect of variations in the number of periods (subdivisions of the time horizon) and the number of consumers. As a consequence, we are able to find an appropriate company-to-consumer ratio when the number of consumers participating in demand response exceeds some threshold. Furthermore, we show, both analytically and numerically, that the multi-period scheme provides incentives for energy consumers to participate in demand response, compared to the single-period framework studied in the literature. In our framework, we provide a condition for the minimum budgets consumers need, and carry out case studies using real life data to demonstrate the benefits of the approach, which show potential savings of up to $30\%$ and equilibrium prices that have low volatility.


title: T-matrix evaluation of three-dimensional acoustic radiation forces on nonspherical objects in Bessel beams with arbitrary order and location

abstract: Acoustic radiation forces (ARFs) induced by a single Bessel beam with arbitrary order and location on a nonspherical shape are studied using the T-matrix method (TMM) in three dimensions. Based on the radiation stress tensor approach and the multipole expansion method for the arbitrary Bessel beam, the ARF expressions are derived in terms of the incident and scattered beam shape coefficients independently with the corresponding homemade code packages. Several numerical experiments are conducted to verify the versatility of the TMM. The axial acoustic radiation forces (ARFs) of several typical shapes are considered in the analysis with the emphasis on the axial ARF reversal and the corresponding physical mechanism. This study may guide the experimental set-up to find negative axial ARFs quickly and effectively based on the predicted parameters with TMM. Relatively elongated shapes may be helpful for pulling forces in Bessel beams. Furthermore, the lateral ARFs for both convex and concave nonspherical shapes are also investigated with different topological charges, cone angles and offsets of the particle centroid to the beam axis in a broadband frequency regime. A brief theoretical derivation of the incident beam shape coefficients for the standing Bessel beams is also given. The present work could help to design the acoustic tweezers numerical toolbox which provides an acoustical alternative to the optical tweezers toolbox.


title: Excitation of high frequency voices from intermediate-mass-ratio inspirals with large eccentricity

abstract: The coalescence of a stellar-mass compact object together with an intermediate-mass black hole, also known as intermediate-mass-ratio inspiral, is usually not expected to be a viable gravitational wave source for the current ground-based gravitational wave detectors, due to the generally lower frequency of such source. In this paper, we adopt the effective-one-body formalism as the equation of motion, and obtain the accurately calculated gravitational waveforms by solving the Teukolsky equation in frequency-domain. We point out that high frequency modes of gravitational waves can be excited by large eccentricities of intermediate-mass-ratio inspirals. These high frequency modes can extend to more than 10 Hz, and enter the designed sensitive band of Advanced LIGO and Advanced Virgo. We propose that such kind of highly eccentric intermediate-mass-ratio inspirals could be feasible sources and potentially observable by the ground-based gravitational wave detectors, like the Advanced LIGO and Advanced Virgo.


title: Understanding caustic crossings in giant arcs: characteristic scales, event rates, and constraints on compact dark matter

abstract: The recent discovery of fast transient events near critical curves of massive galaxy clusters, which are interpreted as highly magnified individual stars in giant arcs due to caustic crossing, opens up the possibility of using such microlensing events to constrain a range of dark matter models such as primordial black holes and scalar field dark matter. Based on a simple analytic model, we study lensing properties of a point mass lens embedded in a high magnification region, and derive the dependence of the peak brightness, microlensing time scales, and event rates on the mass of the point mass lens as well as the radius of a source star that is magnified. We find that the lens mass and source radius of the first event MACS J1149 Lensed Star 1 (LS1) are constrained, with the lens mass range of $0.1~M_\odot \lesssim M \lesssim 4\times 10^3M_\odot$ and the source radius range of $40~R_\odot \lesssim R \lesssim 260~R_\odot$. In the most plausible case with $M\approx 0.3~M_\odot$ and $R\approx 180~R_\odot$, the source star should have been magnified by a factor of $\approx 4300$ at the peak. The derived lens properties are fully consistent with the interpretation that MACS J1149 LS1 is a microlensing event produced by a star that contributes to the intra-cluster light. We argue that compact dark matter models with high fractional mass densities for the mass range $10^{-5}M_\odot \lesssim M\lesssim 10^2M_\odot$ are inconsistent with the observation of MACS J1149 LS1 because such models predict too low magnifications. Our work demonstrates a potential use of caustic crossing events in giant arcs to constrain compact dark matter.


title: Experiment to Form and Characterize a Section of a Spherically Imploding Plasma Liner

abstract: We describe an experiment to form and characterize a section of a spherically imploding plasma liner by merging six supersonic plasma jets that are launched by newly designed contoured-gap coaxial plasma guns. This experiment is a prelude to forming a fully spherical imploding plasma liner using many dozens of plasma guns, as a standoff driver for plasma-jet-driven magneto-inertial fusion. The objectives of the six-jet experiments are to assess the evolution and scalings of liner Mach number and uniformity, which are important metrics for spherically imploding plasma liners to compress magnetized target plasmas to fusion conditions. This paper describes the design of the coaxial plasma guns, experimental characterization of the plasma jets, six-jet experimental setup and diagnostics, initial diagnostic data from three- and six-jet experiments, and the high-level objectives of associated numerical modeling.


title: Precision measurement of the $e^{+}e^{-}~\rightarrow~\Lambda_{c}^{+} \bar{\Lambda}_{c}^{-}$ cross section near threshold

abstract: The Born cross section of the $e^{+}e^{-}~\rightarrow~\Lambda_{c}^{+} \bar{\Lambda}_{c}^{-}$ process is measured with unprecedented precision using data collected with the BESIII detector at $\sqrt{s}=4574.5$, $45580.0$, $4590.0$ and $4599.5$ $\mathrm{MeV}$. The non-zero cross section near the $\Lambda_{c}^{+} \bar{\Lambda}_{c}^{-}$ production threshold is discerned. At center-of-mass energies $\sqrt{s}=4574.5$ and $4599.5 $ $\mathrm{MeV}$, the higher statistics data enable us to measure the $\Lambda_{c}$ polar angle distributions. From these, the ratio between the $\Lambda_{c}$ electric and magnetic form factors ($|G_{E}/G_{M}|$) is measured for the first time. They are found to be $1.14\pm0.14\pm0.07$ and $1.23\pm0.05\pm0.03$ respectively, where the first uncertainties are statistical and the second are systematic.


title: Two-Way Energy Trading and Online Planning for Fifth-Generation Communications with Renewables

abstract: Future fifth-generation (5G) cellular networks, equipped with energy harvesting devices, are uniquely positioned to closely interoperate with smart grid. New interoperable functionalities are discussed in stochastic two-way energy trading and online planning to improve efficiency and productivity. Challenges lie in the unavailability of a-priori knowledge on future wireless channels, energy pricing and harvesting. Lyapunov optimization techniques are utilized to address the challenges and stochastically optimize energy trading and planning. Particularly, it is able to decouple the optimization of energy trading and planning during individual time slots, hence eliminating the need for joint optimization across a large number of slots.


title: Boson peak dynamics of glassy glucose studied by integrated terahertz-band spectroscopy

abstract: We performed terahertz time-domain spectroscopy, low-frequency Raman scattering, and Brillouin light scattering on vitreous glucose to investigate the boson peak (BP) dynamics. In the spectra of {\alpha}({\nu})/{\nu}2 [{\alpha}({\nu}) is the absorption coefficient], the BP is clearly observed around 1.1 THz. Correspondingly, the complex dielectric constant spectra show a universal resonancelike behavior only below the BP frequency. As an analytical scheme, we propose the relative light-vibration coupling coefficient (RCC), which is obtainable from the combination of the far-infrared and Raman spectra. The RCC reveals that the infrared light-vibration coupling coefficient CIR({\nu}) of the vitreous glucose behaves linearly on frequency which deviates from Taraskin's model of CIR({\nu}) = A + B{\nu}2 [S. N. Taraskin et al., Phys. Rev. Lett. 97, 055504 (2006)]. The linearity of CIR({\nu}) might require modification of the second term of the model. The measured transverse sound velocity shows an apparent discontinuity with the flattened mode observed in the inelastic neutron scattering study [N. Violini et al., Phys. Rev. B 85, 134204 (2012)] and suggests a coupling between the transverse acoustic and flattened modes.


title: A Multi-Resolution Model for Non-Gaussian Random Fields on a Sphere with Application to Ionospheric Electrostatic Potentials

abstract: Gaussian random fields have been one of the most popular tools for analyzing spatial data. However, many geophysical and environmental processes often display non-Gaussian characteristics. In this paper, we propose a new class of spatial models for non-Gaussian random fields on a sphere based on a multi-resolution analysis. Using a special wavelet frame, named spherical needlets, as building blocks, the proposed model is constructed in the form of a sparse random effects model. The spatial localization of needlets, together with carefully chosen random coefficients, ensure the model to be non-Gaussian and isotropic. The model can also be expanded to include a spatially varying variance profile. The special formulation of the model enables us to develop efficient estimation and prediction procedures, in which an adaptive MCMC algorithm is used. We investigate the accuracy of parameter estimation of the proposed model, and compare its predictive performance with that of two Gaussian models by extensive numerical experiments. Practical utility of the proposed model is demonstrated through an application of the methodology to a data set of high-latitude ionospheric electrostatic potentials, generated from the LFM-MIX model of the magnetosphere-ionosphere system.


title: Large time behavior of strong solutions to the compressible magnetohydrodynamic system in the critical $L^{p}$ framework

abstract: In this paper, we are concerned with the compressible viscous magnetohydrodynamic (MHD) system and investigate the large time behavior of strong solutions near constant equilibrium. In the eighties, Umeda, Kawashima and Shizuta initiated the dissipative mechanism for a rather general class of symmetric hyperbolic-parabolic systems. From the point of view of dissipativity, Kawashima in his doctoral dissertation established the optimal time-decay estimates of $L^{q}$-$L^{2}(q\in [1,2)$) type for solutions to the MHD system. Here, we shall improve Kawashima's efforts and give more precise description for the large time asymptotic behavior of solutions, not only in extra Lebesgue spaces but also in a full family of Besov norms with the negative regularity index. Precisely, we show that the $L^{p}$ norm (the slightly stronger $\dot{B}^{0}_{p,1}$ one in fact) of global solutions with the critical regularity, decays like $t^{-d(\frac{1}{p}-\frac{1}{4})}$ for $t\rightarrow\infty$. In particular, taking $p=2$ and $d=3$ goes back to the classical time decay $t^{-\frac{3}{4}}$. We derive new estimates which are used to deal with the strong coupling between the magnetic field and fluid dynamics.


title: Theoretical Investigation of the Magnetoelectric Properties of Bi2NiTiO6

abstract: We report the first principle investigations on the structural, electronic, magnetic and ferroelectric properties of a Pb free double perovskite multiferroic Bi2NiTiO6 using density functional theory within the general gradient approximation (GGA) and GGA+U method. Our results show that Bi2NiTiO6 will be an insulator with G-type magnetic ordering in its ground state with Ni2+ in a high spin state and a spin moment of 1.74\mu_B. The paraelectric phase stabilizes in nonmagnetic state with Ni2+ in low spin configuration showing that spin state transition plays an important role in strong magnetoelectric coupling in Bi2NiTiO6. The bonding characteristics of the constituents are analyzed with the help of partial density of states and Born effective charges. The presence of Ti ions at Ni sites suppresses the disproportionation observed in case of BiNiO3 and results in a noncentrosymmetric crystal structure. The coexistence of Bi 6s lone pair and Ti4+ d0 ions which brings covalency produces a polarization of 32 \muCcm-2.


title: Ultra-Wideband Aided Fast Localization and Mapping System

abstract: This paper proposes an ultra-wideband (UWB) aided localization and mapping system that leverages on inertial sensor and depth camera. Inspired by the fact that visual odometry (VO) system, regardless of its accuracy in the short term, still faces challenges with accumulated errors in the long run or under unfavourable environments, the UWB ranging measurements are fused to remove the visual drift and improve the robustness. A general framework is developed which consists of three parallel threads, two of which carry out the visual-inertial odometry (VIO) and UWB localization respectively. The other mapping thread integrates visual tracking constraints into a pose graph with the proposed smooth and virtual range constraints, such that an optimization is performed to provide robust trajectory estimation. Experiments show that the proposed system is able to create dense drift-free maps in real-time even running on an ultra-low power processor in featureless environments.


title: Enhanced Linear-array Photoacoustic Beamforming using Modified Coherence Factor

abstract: Photoacoustic imaging (PAI) is a promising medical imaging modality providing the spatial resolution of ultrasound (US) imaging and the contrast of pure optical imaging. For linear-array PAI, a beamformer has to be used as the reconstruction algorithm. Delay-and-sum (DAS) is the most prevalent beamforming algorithm in PAI. However, using DAS beamformer leads to low resolution images along with significant effects of the off-axis signals. Coherence factor (CF) is a weighting method in which each pixel of the reconstructed image is weighted, based on the spatial spectrum of the aperture, to improve the contrast. In this paper, it has been shown that the numerator of the formula of CF contains a DAS algebra, and it was proposed to use the delay-multiply-and-sum (DMAS) beamformer instead of the available DAS on the numerator. The proposed weighting technique, modified CF (MCF), has been evaluated numerically and experimentally compared to CF, and it was shown that MCF leads to lower sidelobes and better detectable targets. The quantitative results of the experiment (using wire targets) show that MCF leads to for about 45% and 40% improvement, in comparison with CF, in the terms of signal-to-noise ratio and full-width-half-maximum, respectively.


title: On a reduction of the generalized Darboux-Halphen system

abstract: The equations for the general Darboux-Halphen system obtained as a reduction of the self-dual Yang-Mills can be transformed to a third-order system which resembles the classical Darboux-Halphen system with a common additive terms. It is shown that the transformed system can be further reduced to a constrained non-autonomous, non-homogeneous dynamical system. This dynamical system becomes homogeneous for the classical Darboux-Halphen case, and was studied in the context of self-dual Einstein's equations for Bianchi IX metrics. A Lax pair and Hamiltonian for this reduced system is derived and the solutions for the system are prescribed in terms of hypergeometric functions.


title: Search for low mass vector resonances decaying into quark-antiquark pairs in proton-proton collisions at sqrt(s) = 13 TeV

abstract: A search for narrow vector resonances decaying into quark-antiquark pairs is presented. The analysis is based on data collected in proton-proton collisions at sqrt(s) = 13 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 inverse femtobarns. The hypothetical resonance is produced with sufficiently high transverse momentum that its decay products are merged into a single jet with two-prong substructure. A signal would be identified as a peak over a smoothly falling background in the distribution of the invariant mass of the jet, using novel jet substructure techniques. No evidence for such a resonance is observed within the mass range of 50-300 GeV. Upper limits at 95% confidence level are set on the production cross section, and presented in a mass-coupling parameter space. The limits further constrain simplified models of dark matter production involving a mediator interacting between quarks and dark matter particles through a vector or axial-vector current. In the framework of these models, the results are the most sensitive to date, extending for the first time the search region to masses below 100 GeV.


title: Derivations with Leibniz defect

abstract: The non-Leibniz formalism is introduced in this article. The formalism is based on the generalized differentiation operator (kappa-operator) with a non-zero Leibniz defect. The Leibniz defect of the introduced operator linearly depends on one scaling parameter. In a special case, if the Leibniz defect vanishes, the generalized differentiation operator reduces to the common differentiation operator. The kappa-operator allows the formulation of the variational principles and corresponding Lagrange and Hamiltonian equations. The solutions of some generalized dynamical equations are provided closed form.With a positive Leibniz defect the amplitude of free vibration remains constant with time with the fading frequency (<<red shift>>). The negative Leibniz defect leads the opposite behavior, demonstrating the growing frequency (<<blue shift>>). However, the Hamiltonian remains constant in time in both cases. Thus the introduction of non-zero Leibniz defect leads to an alternative mathematical description of the conservative systems.


title: A Product Integration Method for the Approximation of the Early Exercise Boundary in the American Option Pricing Problem

abstract: In this paper, an integral equation representation for the early exercise boundary of an American option contract is considered. Thus far, a number of different techniques have been proposed in the literature to obtain a variety of integral equation forms for the early exercise boundary, all starting from the Black-Scholes partial differential equation. We first present a coherent categorization of exiting integral equation methodologies in the American option pricing literature. In the reminder and based on the fact that the early exercise boundary satisfies a fully nonlinear weakly singular non-standard Volterra integral equation, we propose a product integration approach based on linear barycentric rational interpolation to solve the problem. The price of the option will then be computed using the obtained approximation of the early exercise boundary and a barycentric rational quadrature. The convergence of the approximation scheme will also be analyzed. Finally, some numerical experiments based on the introduced method are presented and compared to some exiting approaches.


title: Accelerated Directional Search with non-Euclidean prox-structure

abstract: In the paper we propose an accelerated directional search method with non-euclidian prox-structure. We consider convex unconstraint optimization problem in $\mathbb{R}^n$. For simplicity we start from the zero point. We expect in advance that 1-norm of the solution is close enough to its 2-norm. In this case the standard accelerated Nesterov's directional search method can be improved. In the paper we show how to make Nesterov's method $n$-times faster (up to a $\log n$-factor) in this case. The basic idea is to use linear coupling, proposed by Allen-Zhu & Orecchia in 2014, and to make Grad-step in 2-norm, but Mirr-step in 1-norm. We show that for constrained optimization problems this approach stable upon an obstacle.


title: Automated Program Analysis for Novice Programmers

abstract: This paper describes how to adapt a static code analyzer to help novice programmers. Current analyzers have been built to give feedback to experienced programmers who build new applications or systems. The type of feedback and the type of analysis of these tools focusses on mistakes that are relevant within that context, and help with debugging the system. When teaching novice programmers this type of advice is often not particularly useful. It would be instead more useful to use these techniques to find problem in the understanding of students of important programming concepts. This paper first explores in what respect static analyzers support the learning and teaching of programming can be implemented based on existing static analysis technology. It presents an extension to static analyzer PMD was made so that feedback messages appear which are easier to understand for novice programmers. To answer the question if these techniques are able to find conceptual mistakes that are characteristic for novice programmers make, we ran it over a number of student projects, and compared these results with publicly available mature software projects.


title: Speaker Role Contextual Modeling for Language Understanding and Dialogue Policy Learning

abstract: Language understanding (LU) and dialogue policy learning are two essential components in conversational systems. Human-human dialogues are not well-controlled and often random and unpredictable due to their own goals and speaking habits. This paper proposes a role-based contextual model to consider different speaker roles independently based on the various speaking patterns in the multi-turn dialogues. The experiments on the benchmark dataset show that the proposed role-based model successfully learns role-specific behavioral patterns for contextual encoding and then significantly improves language understanding and dialogue policy learning tasks.


title: Dynamic Time-Aware Attention to Speaker Roles and Contexts for Spoken Language Understanding

abstract: Spoken language understanding (SLU) is an essential component in conversational systems. Most SLU component treats each utterance independently, and then the following components aggregate the multi-turn information in the separate phases. In order to avoid error propagation and effectively utilize contexts, prior work leveraged history for contextual SLU. However, the previous model only paid attention to the content in history utterances without considering their temporal information and speaker roles. In the dialogues, the most recent utterances should be more important than the least recent ones. Furthermore, users usually pay attention to 1) self history for reasoning and 2) others' utterances for listening, the speaker of the utterances may provides informative cues to help understanding. Therefore, this paper proposes an attention-based network that additionally leverages temporal information and speaker role for better SLU, where the attention to contexts and speaker roles can be automatically learned in an end-to-end manner. The experiments on the benchmark Dialogue State Tracking Challenge 4 (DSTC4) dataset show that the time-aware dynamic role attention networks significantly improve the understanding performance.


title: PCANet-II: When PCANet Meets the Second Order Pooling

abstract: PCANet, as one noticeable shallow network, employs the histogram representation for feature pooling. However, there are three main problems about this kind of pooling method. First, the histogram-based pooling method binarizes the feature maps and leads to inevitable discriminative information loss. Second, it is difficult to effectively combine other visual cues into a compact representation, because the simple concatenation of various visual cues leads to feature representation inefficiency. Third, the dimensionality of histogram-based output grows exponentially with the number of feature maps used. In order to overcome these problems, we propose a novel shallow network model, named as PCANet-II. Compared with the histogram-based output, the second order pooling not only provides more discriminative information by preserving both the magnitude and sign of convolutional responses, but also dramatically reduces the size of output features. Thus we combine the second order statistical pooling method with the shallow network, i.e., PCANet. Moreover, it is easy to combine other discriminative and robust cues by using the second order pooling. So we introduce the binary feature difference encoding scheme into our PCANet-II to further improve robustness. Experiments demonstrate the effectiveness and robustness of our proposed PCANet-II method.


title: Impulse-induced localized nonlinear modes in an electrical lattice

abstract: Intrinsic localized modes, also called discrete breathers, can exist under certain conditions in one-dimensional nonlinear electrical lattices driven by external harmonic excitations. In this work, we have studied experimentally the efectiveness of generic periodic excitations of variable waveform at generating discrete breathers in such lattices. We have found that this generation phenomenon is optimally controlled by the impulse transmitted by the external excitation (time integral over two consecutive zeros


title: On heat kernel decay for the random conductance model

abstract: We study discrete time random walks in an environment of i.i.d. non-negative bounded conductances in $\mathbb{Z}^d$. We are interested in the anomaly of the heat-kernel decay. We improve recent results and techniques.


title: $\phi$ photoproduction on the proton at $E_{\gamma}$ = 1.5 - 2.9 GeV

abstract: Differential cross sections at $t=t_{\text{min}}$ and decay asymmetries for the $\gamma p\rightarrow\phi p$ reaction have been measured using linearly polarized photons in the range 1.5 to 2.9 GeV. These cross sections were used to determine the Pomeron strength factor. The cross sections and decay asymmetries are consistently described by the $t$-channel Pomeron and pseudoscalar exchange model in the $E_{\gamma}$ region above 2.37 GeV. In the lower energy region, an excess over the model prediction is observed in the energy dependence of the differential cross sections at $t=t_{\text{min}}$. This observation suggests that additional processes or interference effects between Pomeron exchange and other processes appear near the threshold region.


title: Solving the Jaynes-Cummings Model with Shift Operators Constructed by Means of the Matrix-Diagonalizing Technique

abstract: The Jaynes-Cummings model is solved with the raising and lowering (shift) operators by using the matrix-diagonalizing technique. Bell nonlocality is also found present ubiquitously in the excitations states of the model.


title: Confirmation detection in human-agent interaction using non-lexical speech cues

abstract: Even if only the acoustic channel is considered, human communication is highly multi-modal. Non-lexical cues provide a variety of information such as emotion or agreement. The ability to process such cues is highly relevant for spoken dialog systems, especially in assistance systems. In this paper we focus on the recognition of non-lexical confirmations such as "mhm", as they enhance the system's ability to accurately interpret human intent in natural communication. The architecture uses a Support Vector Machine to detect confirmations based on acoustic features. In a systematic comparison, several feature sets were evaluated for their performance on a corpus of human-agent interaction in a setting with naive users including elderly and cognitively impaired people. Our results show that using stacked formants as features yield an accuracy of 84% outperforming regular formants and MFCC or pitch based features for online classification.


title: Group-labeled light dual multinets in the projective plane (with Appendix)

abstract: In this paper we investigate light dual multinets labeled by a finite group in the projective plane $PG(2,\mathbb{K})$ defined over a field $\mathbb{K}$. We present two classes of new examples. Moreover, under some conditions on the characteristic $\mathbb{K}$, we classify group-labeled light dual multinets with lines of length least $9$.


title: Multi-Scale Pipeline for the Search of String-Induced CMB Anisotropies

abstract: We propose a multi-scale edge-detection algorithm to search for the Gott-Kaiser-Stebbins imprints of a cosmic string (CS) network on the Cosmic Microwave Background (CMB) anisotropies. Curvelet decomposition and extended Canny algorithm are used to enhance the string detectability. Various statistical tools are then applied to quantify the deviation of CMB maps having a cosmic string contribution with respect to pure Gaussian anisotropies of inflationary origin. These statistical measures include the one-point probability density function, the weighted two-point correlation function (TPCF) of the anisotropies, the unweighted TPCF of the peaks and of the up-crossing map, as well as their cross-correlation. We use this algorithm on a hundred of simulated Nambu-Goto CMB flat sky maps, covering approximately $10\%$ of the sky, and for different string tensions $G\mu$. On noiseless sky maps with an angular resolution of $0.9'$, we show that our pipeline detects CSs with $G\mu$ as low as $G\mu\gtrsim 4.3\times 10^{-10}$. At the same resolution, but with a noise level typical to a CMB-S4 phase II experiment, the detection threshold would be to $G\mu\gtrsim 1.2 \times 10^{-7}$.

categories: astro-ph.CO astro-ph.IM physics.data-an stat.AP


title: UAV-assisted Cooperative Communications with Wireless Information and Power Transfer

abstract: In this paper, we focus on a typical cooperative communication system with one pair of source and destination, where a unmanned aerial vehicle (UAV) flying from a start location to an end location serves as a mobile relay. To efficiently utilize energy in ambient environment, the UAV's transmission capability is powered exclusively by radio signal transmitted from the source via the power-splitting mechanism. In such a cooperative communication system, we study the end-to-end cooperative throughput maximization problem by optimizing the UAV's power profile, power-splitting ratio profile and trajectory for both amplify-and-forward (AF) and decode-and-forward (DF) protocols. The problem is decomposed into two subproblems: profile optimization given trajectory and trajectory optimization given profile. The former one is solved via the dual decomposition with in-depth analysis and the latter one is solved via successive convex optimization, by which a lower bound is iteratively maximized. Then the end-to-end cooperative throughput maximization problem is solved by alternately solving the two subproblems. The numerical results show that with the proposed optimal solution, choice for the UAV's power profile and power-splitting ratio profile is more long-sighted than the greedy strategy from our previous work and the successive optimization is able to converge in a few rounds of iteration. Moreover, as for the end-to-end cooperative throughput, the proposed optimal solution outperforms both static and greedy strategies, especially for the AF protocol.


title: A Many-Objective Evolutionary Algorithm with Angle-Based Selection and Shift-Based Density Estimation

abstract: Evolutionary many-objective optimization has been gaining increasing attention from the evolutionary computation research community. Much effort has been devoted to addressing this issue by improving the scalability of multiobjective evolutionary algorithms, such as Pareto-based, decomposition-based, and indicator-based approaches. Different from current work, we propose a novel algorithm in this paper called AnD, which consists of an angle-based selection strategy and a shift-based density estimation strategy. These two strategies are employed in the environmental selection to delete the poor individuals one by one. Specifically, the former is devised to find a pair of individuals with the minimum vector angle, which means that these two individuals share the most similar search direction. The latter, which takes both the diversity and convergence into account, is adopted to compare these two individuals and to delete the worse one. AnD has a simple structure, few parameters, and no complicated operators. The performance of AnD is compared with that of seven state-of-the-art many-objective evolutionary algorithms on a variety of benchmark test problems with up to 15 objectives. The experimental results suggest that AnD can achieve highly competitive performance. In addition, we also verify that AnD can be readily extended to solve constrained many-objective optimization problems.


title: Peeling off neutron skins from neutron-rich nuclei: Constraints on the symmetry energy from neutron-removal cross sections

abstract: An experimentally constrained equation of state of neutron-rich matter is fundamental for the physics of nuclei and the astrophysics of neutron stars, mergers, core-collapse supernova explosions and the synthesis of heavy elements. To this end, we investigate the potential of constraining the density dependence of the symmetry energy close to saturation density through measurements of neutron-removal cross sections in high-energy nuclear collisions of a few hundreds MeV/nucleon to GeV/nucleon. We show that the sensitivity of the total neutron-removal cross section is high enough so that the required accuracy can be reached experimentally with the recent developments of new detection techniques. We quantify two crucial points to minimize the model dependence of the approach and to reach the required accuracy: the contribution to the cross section from inelastic scattering has to be measured separately in order to allow a direct comparison of experimental cross sections to theoretical cross sections based on density functional theory and eikonal theory. The accuracy of the reaction model should be investigated and quantified by the energy and target dependence of various nucleon-removal cross sections. Our calculations explore the dependence of neutron-removal cross sections on the neutron skin of medium-heavy neutron-rich nuclei, and we demonstrate that the slope parameter L of the symmetry energy could be constrained down to +/-10 MeV by such a measurement with a 2% accuracy of the measured and calculated cross sections.


title: Full-Duplex Relay Selection in Cognitive Underlay Networks

abstract: In this work, we analyze the performance of full-duplex relay selection (FDRS) in spectrum-sharing networks. Contrary to half-duplex relaying, full-duplex relaying (FDR) enables simultaneous listening/forwarding at the secondary relay(s), thereby allowing for a higher spectral efficiency. However, since the source and relay simultaneously transmit in FDR, their superimposed signal at the primary receiver should now satisfy the existing interference constraint, which can considerably limit the secondary network throughput. In this regard, relay selection can offer an adequate solution to boost the secondary throughput while satisfying the imposed interference limit. We first analyze the performance of opportunistic FDRS with residual self-interference (RSI) by deriving the exact cumulative distribution function of its end-to-end signal-to-interference-plus-noise ratio under Nakagami-m fading. We also evaluate the offered diversity gain of relay selection for different full-duplex cooperation schemes in the presence/absence of a direct source-destination link. When the adopted RSI link gain model is sublinear in the relay power, which agrees with recent research findings, we show that remarkable diversity gain can be recovered even in the presence of an interfering direct link. Second, we evaluate the end-to-end performance of FDRS with interference constraints due to the presence of a primary receiver. Finally, the presented exact theoretical findings are verified by numerical simulations.


title: Multiple critical coupling and sensing in a microresonator-waveguide system

abstract: We study the optical transmission of a waveguide that is side-coupled to a high-$Q$ circular microresonator. The coupling is critical if the intrinsic resonator losses equate the coupling losses to the waveguide. When this happens, the transmittance of the waveguide displays resonance dips with maximal depth as the frequency is swept through the resonators resonances. We show that multiple configurations, parameterised by the minimal distance between the resonator and the waveguide, can lead to critical coupling. Indeed, for a sufficiently large resonator radius, the flow of power between the waveguide and the resonator can change sign several times within a single pass. This leads to an oscillatory coupling parameter as a function of the separation distance. As a result, multiple geometrical configurations can lead to critical coupling, even if the waveguide lies in the equatorial plane of the resonator. These results are explained using coupled-mode theory and full wave numerical simulations. In the vicinity of secondary or higher-order critical coupling, the depth of the transmittance dip is very sensitive to the environment. We discuss how this effect can be exploited for sensing purpose. Alternatively, by actively controlling the environment in the secondary critical configuration, the waveguide/resonator system can be driven as an optical switch.


title: Correlations and hyperuniformity in the avalanche size of the Oslo Model

abstract: Certain random processes display anticorrelations resulting in local Poisson-like disorder and global order, where correlations suppress fluctuations. Such processes are called hyperuniform. Using a map to an interface picture we show via analytic calculations that a sequence of avalanche sizes of the Oslo Model is hyperuniform in the temporal domain with the minimal exponent $\lambda=0$. We identify the conserved quantity in the interface picture that gives rise to the hyperuniformity in the avalanche size. We further discuss the fluctuations of the avalanche size in two variants of the Oslo Model. We support our findings with numerical results.


title: Two states

abstract: D. Bures defined a metric $\beta $ on states of a $C^*$-algebra and this concept has been generalized to unital completely positive maps $\phi : \mathcal A \to \mathcal B$, where $\mathcal B$ is either an injective $C^*$-algebra or a von Neumann algebra. We introduce a new distance $\gamma $ for the same classes of unital completely positive maps. We use in our definition the distance between representations on the same Hilbert $C^*$-module in contrast to the Bures metric which uses one representation and distinct vectors. This metric can be expressed in terms of a class of completely positive maps on free products of $C^*$-algebras and in this setting $\gamma $ looks like Wasserstein metric on probability measures. Surprisingly, when the range algebra $\mathcal B$ is injective, $\gamma $ and $\beta $ are related by the following explicit formula: $\beta ^2= 2-\sqrt{4- \gamma ^2} .$ A deep result of Choi and Li on constrained dilation is the main tool in proving this formula.


title: Molecular packing, hydrogen bonding, and fast dynamics in lysozyme/trehalose/glycerol and trehalose/glycerol glasses at low hydration

abstract: Water and glycerol are well-known to facilitate the structural relaxation of amorphous protein matrices. However, several studies evidenced that they may also limit fast ($\sim$ pico-nanosecond, ps-ns) and small-amplitude ($\sim$ \AA ) motions of proteins, which govern their stability in freeze-dried sugar mixtures. To determine how they interact with proteins and sugars in glassy matrices and, thereby, modulate their fast dynamics, we performed molecular dynamics (MD) simulations of lysozyme/trehalose/glycerol (LTG) and trehalose/glycerol (TG) mixtures at low glycerol and water concentrations. Upon addition of glycerol and/or water, the glass transition temperature, T$_{\textrm{g}}$, of LTG and TG mixtures decreases, the molecular packing of glasses is improved, and the mean-square displacements (MSDs) of lysozyme and trehalose either decrease or increase, depending on the time scale and on the temperature considered. A detailed analysis of the hydrogen bonds (HBs) formed between species reveals that water and glycerol may antiplasticize the fast dynamics of lysozyme and trehalose by increasing the total number and/or the strength of the HBs they form in glassy matrices.


title: A new multi line-cusp magnetic field plasma device (MPD) with variable magnetic field for fundamental plasma studies

abstract: One of the fundamental problems is the understanding of physics of electrostatic and electromagnetic fluctuations in multi-scale plasma turbulence. Especially so, in continuously connected plasma regions with varying degree of magnetization. Examples range from multiscale plasmas in Magnetron-like devices to astrophysical plasmas confined by magnetic dipole structures, solar wind driven collision-less and weakly collisional plasmas around Earth, to mention a few. Such plasmas are dominated by both electron scale and ion scale physics as well as finite beta effects. To investigate such processes in laboratory experiments requires excellent control of continuously connected regions of nearly zero plasma beta with finite beta regions as well as the gradient length scales of mean density and temperatures. To address some of these phenomena at laboratory scale, a new multi-line cusp configured plasma device (MPD) consisting of electromagnets with core material has been constructed with a capability to experimentally control the relative volume fractions of magnetized to unmagntized plasma volume as well as accurate control on the gradient length scales of mean density and temperature profiles. The hot tungsten cathode produced Argon plasma in the MPD has been characterised using single Langmuir probes. Argon plasma has been produced in the device over a wide range of pressure 5 x 10-5 mBar to 1 x 10-3 mBar, achieving plasma density range from 109 to 1011cm-3 and temperature in the range 1eV to 8eV.


title: Forward Private Searchable Symmetric Encryption with Optimized I/O Efficiency

abstract: Recently, several practical attacks raised serious concerns over the security of searchable encryption. The attacks have brought emphasis on forward privacy, which is the key concept behind solutions to the adaptive leakage-exploiting attacks, and will very likely to become mandatory in the design of new searchable encryption schemes. For a long time, forward privacy implies inefficiency and thus most existing searchable encryption schemes do not support it. Very recently, Bost (CCS 2016) showed that forward privacy can be obtained without inducing a large communication overhead. However, Bost's scheme is constructed with a relatively inefficient public key cryptographic primitive, and has a poor I/O performance. Both of the deficiencies significantly hinder the practical efficiency of the scheme, and prevent it from scaling to large data settings. To address the problems, we first present FAST, which achieves forward privacy and the same communication efficiency as Bost's scheme, but uses only symmetric cryptographic primitives. We then present FASTIO, which retains all good properties of FAST, and further improves I/O efficiency. We implemented the two schemes and compared their performance with Bost's scheme. The experiment results show that both our schemes are highly efficient, and FASTIO achieves a much better scalability due to its optimized I/O.


title: Higgs production at future $e^+e^-$ colliders in the Georgi-Machacek model

abstract: We study how the dominant single and double SM-like Higgs ($h$) production at future $e^+e^-$ colliders is modified in the Georgi-Machacek (GM) model. On imposing theoretical, indirect and direct constraints, significant deviations of $h$-couplings from their SM values are still possible; for instance, the Higgs-gauge coupling coupling can be corrected by a factor $\kappa_{hVV}\in[0.93,1.15]$ in the allowed parameter space. For the Higgs-strahlung $e^+e^-\to hZ$ and vector boson fusion processes $e^+e^-\to h\nu\bar{\nu},~he^+e^-$, the cross section could increase by $32\%$ or decrease by $13\%$. In the case of associated production with a top quark pair $e^+e^-\to ht\bar{t}$, the cross section can be enhanced up to several times when the custodial triplet scalar $H_3^0$ is resonantly produced. In the meanwhile, the double Higgs production $e^+e^-\to hhZ~(hh\nu\bar{\nu})$ can be maximally enhanced by one order of magnitude at the resonant $H_{1,3}^0$ production. We also include exclusion limits expected from future LHC runs at higher energy and luminosity and discuss their further constraints on the relevant model parameters. We find that the GM model can result in likely measurable deviations of Higgs production from the SM at future $e^+e^-$ colliders.


title: Project Solaris, a Global Network of Autonomous Observatories: Design, Commissioning, and First Science Results

abstract: We present the design and commissioning of Project Solaris, a global network of autonomous observatories. Solaris is a Polish scientific undertaking aimed at the detection and characterization of circumbinary exoplanets and eclipsing binary stars. To accomplish this, a network of four fully autonomous observatories has been deployed in the Southern Hemisphere: Solaris-1 and Solaris-2 in the South African Astronomical Observatory in South Africa; Solaris-3 in Siding Spring Observatory in Australia; and Solaris-4 in Complejo Astronomico El Leoncito in Argentina. The four stations are nearly identical and are equipped with 0.5-m Ritchey-Cr\'etien (f/15) or Cassegrain (f/9, Solaris-3) optics and high-grade 2 K x 2 K CCD cameras with Johnson and Sloan filter sets. We present the design and implementation of low-level security; data logging and notification systems; weather monitoring components; all-sky vision system, surveillance system; and distributed temperature and humidity sensors. We describe dedicated grounding and lighting protection system design and robust fiber data transfer interfaces in electrically demanding conditions. We discuss the outcomes of our design, as well as the resulting software engineering requirements. We describe our system's engineering approach to achieve the required level of autonomy, the architecture of the custom high-level industry-grade software that has been designed and implemented specifically for the use of the network. We present the actual status of the project and first photometric results; these include data and models of already studied systems for benchmarking purposes (Wasp-4b, Wasp-64b, and Wasp-98b transits, PG 1663-018, an eclipsing binary with a pulsator) as well J024946-3825.6, an interesting low-mass binary system for which a complete model is provided for the first time.


title: Insulating-to-Conducting Behavior and Band Profile Across the La0.9Ba0.1MnO3/Nb:SrTiO3 Epitaxial Interface

abstract: La0.9Ba0.1MnO3 is a ferromagnetic insulator in its bulk form, but exhibits metallicity in thin film form. It has a wide potential in a range of spintronic-related applications, and hence it is critical to understand thickness-dependent electronic structure in thin films as well as substrate/film interface effects. Here, using electrical and in-situ photoemission spectroscopy measurements, we report the electronic structure and interface band profile of high-quality layer-by-layer-grown La0.9Ba0.1MnO3 on single crystal Nb:SrTiO3 substrates. A transition from insulating-to-conducting was observed with increasing La0.9Ba0.1MnO3 thickness, which was explained by the determined interface band diagram of La0.9Ba0.1MnO3/Nb: SrTiO3, where a type II heterojunction was formed.


title: Unsupervised Segmentation of Action Segments in Egocentric Videos using Gaze

abstract: Unsupervised segmentation of action segments in egocentric videos is a desirable feature in tasks such as activity recognition and content-based video retrieval. Reducing the search space into a finite set of action segments facilitates a faster and less noisy matching. However, there exist a substantial gap in machine understanding of natural temporal cuts during a continuous human activity. This work reports on a novel gaze-based approach for segmenting action segments in videos captured using an egocentric camera. Gaze is used to locate the region-of-interest inside a frame. By tracking two simple motion-based parameters inside successive regions-of-interest, we discover a finite set of temporal cuts. We present several results using combinations (of the two parameters) on a dataset, i.e., BRISGAZE-ACTIONS. The dataset contains egocentric videos depicting several daily-living activities. The quality of the temporal cuts is further improved by implementing two entropy measures.


title: Reducing backaction when measuring temporal correlations in quantum systems

abstract: Dynamic correlations of quantum observables are challenging to measure due to measurement backaction incurred at early times. Recent work [P. Uhrich et al., Phys. Rev. A, 96:022127 (2017)] has shown that ancilla-based noninvasive measurements are able to reduce this backaction, allowing for dynamic correlations of single-site spin observables to be measured. We generalise this result to correlations of arbitrary spin observables and extend the measurement protocol to simultaneous noninvasive measurements which allow for real and imaginary parts of correlations to be extracted from a single set of measurements. We use positive operator-valued measures to analyse the dynamics generated by the ancilla-based measurements. Using this framework we prove that special observables exist for which measurement backaction is of no concern, so that dynamic correlations of these can be obtained without making use of ancillas.


title: Unsupervised Classification of Intrusive Igneous Rock Thin Section Images using Edge Detection and Colour Analysis

abstract: Classification of rocks is one of the fundamental tasks in a geological study. The process requires a human expert to examine sampled thin section images under a microscope. In this study, we propose a method that uses microscope automation, digital image acquisition, edge detection and colour analysis (histogram). We collected 60 digital images from 20 standard thin sections using a digital camera mounted on a conventional microscope. Each image is partitioned into a finite number of cells that form a grid structure. Edge and colour profile of pixels inside each cell determine its classification. The individual cells then determine the thin section image classification via a majority voting scheme. Our method yielded successful results as high as 90% to 100% precision.


title: Power analysis for a linear regression model when regressors are matrix sampled

abstract: Multiple matrix sampling is a survey methodology technique that randomly chooses a relatively small subset of items to be presented to survey respondents for the purpose of reducing respondent burden. The data produced are missing completely at random (MCAR), and special missing data techniques should be used in linear regression and other multivariate statistical analysis. We derive asymptotic variances of regression parameter estimates that allow us to conduct power analysis for linear regression models fit to the data obtained via a multiple matrix sampling design. The ideas are demonstrated with a variation of the Big Five Inventory of psychological traits. An exploration of the regression parameter space demonstrates instability of the sample size requirements, and substantial losses of precision with matrix-sampled regressors. A simulation with non-normal data demonstrates the advantages of a semi-parametric multiple imputation scheme.


title: Torsion and K-theory for some free wreath products

abstract: We classify torsion actions of free wreath products of arbitrary compact quantum groups and use this to prove that if $\mathbb{G}$ is a torsion-free compact quantum group satisfying the strong Baum-Connes property, then $\mathbb{G}\wr_{\ast}S_{N}^{+}$ also satisfies the strong Baum-Connes property. We then compute the K-theory of free wreath products of classical and quantum free groups by $SO_{q}(3)$.


title: Interplay of Kondo Effect and RKKY Interaction

abstract: A tutorial introduction to the Kondo effect, the RKKY interaction, and a renormalization group treatment to analyze the competition between the two.


title: A Catalog of Cool Dwarf Targets for the Transiting Exoplanet Survey Satellite

abstract: We present a catalog of cool dwarf targets ($V-J>2.7$, $T_{\rm eff} \lesssim 4000 K$) and their stellar properties for the upcoming Transiting Exoplanet Survey Satellite (TESS), for the purpose of determining which cool dwarfs should be observed using two-minute observations. TESS has the opportunity to search tens of thousands of nearby, cool, late K and M-type dwarfs for transiting exoplanets, an order of magnitude more than current or previous transiting exoplanet surveys, such as {\it Kepler}, K2 and ground-based programs. This necessitates a new approach to choosing cool dwarf targets. Cool dwarfs were chosen by collating parallax and proper motion catalogs from the literature and subjecting them to a variety of selection criteria. We calculate stellar parameters and TESS magnitudes using the best possible relations from the literature while maintaining uniformity of methods for the sake of reproducibility. We estimate the expected planet yield from TESS observations using statistical results from the Kepler Mission, and use these results to choose the best targets for two-minute observations, optimizing for small planets for which masses can conceivably be measured using follow up Doppler spectroscopy by current and future Doppler spectrometers. The catalog is incorporated into the TESS Input Catalog and TESS Candidate Target List until a more complete and accurate cool dwarf catalog identified by ESA's Gaia Mission can be incorporated.


title: Where computer vision can aid physics: dynamic cloud motion forecasting from satellite images

abstract: This paper describes a new algorithm for solar energy forecasting from a sequence of Cloud Optical Depth (COD) images. The algorithm is based on the following simple observation: the dynamics of clouds represented by COD images resembles the motion (transport) of a density in a fluid flow. This suggests that, to forecast the motion of COD images, it is sufficient to forecast the flow. The latter, in turn, can be accomplished by fitting a parametric model of the fluid flow to the COD images observed in the past. Namely, the learning phase of the algorithm is composed of the following steps: (i) given a sequence of COD images, the snapshots of the optical flow are estimated from two consecutive COD images; (ii) these snapshots are then assimilated into a Navier-Stokes Equation (NSE), i.e. an initial velocity field for NSE is selected so that the corresponding NSE' solution is as close as possible to the optical flow snapshots. The prediction phase consists of utilizing a linear transport equation, which describes the propagation of COD images in the fluid flow predicted by NSE, to estimate the future motion of the COD images. The algorithm has been tested on COD images provided by two geostationary operational environmental satellites from NOAA serving the west-hemisphere.

categories: math.OC cs.CV physics.flu-dyn physics.geo-ph


title: Nonequilibrium Green's function theory for nonadiabatic effects in quantum electron transport

abstract: We develop nonequilibribrium Green's function based transport theory, which includes effects of nonadiabatic nuclear motion in the calculation of the electric current in molecular junctions. Our approach is based on the separation of slow and fast timescales in the equations of motion for the Green's functions by means of the Wigner representation. Time derivatives with respect to central time serves as a small parameter in the perturbative expansion enabling the computation of nonadiabatic corrections to molecular Green's functions. Consequently, we produce series of analytic expressions for non-adiabatic electronic Green's functions (up to the second order in the central time derivatives); which depend not solely on instantaneous molecular geometry but likewise on nuclear velocities and accelerations. Extended formula for electric current is derived which accounts for the non-adiabatic corrections. This theory is concisely illustrated by the calculations on a model molecular junction.


title: New binary and ternary LCD codes

abstract: LCD codes are linear codes with important cryptographic applications. Recently, a method has been presented to transform any linear code into an LCD code with the same parameters when it is supported on a finite field with cardinality larger than 3. Hence, the study of LCD codes is mainly open for binary and ternary fields. Subfield-subcodes of $J$-affine variety codes are a generalization of BCH codes which have been successfully used for constructing good quantum codes. We describe binary and ternary LCD codes constructed as subfield-subcodes of $J$-affine variety codes and provide some new and good LCD codes coming from this construction.


title: Matching Anonymized and Obfuscated Time Series to Users' Profiles

abstract: Many popular applications use traces of user data to offer various services to their users. However, even if user data is anonymized and obfuscated, a user's privacy can be compromised through the use of statistical matching techniques that match a user trace to prior user behavior. In this work, we derive the theoretical bounds on the privacy of users in such a scenario. We build on our recent study in the area of location privacy, in which we introduced formal notions of location privacy for anonymization-based location privacy-protection mechanisms. Here we derive the fundamental limits of user privacy when both anonymization and obfuscation-based protection mechanisms are applied to users' time series of data. We investigate the impact of such mechanisms on the trade-off between privacy protection and user utility. We first study achievability results for the case where the time-series of users are governed by an i.i.d. process. The converse results are proved both for the i.i.d. case as well as the more general Markov chain model. We demonstrate that as the number of users in the network grows, the obfuscation-anonymization plane can be divided into two regions: in the first region, all users have perfect privacy; and, in the second region, no user has privacy.


title: Asymptotics for the Heat Kernel on H-Type Groups

abstract: We give sharp asymptotic estimates at infinity of all radial partial derivatives of the heat kernel on H-type groups. As an application, we give a new proof of the discreteness of the spectrum of some natural sub-Riemannian Ornstein-Uhlenbeck operators on these groups.


title: Quasi-periodic solution of quasi-linear fifth-order KdV equation

abstract: In this paper, we prove the existence of quasi-periodic small-amplitude solutions for quasi-linear Hamiltonian perturbation of the fifth-order KdV equation on the torus in presence of a quasi-periodic forcing.


title: Exotic superfluidity and pairing phenomena in atomic Fermi gases in mixed dimensions

abstract: Atomic Fermi gases have been an ideal platform for simulating conventional and engineering exotic physical systems owing to their multiple tunable control parameters. Here we investigate the effects of mixed dimensionality on the superfluid and pairing phenomena of a two-component ultracold atomic Fermi gas with a short-range pairing interaction, while one component is confined on a one-dimensional (1D) optical lattice whereas the other is in a homogeneous 3D continuum. We study the phase diagram and the pseudogap phenomena throughout the entire BCS-BEC crossover, using a pairing fluctuation theory. We find that the effective dimensionality of the non-interacting lattice component can evolve from quasi-3D to quasi-1D, leading to strong Fermi surface mismatch. Upon pairing, the system becomes effectively quasi-two dimensional in the BEC regime. The behavior of $T_c$ bears similarity to that of a regular 3D population imbalanced Fermi gas, but with a more drastic departure from the regular 3D balanced case, featuring both intermediate temperature superfluidity and possible pair density wave ground state. Unlike a simple 1D optical lattice case, $T_c$ in the mixed dimensions has a constant BEC asymptote.


title: Full Szeg\H{o}-type trace asymptotics for ergodic operators on large boxes

abstract: We prove full Szeg\H{o}-type large-box trace asymptotics for selfadjoint $\mathbb{Z}^d$-ergodic operators $\Omega\ni \omega\mapsto H_\omega$ acting on $L^2(\mathbb{R}^d)$. More precisely, let $g$ be a bounded, compactly supported and real-valued function such that the (averaged) operator kernel of $g(H_\omega)$ decays sufficiently fast, and let $h$ be a sufficiently smooth compactly supported function. We then prove a full asymptotic expansion of the averaged trace of the operator $h(g(H_\omega)_{[-L,L]^d})$ in terms of the length-scale $L$.


title: Diagramming the Class Diagram: Toward a Unified Modeling Methodology

abstract: The object-oriented class is, in general, the most utilized element in programming and modeling. It is employed throughout the software development process, from early domain analysis phases to later maintenance phases. A class diagram typically uses elements of graph theory, e.g., boxes, ovals, lines. Many researchers have examined the class diagram layout from different perspectives, including visibility, juxtaposability, and aesthetics. While software systems can be incredibly complex, class diagrams represent a very broad picture of the system as a whole. The key to understanding of such complexity is use of tools such as diagrams at various levels of representation. This paper develops a more elaborate diagrammatic description of the class diagram that includes flows of attributes, thus providing a basic representation for specifying behavior and control instead of merely listing methods.


title: Free energy of a folded polymer under cylindrical confinement

abstract: Monte Carlo computer simulations are used to study the conformational free energy of a folded polymer confined to a long cylindrical tube. The polymer is modeled as a hard-sphere chain. Its conformational free energy $F$ is measured as a function of $\lambda$, the end-to-end distance of the polymer. In the case of a flexible linear polymer, $F(\lambda)$ is a linear function in the folded regime with a gradient that scales as $f\equiv |dF/d\lambda| \sim N^0 D^{-1.20\pm 0.01}$ for a tube of diameter $D$ and a polymer of length $N$. This is close to the prediction $f \sim N^0 D^{-1}$ obtained from simple scaling arguments. The discrepancy is due in part to finite-size effects associated with the de-Gennes blob model. A similar discrepancy was observed for the folding of a single arm of a three-arm star polymer. We also examine backfolding of a semiflexible polymer of persistence length $P$ in the classic Odijk regime. In the overlap regime, the derivative scales $f \sim N^0 D^{-1.72\pm 0.02} P^{-0.35\pm 0.01}$, which is close to the prediction $f \sim N^0 D^{-5/3} P^{-1/3}$ obtained from a scaling argument that treats interactions between deflection segments at the second virial level. In addition, the measured free energy cost of forming a hairpin turn is quantitatively consistent with a recent theoretical calculation. Finally, we examine the scaling of $F(\lambda)$ for a confined semiflexible chain in the presence of an S-loop composed of two hairpins. While the predicted scaling of the free energy gradient is the same as that for a single hairpin, we observe a scaling of $f \sim D^{-1.91\pm 0.03} P^{-0.36\pm 0.01}$. Thus, the quantitative discrepancy between this measurement and the predicted scaling is somewhat greater for S-loops than for single hairpins.


title: Enabling Quality Control for Entity Resolution: A Human and Machine Cooperation Framework

abstract: Even though many machine algorithms have been proposed for entity resolution, it remains very challenging to find a solution with quality guarantees. In this paper, we propose a novel HUman and Machine cOoperation (HUMO) framework for entity resolution (ER), which divides an ER workload between the machine and the human. HUMO enables a mechanism for quality control that can flexibly enforce both precision and recall levels. We introduce the optimization problem of HUMO, minimizing human cost given a quality requirement, and then present three optimization approaches: a conservative baseline one purely based on the monotonicity assumption of precision, a more aggressive one based on sampling and a hybrid one that can take advantage of the strengths of both previous approaches. Finally, we demonstrate by extensive experiments on real and synthetic datasets that HUMO can achieve high-quality results with reasonable return on investment (ROI) in terms of human cost, and it performs considerably better than the state-of-the-art alternatives in quality control.


title: Bag-of-Vector Embeddings of Dependency Graphs for Semantic Induction

abstract: Vector-space models, from word embeddings to neural network parsers, have many advantages for NLP. But how to generalise from fixed-length word vectors to a vector space for arbitrary linguistic structures is still unclear. In this paper we propose bag-of-vector embeddings of arbitrary linguistic graphs. A bag-of-vector space is the minimal nonparametric extension of a vector space, allowing the representation to grow with the size of the graph, but not tying the representation to any specific tree or graph structure. We propose efficient training and inference algorithms based on tensor factorisation for embedding arbitrary graphs in a bag-of-vector space. We demonstrate the usefulness of this representation by training bag-of-vector embeddings of dependency graphs and evaluating them on unsupervised semantic induction for the Semantic Textual Similarity and Natural Language Inference tasks.


title: On new theta identities of fermion correlation functions on genus g Riemann surfaces

abstract: Theta identities on genus g Riemann surfaces which decompose simple products of fermion correlation functions with a constraint on their variables are considered. This type of theta identities is, in a sense, dual to Fay s formula, by which it is possible to sum over spin structures of certain part of superstring amplitudes in NSR formalism without using Fay s formula nor Riemann s theta formula in much simpler, more transparent way. Also, such identities will help to cast correlation functions among arbitrary numbers of Kac-Moody currents in a closed form. As for genus 1, the identities are reported before in ref[1] [2]. Based on some notes on genus 1 case which were not reported in ref[1] [2] and relating those to the results of the Dolan Goddard method ref[3] on describing Kac-Moody currents in a closed form, we propose an idea of generalizing genus 1 identities to the case of genus g surfaces. This is not a complete derivation of the higher genus formula due to difficulties of investigating singular part of derivatives of genus g Weierstrass Pe functions. Mathematical issues remained unsolved for genus g >1 are described in the text.


title: Sensitivity to Thermonuclear Reaction Rates in Modeling the Abundance Anomalies of NGC 2419

abstract: Abundance anomalies in globular clusters provide strong evidence for multiple stellar populations within each cluster. These populations are usually interpreted as distinct generations, with the currently observed second-generation stars having formed in part from the ejecta of massive, first-generation "polluter" stars, giving rise to the anomalous abundance patterns. The precise nature of the polluters and their enrichment mechanism are still unclear. Even so, the chemical abundances measured in second-generation stars within the globular cluster NGC 2419 provide insight into this puzzling process. Previous work used Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that could reproduce the observed abundances, thereby placing robust limits on the origins of the polluter material. The effect of individual reaction rates on these conditions has not been studied, however. Thus, we perform an exhaustive sensitivity study on the nuclear physics input to determine which reactions have the greatest impact on these predictions. We find that the $^{30}$Si(p,$\gamma$)$^{31}$P, $^{37}$Ar(p,$\gamma$)$^{38}$K, $^{38}$Ar(p,$\gamma$)$^{39}$K, and $^{39}$K(p,$\gamma$)$^{40}$Ca reactions are all critical in determining the temperature-density conditions, and ultimately, the origin of the polluter material. We conclude with recommendations for future experiments.


title: Efficient Dynamic Searchable Encryption with Forward Privacy

abstract: Searchable symmetric encryption (SSE) enables a client to perform searches over its outsourced encrypted files while preserving privacy of the files and queries. Dynamic schemes, where files can be added or removed, leak more information than static schemes. For dynamic schemes, forward privacy requires that a newly added file cannot be linked to previous searches. We present a new dynamic SSE scheme that achieves forward privacy by replacing the keys revealed to the server on each search. Our scheme is efficient and parallelizable and outperforms the best previous schemes providing forward privacy, and achieves competitive performance with dynamic schemes without forward privacy. We provide a full security proof in the random oracle model. In our experiments on the Wikipedia archive of about four million pages, the server takes one second to perform a search with 100,000 results.


title: Improved Training for Self-Training by Confidence Assessments

abstract: It is well known that for some tasks, labeled data sets may be hard to gather. Therefore, we wished to tackle here the problem of having insufficient training data. We examined learning methods from unlabeled data after an initial training on a limited labeled data set. The suggested approach can be used as an online learning method on the unlabeled test set. In the general classification task, whenever we predict a label with high enough confidence, we treat it as a true label and train the data accordingly. For the semantic segmentation task, a classic example for an expensive data labeling process, we do so pixel-wise. Our suggested approaches were applied on the MNIST data-set as a proof of concept for a vision classification task and on the ADE20K data-set in order to tackle the semi-supervised semantic segmentation problem.


title: Testing for Feature Relevance: The HARVEST Algorithm

abstract: Feature selection with high-dimensional data and a very small proportion of relevant features poses a severe challenge to standard statistical methods. We have developed a new approach (HARVEST) that is straightforward to apply, albeit somewhat computer-intensive. This algorithm can be used to pre-screen a large number of features to identify those that are potentially useful. The basic idea is to evaluate each feature in the context of many random subsets of other features. HARVEST is predicated on the assumption that an irrelevant feature can add no real predictive value, regardless of which other features are included in the subset. Motivated by this idea, we have derived a simple statistical test for feature relevance. Empirical analyses and simulations produced so far indicate that the HARVEST algorithm is highly effective in predictive analytics, both in science and business.


title: The Deep Ritz method: A deep learning-based numerical algorithm for solving variational problems

abstract: We propose a deep learning based method, the Deep Ritz Method, for numerically solving variational problems, particularly the ones that arise from partial differential equations. The Deep Ritz method is naturally nonlinear, naturally adaptive and has the potential to work in rather high dimensions. The framework is quite simple and fits well with the stochastic gradient descent method used in deep learning. We illustrate the method on several problems including some eigenvalue problems.


title: Spherical embeddings of symmetric association schemes in 3-dimensional Euclidean space

abstract: We classify the symmetric association schemes with faithful spherical embedding in 3-dimensional Euclidean space. Our result is based on previous research on primitive association schemes with $m_1 = 3$.


title: Two-way Two-tape Automata

abstract: In this article we consider two-way two-tape (alternating) automata accepting pairs of words and we study some closure properties of this model. Our main result is that such alternating automata are not closed under complementation for non-unary alphabets. This improves a similar result of Kari and Moore for picture languages. We also show that these deterministic, non-deterministic and alternating automata are not closed under composition.


title: An elementary proof of the group law for elliptic curves

abstract: We give an elementary proof of the group law for elliptic curves using explicit formulas.


title: An silicon nanocrystal laser

abstract: Silicon lasers have been the most challenging element for the monolithic integrated Si photonics. Here we report the first successful all-Si laser at room temperature based on silicon nanocrystals (Si NCs) with high optical gains. The active Si NC layer was made from hydrogen silsesquioxane (HSQ) that had undergone a phase separation annealing, followed by high-pressure hydrogen passivation. Using the nanoimprint technique, a second-order Bragg grating was made on the active layer as the distributed feedback (DFB) laser cavity. This Si NC DFB laser device was optically pumped with femtosecond pulses. Criteria for lasing action, i.e., the threshold effect, the polarization dependence and the significant spectral narrowing of stimulated emission were evidently fulfilled, which indicated the realization of an all-Si laser.


title: Maxwell Strata and Cut Locus in Sub-Riemannian Problem on Engel group

abstract: We consider the nilpotent left-invariant sub-Riemannian structure on the Engel group. This structure gives a fundamental local approximation of a generic rank 2 sub-Riemannian structure on a 4-manifold near a generic point (in particular, of the kinematic models of a car with a trailer). On the other hand, this is the simplest sub-Riemannian structure of step three. We describe the global structure of the cut locus (the set of points where geodesics lose their global optimality), the Maxwell set (the set of points that admit more than one minimizer), and the intersection of the cut locus with the caustic (the set of conjugate points along all geodesics). The group of symmetries of the cut locus is described: it is generated by a one-parameter group of dilations $\mathbb{R}_+$ and a discrete group of reflections $\mathbb{Z}_2 \times \mathbb{Z}_2 \times \mathbb{Z}_2$. The cut locus admits a stratification with 6 three-dimensional strata, 12 two-dimensional strata, and 2 one-dimensional strata. Three-dimen-sional strata of the cut locus are Maxwell strata of multiplicity 2 (for each point there are 2 minimizers). Two-dimensional strata of the cut locus consist of conjugate points. Finally, one-dimensional strata are Maxwell strata of infinite multiplicity, they consist of conjugate points as well. Projections of sub-Riemannian geodesics to the 2-dimensional plane of the distribution are Euler elasticae. For each point of the cut locus, we describe the Euler elasticae corresponding to minimizers coming to this point. Finally, we describe the structure of the optimal synthesis, i.e., the set of minimizers for each terminal point in the Engel group.