title: Terahertz Radiation from the Dyakonov-Shur Instability of Hydrodynamic Electrons in a Corbino Geometry

abstract: Hydrodynamic electrons flowing through a two-dimensional channel are predicted to undergo a plasma instability above a critical drift velocity. This Dyakonov-Shur (DS) instability terminates as a coherent nonlinear oscillator which shows promise as a source of radiation that could fill the so-called TeraHertz gap. In this work, we study radial flow in a Corbino disk, and demonstrate how the DS instability is substantially enhanced in this geometry, both in terms of a lower critical drift velocity and a higher generated power. Interestingly, hydrodynamic electron flows were recently reported in a graphene sample of this geometry, and our results are therefore directly relevant to current efforts to detect this experimentally elusive phenomenon. The analysis is based on a hydrodynamic approach and features both linearized calculations as well as full numerical simulations of the Navier-Stokes equation.


title: Microlensing Events in the Galactic Plane Using the Zwicky Transient Facility

abstract: Microlensing is a powerful technique to study the Galactic population of "dark" objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such as the Galactic bulge. Microlensing events in the Galactic plane have the advantage of closer proximity and better constrained relative proper motions, leading to better constrained lens mass estimates at the expense of a lower optical depth compared to events towards the Galactic bulge. We use the Zwicky Transient Facility (ZTF) Data Release 5 (DR5) compiled from 2018--2021 to survey the Galactic plane in the region of $|b| < 20^\circ$. We find a total of 60 candidate microlensing events including three that show a strong microlensing parallax effect. The rate of events traces Galactic structure, decreasing exponentially as a function Galactic longitude with scale length $\ell_0 \sim 37^\circ$. On average, we find Einstein timescales of our microlensing events to be about three times as long ($\sim60$ days) compared to those towards the Galactic bulge ($\sim20$ days). This pilot project demonstrates that microlensing towards the Galactic plane shows strong promise for characterization of dark objects within the Galatic disk.


title: Theoretical predictions for inclusive $B\to X_u \tau \bar\nu$ decay

abstract: With the expected large increase in data sets, previously not measured decays will be studied at Belle II. We derive standard model predictions for the $B\to X_u \tau\bar\nu$ decay rate and distributions. The region in the lepton energy spectrum where higher-dimension operators in the local OPE need to be resummed into the $b$-quark light-cone distribution function is a significantly greater fraction of the phase space than for massless leptons. The finite $\tau$ mass has the novel effect of shifting and squeezing how the distribution function enters the lepton energy spectrum. We also derive new predictions for the $\tau$ polarization.


title: Friction pressure on relativistic bubble walls

abstract: During a cosmological first-order phase transition, particles of the plasma crossing the bubble walls can radiate a gauge boson. The resulting pressure cannot be computed perturbatively for large coupling constant and/or large supercooling. We resum the real and virtual emissions at all leading-log orders, both analytically and numerically using a Monte-Carlo simulation. We find that radiated bosons are dominantly soft and that the resulting retarding pressure on relativistic bubble walls is linear both in the Lorentz boost and in the order parameter, up to a log. We further quantitatively discuss IR cut-offs, wall thickness effects, the impact of various approximations entering the calculation, and comment on the fate of radiated bosons that are reflected.


title: High-precision search for dark photon dark matter with the Parkes Pulsar Timing Array

abstract: The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $\epsilon$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario.


title: Constraints on subleading interactions in beta decay Lagrangian

abstract: We discuss the effective field theory (EFT) for nuclear beta decay. The general quark-level EFT describing charged-current interactions between quarks and leptons is matched to the nucleon-level non-relativistic EFT at the O(MeV) momentum scale characteristic for beta transitions. The matching takes into account, for the first time, the effect of all possible beyond-the-Standard-Model interactions at the subleading order in the recoil momentum. We calculate the impact of all the Wilson coefficients of the leading and subleading EFT Lagrangian on the differential decay width in allowed beta transitions. As an example application, we show how the existing experimental data constrain the subleading Wilson coefficients corresponding to pseudoscalar, weak magnetism, and induced tensor interactions. The data display a 3.5 sigma evidence for nucleon weak magnetism, in agreement with the theory prediction based on isospin symmetry.


title: Peculiar Hydrogen-deficient Carbon Stars: Strontium-Rich Stars and the s-Process

abstract: Context: R Coronae Borealis (RCB) variables and their non-variable counterparts, the dustless Hydrogen-Deficient Carbon (dLHdC) stars have been known to exhibit enhanced s-processed material on their surfaces, especially Sr, Y, and Ba. No comprehensive work has been done to explore the s-process in these types of stars, however one particular RCB star, U Aqr, has been under scrutiny for its extraordinary Sr enhancement. Aims: We aim to identify RCB and dLHdC stars that have significantly enhanced Sr abundances, such as U Aqr, and use stellar evolution models to begin to estimate the type of neutron exposure that occurs in a typical HdC star. Methods: We compare the strength of the Sr II 4077 $\AA$ spectral line to Ca II H to identify the new subclass of Sr-rich HdCs. We additionally use the structural and abundance information from existing RCB MESA models to calculate the neutron exposure parameter, $\tau$ Results: We identify six stars in the Sr-rich class. Two are RCBs, and four are dLHdCs. We additionally find that the preferred RCB MESA model has a neutron exposure $\tau$ ~ 0.1 mb$^{-1}$, which is lower than the estimated $\tau$ between 0.15 and 0.6 mb$^{-1}$ for the Sr-rich star U Aqr found in the literature. We find trends in the neutron exposure corresponding to He-burning shell temperature, metallicity, and assumed s-processing site. Conclusions: We have found a sub-class of 6 HdCs known as the Sr-rich class, which tend to lie in the halo, outside the typical distribution of RCBs and dLHdCs. We find that dLHdC stars are more likely to be Sr-rich than RCBs, with an occurrence rate of ~13\% for dLHdCs and ~2\% for RCBs. This is one of the first potential spectroscopic differences between RCBs and dLHdCs, along with dLHdCs having stronger surface abundances of $^{18}$O.


title: Multi-messenger detection prospects of gamma-ray burst afterglows with optical jumps

abstract: Some afterglow light curves of gamma-ray bursts (GRBs) exhibit very complex temporal and spectral features, such as a sudden intensity jump about one hour after the prompt emission in the optical band. We assume that this feature is due to the late collision of two relativistic shells and investigate the corresponding high-energy neutrino emission within a multi-messenger framework, while contrasting our findings with the ones from the classic afterglow model. For a constant density circumburst medium, the total number of emitted neutrinos can increase by about an order of magnitude when an optical jump occurs with respect to the self-similar afterglow scenario. By exploring the detection prospects with the IceCube Neutrino Observatory and future radio arrays such as IceCube-Gen2 radio, RNO-G and GRAND200k, as well as the POEMMA spacecraft, we conclude that the detection of neutrinos with IceCube-Gen2 radio could enable us to constrain the fraction of GRB afterglows with a jump as well as the properties of the circumburst medium. We also investigate the neutrino signal expected for the afterglows of GRB 100621A and a GRB 130427A-like burst with an optical jump. The detection of neutrinos from GRB afterglows could be crucial to explore the yet-to-be unveiled mechanism powering the optical jumps.


title: HeII Ly$\alpha$ Transmission Spikes and Absorption Troughs in Eight High-resolution Spectra Probing the End of HeII Reionization

abstract: We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=\lambda /\Delta \lambda \simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $\Gamma_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background.


title: R Coronae Borealis and dustless Hydrogen-deficient Carbon stars have different oxygen isotope ratios

abstract: R Coronae Borealis (RCB) and dustless Hydrogen-deficient Carbon (dLHdC) stars are believed to be remnants of low mass white dwarf mergers. These supergiant stars have peculiar hydrogen-deficient carbon-rich chemistries and stark overabundances of $^{18}$O. RCB stars undergo dust formation episodes resulting in large-amplitude photometric variations that are not seen in dLHdC stars. Recently, the sample of known dLHdC stars in the Milky Way has more than quintupled with the discovery of 27 new dLHdC stars. We present medium resolution (R$\approx3000$) near-infrared spectra of 20 newly discovered dLHdC stars. We confirm that unlike RCB stars, dLHdC stars do not show strong blueshifted ($>200$ km s$^{-1}$) He I 1.0833 $\mu$m absorption, suggesting the absence of strong, dust-driven winds around them. We also present medium resolution (R$\approx3000-8000$) $K$-band spectra for 47 RCB stars. We measure the $^{16}$O/$^{18}$O ratios of 7 dLHdC and 31 RCB stars that show $^{12}$C$^{16}$O and $^{12}$C$^{18}$O absorption bands, and present the largest sample of values of $^{16}$O/$^{18}$O for dLHdC and RCB stars to date. We find that six of the seven dLHdC stars have $^{16}$O/$^{18}$O $<0.5$, while 26 of the 31 RCB stars have $^{16}$O/$^{18}$O $>1$. We conclude that most dLHdC stars have lower $^{16}$O/$^{18}$O than most RCB stars. This confirms one of the first spectroscopic differences between RCB and dLHdC stars. Our results rule out the existing picture that RCB stars represent an evolved stage of dLHdC stars. Instead, we suggest that whether the white dwarf merger remnant is a dLHdC or RCB star depends on the mass ratios, masses and compositions of the merging white dwarfs.


title: The dawn of a new era for dustless HdC stars with GAIA eDR3

abstract: Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He white dwarf binary systems, otherwise known as the double-degenerate scenario. We searched for new dustless HdC (dLHdC) stars. We used primarily the 2MASS and GAIA eDR3 catalogues to select candidates that were followed-up spectroscopically. We discovered 27 new dLHdC stars, one new RCB star and two new EHe stars. Surprisingly, 20 of the new dLHdC stars share a characteristic of the known dLHdC star HD 148839, having lower atmospheric hydrogen deficiencies. The uncovered population of dLHdC stars exhibit a Bulge-like distribution, like the RCB stars, but show multiple differences from those that indicate they are a different population of HdC stars following its own evolutionary sequence with a fainter luminosity and also a narrow range of effective temperature. We found indication of a current low dust production activity for four of the new dLHdC stars which could be typical RCB stars passing through a transition time. We have evidence for the first time of a large range of absolute magnitudes in the overall population of HdC stars, spanning over 3 mag. In the favoured formation framework, this is explained by a large range in the initial total WD binary mass which leads to a series of evolutionary sequences with distinct maximum brightness and initial temperature. The cold Galactic RCB stars are also noticeably fainter than the Magellanic ones, possibly due to a difference in metallicity resulting in different WD mass ratio. In our Galaxy, there could be as many dLHdC stars as RCB stars.


title: Renormalized $\epsilon$-finite master integrals and their virtues: the three-loop self energy case

abstract: Loop diagram calculations typically rely on reduction to a finite set of master integrals in $4 - 2\epsilon$ dimensions. It has been shown that for any problem, the masters can be chosen so that their coefficients are finite as $\epsilon \rightarrow 0$. I propose a definition of renormalized $\epsilon$-finite master integrals, which incorporate ultraviolet divergence subtractions in a specific way. A key advantage of this choice is that in expressions for physical observables, expansions to positive powers in $\epsilon$ are never needed. As an example, I provide the subtractions for general three-loop self-energy integrals. The differential equations method is used to compute numerically the renormalized $\epsilon$-finite master integrals for arbitrary external momentum invariant, in special cases with internal masses equal to a single scale or zero. These include the ones necessary for the three-loop QCD corrections to the self-energies of the W, Z, and Higgs bosons. In principle, the same method should provide for numerical computation of general three-loop self energies with any masses.


title: Twist-resilient and robust ferroelectric quantum spin Hall insulators driven by van der Waals interactions

abstract: Quantum spin Hall insulators (QSHI) have been proposed to power a number of applications, many of which rely on the possibility to switch on and off the non-trivial topology. Typically this control is achieved through strain or external electric fields, which require energy consumption to be maintained. On the contrary, a non-volatile mechanism would be highly beneficial and could be realized through ferroelectricity if opposite polarization states are associated with different topological phases. While this is not possible in a single ferroelectric material where the two polarization states are related by inversion, the necessary asymmetry could be introduced by combining a ferroelectric layer with another two-dimensional (2D) trivial insulator. Here, by means of first-principles simulations, not only we propose that this is a promising strategy to engineer non-volatile ferroelectric control of topological order in 2D heterostructures, but also that the effect is robust and can survive up to room temperature, irrespective of the weak van der Waals coupling between the layers. We illustrate the general idea by considering a heterostructure made of a well-known ferroelectric material, In$_2$Se$_3$, and a suitably chosen, easily exfoliable trivial insulator, CuI. In one polarization state the system is trivial, while it becomes a QSHI with a robust band gap upon polarization reversal. Remarkably, the topological band gap is mediated by the interlayer hybridization and allows to maximise the effect of intralayer spin-orbit coupling, promoting a robust ferroelectric topological phase that could not exist in monolayer materials and is resilient against relative orientation and lattice matching between the layers.


title: Chase-escape on the configuration model

abstract: Chase-escape is a competitive growth process in which red particles spread to adjacent empty sites according to a rate-$\lambda$ Poisson process while being chased and consumed by blue particles according to a rate-$1$ Poisson process. Given a growing sequence of finite graphs, the critical rate $\lambda_c$ is the largest value of $\lambda$ for which red fails to reach a positive fraction of the vertices with high probability. We provide a conjecturally sharp lower bound and an implicit upper bound on $\lambda_c$ for supercritical random graphs sampled from the configuration model with independent and identically distributed degrees with finite second moment. We additionally show that the expected number of sites occupied by red undergoes a phase transition and identify the location of this transition.


title: Dirac States in an Inclined Two-Dimensional Su-Schrieffer-Heeger Model

abstract: We propose to realize Dirac states in an inclined two-dimensional Su-Schrieffer-Heeger model on a square lattice. We show that a pair of Dirac points protected by space-time inversion symmetry appear in the semimetal phase. The locations of these Dirac points are not pinned to any high-symmetry points of the Brillouin zone but are tunable through parameter modulations. Interestingly, the merging of two Dirac points undergoes a topological phase transition that leads to either a weak topological insulator or a nodal-line semimetal. We provide a systematic analysis of these topological phases from both bulk and boundary perspectives combined with symmetry arguments. We also discuss feasible experimental platforms to realize our model.


title: The Visibility Graphs of Correlated Time Series Violate the Barthelemy's Conjecture for Degree and Betweenness Centralities

abstract: The problem of betweenness centrality remains a fundamental unsolved problem in complex networks. After a pioneering work by Barthelemy, it has been well-accepted that the maximal betweenness-degree ($b$-$k$) exponent for scale-free (SF) networks is $\eta_{\text{max}}=2$, belonging to scale-free trees (SFTs), based on which one concludes $\delta\ge\frac{\gamma+1}{2}$, where $\gamma$ and $\delta$ are the scaling exponents of the distribution functions of the degree and betweenness centrality, respectively. Here we present evidence for violation of this conjecture for SF visibility graphs (VGs). To this end, we consider the VG of three models: two-dimensional (2D) Bak-Tang-Weisenfeld (BTW) sandpile model, 1D fractional Brownian motion (FBM) and, 1D Levy walks, the two later cases are controlled by the Hurst exponent $H$ and step-index $\alpha$, respectively. Specifically, for the BTW model and FBM with $H\lesssim 0.5$, $\eta$ is greater than $2$, and also $\delta<\frac{\gamma+1}{2}$ for the BTW model, while Barthelemy's conjecture remains valid for the Levy process. We argue that this failure of Barthelemy's conjecture is due to large fluctuations in the scaling $b$-$k$ relation resulting in the violation of hyperscaling relation $\eta=\frac{\gamma-1}{\delta-1}$ and emergent anomalous behaviors for the BTW model and FBM. A super-universal behavior is found for the distribution function for a generalized degree function identical to the Barabasi-Albert network model.


title: The DES view of the Eridanus supervoid and the CMB Cold Spot

abstract: The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $\Lambda$CDM model, only about 10-20$\%$ of the observed temperature depression can be accounted for via its Integrated Sachs-Wolfe imprint. However, $R\gtrsim100~h^{-1}\mathrm{Mpc}$ supervoids elsewhere in the sky have shown ISW imprints $A_{\mathrm{ISW}}\approx5.2\pm1.6$ times stronger than expected from $\Lambda$CDM ($A_{\mathrm{ISW}}=1$), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant under-density in the Cold Spot's direction at $z<0.2$. We also show, with $\mathrm{S/N}\gtrsim5$ significance, that the Eridanus supervoid appears as the most prominent large-scale under-density in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about $30\%$ lower than expected from similar peaks found in N-body simulations based on the standard $\Lambda$CDM model with parameters $\Omega_{\rm m} = 0.279$ and $\sigma_8 = 0.82$. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region.


title: Improving and Maximal Inequalities for Primes in Progressions

abstract: Assume that $ y < N$ are integers, and that $ (b,y) =1$. Define an average along the primes in a progression of diameter $ y$, given by integer $ (b,y)=1 $. \begin{align*} A_{N,y,b} := \frac{\phi (y)}{N} \sum _{\substack{n <N\\n\equiv b\pmod{y}}} \Lambda (n) f(x-n) \end{align*} Above, $\Lambda $ is the von Mangoldt function and $\phi $ is the totient function. We establish improving and maximal inequalities for these averages. These bounds are uniform in the choice of progression. For instance, for $ 1< r < \infty $ there is an integer $N _{y, r}$ so that \begin{align*} \lVert \sup _{N>N _{y,r}} \lvert A_{N,y,b} f \rvert \rVert_{r}\ll \lVert f\rVert_{r}. \end{align*} The implied constant is only a function of $ r$. The uniformity over progressions imposes several novel elements on the proof.


title: Conservative and Adaptive Penalty for Model-Based Safe Reinforcement Learning

abstract: Reinforcement Learning (RL) agents in the real world must satisfy safety constraints in addition to maximizing a reward objective. Model-based RL algorithms hold promise for reducing unsafe real-world actions: they may synthesize policies that obey all constraints using simulated samples from a learned model. However, imperfect models can result in real-world constraint violations even for actions that are predicted to satisfy all constraints. We propose Conservative and Adaptive Penalty (CAP), a model-based safe RL framework that accounts for potential modeling errors by capturing model uncertainty and adaptively exploiting it to balance the reward and the cost objectives. First, CAP inflates predicted costs using an uncertainty-based penalty. Theoretically, we show that policies that satisfy this conservative cost constraint are guaranteed to also be feasible in the true environment. We further show that this guarantees the safety of all intermediate solutions during RL training. Further, CAP adaptively tunes this penalty during training using true cost feedback from the environment. We evaluate this conservative and adaptive penalty-based approach for model-based safe RL extensively on state and image-based environments. Our results demonstrate substantial gains in sample-efficiency while incurring fewer violations than prior safe RL algorithms. Code is available at: https://github.com/Redrew/CAP


title: Quantifying the Excitonic Static Disorder in Organic Semiconductors

abstract: Organic semiconductors are disordered molecular solids and as a result, their internal charge dynamics and ultimately, the performance of the optoelectronic devices they constitute, are governed by energetic disorder. To ascertain how energetic disorder impacts charge generation, exciton transport, charge transport, and the performance of organic semiconductor devices, an accurate approach is first required to measure this critical parameter. In this work, we show that the static disorder has no relation with the so-called Urbach energy in organic semiconductors. Instead, it can be obtained from photovoltaic external quantum efficiency spectra at wavelengths near the absorption onset. We then present a detailed methodology, alongside a computational framework, for quantifying the static energetic disorder associated with singlet excitons. Moreover, the role of optical interference in this analysis is considered to achieve a high-accuracy quantification. Finally, the excitonic static disorder was quantified in several technologically-relevant donor-acceptor blends, including high-efficiency PM6:Y6.


title: EP-PINNs: Cardiac Electrophysiology Characterisation using Physics-Informed Neural Networks

abstract: Accurately inferring underlying electrophysiological (EP) tissue properties from action potential recordings is expected to be clinically useful in the diagnosis and treatment of arrhythmias such as atrial fibrillation, but it is notoriously difficult to perform. We present EP-PINNs (Physics-Informed Neural Networks), a novel tool for accurate action potential simulation and EP parameter estimation, from sparse amounts of EP data. We demonstrate, using 1D and 2D in silico data, how EP-PINNs are able to reconstruct the spatio-temporal evolution of action potentials, whilst predicting parameters related to action potential duration (APD), excitability and diffusion coefficients. EP-PINNs are additionally able to identify heterogeneities in EP properties, making them potentially useful for the detection of fibrosis and other localised pathology linked to arrhythmias. Finally, we show EP-PINNs effectiveness on biological in vitro preparations, by characterising the effect of anti-arrhythmic drugs on APD using optical mapping data. EP-PINNs are a promising clinical tool for the characterisation and potential treatment guidance of arrhythmias.


title: A Tutorial on Bayesian Data Assimilation

abstract: This tutorial provides a broad introduction to Bayesian data assimilation that will be useful to practitioners, in interpreting algorithms and results, and for theoretical studies developing novel schemes with an understanding of the rich history of geophysical data assimilation and its current directions. The simple case of data assimilation in a 'perfect' model is primarily discussed for pedagogical purposes. Some mathematical results are derived at a high-level in order to illustrate key ideas about different estimators. However, the focus of this work is on the intuition behind these methods, where more formal and detailed treatments of the data assimilation problem can be found in the various references. In surveying a variety of widely used data assimilation schemes, the key message of this tutorial is how the Bayesian analysis provides a consistent framework for the estimation problem and how this allows one to formulate its solution in a variety of ways to exploit the operational challenges in the geosciences.


title: Mode solutions to the wave equation on a rotating cosmic string background

abstract: A static rotating cosmic string metric is singular along a timelike line and fails to be globally hyperbolic; these features make it difficult to solve the wave equation by conventional energy methods. Working on a single angular mode at a time, we use microlocal methods to construct forward parametrices for wave and Klein--Gordon equations on such backgrounds.


title: Disorder, Low-Energy Excitations, and Topology in the Kitaev Spin Liquid

abstract: The Kitaev model is a fascinating example of an exactly solvable model displaying a spin-liquid ground state in two dimensions. However, deviations from the original Kitaev model are expected to appear in real materials. In this Letter, we investigate the fate of Kitaev's spin liquid in the presence of disorder -- bond defects or vacancies -- for an extended version of the model. Considering static flux backgrounds, we observe a power-law divergence in the low-energy limit of the density of states with a nonuniversal exponent. We link this power-law distribution of energy scales to weakly coupled droplets inside the bulk, in an uncanny similarity to the Griffiths phase often present in the vicinity of disordered quantum phase transitions. If time-reversal symmetry is broken, we find that power-law singularities are tied to the destruction of the topological phase of the Kitaev model in the presence of bond disorder alone. However, there is a transition from this topologically trivial phase with power-law singularities to a topologically nontrivial one for weak to moderate site dilution. Therefore, diluted Kitaev materials are potential candidates to host Kitaev's chiral spin-liquid phase.


title: Resonant response in tunable metasurface based on crossed all-dielectric gratings

abstract: We theoretically demonstrate mechanical tuning of the spectral response of a Mie-like resonant dielectric metasurface consisting of crossed all-dielectric strip gratings. We use two array of parallel dielectric strips superimposed on each other to mechanically alter the crossing angle between them thereby changing the position of the maximum of the spectral response in a wide frequency range. The metasurface structure is considered as a doubly-periodic system in an oblique mesh. This mechanically controllable metasurface can be perspective in developing of a flexible polarization-sensitive optical devices.


title: Learning to Retrieve Passages without Supervision

abstract: Dense retrievers for open-domain question answering (ODQA) have been shown to achieve impressive performance by training on large datasets of question-passage pairs. In this work we ask whether this dependence on labeled data can be reduced via unsupervised pretraining that is geared towards ODQA. We show this is in fact possible, via a novel pretraining scheme designed for retrieval. Our "recurring span retrieval" approach uses recurring spans across passages in a document to create pseudo examples for contrastive learning. Our pretraining scheme directly controls for term overlap across pseudo queries and relevant passages, thus allowing to model both lexical and semantic relations between them. The resulting model, named Spider, performs surprisingly well without any labeled training examples on a wide range of ODQA datasets. Specifically, it significantly outperforms all other pretrained baselines in a zero-shot setting, and is competitive with BM25, a strong sparse baseline. Moreover, a hybrid retriever over Spider and BM25 improves over both, and is often competitive with DPR models, which are trained on tens of thousands of examples. Last, notable gains are observed when using Spider as an initialization for supervised training.


title: Integrated Neighborhood Colorings of Graphs

abstract: The idea that those different from you are "unfriendly" is captured in the definition of unfriendly 2-colorings in graph theory in a paper by Aharoni, Milner and Prikry, where they prove that every finite graph has an unfriendly coloring. We give a more general definition for all n>1, that we call "integrated" rather than "unfriendly." Then we prove that every finite graph has an integrated n-coloring, n>1. We then give some applications to various graph coloring problems and to some max-cut problems.


title: Aymptotics of eigenvalues of the Neumann-Poincar'e operator in 3D elasticity

abstract: We consider the Neumann-Poincar'e (double layer potential) operator in 3D elasicity on a smooth closed surface. Its essential spectrum consists of 3 points. We find the asymptotics of sequences of eigenvalues converging to these three pounts. They are expressed via the eigenvalue distribution for a matrix eigenvalue problem depending on the points of the cotangent bundle of the surface. For the two-sided asymptotics for the distribution of the the union of sequences of eigenvalues converging to the point of the essential spectrum from both sides, an expression is found via the Euler characteristics and the Willmore energy of the surface.


title: Representing Inferences and their Lexicalization

abstract: We have recently begun a project to develop a more effective and efficient way to marshal inferences from background knowledge to facilitate deep natural language understanding. The meaning of a word is taken to be the entities, predications, presuppositions, and potential inferences that it adds to an ongoing situation. As words compose, the minimal model in the situation evolves to limit and direct inference. At this point we have developed our computational architecture and implemented it on real text. Our focus has been on proving the feasibility of our design.


title: A note on Singular integral

abstract: In this paper, we prove that for $\frac{n}{2}+\frac{1}{4}<\alpha \leq\frac{n+1}{2} $, the convolution operator $$S_{\alpha} f(x)=\int_{|y| \geq 1} f(x-y)\left(|y|^{2}-1\right)^{-\alpha} d y$$ is bounded from $L^p$ to $L^q$ for certain values of $p$ and $q$.


title: Knots, minimal surfaces and J-holomorphic curves

abstract: Let $K$ be a knot in the 3-sphere, viewed as the ideal boundary of hyperbolic 4-space $\mathbb{H}^4$. We prove that the number of minimal discs in $\mathbb{H}^4$ with ideal boundary $K$ is a knot invariant. I.e.\ the number is finite and doesn't change under isotopies of $K$. In fact this gives a family of knot invariants, indexed by an integer describing the extrinsic topology of how the disc sits in $\mathbb{H}^4$. These invariants can be seen as Gromov--Witten invariants counting $J$-holomorphic discs in the twistor space $Z$ of $\mathbb{H}^4$. Whilst Gromov--Witten theory suggests the general scheme for defining the invariants, there are substantial differences in how this must be carried out in our situation. These are due to the fact that the geometry of both $\mathbb{H}^4$ and $Z$ becomes singular at infinity, and so the $J$-holomorphic curve equation is degenerate, rather than elliptic, at the boundary. This means that both the Fredholm and compactness arguments involve completely new features.


title: Numerical optimization of amplitude-modulated pulses in microwave-driven entanglement generation

abstract: Microwave control of trapped ions can provide an implementation of high-fidelity two-qubit gates free from errors induced by photon scattering. Furthermore, microwave conductors may be embedded into a scalable trap structure, providing the chip-level integration of control that is desirable for scaling. Recent developments have demonstrated how amplitude modulation of the gate drive can permit a two-qubit entangling operation to become robust against motional mode noise and other experimental imperfections. Here, we discuss a method for the numerical optimization of the microwave pulse envelope to produce gate pulses with improved resilience, faster operation and higher energy efficiency.


title: Finite subgroups of automorphisms of K3 surfaces

abstract: We give a complete classification of finite subgroups of automorphisms of K3 surfaces up to deformation. The classification is in terms of Hodge theoretic data associated to certain conjugacy classes of finite subgroups of the orthogonal group of the K3 lattice. The moduli theory of K3 surfaces, in particular the surjectivity of the period map and the strong Torelli theorem allow us to interpret this datum geometrically. Our approach is computer aided and involves hermitian lattices over number fields.


title: Dark matter component in hadronic models with short-range correlations

abstract: A relativistic mean field hadronic model with a dark matter (DM) particle coupled to nucleons including short-range correlations (SRC) is applied to study neutron stars (NS). The lightest neutralino is chosen as the dark particle candidate, which interacts with nucleons by the exchange of Higgs bosons. A detailed thermodynamical analysis shows that the contribution of the DM fermions to the energy density of the matter composed by these particles and nucleons is completely dominated by the DM kinetic terms. The model reproduces satisfactorily the constraints on the mass-radius diagram obtained from the analysis of the combined data from the NICER mission, LIGO and Virgo collaborations, and mass measurements from radio observations. We show that the SRC balance the reduction of the neutron star mass due to the DM component, and because of that the model is able to present more massive NS. We also present a study of the effect, in the NS mass-radius profiles, of the uncertainties in some bulk parameters related to the hadronic sector. We find that it is possible to generate parametrizations, with DM content, compatible with the recent astrophysical constraints and with the uncertainty in the symmetry energy slope obtained from the results reported by the updated Lead Radius EXperiment (PREX-2).


title: Deterministic and Stochastic in-host Tuberculosis Models for Bacterium-directed and Host-directed Therapy Combination

abstract: Mycobacterium tuberculosis infection can involve all immune system components and can result in different disease outcomes. The antibiotic TB drugs require strict adherence to prevent both disease relapse and mutation of drug- and multidrug-resistant strains. To overcome the constraints of pathogen-directed therapy, host-directed therapy has attracted more attention in recent years as an adjunct therapy to enhance host immunity to fight against this intractable pathogen. The goal of this paper is to investigate in-host tuberculosis models to provide insights into therapy development. Focusing on therapy-targeting parameters, the parameter regions for different disease outcomes are identified from an established ODE model. Interestingly, the ODE model also demonstrates that the immune responses can both benefit and impede disease progression, depending on the number of bacteria engulfed and released by macrophages. We then develop two It\^{o} SDE models, which consider the impact of demographic variations at the cellular level and environmental variations during therapies along with demographic variations. The SDE model with demographic variation suggests that stochastic fluctuations at the cellular level have significant influences on (1) the T-cell population in all parameter regions, (2) the bacterial population when parameters located in the region with multiple disease outcomes, and (3) the uninfected macrophage population in the parameter region representing active disease. Further, considering environmental variations from therapies, the second SDE model suggests that disease progression can slow down if therapies (1) can have fast return rates and (2) can bring parameter values into the disease clearance regions.


title: Scatterbrained: A flexible and expandable pattern for decentralized machine learning

abstract: Federated machine learning is a technique for training a model across multiple devices without exchanging data between them. Because data remains local to each compute node, federated learning is well-suited for use-cases in fields where data is carefully controlled, such as medicine, or in domains with bandwidth constraints. One weakness of this approach is that most federated learning tools rely upon a central server to perform workload delegation and to produce a single shared model. Here, we suggest a flexible framework for decentralizing the federated learning pattern, and provide an open-source, reference implementation compatible with PyTorch.


title: Decomposing the Deep: Finding Class Specific Filters in Deep CNNs

abstract: Interpretability of Deep Neural Networks has become a major area of exploration. Although these networks have achieved state of the art accuracy in many tasks, it is extremely difficult to interpret and explain their decisions. In this work we analyze the final and penultimate layers of Deep Convolutional Networks and provide an efficient method for identifying subsets of features that contribute most towards the network's decision for a class. We demonstrate that the number of such features per class is much lower in comparison to the dimension of the final layer and therefore the decision surface of Deep CNNs lies on a low dimensional manifold and is proportional to the network depth. Our methods allow to decompose the final layer into separate subspaces which is far more interpretable and has a lower computational cost as compared to the final layer of the full network.


title: Searching for axion-like particles with data scouting at ATLAS and CMS

abstract: We investigate the physics case for a dedicated trigger-level analysis for very low mass diphoton resonances at ATLAS and CMS and introduce a new photon isolation criterion. This greatly increases the acceptance for light particles at the expense of a somewhat larger Standard Model background. We show how such an analysis would likely yield new experimental coverage for axion-like particles for masses below 15 GeV, bridging the gap with the region covered by flavor factories.


title: Polarizable embedding complex polarization propagator in four- and two-component frameworks

abstract: Explicit embedding methods combined with the complex polarization propagator (CPP) enable modeling of spectroscopy for increasingly complex systems with a high density of states. We present the first derivation and implementation of the CPP in four- and exact two-component (X2C) polarizable embedding (PE) frameworks. We denote the developed methods PE-4c-CPP and PE-X2C-CPP, respectively. We illustrate the methods by estimating the solvent effect on UV-vis and X-ray atomic absorption (XAS) spectra of [Rh(H2O)6]3+ and [Ir(H2O)6]3+ immersed in aqueous solution. We moreover estimate solvent effects on UV-vis spectra of a platinum complex that can be photo-chemically activated (in water) to kill cancer cells. Our results clearly show that inclusion of the environment is required: UV-vis and (to a lesser degree) XAS spectra can become qualitatively different from vacuum calculations. Comparison of PE-4c-CPP and PE-X2C-CPP methods shows that X2C essentially reproduces the solvent effect obtained with the 4c methods.


title: Partial-wave approach to the Stark resonance problem of the water molecule

abstract: A partial-wave method is developed to deal with small molecules dominated by a central atom as an extension of earlier single-center methods. In particular, a model potential for the water molecule is expanded over a basis of spherical harmonics. A finite element method is employed to generate local polynomial functions in subintervals over a finite range for the radial variable. The angular parts of the system are represented by spherical harmonics. The problem of Stark resonances is treated with the exterior complex scaling method which incorporates a wavefunction discontinuity at the scaling radius. The resultant non-hermitian matrix eigenvalue problem yields resonance positions and widths (decay rates). We present these DC Stark shifts and exponential decay rates for the valence orbitals $1b_1$, $3a_1$, and the bonding orbital $1b_2$. Furthermore, comparison is made with total molecular decay rates and DC shifts obtained recently within the Hartree-Fock and coupled-cluster approaches.


title: Autonomous Navigation System from Simultaneous Localization and Mapping

abstract: This paper presents the development of a Simultaneous Localization and Mapping (SLAM) based Autonomous Navigation system. The motivation for this study was to find a solution for navigating interior spaces autonomously. Interior navigation is challenging as it can be forever evolving. Solving this issue is necessary for multitude of services, like cleaning, the health industry, and in manufacturing industries. The focus of this paper is the description of the SLAM-based software architecture developed for this proposed autonomous system. A potential application of this system, oriented to a smart wheelchair, was evaluated. Current interior navigation solutions require some sort of guiding line, like a black line on the floor. With this proposed solution, interiors do not require renovation to accommodate this solution. The source code of this application has been made open source so that it could be re-purposed for a similar application. Also, this open-source project is envisioned to be improved by the broad open-source community upon past its current state.


title: Designing a Magnetic Measurement Data Acquisition and Control System with Reuse in Mind: A Rotating Coil System Example

abstract: Accelerator magnet test facilities frequently need to measure different magnets on differently equipped test stands and with different instrumentation. Designing a modular and highly reusable system that combines flexibility built-in at the architectural level as well as on the component level addresses this need. Specification of the backbone of the system, with the interfaces and dataflow for software components and core hardware modules, serves as a basis for building such a system. The design process and implementation of an extensible magnetic measurement data acquisition and control system are described, including techniques for maximizing the reuse of software. The discussion is supported by showing the application of this methodology to constructing two dissimilar systems for rotating coil measurements, both based on the same architecture and sharing core hardware modules and many software components. The first system is for production testing 10 m long cryo-assemblies containing two MQXFA quadrupole magnets for the high-luminosity upgrade of the Large Hadron Collider and the second for testing IQC conventional quadrupole magnets in support of the accelerator system at Fermilab.


title: Limit of connected multigraph with fixed degree sequence

abstract: Motivated by the scaling limits of the connected components of the configuration model, we study uniform connected multigraphs with fixed degree sequence $\mathcal{D}$ and with surplus $k$. We call those random graphs $(\mathcal{D},k)$-graphs. We prove that, for every $k\in \mathbb N$, under natural conditions of convergence of the degree sequence, ($\mathcal{D},k)$-graphs converge toward either $(\mathcal{P},k)$-graphs or $(\Theta,k)$-ICRG (inhomogeneous continuum random graphs). We prove similar results for $(\mathcal{P},k)$-graphs and $(\Theta,k)$-ICRG, which have applications to multiplicative graphs. Our approach relies on two algorithms, the cycle-breaking algorithm, and the stick-breaking construction of $\mathcal{D}$-tree that we introduced in a recent paper arXiv:2110.03378. From those algorithms we deduce a biased construction of $(\mathcal{D},k)$-graph, and we prove our results by studying this bias.


title: Predictive model of fermionic dark matter halos with a quantum core and an isothermal atmosphere

abstract: We develop a thermodynamical model of fermionic dark matter halos at finite temperature. Statistical equilibrium states may be justified by a process of violent collisionless relaxation in the sense of Lynden-Bell or from a collisional relaxation of nongravitational origin if the fermions are self-interacting. The most probable state (maximum entropy state) generically has a "core-halo" structure with a quantum core (fermion ball) surrounded by an isothermal atmosphere. The quantum core is equivalent to a polytrope of index $n=3/2$. The Pauli exclusion principle creates a quantum pressure that prevents gravitational collapse and solves the core-cusp problem of the cold dark matter model. The isothermal atmosphere (which is similar to the NFW profile of cold dark matter) accounts for the flat rotation curves of the galaxies at large distances. We numerically solve the equation of hydrostatic equilibrium with the Fermi-Dirac equation of state and determine the density profiles and rotation curves of fermionic dark matter halos.


title: Cosmological test of local position invariance from the asymmetric galaxy clustering

abstract: The local position invariance (LPI) is one of the three major pillars of Einstein equivalence principle, ensuring the space-time independence on the outcomes of local experiments. The LPI has been tested by measuring the gravitational redshift effect in various depths of gravitational potentials. We propose a new cosmological test of the LPI by observing the asymmetry in the cross-correlation function between different types of galaxies, which predominantly arises from the gravitational redshift effect induced by the gravitational potential of halos at which the galaxies reside. We show that the ongoing and upcoming galaxy surveys can give a fruitful constraint on the LPI-violating parameter, $\alpha$, at distant universes (redshift $z\sim0.1-1.8$) over the cosmological scales (separation $s\sim5-10\, {\rm Mpc}/h$) that have not yet been explored, finding that the expected upper limit on $\alpha$ can reach $0.03$.


title: Fixed points, descents, and inversions in parabolic double cosets of the symmetric group

abstract: We consider statistics on permutations chosen uniformly at random from fixed parabolic double cosets of the symmetric group. We show that the distribution of fixed points is asymptotically Poisson and establish central limit theorems for the distribution of descents and inversions. Our proofs use Stein's method with size-bias coupling and dependency graphs, which also gives convergence rates for our distributional approximations. As applications of our size-bias coupling and dependency graph constructions, we obtain concentration of measure results on the number of fixed points, descents, and inversions.


title: Searching for Diamagnetic Blob Accretion in the 74 day K2 Observation of V2400 Ophiuchi

abstract: Since its discovery in 1995, V2400 Ophiuchi (V2400 Oph) has stood apart from most known intermediate polar cataclysmic variables due to its proposed magnetic field strength (9-27 MG) and disk-less accretion. To date, the exact accretion mechanism of the system is still unknown, and standard accretion models fail to accurately predict the peculiar behavior of its lightcurve. We present the K2 Campaign~11 light curve of V2400 Oph recording 74.19 days of photometric data cadenced at 1 minute. The light curve is dominated by aperiodic flickering and quasi-periodic oscillations, which make the beat and spin signals inconspicuous on short timescales. Notably, a log-log full power spectrum shows a break frequency at $\sim10^2$ cycles~d$^{-1}$ similar to some disk-fed systems. Through power spectral analysis, the beat and spin periods are measured as $1003.4\pm0.2$ seconds and $ 927.7\pm 0.1$ seconds respectively. A power spectrum of the entire K2 observation demonstrates beat period dominance. However, time-resolved power spectra reveals a strong dependence between observation length and the dominant frequency of the light curve. For short observations (2-12 hrs) the beat, spin, or first beat harmonic can be observed as the dominant periodic signal. Such incoherence and variability indicate a dynamical accretion system more complex than current intermediate polar theories can explain. We propose that a diamagnetic blob accretion model may serve as a plausible explanation for the accretion mechanism.


title: Localized big bang stability for the Einstein-scalar field equations

abstract: We prove the nonlinear stability in the contracting direction of Friedmann-Lema\^itre-Robertson-Walker (FLRW) solutions to the Einstein-scalar field equations in $n\geq 3$ spacetime dimensions that are defined on spacetime manifolds of the form $(0,t_0]\times \mathbb{T}^{n-1}$, $t_0>0$. Stability is established under the assumption that the initial data is \textit{synchronized}, which means that on the initial hypersurface $\Sigma= \{t_0\}\times \mathbb{T}^{n-1}$ the scalar field $\tau= \exp\bigl(\sqrt{\frac{2(n-2)}{n-1}}\phi\bigr) $ is constant, that is, $\Sigma=\tau^{-1}(\{t_0\})$. As we show that all initial data sets that are sufficiently close to FRLW ones can be evolved via the Einstein-scalar field equation into new initial data sets that are \textit{synchronized}, no generality is lost by this assumption. By using $\tau$ as a time coordinate, we establish that the perturbed FLRW spacetime manifolds are of the form $M = \bigcup_{t\in (0,t_0]}\tau^{-1}(\{t\})\cong (0,t_0]\times \mathbb{T}^{n-1}$, the perturbed FLRW solutions are asymptotically pointwise Kasner as $\tau \searrow 0$, and a big bang singularity, characterised by the blow up of the scalar curvature, occurs at $\tau=0$. An important aspect of our past stability proof is that we use a hyperbolic gauge reduction of the Einstein-scalar field equations. As a consequence, all of the estimates used in the stability proof can be localized and we employ this property to establish a corresponding localized past stability result for the FLRW solutions.


title: Toward a unified theory for the Hadley cell descending and ascending edges

abstract: We present theories for the latitudinal extents of both Hadley cells throughout the annual cycle by combining our recent scaling for the ascending edge latitude (Hill et al. 2021) with the uniform Rossby number (Ro), baroclinic instability-based theory for the poleward, descending edge latitudes of Kang and Lu 2012. The resulting analytic expressions for all three Hadley cell edges are predictive except for diagnosed values of Ro and two proportionality constants. The theory captures the climatological annual cycle of the ascending and descending edges in an Earth-like simulation in an idealized aquaplanet general circulation model (GCM), provided the descending edge prediction is lagged by one month. In simulations in this and two other idealized GCMs with varied planetary rotation rate ($\Omega$), the winter, descending edge of the solsticial, cross-equatorial Hadley cell scales approximately as $\Omega^{-1/2}$ and the summer, ascending edge as $\Omega^{-2/3}$, both in accordance with our theory.


title: Polarized Structure Function $\sigma_{LT'}$ from $\pi^0 p$ Electroproduction Data in the Resonance Region at $0.4$ GeV$^2 < Q^2 < 1.0$ GeV$^2$

abstract: The first results on the $\sigma_{LT'}$ structure function in exclusive $\pi^0p$ electroproduction at invariant masses of the final state of 1.5 GeV $<$ $W$ $<$ 1.8 GeV and in the range of photon virtualities 0.4 GeV$^2 < Q^2 < 1.0$ GeV$^2$ were obtained from data on beam spin asymmetries and differential cross sections measured with the CLAS detector at Jefferson Lab. The Legendre moments determined from the $\sigma_{LT'}$ structure function have demonstrated sensitivity to the contributions from the nucleon resonances in the second and third resonance regions. These new data on the beam spin asymmetries in $\pi^0p$ electroproduction extend the opportunities for the extraction of the nucleon resonance electroexcitation amplitudes in the mass range above 1.6 GeV.


title: Emergence and Manipulation of non-equilibrium Yu-Shiba-Rusinov states

abstract: The experimental advances in the study of time-dependent phenomena has opened a new path to investigating the complex electronic structure of strongly correlated and topological materials. Yu-Shiba-Rusinov (YSR) states induced by magnetic impurities in $s$-wave superconductors provide an ideal candidate system to study the response of a system to time-dependent manipulations of the magnetic environment. Here, we show that by imposing a time-dependent change in the magnetic exchange coupling, by changing the relative alignment of magnetic moments in an impurity dimer, or through a periodic drive of the impurity moment, one can tune the system through a time-dependent quantum phase transition, in which the system undergoes a transition from a singlet to a doublet ground state. We show that the electronic response of the system to external perturbations can be imaged through the time-dependent differential conductance, $dI(t)/dV$, which, in analogy to the equilibrium case, is proportional to a non-equilibrium local density of states. Our results open the path to visualizing the response of complex quantum systems to time-dependent external perturbations.


title: Multi-wavelength quasi-periodic pulsations in a stellar superflare

abstract: We present the first multi-wavelength simultaneous detection of QPP in a superflare (more than a thousand times stronger than known solar flares) on a cool star, in soft X-rays (SXR, with XMM-Newton) and white light (WL, with Kepler). It allowed for the first-ever analysis of oscillatory processes in a stellar flare simultaneously in thermal and non-thermal emissions, conventionally considered to come from the corona and chromosphere of the star, respectively. The observed QPP have periods $1.5 \pm 0.15$ hours (SXR) and $3 \pm 0.6$ hours (WL), and correlate well with each other. The unique relationship between the observed parameters of QPP in SXR and WL allowed us to link them with oscillations of the electric current in the flare loop, which directly affect the dynamics of non-thermal electrons and indirectly (via Ohmic heating) the thermal plasma. These findings could be considered in favour of the equivalent LCR-contour model of a flare loop, at least in the extreme conditions of a stellar superflare.


title: The Minimal Astration Hypothesis -- a Necessity for Solving the Dust Budget Crisis?

abstract: Assuming that gas and dust separate in the interstellar medium (ISM) so that high-density regions, where stars can form, are almost devoid of dust, the amount of metals being removed from the ISM can be significantly reduced (minimized astration). Here, it is shown by simple analytical models that this may increase the total metal budget of a galaxy considerably. It is suggested that these extra metals may increase the mass of dust such that the "dust budget crisis", i.e., the fact that there seems to be more dust at high redshifts than can be accounted for, can be ameliorated. Reducing the amount of astration, the metal budget can be more than doubled, in particular for systems that evolve under continuous gas accretion.


title: Kerr-Schild-Kundt Metrics in Generic Einstein-Maxwell Theories

abstract: We study Kerr-Schild-Kundt class of metrics in generic gravity theories with Maxwell's field. We prove that these metrics linearize and simplify the field equations of generic gravity theories with Maxwell's field.


title: The Importance of Discussing Assumptions when Teaching Bootstrapping

abstract: Bootstrapping and other resampling methods are increasingly appearing in the textbooks and curricula of courses that introduce undergraduate students to statistical methods. In order to teach the bootstrap well, students and instructors need to be aware of the assumptions behind these intervals. In this article we discuss important assumptions about simple non-parametric bootstrap intervals and their corresponding hypothesis tests. We present simulations that instructors can use to help students understand some of the assumptions behind these methods. The simulations will be especially relevant to instructors who desire to increase accessibility for students from non-mathematical backgrounds, including those with math anxiety.


title: Broadband, efficient extraction of quantum light by a photonic device comprised of a metallic nano-ring and a gold back reflector

abstract: To implement quantum light sources based on quantum emitters in applications, it is desirable to improve the extraction efficiency of single photons. In particular controlling the directionality and solid angle of the emission are key parameters, for instance, to couple single photons into optical fibers and send the information encoded in quantum light over long distances, for quantum communication applications. In addition, fundamental studies of the radiative behavior of quantum emitters, including studies of coherence and blinking, benefit from such improved photon collection. Quantum dots grown via Stranski-Krastanov technique have shown to be good candidates for bright, coherent, indistinguishable quantum light emission. However, one of the challenges associated with these quantum light sources arises from the fact that the emission wavelengths can vary from one emitter to the other. To this end, broadband light extractors that do not rely on high-quality factor optical cavities would be desirable, so that no tuning between the quantum dot emission wavelength and the resonator used to increase the light extraction is needed. Here, we show that metallic nano-rings combined with gold back reflectors increase the collection efficiency of single photons and we study the statistics of this effect when quantum dots are spatially randomly distributed within the nano-rings. We show an average increase in the brightness of about a factor 7.5, when comparing emitters within and outside the nano-rings in devices with a gold back reflector, we measure count rates exceeding 7 x 10^6 photons per second and single photon purities as high as 85% +/- 1%. These results are important steps towards the realisation of scalable, broadband, easy to fabricate sources of quantum light for quantum communication applications.

categories: physics.optics cond-mat.mes-hall physics.app-ph quant-ph


title: Trees with power-like height dependent weight

abstract: We consider planar rooted random trees whose distribution is even for fixed height $h$ and size $N$ and whose height dependence is given by a power function $h^\alpha$. Defining the total weight for such trees of fixed size to be $Z_N$, a detailed analysis of the analyticity properties of the corresponding generating function is provided. Based on this, we determine the asymptotic form of $Z_N$ and show that the local limit at large size is identical to the Uniform Infinite Planar Tree, independent of the exponent $\alpha$ of the height distribution function.


title: A Mathematical Model for the Origin of Name Brands and Generics

abstract: Firms in the U.S. spend over 200 billion dollars each year advertising their products to consumers, around one percent of the country's gross domestic product. It is of great interest to understand how that aggregate expenditure affects prices, market efficiency, and overall welfare. Here, we present a mathematical model for the dynamics of competition through advertising and find a surprising prediction: when advertising is relatively cheap compared to the maximum benefit advertising offers, rational firms split into two groups, one with significantly less advertising (a "generic" group) and one with significantly more advertising (a "name brand" group). Our model predicts that this segmentation will also be reflected in price distributions; we use large consumer data sets to test this prediction and find good qualitative agreement.


title: Constructive nonlocal games with very small classical values

abstract: There are few explicit examples of two player nonlocal games with a large gap between classical and quantum value. One of the reasons is that estimating the classical value is usually a hard computational task. This paper is devoted to analyzing classical values of the so-called linear games (generalization of XOR games to a larger number of outputs). We employ nontrivial results from graph theory and combine them with number theoretic results used earlier in the context of harmonic analysis to obtain a novel tool -- {\it the girth method} -- allowing to provide explicit examples of linear games with prescribed low classical value. In particular, we provide games with minimal possible classical value. We then speculate on the potential unbounded violation, by comparing the obtained classical values with a known upper bound for the quantum value. If this bound can be even asymptotically saturated, our games would have the best ratio of quantum to classical value as a function of the product of the number of inputs and outputs when compared to other explicit (i.e. non-random) constructions.


title: Classifying Emails into Human vs Machine Category

abstract: It is an essential product requirement of Yahoo Mail to distinguish between personal and machine-generated emails. The old production classifier in Yahoo Mail was based on a simple logistic regression model. That model was trained by aggregating features at the SMTP address level. We propose building deep learning models at the message level. We built and trained four individual CNN models: (1) a content model with subject and content as input; (2) a sender model with sender email address and name as input; (3) an action model by analyzing email recipients' action patterns and correspondingly generating target labels based on senders' opening/deleting behaviors; (4) a salutation model by utilizing senders' "explicit salutation" signal as positive labels. Next, we built a final full model after exploring different combinations of the above four models. Experimental results on editorial data show that our full model improves the adjusted-recall from 70.5% to 78.8% compared to the old production model, while at the same time lifts the precision from 94.7% to 96.0%. Our full model also significantly beats the state-of-the-art Bert model at this task. This full model has been deployed into the current production system (Yahoo Mail 6).


title: Neighborhood Random Walk Graph Sampling for Regularized Bayesian Graph Convolutional Neural Networks

abstract: In the modern age of social media and networks, graph representations of real-world phenomena have become an incredibly useful source to mine insights. Often, we are interested in understanding how entities in a graph are interconnected. The Graph Neural Network (GNN) has proven to be a very useful tool in a variety of graph learning tasks including node classification, link prediction, and edge classification. However, in most of these tasks, the graph data we are working with may be noisy and may contain spurious edges. That is, there is a lot of uncertainty associated with the underlying graph structure. Recent approaches to modeling uncertainty have been to use a Bayesian framework and view the graph as a random variable with probabilities associated with model parameters. Introducing the Bayesian paradigm to graph-based models, specifically for semi-supervised node classification, has been shown to yield higher classification accuracies. However, the method of graph inference proposed in recent work does not take into account the structure of the graph. In this paper, we propose a novel algorithm called Bayesian Graph Convolutional Network using Neighborhood Random Walk Sampling (BGCN-NRWS), which uses a Markov Chain Monte Carlo (MCMC) based graph sampling algorithm utilizing graph structure, reduces overfitting by using a variational inference layer, and yields consistently competitive classification results compared to the state-of-the-art in semi-supervised node classification.


title: Toward an automated-algebra framework for high orders in the virial expansion of quantum matter

abstract: The virial expansion provides a non-perturbative view into the thermodynamics of quantum many-body systems in dilute regimes. While powerful, the expansion is challenging as calculating its coefficients at each order $n$ requires analyzing (if not solving) the quantum $n$-body problem. In this work, we present a comprehensive review of automated algebra methods, which we developed to calculate high-order virial coefficients. The methods are computational but non-stochastic, thus avoiding statistical effects; they are also for the most part analytic, not numerical, and amenable to massively parallel computer architectures. We show formalism and results for coefficients characterizing the thermodynamics (pressure, density, energy, static susceptibilities) of homogeneous and harmonically trapped systems, and explain how to generalize them to other observables such as the momentum distribution, Tan's contact, and the structure factor.


title: Learning to track environment state via predictive autoencoding

abstract: This work introduces a neural architecture for learning forward models of stochastic environments. The task is achieved solely through learning from temporal unstructured observations in the form of images. Once trained, the model allows for tracking of the environment state in the presence of noise or with new percepts arriving intermittently. Additionally, the state estimate can be propagated in observation-blind mode, thus allowing for long-term predictions. The network can output both expectation over future observations and samples from belief distribution. The resulting functionalities are similar to those of a Particle Filter (PF). The architecture is evaluated in an environment where we simulate objects moving. As the forward and sensor models are available, we implement a PF to gauge the quality of the models learnt from the data.


title: CEM-GD: Cross-Entropy Method with Gradient Descent Planner for Model-Based Reinforcement Learning

abstract: Current state-of-the-art model-based reinforcement learning algorithms use trajectory sampling methods, such as the Cross-Entropy Method (CEM), for planning in continuous control settings. These zeroth-order optimizers require sampling a large number of trajectory rollouts to select an optimal action, which scales poorly for large prediction horizons or high dimensional action spaces. First-order methods that use the gradients of the rewards with respect to the actions as an update can mitigate this issue, but suffer from local optima due to the non-convex optimization landscape. To overcome these issues and achieve the best of both worlds, we propose a novel planner, Cross-Entropy Method with Gradient Descent (CEM-GD), that combines first-order methods with CEM. At the beginning of execution, CEM-GD uses CEM to sample a significant amount of trajectory rollouts to explore the optimization landscape and avoid poor local minima. It then uses the top trajectories as initialization for gradient descent and applies gradient updates to each of these trajectories to find the optimal action sequence. At each subsequent time step, however, CEM-GD samples much fewer trajectories from CEM before applying gradient updates. We show that as the dimensionality of the planning problem increases, CEM-GD maintains desirable performance with a constant small number of samples by using the gradient information, while avoiding local optima using initially well-sampled trajectories. Furthermore, CEM-GD achieves better performance than CEM on a variety of continuous control benchmarks in MuJoCo with 100x fewer samples per time step, resulting in around 25% less computation time and 10% less memory usage. The implementation of CEM-GD is available at $\href{https://github.com/KevinHuang8/CEM-GD}{\text{https://github.com/KevinHuang8/CEM-GD}}$.


title: Jet Charge: A new tool to probe the anomalous $Zb\bar{b}$ couplings at the EIC

abstract: We propose to utilize the average jet charge weighted single-spin asymmetry $A_e^{bQ}$ in the neutral current DIS scattering with one $b$-tagged jet produced at the EIC to probe the $Zb\bar{b}$ interactions. This novel observable is sensitive to the axial-vector component of the $Zb\bar{b}$ coupling and is complementary to the single-spin asymmetry measurement at the EIC and the $gg\to Zh$ cross section measurement at the HL-LHC in the determination of the $Zb\bar{b}$ couplings. We show that the apparent degeneracy of the $Zb\bar{b}$ couplings, implied by the LEP and SLC precision electroweak data, can be broken by the $A_e^{bQ}$ measurement at the EIC. With a large enough integrated luminosity collected at the EIC, the measurement of $A_e^{bQ}$ could also resolve the long-standing discrepancy of $A_{\rm FB}^b$ at the LEP, strongly dependent on the $b$-tagging efficiency.


title: Flow of (higher) Berry curvature and bulk-boundary correspondence in parametrized quantum systems

abstract: This paper is concerned with the physics of parametrized gapped quantum many-body systems, which can be viewed as a generalization of conventional topological phases of matter. In such systems, rather than considering a single Hamiltonian, one considers a family of Hamiltonians that depend continuously on some parameters. After discussing the notion of phases of parametrized systems, we formulate a bulk-boundary correspondence for an important bulk quantity, the Kapustin-Spodyneiko higher Berry curvature, first in one spatial dimension and then in arbitrary dimension. This clarifies the physical interpretation of the higher Berry curvature, which in one spatial dimension is a flow of (ordinary) Berry curvature. In d dimensions, the higher Berry curvature is a flow of (d-1)-dimensional higher Berry curvature. Based on this, we discuss one-dimensional systems that pump Chern number to/from spatial boundaries, resulting in anomalous boundary modes featuring isolated Weyl points. In higher dimensions, there are pumps of the analogous quantized invariants obtained by integrating the higher Berry curvature. We also discuss the consequences for parametrized systems of Kitaev's proposal that invertible phases are classified by a generalized cohomology theory, and emphasize the role of the suspension isomorphism in generating new examples of parametrized systems from known invertible phases. Finally, we present a pair of general quantum pumping constructions, based on physical pictures introduced by Kitaev, which take as input a d-dimensional parametrized system, and produce new (d+1)-dimensional parametrized systems. These constructions are useful for generating examples, and we conjecture that one of the constructions realizes the suspension isomorphism in a generalized cohomology theory of invertible phases.


title: Entropy stable discontinuous Galerkin methods for the shallow water equations with subcell positivity preservation

abstract: High order schemes are known to be unstable in the presence of shock discontinuities or under-resolved solution features, and have traditionally required additional filtering, limiting, or artificial viscosity to avoid solution blow up. Entropy stable schemes address this instability by ensuring that physically relevant solutions satisfy a semi-discrete entropy inequality independently of discretization parameters. However, additional measures must be taken to ensure that solutions satisfy physical constraints such as positivity. In this work, we present a high order entropy stable discontinuous Galerkin (ESDG) method for the nonlinear shallow water equations (SWE) on two-dimensional (2D) triangular meshes which preserves the positivity of the water heights. The scheme combines a low order positivity preserving method with a high order entropy stable method using convex limiting. This method is entropy stable and well-balanced for fitted meshes with continuous bathymetry profiles.


title: Incomplete RG: Hawking-Page transition, C-theorem and relevant scalar deformations of global AdS

abstract: We discuss relevant scalar deformations of a holographic theory with a compact boundary. An example of such a theory would be the global AdS$_4$ with its spatially compact boundary $S^2$. To introduce a relevant deformation, we choose to turn on a time-independent and spatially homogeneous non-normalizable scalar operator with $m^2 = -2$. The finite size of a compact boundary cuts down the RG flow at a finite length scale leading to an incomplete RG flow to IR. We discuss a version of {\it incomplete} C-theorem and an {\it incomplete} attractor like mechanism. We discuss the implication of our results for entanglement entropy and geometric quantities like scalar curvature, volume and mass scale of fundamental excitation of the how these quantities increase or decrease (often monotonically) with the strength of the deformation. Thermal physics of a holographic theory defined on a compact boundary is more interesting than its non-compact counterpart. It is well known that with a compact boundary, there is a possibility of a first order Hawking-Page transition dual to a de-confinement phase transition. From a gravity perspective, a relevant deformation dumps negative energy inside the bulk, increasing the effective cosmological constant ($\Lambda$) of the AdS. Dumping more negative energy in the bulk would make the HP transition harder and the corresponding HP transition temperature would increase. However, we have found the size of the BH at the transition temperature decreases.


title: Learn bifurcations of nonlinear parametric systems via equation-driven neural networks

abstract: Nonlinear parametric systems have been widely used in modeling nonlinear dynamics in science and engineering. Bifurcation analysis of these nonlinear systems on the parameter space are usually used to study the solution structure such as the number of solutions and the stability. In this paper, we develop a new machine learning approach to compute the bifurcations via so-called equation-driven neural networks (EDNNs). The EDNNs consist of a two-step optimization: the first step is to approximate the solution function of the parameter by training empirical solution data; the second step is to compute bifurcations by using the approximated neural network obtained in the first step. Both theoretical convergence analysis and numerical implementation on several examples have been performed to demonstrate the feasibility of the proposed method.


title: Representation and Invariance in Reinforcement Learning

abstract: Researchers have formalized reinforcement learning (RL) in different ways. If an agent in one RL framework is to run within another RL framework's environments, the agent must first be converted, or mapped, into that other framework. Whether or not this is possible depends on not only the RL frameworks in question and but also how intelligence itself is measured. In this paper, we lay foundations for studying relative-intelligence-preserving mappability between RL frameworks. We define two types of mappings, called weak and strong translations, between RL frameworks and prove that existence of these mappings enables two types of intelligence comparison according to the mappings preserving relative intelligence. We investigate the existence or lack thereof of these mappings between: (i) RL frameworks where agents go first and RL frameworks where environments go first; and (ii) twelve different RL frameworks differing in terms of whether or not agents or environments are required to be deterministic. In the former case, we consider various natural mappings between agent-first and environment-first RL and vice versa; we show some positive results (some such mappings are strong or weak translations) and some negative results (some such mappings are not). In the latter case, we completely characterize which of the twelve RL-framework pairs admit weak translations, under the assumption of integer-valued rewards and some additional mild assumptions.


title: Optically defined cavities in driven-dissipative photonic lattices

abstract: The engineering of localised modes in photonic structures is one of the main targets of modern photonics. An efficient strategy to design these modes is to use the interplay of constructive and destructive interference in periodic photonic lattices. This mechanism is at the origin of defect modes in photonic bandgaps, bound states in the continuum and compact localised states in flat bands. Here we show that in lattices of lossy resonators, the addition of external optical drives with controlled phase enlarges the possibilities of manipulating interference effects and allows designing novel types of localised modes. Using a honeycomb lattice of coupled micropillars resonantly driven with several laser spots at energies within its photonic bands we demonstrate the localisation of light in at-will geometries down to a single site. These localised modes can be seen as fully reconfigurable optical cavities with the potentiality of enhancing nonlinear effects and of controlling light-matter interactions with single site resolution.


title: Squeezing Quantum Many-Body Scars

abstract: We develop an analytical approach for the description of quantum many-body scars in PXP models. We show that the scarred dynamics in the PXP model on a complete bipartite graph can be interpreted as a one-dimensional chiral scattering problem, and solve this problem analytically. The insights from this analysis allow us to predict that dynamical signatures of scars in PXP models can be enhanced by spin squeezing the initial states. We show numerically that this stabilization mechanism applies not only to the complete bipartite graph but also to one- and two-dimensional lattices, which are relevant for Rydberg atom array experiments. Moreover, our findings provide a physical motivation for Hamiltonian deformations reminiscent of those known to produce perfect scars.


title: Separate Exchangeability as Modeling Principle in Bayesian Nonparametrics

abstract: We argue for the use of separate exchangeability as a modeling principle in Bayesian nonparametric (BNP) inference. Separate exchangeability is \emph{de facto} widely applied in the Bayesian parametric case, e.g., it naturally arises in simple mixed models. However, while in some areas, such as random graphs, separate and (closely related) joint exchangeability are widely used, it is curiously underused for several other applications in BNP. We briefly review the definition of separate exchangeability focusing on the implications of such a definition in Bayesian modeling. We then discuss two tractable classes of models that implement separate exchangeability that are the natural counterparts of familiar partially exchangeable BNP models. The first is nested random partitions for a data matrix, defining a partition of columns and nested partitions of rows, nested within column clusters. Many recent models for nested partitions implement partially exchangeable models related to variations of the well-known nested Dirichlet process. We argue that inference under such models in some cases ignores important features of the experimental setup. We obtain the separately exchangeable counterpart of such partially exchangeable partition structures. The second class is about setting up separately exchangeable priors for a nonparametric regression model when multiple sets of experimental units are involved. We highlight how a Dirichlet process mixture of linear models known as ANOVA DDP can naturally implement separate exchangeability in such regression problems. Finally, we illustrate how to perform inference under such models in two real data examples.


title: Quantitatively improved finite-size criteria for spectral gaps

abstract: Finite-size criteria have emerged as an effective tool for deriving spectral gaps in higher-dimensional frustration-free quantum spin systems. We quantitatively improve the existing finite-size criteria by introducing a novel subsystem weighting scheme. The approach applies to Euclidean lattices of any dimension, the honeycomb lattice, and the triangular lattice.

categories: math-ph cond-mat.stat-mech math.MP quant-ph


title: Herding of proteins by the ends of shrinking polymers

abstract: The control of biopolymer length is mediated by proteins that localize to polymer ends and regulate polymerization dynamics. Several mechanisms have been proposed to achieve end localization. Here, we propose a novel mechanism by which a protein that binds to a shrinking polymer and slows its shrinkage will be spontaneously enriched at the shrinking end through a "herding" effect. We formalize this process using both a lattice-gas model and a continuum description, and we present experimental evidence that the microtubule regulator spastin employs this mechanism. Our findings extend to more general problems involving diffusion within shrinking domains.


title: Observations with the Southern Connecticut Stellar Interferometer. I. Instrument Description and First Results

abstract: We discuss the design, construction, and operation of a new intensity interferometer, based on the campus of Southern Connecticut State University in New Haven, Connecticut. While this paper will focus on observations taken with an original two-telescope configuration, the current instrumentation consists of three portable 0.6-m Dobsonian telescopes with single-photon avalanche diode (SPAD) detectors located at the Newtonian focus of each telescope. Photons detected at each station are time-stamped and read out with timing correlators that can give cross-correlations in timing to a precision of 48 ps. We detail our observations to date with the system, which has now been successfully used at our university in 16 nights of observing. Components of the instrument were also deployed on one occasion at Lowell Observatory, where the Perkins and Hall telescopes were made to function as an intensity interferometer. We characterize the performance of the instrument in detail. In total, the observations indicate the detection of a correlation peak at the level of 6.76-sigma when observing unresolved stars, and consistency with partial or no detection when observing at a baseline sufficient to resolve the star. Using these measurements we conclude that the angular diameter of Arcturus is larger than 15 mas, and that of Vega is between 0.8 and 17 mas. While the uncertainties are large at this point, both results are consistent with measures from amplitude-based long baseline optical interferometers.


title: Global energetics of solar flares. XIII. The Neupert effect and acceleration of coronal mass ejections

abstract: Our major aim is a height-time model $r(t)$ of the propagation of {\sl Coronal Mass Ejections (CMEs)}, where the lower corona is self-consistently connected to the heliospheric path. We accomplish this task by using the Neupert effect to derive the peak time, duration, and rate of the CME acceleration phase, as obtained from the time derivative of the {\sl soft X-ray (SXR)} light curve. This novel approach offers the advantage to obtain the kinematics of the CME height-time profile $r(t)$, the CME velocity profile $v(t)=dr(t)/dt$, and the CME acceleration profile $a(t)=dv(t)/dt$ from {\sl Geostationary Orbiting Earth Satellite (GOES)} and white-light data, without the need of {\sl hard X-ray (HXR)} data. We apply this technique to a data set of 576 (GOES X and M-class) flare events observed with GOES and the {\sl Large Angle Solar Coronagraph (LASCO)}. Our analysis yields acceleration rates in the range of $a_A = 0.1-13$ km s$^{-2}$, acceleration durations of $\tau_A = 1.2-45$ min, and acceleration distances in the range of $d_A = 3-1063$ Mm, with a median of $d_A=39$ Mm, which corresponds to the hydrostatic scale height of a corona with a temperature of $T_e \approx 0.8$ MK. The results are consistent with standard flare/CME models that predict magnetic reconnection and synchronized (primary) acceleration of CMEs in the low corona (at a height of ~0.1 R_sun), while secondary (weaker) acceleration may occur further out at heliospheric distances.


title: The impact of aging on human brain network target controllability

abstract: Understanding how few distributed areas can steer large-scale brain activity is a fundamental question that has practical implications, which range from inducing specific patterns of behavior to counteracting disease. Recent endeavors based on network controllability provided fresh insights into the potential ability of single regions to influence whole brain dynamics through the underlying structural connectome. However, controlling the entire brain activity is often unfeasible and might not always be necessary. The question whether single areas can control specific target subsystems remains crucial, albeit still poorly explored. Furthermore, the structure of the brain network exhibits progressive changes across the lifespan, but little is known about the possible consequences in the controllability properties. To address these questions, we adopted a novel target controllability approach that quantifies the centrality of brain nodes in controlling specific target anatomo-functional systems. We then studied such target control centrality in human connectomes obtained from healthy individuals aged from 5 to 85. Main results showed that the sensorimotor system has a high influencing capacity, but it is difficult for other areas to influence it. Furthermore, we reported that target control centrality varies with age and that temporal-parietal regions, whose cortical thinning is crucial in dementia-related diseases, exhibit lower values in older people. By simulating targeted attacks, such as those 19 occurring in focal stroke, we showed that the ipsilesional hemisphere is the most affected one regardless of the damaged area. Notably, such degradation in target control centrality was more evident in younger people, thus supporting early-vulnerability hypotheses after stroke.


title: Split Moves for Monte-Carlo Tree Search

abstract: In many games, moves consist of several decisions made by the player. These decisions can be viewed as separate moves, which is already a common practice in multi-action games for efficiency reasons. Such division of a player move into a sequence of simpler / lower level moves is called \emph{splitting}. So far, split moves have been applied only in forementioned straightforward cases, and furthermore, there was almost no study revealing its impact on agents' playing strength. Taking the knowledge-free perspective, we aim to answer how to effectively use split moves within Monte-Carlo Tree Search (MCTS) and what is the practical impact of split design on agents' strength. This paper proposes a generalization of MCTS that works with arbitrarily split moves. We design several variations of the algorithm and try to measure the impact of split moves separately on efficiency, quality of MCTS, simulations, and action-based heuristics. The tests are carried out on a set of board games and performed using the Regular Boardgames General Game Playing formalism, where split strategies of different granularity can be automatically derived based on an abstract description of the game. The results give an overview of the behavior of agents using split design in different ways. We conclude that split design can be greatly beneficial for single- as well as multi-action games.


title: Proofs of Two Formulas of Vladeta Jovovic

abstract: In this paper, we first provide an analytic and a bijective proof of a formula stated by Vladeta Jovovic in the OEIS sequence A117989. We also provide a bijective proof of another interesting result stated by him on the same page concerning integer partitions with fixed differences between the largest and smallest parts.


title: Phase Transitions in High Purity Zr Under Dynamic Compression

abstract: We present results from ramp compression experiments on high-purity Zr that show the $\alpha \rightarrow \omega$, $\omega \rightarrow \beta$, as well as reverse $\beta \rightarrow \omega$ phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic $\alpha \rightarrow \omega$ transition occurs $\sim 9$ GPa above the equilibrium phase boundary, the $\omega \rightarrow \beta$ transition occurs within 0.9~GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium $\omega - \beta$ line at $P= 29.2$ GPa, and $T = 490$ K. The thermodynamic path in the interior of the sample lies $\sim 100$ K above the isentrope at the point of the $\omega \rightarrow \beta$ transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium $\alpha \rightarrow \omega$ transition. The inferred rate of the $\alpha \rightarrow \omega$ transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. The small fractional volume change $\Delta V/V \approx 0.1$ at the $\alpha \rightarrow \omega$ transition amplifies the effect of shear stress, and we estimate that for this case shear stress is equivalent to a pressure increase in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.


title: Pluto's atmosphere in plateau phase since 2015 from a stellar occultation at Devasthal

abstract: A stellar occultation by Pluto was observed on 6 June 2020 with the 1.3-m and 3.6-m telescopes located at Devasthal, Nainital, India, using imaging systems in the I and H bands, respectively. From this event, we derive a surface pressure for Pluto's atmosphere of $p_{\rm surf}= 12.23^{+0.65}_{-0.38} $~$\mu$bar. This shows that Pluto's atmosphere is in a plateau phase since mid-2015, a result which is in excellent agreement with the Pluto volatile transport model of Meza et al. (2019}. This value does not support the pressure decrease reported by independent teams, based on occultations observed in 2018 and 2019, see Young et al. (2021} and Arimatu et al. (2020), respectively.


title: Nonlinear Discrete-time System Identification without Persistence of Excitation: Finite-time Concurrent Learning Methods

abstract: This paper deals with the problem of finite-time learning for unknown discrete-time nonlinear systems' dynamics, without the requirement of the persistence of excitation. Two finite-time concurrent learning methods are presented to approximate the uncertainties of the discrete-time nonlinear systems in an online fashion by employing current data along with recorded experienced data satisfying an easy-to-check rank condition on the richness of the recorded data which is less restrictive in comparison with persistence of excitation condition. For the proposed finite-time concurrent learning methods, rigorous proofs guarantee the finite-time convergence of the estimated parameters to their optimal values based on the discrete-time Lyapunov analysis. Compared with the existing work in the literature, simulation results illustrate that the proposed methods can timely and precisely approximate the uncertainties.


title: The Effect of Parameter Variations on the Performance of the Josephson Travelling Wave Parametric Amplifiers

abstract: We have simulated the performance of the Josephson Travelling Wave Parametric Amplifier (JTWPA) based on the one-dimensional array of RF SQUIDs. Unlike the ideal model in which all SQUIDs are assumed to be identical, we allowed variation of the device parameters such as the geometric inductance of the SQUID loop, capacitance to ground, Josephson junction capacitance and critical current. Our simulations confirm the negative effects of variation of the device parameters leading to microwave reflections between individual cells and the shift of the flux bias from the optimal point. The strongest effect is caused by the variation of the geometric inductance as it varies both the wave impedance and the flux bias. The most detrimental, however, are point defects, such as shorts to ground making the circuit opaque to microwaves. This imposes stringent requirements on the fabrication process making it extremely challenging. We highlight the strict limitations on parameter spread in these devices while also discussing the robustness of the scheme to variation


title: Strong coupling in a Au plasmonic antenna-SiO$_{2}$ layer system: a hybrid mode analysis

abstract: A detailed analysis of the optical response of a system accommodating several coupled modes is needed for the complete understanding of the strong coupling effect. In this paper, we report on the analysis of scattering cross section spectra of Au antennas on a SiO$_{2}$ layer on a Si substrate in the IR region. A classical model of coupled oscillators is used for determining the resonant energies, damping rates and coupling strengths of four phonon polariton modes in the SiO$_{2}$ layer coupled to a localized surface plasmon mode in a Au antenna. The calculated Hopfield mixing coefficients then show the contribution of the individual uncoupled modes to the hybrid modes of the coupled system.