title: Finite-Size Effects On The Self-Force

abstract: Electromagnetic and linear gravitational radiation do not solely propagate on the null cone in 3+1 dimensions in curved spacetimes, contrary to their well-known behavior in flat spacetime. Their additional propagation inside the null cone is known as the tail effect. A compact body will produce a signal whose tail will interact with its future worldline, thus producing a tail-induced self-force. We present new results for the tail-induced scalar, electromagnetic and gravitational self-force for a test mass in orbit around a central mass, including effects from the internal structure of that body.


title: Edge theories of 2D fermionic symmetry protected topological phases protected by unitary Abelian symmetries

abstract: Abelian Chern-Simons theory, characterized by the so-called $K$ matrix, has been quite successful in characterizing and classifying Abelian fractional quantum hall effect (FQHE) as well as symmetry protected topological (SPT) phases, especially for bosonic SPT phases. However, there are still some puzzles in dealing with fermionic SPT(fSPT) phases. In this paper, we utilize the Abelian Chern-Simons theory to study the fSPT phases protected by arbitrary Abelian total symmetry $G_f$. Comparing to the bosonic SPT phases, fSPT phases with Abelian total symmetry $G_f$ has three new features: (1) it may support gapless majorana fermion edge modes, (2) some nontrivial bosonic SPT phases may be trivialized if $G_f$ is a nontrivial extention of bosonic symmetry $G_b$ over $\mathbb{Z}_2^f$, (3) certain intrinsic fSPT phases can only be realized in interacting fermionic system. We obtain edge theories for various fSPT phases, which can also be regarded as conformal field theories (CFT) with proper symmetry anomaly. In particular, we discover the construction of Luttinger liquid edge theories with central charge $n-1$ for Type-III bosonic SPT phases protected by $(\mathbb{Z}_n)^3$ symmetry and the Luttinger liquid edge theories for intrinsically interacting fSPT protected by unitary Abelian symmetry. The ideas and methods used here might be generalized to derive the edge theories of fSPT phases with arbitrary unitary finite Abelian total symmetry $G_f$.


title: Cubic Stylization

abstract: We present a 3D stylization algorithm that can turn an input shape into the style of a cube while maintaining the content of the original shape. The key insight is that cubic style sculptures can be captured by the as-rigid-as-possible energy with an l1-regularization on rotated surface normals. Minimizing this energy naturally leads to a detail-preserving, cubic geometry. Our optimization can be solved efficiently without any mesh surgery. Our method serves as a non-realistic modeling tool where one can incorporate many artistic controls to create stylized geometries.


title: Optical and photoemission investigation of structural and magnetic transitions in the iron-based superconductor Sr$_\mathbf{0.67}$Na$_\mathbf{0.33}$Fe$_\mathbf{2}$As$_\mathbf{2}$

abstract: We report the temperature-dependent optical conductivity and ARPES studies of the iron-based superconductor (SC) Sr$_{0.67}$Na$_{0.33}$Fe$_2$As$_2$ in the high-temperature tetragonal paramagnetic phase; below the structural and magnetic transitions at $T_{\rm N}\simeq$125 K in the orthorhombic spin-density-wave (SDW)-like phase, and $T_r\simeq$42 K in the reentrant tetragonal double-Q magnetic phase where both charge and SDW order exist; and below the SC transition at $T_c\simeq$10 K. The free-carrier component in the optical conductivity is described by two Drude contributions; one strong and broad, the other weak and narrow. The broad Drude component decreases dramatically below $T_{\rm N}$ and $T_r$, with much of its strength being transferred to a bound excitation in the mid-infrared, while the narrow Drude component shows no anomalies at either of the transitions, actually increasing in strength at low temperature while narrowing dramatically. The behavior of an infrared-active mode suggests zone-folding below $T_r$. Below $T_c$ the dramatic decrease in the low-frequency optical conductivity signals the formation of a SC energy gap. ARPES reveals hole-like bands at the center of the Brillouin zone (BZ), with both electron- and hole-like bands at the corners. Below $T_{\rm N}$, the hole pockets at the center of the BZ decrease in size, consistent with the behavior of the broad Drude component; while below $T_r$ the electron-like bands shift and split, giving rise to a low-energy excitation in the optical conductivity at ~20 meV. The magnetic states, with resulting SDW and charge-SDW order, respectively, lead to a significant reconstruction of the Fermi surface that has profound implications for the transport originating from the electron and hole pockets, but appears to have relatively little impact on the SC in this material.


title: MOBSTER -- III. HD 62658: a magnetic Bp star in an eclipsing binary with a non-magnetic 'identical twin'

abstract: HD 62658 (B9p V) is a little-studied chemically peculiar star. Light curves obtained by the Kilodegree Extremely Little Telescope (KELT) and Transiting Exoplanet Survey Satellite (TESS) show clear eclipses with a period of about 4.75 d, as well as out-of-eclipse brightness modulation with the same 4.75 d period, consistent with synchronized rotational modulation of surface chemical spots. High-resolution ESPaDOnS circular spectropolarimetry shows a clear Zeeman signature in the line profile of the primary; there is no indication of a magnetic field in the secondary. PHOEBE modelling of the light curve and radial velocities indicates that the two components have almost identical masses of about 3 M$_\odot$. The primary's longitudinal magnetic field $\langle B_z \rangle$ varies between about $+100$ and $-250$ G, suggesting a surface magnetic dipole strength $B_{\rm d} = 850$~G. Bayesian analysis of the Stokes $V$ profiles indicates $B_{\rm d} = 650$~G for the primary and $B_{\rm d} < 110$ G for the secondary. The primary's line profiles are highly variable, consistent with the hypothesis that the out-of-eclipse brightness modulation is a consequence of rotational modulation of that star's chemical spots. We also detect a residual signal in the light curve after removal of the orbital and rotational modulations, which might be pulsational in origin; this could be consistent with the weak line profile variability of the secondary. This system represents an excellent opportunity to examine the consequences of magnetic fields for stellar structure via comparison of two stars that are essentially identical with the exception that one is magnetic. The existence of such a system furthermore suggests that purely environmental explanations for the origin of fossil magnetic fields are incomplete.


title: Thermal conductivity of graphene polymorphs and compounds: from C3N to graphdiyne lattices

abstract: Tremendous experimental and theoretical attempts to find carbon based two-dimensional semiconductors have yielded a wide variety of graphene polymorphs, such as carbon-nitride, carbonboride, graphyne and graphdiyne 2D materials with highly attractive physical and chemical properties. In this study, by conducting extensive non-equilibrium molecular dynamics simulations, we have calculated and compared the thermal conductivity of thirteen prominent carbon-based structures at different lengths and two main chirality directions. Acquired results show that the structures of C3N, C3B and C2N exhibit the highest thermal conductivity, respectively, which suggest them as suitable candidates for thermal management systems in order to enhance the heat dissipation rates. In contrast, generally graphdiyne lattices and in particular 18-6-Gdy graphdiyne yields the lowest thermal conductivity, which can be a promising feature for thermoelectric applications. As a remarkable finding, we could establish connections between the thermal conductivity and density or Young's modulus of carbon based 2D systems, which can be employed to estimate the thermal conductivity of other polymorphs. Those results can provide a comprehensive viewpoint on the thermal transport properties of the nonporous and exceedingly porous carbon based 2D materials and may be used as useful guides for future designs in thermal management.


title: A Case Study on Combining ASR and Visual Features for Generating Instructional Video Captions

abstract: Instructional videos get high-traffic on video sharing platforms, and prior work suggests that providing time-stamped, subtask annotations (e.g., "heat the oil in the pan") improves user experiences. However, current automatic annotation methods based on visual features alone perform only slightly better than constant prediction. Taking cues from prior work, we show that we can improve performance significantly by considering automatic speech recognition (ASR) tokens as input. Furthermore, jointly modeling ASR tokens and visual features results in higher performance compared to training individually on either modality. We find that unstated background information is better explained by visual features, whereas fine-grained distinctions (e.g., "add oil" vs. "add olive oil") are disambiguated more easily via ASR tokens.


title: On the shape of the $\bar d-\bar u$ asymmetry

abstract: Using data from a recent reanalysis of neutron structure functions extracted from inclusive proton and deuteron deep-inelastic scattering (DIS), we re-examine the constraints on the shape of the $\bar d-\bar u$ asymmetry in the proton at large parton momentum fractions $x$. A global analysis of the proton-neutron structure function difference from BCDMS, NMC, SLAC and Jefferson Lab DIS measurements, and of Fermilab Drell-Yan lepton-pair production cross sections, suggests that existing data can be well described with $\bar d > \bar u$ for all values of $x$ currently accessible. We compare the shape of the fitted $\bar d-\bar u$ distributions with expectations from nonperturbative models based on chiral symmetry breaking, which can be tested by upcoming Drell-Yan data from the SeaQuest experiment at larger values of $x$.


title: Multi-Modal Machine Learning for Flood Detection in News, Social Media and Satellite Sequences

abstract: In this paper we present our methods for the MediaEval 2019 Mul-timedia Satellite Task, which is aiming to extract complementaryinformation associated with adverse events from Social Media andsatellites. For the first challenge, we propose a framework jointly uti-lizing colour, object and scene-level information to predict whetherthe topic of an article containing an image is a flood event or not.Visual features are combined using early and late fusion techniquesachieving an average F1-score of82.63,82.40,81.40and76.77. Forthe multi-modal flood level estimation, we rely on both visualand textual information achieving an average F1-score of58.48and46.03, respectively. Finally, for the flooding detection in time-based satellite image sequences we used a combination of classicalcomputer-vision and machine learning approaches achieving anaverage F1-score of58.82%


title: Gordian Adjacency for Positive Braid Knots

abstract: A knot $K_1$ is said to be Gordian adjacent to a knot $K_2$ if $K_1$ is an intermediate knot on an unknotting sequence of $K_2$. We extend previous results on Gordian adjacency by showing sufficient conditions for Gordian adjacency between classes of positive braid knots through manipulations of braid words. In addition, we explore unknotting sequences of positive braid knots and give a proof that there are only finitely many positive braid knots for a given unknotting number.


title: Algorithm-Dependent Generalization Bounds for Overparameterized Deep Residual Networks

abstract: The skip-connections used in residual networks have become a standard architecture choice in deep learning due to the increased training stability and generalization performance with this architecture, although there has been limited theoretical understanding for this improvement. In this work, we analyze overparameterized deep residual networks trained by gradient descent following random initialization, and demonstrate that (i) the class of networks learned by gradient descent constitutes a small subset of the entire neural network function class, and (ii) this subclass of networks is sufficiently large to guarantee small training error. By showing (i) we are able to demonstrate that deep residual networks trained with gradient descent have a small generalization gap between training and test error, and together with (ii) this guarantees that the test error will be small. Our optimization and generalization guarantees require overparameterization that is only logarithmic in the depth of the network, while all known generalization bounds for deep non-residual networks have overparameterization requirements that are at least polynomial in the depth. This provides an explanation for why residual networks are preferable to non-residual ones.


title: Automated Enriched Medical Concept Generation for Chest X-ray Images

abstract: Decision support tools that rely on supervised learning require large amounts of expert annotations. Using past radiological reports obtained from hospital archiving systems has many advantages as training data above manual single-class labels: they are expert annotations available in large quantities, covering a population-representative variety of pathologies, and they provide additional context to pathology diagnoses, such as anatomical location and severity. Learning to auto-generate such reports from images present many challenges such as the difficulty in representing and generating long, unstructured textual information, accounting for spelling errors and repetition/redundancy, and the inconsistency across different annotators. We therefore propose to first learn visually-informative medical concepts from raw reports, and, using the concept predictions as image annotations, learn to auto-generate structured reports directly from images. We validate our approach on the OpenI [2] chest x-ray dataset, which consists of frontal and lateral views of chest x-ray images, their corresponding raw textual reports and manual medical subject heading (MeSH ) annotations made by radiologists.


title: eV-scale sterile neutrinos from an extra dimension

abstract: Motivated by the short-baseline neutrino oscillation anomalies that suggest the existence of sterile neutrinos at the eV scale, we construct a scenario of a seesaw mechanism for 3+1 light neutrinos implemented by warped compactification of an extra dimension. As the seesaw mechanism necessitates at least two right-handed neutrinos at mass scales much larger than eV, incorporating an eV-scale sterile neutrino into a seesaw entails large mass hierarchies among the singlet neutrinos. We show that such hierarchies can be naturally explained by moderate fluctuations of the five-dimensional fermion mass parameters.


title: Explicit and implicit error inhibiting schemes with post-processing

abstract: Efficient high order numerical methods for evolving the solution of an ordinary differential equation are widely used. The popular Runge--Kutta methods, linear multi-step methods, and more broadly general linear methods, all have a global error that is completely determined by analysis of the local truncation error. In prior work in we investigated the interplay between the local truncation error and the global error to construct {\em error inhibiting schemes} that control the accumulation of the local truncation error over time, resulting in a global error that is one order higher than expected from the local truncation error. In this work we extend our error inhibiting framework to include a broader class of time-discretization methods that allows an exact computation of the leading error term, which can then be post-processed to obtain a solution that is two orders higher than expected from truncation error analysis. We define sufficient conditions that result in a desired form of the error and describe the construction of the post-processor. A number of new explicit and implicit methods that have this property are given and tested on a variety of ordinary and partial differential equation. We show that these methods provide a solution that is two orders higher than expected from truncation error analysis alone.


title: Transfer operator approach to ray-tracing in circular domains

abstract: The computation of wave-energy distributions in the mid-to-high frequency regime can be reduced to ray-tracing calculations. Solving the ray-tracing problem in terms of an operator equation for the energy density leads to an inhomogeneous equation which involves a Perron-Frobenius operator defined on a suitable Sobolev space. Even for fairly simple geometries, let alone realistic scenarios such as typical boundary value problems in room acoustics or for mechanical vibrations, numerical approximations are necessary. Here we study the convergence of approximation schemes by rigorous methods. For circular billiards we prove that convergence of finite-rank approximations using a Fourier basis follows a power law where the power depends on the smoothness of the source distribution driving the system. The relevance of our studies for more general geometries is illustrated by numerical examples.


title: Arithmetic {B}reuil-{K}isin-{F}argues modules and comparison of integral {p}-adic Hodge theories

abstract: Let $K$ be a discrete valuation field with perfect residue field, we study the functor from weakly admissible filtered $(\varphi,N,G_K)$-modules over $K$ to the isogeny category of Breuil-Kisin-Fargues $G_K$-modules. This functor is the composition of a functor defined by Fargues-Fontaine from weakly admissible filtered $(\varphi,N,G_K)$-modules to $G_K$-equivariant modifications of vector bundles over the Fargues-Fontaine curve $X_{FF}$, with the functor of Fargues-Scholze that between the category of admissible modifications of vector bundles over $X_{FF}$ and the isogeny category of Breuil-Kisin-Fargues modules. We study those objects appear in the essential image of the above functor and call them arithmetic BKF modules. We show certain rigidity result of arithmetic BKF modules and use it to compare existing $p$-adic Hodge theories at ${A}_{\mathrm{inf}}$ level.


title: Deformable Kernels: Adapting Effective Receptive Fields for Object Deformation

abstract: Convolutional networks are not aware of an object's geometric variations, which leads to inefficient utilization of model and data capacity. To overcome this issue, recent works on deformation modeling seek to spatially reconfigure the data towards a common arrangement such that semantic recognition suffers less from deformation. This is typically done by augmenting static operators with learned free-form sampling grids in the image space, dynamically tuned to the data and task for adapting the receptive field. Yet adapting the receptive field does not quite reach the actual goal -- what really matters to the network is the "effective" receptive field (ERF), which reflects how much each pixel contributes. It is thus natural to design other approaches to adapt the ERF directly during runtime. In this work, we instantiate one possible solution as Deformable Kernels (DKs), a family of novel and generic convolutional operators for handling object deformations by directly adapting the ERF while leaving the receptive field untouched. At the heart of our method is the ability to resample the original kernel space towards recovering the deformation of objects. This approach is justified with theoretical insights that the ERF is strictly determined by data sampling locations and kernel values. We implement DKs as generic drop-in replacements of rigid kernels and conduct a series of empirical studies whose results conform with our theories. Over several tasks and standard base models, our approach compares favorably against prior works that adapt during runtime. In addition, further experiments suggest a working mechanism orthogonal and complementary to previous works.


title: The aeolian-erosion barrier for the growth of metre-size objects in protoplanetary-discs

abstract: Aeolian-erosion is a destructive process which can erode small-size planetary objects through their interaction with a gaseous environment. Aeolian-erosion operates in a wide range of environments and under various conditions. Aeolian-erosion has been extensively explored in the context of geophysics in terrestrial planets. Here we show that aeolian-erosion of cobbles, boulders and small planetesimals in protoplanetary-discs can constitute a significant barrier for the early stages of planet formation. We use analytic calculations to show that under the conditions prevailing in protoplanetary-discs small bodies ($10-10^4 \rm{m}$) are highly susceptible to gas-drag aeolian-erosion. At this size-range aeolian-erosion can efficiently erode the planetesimals down to tens-cm size and quench any further growth of such small bodies. It thereby raises potential difficulties for channels suggested to alleviate the metre-size barrier. Nevertheless, the population of $\sim$decimetre-size cobbles resulting from aeolian-erosion might boost the growth of larger (>km size) planetesimals and planetary embryos through increasing the efficiency of pebble-accretion, once/if such large planetesimals and planetary embryos exist in the disc.


title: Solving the non-relativistic electronic Schrodinger equation with manipulating the coupling strength parameter over the electron-electron Coulomb integrals

abstract: The non-relativistic electronic Hamiltonian, H(a)= Hkin + Hne + aHee, extended with coupling strength parameter (a), allows to switch the electron-electron repulsion energy off and on. First, the easier a=0 case is solved and the solution of real (physical) a=1 case is generated thereafter from it to calculate the total electronic energy (Etotal electr,K) mainly for ground state (K=0). This strategy is worked out with utilizing generalized Moller-Plesset (MP), square of Hamiltonian (H2) and Configuration interactions (CI) devices. Applying standard eigensolver for Hamiltonian matrices (one or two times) buys off the needs of self-consistent field (SCF) convergence in this algorithm, along with providing the correction for basis set error and correlation effect. (SCF convergence is typically performed in the standard HF-SCF/basis/a=1 routine in today practice.)


title: Strain relaxation during the layer by layer growth of cubic CdSe onto ZnSe

abstract: A detailed reflection high-energy electron diffraction analysis shows relevant features of the lattice parameter relaxation of CdSe thin films grown in a layer-by-layer mode onto ZnSe. In situ investigations of different azimuths show a clear lattice parameter oscillation in the 110 azimuth. The lattice parameter has a minimum value ~similar to that of ZnSe! during Se exposure steps, and a higher and increasing lattice parameter during Cd exposure steps. The behavior is ascribed to the formation of CdSe islands during Cd exposure steps. The cumulative effect in CdSe exposure steps is considered to be a consequence of a decrease in the island size with the number of cycles. Actual plastic deformation does occur after 5 ML.


title: Adsorption, desorption, and interdiffusion in atomic layer epitaxy of CdTe and CdZnTe

abstract: The mechanisms controlling the growth rate and composition of epitaxial CdTe and CdZnTe films were studied. The films were grown by isothermal closed space configuration technique. A GaAs 100 substrate was exposed sequentially to the elemental sources, Zn, Te, and Cd, in isothermal conditions. While growth of ZnTe followed an atomic layer epitaxy, ALE, regime self regulated at one monolayer per cycle; the CdTe films revealed different growth rates in dependence of the growth parameters,exposure and purge times. Combination of short purge times and larger Cd exposure times led to not self regulated growth regime for CdTe. This is ascribed to large Cd coverages that were dependent on Cd exposure times, following a Brunauer-Emmett and Teller-type adsorption. However, for longer purge times and or short Cd exposure times, an ALE self regulated regime was achieved with 2 ML per cycle. In this sense, the self-regulation of the growth is limited by desorption, instead of absorption, as in the traditional growth technique. Cd atoms substitution by Zn atoms and subsequent evaporation of surface Cd atoms during Zn exposure has been proved. The influence of these facts on the growth and composition of the alloy is discussed.


title: V-Gas: Generating High Gas Consumption Inputs to Avoid Out-of-Gas Vulnerability

abstract: The out-of-gas error occurs when smart contract programs are provided with inputs that cause excessive gas consumption, and would be easily exploited to make the DoS attack. Multiple approaches have been proposed to estimate the gas limit of a function in smart contracts to avoid such error. However, under estimation often happens when the contract is complicated. In this work, we propose V-Gas, which could automatically generate inputs that maximizes the gas cost and reduce the under estimation cases. V-Gas is designed based on feedback-directed mutational fuzz testing. First, V-Gas builds the gas weighted control flow graph (CFG) of functions in smart contracts. Then, V-Gas develops gas consumption guided selection and mutation strategies to generate the input that maximize the gas consumption. For evaluation, we implement V-Gas based on js-evm, a widely used ethereum virtual machine written in javascript, and conduct experiments on 736 real-world transactions recorded on Ethereum. 44.02\% of the transactions would have out-of-gas errors under the estimation results given by solc, means that the recorded real gas consumption for those recorded transactions is larger than the gas limit value estimated by solc. While V-Gas could reduce the under estimation ratio to 13.86\%. Furthermore, V-Gas has exposed 25 previously unknown out-of-gas vulnerabilities in those widely-used smart contracts, 5 of which have been assigned unique CVE identifiers in the US National Vulnerability Database.


title: Fractional integral equations tell us how to impose initial values in fractional differential equations

abstract: The goal of this work is to discuss how should we impose initial values in fractional problems to ensure that they have exactly one smooth unique solution, where smooth simply means that the solution lies in a certain suitable space of fractional differentiability. For the sake of simplicity and to show the fundamental ideas behind our arguments, we will do this only for the Riemann-Liouville case of linear equations with constant coefficients. In a few words, we study the natural consequences in fractional differential equations of the already existing results involving existence and uniqueness for their integral analogues, in terms of the Riemann-Liouville fractional integral. Under this scope, we derive naturally several interesting results. One of the most astonishing ones is that a fractional differential equation of order $\beta>0$ with Riemann-Liouville derivatives can demand, in principle, less initial values than $\lceil \beta \rceil$ to have a uniquely determined solution. In fact, if not all the involved derivatives have the same decimal part, the amount of conditions is given by $\lceil \beta - \beta_* \rceil$ where $\beta_*$ is the highest order in the differential equation such that $\beta-\beta_*$ is not an integer.


title: Hadron production in pp and p-Pb collisions: A mass dependent phenomenon

abstract: The mass dependence plays a significant role in the yield enhancement or suppression of hadrons in pp and p-Pb collisions at the LHC energies. This has been observed by parameterizing the variation of yield ratios between any two hadrons with event charged-particle multiplicity using a single empirical function. We notice that this variation is independent of all quantum numbers and solely depends on masses of hadrons and masses of their valence quarks. The function shows that the amount of quark deconfinement increases with event multiplicity, and the quark coalescence favours more the production of heavier hadrons compared to lighter ones.


title: The false myth of the rise in self-citations, and the impressively positive effect of bibliometric evaluations on the increase of the impact of Italian research

abstract: It has recently been claimed by Baccini and coauthors that due to ANVUR's bibliometric evaluations of individuals, departments, and universities, in Italy there has been a surge in self-citations in the last ten years, thus increasing the "inwardness" of Italian research more than has happened abroad. We have studied the database of Ioannidis et al. published on 12 August 2019 of the one hundred thousand most "highly cited" scientists, including about two thousand Italians, and we found that the problem of self-citations in relation to this scientific elite is not significant in Italy, while perhaps observing a small deviation in the low scores in the rankings. The effect indicated by Baccini et al. consequently, does not seem worrying for the scientific elite (we quantified it in 2 percent of the total of scientists of the "best" one hundred thousand), and is probably largely concentrated in the further less cited scientists. Evaluation agencies like ANVUR should probably exclude self-citations in future evaluations, for the noise introduced by the young researchers. The overall state of health of the Italian research system and the positive effect of the ANVUR assessments are demonstrated by the number of Italian researchers in the top one hundred thousand, which has increased by comparing the "career" databased of 22 years, with that of the "young" researchers in the "2017" database. Italy, looking at the elite researchers, not only is not the most indulgent in self-citations, but has shown the best improvements, proving that the introduction of ANVUR had a positive effect. Indeed, all countries apart from Italy have suffered a decline, even substantial (-20 percent on a national Japan scale), of the number of researchers present in the 2017 data sets compared to career data. Italy instead shows a +0.2 percent on a global basis and an impressive +11.53 percent on a national basis.


title: Simulation of the multiphase configuration and phase transitions with quantum walks utilizing a step-dependent coin

abstract: Quantum walks are versatile simulators of topological phases and phase transitions as observed in condensed matter physics. Here, we utilize a step dependent coin in quantum walks and investigate what topological phases we can simulate with it, their topological invariants, bound states and possibility of phase transitions. These quantum walks simulate non-trivial phases characterized by topological invariants (winding number) $\pm 1$ which are similar to the ones observed in topological insulators and polyacetylene. We confirm that the number of phases and their corresponding bound states increase step dependently. In contrast, the size of topological phase and distance between two bound states are decreasing functions of steps resulting into formation of multiple phases as quantum walks proceed (multiphase configuration). We show that, in the bound states, the winding number and group velocity are ill-defined, and the second moment of the probability density distribution in position space undergoes an abrupt change. Therefore, there are phase transitions taking place over the bound states and between two topological phases with different winding numbers.


title: Enumeration of Sets of Mutually Orthogonal Latin Rectangles

abstract: We study sets of mutually orthogonal Latin rectangles (MOLR), and a natural variation of the concept of self-orthogonal Latin squares which is applicable on larger sets of mutually orthogonal Latin squares and MOLR, namely that each Latin rectangle in a set of MOLR is isotopic to each other rectangle in the set. We call such a set of MOLR \emph{homogeneous}. In the course of doing this, we perform a complete enumeration of non-isotopic sets of $t$ mutually orthogonal $k\times n$ Latin rectangles for $k\leq n \leq 7$, for all $t < n$. Specifically, we keep track of homogeneous sets of MOLR, as well as sets of MOLR where the autotopism group acts transitively on the rectangles, and we call such sets of MOLR \emph{transitive}. We build the sets of MOLR row by row, and in this process we also keep track of which of the MOLR are homogeneous and/or transitive in each step of the construction process. We use the prefix \emph{stepwise} to refer to sets of MOLR with this property. Sets of MOLR are connected to other discrete objects, notably finite geometries and certain regular graphs. Here we observe that all projective planes of order at most 9 except the Hughes plane can be constructed from a stepwise transitive MOLR.


title: Joint analysis of clinical risk factors and 4D cardiac motion for survival prediction using a hybrid deep learning network

abstract: In this work, a novel approach is proposed for joint analysis of high dimensional time-resolved cardiac motion features obtained from segmented cardiac MRI and low dimensional clinical risk factors to improve survival prediction in heart failure. Different methods are evaluated to find the optimal way to insert conventional covariates into deep prediction networks. Correlation analysis between autoencoder latent codes and covariate features is used to examine how these predictors interact. We believe that similar approaches could also be used to introduce knowledge of genetic variants to such survival networks to improve outcome prediction by jointly analysing cardiac motion traits with inheritable risk factors.


title: The Magnet of the Scattering and Neutrino Detector for the SHiP experiment at CERN

abstract: The Search for Hidden Particles (SHiP) experiment proposal at CERN demands a dedicated dipole magnet for its scattering and neutrino detector. This requires a very large volume to be uniformly magnetized at B > 1.2 T, with constraints regarding the inner instrumented volume as well as the external region, where no massive structures are allowed and only an extremely low stray field is admitted. In this paper we report the main technical challenges and the relevant design options providing a comprehensive design for the magnet of the SHiP Scattering and Neutrino Detector.


title: DNN-Aided Block Sparse Bayesian Learning for User Activity Detection and Channel Estimation in Grant-Free Non-Orthogonal Random Access

abstract: In the upcoming Internet-of-Things (IoT) era, the communication is often featured by massive connection, sporadic transmission, and small-sized data packets, which poses new requirements on the delay expectation and resource allocation efficiency of the Random Access (RA) mechanisms of the IoT communication stack. A grant-free non-orthogonal random access (NORA) system is considered in this paper, which could simultaneously reduce the access delay and support more Machine Type Communication (MTC) devices with limited resources. In order to address the joint user activity detection (UAD) and channel estimation (CE) problem in the grant-free NORA system, we propose a deep neural network-aided message passing-based block sparse Bayesian learning (DNN-MP-BSBL) algorithm. In the DNN-MP-BSBL algorithm, the iterative message passing process is transferred from a factor graph to a deep neural network (DNN). Weights are imposed on the messages in the DNN and trained to minimize the estimation error. It is shown that the trained weights could alleviate the convergence problem of the MP-BSBL algorithm, especially on crowded RA scenarios. Simulation results show that the proposed DNN-MP-BSBL algorithm could improve the UAD and CE accuracy with a smaller number of iterations, indicating its advantages for low-latency grant-free NORA systems.


title: On the prime zeta function and the Riemann hypothesis

abstract: By considering the prime zeta function, the author intended to demonstrate in that the Riemann zeta function zeta(s) does not vanish for Re(s)>1/2, which would have proven the Riemann hypothesis. However, he later realised that the proof of "Theorem 3" is fundamentally flawed. The main tools of our argument are: bounds and oscillation theorems for the prime counting function, classical properties of Dirichlet series and the identity theorem for real-analytic functions.


title: Approximate evolution for a system composed by two coupled Jaynes-Cummings Hamiltonians

abstract: In this work we construct an approximate time evolution operator for a system composed by two coupled Jaynes-Cummings Hamiltonians. We express the full time evolution operator as a product of exponentials and we analyze the validity of our approximations contrasting our analytical results with those obtained by purely numerical methods.


title: Hard X-ray Excess from the Magnificent Seven Neutron Stars

abstract: We report significant hard X-ray excesses in the energy range 2-8 keV for two nearby isolated neutron stars RX J1856.6-3754 and RX J0420.0-5022. These neutron stars have previously been observed in soft X-rays to have nearly thermal spectra at temperatures ~100 eV, which are thought to arise from the warm neutron star surfaces. We find nontrivial hard X-ray spectra well above the thermal surface predictions with archival data from the XMM-Newton and Chandra X-ray telescopes. We analyze possible systematic effects that could generate such spurious signals, such as nearby X-ray point sources and pileup of soft X-rays, but we find that the hard X-ray excesses are robust to these systematics. We also investigate possible sources of hard X-ray emission from the neutron stars and find no satisfactory explanation with known mechanisms, suggesting that a novel source of X-ray emission is at play. We do not find high-significance hard X-ray excesses from the other five Magnificent Seven isolated neutron stars.


title: The Shadow of a Spherically Accreting Black Hole

abstract: We explore a simple spherical model of optically thin accretion on a Schwarzschild black hole, and study the properties of the image as seen by a distant observer. We show that a dark circular region in the center --- a shadow --- is always present. The outer edge of the shadow is located at the photon ring radius $b_{\rm ph} \equiv \sqrt{27}r_g$, where $r_g=GM/c^2$ is the gravitational radius of the accreting mass $M$. The location of the shadow edge is independent of the inner radius at which the accreting gas stops radiating. The size of the observed shadow is thus a signature of the spacetime geometry and it is hardly influenced by accretion details. We briefly discuss the relevance of these results for the Event Horizon Telescope image of the supermassive black hole in M87.


title: LEO-Py: Estimating likelihoods for correlated, censored, and uncertain data with given marginal distributions

abstract: Data with uncertain, missing, censored, and correlated values are commonplace in many research fields including astronomy. Unfortunately, such data are often treated in an ad hoc way in the astronomical literature potentially resulting in inconsistent parameter estimates. Furthermore, in a realistic setting, the variables of interest or their errors may have non-normal distributions which complicates the modeling. I present a novel approach to compute the likelihood function for such data sets. This approach employs Gaussian copulas to decouple the correlation structure of variables and their marginal distributions resulting in a flexible method to compute likelihood functions of data in the presence of measurement uncertainty, censoring, and missing data. I demonstrate its use by determining the slope and intrinsic scatter of the star forming sequence of nearby galaxies from observational data. The outlined algorithm is implemented as the flexible, easy-to-use, open-source Python package LEO-Py.


title: The evolution of rest-frame UV properties, Lya EWs and the SFR-Stellar mass relation at z~2-6 for SC4K LAEs

abstract: We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $\sigma$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.


title: High angular resolution ALMA images of dust and molecules in the SN 1987A ejecta

abstract: We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.


title: Dirac neutrinos from Peccei-Quinn symmetry: a fresh look at the axion

abstract: We show that a very simple solution to the strong CP problem naturally leads to Dirac neutrinos. Small effective neutrino masses emerge from a type I Dirac seesaw mechanism. Neutrino mass limits probe the axion parameters in regions currently inaccessible to conventional searches.


title: 3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$\alpha$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

abstract: We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.


title: Classifying Calabi-Yau threefolds using infinite distance limits

abstract: We present a novel way to classify Calabi-Yau threefolds by systematically studying their infinite volume limits. Each such limit is at infinite distance in Kahler moduli space and can be classified by an associated limiting mixed Hodge structure. We then argue that the such structures are labeled by a finite number of degeneration types that combine into a characteristic degeneration pattern associated to the underlying Calabi-Yau threefold. These patterns provide a new invariant way to present crucial information encoded in the intersection numbers of Calabi-Yau threefolds. For each pattern, we also introduce a Hasse diagram with vertices representing each, possibly multi-parameter, decompactification limit and explain how to read off properties of the Calabi-Yau manifold from this graphical representation. In particular, we show how it can be used to count elliptic, K3, and nested fibrations and determine relations of elliptic fibrations under birational equivalence. We exemplify this for hypersurfaces in toric ambient spaces as well as for complete intersections in products of projective spaces.


title: Understanding Lyman-alpha Nebulae at Low Redshift I: The Sizes, Powering, and Kinematics of "Green Bean" Galaxies

abstract: A new but rare sample of spatially extended emission line nebulae, nicknamed "Green Beans", was discovered at z~0.3 thanks to strong [OIII] emission, and subsequently shown to be local cousins of the Lyman-alpha (Lya) nebulae found at high redshift. Here we use follow-up APO/DIS spectroscopy to better understand how these low redshift Lya nebulae compare to other populations of strong emission line sources. Our spectroscopic data show that low-z Lya nebulae have AGN-like emission line ratios, relatively narrow line widths (FWHM<1000 km s^{-1}), and emission line kinematics resembling those of Type 2 AGN at the same redshift, confirming that they are powered by Type 2 AGN with typical ionizing continua. While low-z Lya nebulae are larger and less concentrated than compact, star-forming Green Pea galaxies, we find that they resemble typical Type 2 AGN in terms of r-band concentration and size. Based on this pilot study, low-z Lya nebulae appear to be a subset of Type 2 AGN with bluer optical continua and high [OIII] equivalent widths but with comparable sizes and similar [OIII] kinematics. These characteristics may simply reflect the fact that low-z Lya nebulae are drawn from the high luminosity end of the Type 2 AGN distribution, with higher nuclear activity driving higher [OIII] equivalent widths and more central star formation leading to bluer optical continua. Deeper spectroscopic follow-up of the full sample will shed further light on these issues and on the relationship between these low-z Lya nebulae and the Lya nebula population at high redshift.


title: The Possible Astrometric Signature of a Planetary-mass Companion to the Nearby Young Star TW Piscis Austrini (Fomalhaut B): Constraints from Astrometry, Radial Velocities, and Direct Imaging

abstract: We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.


title: Repeated Gravitational Lensing of Gravitational Waves in Hierarchical Black Hole Triples

abstract: We consider binary black holes (BBHs) in a hierarchical triple system where a more compact, less-massive binary is emitting detectable gravitational waves (GWs), and the tertiary is a supermassive BH at the center of a nuclear star cluster. As previous works have shown, the orbital motion of the outer binary can generate a detectable relativistic Doppler boost of the GWs emitted by the orbiting inner binary. We show here that for outer-binary orbits with a period of order one year, there can be a non-negligible probability for repeated gravitational lensing of the GWs emitted by the inner binary. Repeating gravitational lensing events could be detected by the LISA observatory as periodic GW amplitude spikes before the BBH enters the LIGO band. Such a detection would confirm the origin of some BBH mergers in nuclear star clusters. GW lensing also offers new testing grounds for strong gravity.


title: ZTF Early Observations of Type Ia Supernovae I: Properties of the 2018 Sample

abstract: Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).


title: Four (Super)luminous Supernovae from the First Months of the ZTF Survey

abstract: We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population.


title: Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): How Well-centered Are Low-z AGN?

abstract: Off-nucleus active galactic nuclei (AGN) can be signposts of inspiraling supermassive black holes (SMBHs) on galactic scales, or accreting SMBHs recoiling after the coalescence of a SMBH binary or slingshot from three-body interactions. Because of the stochastic variability of AGN, the measured photocenter of an unresolved AGN-host system will display astrometric jitter that depends on the off-nucleus distance of the AGN, the total photometric variability of the system, and the AGN-host contrast. Here we use the precision astrometry from Gaia DR2 to constrain the off-nucleus population of a low-redshift (0.3<z<0.8) sample of unobscured broad-line AGN drawn from the SDSS with significant host contribution and photometric variability. We find that Gaia DR2 already provides strong constraints on the projected off-nucleus distance in the sub-kpc regime at these redshifts: 99%, 90% and 40% of AGN must be well-centered to <1 kpc, <500 pc and <100 pc, respectively. Limiting the sample to the most variable subset constrains >99% of AGN to be well-centered below 500 parsec. These results suggest that genuine off-nucleus AGN (offset by > a few hundred pc) must be rare at low redshift. Future Gaia releases of time series of photocenter and flux measurements, improved treatments for extended sources and longer baselines will further tighten these constraints, and enable a systematic full-sky search for rare off-nucleus AGN on ~ 10-1000 pc scales.


title: Probing quantum criticality using nonlinear Hall effect in a metallic Dirac system

abstract: Interaction driven symmetry breaking in a metallic (doped) Dirac system can manifest in the spontaneous gap generation at the nodal point buried below the Fermi level. Across this transition linear conductivity remains finite making its direct observation difficult in linear transport. We propose the nonlinear Hall effect as a direct probe of this transition when inversion symmetry is broken. Specifically, for a two-dimensional Dirac material with a tilted low-energy dispersion, we first predict a transformation of the characteristic inter-band resonance peak into a non-Lorentzian form in the collisionless regime. Furthermore, we show that inversion-symmetry breaking quantum phase transition is controlled by an exotic tilt-dependent line of critical points. As this line is approached from the ordered side, the nonlinear Hall conductivity is suppressed owing to the scattering between the strongly coupled incoherent fermionic and bosonic excitations. Our results should motivate further studies of nonlinear responses in strongly interacting Dirac materials.


title: Consistent QFT description of non-standard neutrino interactions

abstract: Neutrino oscillations are precision probes of new physics beyond the Standard Model. Apart from neutrino masses and mixings, they are also sensitive to possible deviations of low-energy interactions between quarks and leptons from the Standard Model predictions. In this paper we develop a systematic description of such non-standard interactions (NSI) in oscillation experiments within the quantum field theory framework. We calculate the event rate and oscillation probability in the presence of general NSI, starting from the effective field theory (EFT) in which new physics modifies the flavor or Lorentz structure of charged-current interactions between leptons and quarks. We also provide the matching between the EFT Wilson coefficients and the widely used simplified quantum-mechanical approach, where new physics is encoded in a set of production and detection NSI parameters. Finally, we discuss the consistency conditions for the standard NSI approach to correctly reproduce the quantum field theory result.


title: Exotic Spin-Dependent Forces from a Hidden Sector

abstract: New dynamics from hidden sectors may manifest as long-range forces between visible matter particles. The well-known case of Yukawa-like potentials occurs via the exchange of a single virtual particle. However, more exotic behavior is also possible. We present three classes of exotic potentials that are generated by relativistic theories: (i) quantum forces from the loop-level exchange of two virtual particles, (ii) conformal forces from a conformal sector, and (iii) emergent forces from degrees of freedom that only exist in the infrared regime of the theory. We discuss the complementarity of spin-dependent force searches in an effective field theory framework. We identify well-motivated directions to search for exotic spin-dependent forces.


title: Classification of RCFT from Holomorphic Modular Bootstrap: A Status Report

abstract: Following the initial proposal in 1988, there has been much progress in classifying Rational Conformal Field Theories in 2 dimensions from the Holomorphic Bootstrap approach. This method starts by postulating a generic holomorphic Modular Linear Differential Equation of a given order and imposing the requirement of non-negative integrality of the coefficients in the series expansion of the solutions, which are then identified as admissible characters, from which a modular-invariant partition function is constructed. In this short note, the status of this project is summarised.


title: SMArT: Training Shallow Memory-aware Transformers for Robotic Explainability

abstract: The ability to generate natural language explanations conditioned on the visual perception is a crucial step towards autonomous agents which can explain themselves and communicate with humans. While the research efforts in image and video captioning are giving promising results, this is often done at the expense of the computational requirements of the approaches, limiting their applicability to real contexts. In this paper, we propose a fully-attentive captioning algorithm which can provide state-of-the-art performances on language generation while restricting its computational demands. Our model is inspired by the Transformer model and employs only two Transformer layers in the encoding and decoding stages. Further, it incorporates a novel memory-aware encoding of image regions. Experiments demonstrate that our approach achieves competitive results in terms of caption quality while featuring reduced computational demands. Further, to evaluate its applicability on autonomous agents, we conduct experiments on simulated scenes taken from the perspective of domestic robots.


title: Fermionic pole-skipping in holography

abstract: We examine thermal Green's functions of fermionic operators in quantum field theories with gravity duals. The calculations are performed on the gravity side using ingoing Eddington-Finkelstein coordinates. We find that at negative imaginary Matsubara frequencies and special values of the wavenumber, there are multiple solutions to the bulk equations of motion that are ingoing at the horizon and thus the boundary Green's function is not uniquely defined. At these points in Fourier space a line of poles and a line of zeros of the correlator intersect. We analyze these `pole-skipping' points in three-dimensional asymptotically anti-de Sitter spacetimes where exact Green's functions are known. We then generalize the procedure to higher-dimensional spacetimes. We also discuss the special case of a fermion with half-integer mass in the BTZ background. We discuss the implications and possible generalizations of the results.


title: Indigenizing the next decade of astronomy in Canada

abstract: (Abridged) The Truth and Reconciliation Commission of Canada published its calls to action in 2015 with 94 recommendations. Many of these 94 recommendations are directly related to education, language, and culture, some of which the Canadian Astronomy community can address and contribute to as part of reconciliation. The Canadian Astronomy community has an additional obligation since it benefits from facilities on Indigenous territories across Canada and the world. Furthermore, Indigenous people are still underrepresented at all levels in Canadian astronomy. The purpose of this Community Paper is to develop recommendations for the Canadian astronomy community to support Indigenous inclusion in the science community, support Indigenous learning by developing Indigenous-based learning materials and facilitate access to professionals and science activities, and to recognize and acknowledge the great contributions of Indigenous communities to our science activities. As part of this work we propose the ten following recommendations for CASCA as an organization and throughout this Community Paper we will include additional recommendations for individuals: astronomers, students and academics.


title: Colored Multipermutations and a Combinatorial Generalization of Worpitzky's Identity

abstract: Worpitzky's identity expresses $n^p$ in terms of the Eulerian numbers and binomial coefficients: $$n^p = \sum_{i=0}^{p-1} \genfrac<>{0pt}{}{p}{i} \binom{n+i}{p}.$$ Pita-Ruiz recently defined numbers $A_{a,b,r}(p,i)$ implicitly to satisfy a generalized Worpitzky identity $$\binom{an+b}{r}^p = \sum_{i=0}^{rp} A_{a,b,r}(p,i) \binom{n+rp-i}{rp},$$ and asked whether there is a combinatorial interpretation of the numbers $A_{a,b,r}(p,i)$. We provide such a combinatorial interpretation by defining a notion of descents in colored multipermutations, and then proving that $A_{a,b,r}(p,i)$ is equal to the number of colored multipermutations of $\{1^r, 2^r, \ldots, p^r\}$ with $a$ colors and $i$ weak descents. We use this to give combinatorial proofs of several identities involving $A_{a,b,r}(p,i)$, including the aforementioned generalized Worpitzky identity.


title: Investigating the Hubble Constant Tension -- Two Numbers in the Standard Cosmological Model

abstract: The current Hubble constant tension is usually presented by comparing constraints on $H_0$ only. However, the post-recombination background cosmic evolution is determined by two parameters in the standard $\Lambda$CDM model, the Hubble constant ($H_0$) and today's matter energy fraction ($\Omega_{\rm{m}}$). If we therefore compare all constraints individually in the $H_0$-$\Omega_{\rm{m}}$ plane, (1) various constraints can be treated as independently as possible, (2) single-sided constraints are easier to consider, (3) compatibility among different constraints can be viewed in a more robust way, (4) the model dependence of each constraint is clear, and (5) whether or not a nonstandard model is able to reconcile all constraints in tension can be seen more effectively. We perform a systematic comparison of different constraints in the $H_0$-$\Omega_{\rm{m}}$ space based on a flat $\Lambda$CDM model, treating them as separately as possible. Constraints along different degeneracy directions consistently overlap in one region of the space, with the local measurement from Cepheid variable-calibrated supernovae being the most outlying, followed by the time-delay strong-lensing result. Considering the possibility that some nonstandard physics may reconcile the constraints, we provide a general discussion on nonstandard models with modifications at high, mid, or low redshifts, and the effect of local environmental factors. Due to the different responses of individual constraints to a modified model, it is not easy for nonstandard models to reconcile all constraints if none of them have unaccounted-for systematic effects.


title: Universal Upper Bound on the Inflationary Energy Scale from the Trans-Planckian Censorship Conjecture

abstract: We study the constraint on the inflationary energy scale from the recently proposed Trans-Planckian Censorship Conjecture (TCC). We find a universal upper bound on the inflationary Hubble expansion rate $H_\text{inf}$ which is solely determined by the reheating temperature $T_\text{rh}$: $H_\text{inf}/M_\text{Pl}\lesssim T_0/T_\text{rh}$, where $T_0$ is the photon temperature today. The upper limit can be saturated by a post-inflationary oscillatory stage with the critical equation-of-state parameter $w\approx-1/3$, or by an inflation model with multiple stages. In the lowest reheating temperature required for big bang nucleosynthesis, the upper bound on the tensor-to-scalar ratio $r$ at the CMB scales is $r\lesssim10^{-8}$, which can be realized in many string-inspired inflation models.


title: Sequential Phase Transitions and Transient Structured Fluctuations in Two-Dimensional Systems with a High-Density Kagome Lattice Phase

abstract: We report studies of the phase diagrams of two two-dimensional systems composed of particles with repulsive pair potentials, each of which supports a high-density Kagome lattice phase. The commonalities in the phase diagrams of the systems we have studied and other 2D systems suggests the existence of a universal mechanism driving all to favor a similar series of packing arrangements as the density is increased, but the simulations considered show that the only such general rule proposed, namely the S\"uto theorem, is not a necessary condition for the support of multiple distinct lattice structures by a particular pair potential. We also find that close to the liquid-to-Kagome phase transition the transient structured fluctuations in the liquid have Kagome symmetry whereas deeper in the liquid phase the fluctuations have hexagonal symmetry. As the deviation of the liquid density from the transition density decreases fluctuations with hexagonal symmetry are replaced with those with Kagome symmetry. When the transition is string phase-to-Kagome phase the transient structured fluctuations in the string phase near the transition have both six-fold and other than six-fold symmetries, with stronger preference for six-fold symmetry in the Truskett system than in the Torquato system. The path of the string-to-Kagome transition in the Truskett system involves intermediate honeycomb configurations that subsequently buckle to form a Kagome lattice. The path of the string-to-Kagome transition in the Torquato system suggests that the Kagome phase is formed by coiled strings merging together; increasing density generates a Kagome phase with imperfections such as 8-particle rings.


title: On the size of the CO-depletion radius in the IRDC G351.77-0.51

abstract: An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is $<20$ K and the number density of molecular hydrogen exceeds a few $\times$10$^{4}$ cm$^{-3}$. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77-0.51. Continuum observations performed with the $Herschel$ $Space$ $Observatory$ and the $LArge$ $APEX$ $BOlometer$ $CAmera$, together with APEX C$^{18}$O and C$^{17}$O J=2$\rightarrow$1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ cm$^{-3}$. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs.


title: Physical properties of terrestrial planets and water delivery in the habitable zone using N-body simulations with fragmentation

abstract: The goal of this research is to study how the fragmentation of planetary embryos can affect the physical and dynamical properties of terrestrial planets around solar-type stars. Our work focuses on the formation and evolution of planets and water delivery in the Habitable Zone. We distinguish class A and class B HZ planets, which have an accretion seed initially located inside and beyond the snowline, respectively. We develop an Nbody integrator that incorporates fragmentation and hit-and-run collisions, which is called D3 N-body code. From this, we perform 46 numerical simulations of planetary accretion in systems that host two gaseous giants like Jupiter and Saturn. We compare two sets of 23 N-body simulations, one of which includes a realistic collisional treatment and the other one models all impacts as perfect mergers. The final masses of the HZ planets formed in runs with fragmentation are about 15\%-20\% smaller than those obtained without fragmentation. As for the class A HZ planets, those formed in simulations without fragmentation experience very significant increases in mass respect to their initial values, while the growth of those produced in runs with fragmentation is less relevant. We remark that the fragments play a secondary role in the masses of the class A HZ planets, providing less than 30\% of their final values. In runs without fragmentation, the final fraction of water of the class A HZ planets keeps the initial value since they do not accrete water-rich embryos. In runs with fragmentation, the final fraction of water of such planets strongly depends on the model used to distribute the water after each collision. The class B HZ planets do not show significant differences concerning their final water contents in runs with and without fragmentation. From this, we find that the collisional fragmentation is not a barrier to the survival of water worlds in the HZ.


title: H$\alpha$ Morphologies of Star Clusters: A LEGUS study of HII region evolution timescales and stochasticity in low mass clusters

abstract: The morphology of HII regions around young star clusters provides insight into the timescales and physical processes that clear a cluster's natal gas. We study ~700 young clusters (<10Myr) in three nearby spiral galaxies (NGC 7793, NGC 4395, and NGC 1313) using Hubble Space Telescope (HST) imaging from LEGUS (Legacy ExtraGalactic Ultraviolet Survey). Clusters are classified by their H$\alpha$ morphology (concentrated, partially exposed, no-emission) and whether they have neighboring clusters (which could affect the clearing timescales). Through visual inspection of the HST images, and analysis of ages, reddenings, and stellar masses from spectral energy distributions fitting, together with the (U-B), (V-I) colors, we find: 1) the median ages indicate a progression from concentrated (~3 Myr), to partially exposed (~4 Myr), to no H$\alpha$ emission (>5Myr), consistent with the expected temporal evolution of HII regions and previous results. However, 2) similarities in the age distributions for clusters with concentrated and partially exposed H$\alpha$ morphologies imply a short timescale for gas clearing (<1Myr). 3) our cluster sample's median mass is ~1000 M, and a significant fraction (~20%) contain one or more bright red sources (presumably supergiants), which can mimic reddening effects. Finally, 4) the median E(B-V) values for clusters with concentrated H$\alpha$ and those without H$\alpha$ emission appear to be more similar than expected (~0.18 vs. ~0.14, respectively), but when accounting for stochastic effects, clusters without H$\alpha$ emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H$\alpha$ morphological class. Composite isolated clusters also reveal a color and age progression for H$\alpha$ morphological classes, consistent with analysis of the individual clusters.


title: Accurate high-throughput screening of I-II-V 8-electron Half-Heusler compounds for renewable-energy applications

abstract: Renewable energy resources have emerged as the best alternatives to fossil fuel energy which are rapidly declining with time. Here, eight valence-electron count Half-Heusler(HH) alloys have been studied using reliable first principles calculations in the search of potential candidates for renewable energy applications like thermoelectric (TE), solar harvesting, topological insulator (TI) and transparent conductor (TC) applications. The initial screening parameters used for our study are chemical and thermal stability, band gap, nature of bandgap and band inversion strength. We have performed quasistatic G0W0 calculation starting from HSE groundstate wavefunction to predict the most accurate estimation of bandgap for these class of compounds. A total of 960 compounds were simulated. 121 out of 960 compounds were found to be thermally and chemically stable. 31 compounds with bandgap less than 1.5 eV were studied for thermoelectric application out of which 13 compounds were found to show thermoelectric figure of merit ZT > 0.7 for both p-type and n-type conduction. 30 compounds with band gap 1-1.8 eV were studied for optoelectronic application out of which 13 compounds were found to show Spectroscopic Limited Maximum Efficiency (SLME) more than 20%, comparable to existing state of the art materials. 21 compounds were found to show band inversion at ambient conditions which is a necessary condition for topological insulators. The surface band structure calculations for one of the promising candidate was done to check robustness of the topological behaviour. 29 compounds were found to have bandgap more than 2 eV which are promoted for transparent conductor applications with further band engineering. We strongly believe that our calculations will give useful insights to experimentalists for synthesizing and investigating proposed compounds for different energy applications.


title: The Effects of the Problem Hamiltonian Parameters on the Minimum Spectral Gap in Adiabatic Quantum Optimization

abstract: We study the relation between the Ising problem Hamiltonian parameters and the minimum spectral gap (min-gap) of the system Hamiltonian in the Ising-based quantum annealer. The main argument we use in this paper to assess the performance of a QA algorithm is the presence or absence of an anti-crossing during quantum evolution. For this purpose, we introduce a new parametrization definition of the anti-crossing. Using the Maximum-weighted Independent Set (MIS) problem in which there are flexible parameters (energy penalties J between pairs of edges) in an Ising formulation as the model problem, we construct examples to show that by changing the value of J, we can change the quantum evolution from one that has an anti-crossing (that results in an exponential small min-gap) to one that does not have, or the other way around, and thus drastically change (increase or decrease) the min-gap. However, we also show that by changing the value of $J$ alone, one can not avoid the anti-crossing. We recall a polynomial reduction from an Ising problem to an MIS problem to show that the flexibility of changing parameters without changing the problem to be solved can be applied to any Ising problem. As an example, we show that by such a reduction alone, it is possible to remove the anti-crossing and thus increase the min-gap. Our anti-crossing definition is necessarily scaling invariant as scaling the problem Hamiltonian does not change the nature (i.e. presence or absence) of an anti-crossing. As a side note, we show exactly how the min-gap is scaled if we scale the problem Hamiltonian by a constant factor.


title: A Distributed and Integrated Method of Moments for High-Dimensional Correlated Data Analysis

abstract: This paper is motivated by a regression analysis of electroencephalography (EEG) neuroimaging data with high-dimensional correlated responses with multi-level nested correlations. We develop a divide-and-conquer procedure implemented in a fully distributed and parallelized computational scheme for statistical estimation and inference of regression parameters. Despite significant efforts in the literature, the computational bottleneck associated with high-dimensional likelihoods prevents the scalability of existing methods. The proposed method addresses this challenge by dividing responses into subvectors to be analyzed separately and in parallel on a distributed platform using pairwise composite likelihood. Theoretical challenges related to combining results from dependent data are overcome in a statistically efficient way using a meta-estimator derived from Hansen's generalized method of moments. We provide a rigorous theoretical framework for efficient estimation, inference, and goodness-of-fit tests. We develop an R package for ease of implementation. We illustrate our method's performance with simulations and the analysis of the EEG data, and find that iron deficiency is significantly associated with two auditory recognition memory related potentials in the left parietal-occipital region of the brain.


title: Semi-analytic Evaluation of 1, 2 and 3-Electron Coulomb Integrals with Gaussian expansion of Distance Operators W= R$_{C1}^{-n}$R$_{D1}^{-m}$, R$_{C1}^{-n}$r$_{12}^{-m}$, r$_{12}^{-n}$r$_{13}^{-m}$

abstract: The equations derived help to evaluate semi-analytically (mostly for k=1,2 or 3) the important Coulomb integrals Int rho(r1)...rho(rk) W(r1,...,rk) dr1...drk, where the one-electron density, rho(r1), is a linear combination (LC) of Gaussian functions of position vector variable r1. It is capable to describe the electron clouds in molecules, solids or any media/ensemble of materials, weight W is the distance operator indicated in the title. R stands for nucleus-electron and r for electron-electron distances. The n=m=0 case is trivial, the (n,m)=(1,0) and (0,1) cases, for which analytical expressions are well known, are widely used in the practice of computation chemistry (CC) or physics, and analytical expressions are also known for the cases n,m=0,1,2. The rest of the cases - mainly with any real (integer, non-integer, positive or negative) n and m - needs evaluation. We base this on the Gaussian expansion of |r|^-u, of which only the u=1 is the physical Coulomb potential, but the u.ne.1 cases are useful for (certain series based) correction for (the different) approximate solutions of Schrodinger equation, for example, in its wave-function corrections or correlation calculations. Solving the related linear equation system (LES), the expansion |r|^-u about equal to SUM(k=0toL)SUM(i=1toM) Cik r^2k exp(-Aik r^2) is analyzed for |r| = r12 or RC1 with least square fit (LSF) and modified Taylor expansion. These evaluated analytic expressions for Coulomb integrals (up to Gaussian function integrand and the Gaussian expansion of |r|^-u) are useful for the manipulation with higher moments of inter-electronic distances via W, even for approximating Hamiltonian.


title: Walking Dynamics Guaranteed

abstract: We report evidence for a continuous transition from an infrared conformal phase to a chirally broken one in four dimensions. We study a model with two Dirac fermions in the adjoint representation of an SU(2) gauge interaction and a chirally symmetric four-fermion interaction. At large four-fermion coupling, the model goes through a transition into a chirally broken phase and infrared conformality is lost. We show strong evidence that this transition is continuous, which would guarantee walking dynamics within the scaling region in the chirally broken phase.


title: Leveraging Vision Reconstruction Pipelines for Satellite Imagery

abstract: Reconstructing 3D geometry from satellite imagery is an important topic of research. However, disparities exist between how this 3D reconstruction problem is handled in the remote sensing context and how multi-view reconstruction pipelines have been developed in the computer vision community. In this paper, we explore whether state-of-the-art reconstruction pipelines from the vision community can be applied to the satellite imagery. Along the way, we address several challenges adapting vision-based structure from motion and multi-view stereo methods. We show that vision pipelines can offer competitive speed and accuracy in the satellite context.


title: Gravitational Fluctuations as an Alternative to Inflation II. CMB Angular Power Spectrum

abstract: Power spectra always play an important role in the theory of inflation. In particular, the ability to reproduce the galaxy matter power spectrum and the CMB temperature angular power spectrum coefficients to high accuracy is often considered a triumph of inflation. In our previous work, we presented an alternative explanation for the matter power spectrum based on nonperturbative quantum field-theoretical methods applied to Einstein's gravity, instead of inflation models based on scalar fields. In this work, we review the basic concepts and provide further in-depth investigations. We first update the analysis with more recent data sets and error analysis, and then extend our predictions to the CMB angular spectrum coefficients, which we did not consider previously. Then we investigate further the potential freedoms and uncertainties associated with the fundamental parameters that are part of this picture, and show how recent cosmological data provides significant constraints on these quantities. Overall, we find good general consistency between theory and data, even potentially favoring the gravitationally-motivated picture at the largest scales. We summarize our results by outlining how this picture can be tested in the near future with increasingly accurate astrophysical measurements.


title: Spectroscopic signatures of phonons in high pressure superconducting hydrides

abstract: The discovery of superconductivity at 203K in SH$_3$ is an important step toward higher values of $T_c$. Predictions based on state-of-the-art DFT for the electronic structure, including one preceding experimental confirmation, showed the mechanism to be the electron-phonon interaction. This was confirmed in optical spectroscopy measurements. For photon energies between $\sim 450$ and 600 meV in SH$_3$, the reflectance in the superconducting state is below that in its normal state. This difference decreases as $T$ approaches $T_c$. Decreasing absorption with increasing $T$ is opposite to what is expected in ordinary metals. Such an anomalous behavior can be traced back to the energy dependence of the superconducting density of states which is highly peaked at the energy gap value $\Delta$ but decays back to the constant normal state value as energy is increased, on a scale of a few $\Delta$, or by increasing $T$ towards $T=T_c$. The process of phonon-assisted optical absorption is encoded with a knowledge of the $T$-dependence of $\Delta$, the order parameter of the superconducting state. Should the energy of the phonon involved be very large, of order 200 meV or more, this process offers the possibility of observing the closing of the superconducting order parameter with $T$ at correspondingly very large energies. The very recent experimental observation of a $T_c\simeq 250$ K in LaH$_{10}$ has further heightened interest in the hydrides. We compare the relevant phonon structure seen in optics with related features in the real and imaginary part of the frequency dependent gap, quasiparticle density of states, reflectance, absorption, and optical scattering rate. The phonon structures all carry information on the $T_c$ value and the $T$-dependence of the order parameter, and can be used to confirm that the mechanism involved in superconductivity is the electron-phonon interaction.


title: A simple quantum test for smooth horizons

abstract: We develop a new test that provides a necessary condition for a quantum state to be smooth in the vicinity of a null surface: near-horizon modes that can be defined locally near any patch of the null surface must be correctly entangled with each other and with their counterparts across the surface. This test is considerably simpler to implement than a full computation of the renormalized stress-energy tensor. We apply this test to Reissner-Nordstrom black holes in asymptotically anti-de Sitter space and provide numerical evidence that the inner horizon of such black holes is singular in the Hartle-Hawking state. We then consider BTZ black holes, where we show that our criterion for smoothness is satisfied as one approaches the inner horizon from outside. This results from a remarkable conspiracy between the properties of mode-functions outside the outer horizon and between the inner and outer horizon. Moreover, we consider the extension of spacetime across the inner horizon of BTZ black holes and show that it is possible to define modes behind the inner horizon that are correctly entangled with modes in front of the inner horizon. Although this provides additional suggestions for the failure of strong cosmic censorship, we lay out several puzzles that must be resolved before concluding that the inner horizon will be traversable.


title: Rekall: Specifying Video Events using Compositions of Spatiotemporal Labels

abstract: Many real-world video analysis applications require the ability to identify domain-specific events in video, such as interviews and commercials in TV news broadcasts, or action sequences in film. Unfortunately, pre-trained models to detect all the events of interest in video may not exist, and training new models from scratch can be costly and labor-intensive. In this paper, we explore the utility of specifying new events in video in a more traditional manner: by writing queries that compose outputs of existing, pre-trained models. To write these queries, we have developed Rekall, a library that exposes a data model and programming model for compositional video event specification. Rekall represents video annotations from different sources (object detectors, transcripts, etc.) as spatiotemporal labels associated with continuous volumes of spacetime in a video, and provides operators for composing labels into queries that model new video events. We demonstrate the use of Rekall in analyzing video from cable TV news broadcasts, films, static-camera vehicular video streams, and commercial autonomous vehicle logs. In these efforts, domain experts were able to quickly (in a few hours to a day) author queries that enabled the accurate detection of new events (on par with, and in some cases much more accurate than, learned approaches) and to rapidly retrieve video clips for human-in-the-loop tasks such as video content curation and training data curation. Finally, in a user study, novice users of Rekall were able to author queries to retrieve new events in video given just one hour of query development time.


title: Stochastic Model Predictive Control of Autonomous Systems with Non-Gaussian Correlated Uncertainty

abstract: Many systems such as autonomous vehicles and quadrotors are subject to parametric uncertainties and external disturbances. These uncertainties can lead to undesired performance degradation and safety issues. Therefore, it is important to design robust control strategies to safely regulate the dynamics of a system. This paper presents a novel framework for chance-constrained stochastic model predictive control of dynamic systems with non-Gaussian correlated probabilistic uncertainties. We develop a new stochastic Galerkin method to propagate the uncertainties using a new type of basis functions and an optimization-based quadrature rule. This formulation can easily handle non-Gaussian correlated uncertainties that are beyond the capability of generalized polynomial chaos expansions. The new stochastic Galerkin formulation enables us to convert a chance-constraint stochastic model predictive control problem into a deterministic one. We verify our approach by several stochastic control tasks, including obstacle avoidance, vehicle path following, and quadrotor reference tracking.


title: A Locally Adaptive Bayesian Cubature Method

abstract: Bayesian cubature (BC) is a popular inferential perspective on the cubature of expensive integrands, wherein the integrand is emulated using a stochastic process model. Several approaches have been put forward to encode sequential adaptation (i.e. dependence on previous integrand evaluations) into this framework. However, these proposals have been limited to either estimating the parameters of a stationary covariance model or focusing computational resources on regions where large values are taken by the integrand. In contrast, many classical adaptive cubature methods focus computational resources on spatial regions in which local error estimates are largest. The contributions of this work are three-fold: First, we present a theoretical result that suggests there does not exist a direct Bayesian analogue of the classical adaptive trapezoidal method. Then we put forward a novel BC method that has empirically similar behaviour to the adaptive trapezoidal method. Finally we present evidence that the novel method provides improved cubature performance, relative to standard BC, in a detailed empirical assessment.


title: Variation of shear moduli across superconducting phase transitions

abstract: We study how shear moduli of a correlated metal change across superconducting phase transitions. Using a microscopic theory we explain why for most classes of superconductors this change is small. The Fe-based and the A15 systems are notable exceptions where the change is boosted by five orders of magnitude. We show that this boost is a consequence of enhanced nematic correlation. The theory explains the unusual temperature dependence of the orthorhombic shear and the back-bending of the nematic transition line in the superconducting phase of the Fe-based systems.


title: Identifying causal effects in maximally oriented partially directed acyclic graphs

abstract: We develop a necessary and sufficient causal identification criterion for maximally oriented partially directed acyclic graphs (MPDAGs). MPDAGs as a class of graphs include directed acyclic graphs (DAGs), completed partially directed acyclic graphs (CPDAGs), and CPDAGs with added background knowledge. As such, they represent the type of graph that can be learned from observational data and background knowledge under the assumption of no latent variables. Our identification criterion can be seen as a generalization of the g-formula of Robins (1986). We further obtain a generalization of the truncated factorization formula (Pearl, 2009) and compare our criterion to the generalized adjustment criterion of Perkovi\'c et al. (2017) which is sufficient, but not necessary for causal identification.


title: Construction of quantum Dirac observables and the emergence of WKB time

abstract: We describe a method of construction of gauge-invariant operators (Dirac observables or ``evolving constants of motion'') from the knowledge of the eigenstates of the gauge generator in time-reparametrization invariant mechanical systems. These invariant operators evolve unitarily with respect to an arbitrarily chosen time variable. We emphasize that the dynamics is relational, both in the classical and quantum theories. In this framework, we show how the ``emergent Wentzel-Kramers-Brillouin time'' often employed in quantum cosmology arises from a weak-coupling expansion of invariant transition amplitudes, and we illustrate an example of singularity avoidance in a vacuum Bianchi I (Kasner) model.


title: Universality for 1 d random band matrices

abstract: We consider 1d random Hermitian $N\times N$ block band matrices consisting of $W\times W$ random Gaussian blocks (parametrized by $j,k \in\Lambda=[1,n]\cap \mathbb{Z}$, $N=nW$) with a fixed entry's variance $J_{jk}=W^{-1}(\delta_{j,k}+\beta\Delta_{j,k})$ in each block. Considering the limit $W, n\to\infty$, we prove that the behaviour of the second correlation function of such matrices in the bulk of the spectrum, as $W\gg \sqrt{N}$, is determined by the Wigner -- Dyson statistics. The method of the proof is based on the rigorous application of supersymmetric transfer matrix approach developed in [Shcherbina, M., Shcherbina, T.:Universality for 1d random band matrices: sigma-model approximation, J.Stat.Phys. 172, p. 627 -- 664 (2018)]


title: 20Ne + 76Ge elastic and inelastic scattering at 306 MeV

abstract: Background: Double charge exchange (DCE) nuclear reactions have recently attracted much interest as tools to provide experimentally driven information about nuclear matrix elements of interest in the context of neutrinoless double-beta decay. In this framework, a good description of the reaction mechanism and a complete knowledge of the initial and final-state interactions are mandatory. Presently, not enough is known about the details of the optical potentials and nuclear response to isospin operators for many of the projectile-target systems proposed for future DCE studies. Among these, the 20Ne + 76Ge DCE reaction is particularly relevant due to its connection with 76Ge double-beta decay. Purpose: We intend to characterize the initial-state interaction for the 20Ne + 76Ge reactions at 306 MeV bombarding energy and determine the optical potential and the role of the couplings between elastic channel and inelastic transitions to the first low-lying excited states. Methods: We determine the experimental elastic and inelastic scattering cross-section angular distributions, compare the theoretical predictions by adopting different models of optical potentials with the experimental data, and evaluate the coupling effect through the comparison of the distorted-wave Born approximation calculations with the coupled channels ones. Results: Optical models fail to describe the elastic angular distribution above the grazing angle (9.4{\deg}). A correction in the geometry to effectively account for deformation of the involved nuclear systems improves the agreement up to about 14{\deg}. Coupled channels effects are crucial to obtain good agreement at large angles in the elastic scattering cross section.


title: Trading off Complexity for Expressiveness in Programming Languages: Visions and Preliminary Experiences

abstract: When programming resource-scarce embedded smart devices, the designer often requires both the low-level system programming features of a language such as C and higher level capability typical of a language like Java. The choice of a particular language typically implies trade offs between conflicting design goals such as performance, costs, and overheads. The large variety of languages, virtual machines, and translators provides the designer with a dense trade off space, ranging from minimalistic to rich full-fledged approaches, but once a choice is made it is often difficult for the designer to revise it. In this work we propose a system of light-weighted and modular extensions as a method to flexibly reshape the target programming language as needed, adding only those application layer features that match the current design goals. In so doing complexity is made transparent, but not hidden: While the programmer can benefit of higher level constructs, the designer can deal with modular building blocks each characterized by a certain algorithmic complexity and therefore each accountable for a given share of the overhead. As a result the designer is given a finer control on the amount of resources that are consumed by the run-time executive of the chosen programming language.