title: Degradation of domains with sequential field application

abstract: Recent experiments show striking unexpected features when alternating square magnetic field pulses are applied to ferromagnetic samples: domains show area reduction and domains walls change their roughness. We explain these phenomena with a simple scalar-field model, using a numerical protocol that mimics the experimental one. For a bubble and a stripe domain, we reproduce the experimental findings: The domains shrink by a combination of linear and exponential behavior. We also reproduce the roughness exponents found in the experiments. Our results suggest that the observed effects are due to a change in the disorder correlation length when the domain walls are subject to alternating fields during the first cycles, where the initial state of the interface plays a crucial role. Finally, our simulations explain the area loss by the interplay between disorder effects and effective fields induced by the local domain curvature.


title: Jet momentum broadening in the pre-equilibrium Glasma

abstract: Jets are important probes of heavy ion collisions as they can provide information on the interaction of a highly energetic parton with the medium it traverses. In the hydrodynamic stage of dense, strongly interacting matter, these interactions can be explained in terms of scattering processes, soft gluon emission and collinear parton splittings. However, jets originate even before the hydrodynamic stage. Here we report on the first numerical simulation of transverse momentum broadening of jets stemming from the interaction of partons with boost-invariantly expanding Glasma flux tubes. The Glasma stage is a pre-hydrodynamic stage based on the Color Glass Condensate framework. Our calculation shows strong time-dependence and an intrinsic anisotropy of momentum broadening in the directions transverse to the jet propagation direction.


title: Exploring the origin of low-metallicity stars in Milky Way-like galaxies with the NIHAO-UHD simulations

abstract: The kinematics of the most metal-poor stars provide a window into the early formation and accretion history of the Milky Way. Here, we use 5~high-resolution cosmological zoom-in simulations ($\sim~5\times10^6$ star particles) of Milky Way-like galaxies taken from the NIHAO-UHD project, to investigate the origin of low-metallicity stars ([Fe/H]$\leq-2.5$). The simulations show a prominent population of low-metallicity stars confined to the disk plane, as recently discovered in the Milky Way. The ubiquity of this finding suggests that the Milky Way is not unique in this respect. Independently of the accretion history, we find that $\gtrsim~90$ per cent of the retrograde stars in this population are brought in during the initial build-up of the galaxies during the first few Gyrs after the Big Bang. Our results therefore highlight the great potential of the retrograde population as a tracer of the early build-up of the Milky Way. The prograde planar population, on the other hand, is accreted during the later assembly phase and samples the full galactic accretion history. In case of a quiet accretion history, this prograde population is mainly brought in during the first half of cosmic evolution ($t\lesssim7$~Gyr), while, in the case of an on-going active accretion history, later mergers on prograde orbits are also able to contribute to this population. Finally, we note that the Milky Way shows a rather large population of eccentric, very metal-poor planar stars. This is a feature not seen in most of our simulations, with the exception of one simulation with an exceptionally active early building phase.


title: Exclusion of cosmic rays from molecular clouds by self-generated electric fields

abstract: It was recently discovered that in some regions of the Galaxy, the cosmic ray (CR) abundance is several orders of magnitude higher than previously thought. Additionally, there is evidence that in molecular cloud envelopes, the CR ionization may be dominated by electrons. We show that for regions with high, electron-dominated ionization, the penetration of CR electrons into molecular clouds is modulated by the electric field that develops as a result of the charge they deposit. We evaluate the significance of this novel mechanism of self-modulation and show that the CR penetration can be reduced by a factor of a few to a few hundred in high-ionization environments, such as those found near the Galactic center.


title: Precise predictions for $\mathrm{WH}$+jet production at the LHC

abstract: We present precise predictions for the production of a Higgs boson in association with a hadronic jet and a $\mathrm{W}$ boson at hadron colliders. The behaviour of QCD corrections are studied for fiducial cross sections and distributions of the charged gauge boson and jet-related observables. The inclusive process (at least one resolved jet) and the exclusive process (exactly one resolved jet) are contrasted and discussed. The inclusion of QCD corrections up to $\mathcal{O}(\alpha_{\text{s}}^3)$ leads to a clear stabilisation of the predictions and contributes substantially to a reduction of remaining theoretical uncertainties.


title: Scalar Fields in 3D Asymptotically Flat Higher-Spin Gravity

abstract: In this work we construct a novel associative algebra and use it to define a theory of higher-spin gravity in (2+1)-dimensional asymptotically flat spacetimes. Our construction is based on a quotient of the universal enveloping algebra (UEA) of $\mathfrak{isl}(2,\mathbb{R})$ with respect to the ideal generated by its Casimir elements, the mass squared $\mathcal{M}^2$ and the three-dimensional analogue of the square of the Pauli-Lubanski vector $\mathcal{S}$ and propose to call the resulting associative algebra $\mathfrak{ihs}(\mathcal{M}^2,\mathcal{S})$. We provide a definition of its generators and even though we are not yet able to provide the complete set of multiplication rules of this algebra our analysis allows us to study many interesting and relevant sub-structures of $\mathfrak{ihs}(\mathcal{M}^2,\mathcal{S})$. We then show how to consistently couple a scalar field to an $\mathfrak{ihs}(\mathcal{M}^2,\mathcal{S})$ higher-spin gauge theory.


title: Local and global gas metallicity versus stellar age relation in MaNGA galaxies

abstract: The search for new global scaling relations linking physical properties of galaxies has a fundamental interest. Furthermore, their recovery from spatially resolved relations has been one of the spotlights of integral field spectroscopy (IFS). In this study we investigate the existence of global and local relations between stellar age ($Age_\star$) and gas-phase metallicity ($Z_g$). To this aim, we analyze IFS data for a sample of 736 star-forming disk galaxies from the MaNGA survey. We report a positive correlation between the global $Z_g$ and $D(4000)$ (an indicator of stellar age), with a slope that decreases with increasing galaxy mass. Locally, a similar trend is found when analyzing the $Z_g$ and $D(4000)$ of the star-forming regions, as well as the residuals resulting from removing the radial gradients of both parameters. The local laws have systematically smaller slopes than the global one. We ascribe this difference to random errors, that make the true slope of the $Age_\star-Z_g$ relation to be systematically underestimated when performing a least square fitting. The explored relation is intimately linked with the already known relation between gas metallicity and star formation rate at fixed mass, both presenting a common physical origin.


title: Radiative three-body D-meson decays in and beyond the standard model

abstract: We study radiative charm decays $D \to P_1 P_2 \gamma$, $P_{1,2}=\pi,K$ in QCD factorization at leading order and within heavy hadron chiral perturbation theory. Branching ratios including resonance contributions are around $\sim 10^{-3}$ for the Cabibbo-favored modes into $K \pi \gamma$ and $\sim 10^{-5}$ for the singly Cabibbo-suppressed modes into $\pi^+ \pi^- \gamma, K^+ K^- \gamma$, and thus in reach of the flavor factories BES III and Belle II. Dalitz plots and forward-backward asymmetries reveal significant differences between the two QCD frameworks; such observables are therefore ideally suited for a data-driven identification of relevant decay mechanisms in the standard-model dominated $D \to K \pi \gamma$ decays. This increases the potential to probe new physics with the $D \to \pi^+ \pi^- \gamma$ and $D \to K^+ K^- \gamma$ decays, which are sensitive to enhanced dipole operators. CP asymmetries are useful to test the SM and look for new physics in neutral $|\Delta C|=1$ transitions. Cuts in the Dalitz plot enhance the sensitivity to new physics due to the presence of both $s$- and $t,u$-channel intermediate resonances.


title: Scattering Amplitudes for Monopoles: Pairwise Little Group and Pairwise Helicity

abstract: On-shell methods are particularly suited for exploring the scattering of electrically and magnetically charged objects, for which there is no local and Lorentz invariant Lagrangian description. In this paper we show how to construct a Lorentz-invariant S-matrix for the scattering of electrically and magnetically charged particles, without ever having to refer to a Dirac string. A key ingredient is a revision of our fundamental understanding of multi-particle representations of the Poincar\'e group. Surprisingly, the asymptotic states for electric-magnetic scattering transform with an additional little group phase, associated with pairs of electrically and magnetically charged particles. The corresponding "pairwise helicity" is identified with the quantized "cross product" of charges, $e_1 g_2 - e_2 g_1$, for every charge-monopole pair, and represents the extra angular momentum stored in the asymptotic electromagnetic field. We define a new kind of pairwise spinor-helicity variable, which serves as an additional building block for electric-magnetic scattering amplitudes. We then construct the most general 3-point S-matrix elements, as well as the full partial wave decomposition for the $2\to 2$ fermion-monopole S-matrix. In particular, we derive the famous helicity flip in the lowest partial wave as a simple consequence of a generalized spin-helicity selection rule, as well as the full angular dependence for the higher partial waves. Our construction provides a significant new achievement for the on-shell program, succeeding where the Lagrangian description has so far failed.


title: Fractionalisation and dynamics of anyons at $\nu=n+1/3$ in fractional quantum Hall effect and their experimental signatures

abstract: We show the low-lying excitations at filling factor $\nu=n+1/3$ with realistic interactions are contained completely within the well-defined Hilbert space of "Gaffnian quasiholes". Each Laughlin quasihole can thus be understood as a bound state of two Gaffnian quasiholes, which in the lowest Landau level (LLL) behaves like "partons" with "asymptotic freedom" mediated by neutral excitations acting as "gluons". Near the experimentally observed nematic FQH phase in higher LLs, quasiholes become weakly bounded and fractionalise with rich dynamical properties. By studying the effective interactions between quasiholes, we predict a finite temperature phase transition of the Laughlin quasiholes even when the Laughlin ground state remains incompressible, and derive relevant experimental conditions for its possible observations.


title: High-resolution spectroscopy of SN 2017hcc and its blueshifted line profiles from post-shock dust formation

abstract: SN2017hcc was remarkable for being a nearby and strongly polarized superluminous TypeIIn supernova (SN). We obtained high-resolution echelle spectra that we combine with other spectra to investigate its line profile evolution. All epochs reveal narrow P~Cygni components from pre-shock circumstellar material (CSM), indicating an axisymmetric outflow from the progenitor of 40-50 km/s. Intermediate-width and broad components exhibit the classic evolution seen in luminous SNe~IIn: symmetric Lorentzian profiles from pre-shock CSM lines broadened by electron scattering at early times, transitioning at late times to multi-component, irregular profiles coming from the SN ejecta and post-shock shell. As in many SNe~IIn, profiles show a progressively increasing blueshift, with a clear flux deficit in red wings of the intermediate and broad velocity components after day 200. This blueshift develops after the continuum luminosity fades, and in the intermediate-width component, persists at late times even after the SN ejecta fade. In SN2017hcc, the blueshift cannot be explained as occultation by the SN photosphere, pre-shock acceleration of CSM, or a lopsided explosion or CSM. Instead, the blueshift arises from dust formation in the post-shock shell and in the SN ejecta. The effect has a wavelength dependence characteristic of dust, exhibiting an extinction law consistent with large grains. Thus, SN2017hcc experienced post-shock dust formation and had a mildly bipolar CSM shell, similar to SN2010jl. Like other superluminous SNeIIn, the progenitor lost around 10Msun due to extreme eruptive mass loss in the decade before exploding.


title: ALMA and VLA Observations of EX Lupi in its Quiescent State

abstract: Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust.


title: On the reliability of photometric and spectroscopic tracers of halo relaxation

abstract: We characterize the relaxation state of galaxy systems by providing an assessment of the reliability of the photometric and spectroscopic probe via the semi-analytic galaxy evolution model. We quantify the correlations between the dynamical age of simulated galaxy groups and popular proxies of halo relaxation in observation, which are mainly either spectroscopic or photometric. We find the photometric indicators demonstrate a stronger correlation with the dynamical relaxation of galaxy groups compared to the spectroscopic probes. We take advantage of the Anderson Darling statistic ($A^2$) and the velocity-segregation ($\bar{\Delta V}$) as our spectroscopic indicators, and use the luminosity gap ($\Delta m_{12}$) and the luminosity de-centring ($D_{\mathrm{off-set}}$) as photometric ones. Firstly, we find that a combination of $\Delta m_{12}$ and $D_{\mathrm{off-set}}$ evaluated by a bivariant relation ($\mathrm{B} = 0.04 \times \Delta m_{12} - 0.11 \times Log(D_{off-set}) + 0.28$), shows a good correlation with the dynamical age compared to all other indicators. Secondly, by using the observational X-ray surface brightness map, we show that the bivariant relation brings about some acceptable correlations with X-ray proxies. These correlations are as well as the correlations between $A^2$ and X-ray proxies, offering reliable yet fast and economical method of quantifying the relaxation of galaxy systems. This study demonstrates that using photometric data to determine the relaxation status of a group will lead to some promising results that are comparable with the more expensive spectroscopic counterpart.


title: Large Density Perturbations from Reheating to Standard Model particles due to the Dynamics of the Higgs Boson during Inflation

abstract: Cosmic Microwave Background (CMB) observations are used to constrain reheating to Standard Model (SM) particles after a period of inflation. As a light spectator field, the SM Higgs boson acquires large field values from its quantum fluctuations during inflation, gives masses to SM particles that vary from one Hubble patch to another, and thereby produces large density fluctuations. We consider both perturbative and resonant decay of the inflaton to SM particles. For the case of perturbative decay from coherent oscillations of the inflaton after high scale inflation, we find strong constraints on the reheat temperature for the inflaton decay into heavy SM particles. For the case of resonant particle production (preheating) to (Higgsed) SM gauge bosons, we find temperature fluctuations larger than observed in the CMB for a range of gauge coupling that includes those found in the SM and conclude that such preheating cannot be the main source of reheating the Universe after inflation.


title: MUSE Analysis of Gas around Galaxies (MAGG) -- II: Metal-enriched halo gas around z~1 galaxies

abstract: We present a study of the metal-enriched cool halo gas traced by MgII absorption around 228 galaxies at z~0.8-1.5 within 28 quasar fields from the MUSE Analysis of Gas around Galaxies (MAGG) survey. We observe no significant evolution in the MgII equivalent width versus impact parameter relation and in the MgII covering fraction compared to surveys at z<~0.5. The stellar mass, along with distance from galaxy centre, appears to be the dominant factor influencing the MgII absorption around galaxies. With a sample that is 90% complete down to a star formation rate of ~0.1 Msun/yr and up to impact parameters ~250-350 kpc from quasars, we find that the majority (67^{+12}_{-15}% or 14/21) of the MgII absorption systems are associated with more than one galaxy. The complex distribution of metals in these richer environments adds substantial scatter to previously-reported correlations. Multiple galaxy associations show on average five times stronger absorption and three times higher covering fraction within twice the virial radius than isolated galaxies. The dependence of MgII absorption on galaxy properties disfavours the scenario in which a widespread intra-group medium dominates the observed absorption. This leaves instead gravitational interactions among group members or hydrodynamic interactions of the galaxy haloes with the intra-group medium as favoured mechanisms to explain the observed enhancement in the MgII absorption strength and cross section in rich environments.


title: Neutrino Oscillation Constraints on U(1)' Models: from Non-Standard Interactions to Long-Range Forces

abstract: We quantify the effect of gauge bosons from a weakly coupled lepton flavor dependent $U(1)'$ interaction on the matter background in the evolution of solar, atmospheric, reactor and long-baseline accelerator neutrinos in the global analysis of oscillation data. The analysis is performed for interaction lengths ranging from the Sun-Earth distance to effective contact neutrino interactions. We survey $\sim 10000$ set of models characterized by the six relevant fermion $U(1)'$ charges and find that in all cases, constraints on the coupling and mass of the $Z'$ can be derived. We also find that about 5% of the $U(1)'$ model charges lead to a viable LMA-D solution but this is only possible in the contact interaction limit. We explicitly quantify the constraints for a variety of models including $U(1)_{B-3L_e}$, $U(1)_{B-3L_\mu}$, $U(1)_{B-3L_\tau}$, $U(1)_{B-\frac{3}{2}(L_\mu+L_\tau)}$, $U(1)_{L_e-L_\mu}$, $U(1)_{L_e-L_\tau}$, $U(1)_{L_e-\frac{1}{2}(L_\mu+L_\tau)}$. We compare the constraints imposed by our oscillation analysis with the strongest bounds from fifth force searches, violation of equivalence principle as well as bounds from scattering experiments and white dwarf cooling. Our results show that generically, the oscillation analysis improves over the existing bounds from gravity tests for $Z'$ lighter than $\sim 10^{-11}$ eV. In the contact interaction limit, we find that for most models listed above there are values of $g'$ and $M_{Z'}$ for which the oscillation analysis provides constraints beyond those imposed by laboratory experiments. Finally we illustrate the range of $Z'$ and couplings leading to a viable LMA-D solution for two sets of models.


title: Topological superconductivity in proximity to type-II superconductors

abstract: One-dimensional systems proximity-coupled to a superconductor can be driven into a topological superconducting phase by an external magnetic field. Here, we investigate the effect of vortices created by the magnetic field in a type-II superconductor providing the proximity effect. We identify different ways in which the topological protection of Majorana modes can be compromised and discuss strategies to circumvent these detrimental effects. Our findings are also relevant to topological phases of proximitized quantum Hall edge states.


title: The breakBRD Breakdown: Using IllustrisTNG to Track the Quenching of an Observationally-Motivated Sample of Centrally Star-Forming Galaxies

abstract: The observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D$_n$4000 spectral index (Tuttle & Tonnesen 2020). In this paper, we search for breakBRD analogues in the IllustrisTNG simulation and describe their history and future. We find that a small fraction (${\sim}4\%$ at $z=0$; ${\sim}1\%$ at $z=0.5$) of galaxies fulfill the breakBRD criteria, in agreement with observations. In comparison with the mass-weighted parent IllustrisTNG sample, these galaxies tend to consist of a higher fraction of satellite and splashback galaxies. However, the central, non-splashback breakBRD galaxies show similar environments, black hole masses, and merger rates, indicating that there is not a single formation trigger for inner star formation and outer quenching. We determine that breakBRD analogue galaxies as a whole are in the process of quenching. The breakBRD state - with its highly centrally concentrated star formation - is uncommon in the history of either currently quiescent or star-forming galaxies; however, approximately 10% of $10^{10} < M_\ast/M_{\odot} < 10^{11}$ quiescent galaxies at $z=0$ have experienced SFR concentrations comparable to those of the breakBRDs in their past. Additionally, the breakBRD state is short-lived, lasting a few hundred Myr up to ${\sim}2$ Gyr. The observed breakBRD galaxies may therefore be a unique sample of outside-in quenching galaxies.


title: Does locality plus perfect correlation imply determinism?

abstract: A 1964 paper by John Bell gave the first demonstration that quantum mechanics is incompatible with local hidden variables. There is an ongoing and vigorous debate on whether he relied on an assumption of determinism, or instead, as he later claimed, derived determinism from assumptions of locality and perfect correlation. This paper aims to bring clarity to the debate via simple examples and rigorous results. It is shown that the weak form of locality used in Bell's 1964 paper (parameter independence) is insufficient for such a derivation, whereas an independent form called outcome independence is sufficient even when weak locality does not hold. It further follows that outcome independence, by itself, implies that standard quantum mechanics is incomplete. It is also shown that an appeal by Bell to the Einstein-Rosen-Podolsky argument to support his claim fails, via examples that expose logical gaps in this argument. However, replacing the reality criterion underpinning the argument by a stronger criterion enables a rigorous derivation of both weak locality and determinism, as required for Bell's 1964 paper. Consequences for quantum interpretations, locality, and classical common causes are briefly discussed, with reference to an example of local classical indeterminism.


title: Electronic structures, charge transfer and charge orders in twisted transition metal dichalcogenide bilayers

abstract: Moire superlattices of transition metal dichalcogenide (TMD) bilayers have been shown to host correlated electronic states, which arises from the interplay of emergent moire potential and long-range Coulomb interactions. Here we theoretically investigate structural relaxation and single-particle electronic properties in moire superlattices of transition metal dichalcogenide homobilayer and study the ground state charge orders in the effective honeycomb lattice of MX and XM region. From the large-scale density functional theory calculation and continuum model with layer degrees of freedom, we find that the out of plane gating field creates a tunable charge transfer gap and introduces a mass term in the Dirac spectrum. At the flat band limit, we observe a series of charge-ordered insulating states at various filling n = 1/4, 1/3, 1/2, 2/3, 1 with Monte Carlo simulations, and predict that gating field induces a phase transition between different electron crystals at fixed filling n = 1/2, 2/3. Our work demonstrates that transition metal dichalcogenide homobilayer provides a powerful platform for the investigation of tunable charge transfer insulator and charge orders.


title: Simultaneous modelling of matter power spectrum and bispectrum in the presence of baryons

abstract: We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and equilateral bispectrum on scales up to k < 5 h/Mpc and redshifts z<2, as measured in six different cosmological hydrodynamical simulations. The accuracy of our fits are typically 1% for the power spectrum, and for the equilateral and squeezed bispectra, which somewhat degrades to 3% for simulations with extreme feedback prescriptions. Our results support the physical assumptions underlying baryonification approaches, and encourage their use in interpreting weak gravitational lensing and other cosmological observables.


title: Union-Find Decoders For Homological Product Codes

abstract: Homological product codes are a class of codes that can have improved distance while retaining relatively low stabilizer weight. We show how to build union-find decoders for these codes, using a union-find decoder for one of the codes in the product and a brute force decoder for the other code. We apply this construction to the specific case of the product of a surface code with a small code such as a $[[4,2,2]]$ code, which we call an augmented surface code. The distance of the augmented surface code is the product of the distance of the surface code with that of the small code, and the union-find decoder, with slight modifications, can decode errors up to half the distance. We present numerical simulations, showing that while the threshold of these augmented codes is lower than that of the surface code, the low noise performance is improved.


title: Spectral and Imaging properties of Sgr A* from High-Resolution 3D GRMHD Simulations with Radiative Cooling

abstract: The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-relativistic magneto-hydrodynamical simulations, we study the impact of radiative cooling on the dynamical evolution of the accreting plasma, presenting spectral energy distributions and synthetic sub-millimeter images generated from the accretion flow around Sgr A*. These simulations solve the approximated equations for radiative cooling processes self-consistently, including synchrotron, bremsstrahlung, and inverse Compton processes. We find that radiative cooling plays an increasingly important role in the dynamics of the accretion flow as the accretion rate increases: the mid-plane density grows and the infalling gas is less turbulent as cooling becomes stronger. The changes in the dynamical evolution become important when the accretion rate is larger than $10^{-8}\,M_{\odot}~{\rm yr}^{-1}$ ($\gtrsim 10^{-7} \dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}$ is the Eddington accretion rate). The resulting spectra in the cooled models also differ from those in the non-cooled models: the overall flux, including the peak values at the sub-mm and the far-UV, is slightly lower as a consequence of a decrease in the electron temperature. Our results suggest that radiative cooling should be carefully taken into account in modelling Sgr A* and other low-luminosity active galactic nuclei that have a mass accretion rate of $\dot{M} > 10^{-7}\,\dot{M}_{\rm Edd}$.


title: Planet formation and disk mass dependence in a pebble-driven scenario for low mass stars

abstract: Measured disk masses seem to be too low to form the observed population of planetary systems. In this context, we develop a population synthesis code in the pebble accretion scenario, to analyse the disk mass dependence on planet formation around low mass stars. We base our model on the analytical sequential model presented in Ormel et al. 2017 and analyse the populations resulting from varying initial disk mass distributions. Starting out with seeds the mass of Ceres near the ice-line formed by streaming instability, we grow the planets using the Pebble Accretion process and migrate them inwards using Type-I migration. The next planets are formed sequentially after the previous planet crosses the ice-line. We explore different initial distributions of disk masses to show the dependence of this parameter with the final planetary population. Our results show that compact close-in resonant systems can be pretty common around M-dwarfs between $0.09-0.2$ $M_{\odot}$ only when the disks considered are more massive than what is being observed by sub-mm disk surveys. The minimum disk mass to form a Mars-like planet is found to be about $2 \times 10^{-3}$ $M_{\odot}$. Small variation in the disk mass distribution also manifest in the simulated planet distribution. The paradox of disk masses might be caused by an underestimation of the disk masses in observations, by a rapid depletion of mass in disks by planets growing within a million years or by deficiencies in our current planet formation picture.


title: ParaMonte: A high-performance serial/parallel Monte Carlo simulation library for C, C++, Fortran

abstract: ParaMonte (standing for Parallel Monte Carlo) is a serial and MPI/Coarray-parallelized library of Monte Carlo routines for sampling mathematical objective functions of arbitrary-dimensions, in particular, the posterior distributions of Bayesian models in data science, Machine Learning, and scientific inference. The ParaMonte library has been developed with the design goal of unifying the **automation**, **accessibility**, **high-performance**, **scalability**, and **reproducibility** of Monte Carlo simulations. The current implementation of the library includes **ParaDRAM**, a **Para**llel **D**elyaed-**R**ejection **A**daptive **M**etropolis Markov Chain Monte Carlo sampler, accessible from a wide range of programming languages including C, C++, Fortran, with a unified Application Programming Interface and simulation environment across all supported programming languages. The ParaMonte library is MIT-licensed and is permanently located and maintained at [https://github.com/cdslaborg/paramonte](https://github.com/cdslaborg/paramonte).

categories: cs.MS astro-ph.IM physics.data-an q-bio.QM stat.ML


title: On theorems of Chernoff and Ingham on the Heisenberg group

abstract: We prove an analogue of Chernoff's theorem for the sublaplacian on the Heisenberg group and use it prove a version of Ingham's theorem for the Fourier transform on the same group.


title: Quantitative and Interpretable Order Parameters for Phase Transitions from Persistent Homology

abstract: We apply modern methods in computational topology to the task of discovering and characterizing phase transitions. As illustrations, we apply our method to four two-dimensional lattice spin models: the Ising, square ice, XY, and fully-frustrated XY models. In particular, we use persistent homology, which computes the births and deaths of individual topological features as a coarse-graining scale or sublevel threshold is increased, to summarize multiscale and high-point correlations in a spin configuration. We employ vector representations of this information called persistence images to formulate and perform the statistical task of distinguishing phases. For the models we consider, a simple logistic regression on these images is sufficient to identify the phase transition. Interpretable order parameters are then read from the weights of the regression. This method suffices to identify magnetization, frustration, and vortex-antivortex structure as relevant features for phase transitions in our models. We also define "persistence" critical exponents and study how they are related to those critical exponents usually considered.


title: Evidence for Late-Time Feedback from the Discovery of Multiphase Gas in a Massive Elliptical at $z=0.4$

abstract: We report the first detection of multiphase gas within a quiescent galaxy beyond $z\approx0$. The observations use the brighter image of doubly lensed QSO HE 0047$-$1756 to probe the ISM of the massive ($M_{\rm star}\approx 10^{11} \mathrm{M_\odot}$) elliptical lens galaxy at $z_\mathrm{gal}=0.408$. Using Hubble Space Telescope's Cosmic Origins Spectrograph (COS), we obtain a medium-resolution FUV spectrum of the lensed QSO and identify numerous absorption features from $\mathrm{H_2}$ in the lens ISM at projected distance $d=4.6$ kpc. The $\mathrm{H_2}$ column density is $\log N(\mathrm{H_2})/\mathrm{cm^{-2}}=17.8^{+0.1}_{-0.3}$ with a molecular gas fraction of $f_\mathrm{H_2}=2-5\%$, roughly consistent with some local quiescent galaxies. The new COS spectrum also reveals kinematically complex absorption features from highly ionized species O VI and N V with column densities log $N(\mathrm{O VI})/\mathrm{cm^{-2}} =15.2\pm0.1$ and log $N(\mathrm{N V})/\mathrm{cm^{-2}} =14.6\pm0.1$, among the highest known in external galaxies. Assuming the high-ionization absorption features originate in a transient warm ($T\sim10^5\,$K) phase undergoing radiative cooling from a hot halo surrounding the galaxy, we infer a mass accretion rate of $\sim 0.5-1.5\,\mathrm{M_\odot\,yr^{-1}}$. The lack of star formation in the lens suggests the bulk of this flow is returned to the hot halo, implying a heating rate of $\sim10^{48}\,\mathrm{erg\,yr^{-1}}$. Continuous heating from evolved stellar populations (primarily SNe Ia but also winds from AGB stars) may suffice to prevent a large accumulation of cold gas in the ISM, even in the absence of strong feedback from an active nucleus.


title: Lip-reading with Densely Connected Temporal Convolutional Networks

abstract: In this work, we present the Densely Connected Temporal Convolutional Network (DC-TCN) for lip-reading of isolated words. Although Temporal Convolutional Networks (TCN) have recently demonstrated great potential in many vision tasks, its receptive fields are not dense enough to model the complex temporal dynamics in lip-reading scenarios. To address this problem, we introduce dense connections into the network to capture more robust temporal features. Moreover, our approach utilises the Squeeze-and-Excitation block, a light-weight attention mechanism, to further enhance the model's classification power. Without bells and whistles, our DC-TCN method has achieved 88.36% accuracy on the Lip Reading in the Wild (LRW) dataset and 43.65% on the LRW-1000 dataset, which has surpassed all the baseline methods and is the new state-of-the-art on both datasets.


title: On the possibility of through passage of asteroid bodies across the Earth's atmosphere

abstract: We have studied the conditions of through passage of asteroids with diameters 200, 100 and 50 m, consisting of three types of materials -- iron, stone and water ice across the Earth's atmosphere with the minimum trajectory altitude 10--15 km. The conditions of this passage with subsequent exit into outer space with the preservation of a substantial fraction of the initial mass have been found. The results obtained support our idea explaining one of the long-standing problems of astronomy -- the Tunguska phenomenon which has not received reasonable and comprehensive interpretations to date. We argue that Tunguska event was caused by an iron asteroid body, which passed through the Earth's atmosphere and continued to the near-solar orbit.


title: Orbital structure of the effective pairing interaction in the high-temperature superconducting cuprates

abstract: The nature of the effective interaction responsible for pairing in the high-temperature superconducting cuprates remains unsettled. This question has been studied extensively using the simplified single-band Hubbard model, which does not explicitly consider the orbital degrees of freedom of the relevant CuO$_2$ planes. Here, we use a dynamic cluster quantum Monte Carlo approximation to study the orbital structure of the pairing interaction in the three-band Hubbard model, which treats the orbital degrees of freedom explicitly. We find that the interaction predominately acts between neighboring copper orbitals, but with significant additional weight appearing on the surrounding bonding molecular oxygen orbitals. By explicitly comparing these results to those from the simpler single-band Hubbard model, our study provides strong support for the single-band framework for describing superconductivity in the cuprates.


title: Smith theory and cyclic base change functoriality

abstract: Lafforgue and Genestier-Lafforgue have constructed the global and (semisimplified) local Langlands correspondences for arbitrary reductive groups over function fields. We establish various properties of these correspondences regarding functoriality for cyclic base change: For $\mathbb{Z}/p\mathbb{Z}$-extensions of global function fields, we prove the existence of base change for mod $p$ automorphic forms on arbitrary reductive groups. For $\mathbb{Z}/p\mathbb{Z}$-extensions of local function fields, we construct a base change homomorphism for the mod $p$ Bernstein center of any reductive group. We then use this to prove existence of local base change for mod $p$ irreducible representation along $\mathbb{Z}/p\mathbb{Z}$-extensions for all large enough $p$, and that Tate cohomology realizes descent along base change, verifying a function field version of a conjecture of Treumann-Venkatesh. The proofs are based on equivariant localization arguments for the moduli spaces of shtukas. They also draw upon new tools from representation theory, including parity sheaves and Smith-Treumann theory. In particular, we use these to establish a categorification of the base change homomorphism for mod $p$ spherical Hecke algebras, in a joint appendix with Gus Lonergan.


title: Augmenting Scientific Papers with Just-in-Time, Position-Sensitive Definitions of Terms and Symbols

abstract: Despite the central importance of research papers to scientific progress, they can be difficult to read. Comprehension is often stymied when the information needed to understand a passage resides somewhere else: in another section, or in another paper. In this work, we envision how interfaces can bring definitions of technical terms and symbols to readers when and where they need them most. We introduce ScholarPhi, an augmented reading interface with four novel features: (1) tooltips that surface position-sensitive definitions from elsewhere in a paper, (2) a filter over the paper that "declutters" it to reveal how the term or symbol is used across the paper, (3) automatic equation diagrams that expose multiple definitions in parallel, and (4) an automatically generated glossary of important terms and symbols. A usability study showed that the tool helps researchers of all experience levels read papers. Furthermore, researchers were eager to have ScholarPhi's definitions available to support their everyday reading.


title: ATLASGAL-selected massive clumps in the inner Galaxy: VIII. Chemistry of photodissociation regions

abstract: We study ten molecular transitions obtained from an unbiased 3 mm molecular line survey using the IRAM 30 m telescope toward 409 compact dust clumps identified by the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) to understand photodissociation regions (PDRs) associated with the clumps. The main goal of this study is to investigate whether the abundances of the selected molecules show any variations resulting from the PDR chemistry in different clump environments. We selected HCO, HOC$^+$, C$_{2}$H, c-C$_{3}$H$_{2}$, CN, H$^{13}$CN, HC$^{15}$N, and HN$^{13}$C as PDR tracers, and H$^{13}$CO$^{+}$ and C$^{18}$O as dense gas tracers. We find that the abundances of HCO, CN, C$_{2}$H and c-C$_{3}$H$_{2}$ decrease as the H$_{2}$ column density increase, indicating high visual extinction, while those of high density tracers (i.e., H$^{13}$CO$^{+}$ and HC$^{15}$N) are constant. In addition, $N$(HCO)/$N$(H$^{13}$CO$^{+}$) ratios significantly decrease as H$_{2}$ column density increase, and in particular, 82 clumps have $X$(HCO) $\gtrsim 10^{-10}$ and $N$(HCO)/$N$(H$^{13}$CO$^{+}$) $\gtrsim 1$, which are the indication of far-ultraviolet (FUV) chemistry. This suggests the observed HCO abundances are likely associated with FUV radiation illuminating the PDRs. We also find that high $N$(c-C$_{3}$H$_{2}$)/$N$(C$_{2}$H) ratios found for HII regions having high HCO abundances ($\gtrsim\,10^{-10}$) are associated with more evolved clumps with high $L_{\rm bol}$/$M_{\rm clump}$. This trend might be associated with gain-surface processes, which determine initial abundances of these molecules, and time-dependent effects in the clumps corresponding to the envelopes around dense PDRs and HII regions. In addition, some fraction of the measured abundances of the small hydrocarbons of the HII sources can be the result of the photodissociation of PAH molecules.


title: Couplings for Andersen Dynamics

abstract: Andersen dynamics is a standard method for molecular simulations, and a precursor of the Hamiltonian Monte Carlo algorithm used in MCMC inference. The stochastic process corresponding to Andersen dynamics is a PDMP (piecewise deterministic Markov process) that iterates between Hamiltonian flows and velocity randomizations of randomly selected particles. Both from the viewpoint of molecular dynamics and MCMC inference, a basic question is to understand the convergence to equilibrium of this PDMP particularly in high dimension. Here we present couplings to obtain sharp convergence bounds in the Wasserstein sense that do not require global convexity of the underlying potential energy.

categories: math.PR physics.chem-ph physics.comp-ph stat.ML


title: Hybrid European MV-LV Network Models for Smart Distribution Network Modelling

abstract: A pair of European-style, integrated MV-LV circuits are presented, created by combining generic MV and real LV networks. The two models have 86,000 and 113,000 nodes, and are made readily available for download in the OpenDSS file format. Primary substation tap change controls and MV-LV feeders are represented as three-phase unbalanced distribution network models, capturing the coupling of voltages at the MV level. The assumptions made in constructing the models are outlined, including a preconditioning step that reduces the number of nodes by more than five times without affecting the solution. Two flexibility-based case studies are presented, with TSO-DSO and peer-peer-based smart controls considered. The demonstration of the heterogeneous nature of these systems is corroborated by the analysis of measured LV voltage data. The models are intended to aid the development of algorithms for maximising the benefits of smart devices within the context of whole energy systems.


title: Interpreting Internal Consistency of DES Measurements

abstract: Bayesian evidence ratios are widely used to quantify the statistical consistency between different experiments. However, since the evidence ratio is prior dependent, the precise translation between its value and the degree of concordance/discordance requires additional information. The most commonly adopted metric, the Jeffreys scale, can falsely suggest agreement between datasets when priors are chosen to be sufficiently wide. In this work, we examine evidence ratios in a DES-Y1 simulated analysis, focusing on the internal consistency between weak lensing and galaxy clustering. We study two scenarios using simulated data in controlled experiments. First, we calibrate the expected evidence ratio distribution given noise realizations around the best fit DES-Y1 $\Lambda$CDM cosmology. Second, we show the behavior of evidence ratios for noiseless fiducial data vectors simulated using a modified gravity model, which generates internal tension in the $\Lambda$CDM analysis. We show that the choice of prior could conceal the discrepancies between weak lensing and galaxy clustering induced by such models and that the evidence ratio in a DES-Y1 study is, indeed, biased towards agreement.


title: Modular invariant flavor model of $\rm A_4$ and hierarchical structures at nearby fixed points

abstract: In the modular invariant flavor model of $\rm A_4$, we study the hierarchical structure of lepton/quark flavors at nearby fixed points of $\tau=i$ and $\tau=\omega$ of the modulus, which are in the fundamental domain of $\rm PSL(2,\mathbb{Z})$. These fixed points correspond to the residual symmetries $\mathbb{Z}_2^{S}=\{I, S \}$ and $\mathbb{Z}_3^{S}=\{ I, ST, (ST)^2 \}$ of $\rm A_4$, where $S$ and $T$ are generators of the $A_4$ group. The infinite $\tau= i \infty$ also preserves the residual symmetry of the subgroup $\mathbb{Z}^T_3=\{ I,T,T^2 \}$ of $\rm A_4$. We study typical two-type mass matrices for charged leptons and quarks in terms of modular forms of weights $2$, $4$ and $6$ while the neutrino mass matrix with the modular forms of weight $4$ through the Weinberg operator. Linear modular forms are obtained approximately by performing Taylor expansion of modular forms around fixed points. By using them, the flavor structure of the lepton and quark mass matrices are examined at nearby fixed points. The hierarchical structure of these mass matrices is clearly shown in the diagonal base of $S$, $T$ and $ST$. The observed PMNS and CKM mixing matrices can be reproduced at nearby fixed points in some cases of mass matrices. By scanning model parameters numerically at nearby fixed points, our discussion are confirmed for both the normal hierarchy and inverted one of neutrino masses. Predictions are given for the sum of neutrino masses and the CP violating Dirac phase of leptons at each nearby fixed point.


title: Temporal State Machines: Using temporal memory to stitch time-based graph computations

abstract: Race logic, an arrival-time-coded logic family, has demonstrated energy and performance improvements for applications ranging from dynamic programming to machine learning. However, the ad hoc mappings of algorithms into hardware result in custom architectures making them difficult to generalize. We systematize the development of race logic by associating it with the mathematical field called tropical algebra. This association between the mathematical primitives of tropical algebra and generalized race logic computations guides the design of temporally coded tropical circuits. It also serves as a framework for expressing high level timing-based algorithms. This abstraction, when combined with temporal memory, allows for the systematic generalization of race logic by making it possible to partition feed-forward computations into stages and organizing them into a state machine. We leverage analog memristor-based temporal memories to design a such a state machine that operates purely on time-coded wavefronts. We implement a version of Dijkstra's algorithm to evaluate this temporal state machine. This demonstration shows the promise of expanding the expressibility of temporal computing to enable it to deliver significant energy and throughput advantages.


title: Acceleration of Large Margin Metric Learning for Nearest Neighbor Classification Using Triplet Mining and Stratified Sampling

abstract: Metric learning is one of the techniques in manifold learning with the goal of finding a projection subspace for increasing and decreasing the inter- and intra-class variances, respectively. Some of the metric learning methods are based on triplet learning with anchor-positive-negative triplets. Large margin metric learning for nearest neighbor classification is one of the fundamental methods to do this. Recently, Siamese networks have been introduced with the triplet loss. Many triplet mining methods have been developed for Siamese networks; however, these techniques have not been applied on the triplets of large margin metric learning for nearest neighbor classification. In this work, inspired by the mining methods for Siamese networks, we propose several triplet mining techniques for large margin metric learning. Moreover, a hierarchical approach is proposed, for acceleration and scalability of optimization, where triplets are selected by stratified sampling in hierarchical hyper-spheres. We analyze the proposed methods on three publicly available datasets, i.e., Fisher Iris, ORL faces, and MNIST datasets.


title: Compactness versus hugeness at successor cardinals

abstract: If $\kappa$ is regular and $2^{<\kappa}\leq\kappa^+$, then the existence of a weakly presaturated ideal on $\kappa^+$ implies $\square^*_\kappa$. This partially answers a question of Foreman and Magidor about the approachability ideal on $\omega_2$. As a corollary, we show that if there is a presaturated ideal $I$ on $\omega_2$ such that $\mathcal{P}(\omega_2)/I$ is semiproper, then CH holds. We also show some barriers to getting the tree property and a saturated ideal simultaneously on a successor cardinal from conventional forcing methods.


title: Influence of distinct kinds of temporal disorder in discontinuous phase transitions

abstract: Based on the MFT arguments, a general description for discontinuous phase transitions in the presence temporal disorder is considered. Our analysis extends the recent findings [Phys. Rev. E {\bf 98}, 032129 (2018)] by considering other kinds of phase transitions beyond the absorbing ones. The theory is exemplified in one of the simplest (nonequilibrium) order disorder (discontinuous) phase transition with "up-down" $Z_2$ symmetry: the inertial majority vote (IMV) model for two kinds of temporal disorder. As for the APT case, the temporal disorder does not suppress the occurrence of discontinuous phase transitions, but remarkable differences emerge when compared with the pure case. A comparison between the distinct kinds of temporal disorder is also performed beyond the MFT for random-regular (RR) complex topologies.


title: Hubble Constant Measurement with GW190521 as an Eccentric Black Hole Merger

abstract: Gravitational wave observations can be used to accurately measure the Hubble constant $H_0$ and could help understand the present discrepancy between constraints from Type Ia supernovae and the cosmic microwave background. Neutron star mergers are primarily used for this purpose as their electromagnetic emission can be used to greatly reduce measurement uncertainties. Here we estimate $H_0$ using the recently observed black hole merger GW190521 and its candidate electromagnetic counterpart found by ZTF using a highly eccentric explanation of the properties of GW190521. We find that the reconstructed distance of GW190521 and the redshift of the candidate host galaxy are more consistent with standard cosmology for our eccentric model than if we reconstruct the source parameters assuming no eccentricity. We obtain $H_0=88.6^{+17.1}_{-34.3}$\,km\,s$^{-1}$Mpc$^{-1}$ for GW190521, and $H_0=73.4^{+6.9}_{-10.7}$\,km\,s$^{-1}$Mpc$^{-1}$ in combination with the results of the neutron star merger GW170817. Our results indicate that future $H_0$ computations using black hole mergers will need to account for possible eccentricity. For extreme cases, the orbital velocity of binaries in AGN disks can represent a significant systematic uncertainty.


title: Ensemble Multi-Source Domain Adaptation with Pseudolabels

abstract: Given multiple source datasets with labels, how can we train a target model with no labeled data? Multi-source domain adaptation (MSDA) aims to train a model using multiple source datasets different from a target dataset in the absence of target data labels. MSDA is a crucial problem applicable to many practical cases where labels for the target data are unavailable due to privacy issues. Existing MSDA frameworks are limited since they align data without considering conditional distributions p(x|y) of each domain. They also miss a lot of target label information by not considering the target label at all and relying on only one feature extractor. In this paper, we propose Ensemble Multi-source Domain Adaptation with Pseudolabels (EnMDAP), a novel method for multi-source domain adaptation. EnMDAP exploits label-wise moment matching to align conditional distributions p(x|y), using pseudolabels for the unavailable target labels, and introduces ensemble learning theme by using multiple feature extractors for accurate domain adaptation. Extensive experiments show that EnMDAP provides the state-of-the-art performance for multi-source domain adaptation tasks in both of image domains and text domains.


title: Enhanced Bayesian Model Updating with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

abstract: Finite element model updating is challenging because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield responses that are similar with respect to actual measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic finite element model updating using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element model updating. The validity of this framework is systematically demonstrated through case studies.


title: A Framework of Learning Through Empirical Gain Maximization

abstract: We develop in this paper a framework of empirical gain maximization (EGM) to address the robust regression problem where heavy-tailed noise or outliers may present in the response variable. The idea of EGM is to approximate the density function of the noise distribution instead of approximating the truth function directly as usual. Unlike the classical maximum likelihood estimation that encourages equal importance of all observations and could be problematic in the presence of abnormal observations, EGM schemes can be interpreted from a minimum distance estimation viewpoint and allow the ignorance of those observations. Furthermore, it is shown that several well-known robust nonconvex regression paradigms, such as Tukey regression and truncated least square regression, can be reformulated into this new framework. We then develop a learning theory for EGM, by means of which a unified analysis can be conducted for these well-established but not fully-understood regression approaches. Resulting from the new framework, a novel interpretation of existing bounded nonconvex loss functions can be concluded. Within this new framework, the two seemingly irrelevant terminologies, the well-known Tukey's biweight loss for robust regression and the triweight kernel for nonparametric smoothing, are closely related. More precisely, it is shown that the Tukey's biweight loss can be derived from the triweight kernel. Similarly, other frequently employed bounded nonconvex loss functions in machine learning such as the truncated square loss, the Geman-McClure loss, and the exponential squared loss can also be reformulated from certain smoothing kernels in statistics. In addition, the new framework enables us to devise new bounded nonconvex loss functions for robust learning.


title: Detection of Pb II in the Ultraviolet Spectra of Three Metal-Poor Stars

abstract: We report the first detection of the Pb II line at 2203.534 Angstroms in three metal-poor stars, using ultraviolet spectra obtained with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. We perform a standard abundance analysis assuming local thermodynamic equilibrium (LTE) to derive lead (Pb, Z=82) abundances. The Pb II line yields a higher abundance than Pb I lines by +0.36 +/- 0.34 dex and +0.49 +/- 0.28 dex in the stars HD 94028 and HD 196944, where Pb I lines had been detected previously. The Pb II line is likely formed in LTE, and these offsets affirm previous calculations showing that Pb I lines commonly used as abundance indicators underestimate the Pb abundance in LTE. Pb is enhanced in the s-process-enriched stars HD 94028 ([Pb/Fe] = +0.95 +/- 0.14) and HD 196944 ([Pb/Fe] = +2.28 +/- 0.23), and we show that Pb-208 is the dominant Pb isotope in these two stars. The log epsilon(Pb/Eu) ratio in the r-process-enhanced star HD 222925 is 0.76 +/- 0.14, which matches the Solar System r-process ratio and indicates that the Solar System r-process residuals for Pb are, in aggregate, correct. The Th/Pb chronometer in HD 222925 yields an age of 8.2 +/- 5.8 Gyr, and we highlight the potential of the Th/Pb chronometer as a relatively model-insensitive age indicator in r-process-enhanced stars.


title: Covariance Steering of Discrete-Time Stochastic Linear Systems Based on Distribution Distance Terminal Costs

abstract: We consider a class of stochastic optimal control problems for discrete-time stochastic linear systems which seek for control policies that will steer the probability distribution of the terminal state of the system close to a desired Gaussian distribution. In our problem formulation, the closeness between the terminal state distribution and the desired (goal) distribution is measured in terms of the squared Wasserstein distance which is associated with a corresponding terminal cost term. We recast the stochastic optimal control problem as a finite-dimensional nonlinear program and we show that its performance index can be expressed as the difference of two convex functions. This representation of the performance index allows us to find local minimizers of the original nonlinear program via the so-called convex-concave procedure [1]. Subsequently, we consider a similar problem but this time we use a terminal cost that corresponds to the KL divergence. Finally, we present non-trivial numerical simulations to demonstrate the proposed techniques and compare them in terms of computation time.


title: Integrability of local and nonlocal non-commutative fourth order quintic nonlinear Schrodinger equations

abstract: We prove integrability of a generalised non-commutative fourth order quintic nonlinear Schrodinger equation. The proof is relatively succinct and rooted in the linearisation method pioneered by Ch. Poppe. It is based on solving the corresponding linearised partial differential system to generate an evolutionary Hankel operator for the `scattering data'. The time-evolutionary solution to the non-commutative nonlinear partial differential system is then generated by solving a linear Fredholm equation which corresponds to the Marchenko equation. The integrability of reverse space-time and reverse time nonlocal versions, in the sense of Ablowitz and Musslimani, of the fourth order quintic nonlinear Schrodinger equation are proved contiguously by the approach adopted. Further, we implement a numerical integration scheme based on the analytical approach above which involves solving the linearised partial differential system followed by numerically solving the linear Fredholm equation to generate the solution at any given time.


title: In Situ characterization of the proton enhanced conductivity of 20nm TiO2 thin-films obtained on the surface of optical fiber

abstract: We report in situ characterized TiO2 thin-films deposited on optical fiber, having thicknesses in the 20-100nm range, and having enhanced conductivity values of 700S/cm upon interacting with hydrogen. This conductivity was achieved in pure hydrogen at 800-900C, having a measured activation energy of 0.26eV of the hopping type. Given the variability in the observed results, it is postulated that the highest conductivity achievable may be much greater than what is currently demonstrated. The conductivity is retained after cooling to ambient temperatures as confirmed by Hall measurements, and subsequent grazing-incidence x-ray diffraction and TEM measurements show the films to be in the rutile phase. The exceptional conductivity in these films is hypothesized to result from direct proton incorporation into the lattice populating the conduction band with excess electrons, or from altering the Titania lattice to form conductive Magneli phases. The films did not display any evidence of transformations, however formation of Magneli phases was confirmed for powders. These interesting results, observed by examining 20nm films on the surface of optical fiber in combination with the first impedance spectroscopy performed on films on optical fiber in high temperature Fuel Cell type environments, confirm hypotheses arrived at in prior publications where thin-films of Titania had optical properties which could only be explained by the current claim. Titania thin-films on optical-fiber are being explored for high temperature hydrogen derived energy generation, thermo-photonic energy conversion, and associated sensors due to their unique interactions with hydrogen.

categories: cond-mat.mtrl-sci cond-mat.mes-hall physics.app-ph physics.ins-det physics.optics


title: A Note on Measure-Valued Solutions to the Full Euler System

abstract: We construct two particular solutions of the full Euler system which emanate from the same initial data. Our aim is to show that the convex combination of these two solutions form a measure-valued solution which may not be approximated by a sequence of weak solutions. As a result, the weak* closure of the set of all weak solutions, considered as parametrized measures, is not equal to the space of all measure-valued solutions. This is in stark contrast with the incompressible Euler equations.


title: Nb Implanted BaO as a Support for Gold Single Atoms

abstract: Using first-principles modelling based on density functional theory we show that oxides implanted with transition metal can act as support for Au single atoms, which are stable against agglomeration. In our previous work we have shown that implanted transition metal, doped in BaO is stable as interstitial in various charge states by transferring the excess charge to an acceptor level close to VBM. Taking Nb as an example we show that single atom Au has its Fermi level close to the VBM of BaO and hence is able to accept charge from the dopant. This charge transfer process between Nb and Au helps Au atoms bind strongly on the doped BaO support. We also show that these charged Au atoms repel each other and prefer to remain atomically dispersed preventing cluster formation. Substitutional doping of transition metals have earlier been reported to bind Au atoms. However, if doped at interstitial sites, they can bind more Au atoms; for example, 5 Au atoms can be anchored per Nb dopant present in BaO interstitial, compared to 3 Au atoms when Nb is doped at substitutional site. This work paves the way for an altogether new technique of stabilizing noble metal single atoms on transition metal doped oxides.


title: Finite-Time Analysis for Double Q-learning

abstract: Although Q-learning is one of the most successful algorithms for finding the best action-value function (and thus the optimal policy) in reinforcement learning, its implementation often suffers from large overestimation of Q-function values incurred by random sampling. The double Q-learning algorithm proposed in~\citet{hasselt2010double} overcomes such an overestimation issue by randomly switching the update between two Q-estimators, and has thus gained significant popularity in practice. However, the theoretical understanding of double Q-learning is rather limited. So far only the asymptotic convergence has been established, which does not characterize how fast the algorithm converges. In this paper, we provide the first non-asymptotic (i.e., finite-time) analysis for double Q-learning. We show that both synchronous and asynchronous double Q-learning are guaranteed to converge to an $\epsilon$-accurate neighborhood of the global optimum by taking $\tilde{\Omega}\left(\left( \frac{1}{(1-\gamma)^6\epsilon^2}\right)^{\frac{1}{\omega}} +\left(\frac{1}{1-\gamma}\right)^{\frac{1}{1-\omega}}\right)$ iterations, where $\omega\in(0,1)$ is the decay parameter of the learning rate, and $\gamma$ is the discount factor. Our analysis develops novel techniques to derive finite-time bounds on the difference between two inter-connected stochastic processes, which is new to the literature of stochastic approximation.


title: Towards decolonising computational sciences

abstract: This article sets out our perspective on how to begin the journey of decolonising computational fields, such as data and cognitive sciences. We see this struggle as requiring two basic steps: a) realisation that the present-day system has inherited, and still enacts, hostile, conservative, and oppressive behaviours and principles towards women of colour (WoC); and b) rejection of the idea that centering individual people is a solution to system-level problems. The longer we ignore these two steps, the more "our" academic system maintains its toxic structure, excludes, and harms Black women and other minoritised groups. This also keeps the door open to discredited pseudoscience, like eugenics and physiognomy. We propose that grappling with our fields' histories and heritage holds the key to avoiding mistakes of the past. For example, initiatives such as "diversity boards" can still be harmful because they superficially appear reformatory but nonetheless center whiteness and maintain the status quo. Building on the shoulders of many WoC's work, who have been paving the way, we hope to advance the dialogue required to build both a grass-roots and a top-down re-imagining of computational sciences -- including but not limited to psychology, neuroscience, cognitive science, computer science, data science, statistics, machine learning, and artificial intelligence. We aspire for these fields to progress away from their stagnant, sexist, and racist shared past into carving and maintaining an ecosystem where both a diverse demographics of researchers and scientific ideas that critically challenge the status quo are welcomed.


title: Visually-Grounded Planning without Vision: Language Models Infer Detailed Plans from High-level Instructions

abstract: The recently proposed ALFRED challenge task aims for a virtual robotic agent to complete complex multi-step everyday tasks in a virtual home environment from high-level natural language directives, such as "put a hot piece of bread on a plate". Currently, the best-performing models are able to complete less than 5% of these tasks successfully. In this work we focus on modeling the translation problem of converting natural language directives into detailed multi-step sequences of actions that accomplish those goals in the virtual environment. We empirically demonstrate that it is possible to generate gold multi-step plans from language directives alone without any visual input in 26% of unseen cases. When a small amount of visual information is incorporated, namely the starting location in the virtual environment, our best-performing GPT-2 model successfully generates gold command sequences in 58% of cases. Our results suggest that contextualized language models may provide strong visual semantic planning modules for grounded virtual agents.


title: Trustworthy Convolutional Neural Networks: A Gradient Penalized-based Approach

abstract: Convolutional neural networks (CNNs) are commonly used for image classification. Saliency methods are examples of approaches that can be used to interpret CNNs post hoc, identifying the most relevant pixels for a prediction following the gradients flow. Even though CNNs can correctly classify images, the underlying saliency maps could be erroneous in many cases. This can result in skepticism as to the validity of the model or its interpretation. We propose a novel approach for training trustworthy CNNs by penalizing parameter choices that result in inaccurate saliency maps generated during training. We add a penalty term for inaccurate saliency maps produced when the predicted label is correct, a penalty term for accurate saliency maps produced when the predicted label is incorrect, and a regularization term penalizing overly confident saliency maps. Experiments show increased classification performance, user engagement, and trust.


title: Abusive Language Detection and Characterization of Twitter Behavior

abstract: In this work, abusive language detection in online content is performed using Bidirectional Recurrent Neural Network (BiRNN) method. Here the main objective is to focus on various forms of abusive behaviors on Twitter and to detect whether a speech is abusive or not. The results are compared for various abusive behaviors in social media, with Convolutional Neural Netwrok (CNN) and Recurrent Neural Network (RNN) methods and proved that the proposed BiRNN is a better deep learning model for automatic abusive speech detection.


title: TEST_POSITIVE at W-NUT 2020 Shared Task-3: Joint Event Multi-task Learning for Slot Filling in Noisy Text

abstract: The competition of extracting COVID-19 events from Twitter is to develop systems that can automatically extract related events from tweets. The built system should identify different pre-defined slots for each event, in order to answer important questions (e.g., Who is tested positive? What is the age of the person? Where is he/she?). To tackle these challenges, we propose the Joint Event Multi-task Learning (JOELIN) model. Through a unified global learning framework, we make use of all the training data across different events to learn and fine-tune the language model. Moreover, we implement a type-aware post-processing procedure using named entity recognition (NER) to further filter the predictions. JOELIN outperforms the BERT baseline by 17.2% in micro F1.


title: An Iterative Least Squares Method for Proton CT Image Reconstruction

abstract: Clinically useful proton Computed Tomography images will rely on algorithms to find the three-dimensional proton stopping power distribution that optimally fits the measured proton data. We present a least squares iterative method with many features to put proton imaging into a more quantitative framework. These include the definition of a unique solution that optimally fits the protons, the definition of an iteration vector that takes into account proton measurement uncertainties, the definition of an optimal step size for each iteration individually, the ability to simultaneously optimize the step sizes of many iterations, the ability to divide the proton data into arbitrary numbers of blocks for parallel processing and use of graphical processing units, and the definition of stopping criteria to determine when to stop iterating. We find that it is possible, for any object being imaged, to provide assurance that the image is quantifiably close to an optimal solution, and the optimization of step sizes reduces the total number of iterations required for convergence. We demonstrate the use of these algorithms on real data.


title: Reduced neural activity during volatile anesthesia compared to TIVA: evidence from a novel EEG signal processing analysis

abstract: Post-operative cognitive decline is a well-known phenomenon and of crucial importance especially in the elderly. General anesthesia can be accomplished by inhalation-based (volatile) or total intravenous anesthesia (TIVA). While their effects on post-operative symptoms have been investigated, little is known about their influence on brain functionalities during the surgery itself. To assess differences 17 patients were divided to receive either volatile anesthesia (n=9), or TIVA (n=8). The level of anesthesia was kept to be equal in both groups. A single bipolar EEG electrode (Neurosteer system) was placed on the participants foreheads. It presented real-time activity and collected their data during the surgery. The dependent variables included frequency bands (delta, theta, alpha, and beta), and three features (VC9, ST4, and A0) previously extracted with the device and provided by Neurosteer. All surgeries were uneventful, and all patients showed bispectral index (BIS) score less than 60. Feature activity under volatile anesthesia (in comparison to TIVA) was significantly lower for the delta, theta and alpha frequency bands and for the three features. Further analysis showed that the largest difference between anesthesia types was for feature A0. The EEG frequency bands and novel brain activity features provide evidence that volatile anesthesia further reduces components of brain activity in comparison to TIVA anesthesia. Specifically, A0, which previously showed a correlation with cognitive decline severity and cognitive load, exhibited the most prominent difference between anesthesia types. Together, this study suggests that measuring brain activity during anesthesia using sensitive features, enables revealing that different anesthesia types may affect brain activity differently, which could affect the recovery from anesthesia, and consequently reduce post-operative cognitive decline


title: CrowdMOT: Crowdsourcing Strategies for Tracking Multiple Objects in Videos

abstract: Crowdsourcing is a valuable approach for tracking objects in videos in a more scalable manner than possible with domain experts. However, existing frameworks do not produce high quality results with non-expert crowdworkers, especially for scenarios where objects split. To address this shortcoming, we introduce a crowdsourcing platform called CrowdMOT, and investigate two micro-task design decisions: (1) whether to decompose the task so that each worker is in charge of annotating all objects in a sub-segment of the video versus annotating a single object across the entire video, and (2) whether to show annotations from previous workers to the next individuals working on the task. We conduct experiments on a diversity of videos which show both familiar objects (aka - people) and unfamiliar objects (aka - cells). Our results highlight strategies for efficiently collecting higher quality annotations than observed when using strategies employed by today's state-of-art crowdsourcing system.


title: There are no exotic ladder surfaces

abstract: It is an open problem to provide a characterization of quasiconformally homogeneous Riemann surfaces. We show that given the current literature, this problem can be broken into four open cases with respect to the topology of the underlying surface. The main result is a characterization in one of the these open cases; in particular, we prove that every quasiconformally homogeneous ladder surface is quasiconformally equivalent to a regular cover of a closed surface (or, in other words, there are no exotic ladder surfaces).


title: Loops in SU(2) and Factorization, II

abstract: In the prequel to this paper, we proved that for a $SU(2,\mathbb C)$ valued loop having the critical degree of smoothness (one half of a derivative in the $L^2$ Sobolev sense), the following statements are equivalent: (1) the Toeplitz and shifted Toeplitz operators associated to the loop are invertible, (2) the loop has a unique triangular factorization, and (3) the loop has a unique root subgroup factorization. This hinges on some Plancherel-esque formulas for determinants of Toeplitz operators. The main point of this report is is to outline a generalization of this result to loops of vanishing mean oscillation, and to discuss some consequences. This generalization hinges on an operator-theoretic factorization of the Toeplitz operators (not simply their determinants).


title: Sonispace: a simulated-space interface for sound design and experimentation

abstract: The world of audio production and design has long been a difficult one to break into, requiring expertise and a working knowledge of the standard digital audio paradigms. This paper describes a novel interface that makes audio production and design more intuitive for novices, using sound-to-space relations that people have learned throughout daily life, such as the roles of barriers and distance in sound perception. The spatial interface for Sonispace allows users to quickly see the relationships between sound-emitting and sound-effecting objects, and to receive audio feedback as they make changes to the space. Algorithms were developed to resemble real-world sonic physics while being efficient enough to provide a user with immediate audio feedback. A prototype of the interface was tested by a group of participants, who confirmed that the software is accessible by novices and that the spatial interface is an engaging way of mixing audio.


title: On the Bieri-Neumann-Strebel-Renz $\Sigma^1$-invariant of even Artin groups

abstract: We calculate the Bieri-Neumann-Strebel-Renz invariant $\Sigma^1(G)$ for even Artin groups $G$ with underlying graph $\Gamma$ such that if there is a closed reduced path in $\Gamma$ with all labels bigger than 2 then the length of such path is always odd. We show that $\Sigma^1(G)^c$ is a rationally defined spherical polyhedron.


title: Improved Battery State Estimation Under Parameter Uncertainty Caused by Aging Using Expansion Measurements

abstract: Accurate tracking of the internal electrochemical states of lithium-ion battery during cycling enables advanced battery management systems to operate the battery safely and maintain high performance while minimizing battery degradation. To this end, techniques based on voltage measurement have shown promise for estimating the lithium surface concentration of active material particles, which is an important state for avoiding aging mechanisms such as lithium plating. However, methods relying on voltage often lead to large estimation errors when the model parameters change during aging. In this paper, we utilize the in-situ measurement of the battery expansion to augment the voltage and develop an observer to estimate the lithium surface concentration distribution in each electrode particle. We demonstrate that the addition of the expansion signal enables us to correct the negative electrode concentration states in addition to the positive electrode. As a result, compared to a voltage only observer, the proposed observer can successfully recover the surface concentration when the electrodes' stoichiometric window changes, which is a common occurrence under aging by loss of lithium inventory. With a 5% shift in the electrodes' stoichiometric window, the results indicate a reduction in state estimation error for the negative electrode surface concentration. Under this simulated aged condition, the voltage based observer had 9.3% error as compared to the proposed voltage and expansion observer which had 0.1% error in negative electrode surface concentration.


title: Using math in physics -- Overview

abstract: The key difference between math as math and math in science is that in science we blend our physical knowledge with our knowledge of math. This blending changes the way we put meaning to math and even to the way we interpret mathematical equations. Learning to think about physics with math instead of just calculating involves a number of general scientific thinking skills that are often taken for granted (and rarely taught) in physics classes. In this paper, I give an overview of my analysis of these additional skills. I propose specific tools for helping students develop these skills in subsequent papers.


title: Methodological reconstruction of historical seismic events from anecdotal accounts of destructive tsunamis: a case study for the great 1852 Banda arc mega-thrust earthquake and tsunami

abstract: We demonstrate the efficacy of a Bayesian statistical inversion framework for reconstructing the likely characteristics of large pre-instrumentation earthquakes from historical records of tsunami observations. Our framework is designed and implemented for the estimation of the location and magnitude of seismic events from anecdotal accounts of tsunamis including shoreline wave arrival times, heights, and inundation lengths over a variety of spatially separated observation locations. As an initial test case we use our framework to reconstruct the great 1852 earthquake and tsunami of eastern Indonesia. Relying on the assumption that these observations were produced by a subducting thrust event, the posterior distribution indicates that the observables were the result of a massive mega-thrust event with magnitude near 8.8 Mw and a likely rupture zone in the north-eastern Banda arc. The distribution of predicted epicentral locations overlaps with the largest major seismic gap in the region as indicated by instrumentally recorded seismic events. These results provide a geologic and seismic context for hazard risk assessment in coastal communities experiencing growing population and urbanization in Indonesia. In addition, the methodology demonstrated here highlights the potential for applying a Bayesian approach to enhance understanding of the seismic history of other subduction zones around the world.


title: Quantum SU$(2|1)$ supersymmetric $\mathbb{C}^N$ Smorodinsky--Winternitz system

abstract: We study quantum properties of SU$(2|1)$ supersymmetric (deformed ${\cal N}=4$, $d=1$ supersymmetric) extension of the superintegrable Smorodinsky--Winternitz system on a complex Euclidian space $\mathbb{C}^N$. The full set of wave functions is constructed and the energy spectrum is calculated. It is shown that SU$(2|1)$ supersymmetry implies the bosonic and fermionic states to belong to separate energy levels, thus exhibiting the "even-odd" splitting of the spectra. The superextended hidden symmetry operators are also defined and their action on SU$(2|1)$ multiplets of the wave functions is given. An equivalent description of the same system in terms of superconformal SU$(2|1,1)$ quantum mechanics is considered and a new representation of the hidden symmetry generators in terms of the SU$(2|1,1)$ ones is found.


title: Hybrid stars in the light of the merging event GW170817

abstract: We study quark-hadron hybrid stars with sharp phase transitions assuming that phase conversions at the interface are slow. Hadronic matter is described by a set of equations of state (EoS) based on the chiral effective field theory and quark matter by a generic bag model. Due to slow conversions at the interface, there is an extended region of stable hybrid stars with central densities above the density of the maximum mass star. We explore systematically the role of the transition pressure and the energy-density jump $\Delta \epsilon$ at the interface on some global properties of hybrid stars. We find that for a given transition pressure, the radius of the last stable hybrid star decreases as $\Delta \epsilon$ raises resulting in a larger extended branch of stable hybrid stars. Contrary to purely hadronic stars, the tidal deformability $\Lambda$ can be either a decreasing or an increasing function of the stellar mass $M$ and for large values of the transition pressure has a very weak dependence on $M$. Finally, we analyze the tidal deformabilities $\Lambda_1$ and $\Lambda_2$ for a binary system with the same chirp mass as GW170817. In the scenario where at least one of the stars in the binary is hybrid, we find that models with low enough transition pressure are inside the $90 \%$ credible region of GW170817. However, these models have maximum masses below $2 \, M_{\odot}$, in disagreement with observations. We also find that the LIGO/Virgo constrain (at $90\%$ level) and the $2 \, M_{\odot}$ requirement can be simultaneously fulfilled in a scenario where all hybrid configurations have masses larger than $1.6 \, M_{\odot}$ and the hadronic EoS is not too stiff, such as several of our hybrid models involving a hadronic EoS of intermediate stiffness. In such scenario hybrid stars may exist in Nature but both objects in GW170817 were hadronic stars.


title: Rectangles, integer vectors and hyperplanes of the hypercube

abstract: We introduce a family of nonnegative integer vectors - primitive vectors - defining hyperplanes of the real affine cube over $C^n:=\{-1,1\}^n$ and study their properties with respect to the rectangles of the cube. As a consequence we give a short proof that, for small dimensions ($n\leq 7$), the real affine cube can be recovered from its signed rectangles and its signed cocircuits complementary of its facets and skew-facets.


title: Compact 200 line MATLAB code for inverse design in photonics by topology optimization: tutorial

abstract: We provide a compact 200 line MATLAB code demonstrating how topology optimization (TopOpt) as an inverse design tool may be used in photonics, targeting the design of two-dimensional dielectric metalenses and a metallic reflector as examples. The physics model is solved using the finite element method, and the code utilizes MATLAB's fmincon algorithm to solve the optimization problem. In addition to presenting the code itself, we briefly discuss a number of extensions and provide the code required to implement some of these. Finally, we demonstrate the superiority of using a gradient-based method compared to a genetic-algorithm-based method (using MATLAB's ga algorithm) for solving inverse design problems in photonics. The MATLAB software is freely available in the paper and may be downloaded from https://www.topopt.mek.dtu.dk.


title: Flow-driven robotic navigation of microengineered endovascular probes

abstract: Minimally invasive medical procedures, such as endovascular catheterization, have drastically reduced procedure time and associated complications. However, many regions inside the body, such as in the brain vasculature, still remain inaccessible due to the lack of appropriate guidance technologies. Here, experimentally and through numerical simulations, we show that tethered ultra-flexible endovascular microscopic probes can be transported through tortuous vascular networks almost effortlessly by harnessing hydrokinetic energy. Dynamic steering at bifurcations is performed by deformation of the probe head using magnetic actuation. We developed an endovascular microrobotic toolkit with a cross-sectional area that is approximately three orders of magnitude smaller than the smallest microcatheter currently available. Our technology has the potential to revolutionize the state-of-the-art practices as it enhances the reachability, reduces the risk of iatrogenic damage, significantly increases the speed of robot-assisted interventions, and enables the deployment of multiple leads simultaneously through a standard needle injection.


title: Price, Volatility and the Second-Order Economic Theory

abstract: We introduce the price probability measure {\eta}(p;t) that defines the mean price p(1;t), mean square price p(2;t), price volatility {\sigma}p2(t)and all price n-th statistical moments p(n;t) as ratio of sums of n-th degree values C(n;t) and volumes U(n;t) of market trades aggregated during certain time interval {\Delta}. The definition of the mean price p(1;t) coincides with definition of the volume weighted average price (VWAP) introduced at least 30 years ago. We show that price volatility {\sigma}p2(t) forecasting requires modeling evolution of the sums of second-degree values C(2;t) and volumes U(2;t). We call this model as second-order economic theory. We use numerical continuous risk ratings as ground for risk assessment of economic agents and distribute agents by risk ratings as coordinates. We introduce continuous economic media approximation of squares of values and volumes of agents trades and their flows aggregated during time interval {\Delta}. We take into account expectations that govern agents trades and introduce aggregated expectations alike to aggregated trades. We derive equations for continuous economic media approximation on the second-degree trades. In the linear approximation we derive mean square price p(2;t) and volatility {\sigma}p2(t) disturbances as functions of the first and second-degree trades disturbances. Description of each next n-th price statistical moment p(n;t) with respect to the unit price measure {\eta}(p;t) depends on sums of n-th degree values C(n;t) and volumes U(n;t) of market trades and hence requires development of the corresponding n-th order economic theory.


title: An Overview on Optimal Flocking

abstract: The study of robotic flocking has received considerable attention in the past twenty years. As we begin to deploy flocking control algorithms on physical multi-agent and swarm systems, there is an increasing necessity for rigorous promises on safety and performance. In this paper, we present an overview the literature focusing on optimization approaches to achieve flocking behavior that provide strong safety guarantees. We separate the literature into cluster and line flocking, and categorize cluster flocking with respect to the system-level objective, which may be realized by a reactive or planning control algorithm. We also categorize the line flocking literature by the energy-saving mechanism that is exploited by the agents. We present several approaches aimed at minimizing the communication and computational requirements in real systems via neighbor filtering and event-driven planning, and conclude with our perspective on the outlook and future research direction of optimal flocking as a field.


title: Learning to swim in potential flow

abstract: Fish swim by undulating their bodies. These propulsive motions require coordinated shape changes of a body that interacts with its fluid environment, but the specific shape coordination that leads to robust turning and swimming motions remains unclear. To address the problem of underwater motion planning, we propose a simple model of a three-link fish swimming in a potential flow environment and we use model-free reinforcement learning for shape control. We arrive at optimal shape changes for two swimming tasks: swimming in a desired direction and swimming towards a known target. This fish model belongs to a class of problems in geometric mechanics, known as driftless dynamical systems, which allow us to analyze the swimming behavior in terms of geometric phases over the shape space of the fish. These geometric methods are less intuitive in the presence of drift. Here, we use the shape space analysis as a tool for assessing, visualizing, and interpreting the control policies obtained via reinforcement learning in the absence of drift. We then examine the robustness of these policies to drift-related perturbations. Although the fish has no direct control over the drift itself, it learns to take advantage of the presence of moderate drift to reach its target.


title: Macroeconomic forecasting through news, emotions and narrative

abstract: This study proposes a new method of incorporating emotions from newspaper articles into macroeconomic forecasts, attempting to forecast industrial production and consumer prices leveraging narrative and sentiment from global newspapers. For the most part, existing research includes positive and negative tone only to improve macroeconomic forecasts, focusing predominantly on large economies such as the US. These works use mainly anglophone sources of narrative, thus not capturing the entire complexity of the multitude of emotions contained in global news articles. This study expands the existing body of research by incorporating a wide array of emotions from newspapers around the world - extracted from the Global Database of Events, Language and Tone (GDELT) - into macroeconomic forecasts. We present a thematic data filtering methodology based on a bi-directional long short term memory neural network (Bi-LSTM) for extracting emotion scores from GDELT and demonstrate its effectiveness by comparing results for filtered and unfiltered data. We model industrial production and consumer prices across a diverse range of economies using an autoregressive framework, and find that including emotions from global newspapers significantly improves forecasts compared to three autoregressive benchmark models. We complement our forecasts with an interpretability analysis on distinct groups of emotions and find that emotions associated with happiness and anger have the strongest predictive power for the variables we predict.


title: Predicting Non Farm Employment

abstract: U.S. Nonfarm employment is considered one of the key indicators for assessing the state of the labor market. Considerable deviations from the expectations can cause market moving impacts. In this paper, the total U.S. nonfarm payroll employment is predicted before the release of the BLS employment report. The content herein outlines the process for extracting predictive features from the aggregated payroll data and training machine learning models to make accurate predictions. Publically available revised employment report by BLS is used as a benchmark. Trained models show excellent behaviour with R2 of 0.9985 and 99.99% directional accuracy on out of sample periods from January 2012 to March 2020. Keywords Machine Learning; Economic Indicators; Ensembling; Regression, Total Nonfarm Payroll


title: Go with the FLOW: Visualizing spatiotemporal dynamics in optical widefield calcium imaging

abstract: Widefield calcium imaging has recently emerged as a powerful experimental technique to record coordinated large-scale brain activity. These measurements present a unique opportunity to characterize spatiotemporal coherent structures that underlie neural activity across many regions of the brain. In this work, we leverage analytic techniques from fluid dynamics to develop a visualization framework that highlights features of flow across the cortex, mapping wave fronts that may be correlated with behavioral events. First, we transform the time series of widefield calcium images into time-varying vector fields using optic flow. Next, we extract concise diagrams summarizing the dynamics, which we refer to as FLOW (flow lines in optical widefield imaging) portraits. These FLOW portraits provide an intuitive map of dynamic calcium activity, including regions of initiation and termination, as well as the direction and extent of activity spread. To extract these structures, we use the finite-time Lyapunov exponent (FTLE) technique developed to analyze time-varying manifolds in unsteady fluids. Importantly, our approach captures coherent structures that are poorly represented by traditional modal decomposition techniques. We demonstrate the application of FLOW portraits on three simple synthetic datasets and two widefield calcium imaging datasets, including cortical waves in the developing mouse and spontaneous cortical activity in an adult mouse.


title: The effects of nonlinear damping on degenerate parametric amplification

abstract: This paper considers the dynamic response of a single degree of freedom system with nonlinear stiffness and nonlinear damping that is subjected to both resonant direct excitation and resonant parametric excitation, with a general phase between the two. This generalizes and expands on previous studies of nonlinear effects on parametric amplification, notably by including the effects of nonlinear damping, which is commonly observed in a large variety of systems, including micro- and nano-scale resonators. Using the method of averaging, a thorough parameter study is carried out that describes the effects of the amplitudes and relative phase of the two forms of excitation. The effects of nonlinear damping on the parametric gain are first derived. The transitions among various topological forms of the frequency response curves, which can include isolae, dual peaks, and loops, are determined, and bifurcation analyses in parameter spaces of interest are carried out. In general, these results provide a complete picture of the system response and allow one to select drive conditions of interest that avoid bistability while providing maximum amplitude gain, maximum phase sensitivity, or a flat resonant peak, in systems with nonlinear damping.