title: Global Optimality in Separable Dictionary Learning with Applications to the Analysis of Diffusion MRI

abstract: Sparse dictionary learning is a popular method for representing signals as linear combinations of a few elements from a dictionary that is learned from the data. In the classical setting, signals are represented as vectors and the dictionary learning problem is posed as a matrix factorization problem where the data matrix is approximately factorized into a dictionary matrix and a sparse matrix of coefficients. However, in many applications in computer vision and medical imaging, signals are better represented as matrices or tensors (e.g. images or videos), where it may be beneficial to exploit the multi-dimensional structure of the data to learn a more compact representation. One such approach is separable dictionary learning, where one learns separate dictionaries for different dimensions of the data. However, typical formulations involve solving a non-convex optimization problem; thus guaranteeing global optimality remains a challenge. In this work, we propose a framework that builds upon recent developments in matrix factorization to provide theoretical and numerical guarantees of global optimality for separable dictionary learning. We propose an algorithm to find such a globally optimal solution, which alternates between following local descent steps and checking a certificate for global optimality. We illustrate our approach on diffusion magnetic resonance imaging (dMRI) data, a medical imaging modality that measures water diffusion along multiple angular directions in every voxel of an MRI volume. State-of-the-art methods in dMRI either learn dictionaries only for the angular domain of the signals or in some cases learn spatial and angular dictionaries independently. In this work, we apply the proposed separable dictionary learning framework to learn spatial and angular dMRI dictionaries jointly and provide preliminary validation on denoising phantom and real dMRI brain data.


title: Asymptotic behavior of least energy solutions to the Lane-Emden system near the critical hyperbola

abstract: The Lane-Emden system is written as \begin{equation*} \begin{cases} -\Delta u = v^p &\text{in } \Omega,\\ -\Delta v = u^q &\text{in } \Omega,\\ u, v > 0 &\text{in } \Omega,\\ u = v = 0 &\text{on } \partial \Omega \end{cases} \end{equation*} where $\Omega$ is a smooth bounded domain in the Euclidean space $\mathbb{R}^n$ for $n \ge 3$ and $0< p< q <\infty$. The asymptotic behavior of least energy solutions near the critical hyperbola was studied by Guerra \cite{G} when $p \geq 1$ and the domain is convex. In this paper, we cover all the remaining cases $p < 1$ and extend the results to any smooth bounded domain.


title: Deep Learning for Semantic Segmentation on Minimal Hardware

abstract: Deep learning has revolutionised many fields, but it is still challenging to transfer its success to small mobile robots with minimal hardware. Specifically, some work has been done to this effect in the RoboCup humanoid football domain, but results that are performant and efficient and still generally applicable outside of this domain are lacking. We propose an approach conceptually different from those taken previously. It is based on semantic segmentation and does achieve these desired properties. In detail, it is being able to process full VGA images in real-time on a low-power mobile processor. It can further handle multiple image dimensions without retraining, it does not require specific domain knowledge for achieving a high frame rate and it is applicable on a minimal mobile hardware.


title: Bound states and QCD

abstract: The similarities of hadrons and atoms motivate a study of the principles of QED bound states and of their applicability to QCD. The power series in $\alpha$ and $\log\alpha$ of the binding energy is reflected in the Fock expansion of the bound state in temporal gauge ($A^0=0$). Gauss' constraint on physical states fixes the gauge for time independent transformations and determines the instantaneous interaction within each Fock state. Positronium atoms generate a classical (dipole) electric field, whereas there can be no color octet gluon field for color singlet hadrons. Hence the gluon field generated by each color component of a hadron need not vanish at spatial infinity. Gauss' constraint has a homogeneous solution with a single parameter $\Lambda$ that is compatible with Poincar\'e invariance. The corresponding potential is linear for $q\bar q$ and $gg$ Fock states, and confining also for other states ($q\bar qg,\,qqq$). This approach is consistent with the quarkonium phenomenology based on the Cornell potential at lowest order. The relativistic meson and glueball eigenstates of the QCD Hamiltonian with the $O(\alpha_s^0)$ linear potential are determined. The states lie on linear Regge trajectories and their daughters. There are also massless bound states which allow to include a $J^{PC}=0^{++}$ condensate in the perturbative vacuum, thus breaking chiral symmetry spontaneously.


title: Inequalities for $L^p$-norms that sharpen the triangle inequality and complement Hanner's Inequality

abstract: In 2006 Carbery raised a question about an improvement on the na\"ive norm inequality $\|f+g\|_p^p \leq 2^{p-1}(\|f\|_p^p + \|g\|_p^p)$ for two functions in $L^p$ of any measure space. When $f=g$ this is an equality, but when the supports of $f$ and $g$ are disjoint the factor $2^{p-1}$ is not needed. Carbery's question concerns a proposed interpolation between the two situations for $p>2$. The interpolation parameter measuring the overlap is $\|fg\|_{p/2}$. We prove an inequality of this type that is stronger than the one Carbery proposed. Moreover, our stronger inequalities are valid for all $p$.


title: Modeling Daily Seasonality of Mexico City Ozone using Nonseparable Covariance Models on Circles Cross Time

abstract: Mexico City tracks ground-level ozone levels to assess compliance with national ambient air quality standards and to prevent environmental health emergencies. Ozone levels show distinct daily patterns, within the city, and over the course of the year. To model these data, we use covariance models over space, circular time, and linear time. We review existing models and develop new classes of nonseparable covariance models of this type, models appropriate for quasi-periodic data collected at many locations. With these covariance models, we use nearest-neighbor Gaussian processes to predict hourly ozone levels at unobserved locations in April and May, the peak ozone season, to infer compliance to Mexican air quality standards and to estimate respiratory health risk associated with ozone. Predicted compliance with air quality standards and estimated respiratory health risk vary greatly over space and time. In some regions, we predict exceedance of national standards for more than a third of the hours in April and May. On many days, we predict that nearly all of Mexico City exceeds nationally legislated ozone thresholds at least once. In peak regions, we estimate respiratory risk for ozone to be 55% higher on average than the annual average risk and as much at 170% higher on some days.


title: Some remarks regarding l elements defined in algebras obtained by the Cayley-Dickson process

abstract: In this paper, we define a special class of elements in the algebras obtained by the Cayley Dickson process, called l elements. We find conditions such that these elements to be invertible. These conditions can be very useful for finding new identities, identities which can help us in the study of the properties of these algebras.


title: Latency-Energy Tradeoff based on Channel Scheduling and Repetitions in NB-IoT Systems

abstract: Narrowband IoT (NB-IoT) is the latest IoT connectivity solution presented by the 3GPP. NB-IoT introduces coverage classes and introduces a significant link budget improvement by allowing repeated transmissions by nodes that experience high path loss. However, those repetitions necessarily increase the energy consumption and the latency in the whole NB-IoT system. The extent to which the whole system is affected depends on the scheduling of the uplink and downlink channels. We address this question, not treated previously, by developing a tractable model of NB-IoT access protocol operation, comprising message exchanges in random-access, control, and data channels, both in the uplink and downlink. The model is then used to analyze the impact of channel scheduling as well as the interaction of coexisting coverage classes, through derivation of the expected latency and battery lifetime for each coverage class. These results are subsequently employed in investigation of latency-energy tradeoff in NB-IoT channel scheduling as well as determining the optimized operation points. Simulations results show validity of the analysis and confirm that there is a significant impact of channel scheduling on latency and lifetime performance of NBIoT devices.


title: Three Blocks Solvable Lattice Models and Birman--Murakami--Wenzl Algebra

abstract: Birman--Murakami--Wenzl (BMW) algebra was introduced in connection with knot theory. We treat here interaction round the face solvable (IRF) lattice models. We assume that the face transfer matrix obeys a cubic polynomial equation, which is called the three block case. We prove that the three block theories all obey the BMW algebra. We exemplify this result by treating in detail the $SU(2)$ $2\times 2$ fused models, and showing explicitly the BMW structure. We use the connection between the construction of solvable lattice models and conformal field theory. This result is important to the solution of IRF lattice models and the development of new models, as well as to knot theory.


title: Chiral Magnetic Effect in Isobaric Collisions from Anomalous-Viscous Fluid Dynamics (AVFD)

abstract: The isobaric collision experiment at RHIC provides the unique opportunity to detect the possible signal of Chiral Magnetic Effect (CME) in heavy ion collisions. The idea is to contrast the correlation observables of the two colliding systems that supposedly have identical flow-driven background contributions while quite different CME signal contributions due to the 10% variation in their nuclear charge and thus magnetic field strength. With the recently developed quantitative simulation tool for computing CME signal, the Anomalous-Viscous Fluid Dynamics (AVFD), we demonstrate that a joint (multiplicity + elliptic-flow) event selection is crucial for this purpose. We further propose to use the absolute difference between RuRu and ZrZr events (after using identical event selection) for detecting CME signal and make predictions for the correlation observables.


title: Characterization of Model-Based Uncertainties in Incompressible Turbulent Flows by Machine Learning

abstract: This work determines the inaccuracy of using Reynolds averaged Navier Stokes (RANS) turbulence models in transition to turbulent flow regimes by predicting the model-based discrepancies between RANS and large eddy simulation (LES) models and then incorporates the capabilities of machine learning algorithms to characterize the discrepancies which are defined as a function of mean flow properties of RANS simulations. First, three-dimensional CFD simulations using k-omega Shear Stress Transport (SST) and dynamic one-equation subgrid-scale models are conducted in a wall-bounded channel containing a cylinder for RANS and LES, respectively, to identify the turbulent kinetic energy discrepancy. Second, several flow features such as viscosity ratio, wall-distance based Reynolds number, and vortex stretching are calculated from the mean flow properties of RANS. Then these flow features are regressed onto the discrepancy using a Random Forests regression algorithm. Finally, the discrepancy of the test flow is predicted using the trained algorithm. The results reveal that a significant discrepancy exists between RANS and LES simulations and ML algorithm successfully predicts the increased model uncertainties caused by the employment of k-omega SST turbulence model for transitional fluid flows.


title: On the Significance of the Quantum Mechanical Covariance Matrix

abstract: The characterization of quantum correlations, being stronger than classical, yet weaker than those appearing in non-signaling models, still poses many riddles. In this work we show that the extent of binary correlations in a general class of nonlocal theories can be characterized by the~existence of a certain covariance matrix. The set of quantum realizable two-point correlators in the~bipartite case then arises from a subtle restriction on the structure of this general covariance matrix. We also identify a class of theories whose covariance does not have neither a quantum nor an "almost quantum" origin, but which nevertheless produce the accessible two-point quantum mechanical correlators. Our approach leads to richer Bell-type inequalities in which the extent of nonlocality is intimately related to a non-additive entropic measure. In particular, it suggests that the Tsallis entropy with parameter q=1/2 is a natural operational measure of non-classicality. Moreover, when~generalizing this covariance matrix we find novel characterizations of the quantum mechanical set of correlators in multipartite scenarios. All these predictions might be experimentally validated when adding weak measurements to the conventional Bell test (without adding postselection).


title: The QCD crossover at zero and non-zero baryon densities from Lattice QCD

abstract: We map out the QCD crossover line $\frac{T_c(\mu_B)}{T_c(0)} = 1 - \kappa_2 \left( \frac{\mu_B}{T_c(0)} \right)^2 - \kappa_4 \left( \frac{\mu_B}{T_c(0)} \right)^4 + \mathcal{O}(\mu_B^6)$ for the first time up to $\mathcal{O}(\mu_B^4)$ for a strangeness neutral system by performing a Taylor expansion of chiral observables in temperature $T$ and chemical potentials $\mu$. At vanishing chemical potential, we report a crossover temperature $T_c(0) = (156.5 \pm 1.5)\;\mathrm{MeV}$ defined by the average of several second-order chiral susceptibilities. For a system with thermal conditions appropriate for a heavy-ion collision, we determined a curvature from the subtracted condensate as $\kappa_2 = 0.0120(20)$ and from the disconnected susceptibility as $\kappa_2 = 0.0123(30)$. The next order $\kappa_4$ is significantly smaller. We also report the crossover temperature as a function of the chemical potentials for: baryon-number, electric charge, strangeness and isospin. Additionally, we find that $T_c(\mu_B)$ is in agreement with lines of constant energy density and constant entropy density. Along this crossover line, we study net baryon-number fluctuations and show that their increase is substantially smaller compared to that obtained in HRG model calculations. Similarly, we analyze chiral susceptibility fluctuations along the crossover line and show that these are constant. We conclude that no signs for a narrowing of the crossover region can be found for baryon chemical potential $\mu_B < 250\;\mathrm{MeV}$.


title: Phonon-Polaritons in Nanoscale Waveguides

abstract: We investigate the phenomena of Brillouin induced opacity in nanoscale linear waveguides and Brillouin induced transparency in nanoscale ring waveguides. The concept of phonon-polariton is required in order to get a deep understanding of such phenomena. Phonon-polaritons are coherent superposition of photons and phonons obtained in the strong coupling regime induced through stimulated Brillouin scattering in the presence of a pump field. The formation of phonon-polaritons can be observed via linear optical spectra that are calculated here using the input-output formalism. For a linear waveguide the signature of phonon-polaritons appears as two peaks in the transmission spectrum with an opacity gap, and for a ring waveguide as two dips in the transmission spectrum with a transparency gap.


title: The Mathematics of Changing one's Mind, via Jeffrey's or via Pearl's update rule

abstract: Evidence in probabilistic reasoning may be 'hard' or 'soft', that is, it may be of yes/no form, or it may involve a strength of belief, in the unit interval [0, 1]. Reasoning with soft, [0, 1]-valued evidence is important in many situations but may lead to different, confusing interpretations. This paper intends to bring more mathematical and conceptual clarity to the field by shifting the existing focus from specification of soft evidence to accomodation of soft evidence. There are two main approaches, known as Jeffrey's rule and Pearl's method; they give different outcomes on soft evidence. This paper argues that they can be understood as correction and as improvement. It describes these two approaches as different ways of updating with soft evidence, highlighting their differences, similarities and applications. This account is based on a novel channel-based approach to Bayesian probability. Proper understanding of these two update mechanisms is highly relevant for inference, decision tools and probabilistic programming languages.


title: Thermodynamic Capacity of Quantum Processes

abstract: Thermodynamics imposes restrictions on what state transformations are possible. In the macroscopic limit of asymptotically many independent copies of a state---as for instance in the case of an ideal gas---the possible transformations become reversible and are fully characterized by the free energy. In this Letter, we present a thermodynamic resource theory for quantum processes that also becomes reversible in the macroscopic limit, a property that is especially rare for a resource theory of quantum channels. We identify a unique single-letter and additive quantity, the thermodynamic capacity, that characterizes the "thermodynamic value" of a quantum channel, in the sense that the work required to simulate many repetitions of a quantum process employing many repetitions of another quantum process becomes equal to the difference of the respective thermodynamic capacities. On a technical level, we provide asymptotically optimal constructions of universal implementations of quantum processes. A challenging aspect of this construction is the apparent necessity to coherently combine thermal engines that would run in different thermodynamic regimes depending on the input state. Our results have applications in quantum Shannon theory by providing a generalized notion of quantum typical subspaces and by giving an operational interpretation to the entropy difference of a channel.


title: Magnetosphere structure of a Kerr black hole: marginally force-free equatorial boundary condition

abstract: The role of equatorial boundary condition in the structure of a force-free black hole magnetosphere was rarely discussed, since previous studies have been focused on the field lines entering the horizon. However, recent high-accuracy force-free electrodynamics (FFE) simulations \cite{East2018} show that there are both field lines entering the horizon and field lines ending up on the equatorial current sheet within the ergosphere for asymptotic uniform field configuration. For the latter field lines, the equatorial boundary condition is well approximated being marginally force-free, i.e., $B^2-E^2\approx 0$, where $B$ and $E$ are the magnetic and electric field strength, respectively. In this paper, we revisit the uniform field solution to the Kerr BH magnetosphere structure and investigate the role of the marginally force-free equatorial boundary condition. We find this boundary condition plays an important role in shaping the BH magnetosphere in various aspects, including the shape of the light surface, the near-horizon field line configuration and the source of the Poynting flux. We also propose an algorithm for numerically solving the Grad-Shafranov equation and self-consistently imposing the marginally force-free equatorial condition. As a result, we find a good agreement between our numerical solutions and the high-accuracy FFE simulations. We also discuss the applicability of the marginally force-free boundary condition and the numerical algorithm proposed in this paper for general magnetic field configurations.


title: LiZn$_{2}$V$_{3}$O$_{8}$: A new geometrically frustrated cluster spin-glass

abstract: We have investigated the structural and magnetic properties of a new cubic spinel LiZn$_{2}$V$_{3}$O$_{8}$ (LZVO) through x-ray diffraction, dc and ac susceptibility, magnetic relaxation, aging, memory effect, heat capacity and $^{7}$Li nuclear magnetic resonance (NMR) measurements. A Curie-Weiss fit of the dc susceptibility $\chi_{\mathrm{dc}}$($\mathit{T}$) yields a Curie-Weiss temperature $\mathrm{\theta}_{\mathrm{CW}}$ = -185 K. This suggests strong antiferromagnetic (AFM) interactions among the magnetic vanadium ions. The dc and ac susceptibility data indicate the spin-glass behavior below a freezing temperature $T_{f}$ $\simeq$ 3 K. The frequency dependence of the $T_{f}$ is characterized by the Vogel-Fulcher law and critical dynamic scaling behavior or power law. From both fitting, we obtained the value of the characteristic angular frequency $\omega_{0}$ $\approx$ 3.56$\times$10$^{6}$ Hz, the dynamic exponent $\mathit{zv}$ $\approx$ 2.65, and the critical time constant $\tau_{0}$ $\approx$ 1.82$\times$10$^{-6}$ s, which falls in the conventional range for typical cluster spin-glass (CSG) systems. The value of relative shift in freezing temperature $\delta T_{f}$ $\simeq$ 0.039 supports a CSG ground states. We also found aging phenomena and memory effects in LZVO. The asymmetric response of the magnetic relaxation below $T_{f}$ supports the hierarchical model. Heat capacity data show no long-range or short-range ordering down to 2 K. Only about 25% magnetic entropy change $(\Delta S_{\mathrm{m}})$ signifies the presence of strong frustration in the system. The $^{7}$Li NMR spectra show a shift and broadening with decreasing temperature. The spin-lattice and spin-spin relaxation rates show anomalies due to spin freezing around 3 K as the bulk magnetization.


title: Numerical solution for Kapitza waves on a thin liquid film

abstract: The flow of a liquid film over an inclined plane is frequently found in nature and industry, and, under some conditions, instabilities in the free surface may appear. These instabilities are initially two-dimensional surface waves, known as Kapitza waves. Surface waves are important to many industrial applications. For example, liquid films with surface waves are employed to remove heat from solid surfaces. The initial phase of the instability is governed by the Orr-Sommerfeld equation and the appropriate boundary conditions; therefore, the fast and accurate solution of this equation is useful for industry. This paper presents a spectral method to solve the Orr-Sommerfeld equation with free surface boundary conditions. Our numerical approach is based on a Galerkin method with Chebyshev polynomials of the first kind, making it possible to express the Orr-Sommerfeld equation and their boundary conditions as a generalized eigenvalue problem. The main advantages of the present spectral method when compared to others such as, for instance, spectral collocation, are its stability and its readiness in including the boundary conditions in the discretized equations. We compare our numerical results with analytical solutions based on Perturbation Methods, which are valid only for long wave instabilities, and show that the results agree in the region of validity of the long-wave hypothesis. Far from this region, our results are still valid. In addition, we compare our results with published experimental results and the agreement is very good. The method is stable, fast, and capable to solve initial instabilities in free surface flows.


title: ANN-Benchmarks: A Benchmarking Tool for Approximate Nearest Neighbor Algorithms

abstract: This paper describes ANN-Benchmarks, a tool for evaluating the performance of in-memory approximate nearest neighbor algorithms. It provides a standard interface for measuring the performance and quality achieved by nearest neighbor algorithms on different standard data sets. It supports several different ways of integrating $k$-NN algorithms, and its configuration system automatically tests a range of parameter settings for each algorithm. Algorithms are compared with respect to many different (approximate) quality measures, and adding more is easy and fast; the included plotting front-ends can visualise these as images, $\LaTeX$ plots, and websites with interactive plots. ANN-Benchmarks aims to provide a constantly updated overview of the current state of the art of $k$-NN algorithms. In the short term, this overview allows users to choose the correct $k$-NN algorithm and parameters for their similarity search task; in the longer term, algorithm designers will be able to use this overview to test and refine automatic parameter tuning. The paper gives an overview of the system, evaluates the results of the benchmark, and points out directions for future work. Interestingly, very different approaches to $k$-NN search yield comparable quality-performance trade-offs. The system is available at http://ann-benchmarks.com .


title: Efficiently Secure Broadcasting in 5G Wireless Fog-Based-Fronthaul Networks

abstract: Enhanced Diversity and Network Coding (eDC-NC), the synergistic combination of Diversity and modified Triangular Network Coding, was introduced recently to provide efficient and ultra-reliable networking with near-instantaneous fault recovery. In this paper it is shown that eDC-NC technology can efficiently and securely broadcast messages in 5G wireless fog-computing-based Radio Access Networks (F-RAN). In particular, this work is directed towards demonstrating the ability of eDC-NC technology to more efficiently provide secure messages broadcasting than standardized methods such as Secure Multicasting using Secret (Shared) Key Cryptography, such that the adversary has no ability to acquire information even if they wiretap the entire F-RAN network (except of course the source and destination nodes). Our results show that using secure eDC-NC technology in F-RAN fronthaul network will enhance secure broadcasting and provide ultra-reliability networking, near-instantaneous fault recovery, and retain the throughput benefits of Network Coding.


title: Network neuroscience for optimizing brain-computer interfaces

abstract: Human-machine interactions are being increasingly explored to create alternative ways of communication and to improve our daily life. Based on a classification of the user's intention from the user's underlying neural activity, brain-computer interfaces (BCIs) allow direct interactions with the external environment while bypassing the traditional effector of the musculoskeletal system. Despite the enormous potential of BCIs, there are still a number of challenges that limit their societal impact, ranging from the correct decoding of a human's thoughts, to the application of effective learning strategies. Despite several important engineering advances, the basic neuroscience behind these challenges remains poorly explored. Indeed, BCIs involve complex dynamic changes related to neural plasticity at a diverse range of spatiotemporal scales. One promising antidote to this complexity lies in network science, which provides a natural language in which to model the organizational principles of brain architecture and function as manifest in its interconnectivity. Here, we briefly review the main limitations currently affecting BCIs, and we offer our perspective on how they can be addressed by means of network theoretic approaches. We posit that the emerging field of network neuroscience will prove to be an effective tool to unlock human-machine interactions.


title: Atom2Vec: learning atoms for materials discovery

abstract: Exciting advances have been made in artificial intelligence (AI) during the past decades. Among them, applications of machine learning (ML) and deep learning techniques brought human-competitive performances in various tasks of fields, including image recognition, speech recognition and natural language understanding. Even in Go, the ancient game of profound complexity, the AI player already beat human world champions convincingly with and without learning from human. In this work, we show that our unsupervised machines (Atom2Vec) can learn the basic properties of atoms by themselves from the extensive database of known compounds and materials. These learned properties are represented in terms of high dimensional vectors, and clustering of atoms in vector space classifies them into meaningful groups in consistent with human knowledge. We use the atom vectors as basic input units for neural networks and other ML models designed and trained to predict materials properties, which demonstrate significant accuracy.


title: Improved Person Re-Identification Based on Saliency and Semantic Parsing with Deep Neural Network Models

abstract: Given a video or an image of a person acquired from a camera, person re-identification is the process of retrieving all instances of the same person from videos or images taken from a different camera with non-overlapping view. This task has applications in various fields, such as surveillance, forensics, robotics, multimedia. In this paper, we present a novel framework, named Saliency-Semantic Parsing Re-Identification (SSP-ReID), for taking advantage of the capabilities of both clues: saliency and semantic parsing maps, to guide a backbone convolutional neural network (CNN) to learn complementary representations that improves the results over the original backbones. The insight of fusing multiple clues is based on specific scenarios in which one response is better than another, thus favoring the combination of them to increase performance. Due to its definition, our framework can be easily applied to a wide variety of networks and, in contrast to other competitive methods, our training process follows simple and standard protocols. We present extensive evaluation of our approach through five backbones and three benchmarks. Experimental results demonstrate the effectiveness of our person re-identification framework. In addition, we combine our framework with re-ranking techniques to achieve state-of-the-art results on three benchmarks.


title: Caching at the Edge with Fountain Codes

abstract: We address the use of linear randon fountain codes caching schemes in a heterogeneous satellite network. We consider a system composed of multiple hubs and a geostationary Earth orbit satellite. Coded content is memorized in hubs' caches in order to serve immediately the user requests and reduce the usage of the satellite backhaul link. We derive the analytical expression of the average backhaul rate, as well as a tight upper bound to it with a simple expression. Furthermore, we derive the optimal caching strategy which minimizes the average backhaul rate and compare the performance of the linear random fountain code scheme to that of a scheme using maximum distance separable codes. Our simulation results indicate that the performance obtained using fountain codes is similar to that of maximum distance separable codes.


title: NEUZZ: Efficient Fuzzing with Neural Program Smoothing

abstract: Fuzzing has become the de facto standard technique for finding software vulnerabilities. However, even state-of-the-art fuzzers are not very efficient at finding hard-to-trigger software bugs. Most popular fuzzers use evolutionary guidance to generate inputs that can trigger different bugs. Such evolutionary algorithms, while fast and simple to implement, often get stuck in fruitless sequences of random mutations. Gradient-guided optimization presents a promising alternative to evolutionary guidance. Gradient-guided techniques have been shown to significantly outperform evolutionary algorithms at solving high-dimensional structured optimization problems in domains like machine learning by efficiently utilizing gradients or higher-order derivatives of the underlying function. However, gradient-guided approaches are not directly applicable to fuzzing as real-world program behaviors contain many discontinuities, plateaus, and ridges where the gradient-based methods often get stuck. We observe that this problem can be addressed by creating a smooth surrogate function approximating the discrete branching behavior of target program. In this paper, we propose a novel program smoothing technique using surrogate neural network models that can incrementally learn smooth approximations of a complex, real-world program's branching behaviors. We further demonstrate that such neural network models can be used together with gradient-guided input generation schemes to significantly improve the fuzzing efficiency. Our extensive evaluations demonstrate that NEUZZ significantly outperforms 10 state-of-the-art graybox fuzzers on 10 real-world programs both at finding new bugs and achieving higher edge coverage. NEUZZ found 31 unknown bugs that other fuzzers failed to find in 10 real world programs and achieved 3X more edge coverage than all of the tested graybox fuzzers for 24 hours running.


title: Constraints on neutrino speed, weak equivalence principle violation, Lorentz invariance violation, and dual lensing from the first high-energy astrophysical neutrino source TXS 0506+056

abstract: We derive stringent constraints on neutrino speed, weak equivalence principle violation, Lorentz invariance violation, and dual lensing from the first high-energy astrophysical neutrino source: TXS 0506+056. Observation of neutrino (IceCube-170922A) and photons in a similar time frame and from the same direction is used to derive these limits. We describe ways in which these constraints can be further improved by orders of magnitude.

categories: astro-ph.HE astro-ph.CO hep-ex hep-ph hep-th


title: Patient-specific and physiological load sustaining synthetic graft substitutes for fusion of critically sized segmental bone defects

abstract: Critically sized defects are currently treated via autologous or allograft bone grafting, distraction osteogenesis, and membrane induction. However, these methods have major weaknesses which encourage the development of synthetic bone graft substitutes. A bioresorbable and physiologically load sustaining graft substitute was fabricated from polylactic acid using 3D printing, a versatile method which enables parameters to be customizable in a patient-specific manner. The internal architecture consists of vertical and horizontal conduits that are organized to achieve gradients in porosity and pore size in the radial direction to emulate cortical bone encapsulating cancellous bone. A diameter of 30mm and height of 10mm was designed to target the repair of critically sized long bone defects. The compressive properties of these graft substitutes were determined to be comparable to adult human femurs. Furthermore, dynamic fatigue testing at twice the human torso weight (800N) indicated that the bone graft substitute achieved stability as quickly as 200,000 cycles, suggesting its ability to support weight bearing loads while it gradually degrades as bone cells proliferate and mineralize. Additionally, the graft substitutes were seeded with human fetal osteoblasts to validate homogenous cellular attachment and cytocompatibility. Graft substitutes were also dynamically cultured for 28 days under various loading regimens using human mesenchymal stem cell differentiated osteoblasts, and these resulting tissue constructs were subjected to fatigue testing. The frequency of loading during dynamic culturing had a profound effect on the rate of mineralization which was particularly high at 2Hz, confirmed by Alizarin Red staining of the calcium deposits. Taken together, these bone graft substitutes have potential for use as orthopaedic implants to facilitate the repair of critically sized bone defects.


title: Accessible computational materials design with high fidelity and high throughput

abstract: Despite multiple successful applications of high-throughput computational materials design from first principles, there is a number of factors that inhibit its future adoption. Of particular importance are limited ability to provide high fidelity in a reliable manner and limited accessibility to non-expert users. We present example applications of a novel approach, where high-fidelity first-principles simulation techniques, Density Functional Theory with Hybrid Screened Exchange (HSE) and GW approximation, are standardized and made available online in an accessible and repeatable setting. We apply this approach to extract electronic band gaps and band structures for a diverse set of 71 materials ranging from pure elements to III-V and II-VI compounds, ternary oxides and alloys. We find that for HSE and G0W0, the average relative error fits within 20%, whereas for conventional Generalized Gradient Approximation the error is 55%. For HSE we find the average calculation time on an up-to-date server centrally available from a public cloud provider to fit within 48 hours. This work provides a cost-effective, accessible and repeatable practical recipe for performing high-fidelity first-principles calculations of electronic materials in a high-throughput manner.

categories: cond-mat.mtrl-sci cond-mat.other cs.CE cs.DC physics.comp-ph


title: The conjecture of Ulam on the invariance of measure on Hilbert cube

abstract: A conjecture of Ulam states that the standard probability measure $\pi$ on the Hilbert cube $I^\omega$ is invariant under the induced metric $d_a$ when the sequence $a = \{ a_i \}$ of positive numbers satisfies the condition $\sum\limits_{i=1}^\infty a_i^2 < \infty$. This conjecture was proved in \cite{jung1} when $E_1$ is a non-degenerate subset of $M_a$. In this paper, we prove the conjecture of Ulam completely by classifying cylinders as non-degenerate and degenerate cylinders and by treating the degenerate case that was overlooked in the previous paper.


title: Convex Bodies Associated to Tensor Norms

abstract: We determine when a convex body in $\mathbb{R}^d$ is the closed unit ball of a reasonable crossnorm on $\mathbb{R}^{d_1}\otimes\cdots\otimes\mathbb{R}^{d_l},$ $d=d_1\cdots d_l.$ We call these convex bodies "tensorial bodies". We prove that, among them, the only ellipsoids are the closed unit balls of Hilbert tensor products of Euclidean spaces. It is also proved that linear isomorphisms on $\mathbb{R}^{d_1}\otimes\cdots \otimes \mathbb{R}^{d_l}$ preserving decomposable vectors map tensorial bodies into tensorial bodies. This leads us to define a Banach-Mazur type distance between them, and to prove that there exists a Banach-Mazur type compactum of tensorial bodies.


title: Consensus Needs Broadcast in Noiseless Models but can be Exponentially Easier in the Presence of Noise

abstract: Consensus and Broadcast are two fundamental problems in distributed computing, whose solutions have several applications. Intuitively, Consensus should be no harder than Broadcast, and this can be rigorously established in several models. Can Consensus be easier than Broadcast? In models that allow noiseless communication, we prove a reduction of (a suitable variant of) Broadcast to binary Consensus, that preserves the communication model and all complexity parameters such as randomness, number of rounds, communication per round, etc., while there is a loss in the success probability of the protocol. Using this reduction, we get, among other applications, the first logarithmic lower bound on the number of rounds needed to achieve Consensus in the uniform GOSSIP model on the complete graph. The lower bound is tight and, in this model, Consensus and Broadcast are equivalent. We then turn to distributed models with noisy communication channels that have been studied in the context of some bio-inspired systems. In such models, only one noisy bit is exchanged when a communication channel is established between two nodes, and so one cannot easily simulate a noiseless protocol by using error-correcting codes. An $\Omega(\epsilon^{-2} n)$ lower bound on the number of rounds needed for Broadcast is proved by Boczkowski et al. [PLOS Comp. Bio. 2018] in one such model (noisy uniform PULL, where $\epsilon$ is a parameter that measures the amount of noise). In such model, we prove a new $\Theta(\epsilon^{-2} n \log n)$ bound for Broadcast and a $\Theta(\epsilon^{-2} \log n)$ bound for binary Consensus, thus establishing an exponential gap between the number of rounds necessary for Consensus versus Broadcast.


title: On triangular paperfolding patterns

abstract: We introduce patterns on a triangular grid generated by paperfolding operations. We show that in case these patterns are defined using a periodic sequence of foldings, they can also be generated using substitution rules and compute eigenvalues and eigenvectors of the corresponding matrices. We also prove that densities of all basic triangles are equal in these patterns.


title: Smart Charging and Parking of Plug-in Hybrid Electric Vehicles in Microgrids Considering Renewable Energy Sources

abstract: With the current trend of transforming a centralized power system into a decentralized one for efficiency, reliability, and environment reasons, the concept of microgrid that integrates a variety of distributed energy resources (DERs) on the distribution network is gaining popularity. In this paper, we investigate the smart charging and parking of plug-in hybrid electric vehicles (PHEVs) in microgrids with renewable energy sources (RES), such as solar panels, in grid-connected mode. To address the uncertainties associated with RES power output and PHEVs charging condition in the microgrid, we propose a two-stage scenario-based stochastic optimization approach with the objective of providing a proper scheduling for parking and charging of PHEVs that minimizes the average total operating cost while maintaining the reliability of the microgrid. A case study is conducted to show the effectiveness of the proposed approach. Extensive simulation results show that the microgrid can minimize the operating cost and ensure its reliability.


title: Very small insects use novel wing flapping and drag principle to generate the weight-supporting vertical force

abstract: The effect of air viscosity on the flow around an insect wing increases as insect size decreases. For the smallest insects (wing length R below 1 mm), the viscous effect is so large that lift-generation mechanisms used by their larger counterparts become ineffective. How the weight-supporting vertical force is generated is unknown. To elucidate the aerodynamic mechanisms responsible, we measure the wing kinematics of the tiny wasp Encarsia formosa (0.6 mm R) in hovering or very slow ascending flight and compute and analyze the aerodynamic forces. We find that the insects perform two unusual wing-motions. One is rowing: the wings move fast downward and backward, like stroking oars; the other is the previously discovered Weis-Fogh fling. The rowing produces 70 percent of the required vertical force and the Weis-Fogh fling the other 30 percent. The oaring wing mainly produces an approximately up-pointing drag, resulting in the vertical force. Because each oaring produces a starting flow, the drag is unsteady in nature and much greater than that in steady motion at the same velocities and angles of attack. Furthermore, our computation shows that if the tiny wasps employed the usual wing kinematics of the larger insects (flapping back and forth in a horizontal plane), vertical force produced would be only one third of that by the real wing kinematics; i.e. they must use the special wing movements to overcome the problem of large viscous effects encountered by the commonly used flapping kinematics. We for the first time observe very small insects using drag to support their weight and explain how a net vertical force is generated when the drag principle is applied.


title: Partially smoothed information measures

abstract: Smooth entropies are a tool for quantifying resource trade-offs in (quantum) information theory and cryptography. In typical bi- and multi-partite problems, however, some of the sub-systems are often left unchanged and this is not reflected by the standard smoothing of information measures over a ball of close states. We propose to smooth instead only over a ball of close states which also have some of the reduced states on the relevant sub-systems fixed. This partial smoothing of information measures naturally allows to give more refined characterizations of various information-theoretic problems in the one-shot setting. In particular, we immediately get asymptotic second-order characterizations for tasks such as privacy amplification against classical side information or classical state splitting. For quantum problems like state merging the general resource trade-off is tightly characterized by partially smoothed information measures as well.


title: Joint Modeling and Optimization of Search and Recommendation

abstract: Despite the somewhat different techniques used in developing search engines and recommender systems, they both follow the same goal: helping people to get the information they need at the right time. Due to this common goal, search and recommendation models can potentially benefit from each other. The recent advances in neural network technologies make them effective and easily extendable for various tasks, including retrieval and recommendation. This raises the possibility of jointly modeling and optimizing search ranking and recommendation algorithms, with potential benefits to both. In this paper, we present theoretical and practical reasons to motivate joint modeling of search and recommendation as a research direction. We propose a general framework that simultaneously learns a retrieval model and a recommendation model by optimizing a joint loss function. Our preliminary results on a dataset of product data indicate that the proposed joint modeling substantially outperforms the retrieval and recommendation models trained independently. We list a number of future directions for this line of research that can potentially lead to development of state-of-the-art search and recommendation models.


title: Hadronic observables in p+p and d+Au collisions at RHIC using CGC+PYTHIA

abstract: The IP-Glasma model of CGC combined with the Lund model of PYTHIA provides a very successful description of hadron production from gluon dominated non-equilibrium matter in various small collision systems. This new CGC+PYTHIA framework, naturally reproduces several characteristic features of the hadronic observables such as the mass ordering of $v_{2}(p_{T})$ and $\langle p_{T}\rangle$, often ascribed to collectivity driven by hydrodynamics at the LHC. In this contribution we extend our work to provide a systematic comparison of particle spectra and multiplicity distributions in p+p and d+Au collisions at the RHIC.


title: A central limit theorem for star-generators of $S_{\infty}$, which relates to the law of a GUE matrix

abstract: It is well-known that, on a purely algebraic level, a simplified algebraic version of the Central Limit Theorem (CLT) can be proved in the framework of a noncommutative probability space, under the hypotheses that the sequence of non-commutative random variables we consider is exchangeable and obeys a certain vanishing condition of some of its joint moments. In this approach (which covers versions for both the classical CLT and the CLT of free probability), the determination of the resulting limit law has to be addressed on a case-by-case basis. In this paper we discuss an instance of the above theorem which takes place in the framework of the group algebra of the infinite symmetric group $S_{\infty}$: the exchangeable sequence that is considered consists of the star-generators of $S_{\infty}$, and the expectation functional used on the group algebra of $S_{\infty}$ depends in a natural way on a parameter $d$, which is a positive integer. We identify precisely the limit distribution $\mu_d$ for this special instance of exchangeable CLT, via a connection that $\mu_d$ turns out to have with the average empirical eigenvalue distribution of a random GUE matrix of size $d \times d$. Moreover, we put into evidence a multi-variate version of this result which follows from the observation that, on the level of calculations with pair-partitions, the (non-centred) star-generators are related to a (centred) exchangeable sequence of GUE matrices with independent entries


title: Stabilized CutDG methods for advection-reaction problems

abstract: We develop novel stabilized cut discontinuous Galerkin (CutDG) methods for advection-reaction problems. The domain of interest is embedded into a structured, unfitted background mesh in $\mathbb{R}^d$ where the domain boundary can cut through the mesh in an arbitrary fashion. To cope with robustness problems caused by small cut elements, we introduce ghost penalties in the vicinity of the embedded boundary to stabilize certain (semi)-norms associated with the advection and reaction operator. A few abstract assumptions on the ghost penalties are identified enabling us to derive geometrically robust and optimal a priori error and condition number estimates for the stationary advection-reaction problem which hold irrespective of the particular cut configuration. Possible realizations of suitable ghost penalties are discussed. The theoretical results are corroborated by a number of computational studies for various approximation orders and for two and three-dimensional test problems.


title: What is the Wigner function closest to a given square integrable function?

abstract: We consider an arbitrary square integrable function $F$ on the phase space and look for the Wigner function closest to it with respect to the $L^2$ norm. It is well known that the minimizing solution is the Wigner function of any eigenvector associated with the largest eigenvalue of the Hilbert-Schmidt operator with Weyl symbol $F$. We solve the particular case of radial functions on the two-dimensional phase space exactly. For more general cases, one has to solve an infinite dimensional eigenvalue problem. To avoid this difficulty, we consider a finite dimensional approximation and estimate the errors for the eigenvalues and eigenvectors. As an application, we address the so-called Wigner approximation suggested by some of us for the propagation of a pulse in a general dispersive medium. We prove that this approximation never leads to a {\it bona fide} Wigner function. This is our prime motivation for our optimization problem. As a by-product of our results we are able to estimate the eigenvalues and Schatten norms of certain Schatten-class operators. The techniques presented here may be potentially interesting for estimating eigenvalues of localization operators in time-frequency analysis and quantum mechanics.


title: Cross Pixel Optical Flow Similarity for Self-Supervised Learning

abstract: We propose a novel method for learning convolutional neural image representations without manual supervision. We use motion cues in the form of optical flow, to supervise representations of static images. The obvious approach of training a network to predict flow from a single image can be needlessly difficult due to intrinsic ambiguities in this prediction task. We instead propose a much simpler learning goal: embed pixels such that the similarity between their embeddings matches that between their optical flow vectors. At test time, the learned deep network can be used without access to video or flow information and transferred to tasks such as image classification, detection, and segmentation. Our method, which significantly simplifies previous attempts at using motion for self-supervision, achieves state-of-the-art results in self-supervision using motion cues, competitive results for self-supervision in general, and is overall state of the art in self-supervised pretraining for semantic image segmentation, as demonstrated on standard benchmarks.


title: Inference of the Dynamic Aging-related Biological Subnetwork via Network Propagation

abstract: Gene expression (GE) data capture valuable condition-specific information ("condition" can mean a biological process, disease stage, age, patient, etc.) However, GE analyses ignore physical interactions between gene products, i.e., proteins. Since proteins function by interacting with each other, and since biological networks (BNs) capture these interactions, BN analyses are promising. However, current BN data fail to capture condition-specific information. Recently, GE and BN data have been integrated using network propagation (NP) to infer condition-specific BNs. However, existing NP-based studies result in a static condition-specific subnetwork, even though cellular processes are dynamic. A dynamic process of our interest is human aging. We use prominent existing NP methods in a new task of inferring a dynamic rather than static condition-specific (aging-related) subnetwork. Then, we study evolution of network structure with age - we identify proteins whose network positions significantly change with age and predict them as new aging-related candidates. We validate the predictions via e.g., functional enrichment analyses and literature search. Dynamic network inference via NP yields higher prediction quality than the only existing method for inferring a dynamic aging-related BN, which does not use NP.


title: Three-dimensional Stable Matching with Cyclic Preferences

abstract: We consider the three-dimensional stable matching problem with cyclic preferences, a problem originally proposed by Knuth. Despite extensive study of the problem by experts from different areas, the question of whether every instance of this problem admits a stable matching remains unanswered. In 2004, Boros, Gurvich, Jaslar and Krasner showed that a stable matching always exists when the number of agents in each of the groups is three. In 2006, Eriksson, Sj\"ostrand and Strimling showed that a stable matching exists also when the number of agents in each group is four. In this paper, we demonstrate that a stable matching exists when each group has five agents. Furthermore, we show that there are at least two distinct stable matchings in that setting.


title: Calibrating magnification bias for the $E_G$ statistic to test general relativity

abstract: We assess the effect of magnification bias on the $E_G$ statistic for probing gravity. $E_G$, a statistic constructed from power spectrum estimates of both weak lensing and redshift space distortions (RSD), directly tests general relativity (GR) while in principle being independent of clustering bias. This property has motivated its recent use in multiple tests of GR. Recent work has suggested that the magnification bias of galaxies due to foreground matter perturbations breaks the bias-independence of the $E_G$ statistic. The magnitude of this effect is very sensitive to the clustering and magnification biases of the galaxy sample. We show that for realistic values of the clustering and magnification biases, the effect for magnification bias is small relative to statistical errors for most spectroscopic galaxy surveys but large for photometric galaxy surveys. For the cases with significant magnification bias, we propose a method to calibrate magnification bias using measurements of the lensing auto-power spectrum. We test this calibration method using simulations, finding that our calibration method can calibrate $E_G$ from 2-4 times the simulation error to well within the errors, although its accuracy is sensitive to the precision of the measured redshift distribution and magnification bias, but not the clustering bias. This work gives strong evidence that this method will work increasingly well in future CMB lensing surveys.


title: Some open problems on holomorphic foliation theory

abstract: We present a list of open questions in the theory of holomorphic foliations, possibly with singularities. Some problems have been around for a while, others are very accessible.


title: The Consistency of Arithmetic

abstract: In 2010, Vladimir Voevodsky gave a lecture on "What If Current Foundations of Mathematics Are Inconsistent?" Among other things he said that he was seriously suspicious that an inconsistency in PA (first-order Peano arithmetic) might someday be found. About a year later, Edward Nelson announced that he had discovered an inconsistency not just in PA, but in a small fragment of primitive recursive arithmetic. Soon, Daniel Tausk and Terence Tao independently found a fatal error, and Nelson withdrew his claim, stating that consistency of PA was an "open problem." Many mathematicians may find such claims bewildering. Is the consistency of PA really an open problem? If so, would the discovery of an inconsistency in PA cause all of mathematics to come crashing down like a house of cards? This expository article attempts to address these questions, by sketching and discussing existing proofs of the consistency of PA (including Gentzen's proof and Friedman's relative consistency proof that appeals to the Bolzano-Weierstrass theorem). Since Nelson was a self-avowed formalist, the article also examines the implications of formalism.


title: LATE Ain'T Earley: A Faster Parallel Earley Parser

abstract: We present the LATE algorithm, an asynchronous variant of the Earley algorithm for parsing context-free grammars. The Earley algorithm is naturally task-based, but is difficult to parallelize because of dependencies between the tasks. We present the LATE algorithm, which uses additional data structures to maintain information about the state of the parse so that work items may be processed in any order. This property allows the LATE algorithm to be sped up using task parallelism. We show that the LATE algorithm can achieve a 120x speedup over the Earley algorithm on a natural language task.


title: Pairing of the P\"oschl-Teller gas

abstract: We report calculations of equation of state of a model system, representative for a cold Fermi gas, of particles interacting via the P\"oschl-Teller interaction. In successively more sophisticated calculations, we determine the importance of correlations and non-localities. We show that the gas displays, at relatively low density, an instability indicated by a divergence of the in-medium scattering length which occurs well before the divergence of the vacuum scattering length and the spinodal density. We also calculate the pairing gap and show that non-local correlations can change the pairing gap by almost an order of magnitude.

categories: cond-mat.quant-gas cond-mat.str-el cond-mat.supr-con nucl-th quant-ph


title: Nonstandard solutions for a perturbed nonlinear Schr\"{o}dinger system with small coupling coefficients\protect\thanks{A perturbed nonlinear Schr\"{o}dinger system

abstract: In this paper, we consider the following weakly coupled nonlinear Schr\"odinger system \begin{equation*} \left\{ \begin{array}{ll} -\epsilon^{2}\Delta u_1 + V_1(x)u_1 = |u_1|^{2p - 2}u_1 + \beta|u_1|^{p - 2}|u_2|^pu_1, & x\in \mathbb{R}^N,\\ -\epsilon^{2}\Delta u_2 + V_2(x)u_2 = |u_2|^{2p - 2}u_2 + \beta|u_2|^{p - 2}|u_1|^pu_2, & x\in \mathbb{R}^N, \end{array} \right. \end{equation*} where $\epsilon>0$, $\beta\in\mathbb{R}$ is a coupling constant, $2p\in (2,2^*)$ with $2^* = \frac{2N}{N - 2}$ if $N\geq 3$ and $+\infty$ if $N = 1,2$, $V_1$ and $V_2$ belong to $C(\mathbb{R}^N,[0,\infty))$. When $p\ge 2$ and $\beta>0$ is suitably small, we show that the problem has a family of nonstandard solutions $\{w_{\epsilon} = (u^1_{\epsilon},u^2_{\epsilon}):0<\epsilon<\epsilon_{0}\}$ concentrating synchronously at the common local minimum of $V_1$ and $V_2$. All decay rates of $V_i(i=1,2)$ are admissible and we can allow that $\beta>0$ is close to $0$ in this paper. Moreover, the location of concentration points is given by local Pohozaev identities. Our proofs are based on variational methods and the penalized technique.


title: Stable noncommutative polynomials and their determinantal representations

abstract: A noncommutative polynomial is stable if it is nonsingular on all tuples of matrices whose imaginary parts are positive definite. In this paper a characterization of stable polynomials is given in terms of strongly stable linear matrix pencils, i.e., pencils of the form $H+iP_0+P_1x_1+\cdots+P_dx_d$, where $H$ is hermitian and $P_j$ are positive semidefinite matrices. Namely, a noncommutative polynomial is stable if and only if it admits a determinantal representation with a strongly stable pencil. More generally, structure certificates for noncommutative stability are given for linear matrix pencils and noncommutative rational functions.


title: Nucleus-nucleus potential from identical-particle interference

abstract: Based on the quantum interference between two-identical-nucleus scattering at energies around the Coulomb barrier, the barrier positions for $^{58}$Ni+$^{58}$Ni and $^{16}$O+$^{16}$O are extracted from Mott oscillations in the angular distributions around 90$^{\circ}$ for the first time. The angle separation of pairs of Mott scattering valleys around 90$^{\circ}$ has a direct relationship with the closest distance between two nuclei in elastic scattering. Together with the barrier height from fusion excitation function, the extracted barrier position provides a sensitive probe to constrain the model predictions for the nucleus-nucleus potential barrier.


title: Control of Ultracold Atoms with a Chiral Ferromagnetic Film

abstract: We show that the magnetic field produced by a chiral ferromagnetic film can be applied to control ultracold atoms. The film will act as a magnetic mirror or a reflection grating for ultracold atoms when it is in the helical phase or the skyrmion crystal phase respectively. By applying a bias magnetic field and a time-dependent magnetic field, one-dimensional or two-dimensional magnetic lattices including honeycomb, Kagome, triangular types can be created to trap the ultracold atoms. We have also discussed the trapping height, potential barrier, trapping frequency, and Majorana loss rate for each lattice. Our results suggest that the chiral ferromagnetic film can be a platform to develop artificial quantum systems with ultracold atoms based on modern spintronics technologies.


title: Convolutional Sparse Kernel Network for Unsupervised Medical Image Analysis

abstract: The availability of large-scale annotated image datasets and recent advances in supervised deep learning methods enable the end-to-end derivation of representative image features that can impact a variety of image analysis problems. Such supervised approaches, however, are difficult to implement in the medical domain where large volumes of labelled data are difficult to obtain due to the complexity of manual annotation and inter- and intra-observer variability in label assignment. We propose a new convolutional sparse kernel network (CSKN), which is a hierarchical unsupervised feature learning framework that addresses the challenge of learning representative visual features in medical image analysis domains where there is a lack of annotated training data. Our framework has three contributions: (i) We extend kernel learning to identify and represent invariant features across image sub-patches in an unsupervised manner. (ii) We initialise our kernel learning with a layer-wise pre-training scheme that leverages the sparsity inherent in medical images to extract initial discriminative features. (iii) We adapt a multi-scale spatial pyramid pooling (SPP) framework to capture subtle geometric differences between learned visual features. We evaluated our framework in medical image retrieval and classification on three public datasets. Our results show that our CSKN had better accuracy when compared to other conventional unsupervised methods and comparable accuracy to methods that used state-of-the-art supervised convolutional neural networks (CNNs). Our findings indicate that our unsupervised CSKN provides an opportunity to leverage unannotated big data in medical imaging repositories.


title: Multiplicative Schr\"odinger problem and the Dirichlet transport

abstract: We consider an optimal transport problem on the unit simplex whose solutions are given by gradients of exponentially concave functions and prove two main results. First, we show that the optimal transport is the large deviation limit of a particle system of Dirichlet processes transporting one probability measure on the unit simplex to another by coordinatewise multiplication and normalizing. The structure of our Lagrangian and the appearance of the Dirichlet process relate our problem closely to the entropic measure on the Wasserstein space as defined by von-Renesse and Sturm in the context of Wasserstein diffusion. The limiting procedure is a triangular limit where we allow simultaneously the number of particles to grow to infinity while the `noise' tends to zero. The method, which generalizes easily to many other cost functions, including the squared Euclidean distance, provides a novel combination of the Schr\"odinger problem approach due to C. L\'eonard and the related Brownian particle systems by Adams et al.which does not require gamma convergence. Second, we analyze the behavior of entropy along the paths of transport. The reference measure on the simplex is taken to be the Dirichlet measure with all zero parameters which relates to the finite-dimensional distributions of the entropic measure. The interpolating curves are not the usual McCann lines. Nevertheless we show that entropy plus a multiple of the transport cost remains convex, which is reminiscent of the semiconvexity of entropy along lines of McCann interpolations in negative curvature spaces. We also obtain, under suitable conditions, dimension-free bounds of the optimal transport cost in terms of entropy.


title: Time Series Deinterleaving of DNS Traffic

abstract: Stream deinterleaving is an important problem with various applications in the cybersecurity domain. In this paper, we consider the specific problem of deinterleaving DNS data streams using machine-learning techniques, with the objective of automating the extraction of malware domain sequences. We first develop a generative model for user request generation and DNS stream interleaving. Based on these we evaluate various inference strategies for deinterleaving including augmented HMMs and LSTMs on synthetic datasets. Our results demonstrate that state-of-the-art LSTMs outperform more traditional augmented HMMs in this application domain.


title: Some model-theoretic results on the 3-valued paraconsistent first-order logic QCiore

abstract: In this paper the 3-valued paraconsistent first-order logic QCiore is studied from the point of view of Model Theory. The semantics for QCiore is given by partial structures, which are first-order structures in which each n-ary predicate R is interpreted as a triple of paiwise disjoint sets of n-uples representing, respectively, the set of tuples which actually belong to R, the set of tuples which actually do not belong to R, and the set of tuples whose status is dubious or contradictory. Partial structures were proposed in 1986 by I. Mikenberg, N. da Costa and R. Chuaqui for the theory of quasi-truth (or pragmatic truth). In 2014, partial structures were studied by M. Coniglio and L. Silvestrini for a 3-valued paraconsistent first-order logic called LPT1, whose 3-valued propositional fragment is equivalent to da Costa-D'Otaviano's logic J3. This approach is adapted in this paper to QCiore, and some important results of classical Model Theory such as Robinson's joint consistency theorem, amalgamation and interpolation are obtained. Although we focus on QCiore, this framework can be adapted to other 3-valued first-order logics.


title: CLEF: Limiting the Damage Caused by Large Flows in the Internet Core (Technical Report)

abstract: The detection of network flows that send excessive amounts of traffic is of increasing importance to enforce QoS and to counter DDoS attacks. Large-flow detection has been previously explored, but the proposed approaches can be used on high-capacity core routers only at the cost of significantly reduced accuracy, due to their otherwise too high memory and processing overhead. We propose CLEF, a new large-flow detection scheme with low memory requirements, which maintains high accuracy under the strict conditions of high-capacity core routers. We compare our scheme with previous proposals through extensive theoretical analysis, and with an evaluation based on worst-case-scenario attack traffic. We show that CLEF outperforms previously proposed systems in settings with limited memory.


title: Learning and Matching Multi-View Descriptors for Registration of Point Clouds

abstract: Critical to the registration of point clouds is the establishment of a set of accurate correspondences between points in 3D space. The correspondence problem is generally addressed by the design of discriminative 3D local descriptors on the one hand, and the development of robust matching strategies on the other hand. In this work, we first propose a multi-view local descriptor, which is learned from the images of multiple views, for the description of 3D keypoints. Then, we develop a robust matching approach, aiming at rejecting outlier matches based on the efficient inference via belief propagation on the defined graphical model. We have demonstrated the boost of our approaches to registration on the public scanning and multi-view stereo datasets. The superior performance has been verified by the intensive comparisons against a variety of descriptors and matching methods.


title: The spherical metric and univalent harmonic mappings

abstract: Let $f=h+\overline{g}$ be a harmonic univalent map in the unit disk $\mathbb{D}$, where $h $ and $g$ are analytic. We obtain an improved estimate for the second coefficient of $h$. This indeed is the first qualitative improvement after the appearance of the papers by Clunie and Sheil-Small in 1984, and by Sheil-Small in 1990. Also, when the sup-norm of the dilatation is less than $1$, it is shown that the spherical area of the covering surface of $h$ is dominated by the spherical area of the covering surface of $f.$


title: Unified framework for heavy flavor and quarkonium production in high multiplicity p+p and p+A collisions at RHIC and LHC

abstract: We discuss the production of $D$-mesons and $J/\psi$ in high multiplicity proton-proton and proton-nucleus collisions within the Color-Glass-Condensate (CGC) framework. We demonstrate that the modification of the LHC data on $D$ and $J/\psi$ yields in high multiplicity events relative to minimum bias events arise from a significant enhancement of the gluon saturation scales of the corresponding rare parton configurations in the colliding protons and nuclei. For a given event multiplicity, we predict these relative yields to be energy independent from $\sqrt{s}=200$ GeV at RHIC to the highest LHC energies.


title: Nonlocal multicontinua upscaling for multicontinua flow problems in fractured porous media

abstract: Our goal of this paper is to develop a new upscaling method for multicontinua flow problems in fractured porous media. We consider a system of equations that describes flow phenomena with multiple flow variables defined on both matrix and fractures. To construct our upscaled model, we will apply the nonlocal multicontinua (NLMC) upscaling technique. The upscaled coefficients are obtained by using some multiscale basis functions, which are solutions of local problems defined on oversampled regions. For each continuum within a target coarse element, we will solve a local problem defined on an oversampling region obtained by extending the target element by few coarse grid layers, with a set of constraints which enforce the local solution to have mean value one on the chosen continuum and zero mean otherwise. The resulting multiscale basis functions have been shown to have good approximation properties. To illustrate the idea of our approach, we will consider a dual continua background model consisting of discrete fractures in two space dimensions, that is, we consider a system with three continua. We will present several numerical examples, and they show that our method is able to capture the interaction between matrix continua and discrete fractures on the coarse grid efficiently.


title: Visualizing the Hidden Features of Galaxy Morphology with Machine Learning

abstract: We train three convolutional neural networks (CNNs) to classify galaxies with Galaxy Zoo 2 dataset and extract the activations from the last fully connected layer or the last average pooling layer of CNNs to study the high-dimensional abstract feature representations of galaxy images. We apply t-Distributed Stochastic Neighbour Embedding (t-SNE), a popular dimensionality reduction technique, to visualize the high-dimensional galaxy feature representations in two-dimensional scatter plots. From the visualization, we try to understand the galaxy images data itself and obtain some highly valuable insights. For instance, the learned galaxy feature representations from networks indicate that the galaxies belonging to the same class tend to group together, i.e. same morphological galaxies are clustered; The cluster of completely round smooth galaxy and the cluster of in-between smooth galaxy (between completely round and cigar-shaped) are moved closer, compared to other clusters; The cluster of cigar-shaped smooth galaxy and the cluster of edge-on galaxy are intertwined surprisingly; A galaxy mislabelled as spiral galaxy in the original dataset falls in the cluster of completely round smooth galaxy, and manual inspection also identifies out the outlier as a completely round smooth galaxy. These findings will facilitate the study of galaxy morphology.


title: On curves intersecting at most once

abstract: We prove that on a closed surface of genus $g$, the cardinality of a set of simple closed curves in which any two are non-homotopic and intersect at most once is $\lesssim g^2 \log(g)$. This bound matches the largest known constructions to within a logarithmic factor. The proof uses a probabilistic argument in graph theory. It generalizes as well to the case of curves that intersect at most $k$ times in pairs.


title: Time Asymmetry of Cosmic Background Evolution in Loop Quantum Cosmology

abstract: We discuss the asymmetry of cosmic background evolution in time with respect to the quantum bounce in the Loop Quantum Cosmology (LQC), employing the value of scalar field at the bounce $\phi_{\rm B}$. We use the Chaotic and the $R^2$ potentials to demonstrate that a possible deflation before the bounce may counteract the inflation that is needed for resolving the cosmological conundrums, so a certain level of time asymmetry is required for the models in LQC. This $\phi_{\rm B}$ is model dependent and closely related to the amounts of deflation and inflation, so we may use observations to confine $\phi_{\rm B}$ and thus the model parameters. With further studies this formalism should be useful in providing an observational testbed for the LQC models.


title: A High-Accuracy Adaptive Beam Training Algorithm for MmWave Communication

abstract: In millimeter wave communications, beam training is an effective way to achieve beam alignment. Traditional beam training method allocates training resources equally to each beam in the pre-designed beam training codebook. The performance of this method is far from satisfactory, because different beams have different beamforming gain, and thus some beams are relatively more difficult to be distinguished from the optimal beam than the others. In this paper, we pro- pose a new beam training algorithm which adaptively allocates training resources to each beam. Specifically, the proposed algorithm allocates more training symbols to the beams with relatively higher beamforming gain, while uses less resources to distinguish the beams with relatively lower beamforming gain. Through theoretical analysis and numerical simulations, we show that in practical situations the proposed adaptive algorithm asymptotically outperforms the traditional method in terms of beam training accuracy. Moreover, simulations also show that this relative performance behavior holds in non-asymptotic regime.


title: Non-Hermitian Hopf-Link Exceptional Line Semimetals

abstract: We study a new class of non-Hermitian topological phases in three dimension in the absence of any symmetry, where the topological robust band degeneracies are Hopf-link exceptional lines. As a concrete example, we investigate the non-Hermitian band structures of nodal line semimetals under non-Hermitian perturbations, where the Fermi surfaces can transit from 1d nodal lines to 2d twisting surfaces with Hopf-link boundaries when the winding number defined along the nodal line is $\pm 1$. The linking numbers of these linked exceptional line phases are also proposed, based on the integral of Chern-Simons form over the Brillouin zone.


title: Data Reduction in Markov model using EM algorithm

abstract: This paper describes a data reduction technique in case of a markov chain of specified order. Instead of observing all the transitions in a markov chain we record only a few of them and treat the remaining part as missing. The decision about which transitions to be filtered is taken before the observation process starts. Based on the filtered chain we try to estimate the parameters of the markov model using EM algorithm. In the first half of the paper we characterize a class of filtering mechanism for which all the parameters remain identifiable. In the later half we explain methods of estimation and testing about the transition probabilities of the markov chain based on the filtered data. The methods are first developed assuming a simple markov model with each probability of transition positive, but then generalized for models with structural zeroes in the transition probability matrix. Further extension is also done for multiple markov chains. The performance of the developed method of estimation is studied using simulated data along with a real life data.


title: Existence and regularity results for minimal sets; Plateau problem

abstract: Solving the Plateau problem means to find the surface with minimal area among all surfaces with a given boundary. Part of the problem actually consists of giving a suitable definition to the notions of 'surface', 'area' and 'boundary'. In our setting the considered objects are sets whose Hausdorff area is locally finite. The sliding boundary condition is given in term of a one parameter family of compact deformations which allows the boundary of the surface to moove along a closed set. The area functional is related to capillarity and free-boundary problems, and is a slight modification of the Hausdorff area. We focused on minimal boundary cones; that is to say tangent cones on boundary points of sliding minimal surfaces. In particular we studied cones contained in an half-space and whose boundary can slide along the bounding hyperplane. After giving a classification of one-dimensional minimal cones in the half-plane we provided four new two-dimensional minimal cones in the three-dimensional half space (which cannot be obtained as the Cartesian product of the real line with one of the previous cones). We employed the technique of paired calibrations and in one case could also generalise it to higher dimension. In order to prove that the provided list of minimal cones is complete, we started the classification of cones satisfying the necessary conditions for the minimality, and with numeric simulations we obtained better competitors for these new candidates


title: Evidence for breathing-type pseudo Jahn-Teller distortions in the charge density wave phase of 1$T$-TiSe$_2$

abstract: The charge density wave (CDW) phase in 1$T$-TiSe$_2$ is investigated using angle resolved photoemission spectroscopy (ARPES) and neutron scattering measurements. Our ARPES results reveal a clear temperature dependence of the chemical potential of the system. They also demonstrate specific changes encountered by the Se 4$p$ valence and Ti 3$d$ conduction bands as the temperature of the system is decreased through $T_{\text{CDW}}$. The valence band undergoes a downward shift, whereas the conduction band remains unaffected. The crystal structure in the CDW state shows a distinct split of the Ti-Se atomic correlations that are reminiscent of Jahn-Teller distortions, manifested in a breathing-type mode. The ARPES data together with the local structure analysis support a direct link between Jahn-Teller-like distortions and the CDW order in 1$T$-TiSe$_2$.


title: Improving the smoothed complexity of FLIP for max cut problems

abstract: Finding locally optimal solutions for max-cut and max-$k$-cut are well-known PLS-complete problems. An instinctive approach to finding such a locally optimum solution is the FLIP method. Even though FLIP requires exponential time in worst-case instances, it tends to terminate quickly in practical instances. To explain this discrepancy, the run-time of FLIP has been studied in the smoothed complexity framework. Etscheid and R\"{o}glin showed that the smoothed complexity of FLIP for max-cut in arbitrary graphs is quasi-polynomial. Angel, Bubeck, Peres, and Wei showed that the smoothed complexity of FLIP for max-cut in complete graphs is $O(\phi^5n^{15.1})$, where $\phi$ is an upper bound on the random edge-weight density and $n$ is the number of vertices in the input graph. While Angel et al.'s result showed the first polynomial smoothed complexity, they also conjectured that their run-time bound is far from optimal. In this work, we make substantial progress towards improving the run-time bound. We prove that the smoothed complexity of FLIP in complete graphs is $O(\phi n^{7.83})$. Our results are based on a carefully chosen matrix whose rank captures the run-time of the method along with improved rank bounds for this matrix and an improved union bound based on this matrix. In addition, our techniques provide a general framework for analyzing FLIP in the smoothed framework. We illustrate this general framework by showing that the smoothed complexity of FLIP for max-$3$-cut in complete graphs is polynomial and for max-$k$-cut in arbitrary graphs is quasi-polynomial. We believe that our techniques should also be of interest towards addressing the smoothed complexity of FLIP for max-$k$-cut in complete graphs for larger constants $k$.


title: Scene Learning: Deep Convolutional Networks For Wind Power Prediction by Embedding Turbines into Grid Space

abstract: Wind power prediction is of vital importance in wind power utilization. There have been a lot of researches based on the time series of the wind power or speed, but In fact, these time series cannot express the temporal and spatial changes of wind, which fundamentally hinders the advance of wind power prediction. In this paper, a new kind of feature that can describe the process of temporal and spatial variation is proposed, namely, Spatio-Temporal Features. We first map the data collected at each moment from the wind turbine to the plane to form the state map, namely, the scene, according to the relative positions. The scene time series over a period of time is a multi-channel image, i.e. the Spatio-Temporal Features. Based on the Spatio-Temporal Features, the deep convolutional network is applied to predict the wind power, achieving a far better accuracy than the existing methods. Compared with the starge-of-the-art method, the mean-square error (MSE) in our method is reduced by 49.83%, and the average time cost for training models can be shortened by a factor of more than 150.


title: A standard siren measurement of the Hubble constant from GW170817 without the electromagnetic counterpart

abstract: We perform a statistical standard siren analysis of GW170817. Our analysis does not utilize knowledge of NGC 4993 as the unique host galaxy of the optical counterpart to GW170817. Instead, we consider each galaxy within the GW170817 localization region as a potential host; combining the redshift from each galaxy with the distance estimate from GW170817 provides an estimate of the Hubble constant, $H_0$. We then combine the $H_0$ values from all the galaxies to provide a final measurement of $H_0$. We explore the dependence of our results on the thresholds by which galaxies are included in our sample, as well as the impact of weighting the galaxies by stellar mass and star-formation rate. Considering all galaxies brighter than $0.01 L^\star_B$ as equally likely to host a BNS merger, we find $H_0= 76^{+48}_{-23}$ km s$^{-1}$ Mpc$^{-1}$ (maximum a posteriori and 68.3% highest density posterior interval; assuming a flat $H_0$ prior in the range $\left[ 10, 220 \right]$ km s$^{-1}$ Mpc$^{-1}$). Restricting only to galaxies brighter than $0.626 L^\star_B$ tightens the measurement to $H_0= 77^{+37}_{-18}$ km s$^{-1}$ Mpc$^{-1}$. We show that weighting the host galaxies by stellar mass or star-formation rate provides entirely consistent results with potentially tighter constraints. While these statistical estimates are inferior to the value from the counterpart standard siren measurement utilizing NGC 4993 as the unique host, $H_0=76^{+19}_{-13}$ km s$^{-1}$ Mpc$^{-1}$ (determined from the same publicly available data), our analysis is a proof-of-principle demonstration of the statistical approach first proposed by Bernard Schutz over 30 years ago.


title: Decompositions of locally finite endo-length modules over a skeletally small category

abstract: Given a skeletally small category $\mathcal{C}$, we show that any locally finite endo-length $\mathcal{C}$-module is the direct sum of indecomposable $\mathcal{C}$-modules, whose endomorphism algebra is local.


title: Quantum patterns of genome size variation in angiosperms

abstract: The nuclear DNA amount in angiosperms is studied from the eigen-value equation of the genome evolution operator H. The operator H is introduced by physical simulation and it is defined as a function of the genome size N and the derivative with respective to the size. The discontinuity of DNA size distribution and its synergetic occurrence in related angiosperms species are successfully deduced from the solution of the equation. The results agree well with the existing experimental data of Aloe, Clarkia, Nicotiana, Lathyrus, Allium and other genera. It may indicate that the evolutionary constrains on angiosperm genome are essentially of quantum origin.


title: Wireless Powered Communication Networks: TDD or FDD?

abstract: In this paper, we compare two common modes of duplexing in wireless powered communication networks (WPCN); namely TDD and FDD. So far, TDD has been the most widely used duplexing technique due to its simplicity. Yet, TDD does not allow the energy transmitter to function continuously, which means to deliver the same amount of energy as that in FDD, the transmitter has to have a higher maximum transmit power. On the other hand, when regulations for power spectral density limits are not restrictive, using FDD may lead to higher throughput than that of TDD by allocating less bandwidth to energy and therefore leaving more bandwidth for data. Hence, the best duplexing technique to choose for a specific problem needs careful examination and evaluation.


title: High-precision gravimeter based on a nano-mechanical resonator hybrid with an electron spin

abstract: We show that the gravitational acceleration can be measured with the matter-wave Ramsey interferometry, by using a nitrogen-vacancy center coupled to a nano-mechanical resonator. We propose two experimental methods to realize the Hamiltonian, by using either a cantilever resonator or a trapped nanoparticle. The scheme is robust against the thermal noise, and could be realized at the temperature much higher than the quantum regime one. The effects of decoherence on the interferometry fringe visibility is caculated, including the mechanical motional decay and dephasing of the nitrogen-vacancy center. In addition, we demonstrate that under the various sources of random and systematic noises, our gravimeter can be made on-chip and achieving a high measurement precision. Under experimental feasible parameters, the proposed gravimeter could achieve $10^{-10}$ relative precision.


title: Experimental evaluation of the non-classical relation between measurement errors using entangled photon pairs as a probe

abstract: We have experimentally evaluated the non-classical relation between measurement errors in a joint measurement of two non-commuting polarizations by using entangled photon pairs as a probe. The joint measurement was realized by filtering polarization directions that are sensitive to both of the two target polarizations and the same filtering procedures were applied independently to the two photons of an entangled pair. Since the statistically independent measurement errors of two identical measurements performed on entangled pairs will reduce an overall visibility of the correlation by the square of the local visibilities, the squared visibilities of the local measurements can be obtained directly from the measurement data. We apply this method to determine the visibility of the correlation between the products of the two non-commuting polarizations, a characteristic of the measurement errors that can not be obtained locally since there is no self-adjoint operator that describes this product of measurement outcomes as a single photon observable. The experimental results clearly show that the square of the product visibility is negative, indicating the non-classical nature of the statistical relation between the measurement errors of non-commuting polarizations.


title: One-body Langevin dynamics in heavy-ion collisions at intermediate energies

abstract: We present a new framework to treat the dissipation and fluctuation dynamics associated with nucleon-nucleon scattering in heavy-ion collisions. The two-body collision processes are effectively described in terms of the diffusion of nucleons in viscous nuclear media, governed by a set of Langevin equations in momentum space. The new framework combined with the usual mean field dynamics can be used to simulate heavy-ion collisions at intermediate energies. As a proof of principle, we simulate Au + Au reactions and obtain results consistent with other existing codes under the same constrained conditions. We also study the formation of fragments in Sn + Sn reactions at 50 MeV/nucleon, and results are discussed and compared with two other models commonly employed for collisions.


title: An asynchronous message-passing distributed algorithm for the global critical section problem

abstract: This paper considers the global $(l,k)$-CS problem which is the problem of controlling the system in such a way that, at least $l$ and at most $k$ processes must be in the CS at a time in the network. In this paper, a distributed solution is proposed in the asynchronous message-passing model. Our solution is a versatile composition method of algorithms for $l$-mutual inclusion and $k$-mutual exclusion. Its message complexity is $O(|Q|)$, where $|Q|$ is the maximum size for the quorum of a coterie used by the algorithm, which is typically $|Q| = \sqrt{n}$.


title: A latent factor approach for prediction from multiple assays

abstract: In many domains such as healthcare or finance, data often come in different assays or measurement modalities, with features in each assay having a common theme. Simply concatenating these assays together and performing prediction can be effective but ignores this structure. In this setting, we propose a model which contains latent factors specific to each assay, as well as a common latent factor across assays. We frame our model-fitting procedure, which we call the "Sparse Factor Method" (SFM), as an optimization problem and present an iterative algorithm to solve it.


title: On the Warming of the Southern Hemisphere since 1955 and Recent Slowing-Down: Role of Sea Ice, Sun Spot and El Nino Variability

abstract: The importance of the sea ice retreat in the polar regions for the global warming and the role of ice-albedo feedback was recognized by various authors [1,2]. Similar to a recent study of the phenomenon in the Arctic [3] we present a semi-quantitative estimate of the mechanism for the Southern Hemisphere (SH). Using a simple model, we estimate the contribution of ice-albedo feedback to the mean temperature increase in the SH to be 0.5 +/- 0.1 K in the years 1955 to 2015, while from the simultaneous growth of the greenhouse gases (GHG) we derive a direct warming of only 0.2 +/- 0.05 K in the same period. These numbers are in nice accordance with the reported mean temperature rise of 0.75 +/- 0.1 K of the SH in 2015 since 1955 (and relative to 1880). Our data also confirm previously noticed correlations between the annual fluctuations of solar intensity and El Nino observations on the one hand and the annual variability of the SH surface temperature on the other hand. Our calculations indicate a slowing down of the temperature increase during the past few years that is likely to persist. Assuming a continuation of the present trends for the southern sea ice and GHG concentration we predict the further temperature rise to decrease by 33 % in 2015 to 2025 as compared to the previous decade.


title: Low-coherent optical diffraction tomography by angle-scanning illumination

abstract: We propose and experimentally demonstrate temporally low-coherent optical diffraction tomography (ODT) based on angle-scanning Mach-Zehnder interferometry. Using a digital micromirror device based on diffractive tilting, we successfully maintain full-field interference of incoherent light during every scan sequence. The ODT reconstruction principles for temporally incoherent illuminations are thoroughly reviewed and developed. Several limitations of incoherent illumination are also discussed, such as the nondispersive assumption, optical sectioning capacity, and illumination angle limitation. Using the proposed setup and reconstruction algorithms, we successfully demonstrate low-coherent ODT imaging of microspheres, human red blood cells, and eukaryotic cells.


title: A Simple and Efficient Estimation of the Average Treatment Effect in the Presence of Unmeasured Confounders

abstract: Wang and Tchetgen Tchetgen (2017) studied identification and estimation of the average treatment effect when some confounders are unmeasured. Under their identification condition, they showed that the semiparametric efficient influence function depends on five unknown functionals. They proposed to parameterize all functionals and estimate the average treatment effect from the efficient influence function by replacing the unknown functionals with estimated functionals. They established that their estimator is consistent when certain functionals are correctly specified and attains the semiparametric efficiency bound when all functionals are correctly specified. In applications, it is likely that those functionals could all be misspecified. Consequently their estimator could be inconsistent or consistent but not efficient. This paper presents an alternative estimator that does not require parameterization of any of the functionals. We establish that the proposed estimator is always consistent and always attains the semiparametric efficiency bound. A simple and intuitive estimator of the asymptotic variance is presented, and a small scale simulation study reveals that the proposed estimation outperforms the existing alternatives in finite samples.


title: Reference-free polarization-sensitive quantitative phase imaging using singe-point optical phase conjugation

abstract: We propose and experimentally demonstrate a method of polarization-sensitive quantitative phase imaging using two photo detectors. Instead of recording wide-field interference patterns, finding the modulation patterns maximizing focused intensities in terms of the polarization states enables polarization-dependent quantitative phase imaging without the need for a reference beam and an image sensor. The feasibility of the present method is experimentally validated by reconstructing Jones matrices of various samples including a polystyrene microsphere, a maize starch granule, and a rat retinal nerve fiber layer. Since the present method is simple and sufficiently general, we expect that it may offer solutions for quantitative phase imaging of birefringent materials.


title: High Rate RPC detector for LHC

abstract: The High Luminosity LHC (HL-LHC) phase is designed to increase by an order of magnitude the amount of data to be collected by the LHC experiments. The foreseen gradual increase of the instantaneous luminosity of up to more than twice its nominal value of $10\times10^{34}\ {\rm cm}^{-1}{\rm s}^{-2}$ during Phase I and Phase II of the LHC running, presents special challenges for the experiments. The region with high pseudo rapidity ($\eta$) region of the forward muon spectrometer ($2.4 > |\eta| > 1.9$) is not equipped with RPC stations. The increase of the expected particles rate up to 2 kHz cm$^{-1}$ ( including a safety factor 3 ) motivates the installation of RPC chambers to guarantee redundancy with the CSC chambers already present. The current CMS RPC technology cannot sustain the expected background level. A new generation of Glass-RPC (GRPC) using low-resistivity glass was proposed to equip the two most far away of the four high $\eta$ muon stations of CMS. In their single-gap version they can stand rates of few kHz cm$^{-1}$. Their time precision of about 1 ns can allow to reduce the noise contribution leading to an improvement of the trigger rate. The proposed design for large size chambers is examined and some preliminary results obtained during beam tests at Gamma Irradiation Facility (GIF++) and Super Proton Synchrotron (SPS) at CERN are shown. They were performed to validate the capability of such detectors to support high irradiation environment with limited consequence on their efficiency.