title: Yukawa vs. Newton: gravitational forces in a cubic cosmological simulation box

abstract: We study the behaviour of Yukawa and Newtonian gravitational forces in a cubic box with fully periodic boundaries commonly encountered in N-body simulations of the structure formation. Placing a single gravitating body at the origin of coordinates, we reveal the scales at which non-negligible deviation from the Yukawa law occurs when the Newtonian approximation is employed. We discuss the results in terms of the corresponding physical distances today as well as earlier, back at the matter-dominated stage. Revisiting the problem for free boundaries, we also compare the periodic and plain gravitational forces for Yukawa-type interactions.


title: Constraining the SMEFT with a differential cross section measurement of $tWZ$ production at the HL-LHC

abstract: A prospective measurement of the differential cross section of \tWZ production with respect to the transverse momentum of the \PZ boson using a general-purpose detector at the High-Luminosity Large Hadron Collider (HL-LHC) is described. The response of a general-purpose detector at the HL-LHC is simulated and used to estimate the uncertainties and covariances of the differential cross section measurement. Constraints on the Standard Model Effective Field Theory (SMEFT) enabled by the measurement are estimated. A parametric model of the differential cross section in the SMEFT is constructed and is used to determine the expected posterior probability function of six SMEFT Wilson coefficients and the expected 95\% Bayesian credible intervals for each coefficient and pair of coefficients. The intervals suggest that for all coefficients, the measurement will provide competitive but weaker constraints than those derived from other \mbox{HL-LHC} measurements involving top quarks and \PZ bosons. However, as the measurement is simultaneously sensitive to to a unique set of SMEFT coefficients, it will provide a useful input to a global SMEFT analysis that considers many operators.


title: The LEGA-C and SAMI Galaxy Surveys: Quiescent Stellar Populations and the Mass-Size Plane across 6 Gyr

abstract: We investigate the change in mean stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift ($0.60\leq z\leq0.76$) using the LEGA-C Survey, to low redshift ($0.014\leq z\leq0.10$) using the SAMI Galaxy Survey. We find that, similarly to their low-redshift counterparts, the stellar metallicity of quiescent galaxies at $0.60\leq z\leq 0.76$ closely correlates with $M_*/R_\mathrm{e}$ (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals, by determining the escape velocity required by metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. On the other hand, we find no correlation between stellar age and $M_*/R_\mathrm{e}^2$ (stellar mass surface density $\Sigma$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$\Sigma$ relation in the context of the redshift evolution of the star-forming and quiescent populations in the mass--size plane, and find our results can be explained as a consequence of galaxies forming more compactly at higher redshifts, and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The $z\sim 0$ age--$\Sigma$ relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and so a range of surface densities.


title: Atomic structure optimization with machine-learning enabled interpolation between chemical elements

abstract: We introduce a computational method for global optimization of structure and ordering in atomic systems. The method relies on interpolation between chemical elements, which is incorporated in a machine learning structural fingerprint. The method is based on Bayesian optimization with Gaussian processes and is applied to the global optimization of Au-Cu bulk systems, Cu-Ni surfaces with CO adsorption, and Cu-Ni clusters. The method consistently identifies low-energy structures, which are likely to be the global minima of the energy. For the investigated systems with 23-66 atoms, the number of required energy and force calculations is in the range 3-75.


title: Structural biology in the clouds: The WeNMR-EOSC Ecosystem

abstract: Structural biology aims at characterizing the structural and dynamic properties of biological macromolecules at atomic details. Gaining insight into three dimensional structures of biomolecules and their interactions is critical for understanding the vast majority of cellular processes, with direct applications in health and food sciences. Since 2010, the WeNMR project (www.wenmr.eu) has implemented numerous web-based services to facilitate the use of advanced computational tools by researchers in the field, using the high throughput computing infrastructure provided by EGI. These services have been further developed in subsequent initiatives under H2020 projects and are now operating as Thematic Services in the European Open Science Cloud (EOSC) portal (www.eosc-portal.eu), sending >12 millions of jobs and using around 4000 CPU-years per year. Here we review 10 years of successful e-infrastructure solutions serving a large worldwide community of over 23,000 users to date, providing them with user-friendly, web-based solutions that run complex workflows in structural biology. The current set of active WeNMR portals are described, together with the complex backend machinery that allows distributed computing resources to be harvested efficiently.


title: Backward-Compatible Prediction Updates: A Probabilistic Approach

abstract: When machine learning systems meet real world applications, accuracy is only one of several requirements. In this paper, we assay a complementary perspective originating from the increasing availability of pre-trained and regularly improving state-of-the-art models. While new improved models develop at a fast pace, downstream tasks vary more slowly or stay constant. Assume that we have a large unlabelled data set for which we want to maintain accurate predictions. Whenever a new and presumably better ML models becomes available, we encounter two problems: (i) given a limited budget, which data points should be re-evaluated using the new model?; and (ii) if the new predictions differ from the current ones, should we update? Problem (i) is about compute cost, which matters for very large data sets and models. Problem (ii) is about maintaining consistency of the predictions, which can be highly relevant for downstream applications; our demand is to avoid negative flips, i.e., changing correct to incorrect predictions. In this paper, we formalize the Prediction Update Problem and present an efficient probabilistic approach as answer to the above questions. In extensive experiments on standard classification benchmark data sets, we show that our method outperforms alternative strategies along key metrics for backward-compatible prediction updates.


title: Einstein-Podolsky-Rosen uncertainty limits for bipartite multimode states

abstract: Certification and quantification of correlations for multipartite states of quantum systems appear to be a central task in quantum information theory. We give here a unitary quantum-mechanical perspective of both entanglement and Einstein-Podolsky-Rosen (EPR) steering of continuous-variable multimode states. This originates in the Heisenberg uncertainty relations for the canonical quadrature operators of the modes. Correlations of two-party $(N\, \text{vs} \,1)$-mode states are examined by using the variances of a pair of suitable EPR-like observables. It turns out that the uncertainty sum of these nonlocal variables is bounded from below by local uncertainties and is strengthened differently for separable states and for each one-way unsteerable ones. The analysis of the minimal properly normalized sums of these variances yields necessary conditions of separability and EPR unsteerability of $(N\, \text{vs} \,1)$-mode states in both possible ways of steering. When the states and the performed measurements are Gaussian, then these conditions are precisely the previously-known criteria of separability and one-way unsteerability.


title: Constraining the NuMI neutrino flux using inverse muon decay reactions in MINERvA

abstract: Inverse muon decay, $\nu_\mu e^-\to\mu^-\nu_e$, is a reaction whose cross-section can be predicted with very small uncertainties. It has a neutrino energy threshold of $\approx 11$ GeV and can be used to constrain the high-energy part of the flux in the NuMI neutrino beam. This reaction is the dominant source of events which only contain high-energy muons nearly parallel to the direction of the neutrino beam. We have isolated a sample of hundreds of such events in neutrino and anti-neutrino enhanced beams, and have constrained the predicted high-energy flux.


title: Projected Least-Squares Quantum Process Tomography

abstract: We propose and investigate a new method of quantum process tomography (QPT) which we call projected least squares (PLS). In short, PLS consists of first computing the least-squares estimator of the Choi matrix of an unknown channel, and subsequently projecting it onto the convex set of Choi matrices. We consider four experimental setups including direct QPT with Pauli eigenvectors as input and Pauli measurements, and ancilla-assisted QPT with mutually unbiased bases (MUB) measurements. In each case, we provide a closed form solution for the least-squares estimator of the Choi matrix. We propose a novel, two-step method for projecting these estimators onto the set of matrices representing physical quantum channels, and a fast numerical implementation in the form of the hyperplane intersection projection algorithm. We provide rigorous, non-asymptotic concentration bounds, sampling complexities and confidence regions for the Frobenius and trace-norm error of the estimators. For the Frobenius error, the bounds are linear in the rank of the Choi matrix, and for low ranks, they improve the error rates of the least squares estimator by a factor $d^2$, where $d$ is the system dimension. We illustrate the method with numerical experiments involving channels on systems with up to 7 qubits, and find that PLS has highly competitive accuracy and computational tractability.


title: Ultradifferentiable extension theorems: a survey

abstract: We survey ultradifferentiable extension theorems, i.e., quantitative versions of Whitney's classical extension theorem, with special emphasis on the existence of continuous linear extension operators. The focus is on Denjoy-Carleman classes for which we develop the theory from scratch and discuss important related concepts such as (non-)quasianalyticity. It allows us to give an efficient and, to a fair extent, elementary introduction to Braun-Meise-Taylor classes based on their representation as intersections and unions of Denjoy-Carleman classes.


title: Two-phase geothermal model with fracture network and multi-branch wells

abstract: This paper focuses on the numerical simulation of geothermal systems in complex geological settings. The physical model is based on two-phase Darcy flows coupling the mass conservation of the water component with the energy conservation and the liquid vapor thermodynamical equilibrium. The discretization exploits the flexibility of unstructured meshes to model complex geology including conductive faults as well as complex wells. The polytopal and essentially nodal Vertex Approximate Gradient scheme is used for the approximation of the Darcy and Fourier fluxes combined with a Control Volume approach for the transport of mass and energy. Particular attention is paid to the faults which are modelled as two-dimensional interfaces defined as collection of faces of the mesh and to the flow inside deviated or multi-branch wells defined as collection of edges of the mesh with rooted tree data structure. By using an explicit pressure drop calculation, the well model reduces to a single equation based on complementarity constraints with only one well implicit unknown. The coupled systems are solved fully implicitely at each time step using efficient nonlinear and linear solvers on parallel distributed architectures. The convergence of the discrete model is investigated numerically on a simple test case with a Cartesian geometry and a single vertical producer well. Then, the ability of our approach to deal efficiently with realistic test cases is assessed on a high energy faulted geothermal reservoir operated using a doublet of two deviated wells.


title: Magnification-independent Histopathological Image Classification with Similarity-based Multi-scale Embeddings

abstract: The classification of histopathological images is of great value in both cancer diagnosis and pathological studies. However, multiple reasons, such as variations caused by magnification factors and class imbalance, make it a challenging task where conventional methods that learn from image-label datasets perform unsatisfactorily in many cases. We observe that tumours of the same class often share common morphological patterns. To exploit this fact, we propose an approach that learns similarity-based multi-scale embeddings (SMSE) for magnification-independent histopathological image classification. In particular, a pair loss and a triplet loss are leveraged to learn similarity-based embeddings from image pairs or image triplets. The learned embeddings provide accurate measurements of similarities between images, which are regarded as a more effective form of representation for histopathological morphology than normal image features. Furthermore, in order to ensure the generated models are magnification-independent, images acquired at different magnification factors are simultaneously fed to networks during training for learning multi-scale embeddings. In addition to the SMSE, to eliminate the impact of class imbalance, instead of using the hard sample mining strategy that intuitively discards some easy samples, we introduce a new reinforced focal loss to simultaneously punish hard misclassified samples while suppressing easy well-classified samples. Experimental results show that the SMSE improves the performance for histopathological image classification tasks for both breast and liver cancers by a large margin compared to previous methods. In particular, the SMSE achieves the best performance on the BreakHis benchmark with an improvement ranging from 5% to 18% compared to previous methods using traditional features.


title: Gamma-ray Observation of the Cygnus Region in the 100 TeV Energy Region

abstract: We report observations of gamma-ray emissions with energies in the 100 TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1 with the pulsar moving away from its original birthplace situated around the centroid of the observed gamma-ray emission. This work would stimulate further studies of particle acceleration mechanisms at these gamma-ray sources.


title: Reverse Sensitivity Analysis for Risk Modelling

abstract: We consider the problem where a modeller conducts sensitivity analysis of a model consisting of random input factors, a corresponding random output of interest, and a baseline probability measure. The modeller seeks to understand how the model (the distribution of the input factors as well as the output) changes under a stress on the output's distribution. Specifically, for a stress on the output random variable, we derive the unique stressed distribution of the output that is closest in the Wasserstein distance to the baseline output's distribution and satisfies the stress. We further derive the stressed model, including the stressed distribution of the inputs, which can be calculated in a numerically efficient way from a set of baseline Monte Carlo samples and which is implemented in the R package SWIM on CRAN. The proposed reverse sensitivity analysis framework is model-free and allows for stresses on the output such as (a) the mean and variance, (b) any distortion risk measure including the Value-at-Risk and Expected-Shortfall, and (c) expected utility type constraints, thus making the reverse sensitivity analysis framework suitable for risk models.


title: Guided flux motion: models and experiment

abstract: The article provides a brief overview of existing models describing guided flux motion. The first model was proposed by Niessen et al. This model works in the single-vortex approximation and provides qualitative explanation of the effect, but in most cases it gives an overestimated guiding angle. The stochastic model explains the experimentally observed effect of decreasing the guiding angle, the so-called slipping effect, by the influence of thermal fluctuations. However, the performed estimates show that thermal fluctuations should not be so significant. An alternative guiding model is the anisotropic pinning model. This model works in the critical state approximation, and the slipping effect is the result of the combined action of anisotropic pinning and vortex interaction. The predictions of all three models are verified by the example of a superconducting Nb-Ti tape containing about 6 vol. % {/alpha}-Ti, which acts as strong pinning centers. Cutting samples at different angles to rolling allows one to control the direction of the driving force. It was found that when the standard criterion of the electric field is reached, the guiding angle varies significantly at different places along the sample, which is explained by the plastic mode of the vortex matter motion. In this case, the anisotropic pinning model semi-quantitatively predicts the direction angle averaged over the sample length. With an increase in the driving force, the vortex matter motion passes into an elastic mode, and the guiding angle becomes uniform along the sample.


title: NB-IoT via LEO satellites: An efficient resource allocation strategy for uplink data transmission

abstract: In this paper, we focus on the use of Low-Eart Orbit (LEO) satellites providing the Narrowband Internet of Things (NB-IoT) connectivity to the on-ground user equipment (UEs). Conventional resource allocation algorithms for the NBIoT systems are particularly designed for terrestrial infrastructures, where devices are under the coverage of a specific base station and the whole system varies very slowly in time. The existing methods in the literature cannot be applied over LEO satellite-based NB-IoT systems for several reasons. First, with the movement of the LEO satellite, the corresponding channel parameters for each user will quickly change over time. Delaying the scheduling of a certain user would result in a resource allocation based on outdated parameters. Second, the differential Doppler shift, which is a typical impairment in communications over LEO, directly depends on the relative distance among users. Scheduling at the same radio frame users that overcome a certain distance would violate the differential Doppler limit supported by the NB-IoT standard. Third, the propagation delay over a LEO satellite channel is around 4-16 times higher compared to a terrestrial system, imposing the need for message exchange minimization between the users and the base station. In this work, we propose a novel uplink resource allocation strategy that jointly incorporates the new design considerations previously mentioned together with the distinct channel conditions, satellite coverage times and data demands of various users on Earth. The novel methodology proposed in this paper can act as a framework for future works in the field.


title: Unsupervised Spoken Utterance Classification

abstract: An intelligent virtual assistant (IVA) enables effortless conversations in call routing through spoken utterance classification (SUC) which is a special form of spoken language understanding (SLU). Building a SUC system requires a large amount of supervised in-domain data that is not always available. In this paper, we introduce an unsupervised spoken utterance classification approach (USUC) that does not require any in-domain data except for the intent labels and a few para-phrases per intent. USUC is consisting of a KNN classifier (K=1) and a complex embedding model trained on a large amount of unsupervised customer service corpus. Among all embedding models, we demonstrate that Elmo works best for USUC. However, an Elmo model is too slow to be used at run-time for call routing. To resolve this issue, first, we compute the uni- and bi-gram embedding vectors offline and we build a lookup table of n-grams and their corresponding embedding vector. Then we use this table to compute sentence embedding vectors at run-time, along with back-off techniques for unseen n-grams. Experiments show that USUC outperforms the traditional utterance classification methods by reducing the classification error rate from 32.9% to 27.0% without requiring supervised data. Moreover, our lookup and back-off technique increases the processing speed from 16 utterances per second to 118 utterances per second.


title: Solving the subset sum problem with a nonideal biological computer

abstract: We consider the solution of the subset sum problem based on a parallel computer consisting of self-propelled biological agents moving in a nanostructured network that encodes the NP-complete task in its geometry. We develop an approximate analytical method to analyze the effects of small errors in the nonideal junctions composing the computing network by using a Gaussian confidence interval approximation of the multinomial distribution. We concretely evaluate the probability distribution for error-induced paths and determine the minimal number of agents required to obtain a proper solution. We finally validate our theoretical results with exact numerical simulations of the subset sum problem for different set sizes and error probabilities.


title: Homogeneity in the instrument-treatment association is not sufficient for the Wald estimand to equal the average causal effect for a binary instrument and a continuous exposure

abstract: Background: Interpreting instrumental variable results often requires further assumptions in addition to the core assumptions of relevance, independence, and the exclusion restriction. Methods: We assess whether instrument-exposure additive homogeneity renders the Wald estimand equal to the average derivative effect (ADE) in the case of a binary instrument and a continuous exposure. Results: Instrument-exposure additive homogeneity is insufficient for ADE identification when the instrument is binary, the exposure is continuous and the effect of the exposure on the outcome is non-linear on the additive scale. For a binary exposure, the exposure-outcome effect is necessarily additive linear, so the homogeneity condition is sufficient. Conclusions: For binary instruments, instrument-exposure additive homogeneity identifies the ADE if the exposure is also binary. Otherwise, additional assumptions (such as additive linearity of the exposure-outcome effect) are required.


title: Asymptotically Safe Gravity-Fermion systems on curved backgrounds

abstract: We set up a consistent background field formalism for studying the renormalization group (RG) flow of gravity coupled to $N_f$ Dirac fermions on maximally symmetric backgrounds. Based on Wetterich's equation we perform a detailed study of the resulting fixed point structure in a projection including the Einstein-Hilbert action, the fermion anomalous dimension, and a specific coupling of the fermion bilinears to the spacetime curvature. The latter constitutes a mass-type term which breaks chiral symmetry explicitly. Our analysis identifies two infinite families of interacting RG fixed points which are viable candidates to provide a high-energy completion through the asymptotic safety mechanism. The fixed points exist for all values of $N_f$ outside of a small window situated at low values $N_f$ and become weakly coupled in the large $N_f$-limit. Symmetry-wise, they correspond to "quasi-chiral" and "non-chiral" fixed points. The former come with enhanced predictive power, fixing one of the couplings via the asymptotic safety condition. Moreover, the interplay of the fixed points allows for cross-overs from the non-chiral to the chiral fixed point, giving a dynamical mechanism for restoring the symmetry approximately at intermediate scales. Our discussion of chiral symmetry breaking effects provides strong indications that the topology of spacetime plays a crucial role when analyzing whether quantum gravity admits light chiral fermions.


title: Mapping Lamb, Stark and Purcell effects at a chromophore-picocavity junction with hyper-resolved fluorescence microscopy

abstract: The interactions between the excited states of a single chromophore with static and dynamic electric fields confined to a plasmonic cavity of picometer dimensions are investigated in a joint experimental and theoretical effort. In this configuration, the spatial extensions of the confined fields are smaller than the one of the molecular exciton, a property that is used to generate fluorescence maps of the chromophores with intra-molecular resolution. Theoretical simulations of the electrostatic and electrodynamic interactions occurring at the chromophore-picocavity junction are able to reproduce and interpret these hyper-resolved fluorescence maps, and reveal the key role played by subtle variations of Purcell, Lamb and Stark effects at the chromophore-picocavity junction.


title: 3d gravity in Bondi-Weyl gauge: charges, corners, and integrability

abstract: We introduce a new gauge and solution space for three-dimensional gravity. As its name Bondi-Weyl suggests, it leads to non-trivial Weyl charges, and uses Bondi-like coordinates to allow for an arbitrary cosmological constant and therefore spacetimes which are asymptotically locally (A)dS or flat. We explain how integrability requires a choice of integrable slicing and also the introduction of a corner term. After discussing the holographic renormalization of the action and of the symplectic potential, we show that the charges are finite, symplectic and integrable, yet not conserved. We find four towers of charges forming an algebroid given by $\mathfrak{vir}\oplus\mathfrak{vir}\oplus\text{Heisenberg}$ with three central extensions, where the base space is parametrized by the retarded time. These four charges generate diffeomorphisms of the boundary cylinder, Weyl rescalings of the boundary metric, and radial translations. We perform this study both in metric and triad variables, and use the triad to explain the covariant origin of the corner terms needed for renormalization and integrability.


title: SPH modelling of wind-companion interactions in eccentric AGB binary systems

abstract: The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularly, there is a lack of studies investigating the structure formation in eccentric systems. We aim to improve our understanding of the wind morphology of eccentric AGB binary systems by investigating the mechanism responsible for the different small-scale structures and global morphologies that arise in a polytropic wind with different velocities. Using the smoothed particle hydrodynamics (SPH) code Phantom, we generate nine different high-resolution, 3D simulations of an AGB star with a solar-mass companion with various wind velocity and eccentricity combinations. The models assume a polytropic gas, with no additional cooling. We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, that is flattened. In the case of a lower wind velocity, the stronger interaction results in the formation of a high-energy region and bow-shock structure that can shape the wind into an irregular morphology if instabilities arise. High-eccentricity models show a complex, phase-dependent interaction leading to wind structures that are irregular in three dimensions. However, the significant interaction with the companion compresses matter into an equatorial density enhancement, irrespective of eccentricity.


title: Quantum speed of evolution in a Markovian bosonic environment

abstract: We present explicit evaluations of quantum speed limit times pertinent to the Markovian dynamics of an open continuous-variable system. Specifically, we consider the standard setting of a cavity mode of the quantum radiation field weakly coupled to a thermal bosonic reservoir. The evolution of the field state is ruled by the quantum optical master equation, which is known to have an exact analytic solution. Starting from a pure input state, we employ two indicators of how different the initial and evolved states are, namely, the fidelity of evolution and the Hilbert-Schmidt distance of evolution. The former was introduced by del Campo {\em et al.} who derived a time-independent speed limit for the evolution of a Markovian open system. We evaluate it for this field-reservoir setting, with an arbitrary input pure state of the field mode. The resultant formula is then specialized to the coherent and Fock states. On the other hand, we exploit an alternative approach that employs both indicators of evolution mentioned above. Their rates of change have the same upper bound, and consequently provide a unique time-dependent quantum speed limit. It turns out that the associate quantum speed limit time built with the Hilbert-Schmidt metric is tighter than the fidelity-based one. As apposite applications, we investigate the damping of the coherent and Fock states by using the characteristic functions of the corresponding evolved states. General expressions of both the fidelity and the Hilbert-Schmidt distance of evolution are obtained and analyzed for these two classes of input states. In the case of a coherent state, we derive accurate formulas for their common speed limit and the pair of associate limit times.


title: DUKweb: Diachronic word representations from the UK Web Archive corpus

abstract: Lexical semantic change (detecting shifts in the meaning and usage of words) is an important task for social and cultural studies as well as for Natural Language Processing applications. Diachronic word embeddings (time-sensitive vector representations of words that preserve their meaning) have become the standard resource for this task. However, given the significant computational resources needed for their generation, very few resources exist that make diachronic word embeddings available to the scientific community. In this paper we present DUKweb, a set of large-scale resources designed for the diachronic analysis of contemporary English. DUKweb was created from the JISC UK Web Domain Dataset (1996-2013), a very large archive which collects resources from the Internet Archive that were hosted on domains ending in `.uk'. DUKweb consists of a series word co-occurrence matrices and two types of word embeddings for each year in the JISC UK Web Domain dataset. We show the reuse potential of DUKweb and its quality standards via a case study on word meaning change detection.


title: Numerical convergence of discrete extensions in a space-time finite element, fictitious domain method for the Navier-Stokes equations

abstract: A key ingredient of our fictitious domain, higher order space-time cut finite element (CutFEM) approach for solving the incompressible Navier--Stokes equations on evolving domains (cf.\ \cite{Bause2021}) is the extension of the physical solution from the time-dependent flow domain $\Omega_f^t$ to the entire, time-independent computational domain $\Omega$. The extension is defined implicitly and, simultaneously, aims at stabilizing the discrete solution in the case of unavoidable irregular small cuts. Here, the convergence properties of the scheme are studied numerically for variations of the combined extension and stabilization.


title: General Board Game Concepts

abstract: Many games often share common ideas or aspects between them, such as their rules, controls, or playing area. However, in the context of General Game Playing (GGP) for board games, this area remains under-explored. We propose to formalise the notion of "game concept", inspired by terms generally used by game players and designers. Through the Ludii General Game System, we describe concepts for several levels of abstraction, such as the game itself, the moves played, or the states reached. This new GGP feature associated with the ludeme representation of games opens many new lines of research. The creation of a hyper-agent selector, the transfer of AI learning between games, or explaining AI techniques using game terms, can all be facilitated by the use of game concepts. Other applications which can benefit from game concepts are also discussed, such as the generation of plausible reconstructed rules for incomplete ancient games, or the implementation of a board game recommender system.


title: Cooperative Training and Latent Space Data Augmentation for Robust Medical Image Segmentation

abstract: Deep learning-based segmentation methods are vulnerable to unforeseen data distribution shifts during deployment, e.g. change of image appearances or contrasts caused by different scanners, unexpected imaging artifacts etc. In this paper, we present a cooperative framework for training image segmentation models and a latent space augmentation method for generating hard examples. Both contributions improve model generalization and robustness with limited data. The cooperative training framework consists of a fast-thinking network (FTN) and a slow-thinking network (STN). The FTN learns decoupled image features and shape features for image reconstruction and segmentation tasks. The STN learns shape priors for segmentation correction and refinement. The two networks are trained in a cooperative manner. The latent space augmentation generates challenging examples for training by masking the decoupled latent space in both channel-wise and spatial-wise manners. We performed extensive experiments on public cardiac imaging datasets. Using only 10 subjects from a single site for training, we demonstrated improved cross-site segmentation performance and increased robustness against various unforeseen imaging artifacts compared to strong baseline methods. Particularly, cooperative training with latent space data augmentation yields 15% improvement in terms of average Dice score when compared to a standard training method.


title: Impact of $B \to K \nu \bar \nu$ measurements on beyond the Standard Model theories

abstract: Semileptonic flavor changing neutral current transitions with a pair of neutrinos in the final state are very accurately determined in the standard model (SM) and thus provide an accurate and sensitive probe for physics beyond the SM. Until recently, the poor tagging efficiency for the $B\to K^{(*)}\nu \bar{\nu}$ modes made them less advantageous as a probe of new physics (NP) compared to the charged lepton counterparts. The most recent Belle II result on $B\to K \nu \bar{\nu}$ uses an innovative inclusive tagging technique resulting in a higher tagging efficiency; this together with previous BaBar and Belle results indicates a possible enhancement in the branching fraction of $B^+\to K^+ \nu \bar{\nu}$. A reanalysis of the full Belle dataset together with upcoming Belle II dataset is expected to result in a much more precise measurement of this mode. If the branching ratio is indeed found to be enhanced with improved measurements, this would provide an unambiguous signal of NP without uncertainties due to long-distance non-factorizable effects or power corrections (in contrast to $B\to K^{(*)} \ell \ell$). We have explored the possibilities of such an enhancement as a signal of NP within several scenarios, which can also explain some of the other tensions observed in neutral as well as charged current $B$-decays. In an effective field theory approach, with the most general dimension-six Hamiltonian including light right-handed neutrinos, we explore the viability of all scalar and vector leptoquarks as well as the parameter space possible with a generic vector gauge boson $Z^\prime$ model assuming minimal new particle content. While being consistent with all data, correlations between the observed intriguing discrepancies in $B$-decays are also obtained, which will discriminate between the various NP scenarios.


title: Neural Network Layer Algebra: A Framework to Measure Capacity and Compression in Deep Learning

abstract: We present a new framework to measure the intrinsic properties of (deep) neural networks. While we focus on convolutional networks, our framework can be extrapolated to any network architecture. In particular, we evaluate two network properties, namely, capacity, which is related to expressivity, and compression, which is related to learnability. Both these properties depend only on the network structure and are independent of the network parameters. To this end, we propose two metrics: the first one, called layer complexity, captures the architectural complexity of any network layer; and, the second one, called layer intrinsic power, encodes how data is compressed along the network. The metrics are based on the concept of layer algebra, which is also introduced in this paper. This concept is based on the idea that the global properties depend on the network topology, and the leaf nodes of any neural network can be approximated using local transfer functions, thereby, allowing a simple computation of the global metrics. We show that our global complexity metric can be calculated and represented more conveniently than the widely-used VC dimension. We also compare the properties of various state-of-the art architectures using our metrics and use the properties to analyze their accuracy on benchmark image classification datasets.


title: Identifying Processes Governing Damage Evolution in Quasi-Static Elasticity Part I -- Analysis

abstract: We present a quasi-static elasticity model that accounts for damage evolution based on the ideas of Kachanov 1958 and and Rabotnov 1968. We analyze the resulting strongly nonlinear system of differential equations in view of well-posedness. The specific feature is that the displacements are connected to the damage evolution equation via a Nemytskii- or superposition- operator. The novelty in this work is that we present an inverse problem in a parameter identification setting, in which we are able to identify the shape of this Nemytskii-operator. From the material modelling point of view, the relation of material damage and displacements is modelled by the shape of this operator. We establish the Fr\'echet-derivative of the forward operator as well as the adjoint of the derivative and characterize both via systems of linear differential equations. We prove ill-posedness of the inverse problem and provide a sufficient condition for the classical nonlinear Landweber method to converge.


title: A covariant approach to relativistic large-eddy simulations: The fibration picture

abstract: Models of turbulent flows require the resolution of a vast range of scales, from large eddies to small-scale features directly associated with dissipation. As the required resolution is not within reach of large scale numerical simulations, standard strategies involve a smoothing of the fluid dynamics, either through time averaging or spatial filtering. These strategies raise formal issues in general relativity, where the split between space and time is observer dependent. To make progress, we develop a new covariant framework for filtering/averaging based on the fibration of spacetime associated with fluid elements and the use of Fermi coordinates to facilitate a meaningful local analysis. We derive the resolved equations of motion, demonstrating how "effective" dissipative terms arise because of the coarse-graining, and paying particular attention to the thermodynamical interpretation of the resolved quantities. Finally, as the smoothing of the fluid-dynamics inevitably leads to a closure problem, we propose a new closure scheme inspired by recent progress in the modelling of dissipative relativistic fluids, and crucially, demonstrate the linear stability of the proposed model.


title: A chemical nano-reactor based on a levitated nanoparticle in vacuum

abstract: A single levitated nanoparticle is used as a nano-reactor for studying surface chemistry at the nanoscale. Optical levitation under controlled pressure, surrounding gas composition, and humidity provides extreme control over the nanoparticle, including dynamics, charge, and surface chemistry. Using a single nanoparticle avoids ensemble averages and allows to study how the presence of silanol groups at its surface affects the adsorption and desorption of water from the background gas with unprecedented real time, spatial, and temporal resolution. Here, we demonstrate the unique potential of this versatile platform by studying the Zhuravlev model in silica particles. In contrast to standard methods, our system allowed the first observation of an abrupt and irreversible change in scattering cross section, mass, and mechanical eigenfrequency during the dehydroxylation process, indicating changes in density, refractive index and volume.


title: Static output feedback stabilization of uncertain rational nonlinear systems with input saturation

abstract: In this paper, the notion of robust strict QSR-dissipativity is applied to solve the static output feedback control problem for a class of continuous-time nonlinear rational systems subject to input saturation and bounded parametric uncertainties. A local dissipativity condition is combined with generalized sector conditions to formulate the synthesis of a stabilizing controller in terms of linear matrix inequalities. The strategy applies to general static output feedback design without any restrictions on the plant output equation. An iterative algorithm based on linear matrix inequalities is proposed in order to compute the feedback gain matrix that maximizes the estimate of the closed-loop region of attraction. Numerical examples are provided to illustrate the applicability of this new approach in examples borrowed from the literature.


title: Comparison of end-to-end neural network architectures and data augmentation methods for automatic infant motility assessment using wearable sensors

abstract: Infant motility assessment using intelligent wearables is a promising new approach for assessment of infant neurophysiological development, and where efficient signal analysis plays a central role. This study investigates the use of different end-to-end neural network architectures for processing infant motility data from wearable sensors. We focus on the performance and computational burden of alternative sensor encoder and time-series modelling modules and their combinations. In addition, we explore the benefits of data augmentation methods in ideal and non-ideal recording conditions. The experiments are conducted using a data-set of multi-sensor movement recordings from 7-month-old infants, as captured by a recently proposed smart jumpsuit for infant motility assessment. Our results indicate that the choice of the encoder module has a major impact on classifier performance. For sensor encoders, the best performance was obtained with parallel 2-dimensional convolutions for intra-sensor channel fusion with shared weights for all sensors. The results also indicate that a relatively compact feature representation is obtainable for within-sensor feature extraction without a drastic loss to classifier performance. Comparison of time-series models revealed that feed-forward dilated convolutions with residual and skip connections outperformed all RNN-based models in performance, training time, and training stability. The experiments also indicate that data augmentation improves model robustness in simulated packet loss or sensor dropout scenarios. In particular, signal- and sensor-dropout-based augmentation strategies provided considerable boosts to performance without negatively affecting the baseline performance. Overall the results provide tangible suggestions on how to optimize end-to-end neural network training for multi-channel movement sensor data.


title: Point sets and functions inducing tangles of set separations

abstract: Tangles, as introduced by Robertson and Seymour, were designed as an indirect way of capturing clusters in graphs and matroids. They have since been shown to capture clusters in much broader discrete structures too. But not all tangles are induced by a set of points, let alone a cluster. We characterise those that are: the tangles that are induced by a subset of or function on the set of data points whose connectivity structure they are meant to capture. We offer two such characterisations. The first is in terms of how many small sides of a tangle's separations it takes to cover the ground set. The second uses a new notion of duality for oriented set separations that becomes possible if these are no longer required to be separations of graph or matroids.


title: Bimolecular chemistry in the ultracold regime

abstract: Advances in atomic, molecular, and optical (AMO) physics techniques allowed the cooling of simple molecules down to the ultracold regime ($\lesssim$ 1 mK), and opened the opportunities to study chemical reactions with unprecedented levels of control. This review covers recent developments in studying bimolecular chemistry at ultralow temperatures. We begin with a brief overview of methods for producing, manipulating, and detecting ultracold molecules. We then survey experimental works that exploit the controllability of ultracold molecules to probe and modify their long-range interactions. Further combining the use of physical chemistry techniques, such as mass spectrometry and ion imaging, significantly improved the detection of ultracold reactions and enabled explorations of their dynamics in the short-range. We discuss a series of studies on the reaction KRb + KRb $\rightarrow$ K$_2$ + Rb$_2$ initiated below 1 $\mu$K, including the direct observation of a long-lived complex, the demonstration of product rotational state control via conserved nuclear spins, and a test of the statistical model using the complete quantum state distribution of the products.


title: $L^p-L^q$ estimates for the circular maximal operators on Heisenberg radial functions

abstract: $L^p$ boundedness of the circular maximal function $\mathcal M_{\mathbb{H}^1}$ on the Heisenberg group $\mathbb{H}^1$ has received considerable attentions. While the problem still remains open, $L^p$ boundedness of $\mathcal M_{\mathbb{H}^1}$ on Heisenberg radial functions was recently shown for $p>2$ by Beltran, Guo, Hickman, and Seeger [2]. In this paper we extend their result considering the local maximal operator $M_{\mathbb{H}^1}$ which is defined by taking supremum over $1<t<2$. We prove $L^p-L^q$ estimates for $M_{\mathbb{H}^1}$ on Heisenberg radial functions on the optimal range of $p,q$ modulo the borderline cases. Our argument also provides a simpler proof of the aforementioned result due to Beltran et al.


title: Looking for astrometric signals below 20 m/s: A Jupiter-mass planet signature in $\epsilon$ Eri

abstract: The USNO ground-based astrometric program URAT-Bright in combination with the Hipparcos mission epoch astrometry provides precise proper motions of a thousand bright stars in the southern hemisphere on a time basis of about 25 years. Small but statistically significant differences between these proper motions and Gaia EDR3 data can reveal long-period exoplanets similar to Jupiter in the nearest star systems. The presence of such a planet orbiting the magnetically active dwarf $\epsilon$ Eri is confirmed from both URAT--Hipparcos--EDR3 data and Hipparcos--EDR3 data with a corresponding projected velocity of $(+5,+8)$ and $(+6,+13)$ m s$^{-1}$, respectively. These signals are formally significant at a 0.989 and 1.0 confidence. We conclude that the newest astrometric results confirm the existence of a long-period exoplanet orbiting $\epsilon$ Eri, which was marginally detected from precision radial velocity measurements some 20 years ago.


title: CrowdSpeech and VoxDIY: Benchmark Datasets for Crowdsourced Audio Transcription

abstract: Domain-specific data is the crux of the successful transfer of machine learning systems from benchmarks to real life. In simple problems such as image classification, crowdsourcing has become one of the standard tools for cheap and time-efficient data collection: thanks in large part to advances in research on aggregation methods. However, the applicability of crowdsourcing to more complex tasks (e.g., speech recognition) remains limited due to the lack of principled aggregation methods for these modalities. The main obstacle towards designing aggregation methods for more advanced applications is the absence of training data, and in this work, we focus on bridging this gap in speech recognition. For this, we collect and release CrowdSpeech -- the first publicly available large-scale dataset of crowdsourced audio transcriptions. Evaluation of existing and novel aggregation methods on our data shows room for improvement, suggesting that our work may entail the design of better algorithms. At a higher level, we also contribute to the more general challenge of developing the methodology for reliable data collection via crowdsourcing. In that, we design a principled pipeline for constructing datasets of crowdsourced audio transcriptions in any novel domain. We show its applicability on an under-resourced language by constructing VoxDIY -- a counterpart of CrowdSpeech for the Russian language. We also release the code that allows a full replication of our data collection pipeline and share various insights on best practices of data collection via crowdsourcing.


title: Import options for chemical energy carriers from renewable sources to Germany

abstract: Import and export of fossil energy carriers are cornerstones of energy systems world-wide. If energy systems are to become climate neutral and sustainable, fossil carriers need to be substituted with carbon neutral alternatives or electrified if possible. We investigate synthetic chemical energy carriers, H2, CH4, MeOH, NH3 and Fischer-Tropsch fuels (FTF), produced using electricity from RES as fossil substitutes. [...] We model the sourcing of feedstock chemicals, synthesis and transport along nine different Energy Supply Chains to Germany (DE) and compare import options for seven locations around the world against each other and with domestically sourced alternatives on the basis of their respective cost per unit of H2 and energy delivered. We find that for each type of chemical energy carrier, there is an import option with lower costs compared to domestic production in DE. No single exporting country or energy carrier has a unique cost advantage, since for each energy carrier and country there are cost-competitive alternatives. This allows exporter and infrastructure decisions to be made based on other criteria than energy and cost. The lowest cost means for importing of energy and H2 are by H2 pipeline from Denmark, Spain and Western Asia and Northern Africa starting at 36 EUR/MWhLHV to 42 EUR/MWh-LHV or 1.0 EUR/kg-H2 to 1.3 EUR/kg-H2 (in 2050, assuming 5 % p.a. capital cost). For complex energy carriers derived from H2 like CH4, NH3, MeOH or FTF, imports from Argentina by ship to DE are lower cost than closer exporters in the European Union or Western Asia and Northern Africa. For meeting H2 demand, direct H2 imports are more attractive than indirect routes using CH4, MeOH or NH3 imports and subsequent decomposition to H2 because of high capital investment costs and energetic losses. We make our model and data available under open licenses for adaptation and reuse.


title: Model Checking C++ Programs

abstract: In the last three decades, memory safety issues in system programming languages such as C or C++ have been one of the significant sources of security vulnerabilities. However, there exist only a few attempts with limited success to cope with the complexity of C++ program verification. Here we describe and evaluate a novel verification approach based on bounded model checking (BMC) and satisfiability modulo theories (SMT) to verify C++ programs formally. Our verification approach analyzes bounded C++ programs by encoding into SMT various sophisticated features that the C++ programming language offers, such as templates, inheritance, polymorphism, exception handling, and the Standard C++ Libraries. We formalize these features within our formal verification framework using a decidable fragment of first-order logic and then show how state-of-the-art SMT solvers can efficiently handle that. We implemented our verification approach on top of ESBMC. We compare ESBMC to LLBMC and DIVINE, which are state-of-the-art verifiers to check C++ programs directly from the LLVM bitcode. Experimental results show that ESBMC can handle a wide range of C++ programs, presenting a higher number of correct verification results. At the same time, it reduces the verification time if compared to LLBMC and DIVINE tools. Additionally, ESBMC has been applied to a commercial C++ application in the telecommunication domain and successfully detected arithmetic overflow errors, potentially leading to security vulnerabilities.


title: Conditional wavefunction theory: a unified treatment of molecular structure and nonadiabatic dynamics

abstract: We demonstrate that a conditional wavefunction theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wavefunction formally recasts the full interacting wavefunction of a closed system as a set of lower dimensional (conditional) coupled `slices'. We formulate a variational wavefunction ansatz based on a set of conditional wavefunction slices, and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wavefunctions, and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser driven proton transfer, and Berry phase effects induced by a conical intersection. This work paves the road for the application of conditional wavefunction theory in equilibrium and out of equilibrium ab-initio molecular simulations of finite and extended systems.


title: Who Votes for Library Bonds? A Principal Component Exploration

abstract: Previous research has shown a relationship between voter characteristics and voter support for tax bonds. These findings, however, are difficult to interpret because of the high degree of collinearity across the measures. From 13 demographic measures of voters in a library bond election, seven independent principal components were extracted which accounted for 95 percent of the variance. Whereas the direct demographic measures showed inconsistent relationships with voting, the principal components of low SES, college experience, female and service job were related to affirmative voting, while high home value was related to negative voting.


title: Constraining particle acceleration in Sgr A* with simultaneous GRAVITY, Spitzer, NuSTAR and Chandra observations

abstract: We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $\nu L_\nu\propto \nu^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $\nu L_\nu \propto \nu^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $\gamma_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $\gamma_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $\gamma_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $\gamma_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare.


title: Spectral flatness and the volume of intersections of $p$-ellipsoids

abstract: Motivated by classical works of Schechtman and Schmuckenschl\"ager on intersections of $\ell_p$-balls and recent ones in information-based complexity relating random sections of ellipsoids and the quality of random information in approximation problems, we study the threshold behavior of the asymptotic volume of intersections of generalized $p$-ellipsoids. The non-critical behavior is determined under a spectral flatness (Wiener entropy) condition on the semi-axes. In order to understand the critical case at the threshold, we prove a central limit theorem for $q$-norms of points sampled uniformly at random from a $p$-ellipsoid, which is obtained under Noether's condition on the semi-axes.


title: Temporal Analysis of Worldwide War

abstract: Analysis of wars and conflicts between regions has been an important topic of interest throughout the history of humankind. In the latter part of the 20th century, in the aftermath of two World Wars and the shadow of nuclear, biological, and chemical holocaust, more was written on the subject than ever before. Wars have a negative impact on a country's economy, social order, infrastructure, and public health. In this paper, we study the wars fought in history and draw conclusions from that. We explore the participation of countries in wars and the nature of relationships between various countries during different timelines. A big part of today's wars is fought against terrorism. Therefore, this study also attempts to shed light on different countries' exposure to terrorist encounters and analyses the impact of wars on a country's economy in terms of change in GDP.


title: Unitary coupled-cluster based self-consistent polarization propagator theory: a quadratic unitary coupled-cluster singles and doubles scheme

abstract: The development of a quadratic unitary coupled-cluster singles and doubles (qUCCSD) based self-consistent polarization propagator method is reported. We present a simple strategy for truncating the commutator expansion of the UCC transformed Hamiltonian $\bar{H}$. The qUCCSD method for the electronic ground-state includes up to double commutators for the amplitude equations and up to cubic commutators for the energy expression. The qUCCSD excited-state eigenvalue equations include up to double commutators for the singles-singles block of $\bar{H}$, single commutators for the singles-doubles and doubles-singles blocks, and the bare Hamiltonian for the doubles-doubles block. Benchmark qUCCSD calculations of the ground-state properties and excitation energies for representative molecules demonstrate significant improvement of the accuracy and robustness over the previous UCC3 scheme derived using M{\o}ller-Plesset perturbation theory.


title: Interplay between singlet and triplet pairings in multi-band two-dimensional oxide superconductors

abstract: We theoretically study the superconducting properties of multi-band two-dimensional transition metal oxide superconductors by analyzing not only the role played by conventional singlet pairings, but also by the triplet order parameters, favored by the spin-orbit couplings present in these ma- terials. In particular, we focus on the two-dimensional electron gas at the (001) interface between LaAlO3 and SrTiO3 band insulators where the low electron densities and the sizeable spin-orbit couplings affect the superconducting features. Our theoretical study is based on an extended su- perconducting mean-field analysis of the typical multi-band tight-binding Hamiltonian, as well as on a parallel analysis of the effective electronic bands in the low-momentum limit, including static on-site and inter-site intra-band attractive potentials under applied magnetic fields. The presence of triplet pairings is able to strongly reduce the singlet order parameters which, as a result, are no longer a monotonic function of the charge density. The interplay between the singlet and the triplet pairings affects the dispersion of quasi-particle excitations in the Brillouin zone and also induces anisotropy in the superconducting behavior under the action of an in-plane and of an out- of-plane magnetic fields. Finally, non-trivial topological superconducting states become stable as a function of the charge density, as well as of the magnitude and of the orientation of the magnetic field. In addition to the chiral, time-reversal breaking, topological superconducting phase, favored by the linear Rashba couplings and by the on-site attractive potentials in the presence of an out- of-plane magnetic field, we find that a time-reversal invariant topological helical superconducting phase is promoted by not-linear spin-orbit couplings and by the inter-site attractive interactions in the absence of magnetic field.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.supr-con quant-ph


title: Network Design with Service Requirements: Scaling-up the Size of Solvable Problems

abstract: Network design, a cornerstone of mathematical optimization, is about defining the main characteristics of a network satisfying requirements on connectivity, capacity, and level-of-service. It finds applications in logistics and transportation, telecommunications, data sharing, energy distribution, and distributed computing. In multi-commodity network design, one is required to design a network minimizing the installation cost of its arcs and the operational cost to serve a set of point-to-point connections. The definition of this prototypical problem was recently enriched by additional constraints imposing that each origin-destination of a connection is served by a single path satisfying one or more level-of-service requirements, thus defining the Network Design with Service Requirements [Balakrishnan, Li, and Mirchandani. Operations Research, 2017]. These constraints are crucial, e.g., in telecommunications and computer networks, in order to ensure reliable and low-latency communication. In this paper we provide a new formulation for the problem, where variables are associated with paths satisfying the end-to-end service requirements. We present a fast algorithm for enumerating all the exponentially-many feasible paths and, when this is not viable, we provide a column generation scheme that is embedded into a branch-and-cut-and-price algorithm. Extensive computational experiments on a large set of instances show that our approach is able to move a step further in the solution of the Network Design with Service Requirements, compared with the current state-of-the-art.


title: Quantum scrambling of observable algebras

abstract: In this paper we describe an algebraic/geometrical approach to quantum scrambling. Generalized quantum subsystems are described by an hermitian-closed unital subalgebra $\cal A$ of operators evolving through a unitary channel. Qualitatively, quantum scrambling is defined by how the associated physical degrees of freedom get mixed up with others by the dynamics. Quantitatively, this is accomplished by introducing a measure, the geometric algebra anti-correlator (GAAC), of the self-orthogonalization of the commutant of $\cal A$ induced by the dynamics. This approach extends and unifies averaged bipartite OTOC, operator entanglement, coherence generating power and Loschmidt echo. Each of these concepts is indeed recovered by a special choice of $\cal A$. We compute typical values of GAAC for random unitaries, we prove upper bounds and characterize their saturation. For generic energy spectrum we find explicit expressions for the infinite-time average of the GAAC which encode the relation between $\cal A$ and the full system of Hamiltonian eigenstates. Finally, a notion of ${\cal A}$-chaoticity is suggested.


title: Generalized Multivariate Signs for Nonparametric Hypothesis Testing in High Dimensions

abstract: High-dimensional data, where the dimension of the feature space is much larger than sample size, arise in a number of statistical applications. In this context, we construct the generalized multivariate sign transformation, defined as a vector divided by its norm. For different choices of the norm function, the resulting transformed vector adapts to certain geometrical features of the data distribution. Building up on this idea, we obtain one-sample and two-sample testing procedures for mean vectors of high-dimensional data using these generalized sign vectors. These tests are based on U-statistics using kernel inner products, do not require prohibitive assumptions, and are amenable to a fast randomization-based implementation. Through experiments in a number of data settings, we show that tests using generalized signs display higher power than existing tests, while maintaining nominal type-I error rates. Finally, we provide example applications on the MNIST and Minnesota Twin Studies genomic data.


title: An Efficient Particle Tracking Algorithm for Large-Scale Parallel Pseudo-Spectral Simulations of Turbulence

abstract: Particle tracking in large-scale numerical simulations of turbulent flows presents one of the major bottlenecks in parallel performance and scaling efficiency. Here, we describe a particle tracking algorithm for large-scale parallel pseudo-spectral simulations of turbulence which scales well up to billions of tracer particles on modern high-performance computing architectures. We summarize the standard parallel methods used to solve the fluid equations in our hybrid MPI/OpenMP implementation. As the main focus, we describe the implementation of the particle tracking algorithm and document its computational performance. To address the extensive inter-process communication required by particle tracking, we introduce a task-based approach to overlap point-to-point communications with computations, thereby enabling improved resource utilization. We characterize the computational cost as a function of the number of particles tracked and compare it with the flow field computation, showing that the cost of particle tracking is very small for typical applications.


title: Memory Efficient Meta-Learning with Large Images

abstract: Meta learning approaches to few-shot classification are computationally efficient at test time, requiring just a few optimization steps or single forward pass to learn a new task, but they remain highly memory-intensive to train. This limitation arises because a task's entire support set, which can contain up to 1000 images, must be processed before an optimization step can be taken. Harnessing the performance gains offered by large images thus requires either parallelizing the meta-learner across multiple GPUs, which may not be available, or trade-offs between task and image size when memory constraints apply. We improve on both options by proposing LITE, a general and memory efficient episodic training scheme that enables meta-training on large tasks composed of large images on a single GPU. We achieve this by observing that the gradients for a task can be decomposed into a sum of gradients over the task's training images. This enables us to perform a forward pass on a task's entire training set but realize significant memory savings by back-propagating only a random subset of these images which we show is an unbiased approximation of the full gradient. We use LITE to train meta-learners and demonstrate new state-of-the-art accuracy on the real-world ORBIT benchmark and 3 of the 4 parts of the challenging VTAB+MD benchmark relative to leading meta-learners. LITE also enables meta-learners to be competitive with transfer learning approaches but at a fraction of the test-time computational cost, thus serving as a counterpoint to the recent narrative that transfer learning is all you need for few-shot classification.


title: Screening for a Reweighted Penalized Conditional Gradient Method

abstract: The conditional gradient method (CGM) is widely used in large-scale sparse convex optimization, having a low per iteration computational cost for structured sparse regularizers and a greedy approach to collecting nonzeros. We explore the sparsity acquiring properties of a general penalized CGM (P-CGM) for convex regularizers and a reweighted penalized CGM (RP-CGM) for nonconvex regularizers, replacing the usual convex constraints with gauge-inspired penalties. This generalization does not increase the per-iteration complexity noticeably. Without assuming bounded iterates or using line search, we show $O(1/t)$ convergence of the gap of each subproblem, which measures distance to a stationary point. We couple this with a screening rule which is safe in the convex case, converging to the true support at a rate $O(1/(\delta^2))$ where $\delta \geq 0$ measures how close the problem is to degeneracy. In the nonconvex case the screening rule converges to the true support in a finite number of iterations, but is not necessarily safe in the intermediate iterates. In our experiments, we verify the consistency of the method and adjust the aggressiveness of the screening rule by tuning the concavity of the regularizer.


title: Modified Kerr-Schild Method applied to Gravity f(R)

abstract: We propose an adaptation of the Kerr-Schild method by implementing the correspondence relations (mapping) between Ricci-based Gravity (RBG) and General Relativity (GR). Basically, we generate GR known solutions from a canonical metric with a well-defined form and then obtain the configuration of this solution for the Gravity f(R). The new method allows us to perturb the metric associated with GR with a null geodetic vector, but instead of taking us to a new solution described by Einstein's gravitational sector, it allows us to configure this solution in an RBG.


title: Lower bounds on mapping content and quantitative factorization through trees

abstract: We give a simple quantitative condition, involving the "mapping content" of Azzam--Schul, that implies that a Lipschitz map from a Euclidean space to a metric space must be close to factoring through a tree. Using results of Azzam--Schul and the present authors, this gives simple checkable conditions for a Lipschitz map to have a large piece of its domain on which it behaves like an orthogonal projection. The proof involves new lower bounds and continuity statements for mapping content, and relies on a "qualitative" version of the main theorem recently proven by Esmayli--Haj{\l}asz.


title: On spectral properties of the Schreier graphs of the Thompson group $F$

abstract: In this article we study spectral properties of the family of Schreier graphs associated to the action of the Thompson group $F$ on the interval [0,1]. In particular, we describe spectra of Laplace type operators associated to these Schreier graphs and calculate certain spectral measures associated to the Schreier graph $\Upsilon$ of the orbit of 1/2. As a byproduct we calculate the asymptotics of the return probabilities of the simple random walk on $\Upsilon$ starting at 1/2. In addition, given a Laplace type operator $L$ on a tree-like graph we study relations between the spectral measures of $L$ associated to delta functions of different vertices and the spectrum of $L$.


title: Decision-Making Technology for Autonomous Vehicles Learning-Based Methods, Applications and Future Outlook

abstract: Autonomous vehicles have a great potential in the application of both civil and military fields, and have become the focus of research with the rapid development of science and economy. This article proposes a brief review on learning-based decision-making technology for autonomous vehicles since it is significant for safer and efficient performance of autonomous vehicles. Firstly, the basic outline of decision-making technology is provided. Secondly, related works about learning-based decision-making methods for autonomous vehicles are mainly reviewed with the comparison to classical decision-making methods. In addition, applications of decision-making methods in existing autonomous vehicles are summarized. Finally, promising research topics in the future study of decision-making technology for autonomous vehicles are prospected.


title: The Interplay between Memory and Potentials of Mean Force: A Discussion on the Structure of Equations of Motion for Coarse Grained Observables

abstract: The underdamped, non-linear, generalized Langevin equation is widely used to model coarse-grained dynamics of soft and biological materials. By means of a projection operator formalism, we show under which approximations this equation can be obtained from the Hamiltonian dynamics of the underlying microscopic system and in which cases it makes sense to introduce a potential of mean force. We discuss shortcomings of previous derivations presented in the literature and demonstrate the implications of our derivation for the structure of memory terms and their connection to generalized fluctuation-dissipation relations. We show, in particular, that the widely used, simple structure which contains a potential of mean force, a memory term which is linear in the observable, and a fluctuating force which is related to the memory term by a fluctuation-dissipation relation, is neither exact nor can it, in general, be derived as a controlled approximation to the exact dynamics.


title: Evidence for cold plasma in planetary nebulae from radio observations with the LOw Frequency ARray (LOFAR)

abstract: We present observations of planetary nebulae with the LOw Frequency ARray (LOFAR) between 120 and 168 MHz. The images show thermal free-free emission from the nebular shells. We have determined the electron temperatures for spatially resolved, optically thick nebulae. These temperatures are 20 to 60% lower than those estimated from collisionally excited optical emission lines. This strongly supports the existence of a cold plasma component, which co-exists with hot plasma in planetary nebulae. This cold plasma does not contribute to the collisionally excited lines, but does contribute to recombination lines and radio flux. Neither of the plasma components are spatially resolved in our images, although we infer that the cold plasma extends to the outer radii of planetary nebulae. However, more cold plasma appears to exist at smaller radii. The presence of cold plasma should be taken into account in modeling of radio emission of planetary nebulae. Modelling of radio emission usually uses electron temperatures calculated from collisionally excited optical and/or infrared lines. This may lead to an underestimate of the ionized mass and an overestimate of the extinction correction from planetary nebulae when derived from the radio flux alone. The correction improves the consistency of extinction derived from the radio fluxes when compared to estimates from the Balmer decrement flux ratios.


title: Measuring Information Leakage in Non-stochastic Brute-Force Guessing

abstract: This paper proposes an operational measure of non-stochastic information leakage to formalize privacy against a brute-force guessing adversary. The information is measured by non-probabilistic uncertainty of uncertain variables, the non-stochastic counterparts of random variables. For $X$ that is related to released data $Y$, the non-stochastic brute-force leakage is measured by the complexity of exhaustively checking all the possibilities of the private attribute $U$ of $X$ by an adversary. The complexity refers to the number of trials to successfully guess $U$. Maximizing this leakage over all possible private attributes $U$ gives rise to the maximal (i.e., worst-case) non-stochastic brute-force guessing leakage. This is proved to be fully determined by the minimal non-stochastic uncertainty of $X$ given $Y$, which also determines the worst-case attribute $U$ indicating the highest privacy risk if $Y$ is disclosed. The maximal non-stochastic brute-force guessing leakage is shown to be proportional to the non-stochastic identifiability of $X$ given $Y$ and upper bounds the existing maximin information. The latter quantifies the information leakage when an adversary must perfectly guess $U$ in one-shot via $Y$. Experiments are used to demonstrate the tradeoff between the maximal non-stochastic brute-force guessing leakage and the data utility (measured by the maximum quantization error) and to illustrate the relationship between maximin information and stochastic one-shot maximal leakage.


title: The single copy of the gravitational holonomy

abstract: The double copy is a well-established relationship between gravity and gauge theories. It relates perturbative scattering amplitudes as well as classical solutions, and recently there has been mounting evidence that it also applies to non-perturbative information. In this paper, we consider the holonomy properties of manifolds in gravity and prescribe a single copy of gravitational holonomy that differs from the holonomy in gauge theory. We discuss specific cases and give examples where the single copy holonomy group is reduced. Our results may prove useful in extending the classical double copy. We also clarify previous misconceptions in the literature regarding gravitational Wilson lines and holonomy.


title: Black Hole Shadows Constrain Extended Gravity

abstract: The first images of black hole shadows open new possibilities to constrain modern extended gravity theories. We present the method of shadow background calculation for black hole solutions in the form of Taylor series where $g_{11} = - g_{00}^{-1}$. The method is extended to general non-rotating case $g_{11} \neq - g_{00}^{-1}$. The results of the analysis are compared with the predictions of General relativity taking into account the Event Horizon Telescope data. The results for the Horndesky model with the Gauss-Bonnet invariant, loop quantum gravity, Bumblebee model and Gauss-Bonnet gravity are in full agreement with the observations of M87*. In conformal gravity, large values of $m_2$ and $Q_s$ must be excluded. In STEGR $f(Q)$ gravity the observational limits on the parameter $\alpha$ are: $-0.025<\alpha<0.04$. For an alternative generalization of the Bumblebee model with the Schwarzschild approximation: $-0.3 < l < 0.45$. These results demonstrate the maximum one can achieve without taking into account of the rotation of a black hole.


title: Optimization of a quantum control sequence for initializing an NV spin register

abstract: Implementation of many quantum information protocols require an efficient initialization of the quantum register. In the present report, we optimize a population trapping protocol for initializing a hybrid spin register associated a single nitrogen vacancy (NV) center in diamond. We initialize the quantum register by polarizing the electronic and the nuclear spins of the NV with a sequence of microwave, radio-frequency and optical pulses. We use a rate equation model to explain the distribution of population under the effect of the optical pulses. The model is compared to the experimental data obtained by performing partial quantum state tomography. To further increase the spin polarization, we propose a recursive protocol with optimized optical pulses.


title: A Weighted and Normalized Gould-Fernandez brokerage measure

abstract: The Gould and Fernandez local brokerage measure defines brokering roles based on the group membership of the nodes from the incoming and outgoing edges. This paper extends on this brokerage measure to account for weighted edges and introduces the Weighted-Normalized Gould-Fernandez measure (WNGF). The value added of this new measure is demonstrated empirically with both a macro level trade network and a micro level organization network. The measure is first applied to the EUREGIO inter-regional trade dataset and then to an organizational network in a research and development group. The results gained from the WNGF measure are compared to those from two dichotomized networks: a threshold and a multiscale backbone network. The results show that the WNGF generates valid results, consistent with those of the dichotomized network. In addition, it provides the following advantages: (i) it ensures information retention, (ii) since no alterations and decisions have to be made on how to dichotomize the network, the WNGF frees the user from the burden of making assumptions, (iii) it provides a nuanced understanding of each node's brokerage role. These advantages are of special importance when the role of less connected nodes is considered. The two empirical networks used here are for illustrative purposes. Possible applications of WNGF span beyond regional and organizational studies, and into all those contexts where retaining weights is important, for example by accounting for persisting or repeating edges compared to one-time interactions. WNGF can also be used to further analyze networks that measure how often people meet, talk, text, like, or retweet. WNGF makes a relevant methodological contribution as it offers a way to analyze brokerage in weighted, directed, and even complete graphs without information loss that can be used across disciplines and different type of networks.


title: Development of a Josephson junction based single photon microwave detector for axion detection experiments

abstract: Josephson junctions, in appropriate configurations, can be excellent candidates for detection of single photons in the microwave frequency band. Such possibility has been recently addressed in the framework of galactic axion detection. Here are reported recent developments in the modelling and simulation of dynamic behaviour of a Josephson junction single microwave photon detector. For a Josephson junction to be enough sensitive, small critical currents and operating temperatures of the order of ten of mK are necessary. Thermal and quantum tunnelling out of the zero-voltage state can also mask the detection process. Axion detection would require dark count rates in the order of 0.001 Hz. It is, therefore, is of paramount importance to identify proper device fabrication parameters and junction operation point.


title: A Micro-Service based Approach for Constructing Distributed Storage System

abstract: This paper presents an approach for constructing distributed storage system based on micro-service architecture. By building storage functionalities using micro services, we can provide new capabilities to a distributed storage system in a flexible way. We take erasure coding and compression as two case studies to show how to build a micro-service based distributed storage system. We also show that by building erasure coding and compression as micro-services, the distributed storage system still achieves reasonable performance compared to the monolithic one.


title: Asymptotic Analysis of Statistical Estimators related to MultiGraphex Processes under Misspecification

abstract: This article studies the asymptotic properties of Bayesian or frequentist estimators of a vector of parameters related to structural properties of sequences of graphs. The estimators studied originate from a particular class of graphex model introduced by Caron and Fox. The analysis is however performed here under very weak assumptions on the underlying data generating process, which may be different from the model of Caron and Fox or from a graphex model. In particular, we consider generic sparse graph models, with unbounded degree, whose degree distribution satisfies some assumptions. We show that one can relate the limit of the estimator of one of the parameters to the sparsity constant of the true graph generating process. When taking a Bayesian approach, we also show that the posterior distribution is asymptotically normal. We discuss situations where classical random graphs models such as configuration models, sparse graphon models, edge exchangeable models or graphon processes satisfy our assumptions.


title: Decarbonising the EU Power Sector: a Technological and Socio-economic Analysis and the Role of Nuclear

abstract: Low-carbon electricity is a key enabler in combating climate change. Decarbonising the power sector is now at the centre of global and European policies. As the IPCC highlights, pathways where the power sector rapidly decarbonises by 2030 have higher chances of keeping global warming below 1.5$^\circ$C. The electricity sector should be fully decarbonised by 2050 to meet either the 1.5$^\circ$C or 2$^\circ$C targets. This means that EU policy efforts should focus on supporting a maximum reduction of emissions per unit of electricity by 2030 and net-zero emissions by 2050. Reaching these targets is one of the most pressing questions EU policymakers face today. In light of the COVID-19 crisis, EU policies should guide a cost-effective, reliable and environmentally sound transition of the power sector, benefiting EU research and innovation and its citizens. This meta-analysis provides a novel view on historical data and compares data from modelling scenarios identified in the literature. It assesses the current and future role of nuclear energy in decarbonizing the EU power sector, while reviewing socio-economic implications that could arise if limited public support nearly excludes nuclear fission electricity from the future EU power mix. This work highlights relevant socio-economic policy implications and actionable policy recommendations.


title: Review of synergic meteor observations: linking the results from cameras, ionosondes, infrasound and seismic detectors

abstract: Joint evaluation of different meteor observation types support the better understanding of both the meteor phenomenon and the terrestrial atmosphere. Two types of examples are presented in this work, linking ionospheric effects to specific meteors, where almost one third of the meteors emerged at high altitudes were simultaneously recorded with an optical camera. Very few such observations have been realized yet. With a daytime fireball, the recorded infrasound effect and the atmospheric blast produced shock wave related small earthquakes were identified by a network of ground stations. The overview of these observational types highlights specific topics where substantial improvements and discoveries are expected in the near future.


title: Constraints on anomalous quartic gauge couplings via $Z\gamma jj$ production at the LHC

abstract: In this paper, we investigate the contributions of anomalous quartic gauge couplings (aQGCs) to $Z\gamma jj$ production at the Large Hadron Collider (LHC) in the context of Standard Model effective theory (SMEFT). When energy scale is large, the validity of SMEFT becomes an important issue. To ensure the validity, the unitarity bound is applied in a model independent approach, which is found to have significant suppressive effects on the signals of $O_{M_i}$ operators. The kinematic and polarization features of the aQGC signals are also studied. The polarization effect is useful to highlight the signals of $O_{T_i}$ operators. The sensitivity estimates on dimension-8 operators with unitarity bounds at $\sqrt{s}=14$ TeV are obtained.


title: Global Structure Identifiability and Reconstructibility of an NDS with Descriptor Subsystems

abstract: This paper investigates requirements on a networked dynamic system (NDS) such that its subsystem interactions can be solely determined from experiment data or reconstructed from its overall model. The NDS is constituted from several subsystems whose dynamics are described through a descriptor form. Except regularity on each subsystem and the whole NDS, no other restrictions are put on either subsystem dynamics or subsystem interactions. A matrix rank based necessary and sufficient condition is derived for the global identifiability of subsystem interactions, which leads to several conclusions about NDS structure identifiability when there is some a priori information. This matrix also gives an explicit description for the set of subsystem interactions that can not be distinguished from experiment data only. In addition, under a well-posedness assumption, a necessary and sufficient condition is obtained for the reconstructibility of subsystem interactions from an NDS descriptor form model. This condition can be verified with each subsystem separately and is therefore attractive in the analysis and synthesis of a large-scale NDS. Simulation results show that rather than increases monotonically with the distance of subsystem interactions to the undifferentiable set, the magnitude of the external output differences between two NDSs with distinct subsystem interactions increases much more rapidly when one of them is close to be unstable. In addition, directions of probing signals are also very important in distinguishing external outputs of distinctive NDSs.These findings are expected to be helpful in identification experiment designs, etc.


title: On Measuring and Controlling the Spectral Bias of the Deep Image Prior

abstract: The deep image prior showed that a randomly initialized network with a suitable architecture can be trained to solve inverse imaging problems by simply optimizing it's parameters to reconstruct a single degraded image. However, it suffers from two practical limitations. First, it remains unclear how to control the prior beyond the choice of the network architecture. Second, training requires an oracle stopping criterion as during the optimization the performance degrades after reaching an optimum value. To address these challenges we introduce a frequency-band correspondence measure to characterize the spectral bias of the deep image prior, where low-frequency image signals are learned faster and better than high-frequency counterparts. Based on our observations, we propose techniques to prevent the eventual performance degradation and accelerate convergence. We introduce a Lipschitz-controlled convolution layer and a Gaussian-controlled upsampling layer as plug-in replacements for layers used in the deep architectures. The experiments show that with these changes the performance does not degrade during optimization, relieving us from the need for an oracle stopping criterion. We further outline a stopping criterion to avoid superfluous computation. Finally, we show that our approach obtains favorable results compared to current approaches across various denoising, deblocking, inpainting, super-resolution and detail enhancement tasks. Code is available at \url{https://github.com/shizenglin/Measure-and-Control-Spectral-Bias}.


title: Interactive Causal Structure Discovery in Earth System Sciences

abstract: Causal structure discovery (CSD) models are making inroads into several domains, including Earth system sciences. Their widespread adaptation is however hampered by the fact that the resulting models often do not take into account the domain knowledge of the experts and that it is often necessary to modify the resulting models iteratively. We present a workflow that is required to take this knowledge into account and to apply CSD algorithms in Earth system sciences. At the same time, we describe open research questions that still need to be addressed. We present a way to interactively modify the outputs of the CSD algorithms and argue that the user interaction can be modelled as a greedy finding of the local maximum-a-posteriori solution of the likelihood function, which is composed of the likelihood of the causal model and the prior distribution representing the knowledge of the expert user. We use a real-world data set for examples constructed in collaboration with our co-authors, who are the domain area experts. We show that finding maximally usable causal models in the Earth system sciences or other similar domains is a difficult task which contains many interesting open research questions. We argue that taking the domain knowledge into account has a substantial effect on the final causal models discovered.


title: A Deep Learning Based Discontinuous Galerkin Method for Hyperbolic Equations with Discontinuous Solutions and Random Uncertainties

abstract: We propose a deep learning based discontinuous Galerkin method (D2GM) to solve hyperbolic equations with discontinuous solutions and random uncertainties. The main computational challenges for such problems include discontinuities of the solutions and the curse of dimensionality due to uncertainties. Deep learning techniques have been favored for high-dimensional problems but face difficulties when the solution is not smooth, thus have so far been mainly used for viscous hyperbolic system that admits only smooth solutions. We alleviate this difficulty by setting up the loss function using discrete shock capturing schemes--the discontinous Galerkin method as an example--since the solutions are smooth in the discrete space. The convergence of D2GM is established via the Lax equivalence theorem kind of argument. The high-dimensional random space is handled by the Monte-Carlo method. Such a setup makes the D2GM approximate high-dimensional functions over the random space with satisfactory accuracy at reasonable cost. The D2GM is found numerically to be first-order and second-order accurate for (stochastic) linear conservation law with smooth solutions using piecewise constant and piecewise linear basis functions, respectively. Numerous examples are given to verify the efficiency and the robustness of D2GM with the dimensionality of random variables up to $200$ for (stochastic) linear conservation law and (stochastic) Burgers' equation.


title: Defect-Dependent Corrugation in Graphene

abstract: Graphene's intrinsically corrugated and wrinkled topology fundamentally influences its electronic, mechanical, and chemical properties. Experimental techniques allow the manipulation of pristine graphene and the controlled production of defects which allows to control the atomic out-of-plane fluctuations and, thus, tune graphene's properties. Here, we perform large scale machine learning-driven molecular dynamics simulations to understand the impact of defects on the structure of graphene. We find that defects cause significantly higher corrugation leading to a strongly wrinkled surface. The magnitude of this structural transformation strongly depends on the defect concentration and specific type of defect. Analysing the atomic neighborhood of the defects reveals that the extent of these morphological changes depends on the preferred geometrical orientation and the interactions between defects. While our work highlights that defects can strongly affect graphene's morphology, it also emphasises the differences between distinct types by linking the global structure to the local environment of the defects.


title: Distribution of the total angular momentum in relativistic configurations

abstract: This paper is devoted to the analysis of the distribution of the total angular momentum in a relativistic configuration. Using cumulants and generating function formalism this analysis can be reduced to the study of individual subshells with $N$ equivalent electrons of momentum $j$. An expression as a nth-derivative is provided for the generating function of the $J$ distribution and efficient recurrence relations are established. It is shown that this distribution may be represented by a Gram-Charlier-like series which is derived from the corresponding series for the magnetic quantum number distribution. The numerical efficiency of this expansion is fair when the configuration consists of several subshells, while the accuracy is less good when only one subshell is involved. An analytical expression is given for the odd-order momenta while the even-order ones are expressed as a series which provides an acceptable accuracy though being not convergent. Such expressions may be used to obtain approximate values for the number of transitions in a spin-orbit split array: it is shown that the approximation is often efficient when few terms are kept, while some complex cases require to include a large number of terms.