title: Parity alternating permutations starting with an odd integer

abstract: A Parity Alternating Permutation of the set $[n] = \{1, 2,\ldots, n\}$ is a permutation with even and odd entries alternatively. We deal with parity alternating permutations having an odd entry in the first position, PAPs. We study the numbers that count the PAPs with even as well as odd parity. We also study a subclass of PAPs being derangements as well, Parity Alternating Derangements (PADs). Moreover, by considering the parity of these PADs we look into their statistical property of excedance.


title: Will Artificial Intelligence supersede Earth System and Climate Models?

abstract: We outline a perspective of an entirely new research branch in Earth and climate sciences, where deep neural networks and Earth system models are dismantled as individual methodological approaches and reassembled as learning, self-validating, and interpretable Earth system model-network hybrids. Following this path, we coin the term "Neural Earth System Modelling" (NESYM) and highlight the necessity of a transdisciplinary discussion platform, bringing together Earth and climate scientists, big data analysts, and AI experts. We examine the concurrent potential and pitfalls of Neural Earth System Modelling and discuss the open question whether artificial intelligence will not only infuse Earth system modelling, but ultimately render them obsolete.


title: A new determination of the primordial helium abundance using the analyses of HII region spectra from SDSS

abstract: The precision measurement of the primordial helium abundance $Y_p$ is a powerful probe of the early Universe. The most common way to determine $Y_p$ is analyses of observations of metal-poor \HII regions found in blue compact dwarf galaxies. We present the spectroscopic sample of 100 \HII regions collected from the Sloan Digital Sky Survey. The final analysed sample consists of our sample and HeBCD database from Izotov et al. 2007. We use a self-consistent procedure to determine physical conditions, current helium abundances, and metallicities of the \HII regions. From a regression to zero metallicity, we have obtained $Y_p = 0.2462 \pm 0.0022$ which is one of the most stringent constraints obtained with such methods up to date and is in a good agreement with the Planck result $Y_{\rm p}^{\it {Planck}} = 0.2471 \pm 0.0003$. Using the determined value of $Y_p$ and the primordial deuterium abundance taken from Particle Data Group (Zyla et al. 2020) we put a constraint on the effective number of neutrino species $N_{\rm eff} = 2.95 \pm 0.16$ which is consistent with the Planck one $N_{\rm eff} = 2.99 \pm 0.17$. Further increase of statistics potentially allows us to achieve Planck accuracy, which in turn will become a powerful tool for studying the self-consistency of the Standard Cosmological Model and/or physics beyond.


title: An Efficient Model for Scaffold-Mediated Bone Regeneration

abstract: We present a three dimensional, time dependent model for bone regeneration in the presence of porous scaffolds to bridge critical size bone defects. Our approach uses homogenized quantities, thus drastically reducing computational cost compared to models resolving the microstructural scale of the scaffold. Using abstract functional relationships instead of concrete effective material properties, our model can incorporate the homogenized material tensors for a large class of scaffold microstructure designs. We prove an existence and uniqueness theorem for solutions based on a fixed point argument. We include the cases of mixed boundary conditions and multiple, interacting signalling molecules, both being important for application. Furthermore we present numerical simulations showing good agreement with experimental findings.


title: Solving the Same-Different Task with Convolutional Neural Networks

abstract: Deep learning demonstrated major abilities in solving many kinds of different real-world problems in computer vision literature. However, they are still strained by simple reasoning tasks that humans consider easy to solve. In this work, we probe current state-of-the-art convolutional neural networks on a difficult set of tasks known as the same-different problems. All the problems require the same prerequisite to be solved correctly: understanding if two random shapes inside the same image are the same or not. With the experiments carried out in this work, we demonstrate that residual connections, and more generally the skip connections, seem to have only a marginal impact on the learning of the proposed problems. In particular, we experiment with DenseNets, and we examine the contribution of residual and recurrent connections in already tested architectures, ResNet-18, and CorNet-S respectively. Our experiments show that older feed-forward networks, AlexNet and VGG, are almost unable to learn the proposed problems, except in some specific scenarios. We show that recently introduced architectures can converge even in the cases where the important parts of their architecture are removed. We finally carry out some zero-shot generalization tests, and we discover that in these scenarios residual and recurrent connections can have a stronger impact on the overall test accuracy. On four difficult problems from the SVRT dataset, we can reach state-of-the-art results with respect to the previous approaches, obtaining super-human performances on three of the four problems.


title: Bridging the gap between collisional and collisionless shock waves

abstract: While the front of a fluid shock is a few mean-free-paths thick, the front of a collisionless shock can be orders of magnitude thinner. By bridging between a collisional and a collisionless formalism, we assess the transition between these two regimes. We consider non-relativistic, un-magnetized, planar shocks in electron/ion plasmas. In addition, our treatment of the collisionless regime is restricted to high Mach number electrostatic shocks. We find that the transition can be parameterized by the upstream plasma parameter $\Lambda$ which measures the coupling of the upstream medium. For $\Lambda \lesssim 1.12$, the upstream is collisional, i.e. strongly coupled, and the strong shock front is about $\mathcal{M}_1 \lambda_{\mathrm{mfp},1}$ thick, where $\lambda_{\mathrm{mfp},1}$ and $\mathcal{M}_1 $ are the upstream mean-free-path and Mach number respectively. A transition occurs for $\Lambda \sim 1.12$ beyond which the front is $\sim \mathcal{M}_1\lambda_{\mathrm{mfp},1}\ln \Lambda/\Lambda$ thick for $\Lambda\gtrsim 1.12$. Considering $\Lambda$ can reach billions in astrophysical settings, this allows to understand how the front of a collisionless shock can be orders of magnitude smaller than the mean-free-path, and how physics transitions continuously between these 2 extremes.

categories: physics.plasm-ph astro-ph.HE astro-ph.SR physics.space-ph


title: Deep Inverse Design of Reconfigurable Metasurfaces for Future Communications

abstract: Reconfigurable intelligent surfaces (RIS) have recently received significant attention as building blocks for smart radio environments and adaptable wireless channels. By altering the space- and time-varying electromagnetic (EM) properties, the RIS transforms the inherently stochastic nature of the wireless environment into a programmable propagation channel. Conventionally, designing RIS to yield the desired EM response requires trial-and-error by iteratively investigating a large possibility of various geometries and materials through thousands of full-wave EM simulations. In this context, deep learning (DL) techniques are proving critical in reducing the computational cost and time of RIS inverse design. Instead of explicitly solving Maxwell's equations, DL models learn physics-based relationships through supervised training data. Further, generative adversarial networks are shown to synthesize novel RIS designs not previously seen in the literature. This article provides a synopsis of DL techniques for inverse RIS design and optimization to yield targeted EM response necessary for future wireless networks.


title: Generalized Newton-Leibniz Formula and the Embedding of the Sobolev Functions with Dominating Mixed Smoothness into H\"{o}lder Spaces

abstract: It is well-known that the embedding of the Sobolev space of weakly differentiable functions into H\"{o}lder spaces holds if the integrability exponent is higher than the space dimension. In this paper, the embedding of the Sobolev functions into the H\"{o}lder spaces is expressed in terms of the minimal weak differentiability requirement independent of the integrability exponent. The proof is based on the generalization of the Newton-Leibniz formula to the $n$-dimensional rectangle and inductive application of the Sobolev trace embedding results. The method is applied to prove the embedding of the Sobolev spaces with dominating mixed smoothness into H\"{o}lder spaces.


title: Isotopic Separation of Helium through Nanoporous Graphene Membranes: A Ring Polymer Molecular Dynamics Study

abstract: Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of $^4$He and $^3$He. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.


title: Boundary conditions for the N\'{e}el order parameter in a chiral antiferromagnetic slab

abstract: Understanding of the interaction of antiferromagnetic solitons including domain walls and skyrmions with boundaries of chiral antiferromagnetic slabs is important for the design of prospective antiferromagnetic spintronic devices. Here, we derive the transition from spin lattice to micromagnetic nonlinear $\sigma$-model with the corresponding boundary conditions for a chiral cubic G-type antiferromagnet and analyze the impact of the slab boundaries and antisymmetric exchange (Dzyaloshinskii--Moriya interaction) on the vector order parameter. We apply this model to evaluate modifications of antiferromagnetic domain walls and skyrmions upon interaction with boundaries for different strengths of the antisymmetric exchange. Due to the presence of the antisymmetric exchange, both types of antiferromagnetic solitons become broader when approaching the boundary and transform to a mixed Bloch--N\'{e}el structure. Both textures feel the boundary at the distance of about 5 magnetic lengths. In this respect, our model provides design rules for antiferromagnetic racetracks, which can support bulk-like properties of solitons.


title: A Survey on the Ricci flow on Singular Spaces

abstract: In this survey we provide an overview of our recent results concerning Ricci de Turck flow on spaces with isolated conical singularities. The crucial characteristic of the flow is that it preserves the conical singularity. Under certain conditions, Ricci flat metrics with isolated conical singularities are stable and positive scalar curvature is preserved under the flow. We also discuss the relation to Perelman's entropies in the singular setting, and outline open questions and future reseach directions.


title: Reconstruction of nanoscale particles from single-shot wide-angle FEL diffractions patterns with physics-informed neural networks

abstract: Single-shot wide-angle diffraction imaging is a widely used method to investigate the structure of non-crystallizing objects such as nanoclusters, large proteins or even viruses. Its main advantage is that information about the three-dimensional structure of the object is already contained in a single image. This makes it useful for the reconstruction of fragile and non-reproducible particles without the need for tomographic measurements. However, currently there is no efficient numerical inversion algorithm available that is capable of determining the object's structure in real time. Neural networks, on the other hand, excel in image processing tasks suited for such purpose. Here we show how a physics-informed deep neural network can be used to reconstruct complete three-dimensional object models of uniform, convex particles on a voxel grid from single two-dimensional wide-angle scattering patterns. We demonstrate its universal reconstruction capabilities for silver nanoclusters, where the network uncovers novel geometric structures that reproduce the experimental scattering data with very high precision.


title: Multi-hop RIS-Empowered Terahertz Communications: A DRL-based Hybrid Beamforming Design

abstract: Wireless communication in the TeraHertz band (0.1--10 THz) is envisioned as one of the key enabling technologies for the future sixth generation (6G) wireless communication systems scaled up beyond massive multiple input multiple output (Massive-MIMO) technology. However, very high propagation attenuations and molecular absorptions of THz frequencies often limit the signal transmission distance and coverage range. Benefited from the recent breakthrough on the reconfigurable intelligent surfaces (RIS) for realizing smart radio propagation environment, we propose a novel hybrid beamforming scheme for the multi-hop RIS-assisted communication networks to improve the coverage range at THz-band frequencies. Particularly, multiple passive and controllable RISs are deployed to assist the transmissions between the base station (BS) and multiple single-antenna users. We investigate the joint design of digital beamforming matrix at the BS and analog beamforming matrices at the RISs, by leveraging the recent advances in deep reinforcement learning (DRL) to combat the propagation loss. To improve the convergence of the proposed DRL-based algorithm, two algorithms are then designed to initialize the digital beamforming and the analog beamforming matrices utilizing the alternating optimization technique. Simulation results show that our proposed scheme is able to improve 50\% more coverage range of THz communications compared with the benchmarks. Furthermore, it is also shown that our proposed DRL-based method is a state-of-the-art method to solve the NP-hard beamforming problem, especially when the signals at RIS-assisted THz communication networks experience multiple hops.


title: LonelyText: A Short Messaging Based Classification of Loneliness

abstract: Loneliness does not only have emotional implications on a person but also on his/her well-being. The study of loneliness has been challenging and largely inconclusive in findings because of the several factors that might correlate to the phenomenon. We present one approach to predicting this event by discovering patterns of language associated with loneliness. Our results show insights and promising directions for mining text from instant messaging to predict loneliness.


title: The Privacy-Utility Tradeoff of Robust Local Differential Privacy

abstract: We consider data release protocols for data $X=(S,U)$, where $S$ is sensitive; the released data $Y$ contains as much information about $X$ as possible, measured as $\operatorname{I}(X;Y)$, without leaking too much about $S$. We introduce the Robust Local Differential Privacy (RLDP) framework to measure privacy. This framework relies on the underlying distribution of the data, which needs to be estimated from available data. Robust privacy guarantees are ensuring privacy for all distributions in a given set $\mathcal{F}$, for which we study two cases: when $\mathcal{F}$ is the set of all distributions, and when $\mathcal{F}$ is a confidence set arising from a $\chi^2$ test on a publicly available dataset. In the former case we introduce a new release protocol which we prove to be optimal in the low privacy regime. In the latter case we present four algorithms that construct RLDP protocols from a given dataset. One of these approximates $\mathcal{F}$ by a polytope and uses results from robust optimisation to yield high utility release protocols. However, this algorithm relies on vertex enumeration and becomes computationally inaccessible for large input spaces. The other three algorithms are low-complexity and build on randomised response. Experiments verify that all four algorithms offer significantly improved utility over regular LDP.


title: Hydrodynamic self-similar cosmological models

abstract: Hydrodynamic self-similar solutions, as obtained by Chi [J. Math. Phys. 24, 2532 (1983)] have been generalized by introducing new variables in place of the old space and time variables. A systematic procedure of obtaining a complete set of solutions has been suggested. The Newtonian analogs of all homogeneous isotropic Friedmann dust universes with spatial curvature $k = 0, \pm 1$ have been given.


title: A Computational Status Update for Exact Rational Mixed Integer Programming

abstract: The last milestone achievement for the roundoff-error-free solution of general mixed integer programs over the rational numbers was a hybrid-precision branch-and-bound algorithm published by Cook, Koch, Steffy, and Wolter in 2013. We describe a substantial revision and extension of this framework that integrates symbolic presolving, features an exact repair step for solutions from primal heuristics, employs a faster rational LP solver based on LP iterative refinement, and is able to produce independently verifiable certificates of optimality. We study the significantly improved performance and give insights into the computational behavior of the new algorithmic components. On the MIPLIB 2017 benchmark set, we observe an average speedup of 6.6x over the original framework and 2.8 times as many instances solved within a time limit of two hours.


title: A Study of Continuous Vector Representationsfor Theorem Proving

abstract: Applying machine learning to mathematical terms and formulas requires a suitable representation of formulas that is adequate for AI methods. In this paper, we develop an encoding that allows for logical properties to be preserved and is additionally reversible. This means that the tree shape of a formula including all symbols can be reconstructed from the dense vector representation. We do that by training two decoders: one that extracts the top symbol of the tree and one that extracts embedding vectors of subtrees. The syntactic and semantic logical properties that we aim to reserve include both structural formula properties, applicability of natural deduction steps, and even more complex operations like unifiability. We propose datasets that can be used to train these syntactic and semantic properties. We evaluate the viability of the developed encoding across the proposed datasets as well as for the practical theorem proving problem of premise selection in the Mizar corpus.


title: Traffic Flow Estimation using LTE Radio Frequency Counters and Machine Learning

abstract: As the demand for vehicles continues to outpace construction of new roads, it becomes imperative we implement strategies that improve utilization of existing transport infrastructure. Traffic sensors form a crucial part of many such strategies, giving us valuable insights into road utilization. However, due to cost and lead time associated with installation and maintenance of traffic sensors, municipalities and traffic authorities look toward cheaper and more scalable alternatives. Due to their ubiquitous nature and wide global deployment, cellular networks offer one such alternative. In this paper we present a novel method for traffic flow estimation using standardized LTE/4G radio frequency performance measurement counters. The problem is cast as a supervised regression task using both classical and deep learning methods. We further apply transfer learning to compensate that many locations lack traffic sensor data that could be used for training. We show that our approach benefits from applying transfer learning to generalize the solution not only in time but also in space (i.e., various parts of the city). The results are very promising and, unlike competing solutions, our approach utilizes aggregate LTE radio frequency counter data that is inherently privacy-preserving, readily available, and scales globally without any additional network impact.


title: A comparison of three heart rate detection algorithms over ballistocardiogram signals

abstract: Heart rate (HR) detection from ballistocardiogram (BCG) signals is challenging because the signal morphology can vary between and within-subjects. Also, it differs from one sensor to another. Hence, it is essential to evaluate HR detection algorithms across several datasets and under different experimental setups. In this paper, we studied the potential of three HR detection algorithms across four independent BCG datasets. The three algorithms are as follows: the multiresolution analysis of the maximal overlap discrete wavelet transform (MODWT-MRA), continuous wavelet transform (CWT), and template matching (TM). The four datasets were obtained using a microbend fiber optic sensor, a fiber Bragg grating sensor, electromechanical films, and load cells, respectively. The datasets were gathered from: a) 10 patients during a polysomnography study, b) 50 subjects in a sitting position, c) 10 subjects in a sleeping position, and d) 40 subjects in a sleeping position. Overall, CWT with derivative of Gaussian provided superior results compared with the MODWT-MRA, CWT (frequency B-spline), and CWT (Shannon). That said, a BCG template was constructed from DataSet1. Then, it was used for HR detection in the other datasets. The TM method achieved satisfactory results for DataSet2 and DataSet3, but it did not detect the HR of two subjects in DataSet4. The proposed methods were implemented on a Raspberry Pi. As a result, the average time required to analyze a 30-second BCG signal was less than one second for all methods. Yet, the MODWT-MRA had the highest performance with an average time of 0.04 seconds.


title: Intelligent Reflecting Surface Enhanced Multi-UAV NOMA Networks

abstract: Intelligent reflecting surface (IRS) enhanced multi-unmanned aerial vehicle (UAV) non-orthogonal multiple access (NOMA) networks are investigated. A new transmission framework is proposed, where multiple UAV-mounted base stations employ NOMA to serve multiple groups of ground users with the aid of an IRS. The three-dimensional (3D) placement and transmit power of UAVs, the reflection matrix of the IRS, and the NOMA decoding orders among users are jointly optimized for maximization of the sum rate of considered networks. To tackle the formulated mixed-integer non-convex optimization problem with coupled variables, a block coordinate descent (BCD)-based iterative algorithm is developed. Specifically, the original problem is decomposed into three subproblems, which are alternatingly solved by exploiting the penalty method and the successive convex approximation technique. The proposed BCD-based algorithm is demonstrated to be able to obtain a stationary point of the original problem with polynomial time complexity. Numerical results show that: 1) the proposed NOMA-IRS scheme for multi-UAV networks achieves a higher sum rate compared to the benchmark schemes, i.e., orthogonal multiple access (OMA)-IRS and NOMA without IRS; 2) the use of IRS is capable of providing performance gain for multi-UAV networks by both enhancing channel qualities of UAVs to their served users and mitigating the inter-UAV interference; and 3) optimizing the UAV placement can make the sum rate gain brought by NOMA more distinct due to the flexible decoding order design.


title: Formation Rate of Extreme Mass Ratio Inspirals in Active Galactic Nuclei

abstract: Extreme Mass Ratio Inspirals (EMRIs) are important sources for space-borne gravitational wave detectors, such as LISA (Laser Interferometer Space Antenna) and TianQin. Previous EMRI rate studies have focused on the "loss cone" scenario, where stellar-mass black holes (sBHs) are scattered into highly eccentric orbits near the central massive black hole (MBH) via multi-body interaction. In this work, we calculate the rate of EMRIs of an alternative formation channel: EMRI formation assisted by the accretion flow around accreting massive black holes. In this scenario, sBHs and stars on inclined orbits are captured by the accretion disk, and then subsequently migrate towards the MBH, under the influence of density wave generation and head wind. By solving the Fokker-Planck equation incorporating both sBH-sBH/sBH-star scatterings and sBH/star-disk interactions, we find that an accretion disk usually boosts the EMRI formation rate per individual MBH by $\mathcal O(10^1-10^3)$ compared with the canonical "loss cone" formation channel. Taking into account that the fraction of active galactic nucleus (AGNs) is $\sim \mathcal O(10^{-2}-10^{-1})$, where the MBHs are expected to be rapidly accreting, we expect EMRI formation assisted by AGN disks to be an important channel for all EMRIs observed by space-borne gravitational wave detectors. These two channels also predict distinct distributions of EMRI eccentricities and orbit inclinations with respect to the MBH spin equatorial plan, which can be tested by future gravitational wave observations.


title: A superstatistical formulation of complexity measures

abstract: It is discussed how the superstatistical formulation of effective Boltzmann factors can be related to the concept of Kolmogorov complexity, generating an infinite set of complexity measures (CMs) for quantifying information. At this level, the information is treated according to its background, which means that the CM depends on the inherent attributes of the information scenario. While the basic Boltzmann factor directly produces the standard complexity measure (SCM), it succeeds in the description of large-scale scenarios where the data components are not interrelated with themselves, thus adopting the behaviour of a gas. What happens in scenarios in which the presence of sources and sinks of information cannot be neglected, needs of a CM other than the one produced by the ordinary Boltzmann factor. We introduce a set of flexible CMs, without free parameters, that converge asymptotically to the Kolmogorov complexity, but also quantify the information in scenarios with a reasonable small density of states. We prove that these CMs are obtained from a generalised relative entropy and we suggest why such measures are the only compatible generalisations of the SCM.


title: New interpretation of the Extended Geometric Deformation in Isotropic Coordinates

abstract: We study the particular case in which Extended Geometric Deformation does consists of consecutive deformations of temporal and spatial components of the metric, in Schwarzschild-like and isotropic coordinates. In the latter, we present two inequivalent ways to perform this 2-steps GD. This was done in such a way that the method may be applied to different seed solutions. As an example, we use Tolman IV as seed solution, in order to obtain two inequivalent physical solutions with anisotropy in the pressures in Schwarzschild-like coordinates. In the isotropic sector, we obtained four different solutions with anisotropy in the pressures that satisfy physical acceptability conditions, using Gold III as seed solution.


title: Streaming Models for Joint Speech Recognition and Translation

abstract: Using end-to-end models for speech translation (ST) has increasingly been the focus of the ST community. These models condense the previously cascaded systems by directly converting sound waves into translated text. However, cascaded models have the advantage of including automatic speech recognition output, useful for a variety of practical ST systems that often display transcripts to the user alongside the translations. To bridge this gap, recent work has shown initial progress into the feasibility for end-to-end models to produce both of these outputs. However, all previous work has only looked at this problem from the consecutive perspective, leaving uncertainty on whether these approaches are effective in the more challenging streaming setting. We develop an end-to-end streaming ST model based on a re-translation approach and compare against standard cascading approaches. We also introduce a novel inference method for the joint case, interleaving both transcript and translation in generation and removing the need to use separate decoders. Our evaluation across a range of metrics capturing accuracy, latency, and consistency shows that our end-to-end models are statistically similar to cascading models, while having half the number of parameters. We also find that both systems provide strong translation quality at low latency, keeping 99% of consecutive quality at a lag of just under a second.


title: Potential Early Markets for Fusion Energy

abstract: We identify potential early markets for fusion energy and their projected cost targets, based on analysis and synthesis of many relevant, recent studies and reports. Because private fusion companies aspire to start commercial deployment before 2040, we examine cost requirements for fusion-generated electricity, process heat, and hydrogen production based on today's market prices but with various adjustments relating to possible scenarios in 2035, such as "business-as-usual," high renewables penetration, and carbon pricing up to 100 \$/tCO$_2$. Key findings are that fusion developers should consider focusing initially on high-priced global electricity markets and including integrated thermal storage in order to maximize revenue and compete in markets with high renewables penetration. Process heat and hydrogen production will be tough early markets for fusion, but may open up to fusion as markets evolve and if fusion's levelized cost of electricity falls below 50 \$/MWh$_\mathrm{e}$. Finally, we discuss potential ways for a fusion plant to increase revenue via cogeneration (e.g., desalination, direct air capture, or district heating) and to lower capital costs (e.g., by minimizing construction times and interest or by retrofitting coal plants).


title: Probability Analysis and Comparison of Well-Known Integer Factorization Algorithms

abstract: Two prominent methods for integer factorization are those based on general integer sieve and elliptic curve. The general integer sieve method can be specialized to quadratic integer sieve method. In this paper, a probability analysis for the success of these methods is described, under some reasonable conditions. The estimates presented are specialized for the elliptic curve factorization. These methods are compared through heuristic estimates. It is shown that the elliptic curve method is a probabilistic polynomial time algorithm under the assumption of uniform probability distribution for the arising group orders and clearly more likely to succeed, faster asymptotically.


title: Semi-discrete and fully discrete mixed finite element methods for Maxwell viscoelastic model of wave propagation

abstract: Semi-discrete and fully discrete mixed finite element methods are considered for Maxwell-model-based problems of wave propagation in linear viscoelastic solid. This mixed finite element framework allows the use of a large class of existing mixed conforming finite elements for elasticity in the spatial discretization. In the fully discrete scheme, a Crank-Nicolson scheme is adopted for the approximation of the temporal derivatives of stress and velocity variables. Error estimates of the semi-discrete and fully discrete schemes, as well as an unconditional stability result for the fully discrete scheme, are derived. Numerical experiments are provided to verify the theoretical results.


title: QUAREP-LiMi: A community-driven initiative to establish guidelines for quality assessment and reproducibility for instruments and images in light microscopy

abstract: In April 2020, the QUality Assessment and REProducibility for Instruments and Images in Light Microscopy (QUAREP-LiMi) initiative was formed. This initiative comprises imaging scientists from academia and industry who share a common interest in achieving a better understanding of the performance and limitations of microscopes and improved quality control (QC) in light microscopy. The ultimate goal of the QUAREP-LiMi initiative is to establish a set of common QC standards, guidelines, metadata models, and tools, including detailed protocols, with the ultimate aim of improving reproducible advances in scientific research. This White Paper 1) summarizes the major obstacles identified in the field that motivated the launch of the QUAREP-LiMi initiative; 2) identifies the urgent need to address these obstacles in a grassroots manner, through a community of stakeholders including, researchers, imaging scientists, bioimage analysts, bioimage informatics developers, corporate partners, funding agencies, standards organizations, scientific publishers, and observers of such; 3) outlines the current actions of the QUAREP-LiMi initiative, and 4) proposes future steps that can be taken to improve the dissemination and acceptance of the proposed guidelines to manage QC. To summarize, the principal goal of the QUAREP-LiMi initiative is to improve the overall quality and reproducibility of light microscope image data by introducing broadly accepted standard practices and accurately captured image data metrics.


title: Virtual laser scanning with HELIOS++: A novel take on ray tracing-based simulation of topographic 3D laser scanning

abstract: Topographic laser scanning is a remote sensing method to create detailed 3D point cloud representations of the Earth's surface. Since data acquisition is expensive, simulations can complement real data given certain premises are available: i) a model of 3D scene and scanner, ii) a model of the beam-scene interaction, simplified to a computationally feasible while physically realistic level, and iii) an application for which simulated data is fit for use. A number of laser scanning simulators for different purposes exist, which we enrich by presenting HELIOS++. HELIOS++ is an open-source simulation framework for terrestrial static, mobile, UAV-based and airborne laser scanning implemented in C++. The HELIOS++ concept provides a flexible solution for the trade-off between physical accuracy (realism) and computational complexity (runtime, memory footprint), as well as ease of use and of configuration. Unique features of HELIOS++ include the availability of Python bindings (pyhelios) for controlling simulations, and a range of model types for 3D scene representation. HELIOS++ further allows the simulation of beam divergence using a subsampling strategy, and is able to create full-waveform outputs as a basis for detailed analysis. As generation and analysis of waveforms can strongly impact runtimes, the user may set the level of detail for the subsampling, or optionally disable full-waveform output altogether. A detailed assessment of computational considerations and a comparison of HELIOS++ to its predecessor, HELIOS, reveal reduced runtimes by up to 83 %. At the same time, memory requirements are reduced by up to 94 %, allowing for much larger (i.e. more complex) 3D scenes to be loaded into memory and hence to be virtually acquired by laser scanning simulation.


title: Some inequalities of the Edmundson-Lah-Ribaric type for 3-convex functions with applications

abstract: In this paper we derive some Edmundson-Lah-Ribari\v{c} type inequalities for positive linear functionals and 3-convex functions. Main results are applied to the generalized f-divergence functional. Examples with Zipf Mandelbrot law are used to illustrate the results. In addition, obtained results are utilized in constructing some families of exponentially convex functions and Stolarsky-type means.


title: Thermodynamics of light emission

abstract: Some interactions between classical or quantum fields and matter are known to be irreversible processes. Here we associate an entropy to the electromagnetic field from well-known notions of statistical quantum mechanics, in particular the notion of diagonal entropy. We base our work on the study of spontaneous emission and light diffusion. We obtain a quantity which allows to quantify irreversibility for a quantum and classical description of the electromagnetic field, that we can study and interpret from a thermodynamical point of view.


title: The Impact of Multiple Parallel Phrase Suggestions on Email Input and Composition Behaviour of Native and Non-Native English Writers

abstract: We present an in-depth analysis of the impact of multi-word suggestion choices from a neural language model on user behaviour regarding input and text composition in email writing. Our study for the first time compares different numbers of parallel suggestions, and use by native and non-native English writers, to explore a trade-off of "efficiency vs ideation", emerging from recent literature. We built a text editor prototype with a neural language model (GPT-2), refined in a prestudy with 30 people. In an online study (N=156), people composed emails in four conditions (0/1/3/6 parallel suggestions). Our results reveal (1) benefits for ideation, and costs for efficiency, when suggesting multiple phrases; (2) that non-native speakers benefit more from more suggestions; and (3) further insights into behaviour patterns. We discuss implications for research, the design of interactive suggestion systems, and the vision of supporting writers with AI instead of replacing them.


title: SUTRA: A Novel Approach to Modelling Pandemics with Applications to COVID-19

abstract: The Covid-19 pandemic has two key properties: (i) asymptomatic cases (both detected and undetected) that can result in new infections, and (ii) time-varying characteristics due to new variants, Non-Pharmaceutical Interventions etc. We develop a model called SUTRA (Susceptible, Undetected though infected, Tested positive, and Removed Analysis) that takes into account both of these two key properties. While applying the model to a region, two parameters of the model can be learnt from the number of daily new cases found in the region. Using the learnt values of the parameters the model can predict the number of daily new cases so long as the learnt parameters do not change substantially. Whenever any of the two parameters changes due to the key property (ii) above, the SUTRA model can detect that the values of one or both of the parameters have changed. Further, the model has the capability to relearn the changed parameter values, and then use these to carry out the prediction of the trajectory of the pandemic for the region of concern. The SUTRA approach can be applied at various levels of granularity, from an entire country to a district, more specifically, to any large enough region for which the data of daily new cases are available. We have applied the SUTRA model to thirty-two countries, covering more than half of the world's population. Our conclusions are: (i) The model is able to capture the past trajectories very well. Moreover, the parameter values, which we can estimate robustly, help quantify the impact of changes in the pandemic characteristics. (ii) Unless the pandemic characteristics change significantly, the model has good predictive capability. (iii) Natural immunity provides significantly better protection against infection than the currently available vaccines.


title: Misplaced trust? The relationship between trust, ability to identify commercially influenced results, and search engine preference

abstract: People have a high level of trust in search engines, especially Google, but only limited knowledge of them, as numerous studies have shown. This leads to the question: To what extent is this trust justified considering the lack of familiarity among users with how search engines work and the business models they are founded on? We assume that trust in Google, search engine preferences, and knowledge of result types are interrelated. To examine this assumption, we conducted a representative online survey with n = 2,012 German internet users. We show that users with little search engine knowledge are more likely to trust and use Google than users with more knowledge. A contradiction revealed itself - users strongly trust Google, yet they are unable to adequately evaluate search results. This may be problematic since it can potentially affect knowledge acquisition. Consequently, there is a need to promote user information literacy to create a more solid foundation for user trust in search engines.


title: Topology and complexity of the hydrogen bond network in classical models of water

abstract: Over the years, plenty of classical interaction potentials for water have been developed and tested against structural, dynamical and thermodynamic properties. On the other hands, it has been recently observed (F. Martelli et. al, \textit{ACS Nano}, \textbf{14}, 8616--8623, 2020) that the topology of the hydrogen bond network (HBN) is a very sensitive measure that should be considered when developing new interaction potentials. Here we report a thorough comparison of 11 popular non polarizable classical water models against their HBN, which is at the root of water properties. We probe the topology of the HBN using the ring statistics and we evaluate the quality of the network inspecting the percentage of broken and intact HBs. For each water model, we assess the tendency to develop hexagonal rings (that promote crystallization at low temperatures) and pentagonal rings (known to frustrate against crystallization at low temperatures). We then introduce the \emph{network complexity index}, a general descriptor to quantify how much the topology of a given network deviates from that of the ground state, namely of hexagonal or cubic ice. Remarkably, we find that the network complexity index allows us to relate, for the first time, the dynamical properties of different water models with their underlying topology of the HBN. Our study provides a benchmark against which the performances of new models should be tested against, and introduces a general way to quantify the complexity of a network which can be transferred to other materials and that links the topology of the HBN with dynamical properties. Finally, our study introduces a new perspective that can help in rationalizing the transformations among the different phases of water and of other materials.


title: Hardhats and Bungaloos: Comparing Crowdsourced Design Feedback with Peer Design Feedback in the Classroom

abstract: Feedback is an important aspect of design education, and crowdsourcing has emerged as a convenient way to obtain feedback at scale. In this paper, we investigate how crowdsourced design feedback compares to peer design feedback within a design-oriented HCI class and across two metrics: perceived quality and perceived fairness. We also examine the perceived monetary value of crowdsourced feedback, which provides an interesting contrast to the typical requester-centric view of the value of labor on crowdsourcing platforms. Our results reveal that the students (N=106) perceived the crowdsourced design feedback as inferior to peer design feedback in multiple ways. However, they also identified various positive aspects of the online crowds that peers cannot provide. We discuss the meaning of the findings and provide suggestions for teachers in HCI and other researchers interested in crowd feedback systems on using crowds as a potential complement to peers.


title: A Robust Blockchain Readiness Index Model

abstract: As the blockchain ecosystem gets more mature many businesses, investors, and entrepreneurs are seeking opportunities on working with blockchain systems and cryptocurrencies. A critical challenge for these actors is to identify the most suitable environment to start or evolve their businesses. In general, the question is to identify which countries are offering the most suitable conditions to host their blockchain-based activities and implement their innovative projects. The Blockchain Readiness Index (BRI) provides a numerical metric (referred to as the blockchain readiness score) in measuring the maturity/readiness levels of a country in adopting blockchain and cryptocurrencies. In doing so, BRI leverages on techniques from information retrieval to algorithmically derive an index ranking for a set of countries. The index considers a range of indicators organized under five pillars: Government Regulation, Research, Technology, Industry, and User Engagement. In this paper, we further extent BRI with the capability of deriving the index - at the country level - even in the presence of missing information for the indicators. In doing so, we are proposing two weighting schemes namely, linear and sigmoid weighting for refining the initial estimates for the indicator values. A classification framework was employed to evaluate the effectiveness of the developed techniques which yielded to a significant classification accuracy.


title: The Next Decade of Telecommunications Artificial Intelligence

abstract: It has been an exciting journey since the mobile communications and artificial intelligence were conceived 37 years and 64 years ago. While both fields evolved independently and profoundly changed communications and computing industries, the rapid convergence of 5G and deep learning is beginning to significantly transform the core communication infrastructure, network management and vertical applications. The paper first outlines the individual roadmaps of mobile communications and artificial intelligence in the early stage, with a concentration to review the era from 3G to 5G when AI and mobile communications started to converge. With regard to telecommunications artificial intelligence, the paper further introduces in detail the progress of artificial intelligence in the ecosystem of mobile communications. The paper then summarizes the classifications of AI in telecom ecosystems along with its evolution paths specified by various international telecommunications standardization bodies. Towards the next decade, the paper forecasts the prospective roadmap of telecommunications artificial intelligence. In line with 3GPP and ITU-R timeline of 5G & 6G, the paper further explores the network intelligence following 3GPP and ORAN routes respectively, experience and intention driven network management and operation, network AI signalling system, intelligent middle-office based BSS, intelligent customer experience management and policy control driven by BSS and OSS convergence, evolution from SLA to ELA, and intelligent private network for verticals. The paper is concluded with the vision that AI will reshape the future B5G or 6G landscape and we need pivot our R&D, standardizations, and ecosystem to fully take the unprecedented opportunities.


title: Gas cooling of test masses for future gravitational-wave observatories

abstract: Recent observations made with Advanced LIGO and Advanced Virgo have initiated the era of gravitational-wave astronomy. The number of events detected by these "2nd Generation" (2G) ground-based observatories is partially limited by noise arising from temperature-induced position fluctuations of the test mass mirror surfaces used for probing spacetime dynamics. The design of next-generation gravitational-wave observatories addresses this limitation by using cryogenically cooled test masses; current approaches for continuously removing heat (resulting from absorbed laser light) rely on heat extraction via black-body radiation or conduction through suspension fibres. As a complementing approach for extracting heat during observational runs, we investigate cooling via helium gas impinging on the test mass in free molecular flow. We establish a relation between cooling power and corresponding displacement noise, based on analytical models, which we compare to numerical simulations. Applying this theoretical framework with regard to the conceptual design of the Einstein Telescope (ET), we find a cooling power of 10 mW at 18 K for a gas pressure that exceeds the ET design strain noise goal by at most a factor of $\sim 3$ in the signal frequency band from 3 to 11 Hz. A cooling power of 100 mW at 18 K corresponds to a gas pressure that exceeds the ET design strain noise goal by at most a factor of $\sim 11$ in the band from 1 to 28 Hz.


title: Recovery of the Order of Derivation for Fractional Diffusion Equations in an Unknown Medium

abstract: In this work, we investigate the recovery of a parameter in a diffusion process given by the order of derivation in time for a class of diffusion type equations, including both classical and time-fractional diffusion equations, from the flux measurement observed at one point on the boundary. The mathematical model for time-fractional diffusion equations involves a Djrbashian-Caputo fractional derivative in time. We prove a uniqueness result in an unknown medium (e.g., diffusion coefficients, obstacle, initial condition and source), i.e., the recovery of the order of derivation in a diffusion process having several pieces of unknown information. The proof relies on the analyticity of the solution at large time, asymptotic decay behavior, strong maximum principle of the elliptic problem and suitable application of the Hopf lemma. Further we provide an easy-to-implement reconstruction algorithm based on a nonlinear least-squares formulation, and several numerical experiments are presented to complement the theoretical analysis.


title: Stability criterion for linear systems of ordinary differential equations

abstract: The Riccati equation method is used to establish a new stability criteria for linear systems of ordinary differential equations. Two examples are presented in which the obtained result is compared with the results obtained by the Lyapunov and Bogdanov methods, by a method involving estimates of solutions in the Lozinskii logarithmic norms and by the freezing method.


title: Improved Sensitivity of Base Layer on the Performance of Rigid Pavement

abstract: The performance of rigid pavement is greatly affected by the properties of base/subbase as well as subgrade layer. However, the performance predicted by the AASHTOWare Pavement ME design shows low sensitivity to the properties of base and subgrade layers. To improve the sensitivity and better reflect the influence of unbound layers a new set of improved models i.e., resilient modulus (MR) and modulus of subgrade reaction (k-value) are adopted in this study. An Artificial Neural Network (ANN) model is developed to predict the modified k-value based on finite element (FE) analysis. The training and validation datasets in the ANN model consist of 27000 simulation cases with different combinations of pavement layer thickness, layer modulus and slab-base interface bond ratio. To examine the sensitivity of modified MR and k-values on pavement response, eight pavement sections data are collected from the Long-Term Pavement performance (LTPP) database and modeled by using the FE software ISLAB2000. The computational results indicate that the modified MR values have higher sensitivity to water content in base layer on critical stress and deflection response of rigid pavements compared to the results using the Pavement ME design model. It is also observed that the k-values using ANN model has the capability of predicting critical pavement response at any partially bonded conditions whereas the Pavement ME design model can only calculate at two extreme bonding conditions (i.e., fully bonding and no bonding).


title: Linear confinement of generalized KG-oscillator with a uniform magnetic field in Kaluza-Klein theory and Aharonov-Bohm effect

abstract: In this paper, we solve generalized KG-oscillator interacts with a uniform magnetic field in five-dimensional space-time background produced by topological defects under a linear confining potential using the Kaluza-Klein theory. We solve this equation and analyze an analogue of the Aharonov-Bohm effect for bound states. We observe that the energy levels for each radial mode depends on the global parameters characterizing the space-time, the confining potential, and the magnetic field which shows a quantum effect


title: The two-loop perturbative correction to the (g-2)$_{\mu}$ HLbL at short distances

abstract: The short-distance behaviour of the hadronic light-by-light (HLbL) contribution to $(g-2)_{\mu}$ has recently been studied by means of an operator product expansion in a background electromagnetic field. The leading term in this expansion has been shown to be given by the massless quark loop, and the non-perturbative corrections are numerically very suppressed. Here, we calculate the perturbative QCD correction to the massless quark loop. The correction is found to be fairly small compared to the quark loop as far as we study energy scales where the perturbative running for the QCD coupling is well-defined, i.e.~for scales $\mu\gtrsim 1\, \mathrm{GeV}$. This should allow to reduce the large systematic uncertainty associated to high-multiplicity hadronic states.


title: Why are there Seven Sisters?

abstract: There are two puzzles surrounding the Pleiades, or Seven Sisters. First, why are the mythological stories surrounding them, typically involving seven young girls being chased by a man associated with the constellation Orion, so similar in vastly separated cultures, such as the Australian Aboriginal cultures and Greek mythology? Second, why do most cultures call them "Seven Sisters" even though most people with good eyesight see only six stars? Here we show that both these puzzles may be explained by a combination of the great antiquity of the stories combined with the proper motion of the stars, and that these stories may predate the departure of most modern humans out of Africa around 100,000 BC.


title: Bit Commitment in Operational Probabilistic Theories

abstract: The aim of this thesis project is to investigate the bit commitment protocol in the framework of operational probabilistic theories. In particular a careful study is carried on the feasibility of bit commitment in the non-local boxes theory. New aspects of the theory are also presented.


title: The $L^{2}$ sequential convergence of a solution to the mass-critical NLS above the ground state

abstract: In this paper we generalize a weak sequential result of \cite{fan20182} to a non-scattering solutions in dimension $d \geq 2$. No symmetry assumptions are required for the initial data. We build on a previous result of \cite{dodson20202} for one dimension.


title: Deconfining Phase Boundary of Rapidly Rotating Hot and Dense Matter and Analysis of Moment of Inertia

abstract: We discuss the effect of rapid rotation on the phase diagram of hadronic matter. The energy dispersion relation is shifted by an effective chemical potential induced by rotation. This suggests that rotation should lower the critical temperature of chiral restoration, but it is still controversial how the deconfinement temperature should change as a function of angular velocity. We adopt the hadron resonance gas model as an approach free from fitting parameters. We identify the deconfinement from the thermodynamic behavior and find that rotation decreases the deconfinement temperature. We also discuss the spatial inhomogeneity of the pressure and give a semi-quantitative estimate of the moment of inertia.


title: Sparsistent filtering of comovement networks from high-dimensional data

abstract: Network filtering is an important form of dimension reduction to isolate the core constituents of large and interconnected complex systems. We introduce a new technique to filter large dimensional networks arising out of dynamical behavior of the constituent nodes, exploiting their spectral properties. As opposed to the well known network filters that rely on preserving key topological properties of the realized network, our method treats the spectrum as the fundamental object and preserves spectral properties. Applying asymptotic theory for high dimensional data for the filter, we show that it can be tuned to interpolate between zero filtering to maximal filtering that induces sparsity and consistency while having the least spectral distance from a linear shrinkage estimator. We apply our proposed filter to covariance networks constructed from financial data, to extract the key subnetwork embedded in the full sample network.


title: "FISTA" in Banach spaces with adaptive discretisations

abstract: FISTA is a popular convex optimisation algorithm which is known to converge at an optimal rate whenever a minimiser is contained in a suitable Hilbert space. We propose a modified algorithm where each iteration is performed in a subset which is allowed to change at every iteration. Sufficient conditions are provided for guaranteed convergence, although at a reduced rate depending on the conditioning of the specific problem. These conditions have a natural interpretation when a minimiser exists in an underlying Banach space. Typical examples are L1-penalised reconstructions where we provide detailed theoretical and numerical analysis.


title: Understanding Visual Saliency in Mobile User Interfaces

abstract: For graphical user interface (UI) design, it is important to understand what attracts visual attention. While previous work on saliency has focused on desktop and web-based UIs, mobile app UIs differ from these in several respects. We present findings from a controlled study with 30 participants and 193 mobile UIs. The results speak to a role of expectations in guiding where users look at. Strong bias toward the top-left corner of the display, text, and images was evident, while bottom-up features such as color or size affected saliency less. Classic, parameter-free saliency models showed a weak fit with the data, and data-driven models improved significantly when trained specifically on this dataset (e.g., NSS rose from 0.66 to 0.84). We also release the first annotated dataset for investigating visual saliency in mobile UIs.


title: X-ray Scatter Estimation Using Deep Splines

abstract: Algorithmic X-ray scatter compensation is a desirable technique in flat-panel X-ray imaging and cone-beam computed tomography. State-of-the-art U-net based image translation approaches yielded promising results. As there are no physics constraints applied to the output of the U-Net, it cannot be ruled out that it yields spurious results. Unfortunately, those may be misleading in the context of medical imaging. To overcome this problem, we propose to embed B-splines as a known operator into neural networks. This inherently limits their predictions to well-behaved and smooth functions. In a study using synthetic head and thorax data as well as real thorax phantom data, we found that our approach performed on par with U-net when comparing both algorithms based on quantitative performance metrics. However, our approach not only reduces runtime and parameter complexity, but we also found it much more robust to unseen noise levels. While the U-net responded with visible artifacts, our approach preserved the X-ray signal's frequency characteristics.


title: Estimating $\alpha$-Rank by Maximizing Information Gain

abstract: Game theory has been increasingly applied in settings where the game is not known outright, but has to be estimated by sampling. For example, meta-games that arise in multi-agent evaluation can only be accessed by running a succession of expensive experiments that may involve simultaneous deployment of several agents. In this paper, we focus on $\alpha$-rank, a popular game-theoretic solution concept designed to perform well in such scenarios. We aim to estimate the $\alpha$-rank of the game using as few samples as possible. Our algorithm maximizes information gain between an epistemic belief over the $\alpha$-ranks and the observed payoff. This approach has two main benefits. First, it allows us to focus our sampling on the entries that matter the most for identifying the $\alpha$-rank. Second, the Bayesian formulation provides a facility to build in modeling assumptions by using a prior over game payoffs. We show the benefits of using information gain as compared to the confidence interval criterion of ResponseGraphUCB (Rowland et al. 2019), and provide theoretical results justifying our method.


title: Contact-mediated cellular communication supplements positional information to regulate spatial patterning during development

abstract: Development in multi-cellular organisms is marked by a high degree of spatial organization of the cells attaining distinct fates in the embryo. We show that receptor-ligand interaction between cells in close physical proximity adaptively regulates the local process of selective gene expression in the presence of a global field set up by a diffusing morphogen that provides positional cues. This allows information from the cellular neighborhood to be incorporated into the emergent thresholds of morphogen concentration that dictate cell fate, consistent with recent experiments.


title: A Newton's Iteration Converges Quadratically to Nonisolated Solutions Too

abstract: The textbook Newton's iteration is practically inapplicable on solutions of nonlinear systems with singular Jacobians. By a simple modification, a novel extension of Newton's iteration regains its local quadratic convergence toward nonisolated solutions that are semiregular as properly defined regardless of whether the system is square, underdetermined or overdetermined while Jacobians can be rank-deficient. Furthermore, the iteration serves as a regularization mechanism for computing singular solutions from empirical data. When a system is perturbed, its nonisolated solutions can be altered substantially or even disappear. The iteration still locally converges to a stationary point that approximates a singular solution of the underlying system with an error bound in the same order of the data accuracy. Geometrically, the iteration approximately approaches the nearest point on the solution manifold. The method simplifies the modeling of nonlinear systems by permitting nonisolated solutions and enables a wide range of applications in algebraic computation.


title: Approximation capability of two hidden layer feedforward neural networks with fixed weights

abstract: We algorithmically construct a two hidden layer feedforward neural network (TLFN) model with the weights fixed as the unit coordinate vectors of the $d$-dimensional Euclidean space and having $3d+2$ number of hidden neurons in total, which can approximate any continuous $d$-variable function with an arbitrary precision. This result, in particular, shows an advantage of the TLFN model over the single hidden layer feedforward neural network (SLFN) model, since SLFNs with fixed weights do not have the capability of approximating multivariate functions.


title: On the negativity of the Wigner function as a measure of entanglement under quantum polarization converter devices

abstract: We study the behaviour of the Negativity of Wigner Function (NWF) as a measure of entanglement in non-Gaussian states under quantum polarisation converter devices. We analyze comparatively this quantity with other measures of entanglement in a system prepared in a superposition of two-mode coherent states. We show that the (WF) can be identified as a quantifier of non-Gaussian entanglement.


title: The extended Bregman divergence and parametric estimation

abstract: Minimization of suitable statistical distances~(between the data and model densities) has proved to be a very useful technique in the field of robust inference. Apart from the class of $\phi$-divergences of \cite{a} and \cite{b}, the Bregman divergence (\cite{c}) has been extensively used for this purpose. However, since the data density must have a linear presence in the cross product term of the Bregman divergence involving both the data and model densities, several useful divergences cannot be captured by the usual Bregman form. In this respect, we provide an extension of the ordinary Bregman divergence by considering an exponent of the density function as the argument rather than the density function itself. We demonstrate that many useful divergence families, which are not ordinarily Bregman divergences, can be accommodated within this extended description. Using this formulation, one can develop many new families of divergences which may be useful in robust inference. In particular, through an application of this extension, we propose the new class of the GSB divergence family. We explore the applicability of the minimum GSB divergence estimator in discrete parametric models. Simulation studies as well as conforming real data examples are given to demonstrate the performance of the estimator and to substantiate the theory developed.


title: Tensor-Train Networks for Learning Predictive Modeling of Multidimensional Data

abstract: In this work, we firstly apply the Train-Tensor (TT) networks to construct a compact representation of the classical Multilayer Perceptron, representing a reduction of up to 95% of the coefficients. A comparative analysis between tensor model and standard multilayer neural networks is also carried out in the context of prediction of the Mackey-Glass noisy chaotic time series and NASDAQ index. We show that the weights of a multidimensional regression model can be learned by means of TT network and the optimization of TT weights is a more robust to the impact of coefficient initialization and hyper-parameter setting. Furthermore, an efficient algorithm based on alternating least squares has been proposed for approximating the weights in TT-format with a reduction of computational calculus, providing a much faster convergence than the well-known adaptive learning-method algorithms, widely applied for optimizing neural networks.


title: Structural and transport properties of 4f electron doped Y1-x(Dy)xPdBi topological semi-metallic thin films

abstract: We report the effect of 4f electron doping on structural, electrical and magneto-transport properties of Dy doped half Heusler Y1-x(Dy)xPdBi (x =0, 0.2, 0.5, 1) thin films grown by pulsed laser deposition. The Dy doping leads to lattice contraction which increases from 0% for the parent x =0 sample to approx 1.3% for x=1 sample with increase in Dy doping. The electrical transport measurements show a typical semi-metallic behaviour in the temperature range 3K to 300K and a sharp drop in resistivity at low temperatures (less than 3K) for all the samples. Magnetotransport measurements and Shubnikov de-Hass oscillations at high magnetic fields demonstrate that for these topologically non-trivial samples, Dy doping induced lattice contraction plays an active role in modifying the Fermi surface, carrier concentration and the effective electron mass. There is an uniform suppression of the onset of superconductivity with increased Dy doping which is possibly related to the increasing local exchange field arising from the 4f electrons in Dy. Our results indicate that we can tune various band structure parameters of YPdBi by f electron doping and strained thin films of Y1-x(Dy)xPdBi show surface dominated relativistic carrier transport at low temperatures.


title: Comment on "Strong Quantum Darwinism and Strong Independence are Equivalent to Spectrum Broadcast Structure"

abstract: In a recent Letter [Phys. Rev. Lett. 122, 010403 (2019)], an equivalence is proposed between the so-called Spectrum Broadcast Structure for a system-multienvironment quantum state, and the conjunction of two information-theory notions: (a) Strong Quantum Darwinism; and (b) Strong Independence. Here, we show that the mathematical formulation of condition (b) by the authors (namely, the pairwise independence of the fragments of the environment, conditioned on the system), is necessary but not sufficient to ensure the equivalence. We propose a simple counter-example, together with a strengthened formulation of condition (b), ensuring the equivalence proposed by the authors.


title: Third-harmonic light polarization control in magnetically-resonant silicon metasurfaces

abstract: Nonlinear metasurfaces have become prominent tools for controlling and engineering light at the nanoscale. Usually, the polarization of the total generated third harmonic is studied. However, diffraction orders may present different polarizations. Here, we design an high quality factor silicon metasurface for third harmonic generation and perform back focal plane imaging of the diffraction orders, which present a rich variety of polarization states. Our results demonstrate the possibility of tailoring the polarization of the generated nonlinear diffraction orders paving the way to a higher degree of wavefront control.


title: Super-Spreaders Out, Super-Spreading In: The Effects of Infectiousness Heterogeneity and Lockdowns on Herd Immunity

abstract: Recently, [8] has proposed that heterogeneity of infectiousness (and susceptibility) across individuals in infectious diseases, plays a major role in affecting the Herd Immunity Threshold (HIT). Such heterogeneity has been observed in COVID-19 and is recognized as overdispersion (or "super-spreading"). The model of [8] suggests that super-spreaders contribute significantly to the effective reproduction factor, R, and that they are likely to get infected and immune early in the process. Consequently, under R_0 = 3 (attributed to COVID-19), the Herd Immunity Threshold (HIT) is as low as 5%, in contrast to 67% according to the traditional models [1, 2, 4, 10]. This work follows up on [8] and proposes that heterogeneity of infectiousness (susceptibility) has two "faces" whose mix affects dramatically the HIT: (1) Personal-Trait-, and (2) Event-Based- Infectiousness (Susceptibility). The former is a personal trait of specific individuals (super-spreaders) and is nullified once those individuals are immune (as in [8]). The latter is event-based (e.g cultural super-spreading events) and remains effective throughout the process, even after the super-spreaders immune. We extend [8]'s model to account for these two factors, analyze it and conclude that the HIT is very sensitive to the mix between (1) and (2), and under R_0 = 3 it can vary between 5% and 67%. Preliminary data from COVID-19 suggests that herd immunity is not reached at 5%. We address operational aspects and analyze the effects of lockdown strategies on the spread of a disease. We find that herd immunity (and HIT) is very sensitive to the lockdown type. While some lockdowns affect positively the disease blocking and increase herd immunity, others have adverse effects and reduce the herd immunity.


title: Analytical Model for the Current Density in the Electrochemical Synthesis of Porous Silicon Structures with a Lateral Gradient

abstract: Layered optical devices with a lateral gradient can be fabricated through electrochemical synthesis of porous silicon (PS) using a position dependent etching current density $\bm j(\bm r_\|)$. Predicting the local value of $\bm j(\bm r_\|)$ and the corresponding porosity $p(\bm r_\|)$ and etching rate $v(\bm r_\|)$ is desirable for their systematic design. We develop a simple analytical model for the calculation of $\bm j(\bm r_\|)$ within a prism shaped cell. Graded single layer PS samples were synthesized and their local calibration curves $p$ vs $\bm j$ and $v$ vs $\bm j$ were obtained from our model and their reflectance spectra. The agreement found between the calibration curves from different samples shows that from one sample we could obtain full calibration curves which may be used to predict, design, and fabricate more complex non-homogeneous multilayered devices with lateral gradients for manifold applications.


title: Correcting ray distortion in tomographic additive manufacturing

abstract: Light-based additive manufacturing techniques enable a rapid transition from object design to production. In these approaches, a 3D object is typically built by successive polymerization of 2D layers in a photocurable resin. A recently demonstrated technique, however, uses tomographic dose patterning to establish a 3D light dose distribution within a cylindrical glass vial of photoresin. Lensing distortion from the cylindrical vial is currently mitigated by either an index matching bath around the print volume or a cylindrical lens. In this work, we show that these hardware approaches to distortion correction are unnecessary. Instead, we demonstrate how the lensing effect can be computationally corrected by resampling the parallel-beam radon transform into an aberrated geometry. We also demonstrate a more general application of our computational approach by correcting for non-telecentricity inherent in most optical projection systems. We expect that our results will underpin a more simple and flexible class of tomographic 3D printers where deviations from the assumed parallel-beam projection geometry are rectified computationally.


title: Discrete phase space and continuous time relativistic quantum mechanics II: Peano circles, hyper-tori phase cells, and fibre bundles

abstract: The discrete phase space and continuous time representation of relativistic quantum mechanics is further investigated here as a continuation of paper I [1]. The main mathematical construct used here will be that of an area-filling Peano curve. We show that the limit of a sequence of a class of Peano curves is a Peano circle denoted as $\bar{S}^{1}_{n}$, a circle of radius $\sqrt{2n+1}$ where $n \in \{0,1,\cdots\}$. We interpret this two-dimensional Peano circle in our framework as a phase cell inside our two-dimensional discrete phase plane. We postulate that a first quantized Planck oscillator, being very light, and small beyond current experimental detection, occupies this phase cell $\bar{S}^{1}_{n}$. The time evolution of this Peano circle sweeps out a two-dimensional vertical cylinder analogous to the world-sheet of string theory. Extending this to three dimensional space, we introduce a $(2+2+2)$-dimensional phase space hyper-tori $\bar{S}^{1}_{n^1} \times \bar{S}^{1}_{n^2} \times \bar{S}^{1}_{n^3}$ as the appropriate phase cell in the physical dimensional discrete phase space. A geometric interpretation of this structure in state space is given in terms of product fibre bundles. We also study free scalar Bosons in the background $[(2+2+2)+1]$-dimensional discrete phase space and continuous time state space using the relativistic partial difference-differential Klein-Gordon equation. The second quantized field quantas of this system can cohabit with the tiny Planck oscillators inside the $\bar{S}^{1}_{n^1} \times \bar{S}^{1}_{n^2} \times \bar{S}^{1}_{n^3}$ phase cells for eternity. Finally, a generalized free second quantized Klein-Gordon equation in a higher $[(2+2+2)N+1]$-dimensional discrete state space is explored. The resulting discrete phase space dimension is compared to the significant spatial dimensions of some of the popular models of string theory.


title: Gravity Optimizer: a Kinematic Approach on Optimization in Deep Learning

abstract: We introduce Gravity, another algorithm for gradient-based optimization. In this paper, we explain how our novel idea change parameters to reduce the deep learning model's loss. It has three intuitive hyper-parameters that the best values for them are proposed. Also, we propose an alternative to moving average. To compare the performance of the Gravity optimizer with two common optimizers, Adam and RMSProp, five standard datasets were trained on two VGGNet models with a batch size of 128 for 100 epochs. Gravity hyper-parameters did not need to be tuned for different models. As will be explained more in the paper, to investigate the direct impact of the optimizer itself on loss reduction no overfitting prevention technique was used. The obtained results show that the Gravity optimizer has more stable performance than Adam and RMSProp and gives greater values of validation accuracy for datasets with more output classes like CIFAR-100 (Fine).


title: Dense outlier detection and open-set recognition based on training with noisy negative images

abstract: Deep convolutional models often produce inadequate predictions for inputs foreign to the training distribution. Consequently, the problem of detecting outlier images has recently been receiving a lot of attention. Unlike most previous work, we address this problem in the dense prediction context in order to be able to locate outlier objects in front of in-distribution background. Our approach is based on two reasonable assumptions. First, we assume that the inlier dataset is related to some narrow application field (e.g.~road driving). Second, we assume that there exists a general-purpose dataset which is much more diverse than the inlier dataset (e.g.~ImageNet-1k). We consider pixels from the general-purpose dataset as noisy negative training samples since most (but not all) of them are outliers. We encourage the model to recognize borders between known and unknown by pasting jittered negative patches over inlier training images. Our experiments target two dense open-set recognition benchmarks (WildDash 1 and Fishyscapes) and one dense open-set recognition dataset (StreetHazard). Extensive performance evaluation indicates competitive potential of the proposed approach.


title: It Takes Two to Tango: Combining Visual and Textual Information for Detecting Duplicate Video-Based Bug Reports

abstract: When a bug manifests in a user-facing application, it is likely to be exposed through the graphical user interface (GUI). Given the importance of visual information to the process of identifying and understanding such bugs, users are increasingly making use of screenshots and screen-recordings as a means to report issues to developers. However, when such information is reported en masse, such as during crowd-sourced testing, managing these artifacts can be a time-consuming process. As the reporting of screen-recordings in particular becomes more popular, developers are likely to face challenges related to manually identifying videos that depict duplicate bugs. Due to their graphical nature, screen-recordings present challenges for automated analysis that preclude the use of current duplicate bug report detection techniques. To overcome these challenges and aid developers in this task, this paper presents Tango, a duplicate detection technique that operates purely on video-based bug reports by leveraging both visual and textual information. Tango combines tailored computer vision techniques, optical character recognition, and text retrieval. We evaluated multiple configurations of Tango in a comprehensive empirical evaluation on 4,860 duplicate detection tasks that involved a total of 180 screen-recordings from six Android apps. Additionally, we conducted a user study investigating the effort required for developers to manually detect duplicate video-based bug reports and compared this to the effort required to use Tango. The results reveal that Tango's optimal configuration is highly effective at detecting duplicate video-based bug reports, accurately ranking target duplicate videos in the top-2 returned results in 83% of the tasks. Additionally, our user study shows that, on average, Tango can reduce developer effort by over 60%, illustrating its practicality.


title: Randomization Inference beyond the Sharp Null: Bounded Null Hypotheses and Quantiles of Individual Treatment Effects

abstract: Randomization inference (RI) is typically interpreted as testing Fisher's "sharp" null hypothesis that all unit-level effects are exactly zero. This hypothesis is often criticized as restrictive and implausible, making its rejection scientifically uninteresting. We show, however, that many randomization tests are also valid for a "bounded" null hypothesis under which the unit-level effects are all non-positive (or all non-negative) but are otherwise heterogeneous. In addition to being more plausible a priori, bounded nulls are closely related to substantively important concepts such as monotonicity and Pareto efficiency. Reinterpreting RI in this way also dramatically expands the range of inferences possible in this framework. We show that exact confidence intervals for the maximum (or minimum) unit-level effect can be obtained by inverting tests for a sequence of bounded nulls. We also generalize RI to cover inference for quantiles of the individual effect distribution as well as for the proportion of individual effects larger (or smaller) than a given threshold. The proposed confidence intervals for all effect quantiles are simultaneously valid, in the sense that no correction for multiple analyses is required, and are thus a "free lunch" added to conventional RI. In sum, our reinterpretation and generalization provide a broader justification for randomization tests and a basis for exact nonparametric inference for effect quantiles. We illustrate our methods with simulations and applications, finding that Stephenson rank statistics can provide more informative results than the more common Wilcoxon rank or difference-in-means statistics. We also provide an R package RIQITE implementing the proposed approach.


title: Sharp $\left( H_{p},L_{p}\right) $ type inequalities of maximal operators of $T$ means with respect to Vilenkin systems with monotone coefficients

abstract: In this paper we prove and discuss some new $\left( H_{p},L_{p}\right)$ type inequalities of maximal operators of $T$ means with respect to the Vilenkin systems with monotone coefficients. We also apply these inequalities to prove strong convergence theorems of such $T$ means. We also show that these results are the best possible in a special sense. As applications, both some well-known and new results are pointed out.


title: Flexible estimation of the state dwell-time distribution in hidden semi-Markov models

abstract: Hidden semi-Markov models generalise hidden Markov models by explicitly modelling the time spent in a given state, the so-called dwell time, using some distribution defined on the natural numbers. While the (shifted) Poisson and negative binomial distribution provide natural choices for such distributions, in practice, parametric distributions can lack the flexibility to adequately model the dwell times. To overcome this problem, a penalised maximum likelihood approach is proposed that allows for a flexible and data-driven estimation of the dwell-time distributions without the need to make any distributional assumption. This approach is suitable for direct modelling purposes or as an exploratory tool to investigate the latent state dynamics. The feasibility and potential of the suggested approach is illustrated by modelling muskox movements in northeast Greenland using GPS tracking data. The proposed method is implemented in the R-package PHSMM which is available on CRAN.


title: Assessment of turbulence-chemistry interaction models in MILD combustion regime

abstract: The present paper reports on the assessment of different turbulence-chemistry interaction closures for modeling turbulent combustion in the Moderate and Intense Low oxygen Dilution combustion regime. 2D RANS simulations have been carried out to model flames issuing from two burners DJHC burner and Adelaide JHC burner which imitate the MILD combustion. In order to model these flames, two different approaches of turbulence-chemistry interaction models are considered; while in the PDF based modeling, two different variants are invoked to understand the applicability of the PDF based models in the MILD regime: one is based on presumed shape PDF approach and the other one is transported PDF approach. A comprehensive study has been carried out by comparing the results obtained from these different models. For the DJHC burner, the computations are carried out for a jet speed corresponding to Reynolds numbers of Re=4100, whereas the Adelaide JHC burner computations are performed for a jet speed corresponding to Reynolds number of Re=10000. The effects of molecular diffusion on the flame characteristics are also studied by using different micro-mixing models. In the case of DJHC burner, it has been observed that the mean axial velocity and the turbulent kinetic energy profiles are in good agreement with the measurements. However, the temperature profiles are over-predicted in the downstream region by both EDC and the PDF based models. In the context of Adelaide JHC burner, the profiles of the temperature and the mass fraction of major species (CH4, H2, O2, H2O, CO2) obtained using LPDF approach are in better agreement with the measurements compared to those obtained using EDC model; although, both the solution approaches fail to capture CO and OH radical profiles.


title: Large Eddy Simulation of MILD combustion using PDF based turbulence-chemistry interaction models

abstract: The present work reports on LES of turbulent reacting flow of the Delft-Jet-in-Hot-Coflow (DJHC) burner, emulating MILD combustion, with transported PDF based combustion models using ANSYS FLUENT 13.0. Two different eddy viscosity models for LES (Dynamic Smagorinsky and Kinetic Energy Transport) along with two solution approaches for PDF transport equation, i.e. Eulerian and Lagrangian, have been used in the present study. Moreover, the effects of chemical kinetics and the micro-mixing models have also been investigated for two different fuel jet Reynolds number (Re = 4,100 and Re = 8,800). The mean velocity and turbulent kinetic energy predicted by the different models are in good agreement with experimental data. Both the composition PDF models predict an early ignition resulting in higher radial mean temperature predictions at burner exit. The models, however, correctly predict the formation mechanism of ignition kernels and the decreasing trend of the lift-off height with increasing jet Reynolds number, as observed experimentally.


title: Large Eddy Simulation of a Premixed Bunsen flame using a modified Thickened-Flame model at two Reynolds number

abstract: A modified Thickened Flame (TF) model based on Large Eddy Simulation methodology is used to investigate premixed combustion and the model predictions are evaluated by comparing with the piloted premixed stoichiometric methane-air flame data (Chen et al., 1996) for Reynolds numbers Re = 24,000 (flame F3) and Re=52,000 (flame F1). The basic idea of Thickened-Flame approach is that the flame front is artificially thickened to resolve on the computational LES grid while keeping the laminar flame speed ( ) constant. The artificially thickening of the flame front is obtained by enhancing the molecular diffusion and decreasing the pre-exponential factor of the Arrhenius law. Since the flame front is artificially thickened, the response of the thickened flame to turbulence is affected and taken care of by incorporating an efficiency function (E) in the governing equations. The efficiency function (E) in the modified TF model is proposed based on the direct numerical simulations (DNS) data set of flame-vortex interactions (Colin et al., 2000). The predicted simulation results are compared with the experimental data and with computations reported using RANS based probability distribution function (PDF) modeling approach (Lindstedt, R. P. and Vaos, E. M., 2006, Transported PDF modeling of high-Reynolds-number premixed turbulent flames, Combustion and Flame, 145, 495) and RANS based G-equation approach (Herrmann, M., 2006). It is shown that the results with the modified TF model are generally in good agreement with the data, with the TF predictions consistently comparable to the PDF model predictions and superior to the results with the G-equation approach.


title: Mass sensing for the advanced fabrication of nanomechanical resonators

abstract: We report on a nanomechanical engineering method to monitor matter growth in real time via e-beam electromechanical coupling. This method relies on the exceptional mass sensing capabilities of nanomechanical resonators. Focused electron beam induced deposition (FEBID) is employed to selectively grow platinum particles at the free end of singly clamped nanotube cantilevers. The electron beam has two functions: it allows both to grow material on the nanotube and to track in real time the deposited mass by probing the noise-driven mechanical resonance of the nanotube. On the one hand, this detection method is highly effective as it can resolve mass deposition with a resolution in the zeptogram range; on the other hand, this method is simple to use and readily available to a wide range of potential users, since it can be operated in existing commercial FEBID systems without making any modification. The presented method allows to engineer hybrid nanomechanical resonators with precisely tailored functionality. It also appears as a new tool for studying growth dynamics of ultra-thin nanostructures, opening new opportunities for investigating so far out-of-reach physics of FEBID and related methods.


title: Joint Identification through Hybrid Models Improved by Correlations

abstract: In mechanical systems coupled with joints, accurate prediction of the joint characteristics is extremely important. Despite years of research, a lot is yet to be learnt about the joints' interface dynamics. The problem becomes even more difficult when the interface Degrees-of-Freedom (DoF) are inaccessible for Frequency Response Function (FRF) measurements. This is, for example, the case of bladed-disk systems with dove-tail or fir-tree type joints. Therefore, an FRF based expansion method called System Equivalent Model Mixing (SEMM) is used to obtain expanded interface dynamics. The method uses numerical and experimental sub-models of each component and their assembly to produce the respective expanded or hybrid sub-models. By applying substructure decoupling to these sub-models, the joint can be identified. However, the joint can be noisy due to expansion and measurement errors which propagate to the hybrid sub-models. In this paper, a correlation based approach is proposed in the SEMM method wherein the quality of the expanded sub-models is improved. In this new approach, several expanded models are generated systematically using different combinations of the experimental FRFs and computing a parameter, Frequency Response Assurance Criteria (FRAC), to evaluate quality of the contribution of the different measurements. The lowest correlated channels or FRFs can be filtered out based on a certain threshold value of FRAC. Using the improved hybrid sub-models, the joint identification also shows a remarkable improvement. The test object for the method is an assembly of disk and one blade with a dove-tail joint.


title: Coordinate conditions and field equations for pure composite gravity

abstract: Whenever an alternative theory of gravity is formulated in a background Minkowski space, the conditions characterizing admissible coordinate systems, in which the alternative theory of gravity may be applied, play an important role. We here propose Lorentz covariant coordinate conditions for the composite theory of pure gravity developed from the Yang-Mills theory based on the Lorentz group, thereby completing this previously proposed higher derivative theory of gravity. The physically relevant static isotropic solutions are determined by various methods, the high-precision predictions of general relativity are reproduced, and an exact black-hole solution with mildly singular behavior is found.


title: Episodic accretion in focus: revealing the environment of FU Orionis-type stars

abstract: The earliest phases of star formation are characterised by intense mass accretion from the circumstellar disk to the central star. One group of young stellar objects, the FU Orionis-type stars exhibit accretion rate peaks accompanied by bright eruptions. The occurrence of these outbursts might solve the luminosity problem of protostars, play a key role in accumulating the final star mass, and have a significant effect on the parameters of the envelope and the disk. In the framework of the Structured Accretion Disks ERC project, we are conducting a systematic investigation of these sources with millimeter interferometry to examine whether they represent normal young stars in exceptional times or they are unusual objects. Our results show that FU Orionis-type stars can be similar to both Class I and Class II systems and may be in a special evolutionary phase between the two classes with their infall-driven episodic eruptions being the main driving force of the transition.


title: Equivariant and supersymmetric localization in QFT

abstract: Equivariant localization theory is a powerful tool that has been extensively used in the past thirty years to elegantly obtain exact integration formulas, in both mathematics and physics. These integration formulas are proved within the mathematical formalism of equivariant cohomology, a variant of standard cohomology theory that incorporates the presence of a symmetry group acting on the space at hand. A suitable infinite-dimensional generalization of this formalism is applicable to a certain class of Quantum Field Theories (QFT) endowed with supersymmetry. In this thesis we review the formalism of equivariant localization and some of its applications in Quantum Mechanics (QM) and QFT. We start from the mathematical description of equivariant cohomology and related localization theorems of finite-dimensional integrals in the case of an Abelian group action, and then we discuss their formal application to infinite-dimensional path integrals in QFT. We summarize some examples from the literature of computations of partition functions and expectation values of supersymmetric operators in various dimensions. For 1-dimensional QFT, that is QM, we review the application of the localization principle to the derivation of the Atiyah-Singer index theorem applied to the Dirac operator on a twisted spinor bundle. In 3 and 4 dimensions, we examine the computation of expectation values of certain Wilson loops in supersymmetric gauge theories and their relation to 0-dimensional theories described by "matrix models". Finally, we review the formalism of non-Abelian localization applied to 2-dimensional Yang-Mills theory and its application in the mapping between the standard "physical" theory and a related "cohomological" formulation.


title: Global Protoplanetary Disk Simulations: Dead Zone Formation and FUor Outbursts

abstract: We conducted global hydrodynamic simulations of protoplanetary disk evolution with an adaptive Shakura-Sunyaev {\alpha} prescription to represent the layered disk structure, and starting with the collapse phase of the molecular cloud. With the canonical values of model parameters, self-consistent dead zones formed at the scale of a few au. The instabilities associated with the dead zone and corresponding outbursts, similar to FUor eruptions, were also observed in the simulations.


title: Differentiable Trust Region Layers for Deep Reinforcement Learning

abstract: Trust region methods are a popular tool in reinforcement learning as they yield robust policy updates in continuous and discrete action spaces. However, enforcing such trust regions in deep reinforcement learning is difficult. Hence, many approaches, such as Trust Region Policy Optimization (TRPO) and Proximal Policy Optimization (PPO), are based on approximations. Due to those approximations, they violate the constraints or fail to find the optimal solution within the trust region. Moreover, they are difficult to implement, often lack sufficient exploration, and have been shown to depend on seemingly unrelated implementation choices. In this work, we propose differentiable neural network layers to enforce trust regions for deep Gaussian policies via closed-form projections. Unlike existing methods, those layers formalize trust regions for each state individually and can complement existing reinforcement learning algorithms. We derive trust region projections based on the Kullback-Leibler divergence, the Wasserstein L2 distance, and the Frobenius norm for Gaussian distributions. We empirically demonstrate that those projection layers achieve similar or better results than existing methods while being almost agnostic to specific implementation choices. The code is available at https://git.io/Jthb0.


title: Some applications of transversality for infinite dimensional manifolds

abstract: We present some transversality results for a category of Fr\'{e}chet manifolds, the so-called $MC^k$-Fr\'{e}chet manifolds. In this context, we apply the obtained transversality results to construct the degree of nonlinear Fredholm mappings by virtue of which we prove a rank theorem, an invariance of domain theorem and a Bursuk-Ulam type theorem.


title: The radial structure of planetary bodies formed by the streaming instability

abstract: Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global and lack the resolution to follow the final collapse stage of a pebble cloud within its Hill radius. We aim to track the collapse of a gravitationally bound pebble cloud subject to mutual collisions and gas drag with the representative particle approach. We determine the radial pebble size distribution of the collapsed core and the impact of mutual pebble collisions on the pebble size distribution. We find that virial equilibrium is never reached during the cloud evolution and that, in general, pebbles with given Stokes number (St) collapse towards an optically thick core in a sequence from aerodynamically largest to aerodynamically smallest. We show that at the location for which the core becomes optically thick, the terminal velocity is well below the fragmentation threshold velocity. While collisional processing is negligible during cloud evolution, the collisions that do occur are sticking. These results support the observations that comets and small planetary bodies are composed of primordial pebbles in the milimeter to centimeter size range