title: Highly strained mixtures of bidimensional soft and rigid grains: an experimental approach from the local scale

abstract: Granular systems are not always homogeneous and can be composed of grains with very different mechanical properties. To improve our understanding of the behavior of real granular systems, in this experimental study, we compress 2D bidisperse systems made of both soft and rigid grains. By means of a recently developed experimental set-up, \md{from the measurement of the displacement field we can} follow all the mechanical observables of this granular medium from the inside of each particle up-to-the whole system scale. We \md{are able to} detect the jamming transition from these observables and study their evolution deep in the jammed state for packing fractions as high as $0.915$. We show the uniqueness of the behavior of such a system, \md{in which way} it is similar to purely soft or rigid systems and how it is different from them. This study constitutes the first step toward a better understanding of the mechanical behavior of granular materials that are polydisperse in terms of grain rheology.


title: Simulation of adiabatic quantum computing for molecular ground states

abstract: Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy, and with only modest computational resources via an adaptive sampling configuration interaction (ASCI) scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error, and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Further, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wavefunction instead of the Hartree-Fock initial state and (ii)~a non-linear interpolation between initial and target Hamiltonians. We find that starting with a CASCI wavefunction with a limited active space yields substantial speedups for many of the systems examined while non-linear interpolation does not. Additionally, we observe interesting trends in the minimum gap location (based on the initial state) as well as how critical time can depend on certain molecular properties such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.


title: Spin-triplet superconductivity from interband effect in doped insulators

abstract: Despite being of fundamental importance and potential interest for topological quantum computing, spin-triplet superconductors remain rare in solid state materials after decades of research. In this work, we present a general mechanism for spin-triplet superconductivity in multi-band systems, where a non-retarded pairing interaction between conduction electrons is produced by their electronic repulsion to a third-electron undergoing a virtual interband transition. Our theory is analytically controlled by an interband hybridization parameter, and explicitly demonstrated in doped band insulators with the example of an extended Hubbard model. In light of this theory, we propose that recently discovered dilute superconductors such as ZrNCl and WTe$_2$ are spin-triplet, and compare the expected consequences of our theory with experimental data.


title: Correspondence between Dicke-model semiclasscial dynamics in the superradiant dipolar phase and the Euler heavy top

abstract: Analytic expression is found for the frequency dependence of transmission coefficient of a transmission line inductively coupled to the microwave cavity with superradiant condensate. Sharp transmission drops reflect condensate's frequencies spectrum. These results pave way to direct detection of emergence of the superradiant condensates in quantum metamaterials. Results are based on the analytic solutions of the nonlinear semiclassical dynamics of superradiant photonic condensate in the Dicke model of an ensemble of two-level atoms dipolar coupled to the electromagnetic field in the microwave cavity. In adiabatic limit with respect to photon degree of freedom the system is approximately integrable, with evolution being expressed via Jacobi elliptic functions of real time. Depending on the coupling strength, the semiclassical coordinate of superradiant condensate in the ground state either oscillates in one of the two degenerate minima of condensate's potential energy or traverses between them over the saddle point. An experimental setup for measuring of the breakdown of the normal phase of the Dicke model via coupling to the transmission line is proposed. A one-to-one mapping of semiclassical motion of superradiant condensate on the nodding of unstable Lagrange "sleeping top" also turns Dicke model into analogue device for modelling dynamics of mechanical systems.


title: Evolution and mass extinctions as lognormal stochastic processes

abstract: In a series of recent papers and in a book, this author put forward a mathematical model capable of embracing the search for extra-terrestrial intelligence (SETI), Darwinian Evolution and Human History into a single, unified statistical picture, concisely called Evo-SETI. The relevant mathematical tools are: (1) Geometric Brownian motion (GBM), the stochastic process representing evolution as the stochastic increase of the number of species living on Earth over the last 3.5 billion years. This GBM is well known in the mathematics of finances (Black-Sholes models). (2) The probability distributions known as b-lognormals, i.e. lognormals starting at a certain positive instant b>0 rather than at the origin. In the framework of Darwinian Evolution, the resulting mathematical construction was shown to be what evolutionary biologists call Cladistics. (3) The (Shannon) entropy of such b-lognormals is then seen to represent the 'degree of progress' reached by each living organism or by each big set of living organisms, like historic human civilizations. (4) All these results also match with SETI in that the statistical Drake equation (generalization of the ordinary Drake equation to encompass statistics) leads just to the lognormal distribution as the probability distribution for the number of extra-terrestrial civilizations existing in the Galaxy. (5) The well-known 'Molecular Clock of Evolution', namely the 'constant rate of Evolution at the molecular level' as shown by Kimura's Neutral Theory of Molecular Evolution, identifies with growth rate of the entropy of our Evo-SETI model. (6) Furthermore, we apply our Evo-SETI model to lognormal stochastic processes other than GBMs. For instance, we provide two models for the mass extinctions that occurred in the past. (7) Finally, we show that the Markov & Korotayev model for Darwinian Evolution identifies with an Evo-SETI model.


title: QUCoughScope: An Artificially Intelligent Mobile Application to Detect Asymptomatic COVID-19 Patients using Cough and Breathing Sounds

abstract: In the break of COVID-19 pandemic, mass testing has become essential to reduce the spread of the virus. Several recent studies suggest that a significant number of COVID-19 patients display no physical symptoms whatsoever. Therefore, it is unlikely that these patients will undergo COVID-19 test, which increases their chances of unintentionally spreading the virus. Currently, the primary diagnostic tool to detect COVID-19 is RT-PCR test on collected respiratory specimens from the suspected case. This requires patients to travel to a laboratory facility to be tested, thereby potentially infecting others along the way.It is evident from recent researches that asymptomatic COVID-19 patients cough and breath in a different way than the healthy people. Several research groups have created mobile and web-platform for crowdsourcing the symptoms, cough and breathing sounds from healthy, COVID-19 and Non-COVID patients. Some of these data repositories were made public. We have received such a repository from Cambridge University team under data-sharing agreement, where we have cough and breathing sound samples for 582 and 141 healthy and COVID-19 patients, respectively. 87 COVID-19 patients were asymptomatic, while rest of them have cough. We have developed an Android application to automatically screen COVID-19 from the comfort of people homes. Test subjects can simply download a mobile application, enter their symptoms, record an audio clip of their cough and breath, and upload the data anonymously to our servers. Our backend server converts the audio clip to spectrogram and then apply our state-of-the-art machine learning model to classify between cough sounds produced by COVID-19 patients, as opposed to healthy subjects or those with other respiratory conditions. The system can detect asymptomatic COVID-19 patients with a sensitivity more than 91%.


title: Grain Boundary Segregation Transitions and Critical Phenomena in Binary Regular Solutions: A Systematics of Complexion Diagrams with Universal Characters

abstract: A systematics of grain boundary (GB) segregation transitions and critical phenomena has been derived to expand the classical GB segregation theory. Using twist GBs as an example, this study uncovers when GB layering vs. prewetting transitions should occur and how they are related to one another. Moreover, a novel descriptor, normalized segregation strength (phi_seg), is introduced. It can represent several factors that control GB segregation, including strain and bond energies, as well as misorientation for small-angle GBs (in a mean-field approximation), which had to be treated separately in prior models. In a strong segregation system with a large phi_seg, first-order layering transitions occur at low temperatures and become continuous above GB roughing temperatures. With reducing phi_seg, the layering transitions gradually merge and finally lump into prewetting transitions without quantized layer numbers, akin to Cahn's critical-point wetting model. Furthermore, GB complexion diagrams with universal characters are constructed as the GB counterpart to the classical exemplar of Pelton-Thompson regular-solution binary bulk phase diagrams.


title: The Concept of an Autonomic Avionics Platform and the Resulting Software Engineering Challenges

abstract: The self-* properties commonly associated with the concept of autonomic computing are capabilities desirable for avionics software platforms. They decrease the configuration effort and inherently provide new fault tolerance and resource savings possibilities. The rigid certification process and the requirements for a static and predetermined system behavior are, however, in contradiction with the adaptive and flexible nature of autonomic computing systems. We propose a partition-based architecture providing autonomic features for avionics software platforms while being compliant to regulations and accepted technologies, such as ARINC 653. The core is a platform consciousness based on a domain-specific model and a novel MAP-QE-K cycle. Moreover, we suggest a planning intelligence, a virtual qualification authority, and a minimized execution unit. For each component we define the required design assurance level and possible realization techniques. We discuss the overall feasibility and point out central challenges in the fields of runtime verification and models at runtime. These challenges need to be solved up to the realization of autonomic avionics, e.g. a virtual security assessment and a qualifiable domain-specific model database.


title: Machine learning based in situ quality estimation by molten pool condition-quality relations modeling using experimental data

abstract: The advancement of machine learning promises the ability to accelerate the adoption of new processes and property designs for metal additive manufacturing. The molten pool geometry and molten pool temperature are the significant indicators for the final part's geometric shape and microstructural properties for the Wire-feed laser direct energy deposition process. Thus, the molten pool condition-property relations are of preliminary importance for in situ quality assurance. To enable in situ quality monitoring of bead geometry and characterization properties, we need to continuously monitor the sensor's data for molten pool dimensions and temperature for the Wire-feed laser additive manufacturing (WLAM) system. We first develop a machine learning convolutional neural network (CNN) model for establishing the correlations from the measurable molten pool image and temperature data directly to the geometric shape and microstructural properties. The multi-modality network receives both the camera image and temperature measurement as inputs, yielding the corresponding characterization properties of the final build part (e.g., fusion zone depth, alpha lath thickness). The performance of the CNN model is compared with the regression model as a baseline. The developed models enable molten pool condition-quality relations mapping for building quantitative and collaborative in situ quality estimation and assurance framework.


title: Understanding performance variability in standard and pipelined parallel Krylov solvers

abstract: In this work, we collect data from runs of Krylov subspace methods and pipelined Krylov algorithms in an effort to understand and model the impact of machine noise and other sources of variability on performance. We find large variability of Krylov iterations between compute nodes for standard methods that is reduced in pipelined algorithms, directly supporting conjecture, as well as large variation between statistical distributions of runtimes across iterations. Based on these results, we improve upon a previously introduced nondeterministic performance model by allowing iterations to fluctuate over time. We present our data from runs of various Krylov algorithms across multiple platforms as well as our updated non-stationary model that provides good agreement with observations. We also suggest how it can be used as a predictive tool.


title: Transfer Learning with Ensembles of Deep Neural Networks for Skin Cancer Detection in Imbalanced Data Sets

abstract: Several machine learning techniques for accurate detection of skin cancer from medical images have been reported. Many of these techniques are based on pre-trained convolutional neural networks (CNNs), which enable training the models based on limited amounts of training data. However, the classification accuracy of these models still tends to be severely limited by the scarcity of representative images from malignant tumours. We propose a novel ensemble-based CNN architecture where multiple CNN models, some of which are pre-trained and some are trained only on the data at hand, along with auxiliary data in the form of metadata associated with the input images, are combined using a meta-learner. The proposed approach improves the model's ability to handle limited and imbalanced data. We demonstrate the benefits of the proposed technique using a dataset with 33126 dermoscopic images from 2056 patients. We evaluate the performance of the proposed technique in terms of the F1-measure, area under the ROC curve (AUC-ROC), and area under the PR-curve (AUC-PR), and compare it with that of seven different benchmark methods, including two recent CNN-based techniques. The proposed technique compares favourably in terms of all the evaluation metrics.


title: Exemplars can Reciprocate Principal Components

abstract: This paper presents a clustering algorithm that is an extension of the Category Trees algorithm. Category Trees is a clustering method that creates tree structures that branch on category type and not feature. The development in this paper is to consider a secondary order of clustering that is not the category to which the data row belongs, but the tree, representing a single classifier, that it is eventually clustered with. Each tree branches to store subsets of other categories, but the rows in those subsets may also be related. This paper is therefore concerned with looking at that second level of clustering between the other category subsets, to try to determine if there is any consistency over it. It is argued that Principal Components may be a related and reciprocal type of structure, and there is an even bigger question about the relation between exemplars and principal components, in general. The theory is demonstrated using the Portugal Forest Fires dataset as a case study. The Category Trees are then combined with other Self-Organising algorithms from the author and it is suggested that they all belong to the same family type, which is an Entropy-style of classifier.


title: Learning to Robustly Negotiate Bi-Directional Lane Usage in High-Conflict Driving Scenarios

abstract: Recently, autonomous driving has made substantial progress in addressing the most common traffic scenarios like intersection navigation and lane changing. However, most of these successes have been limited to scenarios with well-defined traffic rules and require minimal negotiation with other vehicles. In this paper, we introduce a previously unconsidered, yet everyday, high-conflict driving scenario requiring negotiations between agents of equal rights and priorities. There exists no centralized control structure and we do not allow communications. Therefore, it is unknown if other drivers are willing to cooperate, and if so to what extent. We train policies to robustly negotiate with opposing vehicles of an unobservable degree of cooperativeness using multi-agent reinforcement learning (MARL). We propose Discrete Asymmetric Soft Actor-Critic (DASAC), a maximum-entropy off-policy MARL algorithm allowing for centralized training with decentralized execution. We show that using DASAC we are able to successfully negotiate and traverse the scenario considered over 99% of the time. Our agents are robust to an unknown timing of opponent decisions, an unobservable degree of cooperativeness of the opposing vehicle, and previously unencountered policies. Furthermore, they learn to exhibit human-like behaviors such as defensive driving, anticipating solution options and interpreting the behavior of other agents.


title: Selective information exchange in collaborative clustering using regularized Optimal Transport

abstract: Collaborative learning has recently achieved very significant results. It still suffers, however, from several issues, including the type of information that needs to be exchanged, the criteria for stopping and how to choose the right collaborators. We aim in this paper to improve the quality of the collaboration and to resolve these issues via a novel approach inspired by Optimal Transport theory. More specifically, the objective function for the exchange of information is based on the Wasserstein distance, with a bidirectional transport of information between collaborators. This formulation allows to learns a stopping criterion and provide a criterion to choose the best collaborators. Extensive experiments are conducted on multiple data-sets to evaluate the proposed approach.


title: IMF-induced intrinsic uncertainties on measuring galaxy distances based on the number of giant stars: the case of the ultra-diffuse galaxy NGC 1052-DF2

abstract: The surface brightness fluctuation (SBF) technique is one of the distance measurement methods that has been applied on the low surface brightness galaxy NGC 1052-DF2 yielding a distance of about 20 Mpc implying it to be a dark matter deficient galaxy. We assume the number of giant stars above a given luminosity threshold to represent the SBF magnitude. The SBF magnitude depends on the distance, but this is degenerate with the star formation history (SFH). Using a stellar population synthesis model we calculate the number of giant stars for stellar populations with different galaxy-wide stellar initial mass functions (gwIMFs), ages, metallicities and SFHs. If the gwIMF is the invariant canonical IMF, the 1$\sigma$ (3$\sigma$) uncertainty in colour allows a distance as low as 12 Mpc (8 Mpc). If instead the true underlying gwIMF is the integrated galaxy-wide IMF (IGIMF) then overestimating distances for low-mass galaxies would be a natural result, allowing NGC 1052-DF2 to have a distance of 11 Mpc within the 1$\sigma$ colour uncertainty. Finally, we show that our main conclusion on the existence of a bias in the SBF distance estimation is not much affected by changing the luminosity lower limit for counting giant stars.


title: The SIBELIUS Project: E Pluribus Unum

abstract: We introduce "Simulations Beyond The Local Universe" (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provides an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group characteristics emanating from different scales: the total mass, orientation, orbital energy and the angular momentum are largely determined by modes above $\lambda$ = 1.6 comoving Mpc (cMpc) in the primordial density field. Smaller scale variations are mostly manifest as perturbations to the MW-M31 orbit, and we find that the observables commonly used to describe the Local Group -- the MW-M31 separation and radial velocity -- are transient and depend on specifying scales down to 0.2 cMpc in the primordial density field. We further find that the presence of M33/LMC analogues significantly affects the MW-M31 orbit and its sensitivity to small-scale perturbations. We construct initial conditions that lead to the formation of a Local Group whose primary observables precisely match the current observations.


title: $A_{FB}$ in the SMEFT: precision Z physics at the LHC

abstract: We study the forward-backward asymmetry $A_{FB}$ in $pp \to \ell^+\ell^-$ at the Z peak within the Standard Model Effective Field Theory (SMEFT). We find that this observable provides per mille level constraints on the vertex corrections of the Z boson to quarks,which close a flat direction in the electroweak precision SMEFT fit. Moreover, we show that current $A_{FB}$ data is precise enough so that its inclusion in the fit improves significantly LEP bounds even in simple New Physics setups. This demonstrates that the LHC can compete with and complement LEP when it comes to precision measurements of the Z boson properties


title: Fuzzball Shadows: Emergent Horizons from Microstructure

abstract: We study the physical properties of four-dimensional, string-theoretical, horizonless "fuzzball" geometries by imaging their shadows. Their microstructure traps light rays straying near the would-be horizon on long-lived, highly redshifted chaotic orbits. In fuzzballs sufficiently near the scaling limit this creates a shadow much like that of a black hole, while avoiding the paradoxes associated with an event horizon. Observations of the shadow size and residual glow can potentially discriminate between fuzzballs away from the scaling limit and alternative models of black compact objects.


title: SOAR TESS Survey. II: The impact of stellar companions on planetary populations

abstract: We present the results of the second year of exoplanet candidate host speckle observations from the SOAR TESS survey. We find 89 of the 589 newly observed TESS planet candidate hosts have companions within 3\arcsec, resulting in light curve dilution, that if not accounted for leads to underestimated planetary radii. We combined these observations with those from paper I to search for evidence of the impact binary stars have on planetary systems. Removing the quarter of the targets observed identified as false-positive planet detections, we find that transiting planet are suppressed by nearly a factor-of-seven in close solar-type binaries, nearly twice the suppression previously reported. The result on planet occurrence rates that are based on magnitude limited surveys is an overestimation by a factor of two if binary suppression is not taken into account. We also find tentative evidence for similar close binary suppression of planets in M-dwarf systems. Lastly, we find that the high rates of widely separated companions to hot Jupiter hosts previously reported was likely a result of false-positive contamination in our sample.


title: $W^+W^-$ production at NNLO+PS with MiNNLO$_{\rm PS}$

abstract: We consider $W^+W^-$ production in hadronic collisions and present the computation of next-to-next-to-leading order accurate predictions consistently matched to parton showers (NNLO+PS) using the MiNNLO$_{\rm PS}$ method. Spin correlations, interferences and off-shell effects are included by calculating the full process $pp \to e^+\nu_e \mu^-\bar{\nu}_\mu$. This is the first NNLO+PS calculation for $W^+W^-$ production that does not require an a-posteriori multi-differential reweighting. The evaluation time of the two-loop contribution has been reduced by more than one order of magnitude through a four-dimensional cubic spline interpolation. We find good agreement with the inclusive and fiducial cross sections measured by ATLAS and CMS. Both NNLO corrections and matching to parton showers are important for an accurate simulation of the $W^+W^-$ signal, and their matching provides the best description of fully exclusive $W^+W^-$ events to date.


title: Gravity as an ensemble and the moment problem

abstract: If a bulk gravitational path integral can be identified with an average of partition functions over an ensemble of boundary quantum theories, then a corresponding moment problem can be solved. We review existence and uniqueness criteria for the Stieltjes moment problem, which include an infinite set of positivity conditions. The existence criteria are useful to rule out an ensemble interpretation of a theory of gravity, or to indicate incompleteness of the gravitational data. We illustrate this in a particular class of 2D gravities including variants of the CGHS model and JT supergravity. The uniqueness criterium is relevant for an unambiguous determination of quantities such as $\overline{\log Z(\beta)}$ or the quenched free energy. We prove in JT gravity that perturbation theory, both in the coupling which suppresses higher-genus surfaces and in the temperature, fails when the number of boundaries is taken to infinity. Since this asymptotic data is necessary for the uniqueness problem, the question cannot be settled without a nonperturbative completion of the theory.


title: Dynamically Induced Exceptional Phases in Quenched Interacting Semimetals

abstract: We report on the dynamical formation of exceptional degeneracies in basic correlation functions of non-integrable one- and two-dimensional systems quenched to the vicinity of a critical point. Remarkably, fine-tuned semi-metallic points in the phase diagram of the considered systems are thereby promoted to topologically robust non-Hermitian (NH) nodal phases emerging in the coherent long-time evolution of a dynamically equilibrating system. In the framework of non-equilibrium Green's function methods within the conserving second Born approximation, we predict observable signatures of these novel NH nodal phases in simple spectral functions as well as in the time-evolution of momentum distribution functions.

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


title: The 800pc long tidal tails of the Hyades star cluster: Possible discovery of candidate epicyclic overdensities from an open star cluster

abstract: The tidal tails of stellar clusters provide an important tool for studying the birth conditions of the clusters and their evolution, coupling, and interaction with the Galactic potential. We present the N-body evolution of a Hyades-like stellar cluster with backward-integrated initial conditions on a realistic 3D orbit in the Milky Way computed within the AMUSE framework. For the first time, we explore the effect of the initial cluster rotation and the presence of lumps in the Galactic potential on the formation and evolution of tidal tails. We show that the tidal tails are not naturally clustered in any coordinate system. Models with initial rotation result in significant differences in the cluster mass loss and follow different angular momentum time evolution. The orientation of the tidal tails relative to the motion vector of the cluster and the current cluster angular momentum constrain the initial rotation of the cluster. We highlight the use of the convergent point (CP) method in searches for co-moving groups and introduce a new compact CP (CCP) method that accounts for internal kinematics based on an assumed model. Using the CCP method, we are able to recover candidate members of the Hyades tidal tails in the Gaia DR2 and eDR3 reaching a total extent of almost 1kpc. We confirm the previously noted asymmetry in the detected tidal tails. In the eDR3 data we recovered spatial overdensities in the leading and trailing tails that are kinematically consistent with being epicyclic overdensities and thus present candidates for the first such detection in an open star cluster. We show that the epicyclic overdensities are able to provide constraints not only on the cluster properties, but also on the Galactic potential. Finally, based on N-body simulations, a close encounter with a massive Galactic lump can explain the observed asymmetry in the tidal tails of the Hyades.(abriged)


title: Radiative Driving of the AGN Outflows in the Narrow-Line Seyfert 1 Galaxy NGC 4051

abstract: We explore the properties of ionized gas in the nuclear and circumnuclear environment of the narrow-line Seyfert 1 galaxy NGC 4051 using spectroscopic and imaging observations from the Hubble Space Telescope (HST) and Apache Point Observatory (APO)'s ARC 3.5m Telescope. We identify an unresolved moderate-density intermediate width component and a high-density broad component in the optical emission lines from the active nucleus, as well as spatially-resolved emission extending up to $\sim$1 kpc in the AGN ionized narrow-line region (NLR) and $\sim$8 kpc in the stellar ionized host galaxy. The HST narrow-band image reveals a distinct conical structure in [O III] emission towards the NE, and the ionized gas kinematics shows up to two blueshifted velocity components, indicating outflows along the edges of a cone. We introduce an improved model of biconical outflow, with our line of sight passing through the wall of the cone, which suggests that the large number of outflowing UV absorbers seen in NGC 4051 are NLR clouds in absorption. Using the de-projection factors from the biconical geometry, we measure true outflow velocities up to 680 km s$^{-1}$ at a distance of $\sim$350 pc, however, we do not find any rotational signature inside a projected distance $\leq$ 10'' ($\sim$800 pc) from the nucleus. We compare the gas kinematics with analytical models based on a radiation-gravity formalism, which show that most of the observed NLR outflows are launched within $\sim$0.5 pc of the nucleus and can travel up to $\sim$1 kpc from this low-luminosity AGN.


title: Dipole Radiation and Beyond from Axion Stars in Electromagnetic Fields

abstract: We investigate the production of photons from coherently oscillating, spatially localized clumps of axionic fields (oscillons and axion stars) in the presence of external electromagnetic fields. We delineate different qualitative behaviour of the photon luminosity in terms of an effective dimensionless coupling parameter constructed out of the axion-photon coupling, and field amplitude, oscillation frequency and radius of the axion star. For small values of this dimensionless coupling, we provide a general analytic formula for the dipole radiation field and the photon luminosity per solid angle, including a strong dependence on the radius of the configuration. For moderate to large coupling, we report on a non-monotonic behavior of the luminosity with the coupling strength in the presence of external magnetic fields. After an initial rise in luminosity with the coupling strength, we see a suppression (by an order of magnitude or more compared to the dipole radiation approximation) at moderately large coupling. At sufficiently large coupling, we find a transition to a regime of exponential growth of the luminosity due to parametric resonance. We carry out 3+1 dimensional lattice simulations of axion electrodynamics, at small and large coupling, including non-perturbative effects of parametric resonance as well as backreaction effects when necessary. We also discuss medium (plasma) effects that lead to resonant axion to photon conversion, relevance of the coherence of the soliton, and implications of our results in astrophysical and cosmological settings.


title: Large Pauli Limit Violation and Reentrant Superconductivity in Magic-Angle Twisted Trilayer Graphene

abstract: Moir\'e quantum matter has emerged as a novel materials platform where correlated and topological phases can be explored with unprecedented control. Among them, magic-angle systems constructed from two or three layers of graphene have shown robust superconducting phases with unconventional characteristics. However, direct evidence for unconventional pairing remains to be experimentally demonstrated. Here, we show that magic-angle twisted trilayer graphene (MATTG) exhibits superconductivity up to in-plane magnetic fields in excess of 10 T, which represents a large ($2\sim3$ times) violation of the Pauli limit for conventional spin-singlet superconductors. This observation is surprising for a system which is not expected to have strong spin-orbit coupling. Furthermore, the Pauli limit violation is observed over the entire superconducting phase, indicating that it is not related to a possible pseudogap phase with large superconducting amplitude pairing. More strikingly, we observe reentrant superconductivity at large magnetic fields, which is present in a narrower range of carrier density and displacement field. These findings suggest that the superconductivity in MATTG is likely driven by a mechanism resulting in non-spin-singlet Cooper pairs, where the external magnetic field can cause transitions between phases with potentially different order parameters. Our results showcase the richness of moir\'e superconductivity and may pave a new route towards designing next-generation exotic quantum matter.


title: Dynamical dust traps in misaligned circumbinary discs: analytical theory and numerical simulations

abstract: Recent observations have shown that circumbinary discs can be misaligned with respect to the binary orbital plane.The lack of spherical symmetry, together with the non-planar geometry of these systems, causes differential precession which might induce the propagation of warps. While gas dynamics in such environments is well understood, little is known about dusty discs. In this work, we analytically study the problem of dust traps formation in misaligned circumbinary discs. We find that pile-ups may be induced not by pressure maxima, as the usual dust traps, but by a difference in precession rates between the gas and dust. Indeed, this difference makes the radial drift inefficient in two locations, leading to the formation of two dust rings whose position depends on the system parameters. This phenomenon is likely to occur to marginally coupled dust particles $(\text{St}\gtrsim1)$ as both the effect of gravitational and drag force are considerable. We then perform a suite of three-dimensional SPH numerical simulations to compare the results with our theoretical predictions. We explore the parameter space, varying stellar mass ratio, disc thickness, radial extension, and we find a general agreement with the analytical expectations. Such dust pile-up prevents radial drift, fosters dust growth and may thus promote the planet formation in circumbinary discs.


title: Axion Haloscope Array With $\mathcal{PT}$ Symmetry

abstract: We generalize the recently proposed $\mathcal{P T}$-symmetric axion haloscope to a larger array with more $\mathcal{P T}$-symmetric structures. By broadening the response bandwidth of the signal without increasing the readout noise, the optimized scan rate of the axion haloscope is significantly enhanced, as well as is the signal power. Furthermore, we show that the robustness of the detector towards the variations of the array coupling is the strongest when a binary tree structure is introduced which contains a largely enhanced $\mathcal{P T}$ symmetry. The multiple allowed probing sensors can further increase the scan rate by a factor of the sensors' number due to the correlation of the signals. This type of array can strongly boost the search for an axion compared to single-mode resonant detection. The enhancement to the scan rate becomes the most manifest when applied to the proposed detection using a superconducting radio-frequency cavity with an ac magnetic field where most of the parameter space of the QCD axion above kHz can be probed.


title: INSPIRE: INvestigating Stellar Population In RElics II. First Data Release (DR1)

abstract: The INvestigating Stellar Population In RElics is an on-going project targeting 52 ultra-compact massive galaxies at 0.1<z<0.5 with the X-Shooter@VLT spectrograph (XSH). These objects are the perfect candidates to be 'relics', massive red-nuggets formed at high-z (z>2) through a short and intense star formation burst, that evolved passively and undisturbed until the present-day. Relics provide a unique opportunity to study the mechanisms of star formation at high-z. In this paper, we present the first INSPIRE Data Release, comprising 19 systems with observations completed in 2020. We use the methods already presented in the INSPIRE Pilot, but revisiting the 1D spectral extraction. For these 19 systems, we obtain an estimate of the stellar velocity dispersion, fitting separately the two UVB and VIS XSH arms at their original resolution. We estimate [Mg/Fe] abundances via line-index strength and mass-weighted integrated stellar ages and metallicities with full spectral fitting on the combined spectrum. Ages are generally old, in agreement with the photometric ones, and metallicities are almost always super-solar, confirming the mass-metallicity relation. The [Mg/Fe] ratio is also larger than solar for the great majority of the galaxies, as expected. We find that 10 objects have formed more than 75% of their stellar mass (M*) within 3 Gyr from the Big Bang and classify them as relics. Among these, we identify 4 galaxies which had already fully assembled their M* by that time. They are therefore `extreme relics' of the ancient Universe. The INSPIRE DR1 catalogue of 10 known relics to-date augment by a factor of 3.3 the total number of confirmed relics, also enlarging the redshift window. It is therefore the largest publicly available collection. Thanks to the larger number of systems, we can also better quantify the existence of a 'degree of relicness', already hinted at the Pilot Paper.


title: A brief review on primordial black holes as dark matter

abstract: Primordial black holes (PBHs) represent a natural candidate for one of the components of the dark matter (DM) in the Universe. In this review, we shall discuss the basics of their formation, abundance and signatures. Some of their characteristic signals are examined, such as the emission of particles due to Hawking evaporation and the accretion of the surrounding matter, effects which could leave an impact in the evolution of the Universe and the formation of structures. The most relevant probes capable of constraining their masses and population are discussed.


title: The hydrodynamic evolution of binary black holes embedded within the vertically stratified disks of active galactic nuclei

abstract: Stellar-mass black holes can become embedded within the gaseous disks of active galactic nuclei (AGNs). Afterwards, their interactions are mediated by their gaseous surroundings. In this work, we study the evolution of stellar-mass binary black holes (BBHs) embedded within AGN disks using a combination of three-dimensional hydrodynamic simulations and analytic methods, focusing on environments in which the AGN disk scale height $H$ is $\gtrsim$ the BBH sphere of influence. We model the local surroundings of the embedded BBHs using a wind tunnel formalism and characterize different accretion regimes based on the local properties of the disk, which range from wind-dominated to quasi-spherical. We use our simulations to develop prescriptions for mass accretion and drag for embedded BBHs. We use these prescriptions, along with AGN disk models that can represent the Toomre-unstable outer regions of AGN disks, to study the long-term evolution of the BBHs as they migrate through the disk. We find that BBHs typically merge within $\lesssim 5-30\,{\rm Myr}$, increasing their mass significantly in the process, allowing BBHs to enter (or cross) the pair-instability supernova mass gap. The rate at which gas is supplied to these BBHs often exceeds the Eddington limit, sometimes by several orders of magnitude. We conclude that most embedded BBHs will merge before migrating significantly in the disk. Depending on the conditions of the ambient gas and the distance to the system, LISA can detect the transition between the gas-dominated and gravitational wave dominated regime for inspiraling BBHs that are formed sufficiently close to the AGN ($\lesssim$ 0.1 pc). We also discuss possible electromagnetic signatures during and following the inspiral, finding that it is generally unlikely but not inconceivable for the bolometric luminosity of the BBH to exceed that of the host AGN.


title: Multi-TeV Signals of Baryogenesis in Higgs Troika Model

abstract: A modest extension of the Standard Model by two additional Higgs doublets - the Higgs Troika Model - can provide a well-motivated scenario for successful baryogenesis if neutrinos are Dirac fermions. Adapting the "Spontaneous Flavor Violation" framework, we consider a version of the Troika model where light quarks have significant couplings to the new multi-TeV Higgs states. Resonant production of new scalars leading to di-jet or top-pair signals are typical predictions of this setup. The initial and final state quarks relevant to the collider phenomenology also play a key role in baryogenesis, potentially providing direct access to the relevant early Universe physics in high energy experiments. Viable baryogenesis generally prefers some hierarchy of masses between the observed and the postulated Higgs states. We show that there is a complementarity between direct searches at a future 100 TeV $pp$ collider and indirect searches at flavor experiments, with both sensitive to different regions of parameter space relevant for baryogenesis. In particular, measurements of $D-\bar{D}$ mixing at LHCb probe much of the interesting parameter space. Direct and indirect searches can uncover the new Higgs states up to masses of $\mathcal{O}(10)$ TeV, thereby providing an impressive reach to investigate this model.


title: The Balloon-Borne Cryogenic Telescope Testbed Mission: Bulk Cryogen Transfer at 40 km Altitude

abstract: The Balloon-Borne Cryogenic Telescope Testbed (BOBCAT) is a stratospheric balloon payload to develop technology for a future cryogenic suborbital observatory. A series of flights are intended to establish ultra-light dewar performance and open-aperture observing techniques for large (3--5 meter diameter) cryogenic telescopes at infrared wavelengths. An initial flight in 2019 demonstrated bulk transfer of liquid nitrogen and liquid helium at stratospheric altitudes. An 827 kg payload carried 14 liters of liquid nitrogen (LN2) and 268 liters of liquid helium (LHe) in pressurized storage dewars to an altitude of 39.7 km. Once at float altitude, liquid nitrogen transfer cooled a separate, unpressurized bucket dewar to a temperature of 65 K, followed by the transfer of 32 liters of liquid helium from the storage dewar into the bucket dewar. Calorimetric tests measured the total heat leak to the LHe bath within bucket dewar. A subsequent flight will replace the receiving bucket dewar with an ultra-light dewar of similar size to compare the performance of the ultra-light design to conventional superinsulated dewars.


title: Transformer-Based Attention Networks for Continuous Pixel-Wise Prediction

abstract: While convolutional neural networks have shown a tremendous impact on various computer vision tasks, they generally demonstrate limitations in explicitly modeling long-range dependencies due to the intrinsic locality of the convolution operation. Initially designed for natural language processing tasks, Transformers have emerged as alternative architectures with innate global self-attention mechanisms to capture long-range dependencies. In this paper, we propose TransDepth, an architecture that benefits from both convolutional neural networks and transformers. To avoid the network losing its ability to capture local-level details due to the adoption of transformers, we propose a novel decoder that employs attention mechanisms based on gates. Notably, this is the first paper that applies transformers to pixel-wise prediction problems involving continuous labels (i.e., monocular depth prediction and surface normal estimation). Extensive experiments demonstrate that the proposed TransDepth achieves state-of-the-art performance on three challenging datasets. Our code is available at: https://github.com/ygjwd12345/TransDepth.


title: Recursion relations for 5-point conformal blocks

abstract: We consider 5-point functions in conformal field theories in d > 2 dimensions. Using weight-shifting operators, we derive recursion relations which allow for the computation of arbitrary conformal blocks appearing in 5-point functions of scalar operators, reducing them to a linear combination of blocks with scalars exchanged. We additionally derive recursion relations for the conformal blocks which appear when one of the external operators in the 5-point function has spin 1 or 2. Our results allow us to formulate positivity constraints using 5-point functions which describe the expectation value of the energy operator in bilocal states created by two scalars.


title: Contraction criteria for Brownian filtrations samplings

abstract: This paper investigates the problem to determine whether a given stochastic process generates a sampled Brownian filtration. A fairly general sufficient condition is obtained by applying the Frank H. Clarke contraction criteria to a functional whose construction relies on the quasi-invariance properties of a sampled Wiener measure. The latter are investigated from the precise analytic structure underlying this sampled Brownian motion. In particular, the Cameron-Martin space is identified from usual Lebesgue integrals over a Guseinov measure. As an application we obtain sufficient conditions for the existence of a solution to a class of not necessarily Markovian stochastic differential equations driven by a sampled Brownian motion. By taking a sampling times set which coincides with the unit interval, the results apply to continuous time models.


title: Modelling intransitivity in pairwise comparisons with application to baseball data

abstract: The seminal Bradley-Terry model exhibits transitivity, i.e., the property that the probabilities of player A beating B and B beating C give the probability of A beating C, with these probabilities determined by a skill parameter for each player. Such transitive models do not account for different strategies of play between each pair of players, which gives rise to {\it intransitivity}. Various intransitive parametric models have been proposed but they lack the flexibility to cover the different strategies across $n$ players, with the $O(n^2)$ values of intransitivity modelled using O(n) parameters, whilst they are not parsimonious when the intransitivity is simple. We overcome their lack of adaptability by allocating each pair of players to one of a random number of $K$ intransitivity levels, each level representing a different strategy. Our novel approach for the skill parameters involves having the $n$ players allocated to a random number of $A<n$ distinct skill levels, to improve efficiency and avoid false rankings. Although we may have to estimate up to $O(n^2)$ unknown parameters for $(A,K)$ we anticipate that in many practical contexts $A+K < n$. Using a Bayesian hierarchical model, $(A,K)$ are treated as unknown, and inference is conducted via a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm. Our semi-parametric model, which gives the Bradley-Terry model when $(A=n-1, K=0)$, is shown to have an improved fit relative to the Bradley-Terry, and the existing intransitivity models, in out-of-sample testing when applied to simulated and American League baseball data. Supplementary materials for the article are available online.


title: Am I fit for this physical activity? Neural embedding of physical conditioning from inertial sensors

abstract: Inertial Measurement Unit (IMU) sensors are present in everyday devices such as smartphones and fitness watches. As a result, the array of health-related research and applications that tap onto this data has been growing, but little attention has been devoted to the prediction of an individual's heart rate (HR) from IMU data, when undergoing a physical activity. Would that be even possible? If so, this could be used to design personalized sets of aerobic exercises, for instance. In this work, we show that it is viable to obtain accurate HR predictions from IMU data using Recurrent Neural Networks, provided only access to HR and IMU data from a short-lived, previously executed activity. We propose a novel method for initializing an RNN's hidden state vectors, using a specialized network that attempts to extract an embedding of the physical conditioning (PCE) of a subject. We show that using a discriminator in the training phase to help the model learn whether two PCEs belong to the same individual further reduces the prediction error. We evaluate the proposed model when predicting the HR of 23 subjects performing a variety of physical activities from IMU data available in public datasets (PAMAP2, PPG-DaLiA). For comparison, we use as baselines the only model specifically proposed for this task and an adapted state-of-the-art model for Human Activity Recognition (HAR), a closely related task. Our method, PCE-LSTM, yields over 10% lower mean absolute error. We demonstrate empirically that this error reduction is in part due to the use of the PCE. Last, we use the two datasets (PPG-DaLiA, WESAD) to show that PCE-LSTM can also be successfully applied when photoplethysmography (PPG) sensors are available, outperforming the state-of-the-art deep learning baselines by more than 30%.


title: Back to sources -- the role of losses and coherence in super-resolution imaging revisited

abstract: Photon losses are intrinsic for any translationally invariant optical imaging system with a non-trivial Point Spread Function, and the relation between the transmission factor and the coherence properties of an imaged object is universal -- we demonstrate the rigorous proof of this statement, based on the principles of quantum mechanics. The fundamental limit on the precision of estimating separation between two partially coherent sources is then derived. The careful study of the role of photon losses allows to resolve conflicting claims present in previous works. We compute the Quantum Fisher Information for the generic model of optical 4f imaging system, and use prior considerations to validate the result for a general, translationally invariant imaging apparatus. We prove that the spatial-mode demultiplexing (SPADE) measurement, optimal for non-coherent sources, remains optimal for an arbitrary degree of coherence. Moreover, we show that some approximations, omnipresent in theoretical works about optical imaging, inevitably lead to unphysical, zero-transmission models, resulting in misleading claims regarding fundamental resolution limits.


title: The Hubbard Model

abstract: The repulsive Hubbard model has been immensely useful in understanding strongly correlated electron systems, and serves as the paradigmatic model of the field. Despite its simplicity, it exhibits a strikingly rich phenomenology which is reminiscent of that observed in quantum materials. Nevertheless, much of its phase diagram remains controversial. Here, we review a subset of what is known about the Hubbard model, based on exact results or controlled approximate solutions in various limits, for which there is a suitable small parameter. Our primary focus is on the ground state properties of the system on various lattices in two spatial dimensions, although both lower and higher dimensions are discussed as well. Finally, we highlight some of the important outstanding open questions.


title: The Structural Evolution of Isolated Galaxies at Low Redshift in the IllustrisTNG Simulation

abstract: We study the structural evolution of isolated star-forming galaxies in the Illustris TNG100-1 hydrodynamical simulation, with a focus on investigating the growth of the central core density within 2 kpc ($\Sigma_{*,2kpc}$) in relation to total stellar mass ($M_*$) at z < 0.5. First, we show that several observational trends in the $\Sigma_{*,2kpc}$-$M_*$ plane are qualitatively reproduced in IllustrisTNG, including the distributions of AGN, star forming galaxies, quiescent galaxies, and radial profiles of stellar age, sSFR, and metallicity. We find that galaxies with dense cores evolve parallel to the $\Sigma_{*,2kpc}$-$M_*$ relation, while galaxies with diffuse cores evolve along shallower trajectories. We investigate possible drivers of rapid growth in $\Sigma_{*,2kpc}$ compared to $M_*$. Both the current sSFR gradient and the BH accretion rate are indicators of past core growth, but are not predictors of future core growth. Major mergers (although rare in our sample; $\sim$10%) cause steeper core growth, except for high mass ($M_*$ >$\sim$ $10^{10} M_{\odot}$) mergers, which are mostly dry. Disc instabilities, as measured by the fraction of mass with Toomre Q < 2, are not predictive of rapid core growth. Instead, rapid core growth results in more stable discs. The cumulative black hole feedback history sets the maximum rate of core growth, preventing rapid growth in high-mass galaxies ($M_*$ >$\sim$ $10^{9.5} M_{\odot}$). For massive galaxies the total specific angular momentum of accreting gas is the most important predictor of future core growth. Our results suggest that the angular momentum of accreting gas controls the slope, width and zero-point evolution of the $\Sigma_{*,2kpc}$-$M_*$ relation.


title: Detecting the radiative decay of the cosmic neutrino background with line-intensity mapping

abstract: We study the possibility to use line-intensity mapping (LIM) to seek photons from the radiative decay of neutrinos in the cosmic neutrino background. The Standard Model prediction for the rate for these decays is extremely small, but it can be enhanced if new physics increases the neutrino electromagnetic moments. The decay photons will appear as an interloper of astrophysical spectral lines. We propose that the neutrino-decay line can be identified with anisotropies in LIM clustering and also with the voxel intensity distribution. Ongoing and future LIM experiments will have -- depending on the neutrino hierarchy, transition and experiment considered -- a sensitivity to an effective electromagnetic transition moment $\sim 10^{-12}\, -\,10^{-8}\, (m_ic^2/{0.1 \rm eV})^{3/2}\mu_{\rm B}$, where $m_i$ is the mass of the decaying neutrino and $\mu_{\rm B}$ is the Bohr magneton. This will be significantly more sensitive than cosmic microwave background spectral distortions, and it will be competitive with stellar cooling studies. As a byproduct, we also report an analytic form of the one-point probability distribution function for neutrino-density fluctuations, obtained from the Quijote simulations using symbolic regression.


title: Mass-ratio and Magnetic Flux-Dependence of Modulated Accretion from Circumbinary Disks

abstract: Accreting supermassive binary black holes (SMBBHs) are potential multi-messenger sources because they emit both gravitational wave and electromagnetic (EM) radiation. Past work has shown that their EM output may be periodically modulated by an asymmetric density distribution in the circumbinary disk, often called an "overdensity" or "lump;" this modulation could possibly be used to identify a source as a binary. We explore the sensitivity of the overdensity to SMBBH mass ratio and magnetic flux through the accretion disk. We find that the relative amplitude of the overdensity and its associated EM periodic signal both degrade with diminishing mass ratio, vanishing altogether somewhere between 1:2 and 1:5. Greater magnetization also weakens the lump and any modulation of the light output. We develop a model to describe how lump formation results from internal stress degrading faster in the lump region than it can be rejuvenated through accretion inflow, and predicts a threshold value in specific internal stress below which lump formation should occur and which all our lump-forming simulations satisfy. Thus, detection of such a modulation would provide a constraint on both mass-ratio and magnetic flux piercing the accretion flow.


title: Universal dynamics of magnetic monopoles in two-dimensional kagom\'{e} ice

abstract: A magnetic monopole in spin ice is a novel quasiparticle excitation in condensed matter physics, and we found that the ac frequency dependent magnetic susceptibility $\chi(\omega)$ in the two-dimensional (2D) spin ice (so-called kagom\'{e} ice) of Dy$_2$Ti$_2$O$_7$ shows a single scaling form. This behavior can be understood in terms of the dynamical scaling law for 2D Coulomb gas (CG) systems [Phys. Rev. B 90, 144428 (2014)], characterized by the charge correlation length $\xi (\propto1/\sqrt{\omega_1})$, where $\omega_{1}$ is a characteristic frequency proportional to the peak position of the imaginary part of $\chi(\omega)$. It is a generic behavior among a wide variety of models such as the vortex dynamics of 2D superconductors, 2D superfluids, classical XY magnets, and dynamics of melting of Wigner crystals.


title: A Deep Learning Approach for Active Anomaly Detection of Extragalactic Transients

abstract: There is a shortage of multi-wavelength and spectroscopic followup capabilities given the number of transient and variable astrophysical events discovered through wide-field, optical surveys such as the upcoming Vera C. Rubin Observatory. From the haystack of potential science targets, astronomers must allocate scarce resources to study a selection of needles in real time. Here we present a variational recurrent autoencoder neural network to encode simulated Rubin Observatory extragalactic transient events using 1% of the PLAsTiCC dataset to train the autoencoder. Our unsupervised method uniquely works with unlabeled, real time, multivariate and aperiodic data. We rank 1,129,184 events based on an anomaly score estimated using an isolation forest. We find that our pipeline successfully ranks rarer classes of transients as more anomalous. Using simple cuts in anomaly score and uncertainty, we identify a pure (~95% pure) sample of rare transients (i.e., transients other than Type Ia, Type II and Type Ibc supernovae) including superluminous and pair-instability supernovae. Finally, our algorithm is able to identify these transients as anomalous well before peak, enabling real-time follow up studies in the era of the Rubin Observatory.


title: Interplay between powder catchment efficiency and layer height in self-stabilized laser metal deposition

abstract: In laser metal deposition (LMD) the height of the deposited track can vary within and between layers, causing significant deviations during the process evolution. Previous works have shown that in certain conditions a self-stabilizing mechanism occurs, maintaining a regular height growth and a constant standoff distance between the part and the deposition nozzle. Here we analyze the link between the powder catchment efficiency and the deposition height stability. To this purpose, a monitoring system was developed to study the deposition in different process conditions, using inline measurements of the specimen weight in combination with the layer height information obtained with coaxial optical triangulation. An analytical model was used to estimate the deposition efficiency in real-time from the height monitoring and the process parameters, which was verified by the direct mass measurements. The results show that the track height stabilization is associated to a reduction of the powder catchment efficiency, which is governed by the melt pool relative position with respect to the powder cone and the laser beam. For a given set of parameters, the standoff distance can be estimated to achieve the highest powder catchment efficiency and a regular height through the build direction.


title: Discovering features in gravitational-wave data through detector characterization, citizen science and machine learning

abstract: The observation of gravitational waves is hindered by the presence of transient noise (glitches). We study data from the third observing run of the Advanced LIGO detectors, and identify new glitch classes. Using training sets assembled by monitoring of the state of the detector, and by citizen-science volunteers, we update the Gravity Spy machine-learning algorithm for glitch classification. We find that a new glitch class linked to ground motion at the detector sites is especially prevalent, and identify two subclasses of this linked to different types of ground motion. Reclassification of data based on the updated model finds that 27 % of all transient noise at LIGO Livingston belongs to the new glitch class, making it the most frequent source of transient noise at that site. Our results demonstrate both how glitch classification can reveal potential improvements to gravitational-wave detectors, and how, given an appropriate framework, citizen-science volunteers may make discoveries in large data sets.


title: QCD factorization for twist-three axial-vector parton quasidistributions

abstract: The transverse component of the axial-vector correlation function of quark fields is a natural starting object for lattice calculations of twist-3 nucleon parton distribution functions. In this work we derive the corresponding factorization expression in terms of twist-2 and twist-3 collinear distributions to one-loop accuracy. The results are presented both in position space, as the factorization theorem for Ioffe-time distributions, and in momentum space, for the axial-vector quasi- and pseudodistributions.


title: Leveraging Spatial and Photometric Context for Calibrated Non-Lambertian Photometric Stereo

abstract: The problem of estimating a surface shape from its observed reflectance properties still remains a challenging task in computer vision. The presence of global illumination effects such as inter-reflections or cast shadows makes the task particularly difficult for non-convex real-world surfaces. State-of-the-art methods for calibrated photometric stereo address these issues using convolutional neural networks (CNNs) that primarily aim to capture either the spatial context among adjacent pixels or the photometric one formed by illuminating a sample from adjacent directions. In this paper, we bridge these two objectives and introduce an efficient fully-convolutional architecture that can leverage both spatial and photometric context simultaneously. In contrast to existing approaches that rely on standard 2D CNNs and regress directly to surface normals, we argue that using separable 4D convolutions and regressing to 2D Gaussian heat-maps severely reduces the size of the network and makes inference more efficient. Our experimental results on a real-world photometric stereo benchmark show that the proposed approach outperforms the existing methods both in efficiency and accuracy.


title: Dynamical Fragile Topology in Floquet Crystals

abstract: Although fragile topology has been intensely studied in static crystals in terms of Wannier obstruction, it is not clear how to generalize the concept to dynamical systems. In this work, we generalize the concept of fragile topology, and provide a definition of fragile topology for noninteracting Floquet crystals, which we refer to as dynamical fragile topology. In contrast to the static fragile topology defined by Wannier obstruction, dynamical fragile topology is defined for the nontrivial quantum dynamics characterized by the obstruction to static limits (OTSL). Specifically, the OTSL of a Floquet crystal is fragile if and only if it disappears after adding a symmetry-preserving static Hamiltonian in a direct-sum way preserving the relevant gaps (RGs). We further present a concrete 2+1D example for dynamical fragile topology, based on a model that is qualitatively the same as the dynamical model with anomalous chiral edge modes in [Rudner et al., Phys. Rev. X 3, 031005 (2013)]. The fragile OTSL in the 2+1D example exhibits anomalous chiral edge modes for a natural open boundary condition, and does not require any crystalline symmetries besides lattice translations. Our work paves the way to study fragile topology for general quantum dynamics.


title: Quantum field theory and the Bieberbach conjecture

abstract: An intriguing correspondence between ingredients in geometric function theory related to the famous Bieberbach conjecture (de Branges' theorem) and the non-perturbative crossing symmetric representation of 2-2 scattering amplitudes of identical scalars is pointed out. Using the dispersion relation and unitarity, we are able to derive several inequalities, analogous to those which arise in the discussions of the Bieberbach conjecture. We derive new and strong bounds on the ratio of certain Wilson coefficients and demonstrate that these are obeyed in one-loop $\phi^4$ theory, tree level string theory as well as in the S-matrix bootstrap. Further, we find two sided bounds on the magnitude of the scattering amplitude, which are shown to be respected in all the contexts mentioned above. Translated to the usual Mandelstam variables, for large $|s|$, fixed $t$, the upper bound reads $|\mathcal{M}(s,t)|\lesssim |s^2|$. We discuss how Szeg\"{o}'s theorem corresponds to a check of univalence in an EFT expansion, while how the Grunsky inequalities translate into nontrivial, nonlinear inequalities on the Wilson coefficients.


title: Approximating matrix eigenvalues by subspace iteration with repeated random sparsification

abstract: Traditional numerical methods for calculating matrix eigenvalues are prohibitively expensive for high-dimensional problems. Iterative random sparsification methods allow for the estimation of a single dominant eigenvalue at reduced cost by leveraging repeated random sampling and averaging. We present a general approach to extending such methods for the estimation of multiple eigenvalues and demonstrate its performance for several benchmark problems in quantum chemistry.

categories: math.NA cond-mat.str-el cs.NA physics.chem-ph


title: Non-Fermi liquid behavior in the Sachdev-Ye-Kitaev model for a one dimensional incoherent semimetal

abstract: Abstract We study a two-band dispersive Sachdev-Ye-Kitaev (SYK) model in 1 + 1 dimension. We suggest a model that describes a semimetal with quadratic dispersion at half-filling. We compute the Green's function at the saddle point using a combination of analytical and numerical methods. Employing a scaling symmetry of the Schwinger-Dyson equations that becomes transparent in the strongly dispersive limit, we show that the exact solution of the problem yields a distinct type of non-Fermi liquid with sublinear $\rho\propto T^{2/5}$ temperature dependence of the resistivity. A scaling analysis indicates that this state corresponds to the fixed point of the dispersive SYK model for a quadratic band touching semimetal.


title: Approximation of multipartite quantum states and the relative entropy of entanglement

abstract: Special approximation technique for analysis of different characteristics of states of multipartite infinite-dimensional quantum systems is proposed and applied to study of the relative entropy of entanglement and its regularisation. We prove several results about analytical properties of the multipartite relative entropy of entanglement and its regularization (the lower semicontinuity on wide class of states, the uniform continuity under the energy constraints, etc.). We establish a finite-dimensional approximation property for the relative entropy of entanglement and its regularization that allows to generalize to the infinite-dimensional case the results proved in the finite-dimensional settings.


title: The quest for scaling BFT Consensus through Tree-Based Vote Aggregation

abstract: With the growing commercial interest in blockchain, permissioned implementations have received increasing attention. Unfortunately, existing BFT consensus protocols that are the backbone of permissioned blockchains, either scale poorly or offer limited throughput. Most of these algorithms require at least one process to receive and validate the votes from all other processes and then broadcast the result, which is inherently non-scalable. Some algorithms avoid this bottleneck by using aggregation trees to collect and validate votes. However, to the best of our knowledge, such algorithms offer limited throughput and degrade quickly in the presence of faults. In this paper we propose \thesystem, the first BFT communication abstraction that organizes participants in a tree to perform scalable vote aggregation and that, in faulty runs, is able to terminate the protocol within an optimal number of reconfigurations ($f+1$). We define precisely which aggregation trees allow for optimal reconfiguration and show that, unlike previous protocols, when using these configurations, \thesystem scales to large number of processes and outperforms HotStuff's throughput by up to 38x.


title: On transference principle and Nesterenko's linear independence criterion

abstract: We consider the problem of simultaneous approximation of real numbers $\theta_1, \ldots,\theta_n$ with rationals and the dual problem of approximating zero with the values of the linear form $x_0+\theta_1x_1+\ldots+\theta_nx_n$ at integer points. In this setting we analyse two transference inequalities obtained by Schmidt and Summerer. We present a rather simple geometric observation, which proves their result. We also derive several corollaries previously unknown. Particularly, we show that, together with the transference inequalities for uniform exponents, Schmidt and Summerer's inequalities imply the inequalities by Bugeaud and Laurent and "one half" of the inequalities by Marnat and Moshchevitin. Besides that, we show that our main construction provides a rather simple proof of Nesterenko's linear independence criterion.


title: Skew-orthogonal polynomials in the complex plane and their Bergman-like kernels

abstract: Non-Hermitian random matrices with symplectic symmetry provide examples for Pfaffian point processes in the complex plane. These point processes are characterised by a matrix valued kernel of skew-orthogonal polynomials. We develop their theory in providing an explicit construction of skew-orthogonal polynomials in terms of orthogonal polynomials that satisfy a three-term recurrence relation, for general weight functions in the complex plane. New examples for symplectic ensembles are provided, based on recent developments in orthogonal polynomials on planar domains or curves in the complex plane. Furthermore, Bergman-like kernels of skew-orthogonal Hermite and Laguerre polynomials are derived, from which the conjectured universality of the elliptic symplectic Ginibre ensemble and its chiral partner follow in the limit of strong non-Hermiticity at the origin. A Christoffel perturbation of skew-orthogonal polynomials as it appears in applications to quantum field theory is provided.


title: End-to-End Trainable Multi-Instance Pose Estimation with Transformers

abstract: We propose an end-to-end trainable approach for multi-instance pose estimation, called POET (POse Estimation Transformer). Combining a convolutional neural network with a transformer encoder-decoder architecture, we formulate multiinstance pose estimation from images as a direct set prediction problem. Our model is able to directly regress the pose of all individuals, utilizing a bipartite matching scheme. POET is trained using a novel set-based global loss that consists of a keypoint loss, a visibility loss and a class loss. POET reasons about the relations between multiple detected individuals and the full image context to directly predict their poses in parallel. We show that POET achieves high accuracy on the COCO keypoint detection task while having less parameters and higher inference speed than other bottom-up and top-down approaches. Moreover, we show successful transfer learning when applying POET to animal pose estimation. To the best of our knowledge, this model is the first end-to-end trainable multi-instance pose estimation method and we hope it will serve as a simple and promising alternative.


title: Systematic benchmarking of HTTPS third party copy on 100Gbps links using XRootD

abstract: The High Luminosity Large Hadron Collider provides a data challenge. The amount of data recorded from the experiments and transported to hundreds of sites will see a thirty fold increase in annual data volume. A systematic approach to contrast the performance of different Third Party Copy(TPC) transfer protocols arises. Two contenders, XRootD-HTTPS and the GridFTP are evaluated in their performance for transferring files from one server to an-other over 100Gbps interfaces. The benchmarking is done by scheduling pods on the Pacific Research Platform Kubernetes cluster to ensure reproducible and repeatable results. This opens a future pathway for network testing of any TPC transfer protocol.


title: Revisiting Attenuation Curves: the Case of NGC 3351

abstract: Multi-wavelength images from the farUV (~0.15 micron) to the sub-millimeter of the central region of the galaxy NGC 3351 are analyzed to constrain its stellar populations and dust attenuation. Despite hosting a ~1 kpc circumnuclear starburst ring, NGC 3351 deviates from the IRX-beta relation, the relation between the infrared-to-UV luminosity ratio and the UV continuum slope (beta) that other starburst galaxies follow. To understand the reason for the deviation, we leverage the high angular resolution of archival nearUV-to-nearIR HST images to divide the ring into ~60-180 pc size regions and model each individually. We find that the UV slope of the combined intrinsic (dust-free) stellar populations in the central region is redder than what is expected for a young model population. This is due to the region's complex star formation history, which boosts the nearUV emission relative to the farUV. The resulting net attenuation curve has a UV slope that lies between those of the starburst attenuation curve (Calzetti et al. 2000) and the Small Magellanic Cloud extinction curve; the total-to-selective attenuation value, R'(V)=4.93, is larger than both. As found for other star-forming galaxies, the stellar continuum of NGC 3351 is less attenuated than the ionized gas, with E(B-V)_{star}=0.40 E(B-V)_{gas}. The combination of the `red' intrinsic stellar population and the new attenuation curve fully accounts for the location of the central region of NGC 3351 on the IRX-beta diagram. Thus, the observed characteristics result from the complex mixture of stellar populations and dust column densities in the circumnuclear region. Despite being a sample of one, these findings highlight the difficulty of defining attenuation curves of general applicability outside the regime of centrally-concentrated starbursts.


title: E$_9$ exceptional field theory II. The complete dynamics

abstract: We construct the first complete exceptional field theory that is based on an infinite-dimensional symmetry group. E$_9$ exceptional field theory provides a unified description of eleven-dimensional and type IIB supergravity covariant under the affine Kac-Moody symmetry of two-dimensional maximal supergravity. We present two equivalent formulations of the dynamics, which both rely on a pseudo-Lagrangian supplemented by a twisted self-duality equation. One formulation involves a minimal set of fields and gauge symmetries, which uniquely determine the entire dynamics. The other formulation extends $\mathfrak{e}_9$ by half of the Virasoro algebra and makes direct contact with the integrable structure of two-dimensional supergravity. Our results apply directly to other affine Kac-Moody groups, such as the Geroch group of general relativity.


title: Boundary Layers of Accretion Disks: Discovery of Vortex-Driven Modes and Other Waves

abstract: Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on their morphological characteristics in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects.


title: Triangular Lat-Igusa-Todorov algebras

abstract: Recently the authors D. Bravo, M. Lanzilotta, O. Mendoza and J. Vivero gave a generalization of the concept of Igusa-Todorov algebra and proved that those algebras, named Lat-Igusa-Todorov (LIT for short), satisfy the finitistic dimension conjecture. In this paper we explore the scope of that generalization and give conditions for a triangular matrix algebra to be LIT in terms of the algebras and the bimodule used in its definition. As an application we obtain that the tensor product of an LIT $\mathbb{K}$-algebra with a path algebra of a quiver whose underlying graph is a tree, is LIT.


title: On the moduli spaces of framed logarithmic connections on a Riemann surface

abstract: We describe some results on moduli space of logarithmic connections equipped with framings on a $n$-pointed compact Riemann surface.


title: Moving from Linear to Conic Markets for Electricity

abstract: We propose a new forward electricity market framework that admits heterogeneous market participants with second-order cone strategy sets, who accurately express the nonlinearities in their costs and constraints through conic bids, and a network operator facing conic operational constraints. In contrast to the prevalent linear-programming-based electricity markets, we highlight how the inclusion of second-order cone constraints improves uncertainty-, asset- and network-awareness of the market, which is key to the successful transition towards an electricity system based on weather-dependent renewable energy sources. We analyze our general market-clearing proposal using conic duality theory to derive efficient spatially-differentiated prices for the multiple commodities, comprising of energy and flexibility services. Under the assumption of perfect competition, we prove the equivalence of the centrally-solved market-clearing optimization problem to a competitive spatial price equilibrium involving a set of rational and self-interested participants and a price setter. Finally, under common assumptions, we prove that moving towards conic markets does not incur the loss of desirable economic properties of markets, namely market efficiency, cost recovery and revenue adequacy. Our numerical studies focus on the specific use case of uncertainty-aware market design and demonstrate that the proposed conic market brings advantages over existing alternatives within the linear programming market framework.


title: Sparticle and Higgs boson masses from the landscape: dynamical versus spontaneous supersymmetry breaking

abstract: Perturbative supersymmetry breaking on the landscape of string vacua is expected to favor large soft terms as a power-law or log distribution, but tempered by an anthropic veto of inappropriate vacua or vacua leading to too large a value for the derived weak scale -- a violation of the atomic principle. Indeed, scans of such vacua yield a statistical prediction for light Higgs boson mass m_h~ 125 GeV with sparticles (save possibly light higgsinos) typically beyond LHC reach. In contrast, models of dynamical SUSY breaking (DSB) -- with a hidden sector gauge coupling g^2 scanned uniformly -- lead to gaugino condensation and a uniform distribution of soft parameters on a log scale. Then soft terms are expected to be distributed as $m_{\rm soft}^{-1}$ favoring small values. A scan of DSB soft terms generally leads to $m_h\ll 125$ GeV and sparticle masses usually below LHC limits. Thus, the DSB landscape scenario seems excluded from LHC search results. An alternative is that the exponential suppression of the weak scale is set anthropically on the landscape via the atomic principle.


title: Origins and Demographics of Wandering Black Holes

abstract: We characterise the population of wandering black holes, defined as those physically offset from their halo centres, in the Romulus cosmological simulations. Unlike most other currently available cosmological simulations, black holes are seeded based on local gas properties and are permitted to evolve dynamically without being fixed at halo centres. Tracking these black holes allows us to make robust predictions about the offset population. We find that the number of wandering black holes scales roughly linearly with the halo mass, such that we expect thousands of wandering black holes in galaxy cluster halos. Locally, these wanderers account for around 10 per cent of the local black hole mass budget once seed masses are accounted for. Yet for higher redshifts ($z\gtrsim 4$), wandering black holes both outweigh and outshine their central supermassive counterparts. Most wandering black holes, we find, remain close to the seed mass and originate from the centres of previously disrupted satellite galaxies. While most do not retain a resolved stellar counterpart, those that do are situated farther out at larger fractions of the virial radius. Wanderers with higher luminosities are preferentially at lower radius, more massive, and either closer to their host's mid-planes or associated with a stellar overdensity. This analysis shows that our current census of supermassive black holes is incomplete and that a substantial population of off-centre wanderers likely exists.


title: Anomaly Detection Algorithms for Location Security in 5G Scenarios

abstract: Location based services are expected to play a major role in future generation cellular networks, starting from the incoming 5G systems. At the same time, localization technologies may be severely affected by attackers capable to deploy low cost fake base stations and use them to alter localization signals. In this paper, we concretely focus on two classes of threats: noise-like jammers, whose objective is to reduce the signal-to-noise ratio, and spoofing/meaconing attacks, whose objective is to inject false or erroneous information into the receiver. Then, we formulate the detection problems as binary hypothesis tests and solve them resorting to the generalized likelihood ratio test design procedure as well as the Latent Variable Models, which involves the expectation-maximization algorithm to estimate the unknown data distribution parameters. The proposed techniques can be applied to a large class of location data regardless the subsumed network architecture. The performance analysis is conducted over simulated data generated by using measurement models from the literature and highlights the effectiveness of the proposed approaches in detecting the aforementioned classes of attacks.


title: The HETDEX Survey: The Ly$\alpha$ Escape Fraction from 3D-HST Emission Line Galaxies at $z \sim 2$

abstract: We measure the Ly$\alpha$ escape fraction of 935 [OIII]-emitting galaxies between $1.9 < z < 2.35$ by comparing stacked spectra from the Hubble Space Telescope/WFC3's near-IR grism to corresponding stacks from the Hobby Eberly Telescope Dark Energy Experiment's Internal Data Release 2. By measuring the stacks' H$\beta$ to Ly$\alpha$ ratios, we determine the Ly$\alpha$ escape fraction as a function of stellar mass, star formation rate, internal reddening, size, and [OIII]/H$\beta$ ratio. We show that the escape fraction of Ly$\alpha$ correlates with a number of parameters, such as galaxy size, star formation rate, and nebular excitation. However, we also demonstrate that most of these relations are indirect, and the primary variables that control the escape of Ly$\alpha$ are likely stellar mass and internal extinction. Overall, the escape of Ly$\alpha$ declines from $\gtrsim 18\%$ in galaxies with $\log M/M_{\odot} \lesssim 9$ to $\lesssim 1\%$ for systems with $\log M/M_{\odot} \gtrsim 10$, with the sample's mean escape fraction being $6.0^{+0.6\%}_{-0.5\%}$.


title: SRB measures for Anosov actions

abstract: Given a general Anosov $\mathbb{R}^\kappa$ action on a closed manifold, we study properties of certain invariant measures that have recently been introduced in \cite{BGHW20} using the theory of Ruelle-Taylor resonances. We show that these measures share many properties of Sinai-Ruelle-Bowen measures for general Anosov flows such as smooth disintegrations along the unstable foliation, positive Lebesgue measure basins of attraction and a Bowen formula in terms of periodic orbits. Finally we show that if the action in the positive Weyl chamber is transitive, the measure is unique and has full support.


title: Meridional Circulation of Dust and Gas in the Circumstellar Disk: Delivery of Solids onto the Circumplanetary Region

abstract: We carried out 3D dust+gas radiative hydrodynamic simulations of forming planets. We investigated a parameter grid of Neptune-, Saturn-, Jupiter-, and 5 Jupiter-mass planets at 5.2, 30, 50 AU distance from their star. We found that the meridional circulation \citep{Szulagyi14,FC16} drives a strong vertical flow for the dust as well, hence the dust is not settled in the midplane, even for mm-sized grains. The meridional circulation will deliver dust and gas vertically onto the circumplanetary region, efficiently bridging over the gap. The Hill-sphere accretion rates for the dust are $\sim10^{-8}$ to $10^{-10}$ $\rm{M_{Jup}/yr}$, increasing with planet-mass. For the gas component, the gain is $10^{-6}$ to $10^{-8}$ $\rm{M_{Jup}/yr}$. The difference between the dust and gas accretion rates is smaller with decreasing planetary mass. In the vicinity of the planet, the mm-grains can get trapped easier than the gas, which means the circumplanetary disk might be enriched with solids in comparison to the circumstellar disk. We calculated the local dust-to-gas ratio (DTG) everywhere in the circumstellar disk and identified the altitude above the midplane where the DTG is 1, 0.1, 0.01, 0.001. The larger the planetary mass, the higher the mm-sized dust is delivered and a larger fraction of the dust disk is lifted by the planet. The stirring of mm-dust is negligible for Neptune-mass planets or below, but significant above Saturn-mass.


title: Tensor-Train Numerical Integration of Multivariate Functions with Singularities

abstract: Numerical integration is a classical problem emerging in many fields of science. Multivariate integration cannot be approached with classical methods due to the exponential growth of the number of quadrature nodes. We propose a method to overcome this problem. Tensor-train decomposition of a tensor approximating the integrand is constructed and used to evaluate a multivariate quadrature formula. We show how to deal with singularities in the integration domain and conduct theoretical analysis of the integration accuracy. The reference open-source implementation is provided.


title: Instability analysis for spiral arms of local galaxies: M51, NGC3627 and NGC628

abstract: We investigate dynamical states of grand-design spiral arms in three local galaxies: M51, NGC3627 and NGC628. Based on linear perturbation analysis considering multiple components in the galaxies, we compute instability parameters of the spiral arms using their observational data and argue whether the arms will fragment by their self-gravity. Our analysis utilises observations of carbon monoxide (CO), 21-centimetre line emission and multi-band photometric images for molecular gas, atomic gas and stellar components in the arms, respectively. We find that the grand-design arms of these galaxies indicate marginally stable states, and hence they are not on the way to fragment. We consider this to be consistent with the commonness of spiral galaxies and the relative rarity of fragmented discs at low redshifts. In the analysis, molecular gas is the dominant component to determine the (in)stability of the arms, whereas atomic gas and stars are far less important. Therefore, the results of our analysis are sensitive to an assumed CO-to-H$_{\rm 2}$ conversion factor. If we assume a typical scatter of the measurements and admit nearly twice as large a conversion factor as our fiducial value, our analysis results in predicting the instability for the spiral arms. More sophisticated determination of the conversion factor is required for more accurate analysis for the (in)stability of spiral arms.


title: The IoT Exchange

abstract: The IoT ecosystem suffers from a variety of problems around security, identity, access control, data flow and data storage that introduce friction into interactions between various parties. In many respects, the situation is similar to the early days of the Internet, where, prior to the establishment of Internet Exchanges, routing between different BGP autonomous systems was often point to point. We propose a similar solution, the IoT Exchange, where IoT device owners can register their devices and offer data for sale or can upload data into the IoT services of any of the big hyperscale cloud platforms for further processing. The goal of the IoT Exchange is to break down the silos within which device wireless connectivity types and cloud provider IoT systems constrain users to operate. In addition, if the device owner needs to maintain the data close to the edge to reduce access latency, the MillenniumDB service running in an edge data center with minimal latency to the edge device, provides a database with a variety of schema engines (SQL, noSQL, etc). The IoT exchange uses decentralized identifiers for identity management and verifiable credentials for authorizing software updates and to control access to the devices, to avoid dependence on certificate authorities and other centralized identity and authorization management systems. In addition, verifiable credentials provide a way whereby privacy preserving processing can be applied to traffic between a device and an end data or control customer, if some risk of privacy compromise exists.


title: Impulsive fire disturbance in a savanna model: Tree-grass coexistence states, multiple stable system states, and resilience

abstract: Savanna ecosystems are shaped by the frequency and intensity of regular fires. We model savannas via an ordinary differential equation (ODE) encoding a one-sided inhibitory Lotka-Volterra interaction between trees and grass. By applying fire as a discrete disturbance, we create an impulsive dynamical system that allows us to identify the impact of variation in fire frequency and intensity. The model exhibits three different bistability regimes: between savanna and grassland; two savanna states; and savanna and woodland. The impulsive model reveals rich bifurcation structures in response to changes in fire intensity and frequency -- structures that are largely invisible to analogous ODE models with continuous fire. In addition, by using the amount of grass as an example of a socially-valued function of the system state, we examine the resilience of the social value to different disturbance regimes. We find that large transitions ("tipping") in the valued quantity can be triggered by small changes in disturbance regime.


title: Numerical Simulation of Vortex-Induced Vibration With Bistable Springs : Consistency with the Equilibrium Constraint

abstract: We present results from 2-D numerical simulations based on Immersed Boundary Method of a cylinder in uniform fluid flow attached to bistable springs undergoing Vortex-Induced Vibrations (VIV). The elastic spring potential for the bistable springs, consisting of 2 potential wells, is defined by the spacing between the potential minima and the depth of the potential wells. We perform simulations of VIV with linear spring, as well as bistable springs with two different inter-well separations, over a wide range of reduced velocity. As expected, large oscillation amplitudes correspond to lock-in of the lift force with the natural frequency of the spring-mass system. The range of reduced velocity over which lock-in occurs is significantly higher for VIV with bistable springs compared to VIV with linear springs, although the maximum possible amplitude appears to be independent of the spring type. For VIV with bistable springs, the cylinder undergoes double-well oscillations in the lock-in regime. Range of reduced velocity over which lock-in occurs increases when the inter-well distance is reduced. The vortex shedding patterns and amplitude trends look similar at the same equivalent reduced velocity for the different springs. The results here are consistent with our prior theory, in which we propose a new "Equilibrium-Constraint (EC)" based on average kinetic energy budget of the structure. For a given spring potential, the intersection of natural frequency curves with the EC curve yields the possible range of reduced velocities over which lock-in should occur. Our numerical simulations show a collapse of the amplitude-versus-structure frequency data for all the simulations onto roughly the same curve, thus supporting the existence of the EC, and providing an explanation for the trends in the VIV oscillations.


title: A Joint Reinforcement-Learning Enabled Caching and Cross-Layer Network Code for Sum-Rate Maximization in F-RAN with D2D Communications

abstract: In this paper, we leverage reinforcement learning (RL) and cross-layer network coding (CLNC) for efficiently pre-fetching users' contents to the local caches and delivering these contents to users in a downlink fog-radio access network (F-RAN) with device-to-device (D2D) communications. In the considered system, fog access points (F-APs) and cache-enabled D2D (CE-D2D) users are equipped with local caches for alleviating traffic burden at the fronthaul, while users' contents can be easily and quickly accommodated. In CLNC, the coding decisions take users' contents, their rates, and power levels of F-APs and CE-D2D users into account, and RL optimizes caching strategy. Towards this goal, a joint content placement and delivery problem is formulated as an optimization problem with a goal to maximize system sum-rate. For this NP-hard problem, we first develop an innovative decentralized CLNC coalition formation (CLNC-CF) algorithm to obtain a stable solution for the content delivery problem, where F-APs and CE-D2D users utilize CLNC resource allocation. By taking the behavior of F-APs and CE-D2D users into account, we then develop a multi-agent RL (MARL) algorithm for optimizing the content placements at both F-APs and CE-D2D users. Simulation results show that the proposed joint CLNC-CF and RL framework can effectively improve the sum-rate by up to 30\%, 60\%, and 150\%, respectively, compared to: 1) an optimal uncoded algorithm, 2) a standard rate-aware-NC algorithm, and 3) a benchmark classical NC with network-layer optimization.