title: A new approach to the characteristic polynomial of a random unitary matrix

abstract: Since the seminal work of Keating and Snaith, the characteristic polynomial of a random Haar-distributed unitary matrix has seen several of its functional studied or turned into a conjecture; for instance: $ \bullet $ its value in $1$ (Keating-Snaith theorem), $ \bullet $ the truncation of its Fourier series up to any fraction of its degree, $ \bullet $ the computation of the relative volume of the Birkhoff polytope, $ \bullet $ its products and ratios taken in different points, $ \bullet $ the product of its iterated derivatives in different points, $ \bullet $ functionals in relation with sums of divisor functions in $ \mathbb{F}_q[X] $. $ \bullet $ its mid-secular coefficients, $ \bullet $ the "moments of moments", etc. We revisit or compute for the first time the asymptotics of the integer moments of these last functionals and several others. The method we use is a very general one based on reproducing kernels, a symmetric function generalisation of some classical orthogonal polynomials interpreted as the Fourier transform of particular random variables and a local Central Limit Theorem for these random variables. We moreover provide an equivalent paradigm based on a randomisation of the mid-secular coefficients to rederive them all. These methodologies give a new and unified framework for all the considered limits and explain the apparition of Hankel determinants or Wronskians in the limiting functional.


title: Algorithms and Hardness for Linear Algebra on Geometric Graphs

abstract: For a function $\mathsf{K} : \mathbb{R}^{d} \times \mathbb{R}^{d} \to \mathbb{R}_{\geq 0}$, and a set $P = \{ x_1, \ldots, x_n\} \subset \mathbb{R}^d$ of $n$ points, the $\mathsf{K}$ graph $G_P$ of $P$ is the complete graph on $n$ nodes where the weight between nodes $i$ and $j$ is given by $\mathsf{K}(x_i, x_j)$. In this paper, we initiate the study of when efficient spectral graph theory is possible on these graphs. We investigate whether or not it is possible to solve the following problems in $n^{1+o(1)}$ time for a $\mathsf{K}$-graph $G_P$ when $d < n^{o(1)}$: $\bullet$ Multiply a given vector by the adjacency matrix or Laplacian matrix of $G_P$ $\bullet$ Find a spectral sparsifier of $G_P$ $\bullet$ Solve a Laplacian system in $G_P$'s Laplacian matrix For each of these problems, we consider all functions of the form $\mathsf{K}(u,v) = f(\|u-v\|_2^2)$ for a function $f:\mathbb{R} \rightarrow \mathbb{R}$. We provide algorithms and comparable hardness results for many such $\mathsf{K}$, including the Gaussian kernel, Neural tangent kernels, and more. For example, in dimension $d = \Omega(\log n)$, we show that there is a parameter associated with the function $f$ for which low parameter values imply $n^{1+o(1)}$ time algorithms for all three of these problems and high parameter values imply the nonexistence of subquadratic time algorithms assuming Strong Exponential Time Hypothesis ($\mathsf{SETH}$), given natural assumptions on $f$. As part of our results, we also show that the exponential dependence on the dimension $d$ in the celebrated fast multipole method of Greengard and Rokhlin cannot be improved, assuming $\mathsf{SETH}$, for a broad class of functions $f$. To the best of our knowledge, this is the first formal limitation proven about fast multipole methods.


title: On Finsler surfaces with certain flag curvatures

abstract: In the present paper, we find out necessary and sufficient conditions for a Finsler surface $(M,F)$ to be Landsbregian in terms of the Berwald curvature $2$-forms. We study Finsler surfaces which satisfy some flag curvature $K$ conditions, viz., $V(K)=0,\,\,V(K)= -\mathcal{I}/F^2$ and $V(K)=-\mathcal{I}\,K,$ where $\mathcal{I}$ is the Cartan scalar. In order to do so, we investigate some geometric objects associated with the global Berwald distribution $\mathcal{D}:= \operatorname{span}\{S, H, V:=JH\}$ of a $2$-dimensional Finsler metrizable nonflat spray $S$. We obtain some classifications of such surfaces and show that under what hypothesis these surfaces turn to be Riemannian. The existence of a first integral for the geodesic flow in each case has some remarkable consequences concerning rigidity results. We prove that a Finsler surface with $V(K)= -\mathcal{I}/F^2$ and either $S(K)=0$ or $S(\mathcal{J})=0$ is Riemannian. Further, a Finsler surface with $V(K)=-\mathcal{I}\,K$ and $S(K)=0$ is Riemannian.


title: Geometric Solution of Image Degradation by Diffraction in Lensless Sensing and Microscopy

abstract: This paper proposes a non-computational method of counteracting the effect of image degradation introduced by the diffraction phenomenon in lensless microscopy. All the optical images (whether focused by lenses or not) are diffraction patterns, which preserve the visual information upto a certain extent determined by the size of the point spread functions, like airy disks in some cases. A highly diverging beam can be exploited to reduce the spatial extent of these point spread functions relatively in the transformed projective space, which can help us in the spatial unmixing of the visual information. The principle has been experimentally validated by the lensless imaging of red blood cells of diameter ~6-9 micrometers and a photolithography mask with features in micrometer scale. The important advantages of the proposed approach of non-computational shadow microscopy are the improved depth of field and a drastic increase in the sensor to sample working distance. The imaging method can also be used as a projection technique in the multi-angle optical computed tomography (CT).


title: Fixing the vector coupling constant $g_\rho = 5.0$ in the NJL model and final-state interactions in processes $\rho \to e^+e^-[\mu^+\mu^-]$, $\rho \to \pi^+\pi^-$, $\tau \to \pi^-\pi^0\nu_\tau$, and $e^+e^-\to\pi^+\pi^-$

abstract: The possibility to use the width of the decay $\rho \to e^{+}e^{-}$ to fix the input parameter $g_\rho=5.0$ of the $SU(2) \times SU(2)$ chiral-symmetric Nambu--Jona-Lasinio model is discussed. It is shown that for a consistent simultaneous description of the processes $\rho \to e^{+}e^{-}$, $\rho \to \pi^{+}\pi^{-}$, $\tau^{-} \to \pi^{-}\pi^{0} \nu_{\tau}$, and $e^{+}e^{-} \to \pi^{+}\pi^{-}$ can be constructed. Taking into account the interaction of pions in the final state appears to be important. The obtained theoretical results for the considered processes are in a satisfactory agreement with experimental data.


title: MESA models with magnetic braking

abstract: Two magnetic braking models are implemented in MESA for use in the MIST stellar model grids. Stars less than about 1.3 $M_{\odot}$ are observed to spin down over time through interaction with their magnetized stellar winds (i.e., magnetic braking). This is the basis for gyrochronology, and fundamental to the evolution of lower mass stars. The detailed physics behind magnetic braking are uncertain, as are 1D stellar evolution models. Thus, we calibrate our models and compare to data from open clusters. Each braking model tested here is capable of reproducing the data, albeit with some important distinctions. The Matt et al. (2015) prescription matches the slowly rotating stars observed in open clusters, but tends to overestimate the presence of rapidly rotating stars. The Garraffo et al. (2018) prescription often produces too much angular momentum loss to accurately match the observed slow sequence for lower mass stars, but reproduces the bimodal nature of slow and rapidly rotating stars observed in open clusters fairly well. We find additional evidence that some level of mass dependency may be missing in these braking models to match the rotation periods observed in clusters older than 1 Gyr better.


title: Charge Density Wave Transport in Porous Graphene Nanoribbons

abstract: Porous graphene (PG) forms a class of graphene-related materials with nanoporous architectures. Their unique atomic arrangements present interconnected networks with high surface area and high pore volume. Some remarkable properties of PG, such as high mechanical strength and good thermal stability, have been widely studied. However, their electrical conductivity, and most importantly, their charge transport mechanism are still not fully understood. Herein, we employed a numerical approach based on a 2D tight-binding model Hamiltonian to first reveal the nature of the charge transport mechanism in PG nanoribbons. Results showed that the charge transport in these materials is mediated by charge density waves. These carrier species are dynamically stable and present very shallow lattice distortions. The porosity allows for an alternative to the usual arising of polaron-like charge carriers and it can preserve the PG semiconducting character even in broader nanoribbons. The charge density waves move in PG within the optical regime with terminal velocities varying from 0.50 up to 1.15 A/fs. These velocities are lower than the ones for polarons in conventional graphene nanoribbons (2.2-5.1 A/fs).


title: MAXI J1820+070 with NuSTAR II. Flaring during the hard to soft state transition with a long soft lag

abstract: We continue the analysis of NuSTAR data from the recent discovery outburst of MAXI J1820+070 (optical counterpart ASASSN-18ey), focussing on an observation including unusual flaring behaviour during the hard to soft state transition, which is a short phase of outbursts and so comparatively rarely observed. Two plateaus in flux are separated by a variable interval lasting ~10 ks, which shows dipping then flaring stages. The variability is strongest (with fractional variability up to $F_{\rm Var}\sim10\%$) at high energies and reduces as the contribution from disc emission becomes stronger. Flux-resolved spectra show that the variability is primarily due to the power law flux changing. We also find a long soft lag of the thermal behind the power law emission, which is $20_{-1.2}^{+1.6}$ s during the flaring phase. The lag during the dipping stage has a different lag-energy spectrum, which may be due to a wave passing outwards through the disc. Time resolved spectral fitting suggests that the lag during the flaring stage may be due to the disc re-filling after being disrupted to produce the power law flare, perhaps related to the system settling after the jet ejection which occurred around 1 day before. The timescales of these phenomena imply a low viscosity parameter, $\alpha\sim10^{-3}$, for the inner region of the disc.


title: What Skills do IT Companies look for in New Developers? A Study with Stack Overflow Jobs

abstract: Context: There is a growing demand for information on how IT companies look for candidates to their open positions. Objective: This paper investigates which hard and soft skills are more required in IT companies by analyzing the description of 20,000 job opportunities. Method: We applied open card sorting to perform a high-level analysis on which types of hard skills are more requested. Further, we manually analyzed the most mentioned soft skills. Results: Programming languages are the most demanded hard skills. Communication, collaboration, and problem-solving are the most demanded soft skills. Conclusion: We recommend developers to organize their resum\'e according to the positions they are applying. We also highlight the importance of soft skills, as they appear in many job opportunities.


title: Electronic structure correspondence of singlet-triplet scale separation in strained Sr2RuO4

abstract: At a temperature of roughly 1\,K, \ce{Sr2RuO4} undergoes a transition from a normal Fermi liquid to a superconducting phase. Even while the former is relatively simple and well understood, the superconducting state is not even after 25 years of study. More recently it has been found that critical temperatures can be enhanced by application of uniaxial strain, up to a critical strain, after which it falls off. In this work, we take an `instability' approach and seek for divergences in susceptibilities. This provides an unbiased way to distinguish tendencies to competing ground states. We show that in the unstrained compound the singlet and triplet instabilities of the normal Fermi liquid phase are closely spaced. Under uniaxial strain electrons residing on all orbitals contributing to the Fermiology become more coherent while the electrons of Ru-$d_{xy}$ character become heavier and electrons of Ru-$d_{xz,yz}$ characters become lighter. In the process, Im\,$\chi(\mathbf{q},\omega)$ increases rapidly around the incommensurate vector $\mathbf{q}{=}(0.3,0.3,0)2\pi/a$ while it gets suppressed at all other commensurate vectors, in particular at $q{=}0$, which is essential for spin-triplet superconductivity. Thus the triplet superconducting instability remains the lagging instability of the system and the singlet instability enhances under strain, leading to a large energy-scale separation between these competing instabilities. At large strain an instability to a spin density wave overtakes the superconducting one. The analysis relies on a high-fidelity, \emph{ab initio} description of the one-particle properties and two-particle susceptibilities, based on the Quasiparticle Self-Consistent \emph{GW} approximation augmented by Dynamical Mean Field theory. This approach is described and its high fidelity confirmed by comparing to observed one- and two-particle properties.


title: CMB lensing power spectrum estimation without instrument noise bias

abstract: The power spectrum of cosmic microwave background (CMB) lensing will be measured to sub-percent precision with upcoming surveys, enabling tight constraints on the sum of neutrino masses and other cosmological parameters. Measuring the lensing power spectrum involves the estimation of the connected trispectrum of the four-point function of the CMB map, which requires the subtraction of a large Gaussian disconnected noise bias. This reconstruction noise bias receives contributions both from CMB and foreground fluctuations as well as instrument noise (both detector and atmospheric noise for ground-based surveys). The debiasing procedure therefore relies on the quality of simulations of the instrument noise which may be expensive or inaccurate. We propose a new estimator that makes use of at least four splits of the CMB maps with independent instrument noise. This estimator makes the CMB lensing power spectrum completely insensitive to any assumptions made in modeling or simulating the instrument noise. We show that this estimator, in many practical situations, leads to no substantial loss in signal-to-noise. We provide an efficient algorithm for its computation that scales with the number of splits $m$ as $\mathcal{O}(m^2)$ as opposed to a naive $\mathcal{O}(m^4)$ expectation.


title: On Caristi fixed point theorem for set-valued mappings

abstract: The aim of this paper is to discuss Penot's problem on a generalization of Caristi's fixed point theorem. We settle this problem in the negative and we present some new theorems on the existence of fixed points of set-valued mappings in ordered metric spaces.


title: Bosonic Pfaffian State in the Hofstadter-Bose-Hubbard Model

abstract: Topological states of matter, such as fractional quantum Hall states, are an active field of research due to their exotic excitations. In particular, ultracold atoms in optical lattices provide a highly controllable and adaptable platform to study such new types of quantum matter. However, finding a clear route to realize non-Abelian quantum Hall states in these systems remains challenging. Here we use the density-matrix renormalization-group (DMRG) method to study the Hofstadter-Bose-Hubbard model at filling factor $\nu = 1$ and find strong indications that at $\alpha=1/6$ magnetic flux quanta per plaquette the ground state is a lattice analog of the continuum non-Abelian Pfaffian. We study the on-site correlations of the ground state, which indicate its paired nature at $\nu = 1$, and find an incompressible state characterized by a charge gap in the bulk. We argue that the emergence of a charge density wave on thin cylinders and the behavior of the two- and three-particle correlation functions at short distances provide evidence for the state being closely related to the continuum Pfaffian. The signatures discussed in this letter are accessible in current cold atom experiments and we show that the Pfaffian-like state is readily realizable in few-body systems using adiabatic preparation schemes.


title: MARTY -- Modern ARtificial Theoretical phYsicist: A C++ framework automating symbolic calculations Beyond the Standard Model

abstract: Studies Beyond the Standard Model (BSM) will become more and more important in the near future with a rapidly increasing amount of data from different experiments around the world. The full study of BSM models is in general an extremely time-consuming task involving long and difficult calculations. It is in practice not possible to do exhaustive predictions in these models by hand, in particular if one wants to perform a statistical comparison with data and the SM. Here we present MARTY (Modern ARtificial Theoretical phYsicist), a new C++ framework that fully automates calculations from the Lagrangian to physical quantities such as amplitudes or cross-sections. This framework can fully simplify, automatically and symbolically, physical quantities in a very large variety of models. MARTY can also compute Wilson coefficients in effective theories. This will considerably facilitate the study of BSM models in flavor physics. Contrary to the existing public codes in this field MARTY aims to give a unique, free, open-source, powerful and user-friendly tool for high-energy physicists studying predictive BSM models, in effective or full theories up to the 1-loop level, which does not rely on any external package. With a few lines of code one can gather final expressions that may be evaluated numerically for statistical analysis. Features like automatic generation and manual edition of Feynman diagrams, comprehensive manual and documentation, clear and easy to handle user interface are amongst notable features of MARTY.


title: Development and operation of an electrostatic time-of-flight detector for the Rare RI storage Ring

abstract: An electrostatic time-of-flight detector named E-MCP has been developed for quick diagnostics of circulating beam and timing measurement in mass spectrometry at the Rare-RI Ring in RIKEN. The E-MCP detector consists of a conversion foil, potential grids, and a microchannel plate. Secondary electrons are released from the surface of the foil when a heavy ion hits it. The electrons are accelerated and deflected by 90$^\circ$ toward the microchannel plate by electrostatic potentials. A thin carbon foil and a thin aluminum-coated mylar foil were used as conversion foils. We obtained time resolutions of 69(1) ps and 43(1) ps (standard deviation) for a $^{84}$Kr beam at an energy of 170 MeV/u when using the carbon and the aluminum-coated mylar foils, respectively. A detection efficiency of approximately 90% was obtained for both foils. The E-MCP detector equipped with the carbon foil was installed inside the Rare-RI Ring to confirm particle circulation within a demonstration experiment on mass measurements of nuclei around $^{78}$Ge produced by in-flight fission of uranium beam at the RI Beam Factory in RIKEN. Periodic time signals from circulating ions were clearly observed. Revolution times for $^{78}$Ge, $^{77}$Ga, and $^{76}$Zn were obtained. The results confirmed successful circulation of the short-lived nuclei inside the Rare-RI Ring.


title: Disk fragmentation and intermittent accretion onto supermassive stars

abstract: Supermassive stars (SMSs) with $\sim10^{4-5}~\mathrm{M}_{\odot}$ are candidate objects for the origin of supermassive black holes observed at redshift $z$>6. They are supposed to form in primordial-gas clouds that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and its quiescent periods are longer than the Kelvin-Helmholtz (KH) timescales at the stellar surfaces. In this paper, we examine the role of the ionizing radiation feedback based on the accretion history in two possible SMS-forming clouds extracted from cosmological simulations, following their evolution with vertically-integrated two-dimensional hydrodynamic simulations with detailed thermal and chemical models. The consistent treatment of the gas thermal evolution is crucial for obtaining the realistic accretion history, as we demonstrate by performing an additional run with a barotropic equation of state, in which the fluctuation of the accretion rate is artificially suppressed. We find that although the accretion becomes intermittent due to the formation of spiral arms and clumps in gravitationally unstable disks, the quiescent periods are always shorter than the KH timescales, implying that SMSs can form without affected by the ionizing radiation.


title: Anisotropic Fermionic Quasiparticles

abstract: We have carried out a comprehensive investigation of the quasiparticle properties of a two-dimensional electron gas, interacting via the long-range Coulomb interaction, in the presence of bare mass anisotropy (i.e. with an elliptic noninteracting Fermi surface) by calculating the self-energy, the spectral function, the scattering rate, and the effective mass within the leading order dynamical self-energy approximation. We find novel anisotropic features of quasiparticle properties that are not captured by the isotropic approximation where the anisotropic effective mass is replaced by the isotropic averaged density-of-states mass. Some of these interesting results are: (1) the renormalization of the quasiparticle spectrum becomes highly anisotropic as the quasiparticle energy increases away from the Fermi energy; (2) the inelastic scattering rate features a strong anisotropy, exhibiting an abrupt jump at different injected energies depending on the momentum direction of the injected electron; (3) the effective mass anisotropy is reduced by interactions. Our results and analysis show that the unjustified neglect of the mass anisotropy can lead to an incorrect description of quasiparticle properties of the anisotropic system although the use of an equivalent isotropic approximation using the density-of-states effective mass works as a reasonable approximation in many situations. We also provide a theory using the plasmon-pole approximation, commenting on its validity for anisotropic self-energy calculations. We comment also on the interaction effect on the Fermi surface topology, finding that the elliptic shape of the bare Fermi surface is preserved, with suppressed ellipticity, in the interacting system to a high degree of accuracy. Our theory provides a complete generalization of the existing isotropic many-body theory of interacting electrons to the corresponding anisotropic systems.


title: Universality in asymptotic bounds and its saturation in $2$D CFT

abstract: We study asymptotics of three point coefficients (light-light-heavy) and two point correlators in heavy states in unitary, compact $2$D CFTs. We prove an upper and lower bound on such quantities using numerically assisted Tauberian techniques. We obtain an optimal upper bound on the spectrum of operators appearing with fixed spin from the OPE of two identical scalars. While all the CFTs obey this bound, rational CFTs come close to saturating it. This mimics the scenario of bounds on asymptotic density of states and thereby pronounces an universal feature in asymptotics of 2D CFTs. Next, we clarify the role of smearing in interpreting the asymptotic results pertaining to considerations of eigenstate thermalization in 2D CFTs. In the context of light-light-heavy three point coefficients, we find that the order one number in the bound is sensitive to how close the light operators are from the $\frac{c}{32}$ threshold. In context of two point correlator in heavy state, we find the presence of an enigmatic regime which separates the $AdS_3$ thermal physics and the BTZ black hole physics. Furthermore, we present some new numerical results on the behavior of spherical conformal block.


title: Starburst galaxies strike back: a multi-messenger analysis with Fermi-LAT and IceCube data

abstract: Starburst galaxies, which are known as "reservoirs" of high-energy cosmic-rays, can represent an important high-energy neutrino "factory" contributing to the diffuse neutrino flux observed by IceCube. In this paper, we revisit the constraints affecting the neutrino and gamma-ray hadronuclear emissions from this class of astrophysical objects. In particular, we go beyond the standard prototype-based approach leading to a simple power-law neutrino flux, and investigate a more realistic model based on a data-driven blending of spectral indexes, thereby capturing the observed changes in the properties of individual emitters. We then perform a multi-messenger analysis considering the extragalactic gamma-ray background (EGB) measured by Fermi-LAT and different IceCube data samples: the 7.5-year High-Energy Starting Events (HESE) and the 6-year high-energy cascade data. Along with starburst galaxies, we take into account the contributions from blazars and radio galaxies as well as the secondary gamma-rays from electromagnetic cascades. Remarkably, we find that, differently from the highly-constrained prototype scenario, the spectral index blending allows starburst galaxies to account for up to $40\%$ of the HESE events at $95.4\%$ CL, while satisfying the limit on the non-blazar EGB component. Moreover, values of $\mathcal{O}(100~\mathrm{PeV})$ for the maximal energy of accelerated cosmic-rays by supernovae remnants inside the starburst are disfavoured in our scenario. In broad terms, our analysis points out that a better modeling of astrophysical sources could alleviate the tension between neutrino and gamma-ray data interpretation.


title: Functional Prescription for EFT Matching

abstract: We simplify the one-loop functional matching formalism to develop a streamlined prescription. The functional approach is conceptually appealing: all calculations are performed within the UV theory at the matching scale, and no prior determination of an Effective Field Theory (EFT) operator basis is required. Our prescription accommodates any relativistic UV theory that contains generic interactions (including derivative couplings) among scalar, fermion, and vector fields. As an example application, we match the singlet scalar extended Standard Model (SM) onto SMEFT.


title: The GALAH survey: Chemical homogeneity of the Orion complex

abstract: Due to its proximity, the Orion star forming region is often used as a proxy to study processes related to star formation and to observe young stars in the environment they were born in. With the release of Gaia DR2, the distance measurements to the Orion complex are now good enough that the three dimensional structure of the complex can be explored. Here we test the hypothesis that, due to non-trivial structure and dynamics, and age spread in the Orion complex, the chemical enrichment of youngest stars by early core-collapse supernovae can be observed. We obtained spectra of 794 stars of the Orion complex with the HERMES spectrograph at the Anglo Australian telescope as a part of the GALAH and GALAH-related surveys. We use the spectra of $\sim300$ stars to derive precise atmospheric parameters and chemical abundances of 25 elements for 15 stellar clusters in the Orion complex. We demonstrate that the Orion complex is chemically homogeneous and that there was no self-pollution of young clusters by core-collapse supernovae from older clusters; with a precision of 0.02 dex in relative alpha-elements abundance and 0.06 dex in oxygen abundance we would have been able to detect pollution from a single supernova, given a fortunate location of the SN and favourable conditions for ISM mixing. We estimate that the supernova rate in the Orion complex was very low, possibly producing no supernova by the time the youngest stars of the observed population formed (from around 21 to 8 Myr ago).


title: Testing Leptoquark/EFT in $\bar B \to D^{(*)} l\bar\nu$ at the LHC

abstract: We investigate the current LHC bounds on New Physics (NP) that contributes to $\bar B \to D^{(*)} l\bar\nu$ for $l = (e,\mu,\tau)$ by considering both leptoquark (LQ) models and an effective-field-theory (EFT) Hamiltonian. Experimental analyses from $l+\text{missing}$ searches with high $p_T$ are applied to evaluate the NP constraints with respect to the Wilson coefficients. A novel point of this work is to show difference between LQs and EFT for the applicable LHC bound. In particular, we find that the EFT description is not valid to search for LQs with the mass less than $\lesssim 10\,\text{TeV}$ at the LHC and leads to overestimated bounds. We also discuss future prospects of high luminosity LHC searches including the charge asymmetry of background and signal events. Finally, a combined summary for the flavor and LHC bounds is given, and then we see that in several NP scenarios the LHC constraints are comparable with the flavor ones.


title: Intimate Relationship Between Sterile Neutrino Dark Matter and $\Delta N_{\rm eff}$

abstract: The self-interacting neutrino hypothesis is well motivated for addressing the tension between the origin of sterile neutrino dark matter and indirect detection constraints. It can also result in a number of testable signals from the laboratories to the cosmos. We explore a model of neutrino self-interaction mediated by a Majoron-like scalar with sub-MeV mass, and show that explaining the relic density of sterile neutrino dark matter implies a lower bound on the amount of extra radiation in early universe, in particular $\Delta N_{\rm eff}>0.12$ at the CMB epoch. This lower bound will be further strengthened with an improved $X$-ray search at the Athena observatory. Such an intimate relationship will be unambiguously tested by the upcoming CMB-S4 project.


title: Generation of gravitational waves and tidal disruptions in clumpy galaxies

abstract: Obtaining a better understanding of intermediate-mass black holes (IMBHs) is crucial, as their properties could shed light on the origin and growth of their supermassive counterparts. Massive star-forming clumps, which are present in a large fraction of massive galaxies at $z \sim$ 1-3, are amongst the venues wherein IMBHs could reside. We perform a series of Fokker-Planck simulations to explore the occurrence of tidal disruption (TD) and gravitational wave (GW) events about an IMBH in a massive star-forming clump, modelling the latter so that its mass ($10^8 \,{\rm M}_{\odot}$) and effective radius ($100$ pc) are consistent with the properties of both observed and simulated clumps. We find that the TD and GW event rates are in the ranges $10^{-6}$-$10^{-5}$ and $10^{-8}$-$10^{-7}$ yr$^{-1}$, respectively, depending on the assumptions for the initial inner density profile of the system ($\rho \propto r^{-2}$ or $\propto r^{-1}$) and the initial mass of the central IMBH ($10^5$ or $10^3\,{\rm M}_{\odot}$). By integrating the GW event rate over $z$ = 1-3, we expect that the Laser Interferometer Space Antenna will be able to detect $\sim$2 GW events per yr coming from these massive clumps; the intrinsic rate of TD events from these systems amounts instead to a few $10^3$ per yr, a fraction of which will be observable by, e.g. the Square Kilometre Array and the Advanced Telescope for High Energy Astrophysics. In conclusion, our results support the idea that the forthcoming GW and electromagnetic facilities may have the unprecedented opportunity of unveiling the lurking population of IMBHs.


title: Gravitational origin of dark matter and Majorana neutrino mass with non-minimal quartic inflation

abstract: We propose a minimal framework to address successful quartic inflation, dark matter (DM) production in the early universe and non-vanishing tiny Majorana neutrino mass from a common gravitational origin point of view. In this setup, the quartic inflation is revived successfully via non-minimal coupling of inflaton to gravity while the production of DM takes place from purely gravitational effects through a misalignment mechanism. The generation of light Majorana neutrino mass is aided by explicit breaking of global lepton number symmetry through Planck suppressed operators involving non-zero vacuum expectation value (VEV) of the inflaton field. We present a detailed study of the DM yield in presence of non-minimal inflation, considering both the metric and the Palatini formalisms of gravity wherever appropriate. We reach at some interesting and different results in the DM sector compared to the earlier works in the similar direction with minimal inflationary background. Restricting to the light DM regime ($\mathcal{O}(1)$ keV- $\mathcal{O}(100)$ MeV) where classical production is expected to dominate over the quantum production, we numerically predict the DM mass by varying the DM quartic and non-minimal coupling, to be consistent with relic density requirements. We also obtain some non-trivial dependence of DM phenomenology on some of the relevant parameters of the inflation sector {\it e.g.} non minimal coupling and inflaton VEV. To explore the dependence on inflationary parameters further, we also estimate the DM relic using the same mechanism for two other inflationary models consistent with latest data and observe that one of these models predicts different range of DM mass upto hundreds of TeV.


title: Weighing the Galactic disk in sub-regions of the solar neighbourhood using Gaia DR2

abstract: We infer the gravitational potential of the Galactic disk by analysing the phase-space densities of 120 stellar samples in 40 spatially separate sub-regions of the solar neighbourhood, using Gaia's second data release (DR2), in order to quantify spatially dependent systematic effects that bias this type of measurement. The gravitational potential was inferred under the assumption of a steady state in the framework of a Bayesian hierarchical model. We performed a joint fit of our stellar tracers' three-dimensional velocity distribution, while fully accounting for the astrometric uncertainties of all stars. The inferred gravitational potential is compared, post-inference, to a model for the baryonic matter and halo dark matter components. We see an unexpected but clear trend for all 40 spatially separate sub-regions: Compared to the potential derived from the baryonic model, the inferred gravitational potential is significantly steeper close to the Galactic mid-plane (<60 pc), but flattens such that the two agree well at greater distances (~400 pc). The inferred potential implies a total matter density distribution that is highly concentrated to the Galactic mid-plane and decays quickly with height. Apart from this, there are discrepancies between stellar samples, implying spatially dependent systematic effects which are, at least in part, explained by substructures in the phase-space distributions. In terms of the inferred matter density distribution, the very low matter density that is inferred at greater heights is inconsistent with the observed scale height and matter distribution of the stellar disk, which cannot be explained by a misunderstood density of cold gas or other hidden mass. Our interpretation is that these results must be biased by a time-varying phase-space structure, possibly a breathing mode, that is large enough to affect all stellar samples in the same manner.


title: General cosmological perturbations in teleparallel gravity

abstract: We discuss linear perturbations of the most general class of teleparallel spacetimes with cosmological symmetry, and perform a decomposition of these perturbations into irreducible components. We then study their behavior under gauge transformations, i.e., infinitesimal diffeomorphisms of the background spacetime, and derive a full set of gauge-invariant perturbation variables. Comparing our results with perturbations around the diagonal tetrad corresponding to a spatially flat teleparallel background cosmology, we see the appearance of new terms which depend on the particular choice of a cosmologically symmetric background tetrad. Finally, we demonstrate the application of our findings to teleparallel gravity theories. For convenience, results are displayed both in general and conformal coordinates.


title: Quantum Interference in Jet Substructure from Spinning Gluons

abstract: Collimated sprays of hadrons, called jets, are an emergent phenomenon of Quantum Chromodynamics (QCD) at collider experiments, whose detailed internal structure encodes valuable information about the interactions of high energy quarks and gluons, and their confinement into color-neutral hadrons. The flow of energy within jets is characterized by correlation functions of energy flow operators, with the three-point correlator, being the first correlator with non-trivial shape dependence, playing a special role in unravelling the dynamics of QCD. In this Letter we initiate a study of the three-point energy correlator to all orders in the strong coupling constant, in the limit where two of the detectors are squeezed together. We show that by rotating the two squeezed detectors with respect to the third by an angle $\phi$, a $\cos (2\phi)$ dependence arising from the quantum interference between intermediate virtual gluons with $+/-$ helicity is imprinted on the detector. This can be regarded as a double slit experiment performed with jet substructure, and it provides a direct probe of the ultimately quantum nature of the substructure of jets, and of transverse spin physics in QCD. To facilitate our all-orders analysis, we adopt the Operator Product Expansion (OPE) for light-ray operators in conformal field theory and develop it in QCD. Our application of the light-ray OPE in real world QCD establishes it as a powerful theoretical tool with broad applications for the study of jet substructure.


title: Dependence of Dark Matter - Electron Scattering on the Galactic Dark Matter Velocity Distribution

abstract: The rate of dark matter-electron scattering depends on the underlying velocity distribution of the dark matter halo. Importantly, dark matter-electron scattering is particularly sensitive to the high-velocity tail, which differs significantly amongst the various dark matter halo models. In this work, we summarize the leading halo models and discuss the various parameters which enter them. We recommend updated values for these parameters based on recent studies and measurements. Furthermore, we quantify the dependence of the dark matter-electron scattering rate on the choice of halo model and parameters, and demonstrate how these choices propagate to the predicted cross-section limits. The rate is most sensitive to changes in the circular velocity v0; in silicon targets, we find that the changes in the rate predictions can range from O(0.01%) to O(100%) for contact interactions and O(10%) to O(100%) for long-range interactions.


title: Two-Loop anomalous dimensions of QCD operators up to dimension-sixteen and Higgs EFT amplitudes

abstract: We consider two-loop renormalization of high-dimensional Lorentz scalar operators in the gluonic sector of QCD. These operators appear also in the Higgs effective theory obtained by integrating out the top quark loop in the gluon fusion process. We first discuss the classification of operators and how to construct a good set of basis using both off-shell field theory method and on-shell form factor formalism. To study loop corrections, we apply efficient unitarity-IBP strategy and compute the two-loop minimal form factors of length-3 operators up to dimension sixteen. From the UV divergences of form factor results, we extract the renormalization matrices and analyze the operator mixing behavior in detail. The form factors we compute are also equivalent to Higgs plus three-gluon amplitudes that capture high-order top mass corrections in Higgs EFT. We obtain the analytic finite remainder functions which exhibit several universal transcendentality structures.


title: Orbital Foregrounds for Ultra-Short Duration Transients

abstract: Reflections from objects in Earth orbit can produce sub-second, star-like optical flashes similar to astrophysical transients. Reflections have historically caused false alarms for transient surveys, but the population has not been systematically studied. We report event rates for these orbital flashes using the Evryscope Fast Transient Engine, a low-latency transient detection pipeline for the Evryscopes. We select single-epoch detections likely caused by Earth satellites and model the event rate as a function of both magnitude and sky position. We measure a rate of $1800^{+600}_{-280}$ sky$^{-1}$ hour$^{-1}$, peaking at $m_g = 13.0$, for flashes morphologically degenerate with real astrophysical signals in surveys like the Evryscopes. Of these, $340^{+150}_{-85}$ sky$^{-1}$ hour$^{-1}$ are bright enough to be visible to the naked eye in typical suburban skies with a visual limiting magnitude of $V\approx4$. These measurements place the event rate of orbital flashes orders of magnitude higher than the combined rate of public alerts from all active all-sky fast-timescale transient searches, including neutrino, gravitational-wave, gamma-ray, and radio observatories. Short-timescale orbital flashes form a dominating foreground for un-triggered searches for fast transients in low-resolution, wide-angle surveys. However, events like fast radio bursts (FRBs) with arcminute-scale localization have a low probability ($\sim10^{-5}$) of coincidence with an orbital flash, allowing optical surveys to place constraints on their potential optical counterparts in single images. Upcoming satellite internet constellations, like SpaceX Starlink, are unlikely to contribute significantly to the population of orbital flashes in normal operations.


title: Off-shell Noether currents and potentials for first-order general relativity

abstract: We report off-shell Noether currents obtained from off-shell Noether potentials for first-order general relativity described by $n$-dimensional Palatini and Holst Lagrangians including the cosmological constant. These off-shell currents and potentials are achieved by using the corresponding Lagrangian and the off-shell Noether identities satisfied by diffeomorphisms generated by arbitrary vector fields, local $SO(n)$ or $SO(n-1,1)$ transformations, `improved diffeomorphisms', and the `generalization of local translations' of the orthonormal frame and the connection. A remarkable aspect of our approach is that we do {\it not} use Noether's theorem in its direct form. By construction, the currents are off-shell conserved and lead naturally to the definition of off-shell Noether charges. We also study what we call the `half off-shell' case for both Palatini and Holst Lagrangians. In particular, we find that the resulting diffeomorphism and local $SO(3,1)$ or $SO(4)$ off-shell Noether currents and potentials for the Holst Lagrangian generically depend on the Immirzi parameter, which holds even in the `half off-shell' and on-shell cases. We also study Killing vector fields in the `half off-shell' and on-shell cases. The current theoretical framework is illustrated for the `half off-shell' case in static spherically symmetric and Friedmann--Lemaitre--Robertson--Walker spacetimes in four dimensions.


title: The impact of astrophysical dust grains on the confinement of cosmic rays

abstract: We argue that charged dust grains could significantly impact the confinement and transport of galactic cosmic rays. For sub-GeV to ~1000GeV cosmic rays, small-scale parallel Alfv\'en waves, which isotropize cosmic rays through gyro-resonant interactions, are also gyro-resonant with charged grains. If the dust is nearly stationary, as in the bulk of the interstellar medium, Alfv\'en waves are damped by dust. This will reduce the amplitude of Alfv\'en waves produced by the cosmic rays through the streaming instability, thus enhancing cosmic-ray transport. In well-ionized regions, the dust damping rate is larger by a factor of ~10 than other mechanisms that damp parallel Alfv\'en waves at the scales relevant for ~GeV cosmic rays, suggesting that dust could play a key role in regulating cosmic-ray transport. In astrophysical situations in which the dust moves through the gas with super-Alfv\'enic velocities, Alfv\'en waves are rendered unstable, which could directly scatter cosmic rays. This interaction has the potential to create a strong feedback mechanism where dust, driven through the gas by radiation pressure, then strongly enhances the confinement of cosmic rays, increasing their capacity to drive outflows. This mechanism may act in the circumgalactic medium around star-forming galaxies and active galactic nuclei.


title: SOFIA/FIFI-LS Full-disk [CII] Mapping and CO-dark Molecular Gas across the Nearby Spiral Galaxy NGC 6946

abstract: We present SOFIA/FIFI-LS observations of the [CII] 158${\mu}$m cooling line across the nearby spiral galaxy NGC 6946. We combine these with UV, IR, CO, and H I data to compare [CII] emission to dust properties, star formation rate (SFR), H$_2$, and HI at 560pc scales via stacking by environment (spiral arms, interarm, and center), radial profiles, and individual, beam-sized measurements. We attribute $73\%$ of the [CII] luminosity to arms, and $19\%$ and $8\%$ to the center and interarm region, respectively. [CII]/TIR, [CII]/CO, and [CII]/PAH radial profiles are largely constant, but rise at large radii ($\gtrsim$8kpc) and drop in the center ("[CII] deficit"). This increase at large radii and the observed decline with the 70${\mu}$m/100${\mu}$m dust color are likely driven by radiation field hardness. We find a near proportional [CII]-SFR scaling relation for beam-sized regions, though the exact scaling depends on methodology. [CII] also becomes increasingly luminous relative to CO at low SFR (interarm or large radii), likely indicating more efficient photodissociation of CO and emphasizing the importance of [CII] as an H$_2$ and SFR tracer in such regimes. Finally, based on the observed [CII] and CO radial profiles and different models, we find ${\alpha}_{CO}$ to increase with radius, in line with the observed metallicity gradient. The low ${\alpha}_{CO}$ (galaxy average $\lesssim2\,M_{sun}\,pc^{-2}\,(K\,km\,s^{-1})^{-1}$) and low [CII]/CO ratios ($\sim$400 on average) imply little CO-dark gas across NGC 6946, in contrast to estimates in the Milky Way.


title: When Freeze-out occurs due to a non-Boltzmann suppression: A study of degenerate dark sector

abstract: Exponential suppression or commonly known as the Boltzmann suppression in the number density of dark matter is the key ingredient for creating chemical imbalance prior to the usual thermal freeze-out. A degenerate/quasi-degenerate dark sector can experience a different exponential suppression in the number density analogous to the radioactive decay law leading to a delayed freeze-out mechanism of dark matter known as the co-decaying dark matter. In this work, we study the dynamics of a multicomponent dark matter from thermally decoupled degenerate dark sector in a hidden U$(1)_{X}$ extension of the Standard Model. We compute the relic density of dark matter frozen-out through the co-decaying mechanism by solving four coupled Boltzmann equations. We demonstrate how temperature $T^\prime $ of the dark sector changes due to all types of $3\rightarrow 2$ and $2\rightarrow 2$ interactions along with the eternal expansion of the Universe. We find that $3\rightarrow 2$ interactions enhance $T^\prime$ by producing energetic particles in the dark sector while the excess heat is transferred by $2\rightarrow 2$ interactions to the entire dark sector. As the direct detection is possible only through the feeble portal couplings, we investigate the neutrino and $\gamma$-ray signals from dark matter annihilation via one step cascade processes and compare our results with the measured fluxes of atmospheric neutrinos by Super-Kamiokande and diffuse $\gamma$-rays by Fermi-LAT, EGRET, INTEGRAL collaborations. We find that the present scenario easily evades all the existing bounds from atmospheric neutrino and diffuse $\gamma$-ray observations for degenerate dark sector. However, the constraints are significant for quasi degenerate scenario.


title: Highly Tunable Junctions and Nonlocal Josephson Effect in Magic Angle Graphene Tunneling Devices

abstract: Magic-angle twisted bilayer graphene (MATBG) has recently emerged as a highly tunable two-dimensional (2D) material platform exhibiting a wide range of phases, such as metal, insulator, and superconductor states. Local electrostatic control over these phases may enable the creation of versatile quantum devices that were previously not achievable in other single material platforms. Here, we exploit the electrical tunability of MATBG to engineer Josephson junctions and tunneling transistors all within one material, defined solely by electrostatic gates. Our multi-gated device geometry offers complete control over the Josephson junction, with the ability to independently tune the weak link, barriers, and tunneling electrodes. We show that these purely 2D MATBG Josephson junctions exhibit nonlocal electrodynamics in a magnetic field, in agreement with the Pearl theory for ultrathin superconductors. Utilizing the intrinsic bandgaps of MATBG, we also demonstrate monolithic edge tunneling spectroscopy within the same MATBG devices and measure the energy spectrum of MATBG in the superconducting phase. Furthermore, by inducing a double barrier geometry, the devices can be operated as a single-electron transistor, exhibiting Coulomb blockade. These MATBG tunneling devices, with versatile functionality encompassed within a single material, may find applications in graphene-based tunable superconducting qubits, on-chip superconducting circuits, and electromagnetic sensing in next-generation quantum nanoelectronics.


title: Bringing faint active galactic nuclei (AGNs) to light: a view from large-scale cosmological simulations

abstract: The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyze the population of faint AGN (L_x (2-10 keV) < 10^42 erg/s) in the Illustris, TNG100, EAGLE, and SIMBA cosmological simulations, and find that the properties of their host galaxies vary from one simulation to another. In Illustris and EAGLE, faint AGN are powered by low-mass BHs located in low-mass star-forming galaxies. In TNG100 and SIMBA, they are mostly associated with more massive BHs in quenched massive galaxies. We model the X-ray binary populations (XRB) of the simulated galaxies, and find that AGN often dominate the galaxy AGN+XRB hard X-ray luminosity at z>2, while XRBs dominate in some simulations at z<2. Whether the AGN or XRB emission dominates in star-forming and quenched galaxies depends on the simulations. These differences in simulations can be used to discriminate between galaxy formation models with future high-resolution X-ray observations. We compare the luminosity of simulated faint AGN host galaxies to observations of stacked galaxies from Chandra. Our comparison indicates that the simulations post-processed with our X-ray modeling tend to overestimate the AGN+XRB X-ray luminosity; luminosity that can be strongly affected by AGN obscuration. Some simulations reveal clear AGN trends as a function of stellar mass (e.g., galaxy luminosity drop in massive galaxies), which are not apparent in the observations.


title: Automated approach to measure stellar inclinations: validation through large-scale measurements on the red giant branch

abstract: Measuring stellar inclinations is fundamental to understand planetary formation and dynamics as well as physical conditions during star formation. Oscillation spectra of red giant stars exhibit mixed modes that have both a gravity component from the radiative interior and a pressure component from the convective envelope. Gravity-dominated (g-m) mixed modes split by rotation are well separated inside frequency spectra, making possible accurate measurements of stellar inclinations. This work aims at developing an automated and general approach to measure stellar inclinations, that can be applied to any solar-type pulsator for which oscillation modes are identified, and at validating it using red giant branch stars observed by Kepler. We use the mean height-to-background ratio of dipole mixed modes with different azimuthal orders to measure stellar inclinations. The underlying statistical distribution of inclinations is recovered in an unbiased way using a probability density function for the stellar inclination angle. We derive stellar inclination measurements for 1139 stars on the red giant branch, for which Gehan et al. (2018) have identified the azimuthal order of dipole g-m mixed modes. Raw measured inclinations exhibit strong deviation with respect to isotropy which is expected for random inclinations over the sky. When taking uncertainties into account, the reconstructed distribution of inclinations actually follows the expected isotropic distribution of the rotational axis. This work highlights the biases that affect inclination measurements and provides the way to infer their underlying statistical distribution. When the star is seen either pole-on or equator-on, measurements are challenging and result in a biased distribution. Correcting biases that appear at the low- and high inclination regimes allows us to recover the underlying inclination distribution.


title: Two-Gap Time Reversal Symmetry Breaking Superconductivity in Non-Centrosymmetric LaNiC2

abstract: We report a $\mu$SR investigation of a non-centrosymmetric superconductor (LaNiC$_2$) in single crystal form. Compared to previous $\mu$SR studies of non-centrosymmetric superconducting polycrystalline and powder samples, the unambiguous orientation of single crystals enables a simultaneous determination of the absolute value of the magnetic penetration depth and the vortex core size from measurements that probe the magnetic field distribution in the vortex state. The magnetic field dependence of these quantities unambiguously demonstrates the presence of two nodeless superconducting energy gaps. In addition, we detect weak internal magnetic fields in the superconducting phase, confirming earlier $\mu$SR evidence for a time-reversal symmetry breaking superconducting state. Our results suggest that Cooper pairing in LaNiC$_2$ is characterized by the same interorbital equal-spin pairing model introduced to describe the pairing state in the centrosymmetric superconductor LaNiGa$_2$.


title: Structural Properties and $\alpha$-Decay Chains of Transfermium Nuclei (101$\leq$Z$\leq$110)

abstract: Transfermium nuclei (101$\leq$Z$\leq$110) are investigated thoroughly to describe structural properties viz. deformation, radii, shapes, magicity, etc. as well as their probable decay chains. These properties are explored using relativistic mean-field (RMF) approach and compared with other theories along with available experimental data. Neutron numbers N$=$152 and 162 have come forth with a deformed shell gap whereas N$=$184 is ensured as a spherical magic number. The region with N$>$168 bears witness of the phenomenon of shape transition and shape coexistence for all the considered isotopic chains. Experimental $\alpha$-decay half-lives are compared with our theoretical half-lives obtained by using various empirical/semi-empirical formulas. The recent formula proposed by Manjunatha \textit{et al.}, which results best among the considered 10 formulas, is further modified by adding asymmetry dependent terms ($I$ and $I^2$). This modified Manjunatha formula is utilized to predict probable $\alpha$-decay chains that are found in excellent agreement with available experimental data.


title: The GALAH+ Survey: Third Data Release

abstract: The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spectrograph at the Anglo-Australian Telescope. This release (hereafter GALAH+ DR3) includes all observations from GALAH Phase 1 (bright, main, and faint survey, 70%), K2-HERMES (17%), TESS-HERMES (5%), and a subset of ancillary observations (8%) including the bulge and >75 stellar clusters. We derive stellar parameters $T_\text{eff}$, $\log g$, [Fe/H], $v_\text{mic}$, $v_\text{broad}$ & $v_\text{rad}$ using our modified version of the spectrum synthesis code Spectroscopy Made Easy (SME) and 1D MARCS model atmospheres. We break spectroscopic degeneracies in our spectrum analysis with astrometry from $Gaia$ DR2 and photometry from 2MASS. We report abundance ratios [X/Fe] for 30 different elements (11 of which are based on non-LTE computations) covering five nucleosynthetic pathways. We describe validations for accuracy and precision, flagging of peculiar stars/measurements and recommendations for using our results. Our catalogue comprises 65% dwarfs, 34% giants, and 1% other/unclassified stars. Based on unflagged chemical composition and age, we find 62% young low-$\alpha$, 9% young high-$\alpha$, 27% old high-$\alpha$, and 2% stars with $\mathrm{[Fe/H]} \leq -1$. Based on kinematics, 4% are halo stars. Several Value-Added-Catalogues, including stellar ages and dynamics, updated after $Gaia$ eDR3, accompany this release and allow chrono-chemodynamic analyses, as we showcase.


title: The dynamic effect of mechanical losses of actuators on the equations of motion of legged robots

abstract: Industrial manipulators do not collapse under their own weight when powered off due to the friction in their joints. Although these mechanism are effective for stiff position control of pick-and-place, they are inappropriate for legged robots which must rapidly regulate compliant interactions with the environment. However, no metric exists to quantify the robot's perform degradation due to mechanical losses in the actuators. This letter provides a novel formulation which describes how the efficiency of individual actuators propagate to the equations of motion of the whole robot. We quantitatively demonstrate the intuitive fact that the apparent inertia of the robots increase in the presence of joint friction. We also reproduce the empirical result that robots which employ high gearing and low efficiency actuators can statically sustain more substantial external loads. We expect that the framework presented here will provide the foundations to design the next generation of legged robots which can effectively interact with the world.


title: High Eccentricities and High Masses Characterize Gravitational-wave Captures in Galactic Nuclei as Seen by Earth-based Detectors

abstract: The emission of gravitational waves (GWs) during single-single close encounters in galactic nuclei (GNs) leads to the formation and rapid merger of highly eccentric stellar-mass black hole (BH) binaries. The distinct distribution of physical parameters makes it possible to statistically distinguish this source population from others. Previous studies determined the expected binary parameter distribution for this source population in single GNs. Here we take into account the effects of dynamical friction, post-Newtonian corrections, and observational bias to determine the detected sources' parameter-distributions from all GNs in the Universe. We find that the total binary mass distribution of detected mergers is strongly tilted towards higher masses. The distribution of initial peak GW frequency is remarkably high between 1-70 Hz, ~50% of GW capture sources form above 10 Hz with e >~ 0.95. The eccentricity when first entering the LIGO/Virgo/KAGRA band satisfies e_10Hz > 0.1 for over 92% of sources and e_10Hz > 0.8 for more than half of the sources. At the point when the pericenter reaches 10 GM/c^2 the eccentricity satisfies e_10M > 0.1 for over ~70% of the sources, making single-single GW capture events in GNs the most eccentric source population among the currently known stellar-mass binary BH merger channels in the Universe. We identify correlations between total mass, mass ratio, source detection distance, and eccentricities e_10Hz and e_10M. The recently measured source parameters of GW190521 lie close to the peak of the theoretical distributions and the estimated escape speed of the host environment is ~7.5x10^3 km/s - 1.2x10^4 km/s, making this source a candidate for this astrophysical merger channel.


title: A Comparison Between Joint Space and Task Space Mappings for Dynamic Teleoperation of an Anthropomorphic Robotic Arm in Reaction Tests

abstract: Teleoperation (i.e., controlling a robot with human motion) proves promising in enabling a humanoid robot to move as dynamically as a human. But how to map human motion to a humanoid robot matters because a human and a humanoid robot rarely have identical topologies and dimensions. This work presents an experimental study that utilizes reaction tests to compare the proposed joint space mapping and the proposed task space mapping for dynamic teleoperation of an anthropomorphic robotic arm that possesses human-level dynamic motion capabilities. The experimental results suggest that the robot achieved similar and, in some cases, human-level dynamic performances with both mappings for the six participating human subjects. All subjects became proficient at teleoperating the robot with both mappings after practice, despite that the subjects and the robot differed in size and link length ratio and that the teleoperation required the subjects to move unintuitively. Yet, most subjects developed their teleoperation proficiencies more quickly with the task space mapping than with the joint space mapping after similar amounts of practice. This study also indicates the potential values of a three-dimensional task space mapping, a teleoperation training simulator, and force feedback to the human pilot for intuitive and dynamic teleoperation of a humanoid robot's arms.


title: Trading off 1-norm and sparsity against rank for linear models using mathematical optimization: 1-norm minimizing partially reflexive ah-symmetric generalized inverses

abstract: The M-P (Moore-Penrose) pseudoinverse has as a key application the computation of least-squares solutions of inconsistent systems of linear equations. Irrespective of whether a given input matrix is sparse, its M-P pseudoinverse can be dense, potentially leading to high computational burden, especially when we are dealing with high-dimensional matrices. The M-P pseudoinverse is uniquely characterized by four properties, but only two of them need to be satisfied for the computation of least-squares solutions. Fampa and Lee (2018) and Xu, Fampa, Lee, and Ponte (2019) propose local-search procedures to construct sparse block-structured generalized inverses that satisfy the two key M-P properties, plus one more (the so-called reflexive property). That additional M-P property is equivalent to imposing a minimum-rank condition on the generalized inverse. (Vector) 1-norm minimization is used to induce sparsity and, importantly, to keep the magnitudes of entries under control for the generalized-inverses constructed. Here, we investigate the trade-off between low 1-norm and low rank for generalized inverses that can be used in the computation of least-squares solutions. We propose several algorithmic approaches that start from a $1$-norm minimizing generalized inverse that satisfies the two key M-P properties, and gradually decrease its rank, by iteratively imposing the reflexive property. The algorithms iterate until the generalized inverse has the least possible rank. During the iterations, we produce intermediate solutions, trading off low 1-norm (and typically high sparsity) against low rank.


title: Gravitational dark matter production: primordial black holes and UV freeze-in

abstract: Dark matter (DM) interacting only gravitationally with the standard model could have been produced in the early universe by Hawking evaporation of primordial black holes (PBH). This mechanism is viable in a large range of DM mass, spanning up to the Planck scale. However, DM is also unavoidably produced by the irreducible UV gravitational freeze-in. We show that the latter mechanism sets strong bounds, excluding large regions of the parameter space favored by PBH production.


title: Offline Reinforcement Learning from Human Feedback in Real-World Sequence-to-Sequence Tasks

abstract: Large volumes of interaction logs can be collected from NLP systems that are deployed in the real world. How can this wealth of information be leveraged? Using such interaction logs in an offline reinforcement learning (RL) setting is a promising approach. However, due to the nature of NLP tasks and the constraints of production systems, a series of challenges arise. We present a concise overview of these challenges and discuss possible solutions.


title: Robust quantum gates using smooth pulses and physics-informed neural networks

abstract: The presence of decoherence in quantum computers necessitates the suppression of noise. Dynamically corrected gates via specially designed control pulses offer a path forward, but hardware-specific experimental constraints can cause complications. Existing methods to obtain smooth pulses are either restricted to two-level systems, require an optimization over noise realizations or limited to piecewise-continuous pulse sequences. In this work, we present the first general method for obtaining truly smooth pulses that minimizes sensitivity to noise, eliminating the need for sampling over noise realizations and making assumptions regarding the underlying statistics of the experimental noise. We parametrize the Hamiltonian using a neural network, which allows the use of a large number of optimization parameters to adequately explore the functional control space. We demonstrate the capability of our approach by finding smooth shapes which suppress the effects of noise within the logical subspace as well as leakage out of that subspace.


title: Chern-Simons Gravity Dual of BCFT

abstract: In this paper we provide a Chern-Simons gravity dual of a two dimensional conformal field theory on a manifold with boundaries, so called boundary conformal field theory (BCFT). We determine the correct boundary action on the end of the world brane in the Chern-Simons gauge theory. This reproduces known results of the AdS/BCFT for the Einstein gravity. We also give a prescription of calculating holographic entanglement entropy by employing Wilson lines which extend from the AdS boundary to the end of the world brane. We also discuss a higher spin extension of our formulation.


title: Monitoring the Impact of Wildfires on Tree Species with Deep Learning

abstract: One of the impacts of climate change is the difficulty of tree regrowth after wildfires over areas that traditionally were covered by certain tree species. Here a deep learning model is customized to classify land covers from four-band aerial imagery before and after wildfires to study the prolonged consequences of wildfires on tree species. The tree species labels are generated from manually delineated maps for five land cover classes: Conifer, Hardwood, Shrub, ReforestedTree and Barren land. With an accuracy of $92\%$ on the test split, the model is applied to three wildfires on data from 2009 to 2018. The model accurately delineates areas damaged by wildfires, changes in tree species and rebound of burned areas. The result shows clear evidence of wildfires impacting the local ecosystem and the outlined approach can help monitor reforested areas, observe changes in forest composition and track wildfire impact on tree species.


title: Magnetic and Electric Black Holes in the Vector-Tensor Horndeski Theory

abstract: We construct exact solutions of magnetically charged black holes in the vector-tensor Horndeski gravity and discuss their main features. Unlike the analogous electric case, the field equations are linear in a simple (quite standard) parametrization of the metric tensor and they can be solved analytically even when a cosmological constant is added. The solutions are presented in terms of hypergeometric functions which makes the analysis of the black hole properties relatively straightforward. Some of the aspects of these black holes are quite ordinary like the existence of extremal configurations with maximal magnetic charge for a given mass, or the existence of a mass with maximal temperature for a given charge, but others are somewhat unexpected, like the existence of black holes with a repulsive gravitational field. We perform our analysis for both signs of the non-minimal coupling constant and find black hole solutions in both cases but with significant differences between them. The most prominent difference is the fact that the black holes for the negative coupling constant have a spherical surface of curvature singularity rather than a single point. On the other hand, the gravitational field produced around this kind of black holes is always attractive. Also, for small enough magnetic charge and negative coupling constant, extremal black holes do not exist and all magnetic black holes have a single horizon. In addition we study the trajectories around these magnetic black holes for light as well as massive particles either neutral or electrically charged. Finally, we compare the main features of these black holes with their electric counterparts, adding some aspects that have not been discussed before, like temperature, particle trajectories and light deflection by electrically charged Horndesky black holes.


title: On Key Epidemiological Metrics during Infectious Disease Outbreaks

abstract: In this work, we review the figures used to characterize an epidemic outbreak most. Particular attention is drawn to epidemic spreading at time-varying transition rates. A time-varying SIR-like model is used to describe the epidemic dynamics, from which closed equations relating parameters and key quantities like reproduction ratio and doubling time are derived. The definition and computation of these metrics are revisited in the context of the general solution to the time-varying model dynamics, focusing on the similarities and differences with the time-invariant case. Further, the prediction of these metrics, that is of the disease evolution, as response to different scenarios is also investigated.


title: The GALAH survey: effective temperature calibration from the InfraRed Flux Method in the Gaia system

abstract: In order to accurately determine stellar properties, knowledge of the effective temperature of stars is vital. We implement Gaia and 2MASS photometry in the InfraRed Flux Method and apply it to over 360,000 stars across different evolutionary stages in the GALAH DR3 survey. We derive colour-effective temperature relations that take into account the effect of metallicity and surface gravity over the range 4000 to 8000 kelvin, from very metal-poor stars to super solar metallicities. The internal uncertainty of these calibrations is of order 40-80 kelvin depending on the colour combination used. Comparison against solar-twins, Gaia benchmark stars and the latest interferometric measurements validates the precision and accuracy of these calibrations from F to early M spectral types. We assess the impact of various sources of uncertainties, including the assumed extinction law, and provide guidelines to use our relations. Robust solar colours are also derived.


title: Extremum Seeking Control with Attenuated Steady-State Oscillations

abstract: We propose two perturbation-based extremum seeking control (ESC) schemes for general single input single output nonlinear dynamical systems, having structures similar to that of the classical ESC scheme. We propose novel adaptation laws for the excitation signal amplitudes in each scheme that drive the amplitudes to zero. The rates of decay for both the laws are governed by the gradient measures of the unknown reference-to-output equilibrium map. We show that the proposed ESC schemes achieve practical asymptotic convergence to the extremum with a proper tuning of the parameters in the proposed schemes. As the extremum is reached, and the magnitudes of the gradient measures become small, the excitation signal amplitudes converge to zero. Thus, the proposed schemes ensure that the excitation signal is attenuated as the system output converges to a neighborhood of the extremum and the steady-state oscillations about the extremum, typically observed for the classical ESC schemes, are attenuated. Simulation examples are included to illustrate the effectiveness of the proposed schemes.


title: Physics-Informed Neural Network Super Resolution for Advection-Diffusion Models

abstract: Physics-informed neural networks (NN) are an emerging technique to improve spatial resolution and enforce physical consistency of data from physics models or satellite observations. A super-resolution (SR) technique is explored to reconstruct high-resolution images ($4\times$) from lower resolution images in an advection-diffusion model of atmospheric pollution plumes. SR performance is generally increased when the advection-diffusion equation constrains the NN in addition to conventional pixel-based constraints. The ability of SR techniques to also reconstruct missing data is investigated by randomly removing image pixels from the simulations and allowing the system to learn the content of missing data. Improvements in S/N of $11\%$ are demonstrated when physics equations are included in SR with $40\%$ pixel loss. Physics-informed NNs accurately reconstruct corrupted images and generate better results compared to the standard SR approaches.


title: Level attraction and exceptional points in a resonant spin-orbit torque system

abstract: Level attraction can appear in driven systems where instead of repulsion two modes coalesce in a region separated by two exceptional points. This behavior was proposed for optomechanical and optomagnonic systems, and recently observed for dissipative cavity magnon-polaritons. We demonstrate that such a regime exists in a spin-orbit torque system where a magnetic oscillator is resonantly coupled to an electron reservoir. An instability mechanism necessary for mode attraction can be provided by applying an electric field. The field excites interband transitions between spin-orbit split bands leading to an instability of the magnetic oscillator. Two exceptional points then appear in the oscillator energy spectrum and the region of instability. We discuss conditions under which this can occur and estimate the electric field strength necessary for reaching the attraction region for a spin-orbit torque oscillator with Rashba coupling. A proposal for experimental detection is made using magnetic susceptibility measurements.


title: Filtering for Aggregate Hidden Markov Models with Continuous Observations

abstract: We consider a class of filtering problems for large populations where each individual is modeled by the same hidden Markov model (HMM). In this paper, we focus on aggregate inference problems in HMMs with discrete state space and continuous observation space. The continuous observations are aggregated in a way such that the individuals are indistinguishable from measurements. We propose an aggregate inference algorithm called continuous observation collective forward-backward algorithm. It extends the recently proposed collective forward-backward algorithm for aggregate inference in HMMs with discrete observations to the case of continuous observations. The efficacy of this algorithm is illustrated through several numerical experiments.

categories: stat.ML cs.IT cs.LG cs.SY eess.SY math.IT math.OC


title: Gradient-Based Empirical Risk Minimization using Local Polynomial Regression

abstract: In this paper, we consider the problem of empirical risk minimization (ERM) of smooth, strongly convex loss functions using iterative gradient-based methods. A major goal of this literature has been to compare different algorithms, such as gradient descent (GD) or stochastic gradient descent (SGD), by analyzing their rates of convergence to $\epsilon$-approximate solutions. For example, the oracle complexity of GD is $O(n\log(\epsilon^{-1}))$, where $n$ is the number of training samples. When $n$ is large, this can be expensive in practice, and SGD is preferred due to its oracle complexity of $O(\epsilon^{-1})$. Such standard analyses only utilize the smoothness of the loss function in the parameter being optimized. In contrast, we demonstrate that when the loss function is smooth in the data, we can learn the oracle at every iteration and beat the oracle complexities of both GD and SGD in important regimes. Specifically, at every iteration, our proposed algorithm performs local polynomial regression to learn the gradient of the loss function, and then estimates the true gradient of the ERM objective function. We establish that the oracle complexity of our algorithm scales like $\tilde{O}((p \epsilon^{-1})^{d/(2\eta)})$ (neglecting sub-dominant factors), where $d$ and $p$ are the data and parameter space dimensions, respectively, and the gradient of the loss function belongs to a $\eta$-H\"{o}lder class with respect to the data. Our proof extends the analysis of local polynomial regression in non-parametric statistics to provide interpolation guarantees in multivariate settings, and also exploits tools from the inexact GD literature. Unlike GD and SGD, the complexity of our method depends on $d$ and $p$. However, when $d$ is small and the loss function exhibits modest smoothness in the data, our algorithm beats GD and SGD in oracle complexity for a very broad range of $p$ and $\epsilon$.


title: Hypersim: A Photorealistic Synthetic Dataset for Holistic Indoor Scene Understanding

abstract: For many fundamental scene understanding tasks, it is difficult or impossible to obtain per-pixel ground truth labels from real images. We address this challenge by introducing Hypersim, a photorealistic synthetic dataset for holistic indoor scene understanding. To create our dataset, we leverage a large repository of synthetic scenes created by professional artists, and we generate 77,400 images of 461 indoor scenes with detailed per-pixel labels and corresponding ground truth geometry. Our dataset: (1) relies exclusively on publicly available 3D assets; (2) includes complete scene geometry, material information, and lighting information for every scene; (3) includes dense per-pixel semantic instance segmentations and complete camera information for every image; and (4) factors every image into diffuse reflectance, diffuse illumination, and a non-diffuse residual term that captures view-dependent lighting effects. We analyze our dataset at the level of scenes, objects, and pixels, and we analyze costs in terms of money, computation time, and annotation effort. Remarkably, we find that it is possible to generate our entire dataset from scratch, for roughly half the cost of training a popular open-source natural language processing model. We also evaluate sim-to-real transfer performance on two real-world scene understanding tasks - semantic segmentation and 3D shape prediction - where we find that pre-training on our dataset significantly improves performance on both tasks, and achieves state-of-the-art performance on the most challenging Pix3D test set. All of our rendered image data, as well as all the code we used to generate our dataset and perform our experiments, is available online.


title: Study of General Relativity and $f(R)$ Modified Gravity with Cosmological Constant for Einstein and Jordan Frames

abstract: This thesis investigates a toy model for inflation in a class of modified theories of gravity in the metric formalism. Instead of the standard procedure -- assuming a non-linear Lagrangian $f(R)$ in the Jordan frame -- we start from a simple {\phi}^2 potential in the Einstein frame and investigate the corresponding $f(R)$ in the former picture. Such approach yields plenty of new pieces of information, namely a self-terminating inflationary solution with a linear Lagrangian, a robust criterion for stability of such theories, a dynamical effective potential for the Ricci scalar $R$, the addition of an ad-hoc Cosmological Constant in the Einstein frame leads to a Thermodynamical interpretation of this physical system, which allows further insight on its (meta)stability and evolution.


title: Deciphering the Kinematic Structure of the Small Magellanic Cloud through its Red Giant Population

abstract: We present a new kinematic model for the Small Magellanic Cloud (SMC), using data from the \gaia\ Data Release 2 catalog. We identify a sample of astrometrically well-behaved red giant (RG) stars belonging to the SMC and cross-match with publicly available radial velocity (RV) catalogs. We create a 3D spatial model for the RGs, using RR Lyrae for distance distributions, and apply kinematic models with varying rotation properties and a novel tidal expansion prescription to generate mock proper motion (PM) catalogs. When we compare this series of mock catalogs to the observed RG data, we find a combination of moderate rotation (with a magnitude of $\sim10-20$ km s$^{-1}$ at 1 kpc from the SMC center, inclination between $\sim50-80$ degrees, and a predominantly north-to-south line of nodes position angle of $\sim180$ degrees) and tidal expansion (with a scaling of $\sim10$ km s$^{-1}$ kpc$^{-1}$) is required to explain the PM signatures. The exact best-fit parameters depend somewhat on whether we assess only the PMs or include the RVs as a qualitative check, leaving some small tension remaining between the PM and RV conclusions. In either case, the parameter space preferred by our model is different both from previously inferred rotational geometries, including from the SMC H{\small I} gas and from the RG RV-only analyses, and new SMC PM analyses which conclude that a rotation signature is not detectable. Taken together this underscores the need to treat the SMC as a series of different populations with distinct kinematics.


title: Detecting Backdoors in Neural Networks Using Novel Feature-Based Anomaly Detection

abstract: This paper proposes a new defense against neural network backdooring attacks that are maliciously trained to mispredict in the presence of attacker-chosen triggers. Our defense is based on the intuition that the feature extraction layers of a backdoored network embed new features to detect the presence of a trigger and the subsequent classification layers learn to mispredict when triggers are detected. Therefore, to detect backdoors, the proposed defense uses two synergistic anomaly detectors trained on clean validation data: the first is a novelty detector that checks for anomalous features, while the second detects anomalous mappings from features to outputs by comparing with a separate classifier trained on validation data. The approach is evaluated on a wide range of backdoored networks (with multiple variations of triggers) that successfully evade state-of-the-art defenses. Additionally, we evaluate the robustness of our approach on imperceptible perturbations, scalability on large-scale datasets, and effectiveness under domain shift. This paper also shows that the defense can be further improved using data augmentation.


title: Graphene-assisted metagrating: from coherent to angular-asymmetric control of absorption and reflection

abstract: In this paper, we exploit the metagrating paradigm to achieve coherent control of absorption and reflection in a two-port device. Employing graphene ribbon as a tunable element allows us to, for the first time, realize a reconfigurable metagrating that integrates diversified coherent functionalities into a single planar structure. It is illustrated that the suggested design can behave as a coherent perfect absorber at multiple operating incident angles offered by its period-reconfigurability. Besides, our proposed metagrating is also capable of highly-efficient dynamic beam steering based on the coherent interaction of light with light through finely adjusting its chemical potentials, and its compatibility is then investigated to realize linear all-optical logic gates. Moreover, we use the idea of graphene-based metagrating to put forward an extremely asymmetrical device which exhibits high retroreflection upon illumination from one side and a strong absorption under excitation from the opposite side, revealing a great enhancement in design and fabrication simplicity compared to the previous works. The EM response of our suggested device can also be switched from a strongly asymmetric to a symmetric behavior by merely modulating the chemical potentials of its graphene ribbons.


title: Self-motility of an active particle induced by correlations in the surrounding solution

abstract: Current models of phoretic transport rely on molecular forces creating a "diffuse" particle-fluid interface. We investigate theoretically an alternative mechanism, in which a diffuse interface emerges solely due to a non-vanishing correlation length of the surrounding solution. This mechanism can drive self-motility of a chemically active particle. Numerical estimates indicate that the velocity can reach micrometers per second. The predicted phenomenology includes a bilinear dependence of the velocity on the activity and a possible double velocity reversal upon varying the correlation length.


title: CP violating effects in $^{210}$Fr and prospects for new physics beyond the Standard Model

abstract: We report theoretical results of the electric dipole moment (EDM) of $^{210}$Fr which arises from the interaction of the EDM of an electron with the internal electric field in an atom and the scalar-pseudoscalar electron-nucleus interaction; the two dominant sources of CP violation in this atom. Employing the relativistic coupled-cluster theory, we evaluate the enhancement factors for these two CP violating interactions to an accuracy of about 3% and analyze the contributions of the many-body effects. These two quantities in combination with the projected sensitivity of the $^{210}$Fr EDM experiment provide constraints on new physics beyond the Standard Model. Particularly, we demonstrate that their precise values are necessary to account for the effect of the bottom quark in models in which the Higgs sector is augmented by nonstandard Yukawa interactions such as the two-Higgs doublet model.


title: A hundred new eclipsing binary system candidates studied in a near-infrared window in the VVV survey

abstract: We present the first results obtained from an extensive study of eclipsing binary (EB) system candidates recently detected in the VISTA Variables in the V\'ia L\'actea (VVV) near-infrated (NIR) Survey. We analyze the VVV tile d040 in the southern part of the Galactic disc wherein the interstellar reddening is comparatively low, which makes it possible to detect hundreds of new eclipsing binary candidates. We present here the light curves and the determination of the geometric and physical parameters of the best candidates found in this "NIR window", including 37 contact, 50 detached and 13 semi-detached eclipsing binary systems. We infer that the studied systems have an average of the $K_s$ amplitudes of $0.8$ mag and a median period of 1.22 days where, in general, contact binaries have shorter periods. Using the "Physics Of Eclipsing Binaries" (PHOEBE) interactive interface, which is based on the Wilson and Devinney code, we find that the studied systems have low eccentricities. The studied EBs present mean values of about 5700 K and 4900 K for the $T_1$ and $T_2$ components, respectively. The mean mass-ratio ($q$) for the contact EB stars is $\sim$ 0.44. This new Galactic disc sample is a first approach to the massive study of NIR EB systems.


title: Weinberg operator contribution to the nucleon electric dipole moment in the quark model

abstract: We evaluate the contribution of the CP violating gluon chromo-electric dipole moment (the so-called Weinberg operator, denoted as $w$) to the electric dipole moment (EDM) of nucleons in the nonrelativistic quark model. The CP-odd interquark potential is modeled by the perturbative one-loop level gluon exchange generated by the Weinberg operator with massive quarks and gluons. The nucleon EDM is obtained by solving the nonrelativistic Schr\"{o}dinger equation of the three-quark system using the Gaussian expansion method. It is found that the resulting nucleon EDM, which may reasonably be considered as the irreducible contribution, is smaller than the one obtained after $\gamma_5$-rotating the anomalous magnetic moment using the CP-odd mass calculated with QCD sum rules. We estimate the total contribution to be $d_n = w \times 20 \, e \, {\rm MeV}$ and $d_p = - w \times 18 \, e \, {\rm MeV}$ with 60% of theoretical uncertainty.


title: Electron quasi-itinerancy intertwined with quantum order by disorder in pyrochlore iridate magnetism

abstract: We point out the emergence of magnetism from the interplay of electron quasi-itinerancy and quantum order by disorder in pyrochlore iridates. Like other Mott insulating iridates, the Ir$^{4+}$ ion in pyrochlore iridates develops an effective ${J=1/2}$ moment from the on-site spin-orbit coupling. We consider the generic symmetry-allowed exchange between these local moments on a pyrochlore lattice and obtain the mean-field phase diagram. Assuming the superexchange is mediated by direct and/or indirect electron hopping via intermediate oxygens, we derive the exchange interactions in the strong coupling regime from the Hubbard model. This exchange has a degenerate classical ground state manifold and quantum fluctuation selects a non-coplanar ground state, known as quantum order by disorder. Extending to the intermediate coupling regime, the same non-coplanar order is selected from the degenerate manifold by the kinetic energy, which is dubbed "electron quasi-itinerancy". We discuss the experimental relevance of our results and electron quasi-itinerancy among other iridates and $4d/5d$ magnets.


title: The GALAH survey: tracing the Galactic disk with Open Clusters

abstract: Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time.


title: A tunable population timer in multicellular consortia

abstract: Processing time-dependent information requires cells to quantify the duration of past regulatory events and program the time span of future signals. At the single-cell level, timer mechanisms can be implemented with genetic circuits: sets of genes connected to achieve specific tasks. However, such systems are difficult to implement in single cells due to saturation in molecular components and stochasticity in the limited intracellular space. Multicellular implementations, on the other hand, outsource some of the components of information-processing circuits to the extracellular space, and thereby might escape those constraints. Here we develop a theoretical framework, based on a trilinear coordinate representation, to study the collective behavior of a cellular population composed of three cell types under stationary conditions. This framework reveals that distributing different processes (in our case the production, detection and degradation of a time-encoding signal) across distinct cell types enables the robust implementation of a multicellular timer. Our analysis also shows the circuit to be easily tunable by varying the cellular composition of the consortium.


title: How far do Activated Random Walkers spread from a single source?

abstract: Unlike many particle systems, Activated Random Walk has nontrivial behavior even in one spatial dimension. We prove inner and outer bounds on the spread of n activated random walkers from a single source in Z. The inner bound involves a comparison with the stationary distribution of activated random walkers on a finite interval, while the outer bound involves a comparison with the stabilization of an infinite Bernoulli configuration of activated random walkers on Z.


title: The impact of pre-main sequence stellar evolution on midplane snowline locations and C/O in planet forming discs

abstract: We investigate the impact of pre-main sequence stellar luminosity evolution on the thermal and chemical properties of disc midplanes. We create template disc models exemplifying initial conditions for giant planet formation for a variety of stellar masses and ages. These models include the 2D physical structure of gas as well as 1D chemical structure in the disc midplane. The disc temperature profiles are calculated using fully physically consistent radiative transfer models for stars between 0.5 and 3 Msun and ages up to 10 Myr. The resulting temperature profiles are used to determine how the chemical conditions in the mid-plane change over time. We therefore obtain gas and ice-phase abundances of the main carbon and oxygen carrier species. While the temperature profiles produced are not markedly different for the stars of different masses at early stages (<1 Myr), they start to diverge significantly beyond 2 Myr. Discs around stars with mass >1.5 Msun become warmer over time as the stellar luminosity increases, whereas low-mass stars decrease in luminosity leading to cooler discs. This has an observable effect on the location of the CO snowline, which is located >200 au in most models for a 3 Msun star, but is always within 80 au for 0.5 Msun star. The chemical compositions calculated show that a well defined stellar mass and age range exists in which high C/O gas giants can form. In the case of the exoplanet HR8799b, our models show it must have formed before the star was 1 Myr old.


title: Local Coarsening Algorithms on Adaptively Refined Meshes in 2D and Their Efficient Implementation in MATLAB

abstract: Adaptive meshing includes local refinement as well as coarsening of meshes. Typically, coarsening algorithms are based on an explicit refinement history. In this work, we deal with local coarsening algorithms that build on the refinement strategies for triangular and quadrilateral meshes implemented in the ameshref package (Funken and Schmidt 2018, 2019). The ameshref package is a MATLAB-toolbox for research and teaching purposes which offers the user a certain flexibility in the REFINE step of an adaptive finite element method but can also be used in other contexts like computer graphics. This toolbox is now be extended by the coarsening option. In ameshref, no explicit information about the refinement process is stored, but is instead implicit in the data structure. In this work, we present coarsening algorithms that use easy-to-verify criteria to coarsen adaptively generated meshes by exploiting the data structure. Thereby, the desired properties are guaranteed and computational efficiency is maintained. A MATLAB implementation and some numerical examples are discussed in this work and are included in full in the toolbox ameshcoars (Funken and Schmidt 2020).


title: Direction Matters: On the Implicit Bias of Stochastic Gradient Descent with Moderate Learning Rate

abstract: Understanding the algorithmic bias of \emph{stochastic gradient descent} (SGD) is one of the key challenges in modern machine learning and deep learning theory. Most of the existing works, however, focus on \emph{very small or even infinitesimal} learning rate regime, and fail to cover practical scenarios where the learning rate is \emph{moderate and annealing}. In this paper, we make an initial attempt to characterize the particular regularization effect of SGD in the moderate learning rate regime by studying its behavior for optimizing an overparameterized linear regression problem. In this case, SGD and GD are known to converge to the unique minimum-norm solution; however, with the moderate and annealing learning rate, we show that they exhibit different \emph{directional bias}: SGD converges along the large eigenvalue directions of the data matrix, while GD goes after the small eigenvalue directions. Furthermore, we show that such directional bias does matter when early stopping is adopted, where the SGD output is nearly optimal but the GD output is suboptimal. Finally, our theory explains several folk arts in practice used for SGD hyperparameter tuning, such as (1) linearly scaling the initial learning rate with batch size; and (2) overrunning SGD with high learning rate even when the loss stops decreasing.