title: Dynamically Evolving Bond-Dimensions within the one-site Time-Dependent-Variational-Principle method for Matrix Product States: Towards efficient simulation of non-equilibrium open quantum dynamics

abstract: Understanding the emergent system-bath correlations in non-Markovian and non-perturbative open systems is a theoretical challenge that has benefited greatly from the application of Matrix Product State (MPS) methods. Here, we propose an autonmously adapative variant of the one-site Time-Dependent-Variational-Principle (1TDVP) method for many-body MPS wave-functions in which the local bond-dimensions can evolve to capture growing entanglement 'on the fly'. We achieve this by efficiently examining the effect of increasing each MPS bond-dimension in advance of each dynamic timestep, resulting in an MPS that can dynamically and inhomogeneously restructure itself as the complexity of the dynamics grows across time and space. This naturally leads to more efficient simulations, oviates the need for multiple convergence runs, and, as we demonstrate, is ideally suited to the typical, finite-temperature 'impurity' problems that describe open quantum system connected to multiple environments.


title: Automated Verification of Reactive and Concurrent Programs by Calculation

abstract: Reactive programs combine traditional sequential programming constructs with primitives to allow communication with other concurrent agents. They are ubiquitous in modern applications, ranging from components systems and web services, to cyber-physical systems and autonomous robots. In this paper, we present an algebraic verification strategy for concurrent reactive programs, with a large or infinite state space. We define novel operators to characterise interactions and state updates, and an associated equational theory. With this we can calculate a reactive program's denotational semantics, and thereby facilitate automated proof. Of note is our reasoning support for iterative programs with reactive invariants, based on Kleene algebra, and for parallel composition. We illustrate our strategy by verifying a reactive buffer. Our laws and strategy are mechanised in Isabelle/UTP, our implementation of Hoare and He's Unifying Theories of Programming (UTP) framework, to provide soundness guarantees and practical verification support.


title: Short- and long-term forecasting of electricity prices using embedding of calendar information in neural networks

abstract: Electricity prices strongly depend on seasonality of different time scales, therefore any forecasting of electricity prices has to account for it. Neural networks have proven successful in short-term price-forecasting, but complicated architectures like LSTM are used to integrate the seasonal behaviour. This paper shows that simple neural network architectures like DNNs with an embedding layer for seasonality information can generate a competitive forecast. The embedding-based processing of calendar information additionally opens up new applications for neural networks in electricity trading, such as the generation of price forward curves. Besides the theoretical foundation, this paper also provides an empirical multi-year study on the German electricity market for both applications and derives economical insights from the embedding layer. The study shows that in short-term price-forecasting the mean absolute error of the proposed neural networks with an embedding layer is better than the LSTM and time-series benchmark models and even slightly better as our best benchmark model with a sophisticated hyperparameter optimization. The results are supported by a statistical analysis using Friedman and Holm's tests.


title: Learning "What-if" Explanations for Sequential Decision-Making

abstract: Building interpretable parameterizations of real-world decision-making on the basis of demonstrated behavior -- i.e. trajectories of observations and actions made by an expert maximizing some unknown reward function -- is essential for introspecting and auditing policies in different institutions. In this paper, we propose learning explanations of expert decisions by modeling their reward function in terms of preferences with respect to "what if" outcomes: Given the current history of observations, what would happen if we took a particular action? To learn these cost-benefit tradeoffs associated with the expert's actions, we integrate counterfactual reasoning into batch inverse reinforcement learning. This offers a principled way of defining reward functions and explaining expert behavior, and also satisfies the constraints of real-world decision-making -- where active experimentation is often impossible (e.g. in healthcare). Additionally, by estimating the effects of different actions, counterfactuals readily tackle the off-policy nature of policy evaluation in the batch setting, and can naturally accommodate settings where the expert policies depend on histories of observations rather than just current states. Through illustrative experiments in both real and simulated medical environments, we highlight the effectiveness of our batch, counterfactual inverse reinforcement learning approach in recovering accurate and interpretable descriptions of behavior.


title: Second Order PAC-Bayesian Bounds for the Weighted Majority Vote

abstract: We present a novel analysis of the expected risk of weighted majority vote in multiclass classification. The analysis takes correlation of predictions by ensemble members into account and provides a bound that is amenable to efficient minimization, which yields improved weighting for the majority vote. We also provide a specialized version of our bound for binary classification, which allows to exploit additional unlabeled data for tighter risk estimation. In experiments, we apply the bound to improve weighting of trees in random forests and show that, in contrast to the commonly used first order bound, minimization of the new bound typically does not lead to degradation of the test error of the ensemble.


title: Learning Common Harmonic Waves on Stiefel Manifold -- A New Mathematical Approach for Brain Network Analyses

abstract: Converging evidence shows that disease-relevant brain alterations do not appear in random brain locations, instead, its spatial pattern follows large scale brain networks. In this context, a powerful network analysis approach with a mathematical foundation is indispensable to understand the mechanism of neuropathological events spreading throughout the brain. Indeed, the topology of each brain network is governed by its native harmonic waves, which are a set of orthogonal bases derived from the Eigen-system of the underlying Laplacian matrix. To that end, we propose a novel connectome harmonic analysis framework to provide enhanced mathematical insights by detecting frequency-based alterations relevant to brain disorders. The backbone of our framework is a novel manifold algebra appropriate for inference across harmonic waves that overcomes the limitations of using classic Euclidean operations on irregular data structures. The individual harmonic difference is measured by a set of common harmonic waves learned from a population of individual Eigen systems, where each native Eigen-system is regarded as a sample drawn from the Stiefel manifold. Specifically, a manifold optimization scheme is tailored to find the common harmonic waves which reside at the center of Stiefel manifold. To that end, the common harmonic waves constitute the new neuro-biological bases to understand disease progression. Each harmonic wave exhibits a unique propagation pattern of neuro-pathological burdens spreading across brain networks. The statistical power of our novel connectome harmonic analysis approach is evaluated by identifying frequency-based alterations relevant to Alzheimer's disease, where our learning-based manifold approach discovers more significant and reproducible network dysfunction patterns compared to Euclidian methods.


title: Coupling Learning of Complex Interactions

abstract: Complex applications such as big data analytics involve different forms of coupling relationships that reflect interactions between factors related to technical, business (domain-specific) and environmental (including socio-cultural and economic) aspects. There are diverse forms of couplings embedded in poor-structured and ill-structured data. Such couplings are ubiquitous, implicit and/or explicit, objective and/or subjective, heterogeneous and/or homogeneous, presenting complexities to existing learning systems in statistics, mathematics and computer sciences, such as typical dependency, association and correlation relationships. Modeling and learning such couplings thus is fundamental but challenging. This paper discusses the concept of coupling learning, focusing on the involvement of coupling relationships in learning systems. Coupling learning has great potential for building a deep understanding of the essence of business problems and handling challenges that have not been addressed well by existing learning theories and tools. This argument is verified by several case studies on coupling learning, including handling coupling in recommender systems, incorporating couplings into coupled clustering, coupling document clustering, coupled recommender algorithms and coupled behavior analysis for groups.


title: Holographic theories at finite $\theta$-angle, CP-violation, glueball spectra and strong-coupling instabilities

abstract: A general class of holographic theories with a nontrivial $\theta$-angle are analyzed. The instanton density operator is dual to a bulk axion field. We calculate the ground-state solutions with nontrivial source, $a_{UV}$, for the axion, for both steep and soft dilaton potentials in the IR, and both in $d=3$ and $d=4$. We find all cases to be qualitatively similar. We also calculate the spin$=2,0$ glueball spectra and show that the glueball masses monotonically decrease as functions of $a_{UV}$ (or $\theta$-angle). The slopes of glueball masses are different, generically, in different potentials. In the case of steep dilaton potentials, the glueball (masses)$^2$ turn negative before the maximum of $a_{UV}$ is attained. We interpret this as a signal for a favored instanton condensation in the bulk. We also investigate strong CP-violation in the effective glueball action.


title: Effects of asymmetric splayed columnar defects on the anomalous peak effect in Ba$_{0.6}$K$_{0.4}$Fe$_{2}$As$_{2}$

abstract: In order to enhance the critical current density ($J_c$) of superconductors, introduction of columnar defects (CDs) through swift-particle irradiation is effective. By dispersing the direction of CDs (splayed CDs), not only further enhancement of $J_c$ has been confirmed but also an anomalous peak effect (APE) in $J_c$ at a certain magnetic field determined by the irradiation dose was observed. It has been proposed that the APE arises from the suppression of kink motion and/or the effects of vortex entanglement in systems with splayed CDs. In this study, we measure $J_c$ properties of optimally K-doped Ba-122-type iron-based superconductor Ba$_{0.6}$K$_{0.4}$Fe$_{2}$As$_{2}$ single crystals with splayed CDs that are introduced asymmetrically with respect to the $c$-axis by irradiating 2.6 GeV U ions. We discuss the significance of the average direction of splayed CDs to the APE. We also discuss the relationship between the APE and the entanglement of vortices.


title: Characterization of New Silicon Photomultipliers with Low Dark Noise at Low Temperature

abstract: Silicon photomultipliers (SiPMs) have a low radioactivity, compact geometry, low operation voltage, and reasonable photo-detection efficiency for vacuum ultraviolet light (VUV). Therefore it has the potential to replace photomultiplier tubes (PMTs) for future dark matter experiments with liquid xenon (LXe). However, SiPMs have nearly two orders of magnitude higher dark count rate (DCR) compared to that of PMTs at the LXe temperature ($\sim$ 165 K). This type of high DCR mainly originates from the carriers that are generated by band-to-band tunneling effect. To suppress the tunneling effect, we have developed a new SiPM with lowered electric field strength in cooperation with Hamamatsu Photonics K. K. and characterized its performance in a temperature range of 153 K to 298 K. We demonstrated that the newly developed SiPMs had 6--54 times lower DCR at low temperatures compared to that of the conventional SiPMs.


title: A Novel adaptive optimization of Dual-Tree Complex Wavelet Transform for Medical Image Fusion

abstract: In recent years, many research achievements are made in the medical image fusion field. Fusion is basically extraction of best of inputs and conveying it to the output. Medical Image fusion means that several of various modality image information is comprehended together to form one image to express its information. The aim of image fusion is to integrate complementary and redundant information. In this paper, a multimodal image fusion algorithm based on the dual-tree complex wavelet transform (DT-CWT) and adaptive particle swarm optimization (APSO) is proposed. Fusion is achieved through the formation of a fused pyramid using the DTCWT coefficients from the decomposed pyramids of the source images. The coefficients are fused by the weighted average method based on pixels, and the weights are estimated by the APSO to gain optimal fused images. The fused image is obtained through conventional inverse dual-tree complex wavelet transform reconstruction process. Experiment results show that the proposed method based on adaptive particle swarm optimization algorithm is remarkably better than the method based on particle swarm optimization. The resulting fused images are compared visually and through benchmarks such as Entropy (E), Peak Signal to Noise Ratio, (PSNR), Root Mean Square Error (RMSE), Standard deviation (SD) and Structure Similarity Index Metric (SSIM) computations.


title: On the regularization of Cauchy-type integral operators via the density interpolation method and applications

abstract: This paper presents a regularization technique for the high order efficient numerical evaluation of nearly singular, principal-value, and finite-part Cauchy-type integral operators. By relying on the Cauchy formula, the Cauchy-Goursat theorem, and on-curve Taylor interpolations of the input density, the proposed methodology allows to recast the Cauchy and associated integral operators as smooth contour integrals. As such, they can be accurately evaluated everywhere in the complex plane -- including at problematic points near and on the contour -- by means of elementary quadrature rules. Applications of the technique to the evaluation of the Laplace layer potentials and related integral operators, as well as to the computation conformal mappings, are examined in detail. The former application, in particular, amounts to a significant improvement over the recently introduced harmonic density interpolation method. Spectrally accurate discretization approaches for smooth and piecewise smooth contours are presented. A variety of numerical examples, including the solution of weakly singular and hypersingular Laplace boundary integral equations, and the evaluation of challenging conformal mappings, demonstrate the effectiveness and accuracy of the density interpolation method in this context.


title: A machine learning aided global diagnostic and comparative tool to assess effect of quarantine control in Covid-19 spread

abstract: We have developed a globally applicable diagnostic Covid-19 model by augmenting the classical SIR epidemiological model with a neural network module. Our model does not rely upon previous epidemics like SARS/MERS and all parameters are optimized via machine learning algorithms employed on publicly available Covid-19 data. The model decomposes the contributions to the infection timeseries to analyze and compare the role of quarantine control policies employed in highly affected regions of Europe, North America, South America and Asia in controlling the spread of the virus. For all continents considered, our results show a generally strong correlation between strengthening of the quarantine controls as learnt by the model and actions taken by the regions' respective governments. Finally, we have hosted our quarantine diagnosis results for the top 70 affected countries worldwide, on a public platform, which can be used for informed decision making by public health officials and researchers alike.


title: Thermodynamically Consistent Darcy-Brinkman-Forchheimer Framework in Matrix Acidization

abstract: Matrix acidization is an important technique to enhance oil production at the tertiary recovery stage, and its numerical simulation is never concluded. From one of the earliest models, i.e. the two-scale model (Darcy framework), the Darcy-Brinkman-Forchheimer (DBF) framework is developed by adding Brinkman term and Forchheimer term to the momentum conservation equation. However, in the momentum conservation equation of the DBF framework, porosity is put outside of the time derivation term, which cannot describe the change of porosity well. Thus, this work changes the expression so that the modified momentum conservation equation can satisfy Newton's second law. The modified framework is called improved DBF framework. Furthermore, based on the improved DBF framework, the thermal DBF framework is given by introducing the energy balance equation to the improved DBF framework. Both of the frameworks are verified by the former works through numerical experiments and chemical experiments in labs. Parallelization to the codes of the complicated frameworks is also realized, and good scalability can be achieved.


title: Evaluating the reliability of acoustic speech embeddings

abstract: Speech embeddings are fixed-size acoustic representations of variable-length speech sequences. They are increasingly used for a variety of tasks ranging from information retrieval to unsupervised term discovery and speech segmentation. However, there is currently no clear methodology to compare or optimise the quality of these embeddings in a task-neutral way. Here, we systematically compare two popular metrics, ABX discrimination and Mean Average Precision (MAP), on 5 languages across 17 embedding methods, ranging from supervised to fully unsupervised, and using different loss functions (autoencoders, correspondence autoencoders, siamese). Then we use the ABX and MAP to predict performances on a new downstream task: the unsupervised estimation of the frequencies of speech segments in a given corpus. We find that overall, ABX and MAP correlate with one another and with frequency estimation. However, substantial discrepancies appear in the fine-grained distinctions across languages and/or embedding methods. This makes it unrealistic at present to propose a task-independent silver bullet method for computing the intrinsic quality of speech embeddings. There is a need for more detailed analysis of the metrics currently used to evaluate such embeddings.


title: A tumor growth model with autophagy: the reaction-(cross-)diffusion system and its free boundary limit

abstract: In this paper, we propose a tumor growth model to incorporate and investigate the spatial effects of autophagy. The cells are classified into two phases: normal cells and autophagic cells, whose dynamics are also coupled with the nutrients. First, we construct a reaction-(cross-)diffusion system describing the evolution of cell densities, where the drift is determined by the negative gradient of the joint pressure, and the reaction terms manifest the unique mechanism of autophagy. Next, in the incompressible limit, such a cell density model naturally connects to a free boundary system, describing the geometric motion of the tumor region. Analyzing the free boundary model in a special case, we show that the ratio of the two phases of cells exponentially converges to a ``well-mixed" limit. Within this ``well-mixed" limit, we obtain an analytical solution of the free boundary system which indicates the exponential growth of the tumor size in the presence of autophagy in contrast to the linear growth without it. Numerical simulations are also provided to illustrate the analytical properties and to explore more scenarios.


title: Combining Deep Reinforcement Learning and Search for Imperfect-Information Games

abstract: The combination of deep reinforcement learning and search at both training and test time is a powerful paradigm that has led to a number of successes in single-agent settings and perfect-information games, best exemplified by AlphaZero. However, prior algorithms of this form cannot cope with imperfect-information games. This paper presents ReBeL, a general framework for self-play reinforcement learning and search that provably converges to a Nash equilibrium in any two-player zero-sum game. In the simpler setting of perfect-information games, ReBeL reduces to an algorithm similar to AlphaZero. Results in two different imperfect-information games show ReBeL converges to an approximate Nash equilibrium. We also show ReBeL achieves superhuman performance in heads-up no-limit Texas hold'em poker, while using far less domain knowledge than any prior poker AI.


title: Tutorials on Advanced Optimization Methods

abstract: This material provides thorough tutorials on some optimization techniques frequently used in various engineering disciplines, including convex optimization, linearization techniques and mixed-integer linear programming, robust optimization, and equilibrium/game problems. It discusses how to reformulate a difficult (non-convex, multi-agent, min-max) problem to a solver-compatible form (semidefinite program, mixed-integer linear program) via convexification, linearization, and decomposition, so the original problem can be reliably solved by commercial/open-source software. Fundamental algorithms are not the main focus. This material is a good reference for self-learners who have basic knowledge in linear algebra and linear programming. It is one of the main references for an optimization course taught at Tsinghua University.


title: Numerical approximations and error analysis of the Cahn-Hilliard equation with reaction rate dependent dynamic boundary conditions

abstract: We consider numerical approximations and error analysis for the Cahn-Hilliard equation with reaction rate dependent dynamic boundary conditions (P. Knopf et. al., arXiv, 2020). Based on the stabilized linearly implicit approach, a first-order in time, linear and energy stable scheme for solving this model is proposed. The corresponding semi-discretized-in-time error estimates for the scheme are also derived. Numerical experiments, including the comparison with the former work, the convergence results for the relaxation parameter $K\rightarrow0$ and $K\rightarrow\infty$ and the accuracy tests with respect to the time step size, are performed to validate the accuracy of the proposed scheme and the error analysis.


title: Reconstruction Regularized Deep Metric Learning for Multi-label Image Classification

abstract: In this paper, we present a novel deep metric learning method to tackle the multi-label image classification problem. In order to better learn the correlations among images features, as well as labels, we attempt to explore a latent space, where images and labels are embedded via two unique deep neural networks, respectively. To capture the relationships between image features and labels, we aim to learn a \emph{two-way} deep distance metric over the embedding space from two different views, i.e., the distance between one image and its labels is not only smaller than those distances between the image and its labels' nearest neighbors, but also smaller than the distances between the labels and other images corresponding to the labels' nearest neighbors. Moreover, a reconstruction module for recovering correct labels is incorporated into the whole framework as a regularization term, such that the label embedding space is more representative. Our model can be trained in an end-to-end manner. Experimental results on publicly available image datasets corroborate the efficacy of our method compared with the state-of-the-arts.


title: Coexistence of Type-II and Type-IV Dirac Fermions in SrAgBi

abstract: Relativistic massless Weyl and Dirac fermions have isotropic and linear dispersion relations to maintain Poincar\'{e} symmetry, which is the most basic symmetry in high-energy physics. The situation in condensed matter physics is less constrained; only certain subgroups of Poincar\'{e} symmetry -- the 230 space groups that exist in 3D lattices -- need be respected. Then, the free fermionic excitations that have no high energy analogues could exist in solid state systems. Here, We discovered a type of nonlinear Dirac fermion without high-energy analogue in SrAgBi and named it type-IV Dirac fermion. The type-IV Dirac fermion has a nonlinear dispersion relationship and is similar to the type-II Dirac fermion, which has electron pocket and hole pocket. The effective model for the type-IV Dirac fermion is also found. It is worth pointing out that there is a type-II Dirac fermion near this new Dirac fermion. So we used two models to describe the coexistence of these two Dirac fermions. Topological surface states of these two Dirac points are also calculated. We envision that our findings will stimulate researchers to study novel physics of type-IV Dirac fermions, as well as the interplay of type-II and type-IV Dirac fermions.


title: The Wage Premium of Communist Party Membership: Evidence from China

abstract: Social status and political connections could confer large economic benefits to an individual. Previous studies focused on China examine the relationship between Communist party membership and earnings and find a positive correlation. However, this correlation may be partly or totally spurious, thereby generating upwards-biased estimates of the importance of political party membership. Using data from three surveys spanning more than three decades, we estimate the causal effect of Chinese party membership on monthly earnings in in China. We find that, on average, membership in the Communist party of China increases monthly earnings and we find evidence that the wage premium has grown in recent years. We explore for potential mechanisms and we find suggestive evidence that improvements in one's social network, acquisition of job-related qualifications and improvement in one's social rank and life satisfaction likely play an important role. (JEL D31, J31, P2)


title: Cantor sets and fields of reals

abstract: Our main result is a construction of four families C_1,C_2,B_1,B_2 which are equipollent with the power set of the real line R and satisfy the following properties. (i) The members of the families are proper subfields of R whose algebraic closures equal the field C. (ii) Each field in C_1vC_2 contains a Cantor set. (iii) Each field in B_1vB_2 is a Bernstein set. (iv) All fields in C_1vB_1 are isomorphic. (v) If K,L are fields in C_2vB_2 then K is isomorphic to a subfield of L only in the trivial case K=L.


title: Differentiable Manifold Reconstruction for Point Cloud Denoising

abstract: 3D point clouds are often perturbed by noise due to the inherent limitation of acquisition equipments, which obstructs downstream tasks such as surface reconstruction, rendering and so on. Previous works mostly infer the displacement of noisy points from the underlying surface, which however are not designated to recover the surface explicitly and may lead to sub-optimal denoising results. To this end, we propose to learn the underlying manifold of a noisy point cloud from differentiably subsampled points with trivial noise perturbation and their embedded neighborhood feature, aiming to capture intrinsic structures in point clouds. Specifically, we present an autoencoder-like neural network. The encoder learns both local and non-local feature representations of each point, and then samples points with low noise via an adaptive differentiable pooling operation. Afterwards, the decoder infers the underlying manifold by transforming each sampled point along with the embedded feature of its neighborhood to a local surface centered around the point. By resampling on the reconstructed manifold, we obtain a denoised point cloud. Further, we design an unsupervised training loss, so that our network can be trained in either an unsupervised or supervised fashion. Experiments show that our method significantly outperforms state-of-the-art denoising methods under both synthetic noise and real world noise. The code and data are available at https://github.com/luost26/DMRDenoise


title: HeAT -- a Distributed and GPU-accelerated Tensor Framework for Data Analytics

abstract: To cope with the rapid growth in available data, the efficiency of data analysis and machine learning libraries has recently received increased attention. Although great advancements have been made in traditional array-based computations, most are limited by the resources available on a single computation node. Consequently, novel approaches must be made to exploit distributed resources, e.g. distributed memory architectures. To this end, we introduce HeAT, an array-based numerical programming framework for large-scale parallel processing with an easy-to-use NumPy-like API. HeAT utilizes PyTorch as a node-local eager execution engine and distributes the workload on arbitrarily large high-performance computing systems via MPI. It provides both low-level array computations, as well as assorted higher-level algorithms. With HeAT, it is possible for a NumPy user to take full advantage of their available resources, significantly lowering the barrier to distributed data analysis. When compared to similar frameworks, HeAT achieves speedups of up to two orders of magnitude.


title: Sequential design of multi-fidelity computer experiments: maximizing the rate of stepwise uncertainty reduction

abstract: This article deals with the sequential design of experiments for (deterministic or stochastic) multi-fidelity numerical simulators, that is, simulators that offer control over the accuracy of simulation of the physical phenomenon or system under study. Very often, accurate simulations correspond to high computational efforts whereas coarse simulations can be obtained at a smaller cost. In this setting, simulation results obtained at several levels of fidelity can be combined in order to estimate quantities of interest (the optimal value of the output, the probability that the output exceeds a given threshold...) in an efficient manner. To do so, we propose a new Bayesian sequential strategy called Maximal Rate of Stepwise Uncertainty Reduction (MR-SUR), that selects additional simulations to be performed by maximizing the ratio between the expected reduction of uncertainty and the cost of simulation. This generic strategy unifies several existing methods, and provides a principled approach to develop new ones. We assess its performance on several examples, including a computationally intensive problem of fire safety analysis where the quantity of interest is the probability of exceeding a tenability threshold during a building fire.


title: The role of infrared radiation pressure in shaping dusty winds in AGN

abstract: The detection of dusty winds dominating the infrared emission of AGN on parsec scales has revealed the limitations of traditional radiative transfer models based on a toroidal distribution of dusty gas. A new, more complex, dynamical structure is emerging and the physical origin of such dusty winds has to be critically assessed. We present a semi-analytical model to test the hypothesis of radiatively accelerated dusty winds launched by the AGN and by the heated dust itself. The model consists of an AGN and an infrared radiating dusty disk, the latter being the primary mass reservoir for the outflow. We calculate the trajectories of dusty gas clumps in this environment, accounting for both gravity and the AGN radiation as well as the re-radiation by the hot dusty gas clouds themselves. We find that the morphology consists of a disk of material that orbits with sub-Keplerian velocities and a hyperboloid polar wind launched at the inner edge of the dusty disk. This is consistent with high-angular resolution infrared and sub-mm observations of some local Seyfert AGN. The strength of the wind and its orientation depend on the Eddington ratio and the column density of the dusty clumps, which is in agreement with proposed radiation regulated obscuration models developed for the X-ray obscuring material around AGN.


title: Supercontinuum generation assisted by dispersive waves trapping in dispersion-managed integrated silicon waveguides

abstract: Compact chip-scale comb sources are of significant interest for many practical applications. Here, we experimentally study the generation of supercontinuum (SC) in an axially varying integrated waveguide. We show that the local tuning of the dispersion enables the continuous blue shift of dispersive waves thanks to their trapping by the strongly compressed pump pulse. This mechanism provides new insight into supercontinuum generation in a dispersion varying integrated waveguide. Pumped close to 2.2 $\mu$m in the femtosecond regime and at a pulse energy of $\sim$ 4 pJ, the output spectrum extends from 1.1 $\mu$m up to 2.76 $\mu$m and shows good coherence properties. Octave-spanning SC is also observed at input energy as low as $\sim$ 0.9 pJ. We show that the supercontinuum is more robust against variations of the input pulse parameters and is also spectrally flatter in our numerically optimized waveguide than in fixed-width waveguides. This research demonstrates the potential of dispersion varying waveguides for coherent SC generation and paves the way for integrated low power applications, such as chip-scale frequency comb generation, precision spectroscopy, optical frequency metrology, and wide-band wavelength division multiplexing in the near-infrared.


title: Methods for detecting Order-by-Disorder transitions: the example of the Domino model

abstract: Detecting the zero-temperature thermal Order-by-Disorder transition in classical magnetic systems is notably difficult. We propose a method to probe this transition in an indirect way. The idea is to apply adequate magnetic fields to transform the zero temperature transition into a finite temperature sharp crossover, which should be much easier to observe and characterise with usual laboratory methods.


title: Infrared Spectra of Pyroxenes (Crystalline Chain Silicates) at Room Temperature

abstract: Pyroxene crystals are common in meteorites but few compositions have been recognized in astronomical environments. We present quantitative room-temperature spectra of 17 Mg-- Fe-- and Ca--bearing ortho- and clinopyroxenes, and a Ca-pyroxenoid in order to discern trends indicative of crystal structure and a wide range of composition. Data are produced using a Diamond Anvil Cell: our band strengths are up to 6 times higher than those measured in KBr or polyethylene dispersions, which include variations in path length (from grain size) and surface reflections that are not addressed in data processing. Pyroxenes have varied spectra: only two bands, at 10.22~$\mu$m and 15.34~$\mu$m in enstatite (En$_{99}$), are common to all. Peak-wavelengths generally increase as Mg is replaced by Ca or Fe. However, two bands in MgFe-pyroxenes shift to shorter wavelengths as the Fe component increases from 0 to 60 per cent. A high-intensity band shifts from 11.6~$\mu$m to 11.2~$\mu$m and remains at 11.2~$\mu$m as Fe increases to 100~per~cent; it resembles an astronomical feature normally identified with olivine or forsterite. The distinctive pyroxene bands between 13~ and 16~$\mu$m show promise for their identification in MIRI spectra obtained with JWST. The many pyroxene bands between 40 and 80~$\mu$m could be diagnositic of silicate mineralogy if data were obtained with the proposed SPICA telescope. Our data indicate that comparison between room-temperature laboratory bands for enstatite and cold $\sim 10-K$ astronomical dust features at wavelengths $\gtrsim 28~\mu$m can result in the identification of (Mg,Fe)- pyroxenes that contain 7--15 % less Fe-- than their true values because some temperature shifts mimic some compositional shifts. Therefore some astronomical silicates may contain more Fe, and less Mg, than previously thought.


title: Designing a portable radiofrequency source for magnetic hyperthermia applications

abstract: In the 21st century, new and alternative non-radioactive methods for cancer treatment have been developed. One of this is magnetic hyperthermia (MH), which is based on the ability of magnetic nanoparticles (MNP) to transduce electromagnetic energy into heat. If MNP are inside cancer cells, heat induces cells death. The success of HM as cancer therapy depends on two factors, first, the NPM properties and second, the frequency and amplitude of the radio frequency (RF) applied field. Therefore, in this work a design of a portable radio frequency source for magnetic hyperthermia experiments is presented. The design consists of two parts, the first part corresponds to the RF source composed by a RLC parallel circuit, the capacity and inductance values have been chosen in order to achieve the RF range of frequency, and the second part corresponds to the design of a magnetization sensor circuit. While the first circuit produces the RF field, the second circuit serves as a detector of the NPM induced magnetization. The development of this radiofrequency source makes it possible to determine the applicability of an NPM system in the treatment of cancer by magnetic hyperthermia, opening doors for future research focused on the application of nanotechnology in biomedicine.


title: MADGAN: unsupervised Medical Anomaly Detection GAN using multiple adjacent brain MRI slice reconstruction

abstract: Unsupervised learning can discover various unseen abnormalities, relying on large-scale unannotated medical images of healthy subjects. Towards this, unsupervised methods reconstruct a 2D/3D single medical image to detect outliers either in the learned feature space or from high reconstruction loss. However, without considering continuity between multiple adjacent slices, they cannot directly discriminate diseases composed of the accumulation of subtle anatomical anomalies, such as Alzheimer's Disease (AD). Moreover, no study has shown how unsupervised anomaly detection is associated with either disease stages, various (i.e., more than two types of) diseases, or multi-sequence Magnetic Resonance Imaging (MRI) scans. Therefore, we propose unsupervised Medical Anomaly Detection Generative Adversarial Network (MADGAN), a novel two-step method using GAN-based multiple adjacent brain MRI slice reconstruction to detect brain anomalies at different stages on multi-sequence structural MRI: (Reconstruction) Wasserstein loss with Gradient Penalty + 100 L1 loss-trained on 3 healthy brain axial MRI slices to reconstruct the next 3 ones-reconstructs unseen healthy/abnormal scans; (Diagnosis) Average L2 loss per scan discriminates them, comparing the ground truth/reconstructed slices. For training, we use two different datasets composed of 1,133 healthy T1-weighted (T1) and 135 healthy contrast-enhanced T1 (T1c) brain MRI scans for detecting AD and brain metastases/various diseases, respectively. Our Self-Attention MADGAN can detect AD on T1 scans at a very early stage, Mild Cognitive Impairment (MCI), with Area Under the Curve (AUC) 0.727, and AD at a late stage with AUC 0.894, while detecting brain metastases on T1c scans with AUC 0.921.


title: Rate-Splitting Multiple Access: Unifying NOMA and SDMA in MISO VLC Channels

abstract: The increased proliferation of connected devices requires a paradigm shift towards the development of innovative technologies for the next generation of wireless systems. One of the key challenges, however, is the spectrum scarcity, owing to the unprecedented broadband penetration rate in recent years. Visible light communications (VLC) has recently emerged as a possible solution to enable high-speed short-range communications. However, VLC systems suffer from several limitations, including the limited modulation bandwidth of light-emitting diodes. Consequently, several multiple access techniques (MA), e.g., space-division multiple access (SDMA) and non-orthogonal multiple access (NOMA), have been considered for VLC networks. Despite their promising multiplexing gains, their performance is somewhat limited. In this article, we first provide an overview of the key MA technologies used in VLC systems. Then, we introduce rate-splitting multiple access (RSMA), which was initially proposed for RF systems and discuss its potentials in VLC systems. Through system modeling and simulations of an RSMA-based two-use scenario, we illustrate the flexibility of RSMA in generalizing NOMA and SDMA, as well as its superiority in terms of weighted sum rate (WSR) in VLC. Finally, we discuss challenges, open issues, and research directions, which will enable the practical realization of RSMA in VLC.


title: Radio Access Technology Characterisation Through Object Detection

abstract: \ac{RAT} classification and monitoring are essential for efficient coexistence of different communication systems in shared spectrum. Shared spectrum, including operation in license-exempt bands, is envisioned in the \ac{5G} standards (e.g., 3GPP Rel. 16). In this paper, we propose a \ac{ML} approach to characterise the spectrum utilisation and facilitate the dynamic access to it. Recent advances in \acp{CNN} enable us to perform waveform classification by processing spectrograms as images. In contrast to other \ac{ML} methods that can only provide the class of the monitored \acp{RAT}, the solution we propose can recognise not only different \acp{RAT} in shared spectrum, but also identify critical parameters such as inter-frame duration, frame duration, centre frequency, and signal bandwidth by using object detection and a feature extraction module to extract features from spectrograms. We have implemented and evaluated our solution using a dataset of commercial transmissions, as well as in a \ac{SDR} testbed environment. The scenario evaluated was the coexistence of WiFi and LTE transmissions in shared spectrum. Our results show that our approach has an accuracy of 96\% in the classification of \acp{RAT} from a dataset that captures transmissions of regular user communications. It also shows that the extracted features can be precise within a margin of 2\%, %of the size of the image, and is capable of detect above 94\% of objects under a broad range of transmission power levels and interference conditions.


title: Time-dependent atomic magnetometry with a recurrent neural network

abstract: We propose to employ a recurrent neural network to estimate a fluctuating magnetic field from continuous optical Faraday rotation measurement on an atomic ensemble. We show that an encoder-decoder architecture neural network can process measurement data and learn an accurate map between recorded signals and the time-dependent magnetic field. The performance of this method is comparable to Kalman filters while it is free of the theory assumptions that restrict their application to particular measurements and physical systems.


title: Forbidden subspaces for level-1 QAOA and IQP circuits

abstract: We present a thorough investigation of problems that can be solved exactly with the level-1 Quantum Approximate Optimization Algorithm (QAOA). To this end we implicitly define a class of problem Hamiltonians that employed as phase separator in a level-1 QAOA circuit provide unit overlap with a target subspace spanned by a set of computational basis states. For one-dimensional target subspaces we identify instances within the implicitly defined class of Hamiltonians for which Quantum Annealing (QA) and Simulated Annealing (SA) have an exponentially small probability to find the solution. Consequently, our results define a first demarcation line between QAOA, QA and SA, and highlight the fundamental differences between an interference-based search heuristic such as QAOA and heuristics that are based on thermal and quantum fluctuations like SA and QA respectively. Moreover, for two-dimensional solution subspaces we are able to show that the depth of the QAOA circuit grows linearly with the Hamming distance between the two target states. We further show that there are no genuine solutions for target subspaces of dimension higher than $2$ and smaller than $2^n$. We also transfer these results to Instantaneous Quantum Polynomial (IQP) circuits.


title: Lackadaisical quantum walks on triangular and honeycomb 2D grids

abstract: In the typical model, a discrete-time coined quantum walk search has the same running time of $O(\sqrt{N} \log{N})$ for 2D rectangular, triangular and honeycomb grids. It is known that for 2D rectangular grid the running time can be improved to $O(\sqrt{N \log{N}})$ using several different techniques. One of such techniques is adding a self-loop of weight $4/N$ to each vertex (i.e. making the walk lackadaisical). In this paper we apply lackadaisical approach to quantum walk search on triangular and honeycomb 2D grids. We show that for both types of grids adding a self-loop of weight $6/N$ and $3/N$ for triangular and honeycomb grids, respectively, results in $O(\sqrt{N \log{N}})$ running time.


title: Morse-Bott Theory on posets and an homological Lusternik-Schnirelmann Theorem

abstract: We develop Morse-Bott theory on posets, generalizing both discrete Morse-Bott theory for regular complexes and Morse theory on posets. Moreover, we prove a Lusternik-Schnirelmann theorem for general matchings on posets, in particular, for Morse-Bott functions.


title: Are low frequency macroeconomic variables important for high frequency electricity prices?

abstract: Recent research finds that forecasting electricity prices is very relevant. In many applications, it might be interesting to predict daily electricity prices by using their own lags or renewable energy sources. However, the recent turmoil of energy prices and the Russian-Ukrainian war increased attention in evaluating the relevance of industrial production and the Purchasing Managers' Index output survey in forecasting the daily electricity prices. We develop a Bayesian reverse unrestricted MIDAS model which accounts for the mismatch in frequency between the daily prices and the monthly macro variables in Germany and Italy. We find that the inclusion of macroeconomic low frequency variables is more important for short than medium term horizons by means of point and density measures. In particular, accuracy increases by combining hard and soft information, while using only surveys gives less accurate forecasts than using only industrial production data.


title: Non minimal D-type conformal matter compactified on three punctured spheres

abstract: We study compactifications of $6d$ non minimal $(D_{p+3},D_{p+3})$ type conformal matter. These can be described by $N$ M5-branes probing a $D_{p+3}$-type singularity. We derive $4d$ Lagrangians corresponding to compactifications of such $6d$ SCFTs on three punctured spheres (trinions) with two maximal punctures and one minimal puncture. The trinion models are described by simple $\mathcal{N}=1$ quivers with $SU(2N)$ gauge nodes. We derive the trinion Lagrangians using RG flows between the aforementioned $6d$ SCFTs with different values of $p$ and their relations to matching RG flows in their compactifications to $4d$. The suggested trinions are shown to reduce to known models in the minimal case of $N=1$. Additional checks are made to show the new minimal punctures uphold the expected S-duality between models in which we exchange two such punctures. We also show that closing the new minimal puncture leads to expected flux tube models.


title: Jumps and Coalescence in the Continuum: a Numerical Study

abstract: The dynamics is studied of an infinite continuum system of jumping and coalescing point particles. In the course of jumps, the particles repel each other whereas their coalescence is free. As the equation of motion we take a kinetic equation, derived by a scaling procedure from the microscopic Fokker-Planck equation corresponding to this kind of motion. The result of the paper is the numerical study (by the Runge-Kutta method) of the solutions of the kinetic equation revealing a number of interesting peculiarities of the dynamics and clarifying the particular role of the jumps and the coalescence in the system's evolution. Possible nontrivial stationary states are also found and analyzed.


title: Asymptotic expansions for the coefficients of extremal quasimodular forms and a conjecture of Kaneko and Koike

abstract: Extremal quasimodular forms have been introduced by M.~Kaneko and M.Koike as as quasimodular forms which have maximal possible order of vanishing at $i\infty$. We show an asymptotic formula for the Fourier coefficients of such forms. This formula is then used to show that all but finitely many Fourier coefficients of such forms of depth $\leq4$ are positive, which partially solves a conjecture stated by M.~Kaneko and M.Koike. Numerical experiments based on constructive estimates confirm the conjecture for weights $\leq200$ and depths between $1$ and $4$.


title: Understanding controlled EV charging impacts using scenario-based forecasting models

abstract: Electrification of transport is a key strategy in reducing carbon emissions. Many countries have adopted policies of complete but gradual transformation to electric vehicles (EVs). However, mass EV adoption also means a spike in load (kW), which in turn can disrupt existing electricity infrastructure. Smart or controlled charging is widely seen as a potential solution to alleviate this stress on existing networks. Learning from the recent EV trials in the UK and elsewhere we take into account two key aspects which are largely ignored in current research: EVs actually charging at any given time and wide range of EV types, especially battery capacity-wise. Taking a minimalistic scenario-based approach, we study forecasting models for mean number of active chargers and mean EV consumption for distinct scenarios. Focusing on residential charging the models we consider range from simple regression models to more advanced machine and deep learning models such as XGBoost and LSTMs. We then use these models to evaluate the impacts of different levels of future EV penetration on a specimen distribution transformer that captures typical real-world scenarios. In doing so, we also initiate the study of different types of controlled charging when fully controlled charging is not possible. This aligns with the outcomes from recent trials which show that a sizeable proportion of EV owners may not prefer fully controlled centralized charging. We study two possible control regimes and show that one is more beneficial from load-on-transformer point of view, while the other may be preferred for other objectives. We show that a minimum of 60% control is required to ensure that transformers are not overloaded during peak hours.


title: Covert Millimeter-Wave Communication: Design Strategies and Performance Analysis

abstract: In this paper, we investigate covert communication over millimeter-wave (mmWave) frequencies. In particular, a mmWave transmitter, referred to as Alice, attempts to reliably communicate to a receiver, referred to as Bob, while hiding the existence of communication from a warden, referred to as Willie. In this regard, operating over the mmWave bands not only increases the covertness thanks to directional beams, but also increases the transmission data rates given much more available bandwidths and enables ultra-low form factor transceivers due to the lower wavelengths used compared to the conventional radio frequency (RF) counterpart. We first assume that the transmitter Alice employs two independent antenna arrays in which one of the arrays is to form a directive beam for data transmission to Bob. The other antenna array is used by Alice to generate another beam toward Willie as a jamming signal while changing the transmit power independently across the transmission blocks in order to achieve the desired covertness. For this dual-beam setup, we characterize Willie's detection error rate with the optimal detector and the closed-form of its expected value from Alice's perspective. We then derive the closed-form expression for the outage probability of the Alice-Bob link, which enables characterizing the optimal covert rate that can be achieved using the proposed setup. We further obtain tractable forms for the ergodic capacity of the Alice-Bob link involving only one-dimensional integrals that can be computed in closed forms for most ranges of the channel parameters. Finally, we highlight how the results can be extended to more practical scenarios, particularly to the cases where perfect information about the location of the passive warden is not available.


title: High order, semi-implicit, energy stable schemes for gradient flows

abstract: We introduce a class of high order accurate, semi-implicit Runge-Kutta schemes in the general setting of evolution equations that arise as gradient flow for a cost function, possibly with respect to an inner product that depends on the solution, and we establish their energy stability. This class includes as a special case high order, unconditionally stable schemes obtained via convexity splitting. The new schemes are demonstrated on a variety of gradient flows, including partial differential equations that are gradient flow with respect to the Wasserstein (mass transport) distance.


title: Efficient, sub-4-cycle, 1-{\mu}m-pumped optical parametric amplifier at 10 {\mu}m based on BaGa4S7

abstract: We report on a {\mu}J-scale mid-infrared optical parametric amplifier (OPA) based on the recently developed wide-bandgap orthorhombic crystal, BaGa4S7 (BGS), and directly compare its performance to that of LiGaS2 (LGS) in the same OPA setup. The source is based on a single OPA stage amplifying supercontinuum seed pulses with a quantum efficiency of 29{\%} at an idler wavelength of 10 {\mu}m, featuring nominally carrier-envelope phase-stable pulses. As a result of pumping the OPA directly at 1 {\mu}m, the overall conversion efficiency far exceeds that of traditional schemes based on OPA's followed by difference frequency generation. Chirp compensation using bulk germanium resulted in 126-fs pulses covering the 7.6-11.5-{\mu}m spectral range. BGS holds great promise for power scaling due to its availability in larger single-crystal sizes than LGS.


title: Phylogenetic networks as circuits with resistance distance

abstract: Phylogenetic networks are notoriously difficult to reconstruct. Here we suggest that it can be useful to view unknown genetic distance along edges in phylogenetic networks as analogous to unknown resistance in electric circuits. This resistance distance, well known in graph theory, turns out to have nice mathematical properties which allow the precise reconstruction of networks. Specifically we show that the resistance distance for a weighted 1-nested network is Kalmanson, and that the unique associated circular split network fully represents the splits of the original phylogenetic network (or circuit). In fact, this full representation corresponds to a face of the balanced minimal evolution polytope for level-1 networks. Thus the unweighted class of the original network can be reconstructed by either the greedy algorithm neighbor-net or by linear programming over a balanced minimal evolution polytope. We begin study of 2-nested networks with both minimum path and resistance distance, and include some counting results for 2-nested networks.


title: Using Machine Learning to Identify the Most At-Risk Students in Physics Classes

abstract: Machine learning algorithms have recently been used to predict students' performance in an introductory physics class. The prediction model classified students as those likely to receive an A or B or students likely to receive a grade of C, D, F or withdraw from the class. Early prediction could better allow the direction of educational interventions and the allocation of educational resources. However, the performance metrics used in that study become unreliable when used to classify whether a student would receive an A, B or C (the ABC outcome) or if they would receive a D, F or withdraw (W) from the class (the DFW outcome) because the outcome is substantially unbalanced with between 10\% to 20\% of the students receiving a D, F, or W. This work presents techniques to adjust the prediction models and alternate model performance metrics more appropriate for unbalanced outcome variables. These techniques were applied to three samples drawn from introductory mechanics classes at two institutions ($N=7184$, $1683$, and $926$). Applying the same methods as the earlier study produced a classifier that was very inaccurate, classifying only 16\% of the DFW cases correctly; tuning the model increased the DFW classification accuracy to 43\%. Using a combination of institutional and in-class data improved DFW accuracy to 53\% by the second week of class. As in the prior study, demographic variables such as gender, underrepresented minority status, first-generation college student status, and low socioeconomic status were not important variables in the final prediction models.


title: Fluid leakage in metallic seals

abstract: Metallic seals are crucial machine elements in many important applications, e.g., in ultrahigh vacuum systems. Due to the high elastic modulus of metals, and the surface roughness which exists on all solid surfaces, if no plastic deformation would occur one expects in most cases large fluid flow channels between the contacting metallic bodies, and large fluid leakage. However, in most applications plastic deformation occurs, at least at the asperity level, which allows the surfaces to approach each other to such an extent that fluid leakage often can be neglected. In this study, we present an experimental set-up for studying the fluid leakage in metallic seals. We study the water leakage between a steel sphere and a steel body (seat) with a conical surface. The experimental results are found to be in good quantitative agreement with a (fitting-parameter-free) theoretical model. The theory predicts that the plastic deformations reduce the leak-rate by a factor $\approx 8$


title: Giant planet formation models with a self-consistent treatment of the heavy elements

abstract: We present a new numerical framework to model the formation and evolution of giant planets. The code is based on the further development of the stellar evolution toolkit Modules for Experiments in Stellar Astrophysics (MESA). The model includes the dissolution of the accreted planetesimals/pebbles, which are assumed to be made of water ice, in the planetary gaseous envelope, and the effect of envelope enrichment on the planetary growth and internal structure is computed self-consistently. We apply our simulations to Jupiter and investigate the impact of different heavy-element and gas accretion rates on its formation history. We show that the assumed runaway gas accretion rate significantly affect the planetary radius and luminosity. It is confirmed that heavy-element enrichment leads to shorter formation timescales due to more efficient gas accretion. We find that with heavy-element enrichment Jupiter's formation timescale is compatible with typical disks' lifetimes even when assuming a low heavy-element accretion rate (oligarchic regime). Finally, we provide an approximation for the heavy-element profile in the innermost part of the planet, providing a link between the internal structure and the planetary growth history.


title: On common extensions of valued fields

abstract: Given a valuation $v$ on a field $K$, an extension $\bar{v}$ to an algebraic closure and an extension $w$ to $K(X)$. We want to study the common extensions of $\bar{v}$ and $w$ to $\bar{K}(X)$. First we give a detailed link between the minimal pairs notion and the key polynomials notion. Then we prove that in the case when $w$ is a transcendental extension, then any sequence of key polynomials admits a maximal element, and in case this sequence does not contain a limit key polynomial, then any root of the last key polynomial, describe a common extension.


title: Kramers pairs of Majorana corner states in a topological insulator bilayer

abstract: We consider a system consisting of two tunnel-coupled two-dimensional topological insulators proximitized by a top and bottom superconductor with a phase difference of $\pi$ between them. We show that this system exhibits a time-reversal invariant second-order topological superconducting phase characterized by the presence of a Kramers pair of Majorana corner states at all four corners of a rectangular sample. We furthermore investigate the effect of a weak time-reversal symmetry breaking perturbation and show that an in-plane Zeeman field leads to an even richer phase diagram exhibiting two nonequivalent phases with two Majorana corner states per corner as well as an intermediate phase with only one Majorana corner state per corner. We derive our results analytically from continuum models describing our system. In addition, we also provide independent numerical confirmation of the resulting phases using discretized lattice representations of the models, which allows us to demonstrate the robustness of the topological phases and the Majorana corner states against parameter variations and potential disorder.


title: Effect of momentum anisotropy on quark matter in the quark-meson model

abstract: We investigate the chiral phase structure of quark matter with spheroidal momentum-space anisotropy specified by one anisotropy parameter $\xi$ in the 2+1 flavor quark-meson model. We find that the chiral phase diagram and the location of the critical endpoint (CEP) are affected significantly by the value of $\xi$. With the increase of $\xi$, the CEP is shifted to smaller temperatures and larger quark chemical potentials. And the temperature of the CEP is more sensitive to the anisotropy parameter than the corresponding quark chemical potential, which is opposite to the study for finite system volume effect. Furthermore, the effects of momentum anisotropy on the thermodynamic properties and scalar (pseudoscalar) meson masses are also studied at vanishing quark chemical potential. The numerical results show that an increase of $\xi$ can hinder the restoration of chiral symmetry. We also find that shear viscosity and electrical conductivity decrease as $\xi$ grows. However, bulk viscosity exhibits a significant non-trivial behavior with $\xi$ in the entire temperature domain of interest.


title: Fractional diffusion-wave equation: hidden regularity for weak solutions

abstract: We will give some regularity results about fractional diffusion-wave equations.


title: Zeros of Jensen polynomials and asymptotics for the Riemann xi function

abstract: The classical criterion of Jensen for the Riemann hypothesis is that all of the associated Jensen polynomials have only real zeros. We find a new version of this criterion, using linear combinations of Hermite polynomials, and show that this condition holds in many cases. Detailed asymptotic expansions are given for the required Taylor coefficients of the xi function at $1/2$ as well as related quantities. These results build on those in the recent paper of Griffin, Ono, Rolen and Zagier.


title: Examples of abelian surfaces failing the local-global principle for isogenies

abstract: We provide examples of abelian surfaces over number fields $K$ whose reductions at almost all good primes possess an isogeny of prime degree $\ell$ rational over the residue field, but which themselves do not admit a $K$-rational $\ell$-isogeny. This builds on work of Cullinan and Sutherland. When $K=\mathbb{Q}$, we identify certain weight-$2$ newforms $f$ with quadratic Fourier coefficients whose associated modular abelian surfaces $A_f$ exhibit such a failure of a local-global principle for isogenies.


title: Unconventional bulk superconductivity in YFe$_2$Ge$_2$ single crystals

abstract: Using a new horizontal flux growth technique to produce high quality crystals of the unconventional superconductor YFe$_2$Ge$_2$ has led to a seven-fold reduction in disorder scattering, resulting in mm-sized crystals with residual resistivities $\simeq \SI{0.45}{\micro\ohm\centi\meter}$, resistivity ratios $\simeq 430$ and sharp superconducting heat capacity anomalies. This enables searching multi-probe experiments investigating the normal and superconducting states of YFe$_2$Ge$_2$. Low temperature heat capacity measurements suggest a significant residual Sommerfeld coefficient, consistent with in-gap states induced by residual disorder as predicted for a sign-changing order parameter.


title: VFL: A Verifiable Federated Learning with Privacy-Preserving for Big Data in Industrial IoT

abstract: Due to the strong analytical ability of big data, deep learning has been widely applied to train the collected data in industrial IoT. However, for privacy issues, traditional data-gathering centralized learning is not applicable to industrial scenarios sensitive to training sets. Recently, federated learning has received widespread attention, since it trains a model by only relying on gradient aggregation without accessing training sets. But existing researches reveal that the shared gradient still retains the sensitive information of the training set. Even worse, a malicious aggregation server may return forged aggregated gradients. In this paper, we propose the VFL, verifiable federated learning with privacy-preserving for big data in industrial IoT. Specifically, we use Lagrange interpolation to elaborately set interpolation points for verifying the correctness of the aggregated gradients. Compared with existing schemes, the verification overhead of VFL remains constant regardless of the number of participants. Moreover, we employ the blinding technology to protect the privacy of the gradients submitted by the participants. If no more than n-2 of n participants collude with the aggregation server, VFL could guarantee the encrypted gradients of other participants not being inverted. Experimental evaluations corroborate the practical performance of the presented VFL framework with high accuracy and efficiency.


title: Electrowetting-Controlled Dropwise Condensation with Patterned Electrodes: Physical Principles, Modeling, and Application Perspectives for Fog Harvesting and Enhanced Heat Transfer

abstract: Patterning the wettability of solid surfaces is a successful strategy to control the dropwise condensation of vapor onto partially wetting solid surfaces. We followed the condensation of water vapor onto electrowetting-functionalized surfaces with structured co-planar electrodes. A detailed analysis of the experimental distribution of millions of drops reveals that despite the presence of contact angle hysteresis and the occurrence of random drop coalescence events, the preferential drop position closely follows the evolution of the local minima of the numerically calculated drop size-dependent electrostatic energy landscape in two dimensions. Even subtle transitions between competing preferred locations are properly reproduced by the model. Based on this quantitative understanding of the condensation patterns, we discuss a series of important follow-up steps that need to be taken to demonstrate a reliable performance gain in various applications focusing in particular on enhanced heat transfer.


title: A Turing mechanism in order to explain the patchy nature of Crohn's disease

abstract: Crohn's disease is an inflammatory bowel disease (IBD) that is not well understood. In particular, unlike other IBDs, the inflamed parts of the intestine compromise deep layers of the tissue and are not continuous but separated and distributed through the whole gastrointestinal tract, displaying a patchy inflammatory pattern. In the present paper, we introduce a toy-model which might explain the appearance of such patterns. We consider a reaction-diffusion system involving bacteria and phagocyte and prove that, under certain conditions, this system might reproduce an activator-inhibitor dynamic leading to the occurrence of Turing-type instabilities. In other words, we prove the existence of stable stationary solutions that are spatially periodic and do not vanish in time. We also propose a set of parameters for which the system exhibits such phenomena and compare it with realistic parameters found in the literature. This is the first time, as far as we know, that a Turing pattern is investigated in inflammatory models.


title: The Winfree model with non-infinitesimal phase-response curve: Ott-Antonsen theory

abstract: A novel generalization of the Winfree model of globally coupled phase oscillators, representing phase reduction under finite coupling, is studied analytically. We consider interactions through a non-infinitesimal (or finite) phase-response curve (PRC), in contrast to the infinitesimal PRC of the original model. For a family of non-infinitesimal PRCs, the global dynamics is captured by one complex-valued ordinary differential equation resorting to the Ott-Antonsen ansatz. The phase diagrams are thereupon obtained for four illustrative cases of non-infinitesimal PRC. Bistability between collective synchronization and full desynchronization is observed in all cases.


title: Four dimensional biharmonic hypersurfaces in nonzero space form have constant mean curvature

abstract: In this paper, through making careful analysis of Gauss and Codazzi equations, we prove that four dimensional biharmonic hypersurfaces in nonzero space form have constant mean curvature. Our result gives the positive answer to the conjecture proposed by Balmus-Montaldo-Oniciuc in 2008 for four dimensional hypersurfaces.


title: Longitudinal Doppler effect in de Sitter expanding universe

abstract: A new redshift formula is obtained considering the longitudinal Doppler effect in the de Sitter expanding universe where the relative geodesic motion is governed by the Lorenzian isometries of our new de Sitter relativity [I. I. Cot\u aescu, {Eur. Phys. J. C} {\bf 77} (2017) 485]. This formula combines in a non-trivial manner the well-known cosmological contribution given by Lama\' itre's law with that of the relative motion of the source with respect to a fixed observer as in de Sitter relativity instead of special relativity. Other related quantities as the dispersion relation, the propagation time of the photon and the real distance between source and observer at the moment of observation are discussed pointing out their specific features in the de Sitter relativity.


title: DICE: Dynamic Interconnections for the Cellular Ecosystem

abstract: To enable roaming of users, the cellular ecosystem integrates many entities and procedures, including specific infrastructure to connect Mobile Network Operators (MNOs), business partnerships or the use of third-party Data Clearing Houses (DCHs) for billing. Many of these rely on specifications rooted in dated and arcane practices, involving long waiting periods for financial clearing, complex billing models, and disparate mechanisms for dealing with inter-MNO disputes. In this paper, we propose a novel solution, DICE (Dynamic Interconnections for the Cellular Ecosystem), aimed at facilitating dynamic collaboration between MNOs, and sustain fluid interconnection models between the end-users and MNOs. DICE uses distributed ledger technology (DLT) to enable MNOs to interact directly, and offer customizable services to their users through the use of crypto-currencies. We leverage real-world data from a major operational MNO in Europe to support our claims, and to extract the requirements for the DICE system. We introduce the DICE protocol, and discuss real-world implementation considerations.


title: Case Survey Studies in Software Engineering Research

abstract: Background: Given the social aspects of Software Engineering (SE), in the last twenty years, researchers from the field started using research methods common in social sciences such as case study, ethnography, and grounded theory. More recently, case survey, another imported research method, has seen its increasing use in SE studies. It is based on existing case studies reported in the literature and intends to harness the generalizability of survey and the depth of case study. However, little is known on how case survey has been applied in SE research, let alone guidelines on how to employ it properly. Aims: This article aims to provide a better understanding of how case survey has been applied in Software Engineering research. Method: To address this knowledge gap, we performed a systematic mapping study and analyzed 12 Software Engineering studies that used the case survey method. Results: Our findings show that these studies presented a heterogeneous understanding of the approach ranging from secondary studies to primary inquiries focused on a large number of instances of a research phenomenon. They have not applied the case survey method consistently as defined in the seminal methodological papers. Conclusions: We conclude that a set of clearly defined guidelines are needed on how to use case survey in SE research, to ensure the quality of the studies employing this approach and to provide a set of clearly defined criteria to evaluate such work.


title: Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids

abstract: The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models. We build grids of models that sample randomly a large space of physical conditions (unobservable quantities such as gas density, temperature, etc.) and compute the corresponding observables (line intensities, column densities) and the ionization fraction. We estimate the predictive power of each potential tracer by training a Random Forest model to predict the ionization fraction from that tracer, based on these model grids. In both translucent medium and cold dense medium conditions, several observable tracers with very good predictive power for the ionization fraction are found. Several tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.). (abridged)


title: The Complexity of the Distributed Constraint Satisfaction Problem

abstract: We study the complexity of the Distributed Constraint Satisfaction Problem (DCSP) on a synchronous, anonymous network from a theoretical standpoint. In this setting, variables and constraints are controlled by agents which communicate with each other by sending messages through fixed communication channels. Our results endorse the well-known fact from classical CSPs that the complexity of fixed-template computational problems depends on the template's invariance under certain operations. Specifically, we show that DCSP($\Gamma$) is polynomial-time tractable if and only if $\Gamma$ is invariant under symmetric polymorphisms of all arities. Otherwise, there are no algorithms that solve DCSP($\Gamma$) in finite time. We also show that the same condition holds for the search variant of DCSP. Collaterally, our results unveil a feature of the processes' neighbourhood in a distributed network, its iterated degree, which plays a major role in the analysis. We explore this notion establishing a tight connection with the basic linear programming relaxation of a CSP.


title: SparseTrain: Exploiting Dataflow Sparsity for Efficient Convolutional Neural Networks Training

abstract: Training Convolutional Neural Networks (CNNs) usually requires a large number of computational resources. In this paper, \textit{SparseTrain} is proposed to accelerate CNN training by fully exploiting the sparsity. It mainly involves three levels of innovations: activation gradients pruning algorithm, sparse training dataflow, and accelerator architecture. By applying a stochastic pruning algorithm on each layer, the sparsity of back-propagation gradients can be increased dramatically without degrading training accuracy and convergence rate. Moreover, to utilize both \textit{natural sparsity} (resulted from ReLU or Pooling layers) and \textit{artificial sparsity} (brought by pruning algorithm), a sparse-aware architecture is proposed for training acceleration. This architecture supports forward and back-propagation of CNN by adopting 1-Dimensional convolution dataflow. We have built %a simple compiler to map CNNs topology onto \textit{SparseTrain}, and a cycle-accurate architecture simulator to evaluate the performance and efficiency based on the synthesized design with $14nm$ FinFET technologies. Evaluation results on AlexNet/ResNet show that \textit{SparseTrain} could achieve about $2.7 \times$ speedup and $2.2 \times$ energy efficiency improvement on average compared with the original training process.


title: Analyzing the dominant SARS-CoV-2 transmission routes towards an ab initio SEIR model

abstract: Identifying the relative importance of the different transmission routes of the SARS-CoV-2 virus is an urgent research priority. To that end, the different transmission routes, and their role in determining the evolution of the Covid-19 pandemic are analyzed in this work. Probability of infection caused by inhaling virus-laden droplets (initial, ejection diameters between $0.5-750\mu m$) and the corresponding desiccated nuclei that mostly encapsulate the virions post droplet evaporation, are individually calculated. At typical, air-conditioned yet quiescent indoor space, for average viral loading, cough droplets of initial diameter between $10-50 \mu m$ have the highest infection probability. However, by the time they are inhaled, the diameters reduce to about $1/6^{th}$ of their initial diameters. While the initially near unity infection probability due to droplets rapidly decays within the first $25s$, the small yet persistent infection probability of desiccated nuclei decays appreciably only by $\mathcal{O} (1000s)$, assuming the virus sustains equally well within the dried droplet nuclei as in the droplets. Combined with molecular collision theory adapted to calculate frequency of contact between the susceptible population and the droplet/nuclei cloud, infection rate constants are derived ab-initio, leading to a SEIR model applicable for any respiratory event - vector combination. Viral load, minimum infectious dose, sensitivity of the virus half-life to the phase of its vector and dilution of the respiratory jet/puff by the entraining air are shown to mechanistically determine specific physical modes of transmission and variation in the basic reproduction number $\mathcal{R}_0$, from first principle calculations.


title: 3D diffractive imaging of nanoparticle ensembles using an X-ray laser

abstract: We report the 3D structure determination of gold nanoparticles (AuNPs) by X-ray single particle imaging (SPI). Around 10 million diffraction patterns from gold nanoparticles were measured in less than 100 hours of beam time, more than 100 times the amount of data in any single prior SPI experiment, using the new capabilities of the European X-ray free electron laser which allow measurements of 1500 frames per second. A classification and structural sorting method was developed to disentangle the heterogeneity of the particles and to obtain a resolution of better than 3 nm. With these new experimental and analytical developments, we have entered a new era for the SPI method and the path towards close-to-atomic resolution imaging of biomolecules is apparent.


title: Maximal Interaction of Electromagnetic Radiation with Corona-Virions

abstract: Absorption and scattering of the impinging electromagnetic waves are the two fundamental operations describing the energy exchange of any, organic or inorganic, particle with its environment. In the case of virion cells, substantial extinction power, counting both absorbing and scattering effects, is a prerequisite for performing a variety of coupling actions against the viral particles and, thus, a highly sought-after feature. By considering realistic dispersion for the dielectric permittivity of proteins and a core-shell modeling allowing for rigorous formulation via Mie theory, we report optical extinction resonances for corona-virions at mid-infrared range that are not significantly perturbed by changes in the objects size or the background host. Our findings indicate the optimal regime for interaction of photonic radiation with viral particles and may assist towards the development of equipment for thermal damage, disintegration or neutralization of coronavirus cells.

categories: physics.app-ph physics.med-ph physics.optics physics.soc-ph


title: Interlacing properties of system-poles, system-zeros and spectral-zeros in MIMO systems

abstract: SISO passive systems with just one type of memory/storage element (either only inductive or only capacitative) are known to have real poles and zeros, and further, with the zeros interlacing poles (ZIP). Due to a variety of definitions of the notion of a system zero, and due to other reasons described in the paper, results involving ZIP have not been extended to MIMO systems. This paper formulates conditions under which MIMO systems too have interlaced poles and zeros. This paper next focusses on the notion of a `spectral zero' of a system, which has been well-studied in various contexts: for example, spectral factorization, optimal charging/discharging of a dissipative system, and even model order reduction. We formulate conditions under which the spectral zeros of a MIMO system are real, and further, conditions that guarantee that the system-zeros, spectral zeros, and the poles are all interlaced. The techniques used in the proofs involve new results in Algebraic Riccati equations (ARE) and Hamiltonian matrices, and these results help in formulating new notions of positive-real balancing, and inter-relations with the existing notion of positive-real balancing; we also relate the positive-real singular values with the eigenvalues of the extremal ARE solutions in the proposed `quasi-balanced' forms.


title: Shadows around the q-metric

abstract: One crucial problem in relativistic astrophysics is that of the nature of black hole candidates. It is usually assumed that astrophysical black holes are described by the Schwarzschild or Kerr space-times; however, there is no direct evidence to assert this. Moreover, there are various solutions in general relativity that can be alternatives to black holes, usually called black hole mimickers. In this work, we study the shadow produced by a compact object described by the q-metric, which is the simplest static and axially symmetric solution of Einstein equations with a non-vanishing quadrupole moment. This particular spacetime has the property of containing an independent parameter $q$, which is related to the compact object deformation. The solution corresponds to naked singularities for some specific values of this parameter. Additionally, we analyze the eigenvalues of the Riemann tensor using the $SO(3,C)$ representation, which allows us to find, in an invariant way, regions where there may be repulsive effects. Furthermore, we numerically solve the motion equations to show the shadow, the Einstein ring, and the gravitational lensing to establish a possible signature of such repulsive effects. We found that as $q$ is smaller, the Einstein ring decreases, but the shape is the same as the Schwarzschild black hole case. However, for values of $q$ lower or equal than $-0.5$, repulsive gravitational effects appear in the gravitational lensing close to the compact object, where a strong dependence of the system to the initial conditions seems to take place.


title: Continuous Fuzzy Transform as Integral Operator

abstract: The Fuzzy transform is ubiquitous in different research fields and applications, such as image and data compression, data mining, knowledge discovery, and the analysis of linguistic expressions. As a generalisation of the Fuzzy transform, we introduce the continuous Fuzzy transform and its inverse, as an integral operator induced by a kernel function. Through the relation between membership functions and integral kernels, we show that the main properties (e.g., continuity, symmetry) of the membership functions are inherited by the continuous Fuzzy transform. Then, the relation between the continuous Fuzzy transform and integral operators is used to introduce a data-driven Fuzzy transform, which encodes intrinsic information (e.g., structure, geometry, sampling density) about the input data. In this way, we avoid coarse fuzzy partitions, which group data into large clusters that do not adapt to their local behaviour, or a too dense fuzzy partition, which generally has cells that are not covered by the data, thus being redundant and resulting in a higher computational cost. To this end, the data-driven membership functions are defined by properly filtering the spectrum of the Laplace-Beltrami operator associated with the input data. Finally, we introduce the space of continuous Fuzzy transforms, which is useful for the comparison of different continuous Fuzzy transforms and for their efficient computation.


title: The role of the multiple excitation manifold in a driven quantum simulator of an antenna complex

abstract: Biomolecular light-harvesting antennas operate as nanoscale devices in a regime where the coherent interactions of individual light, matter and vibrational quanta are non-perturbatively strong. The complex behaviour arising from this could, if fully understood, be exploited for myriad energy applications. However, non-perturbative dynamics are computationally challenging to simulate, and experiments on biomaterials explore very limited regions of the non-perturbative parameter space. So-called `quantum simulators' of light-harvesting models could provide a solution to this problem, and here we employ the hierarchical equations of motion technique to investigate recent superconducting experiments of Poto{\v{c}}nik $\it{et}$ $\it{al.}$ (Nat. Com. 9, 904 (2018)) used to explore excitonic energy capture. By explicitly including the role of optical driving fields, non-perturbative dephasing noise and the full multi-excitation Hilbert space of a three-qubit quantum circuit, we predict the measureable impact of these factors on transfer efficiency. By analysis of the eigenspectrum of the network, we uncover a structure of energy levels that allows the network to exploit optical `dark' states and excited state absorption for energy transfer. We also confirm that time-resolvable coherent oscillations could be experimentally observed, even under strong, non-additive action of the driving and optical fields.


title: The non-isentropic Einstein-Euler system written in a symmetric hyperbolic form

abstract: We cast the non--isentropic Einstein--Euler system into a symmetric hyperbolic form. Such systems are very suited to treat initial value problems of hyperbolic type. We obtain this form by using the pressure $p$ and not the density $\rho$ as a variable. However, the system becomes degenerate when the pressure $p$ approaches zero, and in these cases we regularise the system by replacing the pressure with an appropriate new matter variable, the Makino variable.


title: Linear and logarithmic entanglement production in an interacting chaotic system

abstract: We investigate entanglement growth for a pair of coupled kicked rotors. For weak coupling, the growth of the entanglement entropy is found to be initially linear followed by a logarithmic growth. We calculate analytically the time after which the entanglement entropy changes its profile, and a good agreement with the numerical result is found. We further show that the different regimes of entanglement growth are associated with different rates of energy growth displayed by a rotor. At a large time, energy grows diffusively, which is preceded by an intermediate dynamical localization. The time-span of intermediate dynamical localization decreases with increasing coupling strength. We argue that the observed diffusive energy growth is the result of one rotor acting as an environment to the other which destroys the coherence. We show that the decay of the coherence is initially exponential followed by a power-law.


title: Alternating proximal-gradient steps for (stochastic) nonconvex-concave minimax problems

abstract: Minimax problems of the form $\min_x \max_y \Psi(x,y)$ have attracted increased interest largely due to advances in machine learning, in particular generative adversarial networks. These are typically trained using variants of stochastic gradient descent for the two players. Although convex-concave problems are well understood with many efficient solution methods to choose from, theoretical guarantees outside of this setting are sometimes lacking even for the simplest algorithms. In particular, this is the case for alternating gradient descent ascent, where the two agents take turns updating their strategies. To partially close this gap in the literature we prove a novel global convergence rate for the stochastic version of this method for finding a critical point of $g(\cdot) := \max_y \Psi(\cdot,y)$ in a setting which is not convex-concave.


title: X-ray quasi-periodic oscillations and estimates of the mass and the spin of the neutron star in 4U 1728-34

abstract: Three of the averaged light curves from the X-ray spectrum of 4U 1728-34 contain triads of quasi-periodic oscillations: two kHz, an upper and a lower one, and one low-frequency. In the present work, with the assumption that the triads occur simultaneously and with the application of the relativistic precession model the mass and the spin of the neutron star in the atoll source 4U 1728-34 is obtained. The consistency of the three independent estimates of the mass and the spin is also tested. According to the obtained results, the studied neutron star has relatively high mass and moderate spin.


title: Finite-size scaling analysis of two-dimensional deformed Affleck-Kennedy-Lieb-Tasaki states

abstract: Using tensor network methods, we perform finite-size scaling analysis to study the parameter-induced phase transitions of two-dimensional deformed Affleck-Kennedy-Lieb-Tasaki states. We use higher-order tensor renormalization group method to evaluate the moments and the correlations. Then, the critical point and critical exponents are determined simultaneously by collapsing the data. Alternatively, the crossing points of the dimensionless ratios are used to determine the critical point, and the scaling at the critical point is used to determine the critical exponents. For the transition between the disordered AKLT phase and the ferromagnetic ordered phase, we demonstrate that both the critical point and the exponents can be determined accurately. Furthermore, the values of the exponents confirm that the AKLT-FM transition belongs to the 2D Ising universality class. We also investigate the Berezinskii-Kosterlitz-Thouless transition from the AKLT phase to the critical XY phase. In this case we show that the critical point can be located by the crossing point of the correlation ratio.


title: Transport under advective trapping

abstract: Advective trapping occurs when solute enters low velocity zones in heterogeneous porous media. Classical local modeling approaches combine the impact of slow advection and diffusion into a hydrodynamic dispersion coefficient and many temporally non-local approaches lump these mechanisms into a single memory function. This joint treatment makes parameterization difficult and thus prediction of large scale transport a challenge. Here we investigate the mechanisms of advective trapping and their impact on transport in media composed of a high conductivity background and isolated low permeability inclusions. Breakthrough curves show that effective transport changes from a streamtube-like behavior to genuine random trapping as the degree of disorder of the inclusion arrangement increases. We upscale this behavior using a Lagrangian view point, in which idealized solute particles transition over a fixed distance at random advection times combined with Poissonian advective trapping events. We discuss the mathematical formulation of the upscaled model in the continuous time random walk and mobile-immobile mass transfer frameworks, and derive a model for large scale solute non-Fickian dispersion. These findings give new insight into transport in highly heterogeneous media.


title: Statistical Bootstrapping for Uncertainty Estimation in Off-Policy Evaluation

abstract: In reinforcement learning, it is typical to use the empirically observed transitions and rewards to estimate the value of a policy via either model-based or Q-fitting approaches. Although straightforward, these techniques in general yield biased estimates of the true value of the policy. In this work, we investigate the potential for statistical bootstrapping to be used as a way to take these biased estimates and produce calibrated confidence intervals for the true value of the policy. We identify conditions - specifically, sufficient data size and sufficient coverage - under which statistical bootstrapping in this setting is guaranteed to yield correct confidence intervals. In practical situations, these conditions often do not hold, and so we discuss and propose mechanisms that can be employed to mitigate their effects. We evaluate our proposed method and show that it can yield accurate confidence intervals in a variety of conditions, including challenging continuous control environments and small data regimes.


title: Analytic continuation of the self-energy via Machine Learning techniques

abstract: We develop a novel analytic continuation method for self-energies on the Matsubara domain as computed by quantum Monte Carlo simulations within dynamical mean field theory (QMC+DMFT). Unlike a maximum entropy (maxEn) procedure employed for the last thirty years, our approach is based on a machine learning (ML) technique in combination with the iterative perturbative theory impurity solver of the dynamical mean field theory self-consistent process (IPT+DMFT). The input and output training datasets for ML are simultaneously obtained from IPT+DMFT calculations on Matsubara and real frequency domains, respectively. The QMC+DMFT self-energy on real frequencies is determined from the -- usually noisy -- input QMC+DMFT self-energy on the Matsubara domain and the trained ML kernel. Our approach is free from both, bias of ML training datasets and from fitting parameters present in the maxEn method. We demonstrate the efficiency of the method on the testbed frustrated Hubbard model on the square lattice.


title: Non-backtracking Spectrum: Unitary Eigenvalues and Diagonalizability

abstract: Much effort has been spent on characterizing the spectrum of the non-backtracking matrix of certain classes of graphs, with special emphasis on the leading eigenvalue or the second eigenvector. Much less attention has been paid to the eigenvalues of small magnitude; here, we fully characterize the eigenvalues with magnitude equal to one. We relate the multiplicities of such eigenvalues to the existence of specific subgraphs. We formulate a conjecture on necessary and sufficient conditions for the diagonalizability of the non backtracking matrix. As an application, we establish an interlacing-type result for the Perron eigenvalue.


title: Non-local SPDE limits of spatially-correlated-noise driven spin systems derived to sample a canonical distribution

abstract: We study the macroscopic behavior of a stochastic spin ensemble driven by a discrete Markov jump process motivated by the Metropolis-Hastings algorithm where the proposal is made with spatially correlated (colored) noise, and hence fails to be symmetric. However, we demonstrate a scenario where the failure of proposal symmetry is a higher order effect. Hence, from these microscopic dynamics we derive as a limit as the proposal size goes to zero and the number of spins to infinity, a non-local stochastic version of the harmonic map heat flow (or overdamped Landau-Lipshitz equation). The equation is both mathematically well-posed and samples the canonical/Gibbs distribution related to the kinetic energy. The failure of proposal symmetry due to interaction between the confining geometry of the spin system and the colored noise is in contrast to the uncorrelated, white-noise, driven system. Specifically, the choice of projection of the noise to conserve the magnitude of the spins is crucial to maintaining the proper equilibrium distribution. Numerical simulations are included to verify convergence properties and demonstrate the dynamics.