title: On members of Lucas sequences which are products of Catalan numbers

abstract: We show that if $\{U_n\}_{n\geq 0}$ is a Lucas sequence, then the largest $n$ such that $|U_n|=C_{m_1}C_{m_2}\cdots C_{m_k}$ with $1\leq m_1\leq m_2\leq \cdots\leq m_k$, where $C_m$ is the $m$th Catalan number satisfies $n<6500$. In case the roots of the Lucas sequence are real, we have $n\in \{1,2, 3, 4, 6, 8, 12\}$. As a consequence, we show that if $\{X_n\}_{n\geq 1}$ is the sequence of the $X$ coordinates of a Pell equation $X^2-dY^2=\pm 1$ with a nonsquare integer $d>1$, then $X_n=C_m$ implies $n=1$.


title: A Randomized Algorithm to Reduce the Support of Discrete Measures

abstract: Given a discrete probability measure supported on $N$ atoms and a set of $n$ real-valued functions, there exists a probability measure that is supported on a subset of $n+1$ of the original $N$ atoms and has the same mean when integrated against each of the $n$ functions. If $ N \gg n$ this results in a huge reduction of complexity. We give a simple geometric characterization of barycenters via negative cones and derive a randomized algorithm that computes this new measure by "greedy geometric sampling". We then study its properties, and benchmark it on synthetic and real-world data to show that it can be very beneficial in the $N\gg n$ regime. A Python implementation is available at \url{https://github.com/FraCose/Recombination_Random_Algos}.


title: The Ritva Blockchain: Enabling Confidential Transactions at Scale

abstract: The distributed ledger technology has been widely hailed as the break-through technology. It has realised a great number of application scenarios, and improved workflow of many domains. Nonetheless, there remain a few major concerns in adopting and deploying the distributed ledger technology at scale. In this white paper, we tackle two of them, namely the throughput scalability and confidentiality protection for transactions. We learn from the existing body of research, and build a scale-out blockchain platform that champions privacy called RVChain. RVChain takes advantage of trusted execution environment to offer confidentiality protection for transactions, and scales the throughput of the network in proportion with the number of network participants by supporting parallel shadow chains.


title: Sparse Perturbations for Improved Convergence in Stochastic Zeroth-Order Optimization

abstract: Interest in stochastic zeroth-order (SZO) methods has recently been revived in black-box optimization scenarios such as adversarial black-box attacks to deep neural networks. SZO methods only require the ability to evaluate the objective function at random input points, however, their weakness is the dependency of their convergence speed on the dimensionality of the function to be evaluated. We present a sparse SZO optimization method that reduces this factor to the expected dimensionality of the random perturbation during learning. We give a proof that justifies this reduction for sparse SZO optimization for non-convex functions without making any assumptions on sparsity of objective function or gradient. Furthermore, we present experimental results for neural networks on MNIST and CIFAR that show faster convergence in training loss and test accuracy, and a smaller distance of the gradient approximation to the true gradient in sparse SZO compared to dense SZO.


title: Modelling of daily reference evapotranspiration using deep neural network in different climates

abstract: Precise and reliable estimation of reference evapotranspiration (ET o ) is an essential for the irrigation and water resources management. ET o is difficult to predict due to its complex processes. This complexity can be solved using machine learning methods. This study investigates the performance of artificial neural network (ANN) and deep neural network (DNN) models for estimating daily ET o . Previously proposed ANN and DNN methods have been realized, and their performances have been compared. Six input data including maximum air temperature (T max ), minimum air temperature (T min ), solar radiation (R n ), maximum relative humidity (RH max ), minimum relative humidity (RH min ) and wind speed (U 2 ) are used from 4 meteorological stations (Adana, Aksaray, Isparta and Ni\u{g}de) during 1999-2018 in Turkey. The results have shown that our proposed DNN models achieves satisfactory accuracy for daily ET o estimation compared to previous ANN and DNN models. The best performance has been observed with the proposed model of DNN with SeLU activation function (P-DNN-SeLU) in Aksaray with coefficient of determination (R 2 ) of 0.9934, root mean square error (RMSE) of 0.2073 and mean absolute error (MAE) of 0.1590, respectively. Therefore, the P-DNN-SeLU model could be recommended for estimation of ET o in other climate zones of the world.


title: The isotropy group of a foliation: the local case

abstract: Given a holomorphic singular foliation $\fa$ of $(\C^n,0)$ we define $Iso(\fa)$ as the group of germs of biholomorphisms on $(\C^n,0)$ preserving $\fa$: $Iso(\fa)=\{\Phi\in Diff(\C^n,0)\,|\,\Phi^*(\fa)=\fa\}$. The normal subgroup of $Iso(\fa)$, of biholomorphisms sending each leaf of $\fa$ into itself, will be denoted as $Fix(\fa)$. The corresponding groups of formal biholomorphisms will be denoted as $\wh{Iso}(\fa)$ and $\wh{Fix}(\fa)$, respectively. The purpose of this paper will be to study the quotients $Iso(\fa)/Fix(\fa)$ and $\wh{Fix}(\fa)/\wh{Fix}(\fa)$, mainly in the case of codimension one foliation.


title: Cosmological constraints on non-adiabatic dark energy perturbations

abstract: The exact nature of dark energy is currently unknown and its cosmological perturbations, when dark energy is assumed not to be the cosmological constant, are usually modeled as adiabatic. Here we explore the possibility that dark energy might have a nonadiabatic component and we examine how it would affect several key cosmological observables. We present analytical solutions for the growth rate and growth index of matter density perturbations and compare them to both numerical solutions of the fluid equations and an implementation in the Boltzmann code CLASS, finding that they all agree to well below one percent. We also perform a Monte Carlo analysis to derive constraints on the parameters of the nonadiabatic component using the latest cosmological data, including the temperature and polarization spectra of the Cosmic Microwave Background as observed by Planck, the Baryon Acoustic Oscillations, the Pantheon type Ia supernovae compilation and lastly, measurements of Redshift Space Distortions (RSD) of the growth rate of matter perturbations. We find that the amplitude of the nonadiabatic pressure perturbation is consistent with zero within $1\sigma$. Finally, we also present a new, publicly available, RSD likelihood for MontePython based on the "Gold 2018" growth rate data compilation.


title: A Tsallisian Universe

abstract: In this paper, we adopt the Verlinde hypothesis on the origin of gravity as the consequence of the tendency of systems to increase their entropy, and employ the Tsallis statistics. Thereinafter, modifications to the Newtonian second law of motion, its gravity and radial velocity profile are studied. In addition and in a classical framework, the corresponding cosmology and also its ability in describing the inflationary phases are investigated.


title: Networks beyond pairwise interactions: structure and dynamics

abstract: The complexity of many biological, social and technological systems stems from the richness of the interactions among their units. Over the past decades, a great variety of complex systems has been successfully described as networks whose interacting pairs of nodes are connected by links. Yet, in face-to-face human communication, chemical reactions and ecological systems, interactions can occur in groups of three or more nodes and cannot be simply described just in terms of simple dyads. Until recently, little attention has been devoted to the higher-order architecture of real complex systems. However, a mounting body of evidence is showing that taking the higher-order structure of these systems into account can greatly enhance our modeling capacities and help us to understand and predict their emerging dynamical behaviors. Here, we present a complete overview of the emerging field of networks beyond pairwise interactions. We first discuss the methods to represent higher-order interactions and give a unified presentation of the different frameworks used to describe higher-order systems, highlighting the links between the existing concepts and representations. We review the measures designed to characterize the structure of these systems and the models proposed in the literature to generate synthetic structures, such as random and growing simplicial complexes, bipartite graphs and hypergraphs. We introduce and discuss the rapidly growing research on higher-order dynamical systems and on dynamical topology. We focus on novel emergent phenomena characterizing landmark dynamical processes, such as diffusion, spreading, synchronization and games, when extended beyond pairwise interactions. We elucidate the relations between higher-order topology and dynamical properties, and conclude with a summary of empirical applications, providing an outlook on current modeling and conceptual frontiers.

categories: physics.soc-ph cond-mat.dis-nn cs.SI nlin.AO q-bio.NC


title: AnalogNet: Convolutional Neural Network Inference on Analog Focal Plane Sensor Processors

abstract: We present a high-speed, energy-efficient Convolutional Neural Network (CNN) architecture utilising the capabilities of a unique class of devices known as analog Focal Plane Sensor Processors (FPSP), in which the sensor and the processor are embedded together on the same silicon chip. Unlike traditional vision systems, where the sensor array sends collected data to a separate processor for processing, FPSPs allow data to be processed on the imaging device itself. This unique architecture enables ultra-fast image processing and high energy efficiency, at the expense of limited processing resources and approximate computations. In this work, we show how to convert standard CNNs to FPSP code, and demonstrate a method of training networks to increase their robustness to analog computation errors. Our proposed architecture, coined AnalogNet, reaches a testing accuracy of 96.9% on the MNIST handwritten digits recognition task, at a speed of 2260 FPS, for a cost of 0.7 mJ per frame.


title: Computing spectral measures of self-adjoint operators

abstract: Using the resolvent operator, we develop an algorithm for computing smoothed approximations of spectral measures associated with self-adjoint operators. The algorithm can achieve arbitrarily high-orders of convergence in terms of a smoothing parameter for computing spectral measures of general differential, integral, and lattice operators. Explicit pointwise and $L^p$-error bounds are derived in terms of the local regularity of the measure. We provide numerical examples, including a partial differential operator, a magnetic tight-binding model of graphene, and compute one thousand eigenvalues of a Dirac operator to near machine precision without spectral pollution. The algorithm is publicly available in $\texttt{SpecSolve}$, which is a software package written in MATLAB.

categories: math.NA cs.NA math-ph math.FA math.MP math.SP


title: Artificial Neural Network Based Breast Cancer Screening: A Comprehensive Review

abstract: Breast cancer is a common fatal disease for women. Early diagnosis and detection is necessary in order to improve the prognosis of breast cancer affected people. For predicting breast cancer, several automated systems are already developed using different medical imaging modalities. This paper provides a systematic review of the literature on artificial neural network (ANN) based models for the diagnosis of breast cancer via mammography. The advantages and limitations of different ANN models including spiking neural network (SNN), deep belief network (DBN), convolutional neural network (CNN), multilayer neural network (MLNN), stacked autoencoders (SAE), and stacked de-noising autoencoders (SDAE) are described in this review. The review also shows that the studies related to breast cancer detection applied different deep learning models to a number of publicly available datasets. For comparing the performance of the models, different metrics such as accuracy, precision, recall, etc. were used in the existing studies. It is found that the best performance was achieved by residual neural network (ResNet)-50 and ResNet-101 models of CNN algorithm.


title: Dynamic scaling in the quenched disordered classical $N$-vector model

abstract: We revisit the effects of short-ranged random quenched disorder on the universal scaling properties of the classical $N$-vector model with cubic anisotropy. We set up the nonconserved relaxational dynamics of the model, and study the universal dynamic scaling near the second order phase transition. We extract the critical exponents and the dynamic exponent in a one-loop dynamic renormalisation group calculation with short-ranged isotropic disorder. We show that the dynamics near a critical point is generically slower when the quenched disorder is relevant than when it is not, independent of whether the pure model is isotropic or cubic anisotropic. We demonstrate the surprising thresholdless instability of the associated universality class due to perturbations from rotational invariance breaking quenched disorder-order parameter coupling, indicating breakdown of dynamic scaling. We speculate that this may imply a novel first order transition in the model, induced by a symmetry-breaking disorder.


title: Establishing a nearly closed cycling transition in a polyatomic molecule

abstract: We study optical cycling in the polar free radical calcium monohydroxide (CaOH) and establish an experimental path towards scattering $\sim$$10^4$ photons. We report rovibronic branching ratio measurements with precision at the $\sim10^{-4}$ level and observe weak symmetry-forbidden decays to bending modes with non-zero vibrational angular momentum. Calculations are in excellent agreement with these measurements and predict additional decay pathways. Additionally, we perform high-resolution spectroscopy of the $\widetilde{\text{X}}\,^2\Sigma^+(12^00)$ and $\widetilde{\text{X}}\,^2\Sigma^+(12^20)$ hybrid vibrational states of CaOH. These advances establish a path towards radiative slowing, 3D magneto-optical trapping, and sub-Doppler cooling of CaOH.


title: What's Sex Got To Do With Fair Machine Learning?

abstract: Debate about fairness in machine learning has largely centered around competing definitions of what fairness or nondiscrimination between groups requires. However, little attention has been paid to what precisely a group is. Many recent approaches to "fairness" require one to specify a causal model of the data generating process. These exercises make an implicit ontological assumption that a racial or sex group is simply a collection of individuals who share a given trait. We show this by exploring the formal assumption of modularity in causal models, which holds that the dependencies captured by one causal pathway are invariant to interventions on any other pathways. Causal models of sex propose two substantive claims: 1) There exists a feature, sex-on-its-own, that is an inherent trait of an individual that causally brings about social phenomena external to it in the world; and 2) the relations between sex and its effects can be modified in whichever ways and the former feature would still retain the meaning that sex has in our world. We argue that this ontological picture is false. Many of the "effects" that sex purportedly "causes" are in fact constitutive features of sex as a social status. They give the social meaning of sex features, meanings that are precisely what make sex discrimination a distinctively morally problematic type of action. Correcting this conceptual error has a number of implications for how models can be used to detect discrimination. Formal diagrams of constitutive relations present an entirely different path toward reasoning about discrimination. Whereas causal diagrams guide the construction of sophisticated modular counterfactuals, constitutive diagrams identify a different kind of counterfactual as central to an inquiry on discrimination: one that asks how the social meaning of a group would be changed if its non-modular features were altered.


title: Database of Power Grid Frequency Measurements

abstract: The transformation of the electrical energy system attracts much attention in diverse research communities. Novel approaches for modeling, control, and power grid architectures are widely proposed. Data from actual power system operation are therefore critical to evaluate models and to analyze real-world scenarios. However, they are rarely available. In the present paper, we introduce a precisely time-stamped data set of power grid frequency measurements. It covers twelve synchronous areas of different sizes in one-second resolution. The contribution includes synchronized measurements within the Continental European synchronous area acquired in Portugal, Germany, and Turkey, maximizing the geographical span. Finally, we provide excerpts of the underlying waveform. Data were collected using a self-developed measurement instrument, the Electrical Data Recorder (EDR), connected mostly to conventional power sockets. We close our description with a discussion on measurement error and data quality.


title: Fast and automated biomarker detection in breath samples with machine learning

abstract: Volatile organic compounds (VOCs) in human breath can reveal a large spectrum of health conditions and can be used for fast, accurate and non-invasive diagnostics. Gas chromatography-mass spectrometry (GC-MS) is used to measure VOCs, but its application is limited by expert-driven data analysis that is time-consuming, subjective and may introduce errors. We propose a system to perform GC-MS data analysis that exploits deep learning pattern recognition ability to learn and automatically detect VOCs directly from raw data, thus bypassing expert-led processing. The new proposed approach showed to outperform the expert-led analysis by detecting a significantly higher number of VOCs in just a fraction of time while maintaining high specificity. These results suggest that the proposed method can help the large-scale deployment of breath-based diagnosis by reducing time and cost, and increasing accuracy and consistency.


title: On Lipschitz Normally Embedded complex surface germs

abstract: We undertake a systematic study of Lipschitz Normally Embedded normal complex surface germs. We prove in particular that the topological type of such a germ determines the combinatorics of its minimal resolution which factors through the blowup of its maximal ideal and through its Nash transform, as well as the polar curve and the discriminant curve of a generic plane projection, thus generalizing results of Spivakovsky and Bondil that were known for minimal surface singularities. In an appendix, we give a new example of a Lipschitz Normally Embedded surface singularity.


title: A $p$-adic approach to rational points on curves

abstract: In 1922, Mordell conjectured the striking statement that for a polynomial equation $f(x,y)=0$, if the topology of the set of complex number solutions is complicated enough, then the set of rational number solutions is finite. This was proved by Faltings in 1983, and again by a different method by Vojta in 1991, but neither proof provided a way to provably find all the rational solutions, so the search for other proofs has continued. Recently, Lawrence and Venkatesh found a third proof, relying on variation in families of $p$-adic Galois representations; this is the subject of the present exposition.


title: A generalized spray-flamelet formulation by means of a monotonic variable

abstract: The external structure of the spray-flamelet can be described using the Schvab-Zel'dovich-Li\~nan formulation. The gaseous mixture-fraction variable as function of the physical space, Z(x_i), typically employed for the description of gaseous diffusion flames leads to non-monotonicity behaviour for spray flames due to the extra fuel supplied by vaporisation of droplets distributed into the flow. As a result, the overall properties of spray flames depend not only on Z and the scalar dissipation rate, but also on the spray source term, S_v. We propose a new general coordinate variable which takes into account the spatial information about the entire mixture fraction due to the gaseous phase and droplet vaporisation. This coordinate variable, Z_C(x_i) is based on the cumulative value of the gaseous mixture fraction Z(x_i), and is shown to be monotonic. For pure gaseous flow, the new cumulative function, Z_C, yields the well-established flamelet structure in Z-space. In the present manuscript, the spray-flamelet structure and the new equations for temperature and mass fractions in terms of Z_C are derived and then applied to the canonical counterflow configuration with potential flow. Numerical results are obtained for ethanol and methanol sprays, and the effect of Lewis and Stokes numbers on the spray-flamelet structure are analyzed. The proposed formulation agrees well when mapping the structure back to physical space thereby confirming our integration methodology.


title: Non-destructive characterisation of dopant spatial distribution in cuprate superconductors

abstract: Proper characterisation of investigated samples is vital when studying superconductivity as impurities and doping inhomogeneities can affect the physical properties of the measured system. We present a method where a polarised neutron imaging setup utilises the precession of spin-polarised neutrons in the presence of a trapped field in the superconducting sample to spatially map out the critical temperature for the phase transition between superconducting and non-superconducting states. We demonstrate this method on a superconducting crystal of the prototypical high-temperature superconductor (La,Sr)$_2$CuO$_4$. The results, which are backed up by complementary magnetic susceptibility measurements, show that the method is able to resolve minor variations in the transition temperature across the length of the LSCO crystal, caused by inhomogeneities in strontium doping.


title: Numerical study of non-gyrotropic electron pressure effects in collisionless magnetic reconnection

abstract: We investigate the time evolution of the six-component electron pressure tensor in a hybrid code studying consequences for the two-dimensional reconnection process in an initially perturbed Harris sheet. We put forward that two tensor components (a diagonal and a non-diagonal one) grow in an unstable way unless an isotropization operator is considered. This isotropization term is physically associated with an electron heat flux. As a consequence, we put forward that an enhanced value of a diagonal component is observed in the very middle of field reversal at sub-ion scale. Because of the increase of the kinetic pressure, the magnetic field is decreased in this electron layer, hence increasing the associated out-of-plane current at its edges and leading to its bifurcation. The bifurcation mechanism is based on the presence of electron pressure anisotropy, related to the gradient of inflow electron bulk velocity. The gradient in the inflow direction of the enhanced diagonal electron pressure tensor component results in the deceleration of the ions entering the X-point region. We suggest that bifurcated current sheets resulting from the anisotropies/agyrotropies of the six-component electron pressure tensor correspond to smaller reconnection rates comparing to non-bifurcated ones.


title: An explicit self-dual construction of complete cotorsion pairs in the relative context

abstract: Let $R\to A$ be a homomorphism of associative rings, and let $(\mathcal F,\mathcal C)$ be a hereditary complete cotorsion pair in $R\mathsf{-Mod}$. Let $(\mathcal F_A,\mathcal C_A)$ be the cotorsion pair in $A\mathsf{-Mod}$ in which $\mathcal F_A$ is the class of all left $A$-modules whose underlying $R$-modules belong to $\mathcal F$. Assuming that the $\mathcal F$-resolution dimension of every left $R$-module is finite and the class $\mathcal F$ is preserved by the coinduction functor $\operatorname{Hom}_R(A,-)$, we show that $\mathcal C_A$ is the class of all direct summands of left $A$-modules finitely (co)filtered by $A$-modules coinduced from $R$-modules from $\mathcal C$. Assuming that the class $\mathcal F$ is closed under countable products and preserved by the functor $\operatorname{Hom}_R(A,-)$, we prove that $\mathcal C_A$ is the class of all direct summands of left $A$-modules cofiltered by $A$-modules coinduced from $R$-modules from $\mathcal C$, with the decreasing filtration indexed by the natural numbers. A combined result, based on the assumption that countable products of modules from $\mathcal F$ have finite $\mathcal F$-resolution dimension bounded by $k$, involves cofiltrations indexed by the ordinal $\omega+k$. The dual results also hold, provable by the same technique going back to the author's monograph on semi-infinite homological algebra arXiv:0708.3398. In addition, we discuss the $n$-cotilting and $n$-tilting cotorsion pairs, for which we obtain better results using a suitable version of a classical Bongartz-Ringel lemma. As an illustration of the main results of the paper, we consider certain cotorsion pairs related to the contraderived and coderived categories of curved DG-modules.


title: 6G White Paper on Localization and Sensing

abstract: This white paper explores future localization and sensing opportunities for beyond 5G wireless communication systems by identifying key technology enablers and discussing their underlying challenges, implementation issues, and identifying potential solutions. In addition, we present exciting new opportunities for localization and sensing applications, which will disrupt traditional design principles and revolutionize the way we live, interact with our environment, and do business. Following the trend initiated in the 5G NR systems, 6G will continue to develop towards even higher frequency ranges, wider bandwidths, and massive antenna arrays. In turn, this will enable sensing solutions with very fine range, Doppler and angular resolutions, as well as localization to cm-level degree of accuracy. Moreover, new materials, device types, and reconfigurable surfaces will allow network operators to reshape and control the electromagnetic response of the environment. At the same time, machine learning and artificial intelligence will leverage the unprecedented availability of data and computing resources to tackle the biggest and hardest problems in wireless communication systems. 6G will be truly intelligent wireless systems that will not only provide ubiquitous communication but also empower high accuracy localization and high-resolution sensing services. They will become the catalyst for this revolution by bringing about a unique new set of features and service capabilities, where localization and sensing will coexist with communication, continuously sharing the available resources in time, frequency and space. This white paper concludes by highlighting foundational research challenges, as well as implications and opportunities related to privacy, security, and trust. Addressing these challenges will undoubtedly require an inter-disciplinary and concerted effort from the research community.


title: A Novel Nudity Detection Algorithm for Web and Mobile Application Development

abstract: In our current web and mobile application development runtime nude image content detection is very important. This paper presents a runtime nudity detection method for web and mobile application development. We use two parameters to detect the nude content of an image. One is the number of skin pixels another is face region. A skin color model based on RGB, HSV color spaces are used to detect skin pixels in an image. Google vision api is used to detect the face region. By the percentage of skin regions and face regions an image is identified nude or not. The success of this algorithm exists in detecting skin regions and face regions. The skin detection algorithm can detect skin 95% accurately with a low false-positive rate and the google vision api for web and mobile applications can detect face 99% accurately with less than 1 second time. From the experimental analysis, we have seen that the proposed algorithm can detect 95% percent accurately the nudity of an image.


title: On the Macroscopic limit of Brownian Particles with local interaction

abstract: An interacting particle system made of diffusion processes with local interaction is considered and the macroscopic limit to a nonlinear PDE is investigated. Few rigorous results exists on this problem and in particular the explicit form of the nonlinearity is not known. The paper reviews this subject, some of the main ideas to get the limit nonlinear PDE and provides both heuristic and numerical informations on the precise form of the nonlinearity which are new with respect to the literature and coherent with the few known informations.


title: Temporally-Extended {\epsilon}-Greedy Exploration

abstract: Recent work on exploration in reinforcement learning (RL) has led to a series of increasingly complex solutions to the problem. This increase in complexity often comes at the expense of generality. Recent empirical studies suggest that, when applied to a broader set of domains, some sophisticated exploration methods are outperformed by simpler counterparts, such as {\epsilon}-greedy. In this paper we propose an exploration algorithm that retains the simplicity of {\epsilon}-greedy while reducing dithering. We build on a simple hypothesis: the main limitation of {\epsilon}-greedy exploration is its lack of temporal persistence, which limits its ability to escape local optima. We propose a temporally extended form of {\epsilon}-greedy that simply repeats the sampled action for a random duration. It turns out that, for many duration distributions, this suffices to improve exploration on a large set of domains. Interestingly, a class of distributions inspired by ecological models of animal foraging behaviour yields particularly strong performance.


title: Chronologically dating the early assembly of the Milky Way

abstract: The standard cosmological model ($\Lambda$-CDM) predicts that galaxies are built through hierarchical assembly on cosmological timescales$^{1,2}$. The Milky Way, like other disc galaxies, underwent violent mergers and accretion of small satellite galaxies in its early history. Thanks to Gaia-DR2$^3$ and spectroscopic surveys$^4$, the stellar remnants of such mergers have been identified$^{5-7}$. The chronological dating of such events is crucial to uncover the formation and evolution of the Galaxy at high redshift, but it has so far been challenging owing to difficulties in obtaining precise ages for these oldest stars. Here we combine asteroseismology -- the study of stellar oscillations -- with kinematics and chemical abundances, to estimate precise stellar ages ($\sim$ 11%) for a sample of stars observed by the $\mathit{Kepler}$ space mission$^8$. Crucially, this sample includes not only some of the oldest stars that were formed inside the Galaxy, but also stars formed externally and subsequently accreted onto the Milky Way. Leveraging this resolution in age, we provide compelling evidence in favour of models in which the Galaxy had already formed a substantial population of its stars (which now reside mainly in its thick disc) before the in-fall of the satellite galaxy Gaia-Enceladus/Sausage$^{5,6}$ around 10 billions years ago


title: Coordinating Multiagent Industrial Symbiosis

abstract: We present a formal multiagent framework for coordinating a class of collaborative industrial practices called Industrial Symbiotic Networks (ISNs) as cooperative games. The game-theoretic formulation of ISNs enables systematic reasoning about what we call the ISN implementation problem. Specifically, the characteristics of ISNs may lead to the inapplicability of standard fair and stable benefit allocation methods. Inspired by realistic ISN scenarios and following the literature on normative multiagent systems, we consider regulations and normative socio-economic policies as coordination instruments that in combination with ISN games resolve the situation. In this multiagent system, employing Marginal Contribution Nets (MC-Nets) as rule-based cooperative game representations foster the combination of regulations and ISN games with no loss in expressiveness. We develop algorithmic methods for generating regulations that ensure the implementability of ISNs and as a policy support, present the policy requirements that guarantee the implementability of all the desired ISNs in a balanced-budget way.


title: Geometric Graph Representations and Geometric Graph Convolutions for Deep Learning on Three-Dimensional (3D) Graphs

abstract: The geometry of three-dimensional (3D) graphs, consisting of nodes and edges, plays a crucial role in many important applications. An excellent example is molecular graphs, whose geometry influences important properties of a molecule including its reactivity and biological activity. To facilitate the incorporation of geometry in deep learning on 3D graphs, we define three types of geometric graph representations: positional, angle-geometric and distance-geometric. For proof of concept, we use the distance-geometric graph representation for geometric graph convolutions. Further, to utilize standard graph convolution networks, we employ a simple edge weight / edge distance correlation scheme, whose parameters can be fixed using reference values or determined through Bayesian hyperparameter optimization. The results of geometric graph convolutions, for the ESOL and Freesol datasets, show significant improvement over those of standard graph convolutions. Our work demonstrates the feasibility and promise of incorporating geometry, using the distance-geometric graph representation, in deep learning on 3D graphs.


title: Hyperparameter Selection for Subsampling Bootstraps

abstract: Massive data analysis becomes increasingly prevalent, subsampling methods like BLB (Bag of Little Bootstraps) serves as powerful tools for assessing the quality of estimators for massive data. However, the performance of the subsampling methods are highly influenced by the selection of tuning parameters ( e.g., the subset size, number of resamples per subset ). In this article we develop a hyperparameter selection methodology, which can be used to select tuning parameters for subsampling methods. Specifically, by a careful theoretical analysis, we find an analytically simple and elegant relationship between the asymptotic efficiency of various subsampling estimators and their hyperparameters. This leads to an optimal choice of the hyperparameters. More specifically, for an arbitrarily specified hyperparameter set, we can improve it to be a new set of hyperparameters with no extra CPU time cost, but the resulting estimator's statistical efficiency can be much improved. Both simulation studies and real data analysis demonstrate the superior advantage of our method.


title: Global well-posedness for the three dimensional Muskat problem in the critical Sobolev space

abstract: We prove that the 3D stable Muskat problem is globally well-posed in the critical Sobolev space $\dot H^2 \cap \dot W^{1,\infty}$ provided that the semi-norm $\Vert f_0 \Vert_{\dot H^{2}}$ is small enough. Consequently, this allows the Lipschitz semi-norm to be arbitrarily large. The proof is based on a new formulation of the 3D Muskat problem that allows to capture the hidden oscillatory nature of the problem. The latter formulation allows to prove the $\dot H^{2}$ {\emph{a priori}} estimates. In the literature, all the known global existence results for the 3D Muskat problem are for small slopes (less than 1). This is the first arbitrary large slope theorem for the 3D stable Muskat problem.


title: A Mechanism of Spark Motion in Inner Acceleration Region to Investigate Subpulse Drifting in Pulsars

abstract: Coherent radio emission in pulsars is excited due to instabilities in a relativistically streaming non-stationary plasma flow, which is generated from sparking discharges in the inner acceleration region (IAR) near the stellar surface. A number of detailed works have shown the IAR to be a partially screened gap (PSG) dominated by non-dipolar magnetic fields with continuous outflow of ions from the surface. The phenomenon of subpulse drifting is expected to originate due to variable $\mathbf{E}\times\mathbf{B}$ drift of the sparks in PSG, where the sparks lag behind corotation velocity of the pulsar. Detailed observations show a wide variety of subpulse drifting behaviour where subpulses in different components of the profile have different phase trajectories. But the drifting periodicity is seen to be constant, within measurement errors, across all components of the profile. Using the concept of sparks lagging behind corotation speed in PSG as well as the different orientations of the surface non-dipolar magnetic fields we have simulated the expected single pulse behaviour in a representative sample of pulsars. Our results show that the different types of drifting phase behaviour can be reproduced using these simple assumptions of spark dynamics in a non-dipolar IAR.


title: A probabilistic generative model for semi-supervised training of coarse-grained surrogates and enforcing physical constraints through virtual observables

abstract: The data-centric construction of inexpensive surrogates for fine-grained, physical models has been at the forefront of computational physics due to its significant utility in many-query tasks such as uncertainty quantification. Recent efforts have taken advantage of the enabling technologies from the field of machine learning (e.g. deep neural networks) in combination with simulation data. While such strategies have shown promise even in higher-dimensional problems, they generally require large amounts of training data even though the construction of surrogates is by definition a Small Data problem. Rather than employing data-based loss functions, it has been proposed to make use of the governing equations (in the simplest case at collocation points) in order to imbue domain knowledge in the training of the otherwise black-box-like interpolators. The present paper provides a flexible, probabilistic framework that accounts for physical structure and information both in the training objectives as well as in the surrogate model itself. We advocate a probabilistic (Bayesian) model in which equalities that are available from the physics (e.g. residuals, conservation laws) can be introduced as virtual observables and can provide additional information through the likelihood. We further advocate a generative model i.e. one that attempts to learn the joint density of inputs and outputs that is capable of making use of unlabeled data (i.e. only inputs) in a semi-supervised fashion in order to promote the discovery of lower-dimensional embeddings which are nevertheless predictive of the fine-grained model's output.


title: Depth-Optimized Delay-Aware Tree (DO-DAT) for Virtual Network Function Placement

abstract: With the constant increase in demand for data connectivity, network service providers are faced with the task of reducing their capital and operational expenses while ensuring continual improvements to network performance. Although Network Function Virtualization (NFV) has been identified as a solution, several challenges must be addressed to ensure its feasibility. In this paper, we present a machine learning-based solution to the Virtual Network Function (VNF) placement problem. This paper proposes the Depth-Optimized Delay-Aware Tree (DO-DAT) model by using the particle swarm optimization technique to optimize decision tree hyper-parameters. Using the Evolved Packet Core (EPC) as a use case, we evaluate the performance of the model and compare it to a previously proposed model and a heuristic placement strategy.


title: SaliencyMix: A Saliency Guided Data Augmentation Strategy for Better Regularization

abstract: Advanced data augmentation strategies have widely been studied to improve the generalization ability of deep learning models. Regional dropout is one of the popular solutions that guides the model to focus on less discriminative parts by randomly removing image regions, resulting in improved regularization. However, such information removal is undesirable. On the other hand, recent strategies suggest to randomly cut and mix patches and their labels among training images, to enjoy the advantages of regional dropout without having any pointless pixel in the augmented images. We argue that such random selection strategies of the patches may not necessarily represent sufficient information about the corresponding object and thereby mixing the labels according to that uninformative patch enables the model to learn unexpected feature representation. Therefore, we propose SaliencyMix that carefully selects a representative image patch with the help of a saliency map and mixes this indicative patch with the target image, thus leading the model to learn more appropriate feature representation. SaliencyMix achieves the best known top-1 error of 21.26% and 20.09% for ResNet-50 and ResNet-101 architectures on ImageNet classification, respectively, and also improves the model robustness against adversarial perturbations. Furthermore, models that are trained with SaliencyMix help to improve the object detection performance. Source code is available at https://github.com/SaliencyMix/SaliencyMix.


title: Intrinsically weak magnetic anisotropy of cerium in potential hard-magnetic intermetallics

abstract: Cerium-based intermetallics are currently attracting much interest as a possible alternative to existing high-performance magnets containing scarce heavy rare-earth elements. However, the intrinsic magnetic properties of Ce in these systems are poorly understood due to the difficulty of a quantitative description of the Kondo effect, a many-body phenomenon where conduction electrons screen out the Ce-4f moment. Here, we show that the Ce-4f shell in Ce-Fe intermetallics is partially Kondo screened. The Kondo scale is dramatically enhanced by nitrogen interstitials completely suppressing the Ce-4f contribution to the magnetic anisotropy, in striking contrast to the effect of nitrogenation in isostructural intermetallics containing other rare-earth elements. We determine the full temperature dependence of the Ce-4f single-ion anisotropy and show that even unscreened Ce-4f moments contribute little to the room-temperature intrinsic magnetic hardness. Our study thus establishes fundamental constraints on the potential of cerium-based permanent magnet intermetallics.


title: Inhomogeneous spacetimes in Weyl integrable geometry with matter source

abstract: We investigate the existence of inhomogeneous exact solutions in Weyl Integrable theory with a matter source. In particular, we consider the existence of a dust fluid source while for the underlying geometry we assume a line element which belongs to the family of silent universes. We solve explicitly the field equations and we find the Szekeres spacetimes in Weyl Integrable theory. We show that only the isotropic family can describe inhomogeneous solutions where the LTB spacetimes are included. A detailed analysis of the dynamics of the field equations is given where the past and future attractors are determined. It is interesting that the Kasner spacetimes can be seen as past attractors for the gravitation models, while the unique future attractor describes the Milne universe similar with the behaviour of the gravitational model in the case of General Relativity.


title: Number Theory and Axiomatic Geometry in the Diproche System

abstract: Diproche ("Didactical Proof Checking") is an automatic system for supporting the acquistion of elementary proving skills in the initial phase of university education in mathematics. A key feature of Diproche - which is designed by the example of the Naproche system developed by M. Cramer and others - is an automated proof checker for proofs written in a controlled fragment of natural language specifically designed to capture the language of beginners' proving exercises in mathematics. Both the accepted language and proof methods depend on the didactical and mathematical context and vary with the level of education and the topic proposed. An overall presentation of the system in general was given in Carl and Krapf 2019. Here, we briefly recall the basic architecture of Diproche and then focus on explaining key features and the working principles of Diproche in the sample topics of elementary number theory and axiomatic geometry.


title: Shapley Value as Principled Metric for Structured Network Pruning

abstract: Structured pruning is a well-known technique to reduce the storage size and inference cost of neural networks. The usual pruning pipeline consists of ranking the network internal filters and activations with respect to their contributions to the network performance, removing the units with the lowest contribution, and fine-tuning the network to reduce the harm induced by pruning. Recent results showed that random pruning performs on par with other metrics, given enough fine-tuning resources. In this work, we show that this is not true on a low-data regime when fine-tuning is either not possible or not effective. In this case, reducing the harm caused by pruning becomes crucial to retain the performance of the network. First, we analyze the problem of estimating the contribution of hidden units with tools suggested by cooperative game theory and propose Shapley values as a principled ranking metric for this task. We compare with several alternatives proposed in the literature and discuss how Shapley values are theoretically preferable. Finally, we compare all ranking metrics on the challenging scenario of low-data pruning, where we demonstrate how Shapley values outperform other heuristics.


title: Neural Speaker Diarization with Speaker-Wise Chain Rule

abstract: Speaker diarization is an essential step for processing multi-speaker audio. Although an end-to-end neural diarization (EEND) method achieved state-of-the-art performance, it is limited to a fixed number of speakers. In this paper, we solve this fixed number of speaker issue by a novel speaker-wise conditional inference method based on the probabilistic chain rule. In the proposed method, each speaker's speech activity is regarded as a single random variable, and is estimated sequentially conditioned on previously estimated other speakers' speech activities. Similar to other sequence-to-sequence models, the proposed method produces a variable number of speakers with a stop sequence condition. We evaluated the proposed method on multi-speaker audio recordings of a variable number of speakers. Experimental results show that the proposed method can correctly produce diarization results with a variable number of speakers and outperforms the state-of-the-art end-to-end speaker diarization methods in terms of diarization error rate.


title: Object-Independent Human-to-Robot Handovers using Real Time Robotic Vision

abstract: We present an approach for safe and object-independent human-to-robot handovers using real time robotic vision and manipulation. We aim for general applicability with a generic object detector, a fast grasp selection algorithm and by using a single gripper-mounted RGB-D camera, hence not relying on external sensors. The robot is controlled via visual servoing towards the object of interest. Putting a high emphasis on safety, we use two perception modules: human body part segmentation and hand/finger segmentation. Pixels that are deemed to belong to the human are filtered out from candidate grasp poses, hence ensuring that the robot safely picks the object without colliding with the human partner. The grasp selection and perception modules run concurrently in real-time, which allows monitoring of the progress. In experiments with 13 objects, the robot was able to successfully take the object from the human in 81.9% of the trials.


title: Complex-valued (p,q)-harmonic morphisms from Riemannian manifolds

abstract: We introduce the natural notion of (p,q)-harmonic morphisms between Riemannian manifolds. This unifies several theories that have been studied during the last decades. We then study the special case when the maps involved are complex-valued. For these we find a characterisation and provide new non-trivial examples in important cases.


title: The role of exchangeability in causal inference

abstract: Though the notion of exchangeability has been discussed in the causal inference literature under various guises, it has rarely taken its original meaning as a symmetry property of probability distributions. As this property is a standard component of Bayesian inference, we argue that in Bayesian causal inference it is natural to link the causal model, including the notion of confounding and definition of causal contrasts of interest, to the concept of exchangeability. Here we propose a probabilistic between-group exchangeability property as an identifying condition for causal effects, relate it to alternative conditions for unconfounded inferences (commonly stated using potential outcomes) and define causal contrasts in the presence of exchangeability in terms of posterior predictive expectations for further exchangeable units. While our main focus is on a point treatment setting, we also investigate how this reasoning carries over to longitudinal settings.


title: Automatized Evaluation of Formalization Exercises in Mathematics

abstract: We describe two systems for supporting beginner students in acquiring basic skills in expressing statements in the formalism of first-order predicate logic; the first, called "math dictations", presents users with the task of formalizing a given natural-language sentence, while the second, called "Game of Def", challenges users to give a formal description of a set of a geometric pattern displayed to them. In both cases, an automatic checking takes place.


title: Variational optimization of continuous matrix product states

abstract: Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials, and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.


title: Robust Multiple Stopping -- A Pathwise Duality Approach

abstract: We develop a method to solve, theoretically and numerically, general optimal stopping problems. Our general setting allows for multiple exercise rights, i.e., optimal multiple stopping, for a robust evaluation that accounts for model uncertainty, and for general reward processes driven by multi-dimensional jump-diffusions. Our approach relies on first establishing robust martingale dual representation results for the multiple stopping problem that satisfy appealing pathwise optimality (i.e., almost sure) properties. Next, we exploit these theoretical results to develop upper and lower bounds that, as we formally show, not only converge to the true solution asymptotically, but also constitute genuine pre-limiting upper and lower bounds. We illustrate the applicability of our approach in a few examples and analyze the impact of model uncertainty on optimal multiple stopping strategies.


title: Compressed Sensing Tomography for qudits in Hilbert spaces of non-power-of-two dimensions

abstract: The techniques of low-rank matrix recovery were adapted for Quantum State Tomography (QST) previously by D. Gross et al. [Phys. Rev. Lett. 105, 150401 (2010)], where they consider the tomography of $n$ spin-$1/2$ systems. For the density matrix of dimension $d = 2^n$ and rank $r$ with $r \ll 2^n$, it was shown that randomly chosen Pauli measurements of the order $O(dr \log(d)^2)$ are enough to fully reconstruct the density matrix by running a specific convex optimization algorithm. The result utilized the low operator-norm of the Pauli operator basis, which makes it `incoherent' to low-rank matrices. For quantum systems of dimension $d$ not a power of two, Pauli measurements are not available, and one may consider using SU($d$) measurements. Here, we point out that the SU($d$) operators, owing to their high operator norm, do not provide a significant savings in the number of measurement settings required for successful recovery of all rank-$r$ states. We propose an alternative strategy, in which the quantum information is swapped into the subspace of a power-two system using only $\textrm{poly}(\log(d)^2)$ gates at most, with QST being implemented subsequently by performing $O(dr \log(d)^2)$ Pauli measurements. We show that, despite the increased dimensionality, this method is more efficient than the one using SU($d$) measurements.


title: Practical sensorless aberration estimation for 3D microscopy with deep learning

abstract: Estimation of optical aberrations from volumetric intensity images is a key step in sensorless adaptive optics for 3D microscopy. Recent approaches based on deep learning promise accurate results at fast processing speeds. However, collecting ground truth microscopy data for training the network is typically very difficult or even impossible thereby limiting this approach in practice. Here, we demonstrate that neural networks trained only on simulated data yield accurate predictions for real experimental images. We validate our approach on simulated and experimental datasets acquired with two different microscopy modalities, and also compare the results to non-learned methods. Additionally, we study the predictability of individual aberrations with respect to their data requirements and find that the symmetry of the wavefront plays a crucial role. Finally, we make our implementation freely available as open source software in Python.


title: Scalable quantum processor noise characterization

abstract: Measurement fidelity matrices (MFMs) (also called error kernels) are a natural way to characterize state preparation and measurement errors in near-term quantum hardware. They can be employed in post processing to mitigate errors and substantially increase the effective accuracy of quantum hardware. However, the feasibility of using MFMs is currently limited as the experimental cost of determining the MFM for a device grows exponentially with the number of qubits. In this work we present a scalable way to construct approximate MFMs for many-qubit devices based on cumulant expansion. Our method can also be used to characterize various types of correlation error.


title: OSSOS XX: The Meaning of Kuiper Belt Colors

abstract: Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.


title: The Fate of Binary Stars Hosting Planets upon Interaction with Sgr A* Black hole

abstract: Our Galaxy hosts a very massive object at its centre, often referred to as the supermassive black hole Sgr A*. Its gravitational tidal field is so intense that can strip apart a binary star passing its vicinity and accelerate one of the components of the binary as hypervelocity star (HVS) and grab the other star as S-star. Taken into consideration that many binary star systems are known to host planets, in this paper we aim to broaden the study of the close interaction of binary stars and their planetary systems with Sgr A* massive object. Results are obtained via a high precision $N-$body code including post-Newtonian approximation. We quantify the likelihood of capture and ejection of stars and planets after interaction with Sgr A*, finding that the fraction of stars captured around it is about three times that of the planets (~ 49.4% versus ~14.5%) and the fraction of hypervelocity planet (HVP) ejection is about twice HVSs (~21.7% versus ~9.0%). The actual possibility of observational counterparts deserves further investigation.


title: Extractive contest design

abstract: We consider contest success functions (CSFs) that extract contestants' prize values. In the common-value case, there exists a CSF extractive in any equilibrium. In the observable-private-value case, there exists a CSF extractive in some equilibrium; there exists a CSF extractive in any equilibrium if and only if the number of contestants is greater than or equal to three or the values are homogeneous. In the unobservable-private-value case, there exists no CSF extractive in some equilibrium. When extractive CSFs exist, we explicitly present one of them.


title: Ubiquitous hundred-Gauss magnetic fields in solar spicules

abstract: Aims. We use high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the SST to study the magnetic field in solar spicules. Methods. The equations that result from the application of the Weak Field Approximation (WFA) to the radiative transfer equations are used to infer the LOS component of the magnetic field (BLOS). Two restrictive conditions are imposed on the Stokes I and V profiles at each pixel before they can be used in a Bayesian inversion to compute its BLOS. Results. The LOS magnetic field component has been inferred in six data sets totalling 448 spectral scans in the Ca ii 8542 A line and containing both active region and quiet Sun areas, with values of hundreds of G being abundantly inferred. There seems to be no difference, from the statistical point of view, between the magnetic field strength of spicules in the quiet Sun or near an active region. On the other hand, the BLOS distributions present smaller values on the disk than off-limb, a fact that can be explained by the effect of superposition on the chromosphere of on-disk structures. We find that spicules in the vicinity of a sunspot have a magnetic field polarity (i.e. north or south) equal to that of the sunspot. This paper also contains an analysis of the effect of off-limb overlapping structures on the observed Stokes I and V parameters and the BLOS obtained from the WFA. It is found that this value is equal to or smaller than the largest LOS magnetic field components of the two structures. In addition, using random BLOS, Doppler velocities and line intensities of these two structures leads in ~ 50% of the cases to Stokes I and V parameters unsuitable to be used with the WFA. Conclusions. Our results present a scarcity of LOS magnetic field components smaller than some 50 G, which must not be taken as evidence against the existence of these magnetic field strengths in spicules.


title: Motive of the $SL_4$-character variety of torus knots

abstract: In this paper, we compute the motive of the character variety of representations of the fundamental group of the complement of an arbitrary torus knot into $SL_4(k)$, for any algebraically closed field $k$ of zero characteristic. For that purpose, we introduce a stratification of the variety in terms of the type of a canonical filtration attached to any representation. This allows us to reduce the computation of the motive to a combinatorial problem.


title: The effects on CMB power spectra and bispectra from the polarization rotation and its correlations with temperature and E-polarization

abstract: The Chern-Simons term, through which the cosmic Axion-like field couples to the electromagnetic field, has the effect to rotate CMB polarization directions and to break the CPT symmetry. This rotation will change the CMB power spectra, no matter isotropic or anisotropic the rotation angle is. In this paper we revisit this issue by further considering the correlations between the (anisotropic) rotation angle $\alpha$ and the CMB temperature and (unrotated) $E$ polarization fields. These correlations could be generated in the Axion-like models with nonzero potential under the adiabatic initial condition. We first investigate how these correlations contribute further modifications to the CMB power spectra, then calculate the CMB bispectra for the temperature and rotated polarization fields. These bispectra would vanish if the $T\alpha$ and $E\alpha$ correlations are absent. So, they are useful in searching for CPT violation and the $T\alpha$ and $E\alpha$ correlations arisen in the Axion-like models.


title: Identifying groomed jet splittings in heavy-ion collisions

abstract: Measurements of jet substructure in heavy-ion collisions may provide key insight to the nature of jet quenching in the quark-gluon plasma. Jet grooming techniques from high-energy physics have been applied to heavy-ion collisions in order to isolate theoretically controlled jet observables and explore possible modification to the hard substructure of jets. However, the grooming algorithms used have not been tailored to the unique considerations of heavy-ion collisions, in particular to the experimental challenge of reconstructing jets in the presence of a large underlying event. We report a set of simple studies illustrating the impact of the underlying event on identifying groomed jet splittings in heavy-ion collisions, and on associated groomed jet observables. We illustrate the importance of the selection of grooming algorithm, as certain groomers are more robust to these effects, while others, including those commonly used in heavy-ion collisions, are susceptible to large background effects -- which, when uncontrolled, can mimic a jet quenching signal. These experimental considerations, along with appropriate theoretical motivation, provide input to the choice of grooming algorithms employed in heavy-ion collisions.


title: Robust and Precision Satellite Formation Flying Guidance Using Adaptive Optimal Control Techniques

abstract: The main focus of the work presented in this thesis is to develop an optimal control based formation flying control strategy for high precision formation flying of small satellites that have restricted computation and storage capacity. Using the recently developed model predictive static programming (MPSP), and Generalized MPSP algorithm a suboptimal guidance logic is presented for the formation flying of small satellites. The proposed guidance scheme is valid both for high eccentricity chief satellite orbits as well as the large separation distance between the chief and deputy satellites. Comparative study with standard Linear Quadratic Regulator (LQR) solution (which serves as a guess solution for MPSP) and another nonlinear controller, Finite-time State-Dependent Ricatti Equation (SDRE) reveals that MPSP guidance achieves the objective with higher accuracy and with a lesser amount of control usage.


title: Influence of the Galactic gravitational field on the positional accuracy of extragalactic sources. II Observational appearances and detectability

abstract: We consider a possibility of detecting the jitter effect of apparent celestial positions of distant sources due to local fluctuations of the Galaxy gravitational field. It is proposed to observe two samples of extragalactic sources (target and control) in different sky directions using the high-precision radio interferometry. It is shown that on a scale of ~2 years, it is possible to detect a systematic increase in the standard deviation of measured arc lengths of pairs of target sources compared to the control ones at the $3\sigma$-level if the accuracy of differential astrometric observations is around 10 $\mu$as. For the current state-of-the-art accuracy of 30 $\mu$as achieved at the KVN or VERA interferometers, which have shorter baselines in comparison with VLBI, the target and control samples will differ only at the 2$\sigma$-level on the scale of 10 years. To achieve the $3\sigma$-level on this time interval, it is necessary to improve the accuracy up to ~20 $\mu$as. Other possible effects that can also affect the arc length measurements between two sources are discussed, and an observational strategy to minimize them is suggested.


title: Strain-Driven Structure-Ferroelectricity Relationship in hexagonal TbMnO$_3$ Films

abstract: Thin films and heterostructures of hexagonal manganites as promising multiferroic materials have attracted a considerable interest recently. We report structural transformations of high quality epitaxial h-TMO/YSZ(111) films, analyzed by means of various characterization techniques. A phase transition from P63mc to P63mcm structure at TC~800 K was observed by temperature dependent Raman spectroscopy and optical ellipsometry. The latter probing directly electronic system, indicates its modification at the structural phase transition likely due to charge transfer from oxygen to Mn. In situ transmission electron microscopy (TEM) of the lamella samples displayed an irreversible P63mc-P63mcm transformation and vanishing of ferroelectric domains already at 410 K. After the temperature cycling (300K-1300K-300K) the room temperature TEM of h-TMO films revealed an inhomogeneous microstructure, containing ferroelectric and paraelectric nanodomains with P63mc and P63mcm structure, respectively. We point out a strong influence of stress relaxation, induced by temperature and by constrained sample geometry onto the structure and ferroelectricity in strain-stabilized h-TMO thin films.


title: Age-Based Coded Computation for Bias Reduction in Distributed Learning

abstract: Coded computation can be used to speed up distributed learning in the presence of straggling workers. Partial recovery of the gradient vector can further reduce the computation time at each iteration; however, this can result in biased estimators, which may slow down convergence, or even cause divergence. Estimator bias will be particularly prevalent when the straggling behavior is correlated over time, which results in the gradient estimators being dominated by a few fast servers. To mitigate biased estimators, we design a $timely$ dynamic encoding framework for partial recovery that includes an ordering operator that changes the codewords and computation orders at workers over time. To regulate the recovery frequencies, we adopt an $age$ metric in the design of the dynamic encoding scheme. We show through numerical results that the proposed dynamic encoding strategy increases the timeliness of the recovered computations, which as a result, reduces the bias in model updates, and accelerates the convergence compared to the conventional static partial recovery schemes.


title: Path Integral Approach to Uncertainties in SIR-type Systems

abstract: In this paper I show how path integral techniques can be used to put measures on histories in "susceptible-infectious-recovered" (SIR)-type systems. The standard SIR solution emerges as the classical saddle point of the action describing the measure. One can then expand perturbatively around the background solution, and this paper goes on to work out the covariance of fluctuations around the background solution. Using a Green's function type approach, one simply needs to solve additional ordinary differential equations; an explicit matrix inversion is not required. The computed covariance matrix should be useful in the construction of fast likelihoods for fitting the parameters of SIR-type models to data. A comparison of the predictions of the approach to an ensemble of simulations is presented.


title: Securing Your Collaborative Jupyter Notebooks in the Cloud using Container and Load Balancing Services

abstract: Jupyter has become the go-to platform for developing data applications but data and security concerns, especially when dealing with healthcare, have become paramount for many institutions and applications dealing with sensitive information. How then can we continue to enjoy the data analysis and machine learning opportunities provided by Jupyter and the Python ecosystem while guaranteeing auditable compliance with security and privacy concerns? We will describe the architecture and implementation of a cloud based platform based on Jupyter that integrates with Amazon Web Services (AWS) and uses containerized services without exposing the platform to the vulnerabilities present in Kubernetes and JupyterHub. This architecture addresses the HIPAA requirements to ensure both security and privacy of data. The architecture uses an AWS service to provide JSON Web Tokens (JWT) for authentication as well as network control. Furthermore, our architecture enables secure collaboration and sharing of Jupyter notebooks. Even though our platform is focused on Jupyter notebooks and JupyterLab, it also supports R-Studio and bespoke applications that share the same authentication mechanisms. Further, the platform can be extended to other cloud services other than AWS.


title: Carath\'eodory Sampling for Stochastic Gradient Descent

abstract: Many problems require to optimize empirical risk functions over large data sets. Gradient descent methods that calculate the full gradient in every descent step do not scale to such datasets. Various flavours of Stochastic Gradient Descent (SGD) replace the expensive summation that computes the full gradient by approximating it with a small sum over a randomly selected subsample of the data set that in turn suffers from a high variance. We present a different approach that is inspired by classical results of Tchakaloff and Carath\'eodory about measure reduction. These results allow to replace an empirical measure with another, carefully constructed probability measure that has a much smaller support, but can preserve certain statistics such as the expected gradient. To turn this into scalable algorithms we firstly, adaptively select the descent steps where the measure reduction is carried out; secondly, we combine this with Block Coordinate Descent so that measure reduction can be done very cheaply. This makes the resulting methods scalable to high-dimensional spaces. Finally, we provide an experimental validation and comparison.


title: Private 5G: The Future of Industrial Wireless

abstract: High-performance wireless communication is crucial in digital transformation of industrial systems which is driven by Industry 4.0 and the Industrial Internet initiatives. Among the candidate industrial wireless technologies, 5G (cellular/mobile) holds significant potential. Operation of private (non-public) 5G networks in industrial environments is promising to fully unleash this potential. This article provides a technical overview of private 5G networks. It introduces the concept and functional architecture of private 5G while highlighting the key benefits and industrial use-cases. It explores spectrum opportunities for private 5G networks. It also discusses design aspects of private 5G along with the key challenges. Finally, it explores the emerging standardization and open innovation ecosystem for private 5G.


title: Hierarchical Schr\"{o}dinger-type operators: the case of potentials with local singularities

abstract: The goal of this paper is twofold. We prove that the operator $H=L+V$ , a perturbation of the Taibleson-Vladimirov multiplier $L=\mathfrak{D}^{\alpha}$ by a potential $V(x)=b\left\Vert x\right\Vert ^{-\alpha},$ $b\geq b_{\ast},$ is essentially self-adjoint and non-negative definite (the critical value $b_{\ast}$ depends on $\alpha$ and will be specified later). While the operator $H$ is non-negative definite the potential $V(x)$ may well take negative values, e.g. $b_{\ast}<0$ for all $0<\alpha<1$. The equation $Hu=v$ admiits a Green function $g_{H}(x,y)$, the integral kernel of the operator $H^{-1}$. We obtain sharp lower- and upper bounds on the ratio of the functions $g_{H}(x,y)$ and $g_{L}(x,y)$. Examples illustrate our exposition.


title: Inflationary gravitational waves and exotic pre Big Bang Nucleosynthesis cosmology

abstract: According to the most popular scenario, the early Universe should have experienced an accelerated expansion phase, called Cosmological Inflation, after which the standard Big Bang Cosmology would have taken place giving rise to the radiation-dominated epoch. However, the details of the inflationary scenario are far to be completely understood. Thus, in this paper we study if possible additional (exotic) cosmological phases could delay the beginning of the standard Big Bang history and alter some theoretical predictions related to the inflationary cosmological perturbations, like, for instance, the order of magnitude of the tensor-to-scalar ratio $r$.


title: Parallel single-shot measurement and coherent control of solid-state spins below the diffraction limit

abstract: Solid-state spin defects are a promising platform for quantum science and technology, having realized demonstrations of a variety of key components for quantum information processing, particularly in the area of quantum networks. An outstanding challenge for building larger-scale quantum systems with solid-state defects is realizing high-fidelity control over multiple defects with nanoscale separations, which is required to realize strong spin-spin interactions for multi-qubit logic and the creation of entangled states. In this work, we experimentally demonstrate an optical frequency-domain multiplexing technique, allowing high-fidelity initialization and single-shot spin measurement of six rare earth (Er$^{3+}$) ions, within the sub-wavelength volume of a single, silicon photonic crystal cavity. We also demonstrate sub-wavelength control over coherent spin rotations using an optical AC Stark shift. The demonstrated approach may be scaled to large numbers of ions with arbitrarily small separation, and is a significant step towards realizing strongly interacting atomic defect arrays with applications to quantum information processing and fundamental studies of many-body dynamics.


title: Minimal and nearly minimal measure expansions in connected unimodular groups

abstract: Let $G$ be a connected unimodular group equipped with a (left and hence right) Haar measure $\mu_G$, and suppose $A, B \subseteq G$ are nonempty and compact. An inequality by Kemperman gives us $\mu_G(AB)\geq\min\{\mu_G(A)+\mu_G(B),\mu_G(G)\}.$ Our first result determines the conditions for the equality to hold, providing a complete answer to a question asked by Kemperman in 1964. Our second result characterizes compact and connected $G$, $A$, and $B$ that nearly realize equality, with quantitative bounds having the sharp exponent. This can be seen up-to-constant as a $(3k-4)$-theorem for this setting and confirms the connected case of conjectures by Griesmer and by Tao. As an application, we get a measure expansion gap result for connected compact simple Lie groups. The tools developed in our proof include an analysis of the shape of minimally and nearly minimally expanding pairs of sets, a bridge from this to the properties of a certain pseudometric, and a construction of appropriate continuous group homomorphisms to either $\mathbb{R}$ or $\mathbb{T} = \mathbb{R}/\mathbb{Z}$ from the pseudometric.


title: Efficient tree-structured categorical retrieval

abstract: We study a document retrieval problem in the new framework where $D$ text documents are organized in a {\em category tree} with a pre-defined number $h$ of categories. This situation occurs e.g. with taxomonic trees in biology or subject classification systems for scientific literature. Given a string pattern $p$ and a category (level in the category tree), we wish to efficiently retrieve the $t$ \emph{categorical units} containing this pattern and belonging to the category. We propose several efficient solutions for this problem. One of them uses $n(\log\sigma(1+o(1))+\log D+O(h)) + O(\Delta)$ bits of space and $O(|p|+t)$ query time, where $n$ is the total length of the documents, $\sigma$ the size of the alphabet used in the documents and $\Delta$ is the total number of nodes in the category tree. Another solution uses $n(\log\sigma(1+o(1))+O(\log D))+O(\Delta)+O(D\log n)$ bits of space and $O(|p|+t\log D)$ query time. We finally propose other solutions which are more space-efficient at the expense of a slight increase in query time.


title: Observational Support for Massive Black Hole Formation Driven by Runaway Stellar Collisions in Galactic Nuclei

abstract: We explore here an scenario for massive black hole formation driven by stellar collisions in galactic nuclei, proposing a new formation regime of global instability in nuclear stellar clusters triggered by runaway stellar collisions. Using order of magnitude estimations, we show that observed nuclear stellar clusters avoid the regime where stellar collisions are dynamically relevant over the whole system, while resolved detections of massive black holes are well into such collision-dominated regime. We interpret this result in terms of massive black holes and nuclear stellar clusters being different evolutionary paths of a common formation mechanism, unified under the standard terminology of being both central massive objects. We propose a formation scenario where central massive objects more massive than $\rm \sim 10^8 \, Msun$, which also have relaxation times longer that their collision times, will be too dense (in virial equilibrium) to be globally stable against stellar collisions and most of its mass will collapse towards the formation of a massive black hole. Contrarily, this will only be the case at the core of less dense central massive objects leading to the formation of black holes with much lower black hole efficiencies $\rm \epsilon_{BH} = \frac{M_{BH}}{M_{CMO}}$, with these efficiencies $\rm \epsilon_{BH}$ drastically growing for central massive objects more massive than $\rm \sim 10^7 \, Msun$, approaching unity around $\rm M_{CMO} \sim 10^8 \, Msun$. We show that the proposed scenario successfully explains the relative trends observed in the masses, efficiencies, and scaling relations between massive black holes and nuclear stellar clusters.

categories: astro-ph.GA astro-ph.CO astro-ph.HE gr-qc


title: Hydrodynamic Simulations of Classical Novae; CO and ONe White Dwarfs are Supernova Ia Progenitors

abstract: Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of the thermonuclear runaway (TNR) in the accreted material that ends in a Classical Nova explosion. We have followed the evolution of the TNRs on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs. We report on 3 studies in this paper. First, simulations in which we accrete only solar matter using NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar matter and compare the results. Third, we accrete solar matter until the TNR is ongoing and then switch the composition in the accreted layers to a mixed composition: either 25% WD and 75% solar or 50% WD and 50% Solar matter. We find that the amount of accreted material is inversely proportional to the initial 12C abundance (as expected). Thus, accreting solar matter results in a larger amount of accreted material to fuel the outburst; much larger than in earlier studies where a mixed composition was assumed from the beginning of the simulation. Our most important result is that all these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit.


title: Ubiquity of Fourier Transformation in Optical Sciences

abstract: The Fourier transform operation is an important conceptual as well as computational tool in the arsenal of every practitioner of physical and mathematical sciences. We discuss some of its applications in optical science and engineering, with the goal of providing a broad perspective on the intimate relation between the physical and mathematical concepts that are elegantly interwoven within the theory of Fourier transforms.


title: Solution to the fractional equation with left-sided fractional Bessel derivatives of Gerasimov-Caputo type

abstract: In the article we propose and study a method to solve ordinary differential equations with left-sided fractional Bessel derivatives on semi-axes of Gerasimov-Caputo type. We derive explicit solutions to equations with fractional powers of Bessel operator using the Meijer integral transform.


title: Tangles: a structural approach to artificial intelligence in the empirical sciences (Part I)

abstract: Traditional clustering identifies groups of objects that share certain qualities. Tangles do the converse: they identify groups of qualities that often occur together. They can thereby discover, relate, and structure types: of behaviour, political views, texts, or viruses. If desired, tangles can also be used as a new method for traditional clustering. They offer a precise, quantitative paradigm suited particularly to fuzzy clusters, since they do not require any assignment of objects to the clusters which these collectively form. This is the first of four parts of a book with the above title. The book explores applications outside mathematics of the notion and theory of tangles generalised from the graph tangles know from graph minor theory.


title: Crossing bridges with strong Szego limit theorem

abstract: We develop a new technique for computing a class of four-point correlation functions of heavy half-BPS operators in planar N=4 SYM theory which admit factorization into a product of two octagon form factors with an arbitrary bridge length. We show that the octagon can be expressed as the Fredholm determinant of the integrable Bessel operator and demonstrate that this representation is very efficient in finding the octagons both at weak and strong coupling. At weak coupling, in the limit when the four half-BPS operators become null separated in a sequential manner, the octagon obeys the Toda lattice equations and can be found in a closed form. At strong coupling, we exploit the strong Szego limit theorem to derive the leading asymptotic behavior of the octagon and, then, apply the method of differential equations to determine the remaining subleading terms of the strong coupling expansion to any order in the inverse coupling. To achieve this goal, we generalize results available in the literature for the asymptotic behavior of the determinant of the Bessel operator. As a byproduct of our analysis, we formulate a Szego-Akhiezer-Kac formula for the determinant of the Bessel operator with a Fisher-Hartwig singularity and develop a systematic approach to account for subleading power suppressed contributions.


title: Comparing Moment-Based and Monte Carlo Methods of Radiation Transport Modeling for Type II-Plateau Supernova Light Curves

abstract: Time-dependent electromagnetic signatures from core-collapse supernovae are the result of detailed transport of the shock-deposited and radioactively-powered radiation through the stellar ejecta. Due to the complexity of the underlying radiative processes, considerable approximations are made to simplify key aspects of the radiation transport problem. We present a systematic comparison of the moment-based radiation hydrodynamical code STELLA and the Monte Carlo radiation transport code Sedona in the 1D modeling of Type II-Plateau supernovae. Based on explosion models generated from the Modules for Experiments in Stellar Astrophysics (MESA) instrument, we find remarkable agreements in the modeled light curves and the ejecta structure thermal evolution, affirming the fidelity of both radiation transport modeling approaches. The radiative moments computed directly by the Monte Carlo scheme in Sedona also verify the accuracy of the moment-based scheme. We find that the coarse resolutions of the opacity tables and the numerical approximations in STELLA have insignificant impact on the resulting bolometric light curves, making it an efficient tool for the specific task of optical light curve modeling.


title: ASKAP reveals giant radio halos in two merging SPT galaxy clusters -- Making the case for a direction-dependent pipeline --

abstract: Early science observations from the Australian Square Kilometre Array Pathfinder (ASKAP) have revealed clear signals of diffuse radio emission associated with two clusters detected by the South Pole Telescope via their Sunyaev Zel'dovich signal. SPT CLJ0553-3342 (MACSJ0553.4-3342) and SPT CLJ0638-5358 (Abell S0592) are both high-mass lensing clusters that have undergone major mergers. To improve the data products of these ASKAP early science observations and create science-fidelity images of the galaxy clusters, we performed direction-dependent (DD) calibration and imaging using state-of-the-art software {\sc killMS} and {\sc DDFacet}. We find that artefacts in the ASKAP images are greatly reduced after directional calibration. Here we present our DD calibrated ASKAP radio images of both clusters showing unambiguous giant radio halos with largest linear scales of $\sim1$~Mpc. The halo in MACSJ0553.4-3342 was previously detected with GMRT observations at 323 MHz, but appears more extended in our ASKAP image. Although there is a shock detected in the thermal X-ray emission of this cluster, we find that the particle number density in the shocked region is too low to allow for the generation of a radio shock. The radio halo in Abell S0592 is a new discovery, and the Southwest border of the halo coincides with a shock detected in X-rays. We discuss the origins of these halos considering both the hadronic and turbulent re-acceleration models as well as sources of \textit{seed} electrons. This work gives a positive indication of the potential of ASKAP's Evolutionary Map of the Universe (EMU) survey in detecting intracluster medium radio sources, and showcases the improvement in data products after utilising third-generation calibration techniques.


title: Bilocal quantum criticality

abstract: We consider 2+1 dimensional conformal gauge theories coupled to additional degrees of freedom which induce a spatially local but long-range in time $1/(\tau-\tau')^2$ interaction between gauge-neutral local operators. Such theories have been argued to describe the hole-doped cuprates near optimal doping. We focus on a SU(2) gauge theory with $N_h$ flavors of adjoint Higgs fields undergoing a quantum transition between Higgs and confining phases: the $1/(\tau-\tau')^2$ interaction arises from a spectator large Fermi surface of electrons. The large $N_h$ expansion leads to an effective action containing fields which are bilocal in time but local in space. We find a strongly-coupled fixed point at order $1/N_h$, with dynamic critical exponent $z > 1$. We show that the entropy preserves hyperscaling, but nevertheless leads to a linear in temperature specific heat with a co-efficient which has a finite enhancement near the quantum critical point.


title: Gravitational path integral from the $T^2$ deformation

abstract: We study a $T^2$ deformation of large $N$ conformal field theories, a higher dimensional generalization of the $T\bar T$ deformation. The deformed partition function satisfies a flow equation of the diffusion type. We solve this equation by finding its diffusion kernel, which is given by the Euclidean gravitational path integral in $d+1$ dimensions between two boundaries with Dirichlet boundary conditions for the metric. This is natural given the connection between the flow equation and the Wheeler-DeWitt equation, on which we offer a new perspective by giving a gauge-invariant relation between the deformed partition function and the radial WDW wave function. An interesting output of the flow equation is the gravitational path integral measure which is consistent with a constrained phase space quantization. Finally, we comment on the relation between the radial wave function and the Hartle-Hawking wave functions dual to states in the CFT, and propose a way of obtaining the volume of the maximal slice from the $T^2$ deformation.


title: Empirical Bolometric Correction Coefficients for Nearby Main-Sequence Stars in Gaia Era

abstract: Nearby detached double-lined eclipsing binaries with most accurate data were studied and 290 systems were found with at least one main-sequence component having a metallicity $0.008\leq Z\leq 0.040$. Stellar parameters, light ratios, Gaia DR2 trigonometric parallaxes, extinctions and/or reddening were investigated and only 206 systems were selected eligible to calculate empirical bolometric corrections. NASA-IPAC Galactic dust maps were main source of extinctions. Unreliable extinctions at low Galactic latitudes $|b|\leq5^o$ were replaced with individual determinations, if they exist in the literature, else associated systems are discarded. Main-sequence stars of remaining systems were used to calculate bolometric corrections ($BC$) and to calibrate $BC-T_{eff}$ relation, which is valid in the range 3100-36000 K. De-reddened $(B-V)_0$ colours, on the other hand, allowed us to calibrate two intrinsic colour effective temperature relations, where the linear one is valid for $T_{eff}>10000$ K, while the quadratic relation is valid for $T_{eff}<10000$ K, that is, both are valid in the same temperature range $BC-T_{eff}$ relation is valid. New $BC$ computed from $T_{eff}$ and other astrophysical parameters are tabulated, as well.


title: Non-Hermitian Physics

abstract: A review is given on the foundations and applications of non-Hermitian classical and quantum physics. First, key theorems and central concepts in non-Hermitian linear algebra, including Jordan normal form, biorthogonality, exceptional points, pseudo-Hermiticity and parity-time symmetry, are delineated in a pedagogical and mathematically coherent manner. Building on these, we provide an overview of how diverse classical systems, ranging from photonics, mechanics, electrical circuits, acoustics to active matter, can be used to simulate non-Hermitian wave physics. In particular, we discuss rich and unique phenomena found therein, such as unidirectional invisibility, enhanced sensitivity, topological energy transfer, coherent perfect absorption, single-mode lasing, and robust biological transport. We then explain in detail how non-Hermitian operators emerge as an effective description of open quantum systems on the basis of the Feshbach projection approach and the quantum trajectory approach. We discuss their applications to physical systems relevant to a variety of fields, including atomic, molecular and optical physics, mesoscopic physics, and nuclear physics with emphasis on prominent phenomena/subjects in quantum regimes, such as quantum resonances, superradiance, continuous quantum Zeno effect, quantum critical phenomena, Dirac spectra in quantum chromodynamics, and nonunitary conformal field theories. Finally, we introduce the notion of band topology in complex spectra of non-Hermitian systems and present their classifications by providing the proof, firstly given by this review in a complete manner, as well as a number of instructive examples. Other topics related to non-Hermitian physics, including nonreciprocal transport, speed limits, nonunitary quantum walk, are also reviewed.


title: Efficiently estimating mean, uncertainty and unconstrained large scale fraction of local Universe simulations with paired fixed fields

abstract: Provided a random realization of the cosmological model, observations of our cosmic neighborhood now allow us to build simulations of the latter down to the non-linear threshold. The resulting local Universe models are thus accurate up to a given residual cosmic variance. Namely some regions and scales are apparently not constrained by the data and seem purely random. Drawing conclusions together with their uncertainties involves then statistics implying a considerable amount of computing time. By applying the constraining algorithm to paired fixed fields, this paper diverts the original techniques from their first use to efficiently disentangle and estimate uncertainties on local Universe simulations obtained with random fields. Paired fixed fields differ from random realizations in the sense that their Fourier mode amplitudes are fixed and they are exactly out of phase. Constrained paired fixed fields show that only 20% of the power spectrum on large scales (> tens of megaparsecs) is purely random. Namely 80% of it is partly constrained by the large scale / small scale data correlations. Additionally, two realizations of our local environment obtained with paired fixed fields of the same pair constitute an excellent non-biased average or quasi-linear realization of the latter, namely the equivalent of hundreds of constrained simulations. The variance between these two realizations gives the uncertainty on the achievable local Universe simulations. These two simulations will permit enhancing faster our local cosmic web understanding thanks to a drastically reduced required computational time to appreciate its modeling limits and uncertainties.


title: Effects of the Hubble Parameter on the Cosmic Growth of the First Quasars

abstract: Supermassive black holes (SMBHs) play a crucial role in the evolution of galaxies and are currently detected up to $z\sim 7.5$. Theories describing black hole (BH) growth are challenged by how rapidly seeds with initial mass $M_{\bullet} \lesssim 10^5 \, \mathrm{M_{\odot}}$, formed at $z \sim 20-30$, grew to $M_{\bullet} \sim 10^9 \, \mathrm{M_{\odot}}$ by $z\sim 7$. Here we study the effects of the value of the Hubble parameter, $H_0$, on models describing the early growth of BHs. First, we note that the predicted mass of a quasar at $z=6$ changes by $> 300$% if the underlying Hubble parameter used in the model varies from $H_0 = 65$ to $H_0 = 74$ km s$^{-1}$Mpc$^{-1}$, a range encompassing current estimates. Employing an MCMC approach based on priors from $z \gtrsim 6.5$ quasars and on $H_0$, we study the interconnection between $H_0$ and the parameters describing BH growth: seed mass $M_i$ and Eddington ratio $f_{\rm Edd}$. Assuming an Eddington ratio of $f_{\rm Edd} = 0.7$, in agreement with previous estimates, we find $H_0 = 73.6^{+1.2}_{-3.3}$ km s$^{-1}$Mpc$^{-1}$. In a second analysis, allowing all the parameters to vary freely, we find $\log(M_{i}/\mathrm{M_{\odot}}) > 4.5$ (at 95% CL), $H_0 = 74^{+1.5}_{-1.4}$ km s$^{-1}$Mpc$^{-1}$ and $f_{\rm Edd}=0.77^{+0.035}_{-0.026}$ at 68% CL. Our results on the typical Eddington ratio are in agreement with previous estimates. Current values of the Hubble parameter strongly favour heavy seed formation scenarios, with $M_i \gtrsim 10^4 \, \mathrm{M_{\odot}}$. In our model, with the priors on BH masses of quasars used, light seed formation scenarios are rejected at $\sim 3\sigma$.


title: Global symmetries, hidden sectors and emergent (dark) vector interactions

abstract: Hidden theories coupled to the SM may provide emergent (dark) vectors, that are composites/bound-states of the hidden fields. This is motivated by paradigms emerging from the AdS/CFT correspondence but it is a more general phenomenon. We explore the general setup and find that UV interactions among currents or charged fields give rise to emergent vectors in the IR. We study the general properties of such vectors and argue that they can be generically different from fundamental dark photons that have been studied so far.


title: GenetIC -- a new initial conditions generator to support genetically modified zoom simulations

abstract: We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation.