title: Community-Empowered Air Quality Monitoring System

abstract: Developing information technology to democratize scientific knowledge and support citizen empowerment is a challenging task. In our case, a local community suffered from air pollution caused by industrial activity. The residents lacked the technological fluency to gather and curate diverse scientific data to advocate for regulatory change. We collaborated with the community in developing an air quality monitoring system which integrated heterogeneous data over a large spatial and temporal scale. The system afforded strong scientific evidence by using animated smoke images, air quality data, crowdsourced smell reports, and wind data. In our evaluation, we report patterns of sharing smoke images among stakeholders. Our survey study shows that the scientific knowledge provided by the system encourages agonistic discussions with regulators, empowers the community to support policy making, and rebalances the power relationship between stakeholders.


title: A Systematic DNN Weight Pruning Framework using Alternating Direction Method of Multipliers

abstract: Weight pruning methods for deep neural networks (DNNs) have been investigated recently, but prior work in this area is mainly heuristic, iterative pruning, thereby lacking guarantees on the weight reduction ratio and convergence time. To mitigate these limitations, we present a systematic weight pruning framework of DNNs using the alternating direction method of multipliers (ADMM). We first formulate the weight pruning problem of DNNs as a nonconvex optimization problem with combinatorial constraints specifying the sparsity requirements, and then adopt the ADMM framework for systematic weight pruning. By using ADMM, the original nonconvex optimization problem is decomposed into two subproblems that are solved iteratively. One of these subproblems can be solved using stochastic gradient descent, the other can be solved analytically. Besides, our method achieves a fast convergence rate. The weight pruning results are very promising and consistently outperform the prior work. On the LeNet-5 model for the MNIST data set, we achieve 71.2 times weight reduction without accuracy loss. On the AlexNet model for the ImageNet data set, we achieve 21 times weight reduction without accuracy loss. When we focus on the convolutional layer pruning for computation reductions, we can reduce the total computation by five times compared with the prior work (achieving a total of 13.4 times weight reduction in convolutional layers). Our models and codes are released at https://github.com/KaiqiZhang/admm-pruning


title: MmWave MU-MIMO for Aerial Networks

abstract: Millimeter wave offers high bandwidth for air-to-air (A2A) communication. In this paper, we evaluate the rate performance of a multiuser MIMO (MU-MIMO) configuration where several aircraft communicate with a central hub. We consider a hybrid subarray architecture, single path channels, and realistic atmospheric attenuation effects. We propose a mathematical framework for the analysis of millimeter wave (mmWave) MU-MIMO networks. Via Monte Carlo simulation, we demonstrate that mmWave is a promising technology for delivering gigabit connectivity in next-generation aerial networks.


title: Exploring evidence of interaction between dark energy and dark matter

abstract: In this work we use the latest observations on SNIa, $H(z)$, BAO, $f_{gas}$ in clusters and CMBR, to constrain three models showing an explicit interaction between dark matter and dark energy. In particular, we use the BOSS BAO measurements at $z \simeq 0.32$, $0.57$ and $2.34$, using the full 2-dimensional constraints on the angular and line of sight BAO scale. We find that using all five observational probes together, two of the interaction models show positive evidence at more than 3 $ \sigma $. Although significant, further study is needed to establish this statement firmly.


title: A "no-hair" test for binary black holes

abstract: One of the consequences of the black-hole "no-hair" theorem in general relativity (GR) is that gravitational radiation (quasi-normal modes) from a perturbed Kerr black hole is uniquely determined by its mass and spin. Thus, the spectrum of quasi-normal mode frequencies have to be all consistent with the same value of the mass and spin. Similarly, the gravitational radiation from a coalescing binary black hole system is uniquely determined by a small number of parameters (masses and spins of the black holes and orbital parameters). Thus, consistency between different spherical harmonic modes of the radiation is a powerful test that the observed system is a binary black hole predicted by GR. We formulate such a test, develop a Bayesian implementation, demonstrate its performance on simulated data and investigate the possibility of performing such a test using previous and upcoming gravitational wave observations.


title: Mixed H2/H-infinity Data-Driven Control Design for Hard Disk Drives

abstract: A frequency based data-driven control design considering mixed H2/H-infinity control objectives is developed for multiple input-single output systems. The main advantage of the data-driven control over the model-based control is its ability to use the frequency response measurements of the controlled plant directly without the need to identify a model for the plant. In the proposed methodology, multiple sets of measurements can be considered in the design process to accommodate variations in the system dynamics. The controller is obtained by translating the mixed H2/H-infinity control objectives into a convex optimization problem. The H-infinity norm is used to shape closed loop transfer functions and guarantee closed loop stability, while the H2 norm is used to constrain and/or minimize the variance of signals in the time domain. The proposed data-driven design methodology is used to design a track following controller for a dual-stage HDD. The sensitivity decoupling structure[16] is considered as the controller structure. The compensators inside this controller structure are designed and compared by decoupling the system into two single input-single-output systems as well as solving for a single input-double output controller.


title: Efficient illumination angle self-calibration in Fourier ptychography

abstract: Fourier ptychography captures intensity images with varying source patterns (illumination angles) in order to computationally reconstruct large space-bandwidth-product images. Accurate knowledge of the illumination angles is necessary for good image quality; hence, calibration methods are crucial, despite often being impractical or slow. Here, we propose a fast, robust, and accurate self-calibration algorithm that uses only experimentally-collected data and general knowledge of the illumination setup. First, our algorithm makes a direct estimate of the brightfield illumination angles based on image processing. Then, a more computationally-intensive spectral correlation method is used inside the iterative solver to further refine the angle estimates of both brightfield and darkfield images. We demonstrate our method for correcting large and small misalignment artifacts in both 2D and 3D Fourier ptychography with different source types: an LED array, a galvo-steered laser, and a high-NA quasi-dome LED illuminator.


title: Periodic Solutions to Nonlinear Euler-Bernoulli Beam Equations

abstract: Bending vibrations of thin beams and plates may be described by nonlinear Euler-Bernoulli beam equations with $x$-dependent coefficients. In this paper we investigate existence of families of time-periodic solutions to such a model using Lyapunov-Schmidt reduction and a differentiable Nash-Moser iteration scheme. The results hold for all parameters $(\epsilon,\omega)$ in a Cantor set with asymptotically full measure as $\epsilon\rightarrow0$.


title: A Mathematical Framework for Superintelligent Machines

abstract: We describe a class calculus that is expressive enough to describe and improve its own learning process. It can design and debug programs that satisfy given input/output constraints, based on its ontology of previously learned programs. It can improve its own model of the world by checking the actual results of the actions of its robotic activators. For instance, it could check the black box of a car crash to determine if it was probably caused by electric failure, a stuck electronic gate, dark ice, or some other condition that it must add to its ontology in order to meet its sub-goal of preventing such crashes in the future. Class algebra basically defines the eval/eval-1 Galois connection between the residuated Boolean algebras of 1. equivalence classes and super/sub classes of class algebra type expressions, and 2. a residual Boolean algebra of biclique relationships. It distinguishes which formulas are equivalent, entailed, or unrelated, based on a simplification algorithm that may be thought of as producing a unique pair of Karnaugh maps that describe the rough sets of maximal bicliques of relations. Such maps divide the n-dimensional space of up to 2n-1 conjunctions of up to n propositions into clopen (i.e. a closed set of regions and their boundaries) causal sets. This class algebra is generalized to type-2 fuzzy class algebra by using relative frequencies as probabilities. It is also generalized to a class calculus involving assignments that change the states of programs. INDEX TERMS 4-valued Boolean Logic, Artificial Intelligence, causal sets, class algebra, consciousness, intelligent design, IS-A hierarchy, mathematical logic, meta-theory, pointless topological space, residuated lattices, rough sets, type-2 fuzzy sets


title: Purely electronic nanometallic ReRAM

abstract: Resistance switching random access memory (ReRAM), with the ability to repeatedly modulate electrical resistance, has been highlighted as a feasible high-density memory with the potential to replace negative-AND (NAND) flash memory. Such resistance modulation usually involves ion migration and filament formation, which usually lead to relatively low device reliability and yield. Resistance switching can also come from an entirely electronic origin, as in nanometallic memory, by electron trapping and detrapping. Recent research has revealed additional merits of its mechanism, which entails smart, atomic-sized floating gates that can be easily engineered in amorphous Si, oxides, and nitrides. This article addresses the basic ideas of nanometallic ReRAM, which may also be a contender for analogue computing and non-von Neumann-type computation.


title: On top fan-in vs formal degree for depth-$3$ arithmetic circuits

abstract: We show that over the field of complex numbers, \emph{every} homogeneous polynomial of degree $d$ can be approximated (in the border complexity sense) by a depth-$3$ arithmetic circuit of top fan-in at most $d+1$. This is quite surprising since there exist homogeneous polynomials $P$ on $n$ variables of degree $2$, such that any depth-$3$ arithmetic circuit computing $P$ must have top fan-in at least $\Omega(n)$. As an application, we get a new tradeoff between the top fan-in and formal degree in an approximate analog of the celebrated depth reduction result of Gupta, Kamath, Kayal and Saptharishi [GKKS13]. Formally, we show that if a degree $d$ homogeneous polynomial $P$ can be computed by an arithmetic circuit of size $s$, then for every $t \leq d$, $P$ is in the border of a depth-$3$ circuit of top fan-in $s^{O(t)}$ and formal degree $s^{O(d/t)}$. To the best of our knowledge, the upper bound on the top fan-in in the original proof of [GKKS13] is always at least $s^{O(\sqrt{d})}$, regardless of the formal degree.


title: Seeking Open-Ended Evolution in Swarm Chemistry II: Analyzing Long-Term Dynamics via Automated Object Harvesting

abstract: We studied the long-term dynamics of evolutionary Swarm Chemistry by extending the simulation length ten-fold compared to earlier work and by developing and using a new automated object harvesting method. Both macroscopic dynamics and microscopic object features were characterized and tracked using several measures. Results showed that the evolutionary dynamics tended to settle down into a stable state after the initial transient period, and that the extent of environmental perturbations also affected the evolutionary trends substantially. In the meantime, the automated harvesting method successfully produced a huge collection of spontaneously evolved objects, revealing the system's autonomous creativity at an unprecedented scale.


title: Maximally non-Markovian quantum dynamics without environment backflow of information

abstract: The degree of non-Markovianity allows to characterizing quantum evolutions that depart from a Markovian regime in a similar way as Schmidt number measures the degree of entanglement of pure states. Maximally non-Markovian dynamics are the analogous of maximally entangled states [D. Chru\'sci\'nski and S. Maniscalco, Phys. Rev. Lett. 112, 120404 (2014)]. Here, we demonstrate that there exists a class of maximally non-Markovian quantum evolutions where the associated environment (degrees of freedom not belonging to the system) obeys a Markovian (memoryless) dynamics, which in turn is unperturbed by the system state or dynamics. These properties imply the absence of any "physical environment-to-system backflow of information." Non-Markovian features (as usual in quantum systems coupled to dissipative classical degrees of freedom) arise from a unidirectional dependence of the system dynamics on the reservoir states.


title: Demonstration of spatial-light-modulation-based four-wave mixing in cold atoms

abstract: Long-distance quantum optical communications usually require efficient wave-mixing processes to convert the wavelengths of single photons. Many quantum applications based on electromagnetically induced transparency (EIT) have been proposed and demonstrated at the single-photon level, such as quantum memories, all-optical transistors, and cross-phase modulations. However, EIT-based four-wave mixing (FWM) in a resonant double-$\Lambda$ configuration has a maximum conversion efficiency (CE) of 25% because of absorptive loss due to spontaneous emission. An improved scheme using spatially modulated intensities of two control fields has been theoretically proposed to overcome this conversion limit. In this study, we first demonstrate wavelength conversion from 780 to 795 nm with a 43% CE by using this scheme at an optical density (OD) of 19 in cold $^{87}$Rb atoms. According to the theoretical model, the CE in the proposed scheme can further increase to 96% at an OD of 240 under ideal conditions, thereby attaining an identical CE to that of the previous non resonant double-$\Lambda$ scheme at half the OD. This spatial-light-modulation-based FWM scheme can achieve a near-unity CE, thus providing an easy method of implementing an efficient quantum wavelength converter for all-optical quantum information processing.


title: A Web-based Large-scale Timelapse Editor for Creating and Sharing Guided Video Tours and Interactive Slideshows

abstract: Scientists, journalists, and photographers have used advanced camera technology to capture extremely high-resolution timelapse and developed information visualization tools for data exploration and analysis. However, it takes a great deal of effort for professionals to form and tell stories after exploring data, since these tools usually provide little aids in creating visual elements. We present a web-based timelapse editor to support the creation of guided video tours and interactive slideshows from a collection of large-scale spatial and temporal images. Professionals can embed these two visual elements into web pages in conjunction with various forms of digital media to tell multimodal and interactive stories.


title: Adversarial Training Versus Weight Decay

abstract: Performance-critical machine learning models should be robust to input perturbations not seen during training. Adversarial training is a method for improving a model's robustness to some perturbations by including them in the training process, but this tends to exacerbate other vulnerabilities of the model. The adversarial training framework has the effect of translating the data with respect to the cost function, while weight decay has a scaling effect. Although weight decay could be considered a crude regularization technique, it appears superior to adversarial training as it remains stable over a broader range of regimes and reduces all generalization errors. Equipped with these abstractions, we provide key baseline results and methodology for characterizing robustness. The two approaches can be combined to yield one small model that demonstrates good robustness to several white-box attacks associated with different metrics.


title: Thermodynamics and evolutionary biology through optimal control

abstract: We consider a particular instance of the lift of controlled systems recently proposed in the theory of irreversible thermodynamics and show that it leads to a variational principle for an optimal control in the sense of Pontryagin. Then we focus on two important applications: in thermodynamics and in evolutionary biology. In the thermodynamic context, we show that this principle provides a dynamical implementation of the Second Law, which stabilizes the equilibrium manifold of a system. In the evolutionary context, we discuss several interesting features: it provides a robust scheme for the coevolution of the population and its fitness landscape; it has a clear informational interpretation; it recovers Price equation naturally; and finally, it extends standard evolutionary dynamics to include phenomena such as the emergence of cooperation and the coexistence of qualitatively different phases of evolution, which we speculate can be associated with Darwinism and punctuated equilibria.


title: Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units

abstract: We present an implementation of phaseless Auxiliary-Field Quantum Monte Carlo (ph-AFQMC) utilizing graphical processing units (GPUs). The AFQMC method is recast in terms of matrix operations which are spread across thousands of processing cores and are executed in batches using custom Compute Unified Device Architecture kernels and the hardware-optimized cuBLAS matrix library. Algorithmic advances include a batched Sherman-Morrison-Woodbury algorithm to quickly update matrix determinants and inverses, density-fitting of the two-electron integrals, an energy algorithm involving a high-dimensional precomputed tensor, and the use of single-precision floating point arithmetic. These strategies result in dramatic reductions in wall-times for both single- and multi-determinant trial wavefunctions. For typical calculations we find speed-ups of roughly two orders of magnitude using just a single GPU card. Furthermore, we achieve near-unity parallel efficiency using 8 GPU cards on a single node, and can reach moderate system sizes via a local memory-slicing approach. We illustrate the robustness of our implementation on hydrogen chains of increasing length, and through the calculation of all-electron ionization potentials of the first-row transition metal atoms. We compare long imaginary-time calculations utilizing a population control algorithm with our previously published correlated sampling approach, and show that the latter improves not only the efficiency but also the accuracy of the computed ionization potentials. Taken together, the GPU implementation combined with correlated sampling provides a compelling computational method that will broaden the application of ph-AFQMC to the description of realistic correlated electronic systems.


title: A Markov Chain Sampler for Plane Curves

abstract: A plane curve is a knot diagram in which each crossing is replaced by a 4-valent vertex, and so are dual to a subset of planar quadrangulations. The aim of this paper is to introduce a new tool for sampling diagrams via sampling of plane curves. At present the most efficient method for sampling diagrams is rejection sampling, however that method is inefficient at even modest sizes. We introduce Markov chains that sample from the space of plane curves using local moves based on Reidemeister moves. By then mapping vertices on those curves to crossings we produce random knot diagrams. Combining this chain with flat histogram methods we achieve an efficient sampler of plane curves and knot diagrams. By analysing data from this chain we are able to estimate the number of knot diagrams of a given size and also compute knotting probabilities and so investigate their asymptotic behaviour.


title: Crafting a Toolchain for Image Restoration by Deep Reinforcement Learning

abstract: We investigate a novel approach for image restoration by reinforcement learning. Unlike existing studies that mostly train a single large network for a specialized task, we prepare a toolbox consisting of small-scale convolutional networks of different complexities and specialized in different tasks. Our method, RL-Restore, then learns a policy to select appropriate tools from the toolbox to progressively restore the quality of a corrupted image. We formulate a step-wise reward function proportional to how well the image is restored at each step to learn the action policy. We also devise a joint learning scheme to train the agent and tools for better performance in handling uncertainty. In comparison to conventional human-designed networks, RL-Restore is capable of restoring images corrupted with complex and unknown distortions in a more parameter-efficient manner using the dynamically formed toolchain.


title: Cortex Neural Network: learning with Neural Network groups

abstract: Neural Network has been successfully applied to many real-world problems, such as image recognition and machine translation. However, for the current architecture of neural networks, it is hard to perform complex cognitive tasks, for example, to process the image and audio inputs together. Cortex, as an important architecture in the brain, is important for animals to perform the complex cognitive task. We view the architecture of Cortex in the brain as a missing part in the design of the current artificial neural network. In this paper, we purpose Cortex Neural Network (CrtxNN). The Cortex Neural Network is an upper architecture of neural networks which motivated from cerebral cortex in the brain to handle different tasks in the same learning system. It is able to identify different tasks and solve them with different methods. In our implementation, the Cortex Neural Network is able to process different cognitive tasks and perform reflection to get a higher accuracy. We provide a series of experiments to examine the capability of the cortex architecture on traditional neural networks. Our experiments proved its ability on the Cortex Neural Network can reach accuracy by 98.32% on MNIST and 62% on CIFAR10 at the same time, which can promisingly reduce the loss by 40%.


title: Optical solitons in media with focusing and defocusing saturable nonlinearity and a parity-time-symmetric external potential

abstract: We report results for solitons in models of waveguides with focusing or defocusing saturable nonlinearity and a parity-time(PT )-symmetric complex-valued external potential of the Scarf-II type. The model applies to the nonlinear wave propagation in graded-index optical waveguides with balanced gain and loss. We find both fundamental and multipole solitons for both focusing and defocusing signs of the saturable nonlinearity in such PT -symmetric waveguides. The dependence of the propagation constant on the soliton's power is presented for different strengths of the nonlinearity saturation, S. The stability of fundamental, dipole, tripole, and quadrupole solitons is investigated by means of the linear-stability analysis and direct numerical simulations of the corresponding (1+1)-dimensional nonlinear Schrodinger-type equation. The results show that the instability of the stationary solutions can be mitigated or completely suppressed, increasing the value of S.


title: Individual and Group Stability in Neutral Restrictions of Hedonic Games

abstract: We consider a class of coalition formation games called hedonic games, i.e., games in which the utility of a player is completely determined by the coalition that the player belongs to. We first define the class of subset-additive hedonic games and show that they have the same representation power as the class of hedonic games. We then define a restriction of subset-additive hedonic games that we call subset-neutral hedonic games and generalize a result by Bogomolnaia and Jackson (2002) by showing the existence of a Nash stable partition and an individually stable partition in such games. We also consider neutrally anonymous hedonic games and show that they form a subclass of the subset-additive hedonic games. Finally, we show the existence of a core stable partition that is also individually stable in neutrally anonymous hedonic games by exhibiting an algorithm to compute such a partition.


title: Optimal Solution of Vehicle Routing Problems with Fractional Objective Function

abstract: This work proposes a first extensive analysis of the Vehicle Routing Problem with Fractional Objective Function (vrpfof). We investigate how the principal techniques used either in the context of fractional programming or in the context of vehicle routing problems interact. We present new dual and primal bounding procedures which have been incorporated in an exact method. The method proposed allows to extend specific variants of \vrp\ to their counterpart with a fractional objective function. Extensive numerical experiments prove the validity of our approach.


title: Question Answering over Freebase via Attentive RNN with Similarity Matrix based CNN

abstract: With the rapid growth of knowledge bases (KBs), question answering over knowledge base, a.k.a. KBQA has drawn huge attention in recent years. Most of the existing KBQA methods follow so called encoder-compare framework. They map the question and the KB facts to a common embedding space, in which the similarity between the question vector and the fact vectors can be conveniently computed. This, however, inevitably loses original words interaction information. To preserve more original information, we propose an attentive recurrent neural network with similarity matrix based convolutional neural network (AR-SMCNN) model, which is able to capture comprehensive hierarchical information utilizing the advantages of both RNN and CNN. We use RNN to capture semantic-level correlation by its sequential modeling nature, and use an attention mechanism to keep track of the entities and relations simultaneously. Meanwhile, we use a similarity matrix based CNN with two-directions pooling to extract literal-level words interaction matching utilizing CNNs strength of modeling spatial correlation among data. Moreover, we have developed a new heuristic extension method for entity detection, which significantly decreases the effect of noise. Our method has outperformed the state-of-the-arts on SimpleQuestion benchmark in both accuracy and efficiency.


title: Design and Implementation of a Thomson Parabola for Fluence Dependent Energy-Loss Measurements at the Neutralized Drift Compression eXperiment

abstract: The interaction of ion beams with matter includes the investigation of the basic principles of ion stopping in heated materials. An unsolved question is the effect of different, especially higher, ion beam fluences on ion stopping in solid targets. This is relevant in applications such as in fusion sciences. To address this question, a Thomson parabola was built for the Neutralized Drift Compression eXperiment (NDCX-II) for ion energy-loss measurements at different ion beam fluences. The linear induction accelerator NDCX-II delivers 2 ns short, intense ion pulses, up to several tens of nC/pulse, or 10$^{10}$-10$^{11}$ ions, with a peak kinetic energy of ~1.1 MeV and a minimal spot size of 2 mm FWHM. For this particular accelerator the energy determination with conventional beam diagnostics, for example, time of flight measurements, is imprecise due to the non-trivial longitudinal phase space of the beam. In contrast, a Thomson parabola is well suited to reliably determine the beam energy distribution. The Thomson parabola differentiates charged particles by energy and charge-to-mass ratio, through deflection of charged particles by electric and magnetic fields. During first proof-of-principle experiments, we achieved to reproduce the average initial helium beam energy as predicted by computer simulations with a deviation of only 1.4 %. Successful energy-loss measurements with 1 {\mu}m thick Silicon Nitride foils show the suitability of the accelerator for such experiments. The initial ion energy was determined during a primary measurement without a target, while a second measurement, incorporating the target, was used to determine the transmitted energy. The energy-loss was then determined as the difference between the two energies.


title: Uniqueness and convergence on equilibria of the Keller-Segel system with subcritical mass

abstract: This paper is concerned with the uniqueness of solutions to the following nonlocal semi-linear elliptic equation \begin{equation}\label{ellip}\tag{$\ast$} \Delta u-\beta u+\lambda\frac{e^u}{\int_{\Omega}e^u}=0~\mathrm{in}~\Omega, \end{equation} where $\Omega$ is a bounded domain in $\mathbb{R}^2$ and $\beta, \lambda$ are positive parameters. The above equation arises as the stationary problem of the well-known classical Keller-Segel model describing chemotaxis. For equation \eqref{ellip} with Neumann boundary condition, we establish an integral inequality and prove that the solution of (\ref{ellip}) is unique if $0<\lambda \leq 8\pi$ and $u$ satisfies some symmetric properties. While for \eqref{ellip} with Dirichlet boundary condition, the same uniqueness result is obtained without symmetric condition by a different approach inspired by some recent works [19,21]. As an application of the uniqueness results, we prove that the radially symmetric solution of the classical Keller-Segel system with subcritical mass subject to Neumann boundary conditions will converge to the unique constant equilibrium as time tends to infinity if $\Omega$ is a disc in two dimensions. As far as we know, this is the first result that asserts the exact asymptotic behavior of solutions to the classical Keller-Segel system with subcritical mass in two dimensions.


title: Automated Detection Methods for Solar Activities and an Application for Statistic Analysis of Solar Filament

abstract: With the rapid development of telescopes, both temporal cadence and the spatial resolution of observations are increasing. This in turn generates vast amount of data, which can be efficiently searched only with automated detections in order to derive the features of interest in the observations. A number of automated detection methods and algorithms have been developed for solar activities, based on the image processing and machine learning techniques. In this paper, after briefly reviewing some automated detection methods, we describe our efficient and versatile automated detection method for solar filaments. It is able not only to recognize filaments, determine the features such as the position, area, spine, and other relevant parameters, but also to trace the daily evolution of the filaments. It is applied to process the full disk H-alpha data observed in nearly three solar cycles, and some statistic results are presented.


title: Lab::Measurement - A portable and extensible framework for controlling lab equipment and conducting measurements

abstract: Lab::Measurement is a framework for test and measurement automatization using Perl 5. While primarily developed with applications in mesoscopic physics in mind, it is widely adaptable. Internally, a layer model is implemented. Communication protocols such as IEEE 488, USB Test & Measurement, or, e.g., VXI-11 are addressed by the connection layer. The wide range of supported connection backends enables unique cross-platform portability. At the instrument layer, objects correspond to equipment connected to the measurement PC (e.g., voltage sources, magnet power supplies, multimeters, etc.). The high-level sweep layer automates the creation of measurement loops, with simultaneous plotting and data logging. An extensive unit testing framework is used to verify functionality even without connected equipment. Lab::Measurement is distributed as free and open source software.


title: Abelian networks IV. Dynamics of nonhalting networks

abstract: An abelian network is a collection of communicating automata whose state transitions and message passing each satisfy a local commutativity condition. This paper is a continuation of the abelian networks series of Bond and Levine (2016), for which we extend the theory of abelian networks that halt on all inputs to networks that can run forever. A nonhalting abelian network can be realized as a discrete dynamical system in many different ways, depending on the update order. We show that certain features of the dynamics, such as minimal period length, have intrinsic definitions that do not require specifying an update order. We give an intrinsic definition of the \emph{torsion group} of a finite irreducible (halting or nonhalting) abelian network, and show that it coincides with the critical group of Bond and Levine (2016) if the network is halting. We show that the torsion group acts freely on the set of invertible recurrent components of the trajectory digraph, and identify when this action is transitive. This perspective leads to new results even in the classical case of sinkless rotor networks (deterministic analogues of random walks). In Holroyd et. al (2008) it was shown that the recurrent configurations of a sinkless rotor network with just one chip are precisely the unicycles (spanning subgraphs with a unique oriented cycle, with the chip on the cycle). We generalize this result to abelian mobile agent networks with any number of chips. We give formulas for generating series such as \[ \sum_{n \geq 1} r_n z^n = \det (\frac{1}{1-z}D - A ) \] where $r_n$ is the number of recurrent chip-and-rotor configurations with $n$ chips; $D$ is the diagonal matrix of outdegrees, and $A$ is the adjacency matrix. A consequence is that the sequence $(r_n)_{n \geq 1}$ completely determines the spectrum of the simple random walk on the network.


title: Comparative study of discrete Boltzmann model and Navier-Stokes

abstract: Discrete Boltzmann model (DBM) is a type of coarse-grained mesoscale kinetic model derived from the Boltzmann equation. Physically, it is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model for the relevant thermodynamic non-equilibrium (TNE) behaviours. The Navier-Stokes (NS) model is a traditional macroscopic hydrodynamic model based on continuity hypothesis and conservation laws. In this study, the two models are compared from two aspects, physical capability and computational cost, by simulating two kinds of flow problems including the thermal Couette flow and a Mach 3 step problem. In the cases where the TNE effects are weak, both the two models give accurate results for the hydrodynamic behaviour. Besides, DBM can provide more detailed non-equilibrium information, while the NS is more efficient if concern only the density, momentum, energy and their derived quantities. It is concluded that, if the TNE effects are strong or are to be investigated, the NS is insufficient while DBM is a good choice. While in the cases where the TNE effects are weak and only the macro flow fields are to be studied, the NS is more preferable.


title: Multiple solutions for a fractional Schrodinger equation with potentials

abstract: This paper is devoted to study a class of nonlinear fractional Schr\"{o}dinger equations: \begin{equation*} (-\Delta)^{s}u+V(x)u=f(x,u), \quad \text{in}\: \mathbb{R}^{N}, \end{equation*} where $s\in (0,1)$, $\ N>2s$, $(-\Delta)^{s}$ stands for the fractional Laplacian. First, by using a variational approach, we establish the existence of at least one nontrivial solution for the above equation with a general potential $V(x)$ which is allowed to be sign-changing and a sublinear nonlinearity $f(x,u)$. Next, by using variational methods and the Moser iteration technique, we prove the existence of infinitely many solutions with $V(x)$ is a nonnegative potential and the nonlinearity $f(x,u)$ is locally sublinear with respect to $u$.


title: GeV emission in the region of HESS J1809-193 and HESS J1813-178: is HESS J1809-193 a proton pevatron?

abstract: HESS J1809-193 is an unidentified TeV source discovered by the High Energy Stereoscopic System and originally classified as a pulsar wind nebula (PWN) candidate associated with the pulsar PSR J1809-1917. However, a recent study of deep radio observations and the interstellar medium near the source has found evidence for a hadronic scenario for the gamma-rays. Here, a detailed study of the GeV emission in the region using data from the \emph{Fermi} LAT is presented. The GeV emission has an extended morphology in the region of the TeV emission and the overall spectrum can be accounted for by a cosmic ray population having a simple power-law spectrum with energies extending up to 1 PeV. However, the spectrum at tens of TeV should be observed more deeply in the future to confirm its hadronic nature, and other scenarios involving combinations of leptonic and hadronic emission from several of the known supernova remnants in the region cannot be ruled out. The nearby TeV source HESS J1813-178, thought to be a PWN, is also studied in detail at GeV energies and we find a region of significant emission which is much more extended than the TeV emission and whose spectrum is softer than expected from a PWN but similar to those seen in several star forming regions that are believed to accelerate protons. There is marginal evidence for a GeV point source at the location of the X-ray PWN, beside the extended emission.


title: Increasing Parallelism in the ROOT I/O Subsystem

abstract: When processing large amounts of data, the rate at which reading and writing can take place is a critical factor. High energy physics data processing relying on ROOT is no exception. The recent parallelisation of LHC experiments' software frameworks and the analysis of the ever increasing amount of collision data collected by experiments further emphasized this issue underlying the need of increasing the implicit parallelism expressed within the ROOT I/O. In this contribution we highlight the improvements of the ROOT I/O subsystem which targeted a satisfactory scaling behaviour in a multithreaded context. The effect of parallelism on the individual steps which are chained by ROOT to read and write data, namely (de)compression, (de)serialisation, access to storage backend, are discussed. Performance measurements are discussed through real life examples coming from CMS production workflows on traditional server platforms and highly parallel architectures such as Intel Xeon Phi.


title: Implementing Push-Pull Efficiently in GraphBLAS

abstract: We factor Beamer's push-pull, also known as direction-optimized breadth-first-search (DOBFS) into 3 separable optimizations, and analyze them for generalizability, asymptotic speedup, and contribution to overall speedup. We demonstrate that masking is critical for high performance and can be generalized to all graph algorithms where the sparsity pattern of the output is known a priori. We show that these graph algorithm optimizations, which together constitute DOBFS, can be neatly and separably described using linear algebra and can be expressed in the GraphBLAS linear-algebra-based framework. We provide experimental evidence that with these optimizations, a DOBFS expressed in a linear-algebra-based graph framework attains competitive performance with state-of-the-art graph frameworks on the GPU and on a multi-threaded CPU, achieving 101 GTEPS on a Scale 22 RMAT graph.


title: On the abundance of SRB measures

abstract: We prove the abundance of Sinai-Ruelle-Bowen measures for diffeomorphisms away from ones with a homoclinic tangency. This is motivated by conjectures of Palis on the existence of physical (Sinai-Ruelle-Bowen) measures for global dynamics. The main novelty in this paper is that we have to deeply study Gibbs $cu$-states in different levels. Note that we have to use random perturbations to give some upper bound of the level of Gibbs $cu$-states.


title: Representation Tradeoffs for Hyperbolic Embeddings

abstract: Hyperbolic embeddings offer excellent quality with few dimensions when embedding hierarchical data structures like synonym or type hierarchies. Given a tree, we give a combinatorial construction that embeds the tree in hyperbolic space with arbitrarily low distortion without using optimization. On WordNet, our combinatorial embedding obtains a mean-average-precision of 0.989 with only two dimensions, while Nickel et al.'s recent construction obtains 0.87 using 200 dimensions. We provide upper and lower bounds that allow us to characterize the precision-dimensionality tradeoff inherent in any hyperbolic embedding. To embed general metric spaces, we propose a hyperbolic generalization of multidimensional scaling (h-MDS). We show how to perform exact recovery of hyperbolic points from distances, provide a perturbation analysis, and give a recovery result that allows us to reduce dimensionality. The h-MDS approach offers consistently low distortion even with few dimensions across several datasets. Finally, we extract lessons from the algorithms and theory above to design a PyTorch-based implementation that can handle incomplete information and is scalable.


title: Systematic search for tidal features around nearby galaxies: I. Enhanced SDSS imaging of the Local Volume

abstract: In hierarchical models of galaxy formation, stellar tidal streams are expected around most, if not all, galaxies. Although these features may provide useful diagnostics of the $\Lambda$CDM model, their observational properties remain poorly constrained because they are challenging to detect and interpret and have been studied in detail for only a sparse sampling of galaxy population. More quantitative, systematic approaches are required. We advocate statistical analysis of the counts and properties of such features in archival wide-field imaging surveys for a direct comparison against results from numerical simulations. Thus, in this paper we aim to study systematically the frequency of occurrence and other observational properties of tidal features around nearby galaxies. The sample we construct will act as a foundational dataset for statistical comparison with cosmological models of galaxy formation. Our approach is based on a visual classification of diffuse features around a volume-limited sample of nearby galaxies, using a post-processing of Sloan Digital Sky Survey imaging optimized for the detection of stellar structure with low surface brightness. At a limiting surface brightness of $28\ \mathrm{mag~arcsec^{-2}}$, 14% of the galaxies in our sample exhibit evidence of diffuse features likely to have arisen from minor merging events. Our technique recovers all previously known streams in our sample and yields a number of new candidates. Consistent with previous studies, coherent arc-like features and shells are the most common type of tidal structures found in this study. We conclude that although some detections are ambiguous and could be corroborated or refuted with deeper imaging, our technique provides a reliable foundation for the statistical analysis of diffuse circumgalactic features in wide-area imaging surveys, and for the identification of targets for follow-up studies.


title: Gravitational major-axis contraction of Mercury's elliptical orbit

abstract: The local curvature of the space produced by the Sun causes not only the perihelion precession of Mercury's elliptical orbit, but also the variations of the whole orbit, in comparison with those predicted by the Newtonian theory of gravitation. Calculations show that the gravitational major-axis contraction of the Mercury's elliptical orbit is 1.3 kilometers which can be confirmed by the present astronomical distance measurement technology.


title: One-Dimensional Sums and Finitized Characters of $2 \times 2$ Fused RSOS Models

abstract: Tartaglia and Pearce have argued that the nonunitary $n\times n$ fused Forrester-Baxter $\mbox{RSOS}(m,m')$ models are described, in the continuum scaling limit, by the minimal models ${\cal M}(M,M',n)$ constructed as the higher-level conformal cosets $(A^{(1)}_1)_k\otimes (A^{(1)}_1)_n/(A^{(1)}_1)_{k+n}$ at integer fusion level $n\ge 1$ and fractional level $k=nM/(M'\!-\!M)-2$ with $(M,M')=\big(nm-(n\!-\!1)m',m'\big)$. These results rely on Yang-Baxter integrability and are valid in Regime III for models determined by the crossing parameter $\lambda=(m'\!-\!m)\pi/m'$ in the interval $0<\lambda<\pi/n$. Here we consider the $2\times 2$ $\mbox{RSOS}(m,m')$ models in the interval $\tfrac{\pi}{2}<\lambda<\pi$ and investigate the associated one-dimensional sums. In this interval, we verify that the one-dimensional sums produce new finitized Virasoro characters $ch_{r,s}^{(N)}(q)$ of the minimal models ${\cal M}(m,m',1)$ with $m'>2m$. We further conjecture finitized bosonic forms and check that these agree with the ground state one-dimensional sums out to system sizes $N=12$. The $2\times 2$ $\mbox{RSOS}(m,m')$ models thus realize new Yang-Baxter integrable models in the universality classes of the minimal models ${\cal M}(m,m',1)$. For the series ${\cal M}(m,2m+1,1)$ with $m\ge 2$, the spin-1 one-dimensional sums were previously analysed by Jacob and Mathieu without the underlying Yang-Baxter structure. Finitized Kac characters $\chi_{r,s}^{m,m';(N)}(q)$ for the logarithmic minimal models ${\cal LM}(p,p',1)$ are also obtained for $p'\ge 2p$ by taking the logarithmic limit $m,m'\to\infty$ with $m'/m\to p'/p+$.


title: Spinon Fractionalization from Dynamic Structure Factor of Spin-1/2 Heisenberg Antiferromagnet on the Kagome Lattice

abstract: We study dynamical spin structure factor (DSSF) of S = 1/2 Heisenberg model on the kagome lattice (KAFM) by means of density-matrix renormalization group. By comparison with the well-defined magnetic ordered state and chiral spin liquid sitting nearby in the phase diagram, the KAFM with the nearest-neighbor interaction shows distinct dynamical response behaviors. First of all, the DSSF displays important spectral intensity predominantly at low frequency region around Q = M point in momentum space, and shows a broad spectral distribution at high frequency region for momenta along the boundary of the extended Brillouin zone. Secondly, spinon excitation spectrum is identified from momentum and energy resolved DSSF, which shows critical behavior with much reduced spectrum intensity comparing to the neighboring chiral spin liquid. By adding a weak Dzyaloshinkii-Moriya interaction, the DSSF demonstrates a strong sensitivity to the boundary condition smore consistent with a gapless spin liquid. These results capture the main observations in the inelastic neutron scattering measurements of herbertsmithite, and indicate the spin liquid nature of the ground state with fractionalized spinon excitations. By following the DSSF crossing the quantum phase transition between the CSL and the magnetical ordered phase, we identify the spinon condensation driving the quantum phase transition.


title: TIDBD: Adapting Temporal-difference Step-sizes Through Stochastic Meta-descent

abstract: In this paper, we introduce a method for adapting the step-sizes of temporal difference (TD) learning. The performance of TD methods often depends on well chosen step-sizes, yet few algorithms have been developed for setting the step-size automatically for TD learning. An important limitation of current methods is that they adapt a single step-size shared by all the weights of the learning system. A vector step-size enables greater optimization by specifying parameters on a per-feature basis. Furthermore, adapting parameters at different rates has the added benefit of being a simple form of representation learning. We generalize Incremental Delta Bar Delta (IDBD)---a vectorized adaptive step-size method for supervised learning---to TD learning, which we name TIDBD. We demonstrate that TIDBD is able to find appropriate step-sizes in both stationary and non-stationary prediction tasks, outperforming ordinary TD methods and TD methods with scalar step-size adaptation; we demonstrate that it can differentiate between features which are relevant and irrelevant for a given task, performing representation learning; and we show on a real-world robot prediction task that TIDBD is able to outperform ordinary TD methods and TD methods augmented with AlphaBound and RMSprop.


title: An Investigation of intra-cluster light evolution using cosmological hydro-dynamical simulations

abstract: The intra-cluster light (ICL) in observations is usually identified through the surface brightness limit method. In this paper, for the first time we produce the mock images of galaxy groups and clusters using a cosmological hydro- dynamical simulation, to investigate the ICL fraction and focus on its dependence on observational parameters, e.g., the surface brightness limit (SBL), the effects of cosmological redshift dimming, point spread function and CCD pixel size. Detailed analyses suggest that the width of point spread function has a significant effect on the measured ICL fraction, while the relatively small pixel size shows almost no influence. It is found that the measured ICL fraction depends strongly on the SBL. At a fixed SBL and redshift, the measured ICL fraction decreases with increasing halo mass, while with a much faint SBL, it does not depend on halo mass at low redshifts. In our work, the measured ICL fraction shows clear dependence on the cosmological redshift dimming effect. It is found that there are more mass locked in ICL component than light, suggesting that the use of a constant mass-to-light ratio at high surface brightness levels will lead to an underestimate of ICL mass. Furthermore, it is found that the radial profile of ICL shows a characteristic radius which is almost independent of halo mass. The current measurement of ICL from observations has a large dispersion due to different methods, and we emphasize the importance of using the same definition when observational results are compared with the theoretical predictions.


title: Revealing the Micro-Structure of the Giant Component in Random Graph Ensembles

abstract: The micro-structure of the giant component of the Erd{\H o}s-R\'enyi network and other configuration model networks is analyzed using generating function methods. While configuration model networks are uncorrelated, the giant component exhibits a degree distribution which is different from the overall degree distribution of the network and includes degree-degree correlations of all orders. We present exact analytical results for the degree distributions as well as higher order degree-degree correlations on the giant components of configuration model networks. We show that the degree-degree correlations are essential for the integrity of the giant component, in the sense that the degree distribution alone cannot guarantee that it will consist of a single connected component. To demonstrate the importance and broad applicability of these results, we apply them to the study of the distribution of shortest path lengths on the giant component, percolation on the giant component and the spectra of sparse matrices defined on the giant component. We show that by using the degree distribution on the giant component, one obtains high quality results for these properties, which can be further improved by taking the degree-degree correlations into account. This suggests that many existing methods, currently used for the analysis of the whole network, can be adapted in a straightforward fashion to yield results conditioned on the giant component.


title: On the origin of asymmetric fission of actinides

abstract: Nuclear fission of heavy (actinide) nuclei results predominantly in asymmetric mass-splits. Without quantum shells, which can give extra binding energy to these mass-asymmetric shapes, the nuclei would fission symmetrically. The strongest shell effects are in spherical nuclei, so naturally the spherical "doubly-magic" ${^{132}}$Sn nucleus (${Z=50}$ protons), was expected to play a major role. However, a systematic study of fission has shown that the heavy fragments are distributed around ${Z=52}$ to 56, indicating that ${^{132}}$Sn is not the only driver. Reconciling the strong spherical shell effects at ${Z=50}$ with the different ${Z}$ values of fission fragments observed in nature has been a longstanding puzzle. Here, we show that the final mass asymmetry of the fragments is also determined by the extra stability of octupole (pear-shaped) deformations which have been recently confirmed experimentally around $^{144}$Ba (${Z=56}$), one of very few nuclei with shell-stabilized octupole deformation. Using a modern quantum many-body model of superfluid fission dynamics, we found that heavy fission fragments are produced predominantly with ${52-56}$ protons, associated with significant octupole deformation acquired on the way to fission. These octupole shapes favouring asymmetric fission are induced by deformed shells at ${Z=52}$ and 56. In contrast, spherical "magic" nuclei are very resistant to octupole deformation, which hinders their production as fission fragments. These findings may explain surprising observations of asymmetric fission of lighter than lead nuclei.


title: Development of Integration-Type Silicon-On-Insulator Monolithic Pixel Detectors by Using a Float Zone Silicon

abstract: In this paper, we describe the development of monolithic pixel detectors by using a Silicon-on-Insulator (SOI) technology for X-ray and charged particle applications. The detectors are based on a 0.2 {\mu}m CMOS fully depleted SOI process (Lapis Semiconductor Co., Ltd). The SOI wafer consists of a thick high-resistivity substrate for sensing and a thin low resistivity Si layer for CMOS circuits. We developed the integration-type SOI pixel detector, INTPIX4 mainly for X-ray imaging; it is made of a Float Zone (FZ) or Czochralski (Cz) silicon wafer. Since 2005, Cz SOI detectors were used initially. After 2011, FZ SOI detectors were successfully fabricated. In this paper, we state recent progresses and test results of the SOI monolithic pixel detector using a FZ silicon and compare them with the results obtained using the Cz detector.


title: Keldysh-Rutherford model for attoclock

abstract: We demonstrate a clear similarity between attoclock offset angles and Rutherford scattering angles taking the Keldysh tunnelling width as the impact parameter and the vector potential of the driving pulse as the asymptotic velocity. This simple model is tested against the solution of the time-dependent Schr\"odinger equation using hydrogenic and screened (Yukawa) potentials of equal binding energy. We observe a smooth transition from a hydrogenic to 'hard-zero' intensity dependence of the offset angle with variation of the Yukawa screening parameter. Additionally we make comparison with the attoclock offset angles for various noble gases obtained with the classical-trajectory Monte Carlo method. In all cases we find a close correspondence between the model predictions and numerical calculations. This suggests a largely Coulombic origin of the attoclock offset angle and casts further doubt on its interpretation in terms of a finite tunnelling time.


title: A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). I. The Optical Counterparts of FIRST Radio Sources

abstract: We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50% of the FIRST sources in the search field, which substantially exceed previously reported fractions of SDSS counterparts (i ~< 22) of ~30%. Among the matched sample, 9% are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (AGNs) beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.


title: Regimes of wrinkling in an indented floating elastic sheet

abstract: A thin, elastic sheet floating on the surface of a liquid bath wrinkles when poked at its centre. We study the onset of wrinkling as well as the evolution of the pattern as indentation progresses far beyond the wrinkling threshold. We use tension field theory to describe the macroscopic properties of the deformed film and show that the system passes through a host of different regimes, even while the deflections and strains remain small. We show that the effect of the finite size of the sheet ultimately plays a key role in determining the location of the wrinkle pattern, and obtain scaling relations that characterize the number of wrinkles at threshold and its variation as the indentation progresses. Some of our predictions are confirmed by recent experiments of Ripp \emph{et al.} [arxiv: 1804.02421].


title: Toward Formalizing Teleportation of Pedagogical Artificial Agents

abstract: Our paradigm for the use of artificial agents to teach requires among other things that they persist through time in their interaction with human students, in such a way that they "teleport" or "migrate" from an embodiment at one time t to a different embodiment at later time t'. In this short paper, we report on initial steps toward the formalization of such teleportation, in order to enable an overseeing AI system to establish, mechanically, and verifiably, that the human students in question will likely believe that the very same artificial agent has persisted across such times despite the different embodiments.


title: Modular Generative Adversarial Networks

abstract: Existing methods for multi-domain image-to-image translation (or generation) attempt to directly map an input image (or a random vector) to an image in one of the output domains. However, most existing methods have limited scalability and robustness, since they require building independent models for each pair of domains in question. This leads to two significant shortcomings: (1) the need to train exponential number of pairwise models, and (2) the inability to leverage data from other domains when training a particular pairwise mapping. Inspired by recent work on module networks, this paper proposes ModularGAN for multi-domain image generation and image-to-image translation. ModularGAN consists of several reusable and composable modules that carry on different functions (e.g., encoding, decoding, transformations). These modules can be trained simultaneously, leveraging data from all domains, and then combined to construct specific GAN networks at test time, according to the specific image translation task. This leads to ModularGAN's superior flexibility of generating (or translating to) an image in any desired domain. Experimental results demonstrate that our model not only presents compelling perceptual results but also outperforms state-of-the-art methods on multi-domain facial attribute transfer.


title: Quantizations of the classical time of arrival and their dynamics

abstract: The classical time of arrival in the interacting case is quantized by way of quantizing its expansion about the free time of arrival. The quantization is formulated in coordinate representation which represents ordering rules in terms of two variable polynomial functions. This leads to representations of the quantized time of arrival operators as integral operators whose kernels are determined by the chosen ordering rule. The formulation lends itself to generalization which allows construction of time of arrival operators that cannot be obtained by direct quantization using particular ordering rules. Wey, symmetric and Born-Jordan quantizations are specifically studied. The dynamics of the eigenfunctions of the different time of arrival operators are investigated. The eigenfunctions exhibit unitary arrival at the intended arrival point at their respective eigenvalues.


title: Force-detected high-frequency electron spin resonance spectroscopy using magnet-mounted nanomembrane: robust detection of thermal magnetization modulation

abstract: In this study, we report a conceptually novel broadband high-frequency electron spin resonance (HFESR) spectroscopic technique. In contrast to the ordinary force-detected ESR technique, which detects the magnetization change due to the saturation effect, this method measures the magnetization change due to the change of the sample temperature at resonance. To demonstrate its principle, we developed a silicon nitride nanomembrane-based force-detected ESR spectrometer, which can be stably operated even at high magnetic fields. Test measurements were performed for samples with different spin relaxation times. We succeeded in obtaining a seamless ESR spectrum in magnetic fields of 15~T and frequencies of 636~GHz without significant spectral distortion. A high spin sensitivity of $10^{12}$~spins/G$\cdot$s was obtained, which was independent of the spin relaxation time. These results show that this technique can be used as a practical method in research fields where the HFESR technique is applicable.


title: Learning Latent Events from Network Message Logs

abstract: We consider the problem of separating error messages generated in large distributed data center networks into error events. In such networks, each error event leads to a stream of messages generated by hardware and software components affected by the event. These messages are stored in a giant message log. We consider the unsupervised learning problem of identifying the signatures of events that generated these messages; here, the signature of an error event refers to the mixture of messages generated by the event. One of the main contributions of the paper is a novel mapping of our problem which transforms it into a problem of topic discovery in documents. Events in our problem correspond to topics and messages in our problem correspond to words in the topic discovery problem. However, there is no direct analog of documents. Therefore, we use a non-parametric change-point detection algorithm, which has linear computational complexity in the number of messages, to divide the message log into smaller subsets called episodes, which serve as the equivalents of documents. After this mapping has been done, we use a well-known algorithm for topic discovery, called LDA, to solve our problem. We theoretically analyze the change-point detection algorithm, and show that it is consistent and has low sample complexity. We also demonstrate the scalability of our algorithm on a real data set consisting of $97$ million messages collected over a period of $15$ days, from a distributed data center network which supports the operations of a large wireless service provider.


title: Low-Energy Gamow-Teller Transitions in deformed $N = Z$ odd-odd Nuclei

abstract: We have investigated the Gamow-Teller transitions from ${}^{22}\textrm{Ne}$ to ${}^{22}\textrm{Na}$ by applying the isospin projected antisymmetrized molecular dynamics combined with generator coordinate method. We have found that the GT strength from $^AZ(J^\pi T)={}^{22}\textrm{Ne}(0_1^+1)$ is fragmented into two final states ${}^{22}\textrm{Na}(1_{1,2}^+0)$, which belong to $K=0$ and $K=1$ bands constructed by a prolately deformed $^{20}$Ne core with a $S=1$ proton-neutron ($pn$) pair. Coupling of the intrinsic-spin of the $pn$ pair with the core deformation play an important role in the GT fragmentation. The symmetry breaking in the intrinsic-spin rotation leads to the $SU(4)$-symmetry breaking of the $NN$ pair and causes the GT fragmentation. We have compared the features of the GT transitions with those for ${}^{10}\textrm{Be}\to{}^{10}\textrm{B}$ and discuss the link of the $SU(4)$-symmetry and the GT fragmentation in the deformed systems.


title: Mean Field Network based Graph Refinement with application to Airway Tree Extraction

abstract: We present tree extraction in 3D images as a graph refinement task, of obtaining a subgraph from an over-complete input graph. To this end, we formulate an approximate Bayesian inference framework on undirected graphs using mean field approximation (MFA). Mean field networks are used for inference based on the interpretation that iterations of MFA can be seen as feed-forward operations in a neural network. This allows us to learn the model parameters from training data using back-propagation algorithm. We demonstrate usefulness of the model to extract airway trees from 3D chest CT data. We first obtain probability images using a voxel classifier that distinguishes airways from background and use Bayesian smoothing to model individual airway branches. This yields us joint Gaussian density estimates of position, orientation and scale as node features of the input graph. Performance of the method is compared with two methods: the first uses probability images from a trained voxel classifier with region growing, which is similar to one of the best performing methods at EXACT'09 airway challenge, and the second method is based on Bayesian smoothing on these probability images. Using centerline distance as error measure the presented method shows significant improvement compared to these two methods.


title: Inference on Local Average Treatment Effects for Misclassified Treatment

abstract: We develop point-identification for the local average treatment effect when the binary treatment contains a measurement error. The standard instrumental variable estimator is inconsistent for the parameter since the measurement error is non-classical by construction. We correct the problem by identifying the distribution of the measurement error based on the use of an exogenous variable that can even be a binary covariate. The moment conditions derived from the identification lead to generalized method of moments estimation with asymptotically valid inferences. Monte Carlo simulations and an empirical illustration demonstrate the usefulness of the proposed procedure.


title: Information theory for fields

abstract: A physical field has an infinite number of degrees of freedom since it has a field value at each location of a continuous space. Therefore, it is impossible to know a field from finite measurements alone and prior information on the field is essential for field inference. An information theory for fields is needed to join the measurement and prior information into probabilistic statements on field configurations. Such an information field theory (IFT) is built upon the language of mathematical physics, in particular on field theory and statistical mechanics. IFT permits the mathematical derivation of optimal imaging algorithms, data analysis methods, and even computer simulation schemes. The application of IFT algorithms to astronomical datasets provides high fidelity images of the Universe and facilitates the search for subtle statistical signals from the Big Bang. The concepts of IFT might even pave the road to novel computer simulations that are aware of their own uncertainties.


title: The restriction of $\sigma^{*}$ and $\phi$ on the protoneutron stars PSR J0348+0432

abstract: The restriction of $\sigma^{*}$ and $\phi$ on the protoneutron star (PNS) PSR J0348+0432 is described by the relativistic mean field theory (RMFT) through choosing the effective coupling constants. We use an entropy per baryon S=1 to describe thermal effect on PSR J0348+0432 in this work and compare the differences between PNS PSR J0348+0432 with $\sigma^{*}$,$\phi$ and without $\sigma^{*}$,$\phi$. These effects include the particle number distribution, the mass-radius relation, the moment of inertia and the surface gravitational redshift. The PNS PSR J0348+0432 with $\sigma^{*}$ and $\phi$ has more nucleons and will pushed forward the threshold for the appearance of the hyperons. The mass-radius relations are ($2.010M_{\odot}$,12.6520km) with $\sigma^{*}$ and $\phi$ and ($2.010M_{\odot}$,12.6170km) without $\sigma^{*}$ and $\phi$ respectively. The moment of inertia corresponding PNS PSR J0348+0432 are ($2.010M_{\odot}$,1.510$\times 10^{45}$g.cm$^{2}$) and ($2.010M_{\odot}$,1.559$\times 10^{45}$g.cm$^{2}$) respectively, the surface gravitational redshift corresponding PNS PSR J0348+0432 are ($2.010M_{\odot}$,0.3747) and ($2.010M_{\odot}$,0.3701) respectively. With the help of these calculations, we study the restriction of $\sigma^{*}$ and $\phi$ on the interactions between baryons in PNS core.


title: The KMTNet 2016 Data Release

abstract: We present Korea Microlensing Telescope Network (KMTNet) light curves for microlensing-event candidates for the 2016 season, which covers an area of 97 sq.deg observed at cadences ranging from Gamma=0.2/hr to Gamma=8/hr from three southern sites in Chile, South Africa, and Australia. These 2163 light curves are comprised of 1856 "clear microlensing" and 307 "possible microlensing" events (including 265 previously released from the K2 C9 field). The data policy is very similar to the one governing the 2015 release. The changes relative to 2015 in the algorithms to find and vet microlensing events are comprehensively described.


title: Evaluating the role of data quality when sharing information in hierarchical multi-stock assessment models, with an application to Dover Sole

abstract: An emerging approach to data-limited fisheries stock assessment uses hierarchical multi-stock assessment models to group stocks together, sharing information from data-rich to data-poor stocks. In this paper, we simulate data-rich and data-poor fishery and survey data scenarios for a complex of dover sole stocks. Simulated data for individual stocks were used to compare estimation performance for single-stock and hierarchical multi-stock versions of a Schaefer production model. The single-stock and best performing multi-stock models were then used in stock assessments for the real dover sole data. Multi-stock models often had lower estimation errors than single-stock models when assessment data had low statistical power. Relative errors for productivity and relative biomass parameters were lower for multi-stock assessment model configurations. In addition, multi-stock models that estimated hierarchical priors for survey catchability performed the best under data-poor scenarios. We conclude that hierarchical multi-stock assessment models are useful for data-limited stocks and could provide a more flexible alternative to data-pooling and catch only methods; however, these models are subject to non-linear side-effects of parameter shrinkage. Therefore, we recommend testing hierarchical multi-stock models in closed-loop simulations before application to real fishery management systems.


title: Decomposing qubit positive-operator valued measurements into continuous destructive weak measurements

abstract: It has been shown that any generalized measurement can be decomposed into a sequence of weak measurements corresponding to a stochastic process. However, the weak measurements may require almost arbitrary unitaries, which are unlikely to be realized by any real measurement device. Furthermore, many measurement processes are destructive, like photon counting procedures that terminate once all photons are consumed. One cannot expect to have full control of the evolution of a state under such destructive measurements, and the possible unitaries allow only a limited set of weak measurements. In this paper, we consider a qubit model of destructive weak measurements, which is a toy version of an optical cavity, in which the state of an electromagnetic field mode inside the cavity leaks out and is measured destructively while the vacuum state |0> leaks in to the cavity. At long times, the state of the qubit inevitably evolves to be |0>, and the only available control is the choice of measurement on the external ancilla system. Surprisingly, this very limited model can still perform an arbitrary projective measurement on the qubit in any basis, where the probability of getting an outcome satisfies the usual Born rule. Combining this method with probabilistic post processing, the result can be extended to any generalized measurement with commuting POVM elements. This implies, among other results, that any two-outcome POVM on a qubit can be decomposed into a sequence of destructive weak measurements by this restricted measurement device.


title: Discrete maximal regularity of an implicit Euler--Maruyama scheme with non-uniform time discretisation for a class of stochastic partial differential equations

abstract: An implicit Euler--Maruyama method with non-uniform step-size applied to a class of stochastic partial differential equations is studied. A spectral method is used for the spatial discretization and the truncation of the Wiener process. A discrete analogue of maximal $L^2$-regularity of the scheme and the discretised stochastic convolution is established, which has the same form as their continuous counterpart.


title: The Erdos-Moser sum-free set problem

abstract: We show that if A is a finite set of integers then it has a subset S of size \log^{1+c} |A| (c>0 absolute) such that s+s' is never in A when s and s' are distinct elements of S.


title: A Manga-Driven System Requirements Development PBL Exercise

abstract: We conducted a Project-Based Learning (PBL)-type exercise incorporating Japanese cartoon (manga) techniques into Requirements Development (RD) processes. Manga has established techniques, such as those for character setting and story development, that we thought are also valid for RD processes. Using this manga-driven method, students were able to clarify high-level project goals early in the development life-cycle, and succeeded in defining high quality and unique system ideas.


title: Curvilinear Mesh Adaptation using Radial Basis Function Interpolation and Smoothing

abstract: We present a new iterative technique based on radial basis function (RBF) interpolation and smoothing for the generation and smoothing of curvilinear meshes from straight-sided or other curvilinear meshes. Our technique approximates the coordinate deformation maps in both the interior and boundary of the curvilinear output mesh by using only scattered nodes on the boundary of the input mesh as data sites in an interpolation problem. Our technique produces high-quality meshes in the deformed domain even when the deformation maps are singular due to a new iterative algorithm based on modification of the RBF shape parameter. Due to the use of RBF interpolation, our technique is applicable to both 2D and 3D curvilinear mesh generation without significant modification.


title: Vertex algebras and coordinate rings of semi-infinite flags

abstract: The direct sum of irreducible level one integrable representations of affine Kac-Moody Lie algebra of (affine) type $ADE$ carries a structure of $P/Q$-graded vertex operator algebra. There exists a filtration on this direct sum studied by Kato and Loktev such that the corresponding graded vector space is a direct sum of global Weyl modules. The associated graded space with respect to the dual filtration is isomorphic to the homogenous coordinate ring of semi-infinite flag variety. We describe the ring structure in terms of vertex operators and endow the homogenous coordinate ring with a structure of $P/Q$-graded vertex operator algebra. We use the vertex algebra approach to derive semi-infinite Pl\"ucker-type relations in the homogeneous coordinate ring.


title: Reference-Conditioned Super-Resolution by Neural Texture Transfer

abstract: With the recent advancement in deep learning, we have witnessed a great progress in single image super-resolution. However, due to the significant information loss of the image downscaling process, it has become extremely challenging to further advance the state-of-the-art, especially for large upscaling factors. This paper explores a new research direction in super resolution, called reference-conditioned super-resolution, in which a reference image containing desired high-resolution texture details is provided besides the low-resolution image. We focus on transferring the high-resolution texture from reference images to the super-resolution process without the constraint of content similarity between reference and target images, which is a key difference from previous example-based methods. Inspired by recent work on image stylization, we address the problem via neural texture transfer. We design an end-to-end trainable deep model which generates detail enriched results by adaptively fusing the content from the low-resolution image with the texture patterns from the reference image. We create a benchmark dataset for the general research of reference-based super-resolution, which contains reference images paired with low-resolution inputs with varying degrees of similarity. Both objective and subjective evaluations demonstrate the great potential of using reference images as well as the superiority of our results over other state-of-the-art methods.


title: The Effect of a Light Sterile Neutrino at NO$\nu$A and DUNE

abstract: Now that the NOvA experiment is approaching two years run time and has released some preliminary data, some constraints for the oscillation parameters can be inferred. Currently the best fits for NOvA are three separate results, the reason they are indistinct is that they produce almost degenerate probability curves. It has been postulated that these degeneracies can be resolved by running antineutrinos at NOvA and/or combining its data with T2K. However, this degeneracy resolution power can be compromised if sterile neutrinos are present due to additional degrees of freedom that can significantly alter the oscillation probability for any of the current best fits. We aim to investigate this degradation in predictive power and the effect of the DUNE experiment on it. Now that the NOvA experiment has been running for a few years and has released some preliminary data, some constraints for the oscillation parameters can be inferred. The best fits for NOvA include three degenerate results, the reason they are indistinct is that they produce almost degenerate probability curves. It has been postulated that these degeneracies can be resolved by running antineutrinos at NOvA and/or combining its data with T2K. However, this degeneracy resolution power can be compromised if sterile neutrinos are present due to additional degrees of freedom that can significantly alter the oscillation probability for any of these best fits. We aim to investigate this degradation in predictive power and the effect of the DUNE experiment on it. In light of the 2018 NOvA data we also consider the same fits but with theta_{23}=45 to see if the sensitivity results are different.


title: What's my age?: Predicting Twitter User's Age using Influential Friend Network and DBpedia

abstract: Social media is a rich source of user behavior and opinions. Twitter senses nearly 500 million tweets per day from 328 million users.An appropriate machine learning pipeline over this information enables up-to-date and cost-effective data collection for a wide variety of domains such as; social science, public health, the wisdom of the crowd, etc. In many of the domains, users demographic information is key to the identification of segments of the populations being studied. For instance, Which age groups are observed to abuse which drugs?, Which ethnicities are most affected by depression per location?. Twitter in its current state does not require users to provide any demographic information. We propose to create a machine learning system coupled with the DBpedia graph that predicts the most probable age of the Twitter user. In our process to build an age prediction model using social media text and user meta-data, we explore the existing state of the art approaches. Detailing our data collection, feature engineering cycle, model selection and evaluation pipeline, we will exhibit the efficacy of our approach by comparing with the "predict mean" age estimator baseline.


title: Error Estimation of Numerical Solvers for Linear Ordinary Differential Equations

abstract: Solving Linear Ordinary Differential Equations (ODEs) plays an important role in many applications. There are various numerical methods and solvers to obtain approximate solutions. However, few work about global error estimation can be found in the literature. In this paper, we first give a definition of the residual, based on the piecewise Hermit interpolation, which is a kind of the backward-error of ODE solvers. It indicates the reliability and quality of numerical solution. Secondly, the global error between the exact solution and an approximate solution is the forward error and a bound of it can be given by using the backward-error. The examples in the paper show that our estimate works well for a large class of ODE models.


title: High-temperature operation of a silicon qubit

abstract: This study alleviates the low operating temperature constraint of Si qubits. A qubit is a key element for quantum sensors, memories, and computers. Electron spin in Si is a promising qubit, as it allows both long coherence times and potential compatibility with current silicon technology. Si qubits have been implemented using gate-defined quantum dots or shallow impurities. However, operation of Si qubits has been restricted to milli-Kelvin temperatures, thus limiting the application of the quantum technology. In this study, we addressed a single deep impurity, having strong electron confinement of up to 0.3 eV, using single-electron tunnelling transport. We also achieved qubit operation at 5-10 K through a spin-blockade effect based on the tunnelling transport via two impurities. The deep impurity was implemented by tunnel field-effect transistors (TFETs) instead of conventional FETs. With further improvement in fabrication and controllability, this work presents the possibility of operating silicon spin qubits at elevated temperatures.


title: A Quantitative Comparison of Opacities Calculated Using the Distorted- Wave and $\boldsymbol{R}$-Matrix Methods

abstract: The present debate on the reliability of astrophysical opacities has reached a new climax with the recent measurements of Fe opacities on the Z-machine at the Sandia National Laboratory \citep{Bailey2015}. To understand the differences between theoretical results, on the one hand, and experiments on the other, as well as the differences among the various theoretical results, detailed comparisons are needed. Many ingredients are involved in the calculation of opacities; deconstructing the whole process and comparing the differences at each step are necessary to quantify their importance and impact on the final results. We present here such a comparison using the two main approaches to calculate the required atomic data, the $R$-Matrix and distorted-wave methods, as well as sets of configurations and coupling schemes to quantify the effects on the opacities for the $Fe\ XVII$ and $Ni\ XIV$ ions.


title: Testing equality of spectral density operators for functional linear processes

abstract: The problem of testing equality of the entire second order structure of two independent functional linear processes is considered. A fully functional $L^2$-type test is developed which evaluates, over all frequencies, the Hilbert-Schmidt distance between the estimated spectral density operators of the two processes. The asymptotic behavior of the test statistic is investigated and its limiting distribution under the null hypothesis is derived. Furthermore, a novel frequency domain bootstrap method is developed which approximates more accurately the distribution of the test statistic under the null than the large sample Gaussian approximation obtained. Asymptotic validity of the bootstrap procedure is established and consistency of the bootstrap-based test under the alternative is proved. Numerical simulations show that, even for small samples, the bootstrap-based test has very good size and power behavior. An application to meteorological functional time series is also presented.


title: Spectral theory of pseudo-differential operators of degree 0 and application to forced linear waves

abstract: We extend the results of our paper "Attractors for two dimensional waves with homogeneous Hamiltonians of degree 0" written with Laure Saint-Raymond to the case of forced linear wave equations in any dimension. We prove that, in dimension 2,if the foliation on the boundary at infinity of the energy shell is Morse-Smale, we can apply Mourre's theory and hence get the asymptotics of the forced solution. We also characterize the wavefrontsets of the limit Schwartz distribution using radial propagation estimates.


title: Learning Deep Gradient Descent Optimization for Image Deconvolution

abstract: As an integral component of blind image deblurring, non-blind deconvolution removes image blur with a given blur kernel, which is essential but difficult due to the ill-posed nature of the inverse problem. The predominant approach is based on optimization subject to regularization functions that are either manually designed, or learned from examples. Existing learning based methods have shown superior restoration quality but are not practical enough due to their restricted and static model design. They solely focus on learning a prior and require to know the noise level for deconvolution. We address the gap between the optimization-based and learning-based approaches by learning a universal gradient descent optimizer. We propose a Recurrent Gradient Descent Network (RGDN) by systematically incorporating deep neural networks into a fully parameterized gradient descent scheme. A hyper-parameter-free update unit shared across steps is used to generate updates from the current estimates, based on a convolutional neural network. By training on diverse examples, the Recurrent Gradient Descent Network learns an implicit image prior and a universal update rule through recursive supervision. The learned optimizer can be repeatedly used to improve the quality of diverse degenerated observations. The proposed method possesses strong interpretability and high generalization. Extensive experiments on synthetic benchmarks and challenging real-world images demonstrate that the proposed deep optimization method is effective and robust to produce favorable results as well as practical for real-world image deblurring applications.


title: Scale invariance and constants of motion

abstract: Scale invariance in the theory of classical mechanics can be induced from the scale invariance of background fields. In this paper we consider the relation between the scale invariance and the constants of particle motion in a self-similar spacetime, only in which the symmetry is well-defined and is generated by a homothetic vector. Relaxing the usual conservation condition by the Hamiltonian constraint in a particle system, we obtain a conservation law holding only on the constraint surface in the phase space. By the conservation law, we characterize constants of motion associated with the scale invariance not only for massless particles but for massive particles and classify the condition for the existence of the constants of motion. Furthermore, we find the explicit form of the constants of motion by solving the conservation equations.


title: X-ray ghost tomography: denoising, dose fractionation and mask considerations

abstract: Ghost imaging has recently been successfully achieved in the X-ray regime; due to the penetrating power of X-rays this immediately opens up the possibility of X-ray ghost tomography. No research into this topic currently exists in the literature. Here we present adaptations of conventional tomography techniques to this new ghost imaging scheme. Several numerical implementations for tomography through X-ray ghost imaging are considered. Specific attention is paid to schemes for denoising of the resulting tomographic reconstruction, issues related to dose fractionation, and considerations regarding the ensemble of illuminating masks used for ghost imaging. Each theme is explored through a series of numerical simulations, and several suggestions offered for practical realisations of X-ray ghost tomography.


title: A model of the late universe with viscous Zel'ldovich fluid and decaying vacuum

abstract: Many have speculated about the presence of a stiff fluid in very early stage of the universe. Such a stiff fluid was first introduced by Zel'dovich. Recently the late acceleration of the universe was studied by taking bulk viscous stiff fluid as the dominant cosmic component, but the age predicted by such a model is less than the observed value. We consider a flat universe with viscous stiff fluid and decaying vacuum energy as the cosmic components and found that the model predicts a reasonable background evolution of the universe with de Sitter epoch as end phase of expansion. More over the model also predicts a reasonable value for the age of the present universe. We also perform a dynamical system analysis of the model and found that the end de Sitter phase predicted by the model is stable.


title: Realtime Active Sound Source Localization for Unmanned Ground Robots Using a Self-Rotational Bi-Microphone Array

abstract: This work presents a novel technique that performs both orientation and distance localization of a sound source in a three-dimensional (3D) space using only the interaural time difference (ITD) cue, generated by a newly-developed self-rotational bi-microphone robotic platform. The system dynamics is established in the spherical coordinate frame using a state-space model. The observability analysis of the state-space model shows that the system is unobservable when the sound source is placed with elevation angles of $90$ and $0$ degree. The proposed method utilizes the difference between the azimuth estimates resulting from respectively the 3D and the two-dimensional models to check the zero-degree-elevation condition and further estimates the elevation angle using a polynomial curve fitting approach. Also, the proposed method is capable of detecting a $90$-degree elevation by extracting the zero-ITD signal 'buried' in noise. Additionally, a distance localization is performed by first rotating the microphone array to face toward the sound source and then shifting the microphone perpendicular to the source-robot vector by a predefined distance of a fixed number of steps. The integrated rotational and translational motions of the microphone array provide a complete orientation and distance localization using only the ITD cue. A novel robotic platform using a self-rotational bi-microphone array was also developed for unmanned ground robots performing sound source localization. The proposed technique was first tested in simulation and was then verified on the newly-developed robotic platform. Experimental data collected by the microphones installed on a KEMAR dummy head were also used to test the proposed technique. All results show the effectiveness of the proposed technique.


title: A general range-separated double-hybrid density-functional theory

abstract: A range-separated double-hybrid (RSDH) scheme which generalizes the usual range-separated hybrids and double hybrids is developed. This scheme consistently uses a two-parameter Coulomb-attenuating-method (CAM)-like decomposition of the electron-electron interaction for both exchange and correlation in order to combine Hartree-Fock exchange and second-order M{{\o}}ller-Plesset (MP2) correlation with a density functional. The RSDH scheme relies on an exact theory which is presented in some detail. Several semi-local approximations are developed for the short-range exchange-correlation density functional involved in this scheme. After finding optimal values for the two parameters of the CAM-like decomposition, the RSDH scheme is shown to have a relatively small basis dependence and to provide atomization energies, reaction barrier heights, and weak intermolecular interactions globally more accurate or comparable to range-separated MP2 or standard MP2. The RSDH scheme represents a new family of double hybrids with minimal empiricism which could be useful for general chemical applications.


title: Sharp Beckner-type inequalities for Cauchy and spherical distributions

abstract: Using some harmonic extensions on the upper-half plane, and probabilistic representations, and curvature-dimension inequalities with some negative dimensions, we obtain some new opimal functional inequalities of the Beckner type for the Cauchy type distributions on the Euclidean space. These optimal inequalities appear to be equivalent to some non tight optimal Beckner inequalities on the sphere, and the family appear to be a new form of the Sobolev inequality.


title: The imaginary part of the scattering Green function: monochromatic relations to the real part and uniqueness results for inverse problems

abstract: For many wave propagation problems with random sources it has been demonstrated that cross correlations of wave fields are proportional to the imaginary part of the Green function of the underlying wave equation. This leads to the inverse problem to recover coefficients of a wave equation from the imaginary part of the Green function on some measurement manifold. In this paper we prove, in particular, local uniqueness results for the Schr{\"o}dinger equation with one frequency and for the acoustic wave equation with unknown density and sound speed and two frequencies. As the main tool of our analysis, we establish new algebraic identities between the real and the imaginary part of Green's function, which in contrast to the well-known Kramers-Kronig relations involve only one frequency.


title: Local uniqueness of $m$-bubbling sequences for the Gel'fand equation

abstract: We consider the Gel'fand problem, $$ \begin{cases} \Delta w_{\varepsilon}+\varepsilon^2 h e^{w_{\varepsilon}}=0\quad&\mbox{in}\quad\Omega, w_{\varepsilon}=0\quad&\mbox{on}\quad\partial\Omega, \end{cases} $$ where $h$ is a nonnegative function in ${\Omega\subset\mathbb{R}^2}$. Under suitable assumptions on $h$ and $\Omega$, we prove the local uniqueness of $m-$bubbling solutions for any $\varepsilon>0$ small enough.


title: Ground-based detection of G star superflares with NGTS

abstract: We present high cadence detections of two superflares from a bright G8 star (V = 11.56) with the Next Generation Transit Survey (NGTS). We improve upon previous superflare detections by resolving the flare rise and peak, allowing us to fit a solar flare inspired model without the need for arbitrary break points between rise and decay. Our data also enables us to identify substructure in the flares. From changing starspot modulation in the NGTS data we detect a stellar rotation period of 59 hours, along with evidence for differential rotation. We combine this rotation period with the observed \textit{ROSAT} X-ray flux to determine that the star's X-ray activity is saturated. We calculate the flare bolometric energies as $5.4^{+0.8}_{-0.7}\times10^{34}$ and $2.6^{+0.4}_{-0.3}\times10^{34}$ erg and compare our detections with G star superflares detected in the \textit{Kepler} survey. We find our main flare to be one of the largest amplitude superflares detected from a bright G star. With energies more than 100 times greater than the Carrington event, our flare detections demonstrate the role that ground-based instruments such as NGTS can have in assessing the habitability of Earth-like exoplanets, particularly in the era of \textit{PLATO}.


title: Maximum likelihood estimation for Gaussian processes under inequality constraints

abstract: We consider covariance parameter estimation for a Gaussian process under inequality constraints (boundedness, monotonicity or convexity) in fixed-domain asymptotics. We address the estimation of the variance parameter and the estimation of the microergodic parameter of the Mat\'ern and Wendland covariance functions. First, we show that the (unconstrained) maximum likelihood estimator has the same asymptotic distribution, unconditionally and conditionally to the fact that the Gaussian process satisfies the inequality constraints. Then, we study the recently suggested constrained maximum likelihood estimator. We show that it has the same asymptotic distribution as the (unconstrained) maximum likelihood estimator. In addition, we show in simulations that the constrained maximum likelihood estimator is generally more accurate on finite samples. Finally, we provide extensions to prediction and to noisy observations.