title: New Limits for Large $N$ Matrix and Tensor Models: Large $D$, Melons and Applications

abstract: Large $N$ matrix models play an important role in modern theoretical physics, ranging from quantum chromodynamics to string theory and holography. However, they remain a difficult technical challenge because in most cases it is not known how to perform the sum over planar graphs, which dominate the models at large $N$. In this thesis, we study large $D$ matrix models, which provide a framework to build new limits for matrix models in which the sum over planar graphs simplifies when $D$ is large. The basic degrees of freedom are real matrices of size $N\times N$ with $r$ additional indices of range $D$. They can be interpreted as a real tensor of rank $R=r+2$ with indices of different ranges, making a compelling connection with tensor models. We define a new large $D$ limit for the sum over Feynman graphs of fixed genus in matrix models, based on an enhanced large $D$ scaling of the coupling constants. Then, we show that the resulting large $D$ expansion is well-defined and organized according to a half-integer called the index. When $N=D$, the result provides a new large $N$ limit for general $\text{O}(N)^R$ invariant tensor models. In the large $D$ limit, the sum over planar graphs of large $N$ matrix models simplifies to a non-trivial sum over generalized melonic graphs. This class of graphs extends the one obtained in tensor models with standard scaling and allows for a wider class of interactions, including all the maximally single-trace terms. The general classification of generalized melonic graphs remains an open problem. However, in the case of the complete interaction of order $R+1$ for $R$ a prime number, we identify them in detail and demonstrate that they exhibit the same important features as the SYK model with $q=(R+1)$-fold random interactions, including the emergent conformal symmetry in the infrared regime and maximal chaos.


title: Study of thermal SU(3) supersymmetric Yang-Mills theory and near-conformal theories from the gradient flow

abstract: We compute the renormalization group flow of the mass anomalous dimension in adjoint QCD with $N_{f}=1$, $3/2$, and 2 Dirac fermions, using the gradient flow. Preliminary results are in agreement with at least a near-conformal scenario in all cases. At the largest flavor numbers we obtain the strongest indication for an IR conformal fixed point scenario. Moreover, we provide results for the thermal phase transitions in SU(3) supersymmetric Yang-Mills theory. We find hints for a connection between chiral and center symmetries in terms of a single first order phase transition where chiral symmetry is restored and center symmetry gets broken.


title: FusionStitching: Boosting Execution Efficiency of Memory Intensive Computations for DL Workloads

abstract: Performance optimization is the art of continuous seeking a harmonious mapping between the application domain and hardware. Recent years have witnessed a surge of deep learning (DL) applications in industry. Conventional wisdom for optimizing such workloads mainly focus on compute intensive ops (GEMM, Convolution, etc). Yet we show in this work, that the performance of memory intensive computations is vital to E2E performance in practical DL models. We propose \emph{FusionStitching}, a optimization framework capable of fusing memory intensive \emph{elementwise}, \emph{reduction} and fine grained \emph{GEMM/Batched-GEMM} ops, with or without data dependences, into large computation units, then mapping and transforming them into efficient GPU kernels. We formulate the fusion plan optimization as an integer linear programming (ILP) problem, and propose a set of empirical heuristics to reduce the combinatorial search space. In order to map optimized fusion plans to hardware, we propose a technique to effectively compose various groups of computations into a single GPU kernel, by fully leveraging on chip resources like scratchpads or registers. Experimental results on six benchmarks and four industry scale practical models are encouraging. Overall, \emph{FusionStitching} can reach up to 5.7x speedup compared to Tensorflow baseline, and achieves 1.25x to 1.85x performance speedups compared to current state of the art, with 1.4x on average (geometric mean).


title: Scintillating crystals for the Neutral Particle Spectrometer in Hall C at JLab

abstract: This paper discusses the quality and performance of currently available PbWO$_4$ crystals of relevance to high-resolution electromagnetic calorimetry, e.g. detectors for the Neutral Particle Spectrometer at Jefferson Lab or those planned for the Electron-Ion Collider. Since the construction of the Compact Muon Solenoid (CMS) at the Large Hadron Collider (LHC) and early PANDA (The antiProton ANnihilations at DArmstadt) electromagnetic calorimeter (ECAL) the worldwide availability of high quality PbWO$_4$ production has changed dramatically. We report on our studies of crystal samples from SICCAS/China and CRYTUR/Czech Republic that were produced between 2014 and 2019.


title: Dr. Bertlmann's Socks in the Quaternionic World of Ambidextral Reality

abstract: In this pedagogical paper, John S. Bell's amusing example of Dr. Bertlmann's socks is reconsidered, first within a toy model of a two-dimensional one-sided world of a non-orientable M\"obius strip, and then within a real world of three-dimensional quaternionic sphere, S^3, which results from an addition of a single point to R^3 at infinity. In the quaternionic world, which happens to be the spatial part of a solution of Einstein's field equations of general relativity, the singlet correlations between a pair of entangled fermions can be understood as classically as those between Dr. Bertlmann's colorful socks.


title: The Huntsman Telescope

abstract: The Huntsman Telescope, located at Siding Spring Observatory in Australia, is a system of ten telephoto Canon lenses designed for low surface brightness imaging in the Southern sky. Based upon the Dragonfly Telephoto Array, the refractive lens-based system provides an obstruction free optical path, which reduces the number of scattering surfaces and allows easier access to lower surface brightness levels. In this proceeding, we present an analysis of the impact of flat fielding uncertainty on the limiting low surface brightness levels. We show that a fairly standard set of flat-field data can be well-characterised to a $\sim0.1\%$ level. This corresponds to a 5-$\sigma$ lower limit of $\sim33$ magnitude per arcsecond$^2$, which means that flat fielding is not likely going to set Huntsman's low surface brightness limit. We also present early results of an exoplanet transient mode for Huntsman where all lenses work together to detect subtle variations in the luminosity of relatively bright $V=8-12$ magnitude stars. High-precision exoplanet imaging is ultimately limited by systematic uncertainties, so we anticipate multiple lenses will help to mitigate issues related to pixel-to-pixel and intra-pixel sensitivity variations. Our initial results show we can easily get $\sim0.4\%$ photometric precision with a single, defocused lens.


title: Measurements of electron transport in liquid and gas Xenon using a laser-driven photocathode

abstract: Measurements of electron drift properties in liquid and gaseous xenon are reported. The electrons are generated by the photoelectric effect in a semi-transparent gold photocathode driven in transmission mode with a pulsed ultraviolet laser. The charges drift and diffuse in a small chamber at various electric fields and a fixed drift distance of 2.0 cm. At an electric field of 0.5 kV/cm, the measured drift velocities and corresponding temperature coefficients respectively are $1.97 \pm 0.04$ mm/$\mu$s and $(-0.69\pm0.05)$\%/K for liquid xenon, and $1.42 \pm 0.03$ mm/$\mu$s and $(+0.11\pm0.01)$\%/K for gaseous xenon at 1.5 bar. In addition, we measure longitudinal diffusion coefficients of $25.7 \pm 4.6$ cm$^2$/s and $149 \pm 23$ cm$^2$/s, for liquid and gas, respectively. The quantum efficiency of the gold photocathode is studied at the photon energy of 4.73 eV in liquid and gaseous xenon, and vacuum. These charge transport properties and the behavior of photocathodes in a xenon environment are important in designing and calibrating future large scale noble liquid detectors.


title: Histogram Transform Ensembles for Density Estimation

abstract: We investigate an algorithm named histogram transform ensembles (HTE) density estimator whose effectiveness is supported by both solid theoretical analysis and significant experimental performance. On the theoretical side, by decomposing the error term into approximation error and estimation error, we are able to conduct the following analysis: First of all, we establish the universal consistency under $L_1(\mu)$-norm. Secondly, under the assumption that the underlying density function resides in the H\"{o}lder space $C^{0,\alpha}$, we prove almost optimal convergence rates for both single and ensemble density estimators under $L_1(\mu)$-norm and $L_{\infty}(\mu)$-norm for different tail distributions, whereas in contrast, for its subspace $C^{1,\alpha}$ consisting of smoother functions, almost optimal convergence rates can only be established for the ensembles and the lower bound of the single estimators illustrates the benefits of ensembles over single density estimators. In the experiments, we first carry out simulations to illustrate that histogram transform ensembles surpass single histogram transforms, which offers powerful evidence to support the theoretical results in the space $C^{1,\alpha}$. Moreover, to further exert the experimental performances, we propose an adaptive version of HTE and study the parameters by generating several synthetic datasets with diversities in dimensions and distributions. Last but not least, real data experiments with other state-of-the-art density estimators demonstrate the accuracy of the adaptive HTE algorithm.


title: DDNet: Dual-path Decoder Network for Occlusion Relationship Reasoning

abstract: Occlusion relationship reasoning based on convolution neural networks consists of two subtasks: occlusion boundary extraction and occlusion orientation inference. Due to the essential differences between the two subtasks in the feature expression at the higher and lower stages, it is challenging to carry on them simultaneously in one network. To address this issue, we propose a novel Dual-path Decoder Network, which uniformly extracts occlusion information at higher stages and separates into two paths to recover boundary and occlusion orientation respectively in lower stages. Besides, considering the restriction of occlusion orientation presentation to occlusion orientation learning, we design a new orthogonal representation for occlusion orientation and proposed the Orthogonal Orientation Regression loss which can get rid of the unfitness between occlusion representation and learning and further prompt the occlusion orientation learning. Finally, we apply a multi-scale loss together with our proposed orientation regression loss to guide the boundary and orientation path learning respectively. Experiments demonstrate that our proposed method achieves state-of-the-art results on PIOD and BSDS ownership datasets.


title: Near-resonant light transmission in two-dimensional dense cold atomic media with short-range positional correlations

abstract: Light propagation in disordered media is a fundamental and important problem in optics and photonics. In particular, engineering light-matter interaction in disordered cold atomic ensembles is one of the central topics in modern quantum and atomic optics. The collective response of dense atomic gases under light excitation, which crucially depends on the spatial distribution of atoms and the geometry of the ensemble, has important impacts on quantum technologies like quantum sensors, atomic clocks and quantum information storage. Here we analyze near-resonant light transmission in two-dimensional dense ultracold atomic ensembles with short-range positional correlations. Based on the coupled-dipole simulations under different atom number densities and correlation lengths, we show that the collective effects are strongly influenced by those positional correlations, manifested as significant shifts and broadening or narrowing of transmission resonance lines. The results show that mean-field theories like Lorentz-Lorenz relation are not capable of describing such collective effects. We further investigate the statistical distribution of eigenstates, which are significantly affected by the interplay between dipole-dipole interactions and position correlations. This work can provide profound implications on collective and cooperative effects in cold atomic ensembles as well as the study of mesoscopic physics concerning light transport in strongly scattering disordered media.

categories: physics.optics cond-mat.dis-nn cond-mat.quant-gas quant-ph


title: Flexible graph matching and graph edit distance using answer set programming

abstract: The graph isomorphism, subgraph isomorphism, and graph edit distance problems are combinatorial problems with many applications. Heuristic exact and approximate algorithms for each of these problems have been developed for different kinds of graphs: directed, undirected, labeled, etc. However, additional work is often needed to adapt such algorithms to different classes of graphs, for example to accommodate both labels and property annotations on nodes and edges. In this paper, we propose an approach based on answer set programming. We show how each of these problems can be defined for a general class of property graphs with directed edges, and labels and key-value properties annotating both nodes and edges. We evaluate this approach on a variety of synthetic and realistic graphs, demonstrating that it is feasible as a rapid prototyping approach.


title: Photonic realization of quantum resetting

abstract: Contrary to the usual assumption of at least partial control of quantum dynamics, a surprising recent result proved that an arbitrary quantum state can be probabilistically reset to a state in the past by having it interact with probing systems in a consistent, but $uncontrolled$ way. We present a photonic implementation to achieve this resetting process, experimentally verifying that a state can be probabilistically reset to its past with a fidelity of $0.870\pm0.012$. We further demonstrate the preservation of an entangled state, which still violates a Bell inequality, after half of the entangled pair was reset. The ability to reset uncontrolled quantum states has implications in the foundations of quantum physics and applications in areas of quantum technology.


title: Cascade disruptions in asteroid clusters

abstract: We studied asteroid clusters suggesting a possibility of at least two disruption events in their recent history (< 5 Myr). We searched for new members of known asteroid pairs and clusters and we verified their membership. We found four asteroid clusters, namely the clusters of (11842) Kap'bos, (14627) Emilkowalski, (63440) 2001 MD30 and (157123) 2004 NW5 that show at least two secondary separation events that occurred at significantly different times. We considered a possible formation mechanism for these clusters: The parent of an asteroid cluster was spun up to its critical rotation frequency, underwent a rotation fission and was slowed down by escape of the newly formed secondary/ies. Then the YORP effect spun up the primary again and it reached its critical rotation frequency and underwent another fission. We created a simple model to test whether this scenario is possible for the four clusters. We obtained a good agreement between the model and the cluster properties for the clusters of Kap'bos and (63440). For the cluster of Emilkowalski, the time needed for the primary to reach its critical frequency after the first fission event was predicted to be too long by a factor of several. We suspect, considering also its D type taxonomic classification and the existence of a dust band associated with the cluster, that the asteroid Emilkowalski may actually be a cometary nucleus. For the cluster of (157123), the final rotational frequency of the primary after the last fission event predicted by our model is in a good agreement with the observed rotation frequency of (157123). However, a separation of the older secondary is not possible in our model due to the deficiency of free energy needed for an escape of the large secondary. This could be due to an error in the H value of the secondary or the possibility that we did not find the real primary of this cluster.


title: Bi-orbital sheaves and affine Hecke algebras at roots of unity

abstract: Continuing the study of perverse sheaves on the nilpotent cone of a $\mathbb{Z}/m$-graded Lie algebra initiated by Lusztig--Yun, we study in this work the parabolic induction and introduce the notion of supercuspidal sheaves on the nilpotent cone. One of our main results shows that simple perverse sheaves with nilpotent singular support (called bi-orbital sheaves) are produced by parabolic induction from supercuspidal sheaves. As application, we provide a proof for a theorem announced by I. Grojnowski on the parametrisation of simple modules of affine Hecke algebras at roots of unity via Deligne--Langlands--Lusztig parameters with nilpotent singular support.


title: The Diffuse Ionized Gas Halo of the Small Magellanic Cloud

abstract: Observations with the Wisconsin H-alpha Mapper (WHAM) reveal a large, diffuse ionized halo that surrounds the Small Magellanic Cloud (SMC). We present the first kinematic H-alpha survey of an extended region around the galaxy, from (l,b) = (289.5,-35.0) to (315.1,-5.3) and covering +90 <= vLSR <= +210 km s-1. The ionized gas emission extends far beyond the central stellar component of the galaxy, reaching similar distances to that of the diffuse neutral halo traced by 21 cm observations. H-alpha emission extends several degrees beyond the sensitivity of current H I surveys toward smaller Galactic longitudes and more negative Galactic latitudes. The velocity field of the ionized gas near the SMC appears similar to to the neutral halo of the galaxy. Using the observed emission measure as a guide, we estimate the mass of this newly revealed ionized component to be roughly (0.8 - 1.0) x 10^9 M_sun, which is comparable to the total neutral mass in the same region of (0.9 - 1.1) x 10^9 M_sun. We find ratios of the total ionized gas mass divided by the total neutral plus ionized gas mass in three distinct subregions to be: (1) 46%-54% throughout the SMC and its extended halo, (2) 12%-32% in the SMC Tail that extends toward the Magellanic Bridge, and (3) 65%-79% in a filament that extends away from the SMC toward the Magellanic Stream. This newly discovered, coherent H-alpha filament does not appear to have a well-structured neutral component and is also not coincident with two previously identified filaments traced by 21 cm emission within the Stream.


title: Complexity bounds on supermesh construction for quasi-uniform meshes

abstract: Projecting fields between different meshes commonly arises in computational physics. This operation requires a supermesh construction and its computational cost is proportional to the number of cells of the supermesh $n$. Given any two quasi-uniform meshes of $n_A$ and $n_B$ cells respectively, we show under standard assumptions that n is proportional to $n_A + n_B$. This result substantially improves on the best currently available upper bound on $n$ and is fundamental for the analysis of algorithms that use supermeshes.


title: Time Travel Paradoxes and Multiple Histories

abstract: If time travel is possible, it seems to inevitably lead to paradoxes. These include consistency paradoxes, such as the famous grandfather paradox, and bootstrap paradoxes, where something is created out of nothing. One proposed class of resolutions to these paradoxes allows for multiple histories (or timelines), such that any changes to the past occur in a new history, independent of the one where the time traveler originated. We introduce a simple mathematical model for a spacetime with a time machine, and suggest two possible multiple-histories models, making use of branching spacetimes and covering spaces respectively. We use these models to construct novel and concrete examples of multiple-histories resolutions to time travel paradoxes, and we explore questions such as whether one can ever come back to a previously visited history and whether a finite or infinite number of histories is required. Interestingly, we find that the histories may be finite and cyclic under certain assumptions, in a way which extends the Novikov self-consistency conjecture to multiple histories and exhibits hybrid behavior combining the two. Investigating these cyclic histories, we rigorously determine how many histories are needed to fully resolve time travel paradoxes for particular laws of physics. Finally, we discuss how observers may experimentally distinguish between multiple histories and the Hawking and Novikov conjectures.


title: Uniqueness and Stability of Solutions for a Coupled System of Nabla Fractional Difference Boundary Value Problems

abstract: In this article, we obtain sufficient conditions on existence, uniqueness and Ulam--Hyers stability of solutions for a coupled system of two-point nabla fractional difference boundary value problems, using Banach fixed point theorem and Urs's approach. Finally, we illustrate the applicability of established results through an example.


title: Transaction Confirmation Time Prediction in Ethereum Blockchain Using Machine Learning

abstract: Blockchain offers a decentralized, immutable, transparent system of records. It offers a peer-to-peer network of nodes with no centralised governing entity making it unhackable and therefore, more secure than the traditional paper-based or centralised system of records like banks etc. While there are certain advantages to the paper-based recording approach, it does not work well with digital relationships where the data is in constant flux. Unlike traditional channels, governed by centralized entities, blockchain offers its users a certain level of anonymity by providing capabilities to interact without disclosing their personal identities and allows them to build trust without a third-party governing entity. Due to the aforementioned characteristics of blockchain, more and more users around the globe are inclined towards making a digital transaction via blockchain than via rudimentary channels. Therefore, there is a dire need for us to gain insight on how these transactions are processed by the blockchain and how much time it may take for a peer to confirm a transaction and add it to the blockchain network. This paper presents a novel approach that would allow one to estimate the time, in block time or otherwise, it would take for a mining node to accept and confirm a transaction to a block using machine learning. The paper also aims to compare the predictive accuracy of two machine learning regression models- Random Forest Regressor and Multilayer Perceptron against previously proposed statistical regression model under a set evaluation criterion. The objective is to determine whether machine learning offers a more accurate predictive model than conventional statistical models. The proposed model results in improved accuracy in prediction.


title: Quantum Effects in Gravity Waves

abstract: We discuss the quantum mechanical description of a gravitational wave interacting with a cavity electromagnetic field. Quantum fluctuations of the gravitational vacuum induce squeezing in the optical field. Moreover, this squeezing experiences revivals, a purely quantum effect. Measuring these gravitationally induced revivals, although out of reach from experiments, would provide evidence on the quantum nature of gravity. We also discuss the quantum mechanical treatment of the interaction between coherent and squeezed gravitational wave states and a gravity wave detector. In the case of a coherent gravitational wave, we reproduce the result from the classical theory with a quantum mechanical calculation. The case of a squeezed gravity wave is not calculable within the classical theory, and could provide evidence on the quantum nature of gravity.


title: Etude des effets non lin\'eaires observ\'es sur les oscillations d'un pendule simple

abstract: In this paper we present a study of the non-linear effects of anharmonicity of the potential of the simple pendulum. In a theoretical reminder we highlight that anharmonicity of the potential generates additional harmonics and the non-isochronism of oscillations. These phenomena are all the more important as we move away from the oscillations at small angles, which represent the domain of validity of the harmonic approximation. The measurement is apprehended by means of the acquisition box SYSAM-SP5 coupled with the Latis pro software and the Eurosmart pendulum. We show that only a detailed analysis by fitting the recorded curve can provide sufficient accuracy to describe the quadratic evolution of the period as a function of the amplitude of the oscillations. We we can detect the additional harmonics in the oscillations when the amplitude becomes very high.


title: Cohomology and Deformations of n-Hom-Liebniz algebra morphisms

abstract: In this paper, we introduce cohomology of n-Hom-Liebniz algebra morphisms and formal deformation theory of n-Hom-Liebniz algebra morphisms .


title: Distortion and Faults in Machine Learning Software

abstract: Machine learning software, deep neural networks (DNN) software in particular, discerns valuable information from a large dataset, a set of data. Outcomes of such DNN programs are dependent on the quality of both learning programs and datasets. Unfortunately, the quality of datasets is difficult to be defined, because they are just samples. The quality assurance of DNN software is difficult, because resultant trained machine learning models are unknown prior to its development, and the validation is conducted indirectly in terms of prediction performance. This paper introduces a hypothesis that faults in the learning programs manifest themselves as distortions in trained machine learning models. Relative distortion degrees measured with appropriate observer functions may indicate that there are some hidden faults. The proposal is demonstrated with example cases of the MNIST dataset.


title: A Summary of Formal Specification and Verification of Autonomous Robotic Systems

abstract: Autonomous robotic systems are complex, hybrid, and often safety-critical; this makes their formal specification and verification uniquely challenging. Though commonly used, testing and simulation alone are insufficient to ensure the correctness of, or provide sufficient evidence for the certification of, autonomous robotics. Formal methods for autonomous robotics have received some attention in the literature, but no resource provides a current overview. This short paper summarises the contributions of Luckcuck 2019, which surveys the state-of-the-art in formal specification and verification for autonomous robotics.


title: On electron imaging of small molecules, Part II: tomographic reconstruction from defocus series

abstract: A practical method utilising three-dimensional image pattern matching is proposed which, in principle, is capable of unambiguous determination of the types and positions of atoms in small molecules from defocus series collected at only a few different angular orientations of the molecule. Numerical tests of the method are presented using multislice calculations of defocus series of small biological molecules. The proposed technique can in future gain from Bayesian or machine-learning approaches and is likely to be useful in cryogenic electron microscopy.


title: On electron imaging of small molecules, Part I: relative roles of multiple scattering and Fresnel diffraction

abstract: The relative roles of multiple electron scattering and in-molecule free-space propagation in transmission electron microscopy of small molecules are discussed. It is argued that while multiple scattering tends to have only a moderate effect in this case, the in-molecule Fresnel diffraction is likely to be significant due to the shallow depth of focus under the relevant experimental conditions. As a consequence, diffraction tomography based on the first Born or first Rytov approximation represents a more suitable method for the reconstruction of three-dimensional distribution of the electrostatic potential in this context, compared to conventional computed tomography which is intrinsically based on the projection approximation. A simplified method for diffraction tomography is proposed and tested on numerically simulated examples.


title: Fourier restriction above rectangles

abstract: In this article, we study the problem of obtaining Lebesgue space inequalities for the Fourier restriction operator associated to rectangular pieces of the paraboloid and perturbations thereof. We state a conjecture for the dependence of the operator norms in these inequalities on the sidelengths of the rectangles, prove that this conjecture follows from (a slight reformulation of the) restriction conjecture for elliptic hypersurfaces, and prove that, if valid, the conjecture is essentially sharp. Such questions arise naturally in the study of restriction inequalities for degenerate hypersurfaces; we demonstrate this connection by using our positive results to prove new restriction inequalities for a class of hypersurfaces having some additive structure.


title: Cross-lingual Multi-speaker Text-to-speech Synthesis for Voice Cloning without Using Parallel Corpus for Unseen Speakers

abstract: We investigate a novel cross-lingual multi-speaker text-to-speech synthesis approach for generating high-quality native or accented speech for native/foreign seen/unseen speakers in English and Mandarin. The system consists of three separately trained components: an x-vector speaker encoder, a Tacotron-based synthesizer and a WaveNet vocoder. It is conditioned on 3 kinds of embeddings: (1) speaker embedding so that the system can be trained with speech from many speakers will little data from each speaker; (2) language embedding with shared phoneme inputs; (3) stress and tone embedding which improves naturalness of synthesized speech, especially for a tonal language like Mandarin. By adjusting the various embeddings, MOS results show that our method can generate high-quality natural and intelligible native speech for native/foreign seen/unseen speakers. Intelligibility and naturalness of accented speech is low as expected. Speaker similarity is good for native speech from native speakers. Interestingly, speaker similarity is also good for accented speech from foreign speakers. We also find that normalizing speaker embedding x-vectors by L2-norm normalization or whitening improves output quality a lot in many cases, and the WaveNet performance seems to be language-independent: our WaveNet is trained with Cantonese speech and can be used to generate Mandarin and English speech very well.


title: Dissipation induced Luttinger liquid correlations in a one dimensional Fermi gas

abstract: We study a one-dimensional Fermi gas in the presence of dissipative coupling to environment through the Lindblad equation. The dissipation involves energy exchange with the environment and favours the relaxation of electrons to excitations. After switching on the dissipation, the system approaches a steady state, which is described by a generalized Gibbs ensemble. The fermionic single particle density matrix resembles deceivingly to that in a hermitian interaction quench. It decays inversely with the distance for short times due to the fermionic correlations in the initial state, which changes into a non-integer power law decay for late times, representing dissipation induced Luttinger liquid behaviour. However, the crossover between the two regions occurs due to dissipation induced damping, and is unrelated to the propagation of excitations. The velocity of information spreading is set by the dissipative coupling, and differs significantly from the original sound velocity. The thermodynamic entropy grows as $-t\ln t$ initially, and saturates to an extensive value. Our results can be tested experimentally in one-dimensional Dirac systems.


title: Eigenbranes in Jackiw-Teitelboim gravity

abstract: It was proven recently that JT gravity can be defined as an ensemble of L x L Hermitian matrices. We point out that the eigenvalues of the matrix correspond in JT gravity to FZZT-type boundaries on which spacetimes can end. We then investigate an ensemble of matrices with 1<<N<<L eigenvalues held fixed. This corresponds to a version of JT gravity which includes N FZZT type boundaries in the path integral contour and which is found to emulate a discrete quantum chaotic system. In particular this version of JT gravity can capture the behavior of finite-volume holographic correlators at late times, including erratic oscillations.


title: Unirational Differential Curves and Differential Rational Parametrizations

abstract: In this paper, we study unirational differential curves and the corresponding differential rational parametrizations. We first investigate basic properties of proper differential rational parametrizations for unirational differential curves. Then we show that the implicitization problem of proper linear differential rational parametric equations can be solved by means of differential resultants. Furthermore, for linear differential curves, we give an algorithm to determine whether an implicitly given linear differential curve is unirational and, in the affirmative case, to compute a proper differential rational parametrization for the differential curve.


title: Metric-affine bumblebee gravity: classical aspects

abstract: We consider the metric-affine formulation of bumblebee gravity, derive the field equations, and show that the connection can be written as Levi-Civita of a disformally related metric in which the bumblebee field determines the disformal part. As a consequence, the bumblebee field gets coupled to all the other matter fields present in the theory, potentially leading to nontrivial phenomenological effects. To explore this issue we compute the weak-field, post-Minkowskian limit and study the resulting effective theory. In this scenario, we couple scalar and spinorial matter to the effective metric, and then we explore the physical properties of the VEV of the bumblebee field, focusing mainly on the dispersion relations and the stability of the resulting effective theory.


title: The hidden topological structure of flow network functionality

abstract: The ability to reroute and control flow is vital to the function of venation networks across a wide range of organisms. By modifying individual edges in these networks, either by adjusting edge conductances or creating and destroying edges, organisms can robustly control the propagation of inputs to perform specific tasks. However, a fundamental disconnect exists between the structure and function of these networks: networks with different local architectures can perform the same functions. Here we answer the question of how structural changes at the microscopic level are able to collectively create functionality at the scale of an entire network. Using persistent homology, we analyze networks tuned to perform complex multifunctional tasks. We find that the responses of such networks encode a hidden topological structure composed of sectors of uniform pressure. Although these sectors are not apparent in the underlying network architectures, we find that they nonetheless correlate strongly with the tuned function. We conclude that the connectivity of these sectors, rather than that of the individual nodes, provides a quantitative relationship between structure and function in flow networks. Finally, we use this topological description to place a bound on the limits of task complexity.


title: Deep Learning with Gaussian Differential Privacy

abstract: Deep learning models are often trained on datasets that contain sensitive information such as individuals' shopping transactions, personal contacts, and medical records. An increasingly important line of work therefore has sought to train neural networks subject to privacy constraints that are specified by differential privacy or its divergence-based relaxations. These privacy definitions, however, have weaknesses in handling certain important primitives (composition and subsampling), thereby giving loose or complicated privacy analyses of training neural networks. In this paper, we consider a recently proposed privacy definition termed \textit{$f$-differential privacy} [18] for a refined privacy analysis of training neural networks. Leveraging the appealing properties of $f$-differential privacy in handling composition and subsampling, this paper derives analytically tractable expressions for the privacy guarantees of both stochastic gradient descent and Adam used in training deep neural networks, without the need of developing sophisticated techniques as [3] did. Our results demonstrate that the $f$-differential privacy framework allows for a new privacy analysis that improves on the prior analysis~[3], which in turn suggests tuning certain parameters of neural networks for a better prediction accuracy without violating the privacy budget. These theoretically derived improvements are confirmed by our experiments in a range of tasks in image classification, text classification, and recommender systems. Python code to calculate the privacy cost for these experiments is publicly available in the \texttt{TensorFlow Privacy} library.


title: Periodicity in fields of elongating dunes

abstract: Dune fields are commonly associated with periodic patterns that are among the most recognizable landscapes on Earth and other planetary bodies. However, in zones of limited sediment supply, where periodic dunes elongate and align in the direction of the resultant sand flux, there has been no attempt to explain the emergence of such a regular pattern. Here, we show, by means of numerical simulations, that the elongation growth mechanism does not produce a pattern with a specific wavelength. Periodic elongating dunes appear to be a juxtaposition of individual structures, the arrangement of which is due to regular landforms at the border of the field acting as boundary conditions. This includes, among others, dune patterns resulting from bed instability, or the crestline reorganization induced by dune migration. The wavelength selection in fields of elongating dunes therefore reflects the interdependence of dune patterns over the course of their evolution.


title: Poles and branch cuts in free surface hydrodynamics

abstract: We consider the motion of ideal incompressible fluid with free surface. We analyzed the exact fluid dynamics though the time-dependent conformal mapping $z=x+iy=z(w,t)$ of the lower complex half plane of the conformal variable $w$ into the area occupied by fluid. We established the exact results on the existence vs. nonexistence of the pole and power law branch point solutions for $1/z_w$ and the complex velocity. We also proved the nonexistence of the time-dependent rational solution of that problem for the second and the first order moving pole.


title: Improving EEG based Continuous Speech Recognition

abstract: In this paper we introduce various techniques to improve the performance of electroencephalography (EEG) features based continuous speech recognition (CSR) systems. A connectionist temporal classification (CTC) based automatic speech recognition (ASR) system was implemented for performing recognition. We introduce techniques to initialize the weights of the recurrent layers in the encoder of the CTC model with more meaningful weights rather than with random weights and we make use of an external language model to improve the beam search during decoding time. We finally study the problem of predicting articulatory features from EEG features in this paper.


title: Distributed Linear Quadratic Tracking Control for Leader-Follower Multi-Agent Systems: A Suboptimality Approach

abstract: In this paper, we extend the results from Jiao et al. (2019) on distributed linear quadratic control for leaderless multi-agent systems to the case of distributed linear quadratic tracking control for leader-follower multi-agent systems. Given one autonomous leader and a number of homogeneous followers, we introduce an associated global quadratic cost functional. We assume that the leader shares its state information with at least one of the followers and the communication between the followers is represented by a connected simple undirected graph. Our objective is to design distributed control laws such that the controlled network reaches tracking consensus and, moreover, the associated cost is smaller than a given tolerance for all initial states bounded in norm by a given radius. We establish a centralized design method for computing such suboptimal control laws, involving the solution of a single Riccati inequality of dimension equal to the dimension of the local agent dynamics, and the smallest and the largest eigenvalue of a given positive definite matrix involving the underlying graph. The proposed design method is illustrated by a simulation example.


title: On Symbiosis of Attribute Prediction and Semantic Segmentation

abstract: In this paper, we propose to employ semantic segmentation to improve person-related attribute prediction. The core idea lies in the fact that the probability of an attribute to appear in an image is far from being uniform in the spatial domain. We build our attribute prediction model jointly with a deep semantic segmentation network. This harnesses the localization cues learned by the semantic segmentation to guide the attention of the attribute prediction to the regions where different attributes naturally show up. Therefore, in addition to prediction, we are able to localize the attributes despite merely having access to image-level labels (weak supervision) during training. We first propose semantic segmentation-based pooling and gating, respectively denoted as SSP and SSG. In the former, the estimated segmentation masks are used to pool the final activations of the attribute prediction network, from multiple semantically homogeneous regions. In SSG, the same idea is applied to the intermediate layers of the network. SSP and SSG, while effective, impose heavy memory utilization since each channel of the activations is pooled/gated with all the semantic segmentation masks. To circumvent this, we propose Symbiotic Augmentation (SA), where we learn only one mask per activation channel. SA allows the model to either pick one, or combine (weighted superposition) multiple semantic maps, in order to generate the proper mask for each channel. SA simultaneously applies the same mechanism to the reverse problem by leveraging output logits of attribute prediction to guide the semantic segmentation task. We evaluate our proposed methods for facial attributes on CelebA and LFWA datasets, while benchmarking WIDER Attribute and Berkeley Attributes of People for whole body attributes. Our proposed methods achieve superior results compared to the previous works.


title: On some geometric properties of normalized Wright functions

abstract: The main purpose of the present paper is to determine the radii of starlikeness and convexity associated with lemniscate of Bernoulli and the Janowski function, $(1+Az)/(1+Bz)$ for $-1\leq B<A\leq 1,$ of normalized Wright functions. The key tools in the proof of our main results are the infinite product representation of Wright function and properties of real zeros of the Wright function and its derivative.


title: Dirac star with SU(2) Yang-Mills and Proca fields

abstract: We study spherically symmetric strongly gravitating configurations supported by nonlinear spinor fields and non-Abelian SU(2) Yang-Mills/Proca magnetic fields. Regular asymptotically flat solutions describing objects with positive Arnowitt-Deser-Misner masses are obtained numerically. When the mass of the spinor fields is much smaller than the Planck mass, we find approximate solutions that can describe systems with total masses comparable to the Chandrasekhar mass and with effective radii of the order of kilometers. For the values of the system free parameters used here, we show that the SU(2) magnetic field always gives a small contribution to the total energy density and mass of the configurations under investigation. From the astrophysical point of view, one can regard such objects as magnetized Dirac stars.


title: Low-$n$ global ideal MHD instabilities in CFETR baseline scenario

abstract: This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. A good agreement is achieved between two codes in the growth rates of $n=1-10$ ideal MHD modes for various locations of the perfect conducting first-wall. The higher-$n$ modes are dominated by ballooning modes and localized in the pedestal region, while the lower-$n$ modes have more prominent external kink components and broader mode profiles. The influences of plasma-vacuum profile and wall shape are also examined using NIMROD. In presence of resistive wall, the low-$n$ ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the $n = 1$ resistive wall mode (RWM) is found unstable, which could be fully stabilized by uniform toroidal rotation above 2.9\% core Alfv\'en speed.


title: Enhancing Cross-task Black-Box Transferability of Adversarial Examples with Dispersion Reduction

abstract: Neural networks are known to be vulnerable to carefully crafted adversarial examples, and these malicious samples often transfer, i.e., they remain adversarial even against other models. Although great efforts have been delved into the transferability across models, surprisingly, less attention has been paid to the cross-task transferability, which represents the real-world cybercriminal's situation, where an ensemble of different defense/detection mechanisms need to be evaded all at once. In this paper, we investigate the transferability of adversarial examples across a wide range of real-world computer vision tasks, including image classification, object detection, semantic segmentation, explicit content detection, and text detection. Our proposed attack minimizes the ``dispersion'' of the internal feature map, which overcomes existing attacks' limitation of requiring task-specific loss functions and/or probing a target model. We conduct evaluation on open source detection and segmentation models as well as four different computer vision tasks provided by Google Cloud Vision (GCV) APIs, to show how our approach outperforms existing attacks by degrading performance of multiple CV tasks by a large margin with only modest perturbations linf=16.


title: On topological Rudin's lemma, well-filtered spaces and sober spaces

abstract: Based on topological Rudin's Lemma, we investigate two new kinds of sets - Rudin sets and well-filtered determined sets in $T_0$ topological spaces. Using such sets, we formulate and prove some new characterizations for well-filtered spaces and sober spaces. Part of the work was inspired by Xi and Lawson's work on well-filtered spaces. Our study also lead to a new class of spaces - strong $d$-spaces and some problems whose solutions will strengthen our understanding of the related structures.


title: A direct approach to $K$-reflections of $T_0$ spaces

abstract: In this paper, we provide a direct approach to $\mathbf{K}$-reflections of $T_0$ spaces. For a full subcategory $\mathbf{K}$ of the category of all $T_0$ spaces and a $T_0$ space $X$, let $\mathbf{K}(X)=\{A\subseteq X : A$ is closed and for any continuous mapping $f : X\longrightarrow Y$ to a $\mathbf{K}$-space $Y$, there exists a unique $y_A\in Y$ such that $\overline{f(A)}=\overline{\{y_A\}}\}$ and $P_H(\mathbf{K}(X))$ the space of $\mathbf{K}(X)$ endowed with the lower Vietoris topology. It is proved that if $P_H(\mathbf{K}(X))$ is a $\mathbf{K}$-space, then the pair $\langle X^k=P_H(\mathbf{K}(X)), \eta_X\rangle$, where $\eta_X :X\longrightarrow X^k$, $x\mapsto\overline{\{x\}}$, is the $\mathbf{K}$-reflection of $X$. We call $\mathbf{K}$ an adequate category if for any $T_0$ space $X$, $P_H(\mathbf{K}(X))$ is a $\mathbf{K}$-space. Therefore, if $\mathbf{K}$ is adequate, then $\mathbf{K}$ is reflective in $\mathbf{Top}_0$. It is shown that the category of all sober spaces, that of all $d$-spaces, that of all well-filtered spaces and the Keimel and Lawson's category are all adequate, and hence are all reflective in $\mathbf{Top}_0$. Some major properties of $\mathbf{K}$-spaces and $\mathbf{K}$-reflections of $T_0$ spaces are investigated. In particular, it is proved that if $\mathbf{K}$ is adequate, then the $\mathbf{K}$-reflection preserves finite products of $T_0$ spaces. Our study also leads to a number of problems, whose answering will deepen our understanding of the related spaces and their categorical structures.


title: Joint Spatial and Angular Super-Resolution from a Single Image

abstract: Synthesizing a densely sampled light field from a single image is highly beneficial for many applications. Moreover, jointly solving both angular and spatial super-resolution problem also introduces new possibilities in light field imaging. The conventional method relies on physical-based rendering and a secondary network to solve the angular super-resolution problem. In addition, pixel-based loss limits the network capability to infer scene geometry globally. In this paper, we show that both super-resolution problems can be solved jointly from a single image by proposing a single end-to-end deep neural network that does not require a physical-based approach. Two novel loss functions based on known light field domain knowledge are proposed to enable the network to preserve the spatio-angular consistency between sub-aperture images. Experimental results show that the proposed model successfully synthesizes dense high resolution light field and it outperforms the state-of-the-art method in both quantitative and qualitative criteria. The method can be generalized to arbitrary scenes, rather than focusing on a particular subject. The synthesized light field can be used for various applications, such as depth estimation and refocusing.


title: Teaching Perception

abstract: The visual world is very rich and generally too complex to perceive in its entirety. Yet only certain features are typically required to adequately perform some task in a given situation. Rather than hardwire-in decisions about when and what to sense, this paper describes a robotic system whose behavioral policy can be set by verbal instructions it receives. These capabilities are demonstrated in an associated video showing the fully implemented system guiding the perception of a physical robot in simple scenario. The structure and functioning of the underlying natural language based symbolic reasoning system is also discussed.


title: On quantumness in multi-parameter quantum estimation

abstract: In this article we derive a measure of quantumness in quantum multi-parameter estimation problems. We can show that the ratio between the mean Uhlmann Curvature and the Fisher Information provides a figure of merit which estimates the amount of incompatibility arising from the quantum nature of the underlying physical system. This ratio accounts for the discrepancy between the attainable precision in the simultaneous estimation of multiple parameters and the precision predicted by the Cram\'er-Rao bound. As a testbed for this concept, we consider a quantum many-body system in thermal equilibrium, and explore the quantum compatibility of the model across its phase diagram.


title: A discriminative condition-aware backend for speaker verification

abstract: We present a scoring approach for speaker verification that mimics the standard PLDA-based backend process used in most current speaker verification systems. However, unlike the standard backends, all parameters of the model are jointly trained to optimize the binary cross-entropy for the speaker verification task. We further integrate the calibration stage inside the model, making the parameters of this stage depend on metadata vectors that represent the conditions of the signals. We show that the proposed backend has excellent out-of-the-box calibration performance on most of our test sets, making it an ideal approach for cases in which the test conditions are not known and development data is not available for training a domain-specific calibration model.


title: A Vietnamese information retrieval system for product-price

abstract: A price information retrieval (IR) system allows users to search and view differences among prices of specific products. Building product-price driven IR system is a challenging and active research area. Approaches entirely depending products information provided by shops via interface environment encounter limitations of database. While automatic systems specifically require product names and commercial websites for their input. For both paradigms, approaches of building product-price IR system for Vietnamese are still very limited. In this paper, we introduce an automatic Vietnamese IR system for product-price by identifying and storing Xpath patterns to extract prices of products from commercial websites. Experiments of our system show promising results.


title: An analysis of the graph processing landscape

abstract: The value of graph-based big data can be unlocked by exploring the topology and metrics of the networks they represent, and the computational approaches to this exploration take on many forms. The use-case of performing global computations over a graph, it is first ingested into a graph processing system from one of many digital representations. Extracting information from graphs involves processing all their elements globally, and can be done with single-machine systems (with varying approaches to hardware usage), distributed systems (either homogeneous or heterogeneous groups of machines) and systems dedicated to high-performance computing (HPC). We provide an overview of different aspects of the graph processing landscape and describe classes of systems based on a set of dimensions we detail. The dimensions we detail encompass paradigms to express graph processing, different types of systems to use, coordination and communication models in distributed graph processing, partitioning techniques and different definitions related to the potential for a graph to be updated. This survey is aimed at both the experienced software engineer or researcher as well as the newcomer looking for an understanding of the landscape of solutions (and their limitations) for graph processing.


title: Invariants of relatively generic surface singularities II. Images of Abel maps

abstract: In \cite{R} the author investigated invariants of relatively generic structures on surface singularities generalising results of \cite{NNA1} and \cite{NNA2} about generic analytic structures and generic line bundles to the case of the relative setup, where we fix a given analytic type or line bundle on a smaller subgraph or more generally on a smaller cycle and we choose a relatively generic line bundle or analytic type on the large cycle and managed to compute some of it's invariants, like geometric genus or $h^1$ of natural line bundles. In \cite{NNAD} the authors investigated the images of Abel maps $c^{l'}(Z) : \eca^{l'}(Z) \to \pic^{l'}(Z)$, where $l' \in - S'(|Z|)$, especially the dimensions of the images of these maps and gave two algorithms to compute these invariants from cohomology numbers of cycles and from periodic constants of singularities we get from $\tX$ by blowing it up at generic points sequentially. Furthemore in \cite{NNAD} the authors gave explicit combinatorial formulas in the case of generic singularities. In this paper we want to generalise the theorems from \cite{NNAD} to the relatively generic case. In this case we fix a subsingularity $\tX_1$ for a subgraph $\mathcal{T}_1 \subset \mathcal{T}$ and a relatively generic singularity $\tX$ corresponding to $\tX_1$. Furthermore we fix a line bundle $\calL$ on $\tX_1$ and a Chern class $l' \in - S'$, such that $c^1(\calL) = R(l')$. Our main goal in the article is to compute $\dim(c^{l'}(\eca^{ l', \calL}(Z)))$ from invariants of the subsingularity $\tX_1$ and to conclude a few corollaries.


title: Construction of the classical time crystal Lagrangians from Sisyphus dynamics and duality description with the Li\'enard type equation

abstract: We explore the connection between the equations describing Sisyphus dynamics and the generic Li\'{e}nard type or Li\'{e}nard equation from the viewpoint of branched Hamiltonians. The former provides the appropriate setting for classical time crystal being derivable from a higher order Lagrangian. However it appears the equations of Sisyphus dynamics have a close relation with the Li\'{e}nard-II equation when expressed in terms of the `velocity' variable. Another interesting feature of the equations of Sisyphus dynamics is the appearance of velocity dependent "mass function" in contrast to the more commonly encountered position dependent mass. The consequences of such mass functions seem to have connections to cosmological time crystals .


title: Condensation phenomena in gravity

abstract: Gravity can play a role in critical phenomena. Topological singularities induce ground state degeneracy and break the continuum symmetry of the vacuum. They also generate momenta oscillations about an average momentum and a positive gravitational susceptibility. Gravitational analogues of the laws of Curie and Bloch have been found for a one-dimensional model. The critical temperature for a change in phase from unbound to isolated vortices has also been calculated in a $XY$-model.


title: A degenerate elliptic equation for second order controllability of nonlinear systems

abstract: For a general nonlinear control system, we study the problem of small time local attainability of a target which is the closure of an open set. When the target is smooth and locally the sublevel set of a smooth function, we develop second order attainability conditions as explicit pointwise conditions on the vector fields at points where all the available vector fields are contained in the tangent space of its boundary. Our sufficient condition requires the function defining the target to be a strict supersolution of a second order degenerate elliptic equation and if satisfied, it allows to reach the target with a piecewise constant control with at most one switch. For symmetric systems, our sufficient condition is also necessary and can be reformulated as a suitable symmetric matrix having a negative eigenvalue. For nonlinear affine systems to obtain a necessary and sufficient condition we require an additional request on the drift. Our second order pde has the same role of the Hamilton-Jacobi equation for first order sufficient conditions.


title: Logical Interpretations of Autoencoders

abstract: The unification of low-level perception and high-level reasoning is a long-standing problem in artificial intelligence, which has the potential to not only bring the areas of logic and learning closer together but also demonstrate how abstract concepts might emerge from sensory data. Precisely because deep learning methods dominate perception-based learning, including vision, speech, and linguistic grammar, there is fast-growing literature on how to integrate symbolic reasoning and deep learning. Broadly, efforts seem to fall into three camps: those focused on defining a logic whose formulas capture deep learning, ones that integrate symbolic constraints in deep learning, and others that allow neural computations and symbolic reasoning to co-exist separately, to enjoy the strengths of both worlds. In this paper, we identify another dimension to this inquiry: what do the hidden layers really capture, and how can we reason about that logically? In particular, we consider autoencoders that are widely used for dimensionality reduction and inject a symbolic generative framework onto the feature layer. This allows us, among other things, to generate example images for a class to get a sense of what was learned. Moreover, the modular structure of the proposed model makes it possible to learn relations over multiple images at a time, as well as handle noisy labels. Our empirical evaluations show the promise of this inquiry.


title: The three-loop polarized pure singlet operator matrix element with two different masses

abstract: We present the two-mass QCD contributions to the polarized pure singlet operator matrix element at three loop order in $x$-space. These terms are relevant for calculating the polarized structure function $g_1(x,Q^2)$ at $O(\alpha_s^3)$ as well as for the matching relations in the variable flavor number scheme and the polarized heavy quark distribution functions at the same order. The result for the operator matrix element is given in terms of generalized iterated integrals. These integrals depend on the mass ratio through the main argument, and the alphabet includes square--root valued letters.


title: Bridging the Gap between Semantics and Multimedia Processing

abstract: In this paper, we give an overview of the semantic gap problem in multimedia and discuss how machine learning and symbolic AI can be combined to narrow this gap. We describe the gap in terms of a classical architecture for multimedia processing and discuss a structured approach to bridge it. This approach combines machine learning (for mapping signals to objects) and symbolic AI (for linking objects to meanings). Our main goal is to raise awareness and discuss the challenges involved in this structured approach to multimedia understanding, especially in the view of the latest developments in machine learning and symbolic AI.


title: Minimal Linear Codes Constructed from Functions

abstract: In this paper, we consider minimal linear codes in a general construction of linear codes from q-ary functions. First, we give the sufficient and necessary condition for codewords to be minimal. Second, as an application, we present four constructions of minimal linear codes which contained some recent results as special cases.


title: Assessing Impacts of Abnormal Events on Travel Patterns Leveraging Passively Collected Trajectory Data

abstract: Travel patterns can be impacted by abnormal events. Assessing the impacts has important implications for relief operations and improving preparedness or planning for future events. Conventionally, the assessment is done followed by data collection from post-event surveys, which are economically costly, suffering low-response rate, time-consuming and usually delayed for months (or even years) after an event, leading to inefficient and unreliable assessment and creating obstacles for relief organizations to reach people in need. Penetration of smartphones and services enabled by them continuously generate large amount of trajectory data (e.g., Call Records Data, App-based data), containing trajectories of massive users. These trajectory data are passively and timely collected and without additional cost and contain information of travel patterns of the massive number of individuals in a region for a prolonged time period (e.g., months to years). We propose a framework to assessing the impacts on travel patterns using these data. Utilizing the passively collected trajectory data, the proposed framework seeks to capturing and understanding the full spectrum of travel pattern changes, which helps to assess who, when and how people in a certain area were impacted. The proposed framework is applied to a mobile phone trajectory dataset containing about half-year trajectories of a million anonymous users to assess the impacts of Hurricane Harvey (the second-costliest hurricane in US history). The results are validated and show that the proposed framework can provide a comprehensive assessment of impacts of Harvey on travel patterns, which could guide the response to and the recovery from the impacts.


title: GBCNs: Genetic Binary Convolutional Networks for Enhancing the Performance of 1-bit DCNNs

abstract: Training 1-bit deep convolutional neural networks (DCNNs) is one of the most challenging problems in computer vision, because it is much easier to get trapped into local minima than conventional DCNNs. The reason lies in that the binarized kernels and activations of 1-bit DCNNs cause a significant accuracy loss and training inefficiency. To address this problem, we propose Genetic Binary Convolutional Networks (GBCNs) to optimize 1-bit DCNNs, by introducing a new balanced Genetic Algorithm (BGA) to improve the representational ability in an end-to-end framework. The BGA method is proposed to modify the binary process of GBCNs to alleviate the local minima problem, which can significantly improve the performance of 1-bit DCNNs. We develop a new BGA module that is generic and flexible, and can be easily incorporated into existing DCNNs, such asWideResNets and ResNets. Extensive experiments on the object classification tasks (CIFAR, ImageNet) validate the effectiveness of the proposed method. To highlight, our method shows strong generalization on the object recognition task, i.e., face recognition, facial and person re-identification.


title: The study of the gluon distribution function and reduced cross section behavior using the proton structure function

abstract: The behavior of the gluon distribution function and the reduced cross section considered from the proton structure function and its derivatives at low values of $x$. These behaviors studied and compared with the experimental data. These results are augmented by including an additional higher-twist term in the description of the nonlinear correction. This additional term, modified nonlinear correction, improves the description of the reduced cross sections significantly at low values of $Q^{2}$. We discuss, furthermore, how this behavior can be determine the reduced cross section with respect to the proton parameterization at high-y values. The resulting predictions for $\sigma_{r}$ suggest that further corrections are required for $Q^{2}$ less than about $3~\mathrm{GeV}^{2}$.


title: Solving Traveltime Tomography with Deep Learning

abstract: This paper introduces a neural network approach for solving two-dimensional traveltime tomography (TT) problems based on the eikonal equation. The mathematical problem of TT is to recover the slowness field of a medium based on the boundary measurement of the traveltimes of waves going through the medium. This inverse map is high-dimensional and nonlinear. For the circular tomography geometry, a perturbative analysis shows that the forward map can be approximated by a vectorized convolution operator in the angular direction. Motivated by this and filtered back-projection, we propose an effective neural network architecture for the inverse map using the recently proposed BCR-Net, with weights learned from training datasets. Numerical results demonstrate the efficiency of the proposed neural networks.


title: Fast Fibonacci heaps with worst case extensions

abstract: We are concentrating on reducing overhead of heaps based on comparisons with optimal worstcase behaviour. The paper is inspired by Strict Fibonacci Heaps [1], where G. S. Brodal, G. Lagogiannis, and R. E. Tarjan implemented the heap with DecreaseKey and Meld interface in assymptotically optimal worst case times (based on key comparisons). In the paper [2], the ideas were elaborated and it was shown that the same asymptotical times could be achieved with a strategy loosing much less information from previous comparisons. There is big overhead with maintainance of violation lists in these heaps. We propose simple alternative reducing this overhead. It allows us to implement fast amortized Fibonacci heaps, where user could call some methods in variants guaranting worst case time. If he does so, the heaps are not guaranted to be Fibonacci until an amortized version of a method is called. Of course we could call worst case versions all the time, but as there is an overhead with the guarantee, calling amortized versions is prefered choice if we are not concentrated on complexity of the separate operation. We have shown, we could implement full DecreaseKey-Meld interface, but Meld interface is not natural for these heaps, so if Meld is not needed, much simpler implementation suffices. As I don't know application requiring Meld, we would concentrate on noMeld variant, but we will show the changes could be applied on Meld including variant as well. The papers [1], [2] shown the heaps could be implemented on pointer machine model. For fast practical implementations we would rather use arrays. Our goal is to reduce number of pointer manipulations. Maintainance of ranks by pointers to rank lists would be unnecessary overhead.


title: Resolving isotope splitting of boron-related intracenter transitions in diamond by infrared absorption spectroscopy

abstract: Isotopic enrichment offers cutting-edge properties of materials opening exciting research and development opportunities. In semiconductors, reached progress of ultimate control in growth and doping techniques follows nowadays the high level isotopic purification. This requires deep understanding of isotopic disorder effects and techniques of their effective determination. Isotopic content of both crystal lattice and impurity centers cause the effects, which can be examined by different optical techniques. While disorder in the host lattice can be straight forward evaluated by inelastic light scattering or by SIMS measurements, determination of isotopic contributions of many orders less presented impurities remains challenging and usually observed in high-resolution photoluminescence or optical absorption spectra. Boron-doped diamonds exhibit complex infrared absorption spectra while boron-related luminescence remains unobserved. Boron, as a most light element acting as a hydrogen-like dopant in elemental semiconductors, has a largest relative difference in its isotope masses, and by this, cause the largest isotopic disorder in semiconductors, including diamond, an elemental semiconductor with the lightest atomic mass of a host lattice. This enables an access to the isotopic constitution of boron in diamond by infrared absorption spectroscopy. By comparison of low-temperature absorption spectra of a natural (20% of 10B and 80% of 11B isotopes) and 11B enriched (up to 99%) doped diamonds we differentiate the intracenter transitions related to 10B and to 11B isotopes. We have found that the isotopic spectral lines of the same boron intracenter transition are separated with the energy of about 0.7 meV. This is the largest impurity isotopic shift ever observed in semiconductors doped by hydrogen-like impurity centers.


title: Magneto-conductance of topological junctions based on two-dimensional electron gases reveals Majorana phases

abstract: We calculate the linear conductance of a two-dimensional electron gas (2DEG)-based junction between a normal semiconductor section and a hybrid semiconductor-superconductor section, under perpendicular magnetic field. We consider two important terms often neglected in the literature, the magneto-orbital and transverse Rashba spin-orbit. The strong orbital effect due to the magnetic field yields topological phase transitions to nontrivial phases hosting Majorana modes in the hybrid section. The presence of a potential barrier at the junction interface reveals the Majorana phases as quantized plateaus of high conductance, for low values of the chemical potential. In wide junctions (or large chemical potentials) the phase transitions occur at low magnetic fields but the magneto-conductance becomes anomalous and lacks clearly quantized plateaus.


title: A Stochastic Trust-Region Framework for Policy Optimization

abstract: In this paper, we study a few challenging theoretical and numerical issues on the well known trust region policy optimization for deep reinforcement learning. The goal is to find a policy that maximizes the total expected reward when the agent acts according to the policy. The trust region subproblem is constructed with a surrogate function coherent to the total expected reward and a general distance constraint around the latest policy. We solve the subproblem using a preconditioned stochastic gradient method with a line search scheme to ensure that each step promotes the model function and stays in the trust region. To overcome the bias caused by sampling to the function estimations under the random settings, we add the empirical standard deviation of the total expected reward to the predicted increase in a ratio in order to update the trust region radius and decide whether the trial point is accepted. Moreover, for a Gaussian policy which is commonly used for continuous action space, the maximization with respect to the mean and covariance is performed separately to control the entropy loss. Our theoretical analysis shows that the deterministic version of the proposed algorithm tends to generate a monotonic improvement of the total expected reward and the global convergence is guaranteed under moderate assumptions. Comparisons with the state-of-the-art methods demonstrate the effectiveness and robustness of our method over robotic controls and game playings from OpenAI Gym.


title: PIQA: Reasoning about Physical Commonsense in Natural Language

abstract: To apply eyeshadow without a brush, should I use a cotton swab or a toothpick? Questions requiring this kind of physical commonsense pose a challenge to today's natural language understanding systems. While recent pretrained models (such as BERT) have made progress on question answering over more abstract domains - such as news articles and encyclopedia entries, where text is plentiful - in more physical domains, text is inherently limited due to reporting bias. Can AI systems learn to reliably answer physical common-sense questions without experiencing the physical world? In this paper, we introduce the task of physical commonsense reasoning and a corresponding benchmark dataset Physical Interaction: Question Answering or PIQA. Though humans find the dataset easy (95% accuracy), large pretrained models struggle (77%). We provide analysis about the dimensions of knowledge that existing models lack, which offers significant opportunities for future research.


title: System Performance with varying L1 Instruction and Data Cache Sizes: An Empirical Analysis

abstract: In this project, we investigate the fluctuations in performance caused by changing the Instruction (I-cache) size and the Data (D-cache) size in the L1 cache. We employ the Gem5 framework to simulate a system with varying specifications on a single host machine. We utilize the FreqMine benchmark available under the PARSEC suite as the workload program to benchmark our simulated system. The Out-order CPU (O3) with Ruby memory model was simulated in a Full-System X86 environment with Linux OS. The chosen metrics deal with Hit Rate, Misses, Memory Latency, Instruction Rate, and Bus Traffic within the system. Performance observed by varying L1 size within a certain range of values was used to compute Confidence Interval based statistics for relevant metrics. Our expectations, corresponding experimental observations, and discrepancies are also discussed in this report.


title: Polynomial Trace Identities in $SL(2,{\bf C})$, Quaternion Algebras, and Two-generator Kleinian Groups

abstract: We study certain polynomial trace identities in the group $SL(2,\IC)$ and their application in the theory of discrete groups. We obtain canonical representations for two generator groups in \S 4 and then in \S 5 we give a new proof for Gehring and Martin's polynomial trace identities for good words, and extend that result to a larger class which is also closed under a semigroup operation inducing polynomial composition. This new approach is through the use of quaternion algebras over indefinites and an associated group of units. We obtain structure theorems for these quaternion algebras which appear to be of independent interest in \S 8. Using these quaternion algebras and their units, we consider their relation to discrete subgroups of $SL(2,\IC)$ giving necessary and sufficient criteria for discreteness, and another for arithmeticity \S 9. We then show that for the groups $\IZ_p*\IZ_2$, the complement of the closure of roots of the good word polynomials is precisely the moduli space of geometrically finite discrete and faithful representations a result we show holds in greater generality in \S12.


title: Slow light,double directional rainbow trapping and releasing at terahertz frequencies

abstract: Slow light and rainbow trapping attract many attentions in last twenty years, and in most of the previous works, the researchers achieve the slow light and rainbow trapping with complicate configurations or techniques, for example, metamaterial techniques. In this paper, we propose a simple waveguide consisted with semiconductor and dielectric layers, and we find that the slow-light peaks appear in the dispersion curves. Besides, the cutoff frequency with $v_g = 0$ ($v_g$ is the group velocity) in the slow-light peaks is depending on the thicknesses of the semiconductor and dielectric. Then, we design a tapered, horizontally symmetric structure. By using software COMSOL and finite difference time domain method as well, we achieve double directional rainbow trapping and more importantly, releasing the rainbow in such structure. This approach is promising for optical isolator, optical buffer, optical switch and other optical functional devices in miniaturization optical integrated circuit.


title: Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains

abstract: By discretising space into compartments and letting system dynamics be governed by the reaction-diffusion master equation, it is possible to derive and simulate a stochastic model of reaction and diffusion on an arbitrary domain. However, there are many implementation choices involved in this process, such as the choice of discretisation and method of derivation of the diffusive jump rates, and it is not clear a priori how these affect model predictions. To shed light on this issue, in this work we explore how a variety of discretisations and method for derivation of the diffusive jump rates affect the outputs of stochastic simulations of reaction-diffusion models, in particular using Turing's model of pattern formation as a key example. We consider both static and uniformly growing domains and demonstrate that, while only minor differences are observed for simple reaction-diffusion systems, there can be vast differences in model predictions for systems that include complicated reaction kinetics, such as Turing's model of pattern formation. Our work highlights that care must be taken in using the reaction-diffusion master equation to make predictions as to the dynamics of stochastic reaction-diffusion systems.


title: A High-dimensional M-estimator Framework for Bi-level Variable Selection

abstract: In high-dimensional data analysis, bi-level sparsity is often assumed when covariates function group-wisely and sparsity can appear either at the group level or within certain groups. In such cases, an ideal model should be able to encourage the bi-level variable selection consistently. Bi-level variable selection has become even more challenging when data have heavy-tailed distribution or outliers exist in random errors and covariates. In this paper, we study a framework of high-dimensional M-estimation for bi-level variable selection. This framework encourages bi-level sparsity through a computationally efficient two-stage procedure. In theory, we provide sufficient conditions under which our two-stage penalized M-estimator possesses simultaneous local estimation consistency and the bi-level variable selection consistency if certain nonconvex penalty functions are used at the group level. Both our simulation studies and real data analysis demonstrate satisfactory finite sample performance of the proposed estimators under different irregular settings.


title: ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

abstract: To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.


title: On grouops with formational subnormal or self-normalizing subgroups

abstract: We establish the structure of finite groups with $\mathfrak{F}$-subnormal or self-normalizing primary cyclic subgroups in case $\mathfrak{F}$ is a subgroup-closed saturated superradical formation containing all nilpotent groups.


title: BPS states, conserved charges and centres of symmetric group algebras

abstract: In $\mathcal{N}=4$ SYM with $U(N)$ gauge symmetry, the multiplicity of half-BPS states with fixed dimension can be labelled by Young diagrams and can be distinguished using conserved charges corresponding to Casimirs of $U(N)$. The information theoretic study of LLM geometries and superstars in the dual $AdS_5 \times S^5$ background has raised a number of questions about the distinguishability of Young diagrams when a finite set of Casimirs are known. Using Schur-Weyl duality relations between unitary groups and symmetric groups, these questions translate into structural questions about the centres of symmetric group algebras. We obtain analytic and computational results about these structural properties and related Shannon entropies, and generate associated number sequences. A characterization of Young diagrams in terms of content distribution functions relates these number sequences to diophantine equations. These content distribution functions can be visualized as connected, segmented, open strings in content space.


title: A continuation method in Bayesian inference

abstract: We present a continuation method that entails generating a sequence of transition probability density functions from the prior to the posterior in the context of Bayesian inference for parameter estimation problems. The characterization of transition distributions, by tempering the likelihood function, results in a homogeneous nonlinear partial integro-differential equation whose existence and uniqueness of solutions are addressed. The posterior probability distribution comes as the interpretation of the final state of the path of transition distributions. A computationally stable scaling domain for the likelihood is explored for the approximation of the expected deviance, where we manage to hold back all the evaluations of the forward predictive model at the prior stage. It follows the computational tractability of the posterior distribution and opens access to the posterior distribution for direct samplings. To get a solution formulation of the expected deviance, we derive a partial differential equation governing the moments generating function of the log-likelihood. We show also that a spectral formulation of the expected deviance can be obtained for low-dimensional problems under certain conditions. The computational efficiency of the proposed method is demonstrated through three differents numerical examples that focus on analyzing the computational bias generated by the method, assessing the continuation method in the Bayesian inference with non-Gaussian noise, and evaluating its ability to invert a multimodal parameter of interest.


title: Algorithm for automated tuning of a quantum dot into the single-electron regime

abstract: We report an algorithm designed to perform computer-automated tuning of a single quantum dot with a charge sensor. The algorithm performs an adaptive measurement sequence of sub-sized stability diagrams until the single-electron regime is identified and reached. For each measurement, the signal processing module removes the physical background of the charge sensor to generate a binary image of charge transitions. Then, the image analysis module identifies the position and number of lines using two line detection schemes that are robust to noise and missing data.


title: Assessing Supply Chain Cyber Risks

abstract: Risk assessment is a major challenge for supply chain managers, as it potentially affects business factors such as service costs, supplier competition and customer expectations. The increasing interconnectivity between organisations has put into focus methods for supply chain cyber risk management. We introduce a general approach to support such activity taking into account various techniques of attacking an organisation and its suppliers, as well as the impacts of such attacks. Since data is lacking in many respects, we use structured expert judgment methods to facilitate its implementation. We couple a family of forecasting models to enrich risk monitoring. The approach may be used to set up risk alarms, negotiate service level agreements, rank suppliers and identify insurance needs, among other management possibilities.


title: Mapping and Monitoring Pollution Levels of Carbon Monoxide (CO) using Arduino and Location-Based Service

abstract: Air pollution is a serious issue in the world. Around 98% of cities with a population of over 100,000 people in low and middle-income countries do not meet air quality standards, while in high-income countries, the number has decreased by 52%. Industries and motor vehicles are the biggest contributors to carbon monoxide (CO).Air pollution is a severe issue in the world. The object of research to investigate is the detection of carbon monoxide (CO) levels in an area mapped based on the current CO levels along with the location coordinates via GPS and LBS technology with a microcontroller and sensor-based device.The results of testing the CO level detection devices found that high CO levels are in the afternoon with an average CO level of 49.59656 which means classified as dangerous if we are outdoors more than 30 minutes, can interfere with heart function.


title: Dynamic nuclear polarization and ESR hole burning in As doped silicon

abstract: We present an experimental study of the Dynamic Nuclear Polarization (DNP) of \si{} nuclei in silicon crystals of natural abundance doped with As in the temperature range 0.1-1 K and in strong magnetic field of 4.6 T. This ensures very high degree of electron spin polarization, extremely slow nuclear relaxation and optimal conditions for realization of Overhauser and resolved solid effects. We found that the solid effect DNP leads to an appearance of a pattern of holes and peaks in the ESR line, separated by the super-hyperfine interaction between the donor electron and \si{} nuclei closest to the donor. On the contrary, the Overhauser effect DNP mainly affects the remote \si{} nuclei having the weakest interaction with the donor electron. This leads to an appearance of a very narrow ($\approx$ 3 mG wide) hole in the ESR line. We studied relaxation of the holes after burning, which is caused by the nuclear spin diffusion. Analyzing the spin diffusion data with a simple one-dimensional spectral diffusion model leads to a value of the spectral diffusion coefficient $D=8(3)\times 10^{-3}$ mG$^2$/s. Our data indicate that the spin diffusion is not completely prevented even in the frozen core near the donors. The emergence of the narrow hole after the Overhauser DNP may be explained by a partial "softening" of the frozen core caused by Rabi oscillations of the electron spin.


title: Magnetic domain wall motion in SrRuO$_3$ thin films

abstract: Influence of substrate miscut on magnetization dynamics in SrRuO$_3$ (SRO) thin films was studied. Two films were grown on SrTiO$_3$ substrates with high ($\sim1^{\circ}$) and low ($\sim0.1^{\circ}$) miscut angles, respectively. As expected, high miscut angle leads to suppression of multi-variant growth. By means of SQUID magnetometry, comparable relaxation effects were observed in both the multi-variant and the nearly single-variant sample. Differences in the magnetization reversal process were revealed by magnetic force microscopy. It showed that the multi-variant growth leads to higher density of defects acting as pinning or nucleation sites for magnetic domains, which consequently results in deterioration of magnetic properties. It was demonstrated that the use of high miscut substrate is important for fabrication of high quality SRO thin films with low density of crystallographic defects and excellent magnetic properties.


title: A strongly convergent Krasnosel'ski\v{\i}-Mann-type algorithm for finding a common fixed point of a countably infinite family of nonexpansive operators in Hilbert spaces

abstract: In this article, we propose a Krasnosel'ski\v{\i}-Mann-type algorithm for finding a common fixed point of a countably infinite family of nonexpansive operators $(T_n)_{n \geq 0}$ in Hilbert spaces. We formulate an asymptotic property which the family $(T_n)_{n \geq 0}$ has to fulfill such that the sequence generated by the algorithm converges strongly to the element in $\bigcap_{n \geq 0} \operatorname{Fix} T_n$ with minimum norm. Based on this, we derive a forward-backward algorithm that allows variable step sizes and generates a sequence of iterates that converge strongly to the zero with minimum norm of the sum of a maximally monotone operator and a cocoercive one. We demonstrate the superiority of the forward-backward algorithm with variable step sizes over the one with constant step size by means of numerical experiments on variational image reconstruction and split feasibility problems in infinite dimensional Hilbert spaces.


title: Hybrid Text Feature Modeling for Disease Group Prediction using Unstructured Physician Notes

abstract: Existing Clinical Decision Support Systems (CDSSs) largely depend on the availability of structured patient data and Electronic Health Records (EHRs) to aid caregivers. However, in case of hospitals in developing countries, structured patient data formats are not widely adopted, where medical professionals still rely on clinical notes in the form of unstructured text. Such unstructured clinical notes recorded by medical personnel can also be a potential source of rich patient-specific information which can be leveraged to build CDSSs, even for hospitals in developing countries. If such unstructured clinical text can be used, the manual and time-consuming process of EHR generation will no longer be required, with huge person-hours and cost savings. In this paper, we propose a generic ICD9 disease group prediction CDSS built on unstructured physician notes modeled using hybrid word embeddings. These word embeddings are used to train a deep neural network for effectively predicting ICD9 disease groups. Experimental evaluation showed that the proposed approach outperformed the state-of-the-art disease group prediction model built on structured EHRs by 15% in terms of AUROC and 40% in terms of AUPRC, thus proving our hypothesis and eliminating dependency on availability of structured patient data.


title: A User Study of Perceived Carbon Footprint

abstract: We propose a statistical model to understand people's perception of their carbon footprint. Driven by the observation that few people think of CO2 impact in absolute terms, we design a system to probe people's perception from simple pairwise comparisons of the relative carbon footprint of their actions. The formulation of the model enables us to take an active-learning approach to selecting the pairs of actions that are maximally informative about the model parameters. We define a set of 18 actions and collect a dataset of 2183 comparisons from 176 users on a university campus. The early results reveal promising directions to improve climate communication and enhance climate mitigation.


title: Integral geometry of pairs of hyperplanes or lines

abstract: Crofton's formula of integral geometry evaluates the total motion invariant measure of the set of $k$-dimensional planes having nonempty intersection with a given convex body. This note deals with motion invariant measures on sets of pairs of hyperplanes or lines meeting a convex body. Particularly simple results are obtained if, and only if, the given body is of constant width in the first case, and of constant brightness in the second case.


title: {\it Colloquium:} Statistical Mechanics and Thermodynamics at Strong Coupling: Quantum and Classical

abstract: The statistical mechanical description of small systems staying in thermal equilibrium with an environment can be achieved by means of the Hamiltonian of mean force. In contrast to the reduced density matrix of an open quantum system, or the reduced phase-space probability density function of a classical open system, the Hamiltonian of mean force not only characterizes the reduced state but also contains full information about the thermodynamics of the considered open system. The resulting thermodynamic potentials all assume the form as the difference of the potentials for the total system and the bare environment in the absence of the system. In contrast to work as a mechanical notion, one faces several problems with the definition of heat which turns out to be largely ambiguous in the case of strong coupling between system and environment. The general theory of the thermodynamics of open systems, in particular, in view of strong coupling, is reviewed and illustrated with several examples. The vagueness of heat is discussed in the context of the ambiguities in the definitions of a fluctuating internal energy and other fluctuating thermodynamic potentials.


title: Universal spectral features of different classes of random diffusivity processes

abstract: .Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic $1/f^2$-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.


title: Finite Element Simulations of an Elasto-Viscoplastic Model for Clay

abstract: In this paper, we develop an elasto-viscoplastic (EVP) model for clay using the non-associated flow rule. This is accomplished by using a modified form of the Perzyna's overstressed EVP theory, the critical state soil mechanics, and the multi-surface theory. The new model includes six parameters, five of which are identical to those in the critical state soil mechanics model. The other parameter is the generalized nonlinear secondary compression index. The EVP model was implemented in a nonlinear coupled consolidated code using a finite-element numerical algorithm (AFENA). We then tested the model for different clays, such as the Osaka clay, the San Francisco Bay Mud clay, the Kaolin clay, and the Hong Kong Marine Deposit clay. The numerical results show good agreement with the experimental data.


title: Scalable Extreme Deconvolution

abstract: The Extreme Deconvolution method fits a probability density to a dataset where each observation has Gaussian noise added with a known sample-specific covariance, originally intended for use with astronomical datasets. The existing fitting method is batch EM, which would not normally be applied to large datasets such as the Gaia catalog containing noisy observations of a billion stars. We propose two minibatch variants of extreme deconvolution, based on an online variation of the EM algorithm, and direct gradient-based optimisation of the log-likelihood, both of which can run on GPUs. We demonstrate that these methods provide faster fitting, whilst being able to scale to much larger models for use with larger datasets.