title: Asymptotic low-temperature behavior of two-dimensional RP$^{N-1}$ models

abstract: We investigate the low-temperature behavior of two-dimensional (2D) RP$^{N-1}$ models, characterized by a global O($N$) symmetry and a local ${\mathbb Z}_2$ symmetry. For $N=3$ we perform large-scale simulations of four different 2D lattice models: two standard lattice models and two different constrained models. We also consider a constrained mixed O(3)-RP$^2$ model for values of the parameters such that vector correlations are always disordered. We find that all these models show the same finite-size scaling (FSS) behavior, and therefore belong to the same universality class. However, these FSS curves differ from those computed in the 2D O(3) $\sigma$ model, suggesting the existence of a distinct 2D RP$^2$ universality class. We also performed simulations for $N=4$, and the corresponding FSS results also support the existence of an RP$^3$ universality class, different from the O(4) one.


title: Successful implementation of discrete event simulation: the case of an Italian emergency department

abstract: This paper focuses on the study of a practical management problem faced by a healthcare {\it emergency department} (ED) located in the north of Italy. The objective of our study was to propose organisational changes in the selected ED, which admits approximately 7000 patients per month, aiming at improving key performance indicators related to patient satisfaction, such as the waiting time. Our study is based on a design thinking process that adopts a {\it discrete event simulation} (DES) model as the main tool for proposing changes. We used the DES model to propose and evaluate the impact of different improving scenarios. The model is based on historical data, on the observation of the current ED situation, and information obtained from the ED staff. The results obtained by the DES model have been compared with those related to the existing ED setting, and then validated by the ED managers. Based on the results we obtained, one of the tested scenarios was selected by the ED for implementation.


title: Hybrid Session-based News Recommendation using Recurrent Neural Networks

abstract: We describe a hybrid meta-architecture -- the CHAMELEON -- for session-based news recommendation that is able to leverage a variety of information types using Recurrent Neural Networks. We evaluated our approach on two public datasets, using a temporal evaluation protocol that simulates the dynamics of a news portal in a realistic way. Our results confirm the benefits of modeling the sequence of session clicks with RNNs and leveraging side information about users and articles, resulting in significantly higher recommendation accuracy and catalog coverage than other session-based algorithms.


title: Mixed Integer Linear Programming and Constraint Programming Models for the Online Printing Shop Scheduling Problem

abstract: In this work, the online printing shop scheduling problem is considered. This challenging real problem, that appears in the nowadays printing industry, can be seen as a flexible job shop scheduling problem with sequence flexibility in which precedence constraints among operations of a job are given by an arbitrary directed acyclic graph. In addition, several complicating particularities such as periods of unavailability of the machines, resumable operations, sequence-dependent setup times, partial overlapping among operations with precedence constraints, release times, and fixed operations are also present in the addressed problem. In the present work, mixed integer linear programming and constraint programming models for the minimization of the makespan are presented. Modeling the problem is twofold. On the one hand, the problem is precisely defined. On the other hand, the capabilities and limitations of a commercial software for solving the models are analyzed. Extensive numerical experiments with small-, medium-, and large-sized instances are presented. Numerical experiments show that the commercial solver is able to optimally solve only a fraction of the small-sized instances when considering the mixed integer linear programming model; while all small-sized and a fraction of the medium-sized instances are optimally solved when considering the constraint programming formulation of the problem. Moreover, the commercial solver is able to deliver feasible solutions for the large-sized instances that are of the size of the instances that appear in practice.


title: DCNNs: A Transfer Learning comparison of Full Weapon Family threat detection for Dual-Energy X-Ray Baggage Imagery

abstract: Recent advancements in Convolutional Neural Networks have yielded super-human levels of performance in image recognition tasks [13, 25]; however, with increasing volumes of parcels crossing UK borders each year, classification of threats becomes integral to the smooth operation of UK borders. In this work we propose the first pipeline to effectively process Dual-Energy X-Ray scanner output, and perform classification capable of distinguishing between firearm families (Assault Rifle, Revolver, Self-Loading Pistol,Shotgun, and Sub-Machine Gun) from this output. With this pipeline we compare re-cent Convolutional Neural Network architectures against the X-Ray baggage domain via Transfer Learning and show ResNet50 to be most suitable to classification - outlining a number of considerations for operational success within the domain.


title: Four-dimensional complete gradient shrinking Ricci solitons

abstract: In this article, we study four-dimensional complete gradient shrinking Ricci solitons. We prove that a four-dimensional complete gradient shrinking Ricci soliton satisfying a pointwise condition involving either the self-dual or anti-self-dual part of the Weyl tensor is either Einstein, or a finite quotient of either the Gaussian shrinking soliton $\Bbb{R}^4,$ or $\Bbb{S}^{3}\times\Bbb{R}$, or $\Bbb{S}^{2}\times\Bbb{R}^{2}.$ In addition, we provide some curvature estimates for four-dimensional complete gradient Ricci solitons assuming that its scalar curvature is suitable bounded by the potential function.


title: Lightweight Online Noise Reduction on Embedded Devices using Hierarchical Recurrent Neural Networks

abstract: Deep-learning based noise reduction algorithms have proven their success especially for non-stationary noises, which makes it desirable to also use them for embedded devices like hearing aids (HAs). This, however, is currently not possible with state-of-the-art methods. They either require a lot of parameters and computational power and thus are only feasible using modern CPUs. Or they are not suitable for online processing, which requires constraints like low-latency by the filter bank and the algorithm itself. In this work, we propose a mask-based noise reduction approach. Using hierarchical recurrent neural networks, we are able to drastically reduce the number of neurons per layer while including temporal context via hierarchical connections. This allows us to optimize our model towards a minimum number of parameters and floating-point operations (FLOPs), while preserving noise reduction quality compared to previous work. Our smallest network contains only 5k parameters, which makes this algorithm applicable on embedded devices. We evaluate our model on a mixture of EUROM and a real-world noise database and report objective metrics on unseen noise.


title: Diameter bounds for degenerating Calabi-Yau metrics

abstract: We obtain sharp upper and lower bounds for the diameter of Ricci-flat Kahler metrics on polarized Calabi-Yau degeneration families, as conjectured by Kontsevich-Soibelman.


title: An entropy structure preserving space-time formulation for cross-diffusion systems: Analysis and Galerkin discretization

abstract: Cross-diffusion systems are systems of nonlinear parabolic partial differential equations that are used to describe dynamical processes in several application, including chemical concentrations and cell biology. We present a space-time approach to the proof of existence of bounded weak solutions of cross-diffusion systems, making use of the system entropy to examine long-term behavior and to show that the solution is nonnegative, even when a maximum principle is not available. This approach naturally gives rise to a novel space-time Galerkin method for the numerical approximation of cross-diffusion systems that conserves their entropy structure. We prove existence and convergence of the discrete solutions, and present numerical results for the porous medium, the Fisher-KPP, and the Maxwell-Stefan problem.


title: Normalizing Flows Across Dimensions

abstract: Real-world data with underlying structure, such as pictures of faces, are hypothesized to lie on a low-dimensional manifold. This manifold hypothesis has motivated state-of-the-art generative algorithms that learn low-dimensional data representations. Unfortunately, a popular generative model, normalizing flows, cannot take advantage of this. Normalizing flows are based on successive variable transformations that are, by design, incapable of learning lower-dimensional representations. In this paper we introduce noisy injective flows (NIF), a generalization of normalizing flows that can go across dimensions. NIF explicitly map the latent space to a learnable manifold in a high-dimensional data space using injective transformations. We further employ an additive noise model to account for deviations from the manifold and identify a stochastic inverse of the generative process. Empirically, we demonstrate that a simple application of our method to existing flow architectures can significantly improve sample quality and yield separable data embeddings.


title: Domain Adaptation for Semantic Parsing

abstract: Recently, semantic parsing has attracted much attention in the community. Although many neural modeling efforts have greatly improved the performance, it still suffers from the data scarcity issue. In this paper, we propose a novel semantic parser for domain adaptation, where we have much fewer annotated data in the target domain compared to the source domain. Our semantic parser benefits from a two-stage coarse-to-fine framework, thus can provide different and accurate treatments for the two stages, i.e., focusing on domain invariant and domain specific information, respectively. In the coarse stage, our novel domain discrimination component and domain relevance attention encourage the model to learn transferable domain general structures. In the fine stage, the model is guided to concentrate on domain related details. Experiments on a benchmark dataset show that our method consistently outperforms several popular domain adaptation strategies. Additionally, we show that our model can well exploit limited target data to capture the difference between the source and target domain, even when the target domain has far fewer training instances.


title: Wavelet Augmented Regression Profiling (WARP): improved long-term estimation of travel time series with recurrent congestion

abstract: Reliable estimates of typical travel times allow road users to forward plan journeys to minimise travel time, potentially increasing overall system efficiency. On busy highways, however, congestion events can cause large, short-term spikes in travel time. These spikes make direct forecasting of travel time using standard time series models difficult on the timescales of hours to days that are relevant to forward planning. The problem is that some such spikes are caused by unpredictable incidents and should be filtered out, whereas others are caused by recurrent peaks in demand and should be factored into estimates. Here we present the Wavelet Augmented Regression Profiling (WARP) method for long-term estimation of typical travel times. WARP linearly decomposes historical time series of travel times into two components: background and spikes. It then further separates the spikes into contributions from recurrent and residual congestion. This is achieved using a combination of wavelet transforms, spectral filtering and locally weighted regression. The background and recurrent congestion contributions are then used to estimate typical travel times with horizon of one week in an accurate and computationally inexpensive manner. We train and test WARP on the M6 and M11 motorways in the United Kingdom using 12 weeks of link level travel time data obtained from the UK's National Traffic Information Service (NTIS). In out-of-sample validation tests, WARP compares favourably to estimates produced by a simple segmentation method and to the estimates published by NTIS.


title: Reduction From Non-Unique Games To Boolean Unique Games

abstract: We reduce the problem of proving a "Boolean Unique Games Conjecture" (with gap 1-delta vs. 1-C*delta, for any C> 1, and sufficiently small delta>0) to the problem of proving a PCP Theorem for a certain non-unique game. In a previous work, Khot and Moshkovitz suggested an inefficient candidate reduction (i.e., without a proof of soundness). The current work is the first to provide an efficient reduction along with a proof of soundness. The non-unique game we reduce from is similar to non-unique games for which PCP theorems are known. Our proof relies on a new concentration theorem for functions in Gaussian space that are restricted to a random hyperplane. We bound the typical Euclidean distance between the low degree part of the restriction of the function to the hyperplane and the restriction to the hyperplane of the low degree part of the function.


title: New examples of shrinking Laplacian solitons

abstract: We give a one-parameter family of examples of shrinking Laplacian solitons, which are the second known solutions to the closed $G_2$-Laplacian flow with a finite-time singularity. The torsion forms and the Laplacian and Ricci operators of a large family of $G_2$-structures on different Lie groups are also studied. We apply these formulas to prove that, under a suitable extra condition, there is no closed eigenform for the Laplacian on such family.


title: Virtual mutations of weighted surface algebras

abstract: The finite-dimensional symmetric algebras over an algebraically closed field, based on surface triangulations, motivated by the theory of cluster algebras, have been extensively investigated and applied. In particular, the weighted surface algebras and their deformations were introduced and studied in [16]-[20], and it was shown that all these algebras, except few singular cases, are symmetric tame periodic algebras of period $4$. In this article, using the general form of a weighted surface algebra from [19], we introduce and study so called virtual mutations of weighted surface algebras, which constitute a new large class of symmetric tame periodic algebras of period $4$. We prove that all these algebras are derived equivalent but not isomorphic to weighted surface algebras. We associate such algebras to any triangulated surface, first taking blow-ups of a family of edges to $2$-triangle discs, and then virtual mutations of their weighted surface algebras. The results of this paper form an essential step towards a classification of all tame symmetric periodic algebras.


title: Measurement of the $t\bar{t}$ production cross-section in the lepton+jets channel at $\sqrt{s}=13\;$TeV with the ATLAS experiment

abstract: The $t\bar{t}$ production cross-section is measured in the lepton+jets channel using proton$-$proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV collected with the ATLAS detector at the LHC. The dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. Events with exactly one charged lepton and four or more jets in the final state, with at least one jet containing $b$-hadrons, are used to determine the $t\bar{t}$ production cross-section through a profile-likelihood fit. The inclusive cross-section is measured to be ${\sigma_{\text{inc}} = 830 \pm 0.4~ \text{(stat.)}\pm 36~\text{(syst.)}\pm 14~\text{(lumi.)}~\mathrm{pb}}$ with a relative uncertainty of 4.6 %. The result is consistent with theoretical calculations at next-to-next-to-leading order in perturbative QCD. The fiducial $t\bar{t}$ cross-section within the experimental acceptance is also measured.


title: CLC: Complex Linear Coding for the DNS 2020 Challenge

abstract: Complex-valued processing brought deep learning-based speech enhancement and signal extraction to a new level. Typically, the noise reduction process is based on a time-frequency (TF) mask which is applied to a noisy spectrogram. Complex masks (CM) usually outperform real-valued masks due to their ability to modify the phase. Recent work proposed to use a complex linear combination of coefficients called complex linear coding (CLC) instead of a point-wise multiplication with a mask. This allows to incorporate information from previous and optionally future time steps which results in superior performance over mask-based enhancement for certain noise conditions. In fact, the linear combination enables to model quasi-steady properties like the spectrum within a frequency band. In this work, we apply CLC to the Deep Noise Suppression (DNS) challenge and propose CLC as an alternative to traditional mask-based processing, e.g. used by the baseline. We evaluated our models using the provided test set and an additional validation set with real-world stationary and non-stationary noises. Based on the published test set, we outperform the baseline w.r.t. the scale independent signal distortion ratio (SI-SDR) by about 3dB.


title: General Relativity and the AKSZ construction

abstract: In this note the AKSZ construction is applied to the BFV description of the reduced phase space of the Einstein-Hilbert and of the Palatini--Cartan theories in every space-time dimension greater than two. In the former case one obtains a BV theory for the first-order formulation of Einstein--Hilbert theory, in the latter a BV theory for Palatini--Cartan theory with a partial implementation of the torsion-free condition already on the space of fields. All theories described here are BV versions of the same classical system on cylinders. The AKSZ implementations we present have the advantage of yielding a compatible BV-BFV description, which is the required starting point for a quantization in presence of a boundary.


title: Coconut Palm: Static and Streaming Data Series Exploration Now in your Palm

abstract: Many modern applications produce massive streams of data series and maintain them in indexes to be able to explore them through nearest neighbor search. Existing data series indexes, however, are expensive to operate as they issue many random I/Os to storage. To address this problem, we recently proposed Coconut, a new infrastructure that organizes data series based on a new sortable format. In this way, Coconut is able to leverage state-of-the-art indexing techniques that rely on sorting for the first time to build, maintain and query data series indexes using fast sequential I/Os. In this demonstration, we present Coconut Palm, a new exploration tool that allows to interactively combine different indexing techniques from within the Coconut infrastructure and to thereby seamlessly explore data series from across various scientific domains. We highlight the rich indexing design choices that Coconut opens up, and we present a new recommender tool that allows users to intelligently navigate them for both static and streaming data exploration scenarios.


title: Spiral Thermal Waves Generated by Self-Propelled Camphor Boats

abstract: Spiral thermal surface waves arising from self-propulsion of the camphor-driven objects are reported. Spiral thermal waves were registered for the dissolution and evaporation-guided self-propulsion. Soluto-capillarity is accompanied by thermo-capillarity under self-propulsion of camphor boats. The jump in the surface tension due to the soluto-capillarity is much larger than that due to the thermo-capillarity. The observed thermal effect is related to the adsorption of camphor molecules at the water/vapor interface. The observed spirals are shaped as Archimedean ones.


title: Energy currents in the field formed by waves with different frequencies

abstract: The behavior of instantaneous and averaged vectors of the Poynting vector transverse component for the resulting field formed as a superposition of waves with different frequencies and different polarizations is considered. Results of computer simulation and experimental data are obtained for the superposition of frequently differed vortices with different and the same signs of topological charges. It is shown that the captured particle rotates due to the orbital momentum and substantial increase of particle rotation speed in the case of polychromatic radiation as opposed to monochromatic is observed.


title: Microstructure of the near-wall layer of filtration-induced colloidal assembly

abstract: This paper describes an experimental study of filtration of a colloidal suspension using microfluidic devices. A suspension of micrometer-scale colloids flows through parallel slit-shaped pores at fixed pressure drop. Clogs and cakes are systematically observed at pore entrance, for variable applied pressure drop and ionic strength. Based on image analysis of the layer of colloids close to the device wall, global and local studies are performed to analyse in detail the near-wall layer microstructure. Whereas global porosity of this layer does not seem to be affected by ionic strength and applied pressure drop, a local study shows some heterogeneity: clogs are more porous at the vicinity of the pore than far away. An analysis of medium-range order using radial distribution function shows a slightly more organized state at high ionic strength. This is confirmed by a local analysis using two-dimension continuous wavelet decomposition: the typical size of crystals of colloids is larger for low ionic strength, and it increases with distance from the pores. We bring these results together in a phase diagram involving colloid-colloid repulsive interactions and fluid velocity.


title: Rationality of Equivariant Hilbert Series and Asymptotic Properties

abstract: An FI- or an OI-module $\mathbf{M}$ over a corresponding noetherian polynomial algebra $\mathbf{P}$ may be thought of as a sequence of compatible modules $\mathbf{M}_n$ over a polynomial ring $\mathbf{P}_n$ whose number of variables depends linearly on $n$. In order to study invariants of the modules $\mathbf{M}_n$ in dependence of $n$, an equivariant Hilbert series is introduced if $\mathbf{M}$ is graded. If $\mathbf{M}$ is also finitely generated, it is shown that this series is a rational function. Moreover, if this function is written in reduced form rather precise information about the irreducible factors of the denominator is obtained. This is key for applications. It follows that the Krull dimension of the modules $\mathbf{M}_n$ grows eventually linearly in $n$, whereas the multiplicity of $\mathbf{M}_n$ grows eventually exponentially in $n$. Moreover, for any fixed degree $j$, the vector space dimensions of the degree $j$ components of $\mathbf{M}_n$ grow eventually polynomially in $n$. As a consequence, any graded Betti number of $\mathbf{M}_n$ in a fixed homological degree and a fixed internal degree grows eventually polynomially in $n$. Furthermore, evidence is obtained to support a conjecture that the Castelnuovo-Mumford regularity and the projective dimension of $\mathbf{M}_n$ both grow eventually linearly in $n$. It is also shown that modules $\mathbf{M}$ whose width $n$ components $\mathbf{M}_n$ are eventually Artinian can be characterized by their equivariant Hilbert series. Using regular languages and finite automata, an algorithm for computing equivariant Hilbert series is presented.


title: Single-Shot 3D Detection of Vehicles from Monocular RGB Images via Geometry Constrained Keypoints in Real-Time

abstract: In this paper we propose a novel 3D single-shot object detection method for detecting vehicles in monocular RGB images. Our approach lifts 2D detections to 3D space by predicting additional regression and classification parameters and hence keeping the runtime close to pure 2D object detection. The additional parameters are transformed to 3D bounding box keypoints within the network under geometric constraints. Our proposed method features a full 3D description including all three angles of rotation without supervision by any labeled ground truth data for the object's orientation, as it focuses on certain keypoints within the image plane. While our approach can be combined with any modern object detection framework with only little computational overhead, we exemplify the extension of SSD for the prediction of 3D bounding boxes. We test our approach on different datasets for autonomous driving and evaluate it using the challenging KITTI 3D Object Detection as well as the novel nuScenes Object Detection benchmarks. While we achieve competitive results on both benchmarks we outperform current state-of-the-art methods in terms of speed with more than 20 FPS for all tested datasets and image resolutions.


title: Calibration of Shared Equilibria in General Sum Partially Observable Markov Games

abstract: Training multi-agent systems (MAS) to achieve realistic equilibria gives us a useful tool to understand and model real-world systems. We consider a general sum partially observable Markov game where agents of different types share a single policy network, conditioned on agent-specific information. This paper aims at i) formally understanding equilibria reached by such agents, and ii) matching emergent phenomena of such equilibria to real-world targets. Parameter sharing with decentralized execution has been introduced as an efficient way to train multiple agents using a single policy network. However, the nature of resulting equilibria reached by such agents has not been yet studied: we introduce the novel concept of Shared equilibrium as a symmetric pure Nash equilibrium of a certain Functional Form Game (FFG) and prove convergence to the latter for a certain class of games using self-play. In addition, it is important that such equilibria satisfy certain constraints so that MAS are calibrated to real world data for practical use: we solve this problem by introducing a novel dual-Reinforcement Learning based approach that fits emergent behaviors of agents in a Shared equilibrium to externally-specified targets, and apply our methods to a n-player market example. We do so by calibrating parameters governing distributions of agent types rather than individual agents, which allows both behavior differentiation among agents and coherent scaling of the shared policy network to multiple agents.


title: Classifying Network Vendors at Internet Scale

abstract: In this paper, we develop a method to create a large, labeled dataset of visible network device vendors across the Internet by mapping network-visible IP addresses to device vendors. We use Internet-wide scanning, banner grabs of network-visible devices across the IPv4 address space, and clustering techniques to assign labels to more than 160,000 devices. We subsequently probe these devices and use features extracted from the responses to train a classifier that can accurately classify device vendors. Finally, we demonstrate how this method can be used to understand broader trends across the Internet by predicting device vendors in traceroutes from CAIDA's Archipelago measurement system and subsequently examining vendor distributions across these traceroutes.


title: On the Interoperability of Decentralized Exposure Notification Systems

abstract: This report summarizes the requirements and proposes a high-level solution for interoperability across recently proposed COVID-19 exposure notification efforts. Our focus is on interoperability across exposure notification (EN) applications which are based on the decentralized Bluetooth Low Energy (BLE) protocol driven by Google/Apple Exposure Notifications API (including DP3T and similar protocols). We distinguish different interoperability use cases, such as worldwide public EN interoperability, as well as interoperability in the enterprise EN systems. This report also proposes an API and a backend implementation architecture for EN interoperability. Finally, we propose using a permissioned blockchain-based solution for managing EN backend certificates and configurations (without storing any users' data on the blockchain) for helping address EN interoperability challenges across different vendors.


title: Modelling intergalactic low ionisation metal absorption line systems near the epoch of reionization

abstract: We interpret observations of intergalactic low ionisation metal absorption systems at redshifts z $\gtrsim$5 in terms of pressure-confined clouds. We find clouds confined by the expected pressure of galactic haloes with masses $11<\log M_h/h^{-1}M_\odot<12$ provide a good description of the column density ratios between low ionisation metal absorbers. Some of the ratios, however, require extending conventional radiative transfer models of irradiated slabs to spherical (or cylindrical) clouds to allow for lines of sight passing outside the cores of the clouds. Moderate depletion of silicon onto dust grains is also indicated in some systems. The chemical abundances inferred span the range between solar and massive-star dominated stellar populations as may arise in starburst galaxies. The typical HI column densities matching the data correspond to Damped Lyman-$\alpha$ Absorbers (DLAs) or sub-DLAs, with sizes of 40 pc to 3 kpc, gas masses $3.5<\log M_c/M_\odot<8$ and metallicites $0.001-0.01Z_\odot$. Such systems continue to pose a challenge for galaxy-scale numerical simulations to reproduce.


title: Strategies to build functionals of the density, or functionals of Green's functions: what can we learn?

abstract: The many-body problem can in general not be solved exactly, and one of the most prominent approximations is to build perturbation expansions. A huge variety of expansions is possible, which differ by the quantity to be expanded, the expansion variable, the starting point, and ideas how to resum or terminate the series. Although much has been discussed and much has been done, some choices were made for historical reasons, in particular, limited computation or storage capacities. The present work aims at examining the justifications for different choices made in different contexts, by comparing ingredients of functionals based on \GF s on one side, and on the charge density on the other side. Of particular interest will be the question of how to build an optimal starting point for the approximation of non-local quantities, making use of near- or far-sightedness, and daring to consider models beyond the homogeneous electron gas. This will include the use of connector approximations. We will also discuss why it is a good idea to build functionals of the density.


title: Data Augmentation View on Graph Convolutional Network and the Proposal of Monte Carlo Graph Learning

abstract: Today, there are two major understandings for graph convolutional networks, i.e., in the spectral and spatial domain. But both lack transparency. In this work, we introduce a new understanding for it -- data augmentation, which is more transparent than the previous understandings. Inspired by it, we propose a new graph learning paradigm -- Monte Carlo Graph Learning (MCGL). The core idea of MCGL contains: (1) Data augmentation: propagate the labels of the training set through the graph structure and expand the training set; (2) Model training: use the expanded training set to train traditional classifiers. We use synthetic datasets to compare the strengths of MCGL and graph convolutional operation on clean graphs. In addition, we show that MCGL's tolerance to graph structure noise is weaker than GCN on noisy graphs (four real-world datasets). Moreover, inspired by MCGL, we re-analyze the reasons why the performance of GCN becomes worse when deepened too much: rather than the mainstream view of over-smoothing, we argue that the main reason is the graph structure noise, and experimentally verify our view. The code is available at https://github.com/DongHande/MCGL.


title: A de Finetti-type representation of joint hierarchically exchangeable arrays on directed acyclic graphs

abstract: We define joint exchangeability on arrays indexed by a vector of natural numbers with coordinates being the vertices of directed acyclic graphs (DAGs) using local isomorphisms. The notion provides a new version of exchangeability, which is a joint version of hierarchical exchangeability defined in Jung, L., Staton, Yang (2020). We also prove the existence of a generic representation by independent uniform random variables.


title: Multi-Class Uncertainty Calibration via Mutual Information Maximization-based Binning

abstract: Post-hoc multi-class calibration is a common approach for providing high-quality confidence estimates of deep neural network predictions. Recent work has shown that widely used scaling methods underestimate their calibration error, while alternative Histogram Binning (HB) methods often fail to preserve classification accuracy. When classes have small prior probabilities, HB also faces the issue of severe sample-inefficiency after the conversion into K one-vs-rest class-wise calibration problems. The goal of this paper is to resolve the identified issues of HB in order to provide calibrated confidence estimates using only a small holdout calibration dataset for bin optimization while preserving multi-class ranking accuracy. From an information-theoretic perspective, we derive the I-Max concept for binning, which maximizes the mutual information between labels and quantized logits. This concept mitigates potential loss in ranking performance due to lossy quantization, and by disentangling the optimization of bin edges and representatives allows simultaneous improvement of ranking and calibration performance. To improve the sample efficiency and estimates from a small calibration set, we propose a shared class-wise (sCW) calibration strategy, sharing one calibrator among similar classes (e.g., with similar class priors) so that the training sets of their class-wise calibration problems can be merged to train the single calibrator. The combination of sCW and I-Max binning outperforms the state of the art calibration methods on various evaluation metrics across different benchmark datasets and models, using a small calibration set (e.g., 1k samples for ImageNet).


title: A dynamical system approach to a class of radial weighted fully nonlinear equations

abstract: In this paper we study existence, nonexistence and classification of radial positive solutions of some weighted fully nonlinear equations involving Pucci extremal operators. Our results are entirely based on the analysis of the dynamics induced by an autonomous quadratic system which is obtained after a suitable transformation. This method allows to treat both regular and singular solutions in a unified way, without using energy arguments. In particular we recover known results on regular solutions for the fully nonlinear non weighted problem by alternative proofs. We also slightly improve the classification of the solutions for the extremal operator $\mathcal{M}^-$.


title: The role of swabs in modeling the COVID-19 outbreak in the most affected regions of Italy

abstract: The daily fluctuations in the released number of Covid-19 cases played a big role both at the beginning and in the most critical weeks of the outbreak, when local authorities in Italy had to decide whether to impose a lockdown and at which level. Public opinion was focused on this information as well, to understand how quickly the epidemic was spreading. When an increase/decrease was communicated, especially a large one, it was not easy to understand if it was due to a change in the epidemic evolution or if it was a fluctuation due to other reasons, such as an increase in the number of swabs or a delay in the swab processing. In this work, we propose a nonlinear asymmetric diffusion model, which includes information on the daily number of swabs, to describe daily fluctuations in the number of confirmed cases in addition to the main trend of the outbreak evolution. The class of diffusion models has been selected to develop our proposal, as it also allows estimation of the total number of confirmed cases at the end of the outbreak. The proposed model is compared to six existing models, among which the logistic and the SIRD models are used as benchmarks, in the five most affected Italian regions.


title: SCARE: Side Channel Attack on In-Memory Computing for Reverse Engineering

abstract: In-memory computing architectures provide a much needed solution to energy-efficiency barriers posed by Von-Neumann computing due to the movement of data between the processor and the memory. Functions implemented in such in-memory architectures are often proprietary and constitute confidential Intellectual Property. Our studies indicate that IMCs implemented using RRAM are susceptible to Side Channel Attack. Unlike conventional SCAs that are aimed to leak private keys from cryptographic implementations, SCARE can reveal the sensitive IP implemented within the memory. Therefore, the adversary does not need to perform invasive Reverse Engineering to unlock the functionality. We demonstrate SCARE by taking recent IMC architectures such as DCIM and MAGIC as test cases. Simulation results indicate that AND, OR, and NOR gates (building blocks of complex functions) yield distinct power and timing signatures based on the number of inputs making them vulnerable to SCA. Although process variations can obfuscate the signatures due to significant overlap, we show that the adversary can use statistical modeling and analysis to identify the structure of the implemented function. SCARE can find the implemented IP by testing a limited number of patterns. For example, the proposed technique reduces the number of patterns by 64% compared to a brute force attack for a+bc function. Additionally, analysis shows improvement in SCAREs detection model due to adversarial change in supply voltage for both DCIM and MAGIC. We also propose countermeasures such as redundant inputs and expansion of literals. Redundant inputs can mask the IP with 25% area and 20% power overhead. However, functions can be found by greater RE effort. Expansion of literals incurs 36% power overhead. However, it imposes brute force search by the adversary for which the RE effort increases by 3.04X.


title: Compensating for visibility artefacts in photoacoustic imaging with a deep learning approach providing prediction uncertainties

abstract: Conventional photoacoustic imaging may suffer from the limited view and bandwidth of ultrasound transducers. A deep learning approach is proposed to handle these problems and is demonstrated both in simulations and in experiments on a multi-scale model of leaf skeleton. We employed an experimental approach to build the training and the test sets using photographs of the samples as ground truth images. Reconstructions produced by the neural network show a greatly improved image quality as compared to conventional approaches. In addition, this work aimed at quantifying the reliability of the neural network predictions. To achieve this, the dropout Monte-Carlo procedure is applied to estimate a pixel-wise degree of confidence on each predicted picture. Last, we address the possibility to use transfer learning with simulated data in order to drastically limit the size of the experimental dataset.


title: $L$-series values for sextic twists of elliptic curves over $\mathbb{Q}[\sqrt{-3}]$

abstract: We prove a new formula for the central value of the $L$-function $L(E_{D, \alpha}, 1)$ corresponding to the family of sextic twists over $\mathbb{Q}[\sqrt{-3}]$ of elliptic curves $E_{D, \alpha}: y^2=x^3+16D^2\alpha^3$ for $D$ an integer and $\alpha \in \mathbb{Q}[\sqrt{-3}]$. The formula generalizes the result of cubic twists over $\mathbb{Q}$ of Rodriguez-Villegas and Zagier for a prime $D \equiv 1 (9)$ and of Rosu for general $D$. For $\alpha$ prime and all integers $D$, we also show that the expected value from the Birch and Swinnerton-Dyer conjecture of the order of the Tate-Shafarevich group is an integer square in certain cases, and an integer square up to a factor $2^{2a}3^{2b}$ in general.


title: RASSINE: Interactive tool for normalising stellar spectra I. Description and performance of the code

abstract: Aims: We provide an open-source code allowing an easy, intuitive, and robust normalisation of spectra. Methods: We developed RASSINE, a Python code for normalising merged 1D spectra through the concepts of convex hulls. The code uses six parameters that can be easily fine-tuned. The code also provides a complete user-friendly interactive interface, including graphical feedback, that helps the user to choose the parameters as easily as possible. To facilitate the normalisation even further, RASSINE can provide a first guess for the parameters that are derived directly from the merged 1D spectrum based on previously performed calibrations. Results: For HARPS spectra of the Sun that were obtained with the HELIOS solar telescope, a continuum accuracy of 0.20% on line depth can be reached after normalisation with RASSINE. This is three times better than with the commonly used method of polynomial fitting. For HARPS spectra of $\alpha$ Cen B, a continuum accuracy of 2.0% is reached. This rather poor accuracy is mainly due to molecular band absorption and the high density of spectral lines in the bluest part of the merged 1D spectrum. When wavelengths shorter than 4500 \AA are excluded, the continuum accuracy improves by up to 1.2%. The line-depth precision on individual spectrum normalisation is estimated to be 0.15%, which can be reduced to the photon-noise limit (0.10%) when a time series of spectra is given as input for RASSINE. Conclusions: With a continuum accuracy higher than the polynomial fitting method and a line-depth precision compatible with photon noise, RASSINE is a tool that can find applications in numerous cases, for example stellar parameter determination, transmission spectroscopy of exoplanet atmospheres, or activity-sensitive line detection.


title: Bootstrapping $\ell_p$-Statistics in High Dimensions

abstract: This paper considers a new bootstrap procedure to estimate the distribution of high-dimensional $\ell_p$-statistics, i.e. the $\ell_p$-norms of the sum of $n$ independent $d$-dimensional random vectors with $d \gg n$ and $p \in [1, \infty]$. We provide a non-asymptotic characterization of the sampling distribution of $\ell_p$-statistics based on Gaussian approximation and show that the bootstrap procedure is consistent in the Kolmogorov-Smirnov distance under mild conditions on the covariance structure of the data. As an application of the general theory we propose a bootstrap hypothesis test for simultaneous inference on high-dimensional mean vectors. We establish its asymptotic correctness and consistency under high-dimensional alternatives, and discuss the power of the test as well as the size of associated confidence sets. We illustrate the bootstrap and testing procedure numerically on simulated data.


title: A study on quantum gases: bosons in optical lattices and the one-dimensional interacting Bose gas

abstract: Bosonic atoms confined in optical lattices are described by the Bose-Hubbard model and can exist in two different phases, Mott insulator or superfluid, depending on the strength of the system parameters. In the vicinity of the phase boundary, there are degeneracies that occur between every two adjacent Mott lobes. Because of this, nondegenerate perturbation theory fails to give meaningful results for the condensate density: it predicts a phase transition in a point of the phase diagram where no transition occurs. Motivated by this, we develop two different degenerate perturbative methods to solve the degeneracy-related problems. Moreover, we study the one-dimensional repulsively interacting Bose gas under harmonic confinement, with special attention to the asymptotic behavior of the momentum distribution, which is a universal $k^{-4}$ decay characterized by the Tan's contact. The latter constitutes a direct signature of the short-range correlations in such an interacting system and provides valuable insights about the role of the interparticle interactions. We investigate the system constituted of $N$ interacting particles in the strongly interacting limit. In such a regime, the strong interparticle interaction makes the bosons behave similarly to the ideal Fermi gas. Because of the difficulty in analytically solving the system for $N$ particles at finite interaction, the Tonks-Girardeau regime provides a favorable scenario to probe the contact. Therefore, we are able to provide an analytical formula for the Tan's contact. Furthermore, we analyze the scaling properties of the Tan's contact in terms of $N$ in the high-temperature regime as well as in the strongly interacting regime. Finally, we compare our analytical calculations of the Tan's contact to quantum Monte Carlo simulations and discuss some fundamental differences between the canonical and the grand-canonical ensembles.

categories: cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph physics.optics quant-ph


title: A chromosome map to unveil stellar populations with different magnesium abundances. The case of Omega Centauri

abstract: Historically, photometry has been largely used to identify stellar populations (MPs) in Globular Clusters (GCs) by using diagrams that are based on colours and magnitudes that are mostly sensitive to stars with different metallicities or different abundances of helium, carbon, nitrogen and oxygen. In particular, the pseudo two-colour diagram called chromosome map (ChM), allowed the identification and the characterization of MPs in about 70 GCs by using appropriate filters of the Hubble Space Telescope (HST) that are sensitive to the stellar content of He, C, N, O and Fe. We use here high-precision HST photometry from F275W, F280N, F343N, F373N, and F814W images of Omega Centauri to investigate its MPs. We introduce a new ChM whose abscissa and ordinate are mostly sensitive to stellar populations with different magnesium and nitrogen, respectively, in monometallic GCs. This ChM is effective in disentangling the MPs based on their Mg chemical abundances, allowing us to explore, for the first time, possible relations between the production of these elemental species for large samples of stars in GCs. By comparing the colours of the distinct stellar populations with the colours obtained from appropriate synthetic spectra we provide 'photometric-like' estimates of the chemical composition of each population. Our results show that, in addition to first generation (1G) stars, the metal-poor population of Omega Centauri hosts four groups of second-generation stars with different [N/Fe], namely, 2GA--D. 2GA stars share nearly the same [Mg/Fe] as the 1G, whereas 2GB, 2GC and 2GD stars are Mg depleted by ~0.15, ~0.25 and ~0.45 dex, respectively. We provide evidence that the metal-intermediate populations host stars with depleted [Mg/Fe].


title: Injectivity and range description of first $(k+1)$ integral moment transforms over $m$-tensor fields in $\mathbb{R}^n$

abstract: In this work, we prove a new decomposition result for rank $m$ symmetric tensor fields which generalizes the well known solenoidal and potential decomposition of tensor fields. This decomposition is then used to describe the kernel and to prove an injectivity result for first $(k+1)$ integral moment transforms of symmetric $m$-tensor fields in $\mathbb{R}^n$. Additionally, we also present a range characterization for first $(k+1)$ integral moment transforms in terms of the John's equation.


title: Recurrence relations of Li coefficients

abstract: One of equivalents of the Riemann hypothesis is Li's criterion that all Li coefficients are positive. We study recurrence relations of Li coefficients in this note.


title: An improved sample size calculation method for score tests in generalized linear models

abstract: Self and Mauritsen (1988) developed a sample size determination procedure for score tests in generalized linear models under contiguous alternatives. Its performance may deteriorate when the effect size is large. We propose a modification of the Self-Mauritsen method by taking into account of the variance of the score statistic under both the null and alternative hypotheses, and extend the method to noninferiority trials. The modified approach is employed to calculate the sample size for the logistic regression and negative binomial regression in superiority and noninferiority trials. We further explain why the formulae recently derived by Zhu and Lakkis tend to underestimate the required sample size for the negative binomial regression. Numerical examples are used to demonstrate the accuracy of the proposed method.


title: Differentiable Segmentation of Sequences

abstract: Segmented models are widely used to describe non-stationary sequential data with discrete change points. Their estimation usually requires solving a mixed discrete-continuous optimization problem, where the segmentation is the discrete part and all other model parameters are continuous. A number of estimation algorithms have been developed that are highly specialized for their specific model assumptions. The dependence on non-standard algorithms makes it hard to integrate segmented models in state-of-the-art deep learning architectures that critically depend on gradient-based optimization techniques. In this work, we formulate a relaxed variant of segmented models that enables joint estimation of all model parameters, including the segmentation, with gradient descent. We build on recent advances in learning continuous warping functions and propose a novel family of warping functions based on the two-sided power (TSP) distribution. TSP-based warping functions are differentiable, have simple closed-form expressions, and can represent segmentation functions exactly. Our formulation includes the important class of segmented generalized linear models as a special case, which makes it highly versatile. We use our approach to model the spread of COVID-19 with Poisson regression, apply it on a change point detection task, and learn classification models with concept drift. The experiments show that our approach effectively learns all these tasks with standard algorithms for gradient descent.


title: Emergence of Spacetime: From Entanglement to Einstein

abstract: Here I develop the connection between thermodynamics, entanglement, and gravity. I begin by showing that the classical null energy condition (NEC) can arise as a consequence of the second law of thermodynamics applied to local holographic screens. This is accomplished by essentially reversing the steps of Hawking's area theorem, leading to the Ricci convergence condition as an input, from which an application of Einstein's equations yields the NEC -- even in the presence of 1-loop quantum corrections to the Bekenstein-Hawking entropy formula. Then, by attributing thermodynamics to the stretched horizon of future lightcones -- a timelike hypersurface generated by a collection of radially accelerating observers with constant and uniform proper acceleration -- I derive Einstein's equations from the Clausius relation $T\Delta S_{\text{rev}}=Q$, where $\Delta S_{\text{rev}}$ is the reversible entropy change. Based on this derivation I uncover a local first law of gravity, $\Delta E=T\Delta S-W$, connecting gravitational entropy $S$ to matter energy $E$ and work $W$. I then provide an entanglement interpretation of stretched lightcone thermodynamics by extending the entanglement equilibrium proposal. Using the $\text{AdS}_{3}/\text{CFT}_{2}$ correspondence, I then provide a microscopic explanation of the `thermodynamic volume' in extended black hole thermodynamics and reveal the super-entropicity of $\text{AdS}_{3}$ black holes is due to the gravitational entropy overcounting the number of available dual $\text{CFT}_{2}$ states. Finally, I conclude by providing a recent generalization of the extended first law of entanglement, and study its non-trivial 2+1- and 1+1-dimensional limits, including an extended first law for Jackiw-Teitelboim gravity. This thesis is self-contained and pedagogical by including useful background content relevant to emergent gravity.


title: Fast, Optimal, and Targeted Predictions using Parametrized Decision Analysis

abstract: Prediction is critical for decision-making under uncertainty and lends validity to statistical inference. With targeted prediction, the goal is to optimize predictions for specific decision tasks of interest, which we represent via functionals. Although classical decision analysis extracts predictions from a Bayesian model, these predictions are often difficult to interpret and slow to compute. Instead, we design a class of parametrized actions for Bayesian decision analysis that produce optimal, scalable, and simple targeted predictions. For a wide variety of action parametrizations and loss functions--including linear actions with sparsity constraints for targeted variable selection--we derive a convenient representation of the optimal targeted prediction that yields efficient and interpretable solutions. Customized out-of-sample predictive metrics are developed to evaluate and compare among targeted predictors. Through careful use of the posterior predictive distribution, we introduce a procedure that identifies a set of near-optimal, or acceptable targeted predictors, which provide unique insights into the features and level of complexity needed for accurate targeted prediction. Simulations demonstrate excellent prediction, estimation, and variable selection capabilities. Targeted predictions are constructed for physical activity data from the National Health and Nutrition Examination Survey (NHANES) to better predict and understand the characteristics of intraday physical activity.


title: Distilling Object Detectors with Task Adaptive Regularization

abstract: Current state-of-the-art object detectors are at the expense of high computational costs and are hard to deploy to low-end devices. Knowledge distillation, which aims at training a smaller student network by transferring knowledge from a larger teacher model, is one of the promising solutions for model miniaturization. In this paper, we investigate each module of a typical detector in depth, and propose a general distillation framework that adaptively transfers knowledge from teacher to student according to the task specific priors. The intuition is that simply distilling all information from teacher to student is not advisable, instead we should only borrow priors from the teacher model where the student cannot perform well. Towards this goal, we propose a region proposal sharing mechanism to interflow region responses between the teacher and student models. Based on this, we adaptively transfer knowledge at three levels, \emph{i.e.}, feature backbone, classification head, and bounding box regression head, according to which model performs more reasonably. Furthermore, considering that it would introduce optimization dilemma when minimizing distillation loss and detection loss simultaneously, we propose a distillation decay strategy to help improve model generalization via gradually reducing the distillation penalty. Experiments on widely used detection benchmarks demonstrate the effectiveness of our method. In particular, using Faster R-CNN with FPN as an instantiation, we achieve an accuracy of $39.0\%$ with Resnet-50 on COCO dataset, which surpasses the baseline $36.3\%$ by $2.7\%$ points, and even better than the teacher model with $38.5\%$ mAP.


title: An Agent-based Cloud Service Negotiation in Hybrid Cloud Computing

abstract: With the advent of evolution of cloud computing, large organizations have been scaling the on-premise IT infrastructure to the cloud. Although this being a popular practice, it lacks comprehensive efforts to study the aspects of automated negotiation of resources among cloud customers and providers. This paper proposes a full-fledged framework for the multi-party, multi-issue negotiation system for cloud resources. It introduces a robust cloud marketplace system to buy and sell cloud resources. The Belief-Desire-Intention (BDI) model-based cloud customer and provider agents concurrently negotiate on multiple issues, pursuing a hybrid tactic of time and resource-based dynamic deadline algorithms to generate offers and counter-offers. The cloud marketplace-based system is further augmented with the assignment of behavior norm score and reputation index to the agents to establish trust among them.


title: Optimized probes of $CP$-odd effects in the $t \bar{t} h$ process at hadron colliders

abstract: We use machine learning (ML) and non-ML techniques to study optimized $CP$-odd observables, directly and maximally sensitive to the $CP$-odd $i \tilde \kappa \bar t \gamma^5 t h$ interaction at the LHC and prospective future hadron colliders using the final state with a Higgs boson and a top quark pair, $pp\to t\bar t h$, followed by semileptonic $t$ decays. We perform phase-space optimization of manifestly $CP$-odd observables ($\boldsymbol \omega$), sensitive to the sign of $\tilde \kappa$, and constructed from experimentally accessible final state momenta. We identify a simple optimized linear combination $\boldsymbol \alpha\cdot \boldsymbol\omega$ that gives similar sensitivity as the studied fully fledged ML models. Using $\boldsymbol\alpha\cdot \boldsymbol\omega$ we project the expected sensitivities to $\tilde \kappa$ at HL-LHC, HE-LHC, and FCC-hh.


title: Order-chaos-order and invariant manifolds in the bounded planar Earth-Moon system

abstract: In this work, we investigate the Earth-Moon system, as modeled by the planar circular restricted three-body problem, and relate its dynamical properties to the underlying structure associated with specific invariant manifolds. We consider a range of Jacobi constant values for which the neck around the Lagrangian point $L_1$ is always open but the orbits are bounded due to Hill stability. First, we show that the system displays three different dynamical scenarios in the neighborhood of the Moon: two mixed ones, with regular and chaotic orbits, and an almost entirely chaotic one in between. We then analyze the transitions between these scenarios using the Monodromy matrix theory and determine that they are given by two specific types of bifurcations. After that, we illustrate how the phase space configurations, particularly the shapes of stability regions and stickiness, are intrinsically related to the hyperbolic invariant manifolds of the Lyapunov orbits around $L_1$ and also to the ones of some particular unstable periodic orbits. Lastly, we define transit time in a manner that is useful to depict dynamical trapping and show that the traced geometrical structures are also connected to the transport properties of the system.


title: Optimised allgatherv, reduce_scatter and allreduce communication in message-passing systems

abstract: Collective communications, namely the patterns allgatherv, reduce_scatter, and allreduce in message-passing systems are optimised based on measurements at the installation time of the library. The algorithms used are set up in an initialisation phase of the communication, similar to the method used in so-called persistent collective communication introduced in the literature. For allgatherv and reduce_scatter the existing algorithms, recursive multiply/divide and cyclic shift (Bruck's algorithm) are applied with a flexible number of communication ports per node. The algorithms for equal message sizes are used with non-equal message sizes together with a heuristic for rank reordering. The two communication patterns are applied in a plasma physics application that uses a specialised matrix-vector multiplication. For the allreduce pattern the cyclic shift algorithm is applied with a prefix operation. The data is gathered and scattered by the cores within the node and the communication algorithms are applied across the nodes. In general our routines outperform the non-persistent counterparts in established MPI libraries by up to one order of magnitude or show equal performance, with a few exceptions of number of nodes and message sizes.


title: Proton polarons in HxWO3 by synchrotron photoemission and DFT modelling

abstract: The measurement of hydrogen induced changes on the electronic structure of transition metal oxides by X-ray photoelectron spectroscopy is a challenging endeavor, since the origin of the photoelectron cannot be unambiguously assigned to hydrogen. The H-induced electronic structure changes in tungsten trioxide have been known for more than 100 years, but are still being controversially debated. The controversy stems from the difficulty in disentangling effects due to hydrogenation from the effects of oxygen deficiencies. Using a membrane approach to X-ray photoelectron spectroscopy, in combination with tuneable synchrotron radiation we measure simultaneously core levels and valence band up to a hydrogen pressure of 1000 mbar. Upon hydrogenation, the intensities of the W$^{5+}$ core level and a state close to the Fermi level increase following the pressure-composition isotherm curve of bulk H$_x$WO$_3$. Combining experimental data and density-functional theory the description of the hydrogen induced coloration by a proton polaron model is corroborated. Although hydrogen is the origin of the electronic structure changes near the Fermi edge, the valence band edge is now dominated by tungsten orbitals instead of oxygen as is the case for the pristine oxide having wider implication for its use as (photo-electrochemical) catalyst.


title: Fairness without Demographics through Adversarially Reweighted Learning

abstract: Much of the previous machine learning (ML) fairness literature assumes that protected features such as race and sex are present in the dataset, and relies upon them to mitigate fairness concerns. However, in practice factors like privacy and regulation often preclude the collection of protected features, or their use for training or inference, severely limiting the applicability of traditional fairness research. Therefore we ask: How can we train an ML model to improve fairness when we do not even know the protected group memberships? In this work we address this problem by proposing Adversarially Reweighted Learning (ARL). In particular, we hypothesize that non-protected features and task labels are valuable for identifying fairness issues, and can be used to co-train an adversarial reweighting approach for improving fairness. Our results show that {ARL} improves Rawlsian Max-Min fairness, with notable AUC improvements for worst-case protected groups in multiple datasets, outperforming state-of-the-art alternatives.


title: On four families of power series involving harmonic numbers and central binomial coefficients

abstract: We present several sequences involving harmonic numbers and the central binomial coefficients. The calculational technique is consists of a special summation method that allows, based on proper two-valued integer functions, to calculate different families of power series which involve odd harmonic numbers and central binomial coefficients. Furthermore it is shown that based on these series a new type of nonlinear Euler sums that involve odd harmonic numbers can be calculated in terms of zeta functions.


title: Isobar correlations bearing information on the properties of hot disassembling nuclear sources

abstract: Two-particle correlations based on the multiplicity of selected isobars are found to be sensitive to the parameterization of the fragments' binding energies and the breakup volume assumed in the model calculations. The properties of these correlations have been examined in the framework of the Statistical Multifragmentation Model as a function of the breakup temperature. The model calculations suggest that the maxima of these correlation functions occur at well separated temperatures as the breakup volumes used in the model vary from 3 to 6 times that at normal density. These volumes are within the range assumed in most statistical calculations and supported by experiments. Besides their position, the height and width of the maxima are also found to be sensitive to the parameterization of the fragments' binding energy. The magnitude of all these effects also depends on the isobars considered in the correlations. We found that, due to an interplay between the symmetry energy and the volume dependent terms of the Helmholtz free energy, in the case of nearly symmetric sources, correlations involving light mirror nuclei seem to enhance these effects. We suggest that the proposed correlation functions could be used to extract information on the fragments' energies and on the breakup volume of nuclear sources.


title: Applications of 10 MeV Photons and 100 MeV Protons in Radiotherapy

abstract: Geant4 software is used to study 10 MeV photon beams and 100 MeV proton beams. Conclusions are drawn as to the advantages and disadvantages of each therapy, with some additional commentary on the secondary effects of radiation damage in treatment.


title: Active model split hybrid RANS/LES

abstract: Reliably predictive simulation of complex flows requires a level of model sophistication and robustness exceeding the capabilities of current Reynolds-averaged Navier-Stokes (RANS) models. The necessary capability can often be provided by well-resolved large eddy simulation (LES), but, for many flows of interest, such simulations are too computationally intensive to be performed routinely. In principle, hybrid RANS/LES (HRL) models capable of transitioning through arbitrary levels of modeled and resolved turbulence would ameliorate both RANS deficiencies and LES expense. However, these HRL approaches have led to a host of unique complications, in addition to those already present in RANS and LES. This work proposes a modeling approach aimed at overcoming such challenges. The approach presented here relies on splitting the turbulence model into three distinct components: two responsible for the standard subgrid model roles of either providing the unresolved stress or dissipation and a third which reduces the model length scale by creating resolved turbulence. This formulation renders blending functions unnecessary in HRL. Further, the split-model approach both reduces the physics-approximation burden on simple eddy-viscosity-based models and provides convenient flexibility in model selection. In regions where the resolution is adequate to support additional turbulence, fluctuations are generated at the smallest locally resolved scales of motion. This active forcing drives the system towards a balance between RANS and grid-resolved LES for any combination of resolution and flow while the split-model formulation prevents local disruption to the total stress. The model is demonstrated on fully-developed, incompressible channel flow [1] and the periodic hill [2], in which it is shown to produce accurate results and avoid common HRL shortcomings, such as model stress depletion.


title: Oxide Fermi liquid universality revealed by electron spectroscopy

abstract: We present a combined soft x-ray and high-resolution vacuum-ultraviolet angle-resolved photoemission spectroscopy study of the electron-overdoped cuprate Pr$_{1.3-x}$La$_{0.7}$Ce$_{x}$CuO$_4$ (PLCCO). Demonstration of its highly two-dimensional band structure enabled precise determination of the in-plane self-energy dominated by electron-electron scattering. Through analysis of this self-energy and the Fermi-liquid cut-off energy scale, we find -- in contrast to hole-doped cuprates -- a momentum isotropic and comparatively weak electron correlation in PLCCO. Yet, the self-energies extracted from multiple oxide systems combine to demonstrate a logarithmic divergent relation between the quasiparticle scattering rate and mass. This constitutes a spectroscopic version of the Kadowaki-Woods relation with an important merit -- the demonstration of Fermi liquid quasiparticle lifetime and mass being set by a single energy scale.


title: Discrete Few-Shot Learning for Pan Privacy

abstract: In this paper we present the first baseline results for the task of few-shot learning of discrete embedding vectors for image recognition. Few-shot learning is a highly researched task, commonly leveraged by recognition systems that are resource constrained to train on a small number of images per class. Few-shot systems typically store a continuous embedding vector of each class, posing a risk to privacy where system breaches or insider threats are a concern. Using discrete embedding vectors, we devise a simple cryptographic protocol, which uses one-way hash functions in order to build recognition systems that do not store their users' embedding vectors directly, thus providing the guarantee of computational pan privacy in a practical and wide-spread setting.


title: An Improved Transmission Switching Algorithm for Managing Post-(N-1) Contingency Violations in Electricity Networks

abstract: This letter is a proof of concept for an improved transmission switching (TS) performance by moving the search space to load shed buses. Research from the past shows that changing transmission system topology changes the power flows and removes post contingency violations. Hence, TS can reduce the amount of load shed after an N-1 contingency. One of the major challenges is to find the best TS candidate in a suitable time. In this letter, the best TS candidate is determined by using a novel heuristic bi-level method based on linear sensitivity. The proposed bi-level method is easy to implement in the real world, guarantees removal of post contingency violations, and ranks the best TS candidates based on minimum load shedding possible. Moreover, the proposed method is computationally efficient since it does not involve mixed integer programming. The bi-level method is implemented by modifying the topology of transmission system after the N-1 contingency in the IEEE 39-bus test system and results show that TS with generation re-dispatch is the best solution for load shed recovery to prevent cascading failures. Moreover, the bi-level method performs even for the case when the existing methods in literature fail to completely remove post contingency violations.


title: Perelomov type coherent states of SO(D + 1) in all dimensional loop quantum gravity

abstract: A comprehensive study of the application of SO$(D+1)$ coherent states of Perelomov type to loop quantum gravity in general spacetime dimensions $D+1\geq 3$ is given in this paper. We focus on so-called simple representations of SO$(D+1)$ which solve the simplicity constraint and the associated homogeneous harmonic function spaces. With the harmonic function formulation, we study general properties of the coherent states such as the peakedness properties and the inner product. We also discuss the properties of geometric operators evaluated in the coherent states. In particular, we calculate the expectation value of the volume operator, and the results agree with the ones obtained from the classical label of the coherent states up to error terms which vanish in the limit of large representation labels $N$, i.e. the analogue of the large spin limit in standard $3+1$-dimensional loop quantum gravity.


title: Two superconducting components with different symmetries in Nd1-xSrxNiO2 films

abstract: The pairing mechanism in cuprates remains as one of the most challenging issues in the field of condensed matter physics. The unique 3d9 electron orbital of the Cu2+ ionic states in cuprates is supposed to be the major player for the occurrence of superconductivity. Recently, superconductivity at about 9-15 K was discovered in infinite layer thin films of nickelate Nd1-xSrxNiO2 (x=0.1-0.2) which is believed to have the similar 3d9 orbital electrons. The key issue concerned here is about the superconducting gap function. Here we report the first set data of single particle tunneling measurements on the superconducting nickelate thin films. We find predominantly two types of tunneling spectra, one shows a V-shape feature which can be fitted very well by a d-wave gap function with gap maximum of about 3.9 meV, another one exhibits a full gap of about 2.35 meV. Some spectra demonstrate mixed contributions of these two components. Our results suggest that the newly found Ni-based superconductors play as close analogs to cuprates, and thus demonstrate the commonality of unconventional superconductivity.


title: Analysis of Methanol Maser Flares in G107.298+5.63 andS255-NIRS3

abstract: A 3D maser model has been used to perform an inverse problem on the light curves from three high-amplitude maser flares, selected on the basis of contemporaneous infra-red observations. Plots derived from the model recover the size of the maser cloud, and two parameters linked to saturation, from three observational properties of the light curve. Recovered sizes are consistent with independent interferometric measurements. Maser objects transition between weak and moderate saturation during a flare.


title: DeepQTMT: A Deep Learning Approach for Fast QTMT-based CU Partition of Intra-mode VVC

abstract: Versatile Video Coding (VVC), as the latest standard, significantly improves the coding efficiency over its ancestor standard High Efficiency Video Coding (HEVC), but at the expense of sharply increased complexity. In VVC, the quad-tree plus multi-type tree (QTMT) structure of coding unit (CU) partition accounts for over 97% of the encoding time, due to the brute-force search for recursive rate-distortion (RD) optimization. Instead of the brute-force QTMT search, this paper proposes a deep learning approach to predict the QTMT-based CU partition, for drastically accelerating the encoding process of intra-mode VVC. First, we establish a large-scale database containing sufficient CU partition patterns with diverse video content, which can facilitate the data-driven VVC complexity reduction. Next, we propose a multi-stage exit CNN (MSE-CNN) model with an early-exit mechanism to determine the CU partition, in accord with the flexible QTMT structure at multiple stages. Then, we design an adaptive loss function for training the MSE-CNN model, synthesizing both the uncertain number of split modes and the target on minimized RD cost. Finally, a multi-threshold decision scheme is developed, achieving desirable trade-off between complexity and RD performance. Experimental results demonstrate that our approach can reduce the encoding time of VVC by 44.65%-66.88% with the negligible Bj{\o}ntegaard delta bit-rate (BD-BR) of 1.322%-3.188%, which significantly outperforms other state-of-the-art approaches.


title: Fixing Inventory Inaccuracies At Scale

abstract: Inaccurate records of inventory occur frequently, and by some measures cost retailers approximately 4% in annual sales. Detecting inventory inaccuracies manually is cost-prohibitive, and existing algorithmic solutions rely almost exclusively on learning from longitudinal data, which is insufficient in the dynamic environment induced by modern retail operations. Instead, we propose a solution based on cross-sectional data over stores and SKUs, observing that detecting inventory inaccuracies can be viewed as a problem of identifying anomalies in a (low-rank) Poisson matrix. State-of-the-art approaches to anomaly detection in low-rank matrices apparently fall short. Specifically, from a theoretical perspective, recovery guarantees for these approaches require that non-anomalous entries be observed with vanishingly small noise (which is not the case in our problem, and indeed in many applications). So motivated, we propose a conceptually simple entry-wise approach to anomaly detection in low-rank Poisson matrices. Our approach accommodates a general class of probabilistic anomaly models. We show that the cost incurred by our algorithm approaches that of an optimal algorithm at a min-max optimal rate. Using synthetic data and real data from a consumer goods retailer, we show that our approach provides up to a 10x cost reduction over incumbent approaches to anomaly detection. Along the way, we build on recent work that seeks entry-wise error guarantees for matrix completion, establishing such guarantees for sub-exponential matrices, a result of independent interest.


title: Rigorous computer-assisted bounds on the period doubling renormalisation fixed point and eigenfunctions in maps with critical point of degree 4

abstract: We gain tight rigorous bounds on the renormalisation fixed point for period doubling in families of unimodal maps with degree $4$ critical point. We use a contraction mapping argument to bound essential eigenfunctions and eigenvalues for the linearisation of the operator and for the operator controlling the scaling of added noise. Multi-precision arithmetic with rigorous directed rounding is used to bound operations in a space of analytic functions yielding tight bounds on power series coefficients and universal constants to over $320$ significant figures.


title: Three-dimensional Maxwellian Extended Newtonian gravity and flat limit

abstract: In the present work we find novel Newtonian gravity models in three spacetime dimensions. We first present a Maxwellian version of the extended Newtonian gravity, which is obtained as the non-relativistic limit of a particular $U(1)$-enlargement of an enhanced Maxwell Chern-Simons gravity. We show that the extended Newtonian gravity appears as a particular sub-case. Then, the introduction of a cosmological constant to the Maxwellian extended Newtonian theory is also explored. To this purpose, we consider the non-relativistic limit of an enlarged symmetry. An alternative method to obtain our results is presented by applying the semigroup expansion method to the enhanced Nappi-Witten algebra. The advantages of considering the Lie algebra expansion procedure is also discussed.


title: A Catalogue of Game-Specific Software Nuggets

abstract: With the ever-increasing use of games, game developers are expected to write efficient code supporting several qualities such as security, maintainability, and performance. However, the continuous need to update the features of games in less duration might compel the developers to use anti-patterns, code smells and quick-fix solutions that may affect the functional and non-functional requirements of the game. These bad practices may lead to technical debt, poor program comprehension, and can cause several issues during software maintenance. Here, in this paper, we introduce "Software Nuggets" as a concept that affects software quality in a negative way and as a superset of anti-patterns, code smells, bugs, software bad practices. We call these Software Nuggets as "G-Nuggets" in the context of games. While there exists empirical research on games, we are not aware of any work on understanding and cataloguing these G-Nuggets. Thus, we propose a catalogue of G-Nuggets by mining and analyzing 892 commits, 189 issues, and 104 pull requests from 100 open-source GitHub game repositories. We use regular expressions and thematic analysis on this dataset for cataloguing game-specific Software Nuggets. We present a catalogue of ten G-Nuggets and provide examples for them present online at: https://phoebs88.github.io/A-Catalogue-of-Game-Specific-Software-Nuggets. We believe this catalogue might be helpful for researchers for further empirical research in the domain of games as well as for game developers to improve quality of games.


title: Imaging non-collinear antiferromagnetic textures via single spin relaxometry

abstract: Antiferromagnetic materials are promising platforms for next-generation spintronics owing to their fast dynamics and high robustness against parasitic magnetic fields. However, nanoscale imaging of the magnetic order in such materials with zero net magnetization remains a major experimental challenge. Here we show that non-collinear antiferromagnetic spin textures can be imaged by probing the magnetic noise they locally produce via thermal populations of magnons. To this end, we perform nanoscale, all-optical relaxometry with a scanning quantum sensor based on a single nitrogen-vacancy (NV) defect in diamond. Magnetic noise is detected through an increase of the spin relaxation rate of the NV defect, which results in an overall reduction of its photoluminescence signal under continuous laser illumination. As a proof-of-concept, the efficiency of the method is demonstrated by imaging various spin textures in synthetic antiferromagnets, including domain walls, spin spirals and antiferromagnetic skyrmions. This imaging procedure could be extended to a large class of intrinsic antiferromagnets and opens up new opportunities for studying the physics of localized spin wave modes for magnonics.


title: Computational modelling of cancer evolution by multi-type branching processes

abstract: Metastasis, the spread of cancer cells from a primary tumor to secondary location(s) in the human organism, is the ultimate cause of death for the majority of cancer patients. That is why, it is crucial to understand metastases evolution in order to successfully combat the disease. We consider a metastasized cancer cell population after medical treatment (e.g. chemotherapy). Arriving in a different environment the cancer cells may change their lifespan and reproduction, thus they may proliferate into different types. If the treatment is effective, in the context of branching processes it means, the reproduction of cancer cells is such that the mean offspring of each cell is less than one. However, it is possible mutations to occur during cell division cycle. These mutations can produce a new cancer cell type, which is resistant to the treatment. Cancer cells from this new type may lead to the rise of a non-extinction branching process. The above scenario leads us to the choice of a reducible multi-type age-dependent branching process as a relevant framework for studying the asymptotic behavior of such complex structures. Our previous theoretical results are related to the asymptotic behavior of the waiting time until the first occurrence of a mutant starting a non-extinction process and the modified hazard function as a measure of immediate recurrence of cancer disease. In the present paper these asymptotic results are used for developing numerical schemes and algorithms implemented in Python via the NumPy package for approximate calculation of the corresponding quantities. In conclusion, our conjecture is that this methodology can be advantageous in revealing the role of the lifespan distribution of the cancer cells in the context of cancer disease evolution and other complex cell population systems, in general.


title: On Counterfactual Explanations under Predictive Multiplicity

abstract: Counterfactual explanations are usually obtained by identifying the smallest change made to an input to change a prediction made by a fixed model (hereafter called sparse methods). Recent work, however, has revitalized an old insight: there often does not exist one superior solution to a prediction problem with respect to commonly used measures of interest (e.g. error rate). In fact, often multiple different classifiers give almost equal solutions. This phenomenon is known as predictive multiplicity (Breiman, 2001; Marx et al., 2019). In this work, we derive a general upper bound for the costs of counterfactual explanations under predictive multiplicity. Most notably, it depends on a discrepancy notion between two classifiers, which describes how differently they treat negatively predicted individuals. We then compare sparse and data support approaches empirically on real-world data. The results show that data support methods are more robust to multiplicity of different models. At the same time, we show that those methods have provably higher cost of generating counterfactual explanations under one fixed model. In summary, our theoretical and empiricaln results challenge the commonly held view that counterfactual recommendations should be sparse in general.


title: Autonomous and Semi-Autonomous Intersection Management: A Survey

abstract: Intersection is a major source of traffic delays and accidents within modern transportation systems. Compared to signalized intersection management, autonomous intersection management (AIM) coordinates the intersection crossing at an individual vehicle level with additional flexibility. AIM can potentially eliminate stopping in intersection crossing due to traffic lights while maintaining a safe separation among conflicting movements. In this paper, the state-of-the-art AIM research among various disciplines (e.g., traffic engineering, control engineering) is surveyed from the perspective of three hierarchical layers: corridor coordination layer, intersection management layer, and vehicle control layer. The key aspects of AIM designs are discussed in details, including conflict detection schemes, priority rules, control centralization, computation complexity, etc. The potential improvements for AIM evaluation with the emphasis of realistic scenarios are provided. This survey serves as a comprehensive review of AIM design and provides promising directions for future research.


title: Interpreting Holographic Molecular Binding Assays with Effective Medium Theory

abstract: Holographic molecular binding assays use holographic video microscopy to directly detect molecules binding to the surfaces of micrometer-scale colloidal beads by monitoring associated changes in the beads' light-scattering properties. Holograms of individual spheres are analyzed by fitting to a generative model based on the Lorenz-Mie theory of light scattering. Each fit yields an estimate of a probe bead's diameter and refractive index with sufficient precision to watch the beads grow as molecules bind. Rather than modeling the molecular-scale coating, however, these fits use effective medium theory, treating the coated sphere as if it were homogeneous. This effective-sphere analysis is rapid and numerically robust and so is useful for practical implementations of label-free immunoassays. Here, we assess how effective-sphere properties reflect the properties of molecular-scale coatings by modeling coated spheres with the discrete-dipole approximation and analyzing their holograms with the effective-sphere model.


title: Estimation of Causal Effects in the Presence of Unobserved Confounding in the Alzheimer's Continuum

abstract: Studying the relationship between neuroanatomy and cognitive decline due to Alzheimer's has been a major research focus in the last decade. However, to infer cause-effect relationships rather than simple associations from observational data, we need to (i) express the causal relationships leading to cognitive decline in a graphical model, and (ii) ensure the causal effect of interest is identifiable from the collected data. We derive a causal graph from the current clinical knowledge on cause and effect in the Alzheimer's disease continuum, and show that identifiability of the causal effect requires all confounders to be known and measured. However, in complex neuroimaging studies, we neither know all potential confounders nor do we have data on them. To alleviate this requirement, we leverage the dependencies among multiple causes by deriving a substitute confounder via a probabilistic latent factor model. In our theoretical analysis, we prove that using the substitute confounder enables identifiability of the causal effect of neuroanatomy on cognition. We quantitatively evaluate the effectiveness of our approach on semi-synthetic data, where we know the true causal effects, and illustrate its use on real data on the Alzheimer's disease continuum, where it reveals important causes that otherwise would have been missed.


title: Atomistic Mechanism of Phase Transition in Shock Compressed Gold Revealed by Deep Potential

abstract: A detailed understanding of the material response to rapid compression is challenging and demanding. For instance, the element gold under dynamic compression exhibits complex phase transformations where there exist some large discrepancies between experimental and theoretical studies. Here, we combined large-scale molecular dynamics simulations with a deep potential to elucidate the dynamic compression processes of gold from an atomic level. The potential is constructed by accurately reproducing the free energy surfaces of density-functional-theory calculations for gold, from ambient conditions to 15 500 K and 500 GPa. Within this framework, we extend the simulations up to 200 000 atoms size, and found a much lower pressure threshold for phase transitioning from face-centered cubic (FCC) to body-centered (BCC), as compared to previous calculations. Furthermore, the transition pressure is strongly dependent on the shock direction, namely 159 GPa for (100) orientation and 219 GPa for (110) orientation, respectively. Most importantly, the accurate atomistic perspective presents that the shocked BCC structure contains unique features of medium-range and short-range orders, which is named disorders here. We propose a model and demonstrate that the existence of disorders significantly reduces the Gibbs free energies of shocked structures, therefore leading to the lowering of the phase transition pressure. The present study provides a new path to understand the structure dynamics under extreme conditions.


title: Primordial Black Holes from QCD Axion Bubbles

abstract: We propose a scenario in which a strong Peccei-Quinn (PQ) symmetry breaking in the early universe results in large inhomogeneities of the initial QCD axion field value, leading to the formation of very dense axion bubbles. Some of the axion bubbles subsequently collapse into primordial black holes (PBHs). The spatially homogeneous part of the QCD axion explains dark matter of the universe, while the PBHs arising from the axion bubbles can explain the LIGO events or the seed of supermassive black holes. Interestingly, the mass of PBH is determined by the axion decay constant; for $f_a = 10^{17} (10^{16})$ GeV, the PBH mass is heavier than about $10 (10^4) M_\odot$. In addition, axion miniclusters are also formed from the axion bubbles more abundantly than PBHs, and their masses are expected to be heavier than in the usual scenario based on the spontaneous breaking of the PQ symmetry after inflation.


title: hxtorch: PyTorch for BrainScaleS-2 -- Perceptrons on Analog Neuromorphic Hardware

abstract: We present software facilitating the usage of the BrainScaleS-2 analog neuromorphic hardware system as an inference accelerator for artificial neural networks. The accelerator hardware is transparently integrated into the PyTorch machine learning framework using its extension interface. In particular, we provide accelerator support for vector-matrix multiplications and convolutions; corresponding software-based autograd functionality is provided for hardware-in-the-loop training. Automatic partitioning of neural networks onto one or multiple accelerator chips is supported. We analyze implementation runtime overhead during training as well as inference, provide measurements for existing setups and evaluate the results in terms of the accelerator hardware design limitations. As an application of the introduced framework, we present a model that classifies activities of daily living with smartphone sensor data.


title: Phonon-mediated High-voltage Detector with Background Rejection for Low-mass Dark Matter and Reactor Coherent Neutrino Scattering Experiments

abstract: We report the first demonstration of a phonon-mediated silicon detector technology that provides a primary phonon measurement in a low-voltage region, and a simultaneous indirect measurement of the ionization signal through Neganov-Trofimov-Luke amplification in a high voltage region, both in a monolithic crystal. We present characterization of charge and phonon transport between the two stages of the detector and the resulting background discrimination capability at low energies. This new detector technology has the potential to significantly enhance the sensitivity of dark matter and coherent neutrino scattering experiments beyond the capabilities of current technologies that have limited discrimination at low energies.


title: A Mathematical Negotiation Mechanism for Distributed Procurement Problems and a Hybrid Algorithm for its Solution

abstract: In this paper, a mathematical negotiation mechanism is designed to minimize the negotiators' costs in a distributed procurement problem at two echelons of an automotive supply chain. The buyer's costs are procurement cost and shortage penalty in a one-period contract. On the other hand, the suppliers intend to solve a multi-period, multi-product production planning to minimize their costs. Such a mechanism provides an alignment among suppliers' production planning and order allocation, also supports the partnership with the valued suppliers by taking suppliers' capacities into account. Such a circumstance has been modeled via bi-level programming, in which the buyer acts as a leader, and the suppliers individually appear as followers in the lower level. To solve this nonlinear bi-level programming model, a hybrid algorithm by combining the particle swarm optimization (PSO) algorithm with a heuristic algorithm based on A search is proposed. The heuristic A algorithm is embedded to solve the mixed-integer nonlinear programming (MINLP) sub-problems for each supplier according to the received variable values determined by PSO system particles (buyer's request for quotations (RFQs)). The computational analyses have shown that the proposed hybrid algorithm called PSO-A outperforms PSO-SA and PSO-Greedy algorithms.


title: Alveolar mimics with periodic strain and its effect on the cell layer formation

abstract: We report on the development of a new model of alveolar air-tissue interface on a chip. The model consists of an array of suspended hexagonal monolayers of gelatin nanofibers supported by microframes and a microfluidic device for the patch integration. The suspended monolayers are deformed to a central displacement of 40-80 um at the air-liquid interface by application of air pressure in the range of 200-1000 Pa. With respect to the diameter of the monolayers that is 500 um, this displacement corresponds to a linear strain of 2-10% in agreement with the physiological strain range in the lung alveoli. The culture of A549 cells on the monolayers for an incubation time 1-3 days showed viability in the model. We exerted a periodic strain of 5% at a frequency of 0.2 Hz during 1 hour to the cells. We found that the cells were strongly coupled to the nanofibers, but the strain reduced the coupling and induced remodeling of the actin cytoskeleton, which led to a better tissue formation. Our model can serve as a versatile tool in lung investigations such as in inhalation toxicology and therapy.


title: A Deep Neural Network Model of Particle Thermal Radiation in Packed Bed

abstract: Prediction of particle radiative heat transfer flux is an important task in the large discrete granular systems, such as pebble bed in power plants and industrial fluidized beds. For particle motion and packing, discrete element method (DEM) now is widely accepted as the excellent Lagrangian approach. For thermal radiation, traditional methods focus on calculating the obstructed view factor directly by numerical algorithms. The major challenge for the simulation is that the method is proven to be time-consuming and not feasible to be applied in the practical cases. In this work, we propose an analytical model to calculate macroscopic effective conductivity from particle packing structures Then, we develop a deep neural network (DNN) model used as a predictor of the complex view factor function. The DNN model is trained by a large dataset and the computational speed is greatly improved with good accuracy. It is feasible to perform real-time simulation with DNN model for radiative heat transfer in large pebble bed. The trained model also can be coupled with DEM and used to analyze efficiently the directional radiative conductivity, anisotropic factor and wall effect of the particle thermal radiation.


title: Gaussian Continuous Tensor Network States for Simple Bosonic Field Theories

abstract: Tensor networks states allow to find the low energy states of local lattice Hamiltonians through variational optimization. Recently, a construction of such states in the continuum was put forward, providing a first step towards the goal of solving quantum field theories (QFTs) variationally. However, the proposed manifold of continuous tensor network states (CTNSs) is difficult to study in full generality, because the expectation values of local observables cannot be computed analytically. In this paper, we study a tractable subclass of CTNSs, the Gaussian CTNSs (GCTNSs), and benchmark them on simple quadratic and quartic bosonic QFT Hamiltonians. We show that GCTNSs provide arbitrarily accurate approximations to the ground states of quadratic Hamiltonians, and decent estimates for quartic ones at weak coupling. Since they capture the short distance behavior of the theories we consider exactly, GCTNSs even allow to renormalize away simple divergences variationally. In the end, our study makes it plausible that CTNSs are indeed a good manifold to approximate the low energy states of QFTs.