title: Carpet-2 search for PeV gamma rays associated with IceCube high-energy neutrino events

abstract: Carpet-2 is an air-shower array at Baksan Valley, Russia, equipped with a large-area (175 m^2) muon detector, which makes it possible to separate primary photons from hadrons. We report the first results of the search for primary photons with energies E_\gamma>1 PeV, directionally associated with IceCube high-energy neutrino events, in the data obtained in 3080 days of Carpet-2 live time.


title: Search for astrophysical PeV gamma rays from point sources with Carpet-2

abstract: Early results of the search for E_\gamma>1 PeV cosmic photons from point sources with the data of Carpet-2, an air-shower array equipped with a 175 m^2 muon detector, are presented. They include 95% CL upper limits on PeV photon fluxes from stacked directions of high-energy IceCube neutrino events and from four predefined sources, Crab, Cyg X-3, Mrk 421 and Mrk 501. An insignificant excess of events from Mrk 421 will be further monitored. Prospects of the use of the upgraded installation, Carpet-3 (410 m^2 muon detector), scheduled to start data taking in 2019, for searches of E_\gamma>100 TeV photons, are briefly discussed.


title: Recursive Visual Attention in Visual Dialog

abstract: Visual dialog is a challenging vision-language task, which requires the agent to answer multi-round questions about an image. It typically needs to address two major problems: (1) How to answer visually-grounded questions, which is the core challenge in visual question answering (VQA); (2) How to infer the co-reference between questions and the dialog history. An example of visual co-reference is: pronouns (\eg, ``they'') in the question (\eg, ``Are they on or off?'') are linked with nouns (\eg, ``lamps'') appearing in the dialog history (\eg, ``How many lamps are there?'') and the object grounded in the image. In this work, to resolve the visual co-reference for visual dialog, we propose a novel attention mechanism called Recursive Visual Attention (RvA). Specifically, our dialog agent browses the dialog history until the agent has sufficient confidence in the visual co-reference resolution, and refines the visual attention recursively. The quantitative and qualitative experimental results on the large-scale VisDial v0.9 and v1.0 datasets demonstrate that the proposed RvA not only outperforms the state-of-the-art methods, but also achieves reasonable recursion and interpretable attention maps without additional annotations. The code is available at \url{https://github.com/yuleiniu/rva}.


title: Finite Element Method for Solving the Collective Nuclear Model with Tetrahedral Symmetry

abstract: We apply a new calculation scheme of a finite element method (FEM) for solving an elliptic boundary-value problem describing a quadrupole vibration collective nuclear model with tetrahedral symmetry. We use of shape functions constructed with interpolation Lagrange polynomials on a triangle finite element grid and compare the FEM results with obtained early by a finite difference method.


title: On the replica structure of Sachdev-Ye-Kitaev model

abstract: We investigate existence of replica off-diagonal solutions in the field-theoretical description of Sachdev-Ye-Kitaev model. To this end we evaluate a set of local and non-local dynamic correlation functions in the long time limit. We argue that the structure of the soft-mode Schwarzian action is qualitatively different in replica-diagonal vs. replica-off-diagonal scenarios, leading to distinct long-time predictions for the correlation functions. We then evaluate the corresponding correlation functions numerically and compare the simulations with analytical predictions of replica-diagonal and replica-off-diagonal calculations. We conclude that all our numerical results are in a quantitative agreement with the theory based on the replica-diagonal saddle point plus Schwarzian and massive Gaussian fluctuations (the latter do contain replica off-diagonal components). This seems to exclude any contributions from replica-off-diagonal saddle points, at least on the time scales shorter than the inverse many-body level spacing.


title: Hamiltonian analysis of mimetic scalar gravity revisited

abstract: We perform the Hamiltonian analysis of several mimetic gravity models and compare our results with those obtained previously by different authors. We verify that for healthy mimetic scalar-tensor theories the condition for the corresponding part of the Hamiltonian to be bounded from below is the positive value of the mimetic field energy density $\lambda$. We show that for mimetic dark matter possessing a shift symmetry the mimetic energy density remains positive in time, provided appropriate boundary conditions are imposed on its initial value, while in models without shift symmetry the positive energy density can be maintained by simply replacing $\lambda\to e^{\lambda}$. The same result also applies to mimetic $f(R)$ gravity, which is healthy if the usual stability conditions of the standard $f(R)$ gravity are assumed and $\lambda>0$. In contrast, if we add mimetic matter to an unhealthy seed action, the resulting mimetic gravity theory remains, in general, unstable. As an example, we consider a scalar-tensor theory with the higher-derivative term $(\Box \varphi)^2$, which contains an Ostrogradski ghost. We also revisit results regarding stability issues of linear perturbations around the FLRW background of the mimetic dark matter in the presence of ordinary scalar matter. We find that the presence of conventional matter does not revive dynamical ghost modes (at least in the UV limit). The modes, whose Hamiltonian is not positive definite, are non-propagating (have zero sound speed) and are associated with the mimetic matter itself. They are already present in the case in which the ordinary scalar fluid is absent, causing a growth of dust overdensity.


title: Modelling periodic structure formation on 100Cr6 steel after irradiation with femtosecond-pulsed laser beams

abstract: We investigate the periodic structure formation upon intense femtosecond pulsed irradiation of chrome steel (100Cr6) for linearly polarised laser beams. The underlying physical mechanism of the laser induced periodic structures is explored, their spatial frequency is calculated and theoretical results are compared with experimental observations. The proposed theoretical model comprises estimations of electron excitation, heat transfer, relaxation processes, and hydrodynamics-related mass transport. Simulations describe the sequential formation of sub-wavelength ripples and supra-wavelength grooves. In addition, the influence of the laser wavelength on the periodicity of the structures is discussed. The proposed theoretical investigation offers a systematic methodology towards laser processing of steel surfaces with important applications.


title: TASI Lectures on the Strong CP Problem and Axions

abstract: These TASI lectures where given in 2018 and provide an introduction to the Strong CP problem and the axion. I start by introducing the Strong CP problem from both a classical and a quantum mechanical perspective, calculating the neutron eDM and discussing the $\theta$ vacua. Next, I review the various solutions and discuss the active areas of axion model building. Finally, I summarize various experiments proposed to look for these solutions.


title: Associating Host Galaxy Candidates to Massive Black Hole Binaries resolved by Pulsar Timing Arrays

abstract: We propose a novel methodology to select host galaxy candidates of future pulsar timing array (PTA) detections of resolved gravitational waves (GWs) from massive black hole binaries (MBHBs). The method exploits the physical dependence of the GW amplitude on the MBHB chirp mass and distance to the observer, together with empirical MBH mass-host galaxy correlations, to rank potential host galaxies in the mass-redshift plane. This is coupled to a null-stream based likelihood evaluation of the GW amplitude and sky position in a Bayesian framework that assigns to each galaxy a probability of hosting the MBHB generating the GW signal. We test our algorithm on a set of realistic simulations coupling the likely properties of the first PTA resolved GW signal to synthetic all-sky galaxy maps. For a foreseeable PTA sky-localization precision of 100 squared degrees, we find that the GW source is hosted with 50% (90%) probability within a restricted number of <50 (<500) potential hosts. These figures are orders of magnitude smaller than the total number of galaxies within the PTA sky error-box, enabling extensive electromagnetic follow-up campaigns on a limited number of targets.


title: Spectral multipliers and wave equation for sub-Laplacians: lower regularity bounds of Euclidean type

abstract: Let $\mathscr{L}$ be a smooth second-order real differential operator in divergence form on a manifold of dimension $n$. Under a bracket-generating condition, we show that the ranges of validity of spectral multiplier estimates of Mihlin--H\"ormander type and wave propagator estimates of Miyachi--Peral type for $\mathscr{L}$ cannot be wider than the corresponding ranges for the Laplace operator on $\mathbb{R}^n$. The result applies to all sub-Laplacians on Carnot groups and more general sub-Riemannian manifolds, without restrictions on the step. The proof hinges on a Fourier integral representation for the wave propagator associated with $\mathscr{L}$ and nondegeneracy properties of the sub-Riemannian geodesic flow.


title: Probing New Physics in $B\to\pi K$ Decays

abstract: Over the last two decades, $B \to \pi K$ modes have received a lot of attention. They are particularly interesting since the main contribution to these decays comes from QCD penguin topologies. Furthermore, electroweak penguin amplitudes enter at a level comparable to the tree topologies. In the past, a discrepancy was found in the correlation between the CP asymmetries of $B_d^0 \to \pi^0 K_{\rm S}$. We provide the state-of-the-art picture of this situation and consider new constraints, finding that the effect has become larger. An attractive explanation is offered by a modified electroweak penguin sector. We propose a new method to determine the relevant parameters. It employs an isospin relation between the amplitudes of the $B \to \pi K$ decays, and makes only minimal use of the $SU(3)$ flavour symmetry. Using current data as well as future scenarios, we demonstrate how the mixing-induced CP asymmetry of $B_d^0 \to \pi^0 K_{\rm S}$ plays a key role in this analysis. The application of our strategy at the next generation of $B$-physics experiments may establish New Physics and reveal new sources of CP violation in the electroweak penguin sector.


title: Formation of periodic surface structures on dielectrics after irradiation with laser beams of spatially variant polarisation: a comparative study

abstract: A comparative study is performed to explore the periodic structure formation upon intense femtosecond pulsed irradiation of dielectrics with radially and azimuthally polarised beams. Laser conditions have been selected appropriately to produce excited carriers with densities below the optical breakdown threshold in order to highlight the role of phase transitions in surface modification mechanisms. The frequency of the laser induced structures is calculated based on a theoretical model that comprises estimation of electron density excitation, heat transfer, relaxation processes, and hydrodynamics-related mass transport. The influence of the laser wavelength in the periodicity of the structures is also unveiled. The decreased energy absorption for azimuthally polarised beams yields periodic structures with smaller frequencies which are more pronounced as the number of laser pulses applied to the irradiation spot increases. Similar results are obtained for laser pulses of larger photon energy and higher fluences. All induced periodic structures are oriented parallel to the laser beam polarisation.


title: MAUS: The MICE Analysis User Software

abstract: The Muon Ionization Cooling Experiment (MICE) collaboration has developed the MICE Analysis User Software (MAUS) to simulate and analyze experimental data. It serves as the primary codebase for the experiment, providing for offline batch simulation and reconstruction as well as online data quality checks. The software provides both traditional particle-physics functionalities such as track reconstruction and particle identification, and accelerator physics functions, such as calculating transfer matrices and emittances. The code design is object orientated, but has a top-level structure based on the Map-Reduce model. This allows for parallelization to support live data reconstruction during data-taking operations. MAUS allows users to develop in either Python or C++ and provides APIs for both. Various software engineering practices from industry are also used to ensure correct and maintainable code, including style, unit and integration tests, continuous integration and load testing, code reviews, and distributed version control. The software framework and the simulation and reconstruction capabilities are described.


title: Ghost and tachyon free Poincar\'e gauge theories: a systematic approach

abstract: A systematic method is presented for determining the conditions on the parameters in the action of a parity-preserving gauge theory of gravity for it to contain no ghost or tachyon particles. The technique naturally accommodates critical cases in which the parameter values lead to additional gauge invariances. The method is implemented as a computer program, and is used here to investigate the particle content of parity-conserving Poincar\'e gauge theory, which we compare with previous results in the literature. We find 450 critical cases that are free of ghosts and tachyons, and we further identify 10 of these that are also power-counting renormalizable, of which four have only massless tordion propagating particles and the remaining six have only a massive tordion propagating mode.


title: Theoretical Guarantees of Deep Embedding Losses Under Label Noise

abstract: Collecting labeled data to train deep neural networks is costly and even impractical for many tasks. Thus, research effort has been focused in automatically curated datasets or unsupervised and weakly supervised learning. The common problem in these directions is learning with unreliable label information. In this paper, we address the tolerance of deep embedding learning losses against label noise, i.e. when the observed labels are different from the true labels. Specifically, we provide the sufficient conditions to achieve theoretical guarantees for the 2 common loss functions: marginal loss and triplet loss. From these theoretical results, we can estimate how sampling strategies and initialization can affect the level of resistance against label noise. The analysis also helps providing more effective guidelines in unsupervised and weakly supervised deep embedding learning.


title: Phenomenological Extension for Tidal Charge Black Hole

abstract: A simple phenomenological extension of the black hole solution with tidal charge is proposed. Empirical data on the Sgr A* is consistent with the suggested metric which serves as a generalisation of the Reissner-Nordstrom one. Such a generalisation includes the leading effects beyond general relativity so, the discussed metric can explain wider range of gravitational effects. We discuss physical features of an object described by the proposed metric, namely, the size of its shadow and the innermost stable circular orbit radius.


title: Exploring the legacy of big stargazing events

abstract: We assess the impact of annual large-scale stargazing events, following the success of Stargazing Live and based on our experiences running events in Portsmouth, UK


title: Probing New Physics by the Tail of the Off-shell Higgs in $V_LV_L$ Mode

abstract: Off-shell Higgs at the high mass tail may shed light on the underlying mechanism of the electroweak symmetry breaking. Due to the large cancellation in the standard model (SM) between the box and Higgs-mediated triangle diagrams, the $gg\to WW(ZZ)$ process in the SM is dominated by the $V_T V_T$ transverse-mode at the high mass tail. The cancellation does not necessarily hold, when there is a sufficiently large new physics contribution resulting in $V_LV_L$ longitudinal mode, which is commonly the case when the Higgs sector is modified. Thus the $V_LV_L$ final states in the high mass tail can be utilized as a sensitive probe for such models. In the paper we focus on a study of the $gg \to ZZ$ process in the fully leptonic decay modes, proposing to utilize the polarization modes of the off-shell Higgs to probe new physics, whose contribution mainly shows in the longitudinal mode. As examples, we analyze three different Higgs sector new physics cases (Higgs portal with a light scalar, a broad-width scalar that mixes with the Higgs, and quantum critical Higgs models), and demonstrate that the angular information relating to the polarization serves as very sensitive probe for such new physics.


title: The structure of normal Hausdorff operators on Lebesgue spaces

abstract: We consider a generalization of Hausdorff operator and introduce the notion of the symbol of such an operator. Using this notion we describe the structure and investigate important properties (such as invertibility, spectrum, norm, and compactness) of normal generalized Hausdorff operators on Lebesgue spaces over $\mathbb{R}^n.$ The examples of Ces\`{a}ro operators are considered.


title: Harmonic superspace approach to the effective action in six-dimensional supersymmetric gauge theories

abstract: We review the recent progress in studying the quantum structure of $6D$, ${\cal N}=(1,0)$ and ${\cal N}=(1,1)$ supersymmetric gauge theories formulated through unconstrained harmonic superfields. The harmonic superfield approach allows one to carry out the quantization and calculations of the quantum corrections in a manifestly ${\cal N}=(1,0)$ supersymmetric way. The quantum effective action is constructed with the help of the background field method that secures the manifest gauge invariance of the results. Although the theories under consideration are not renormalizable, the extended supersymmetry essentially improves the ultraviolet behavior of the lowest-order loops. The ${\cal N}=(1,1)$ supersymmetric Yang--Mills theory turns out to be finite in the one-loop approximation in the minimal gauge. Also some two-loop divergences are shown to be absent in this theory. Analysis of the divergences is performed both in terms of harmonic supergraphs and by the manifestly gauge covariant superfield proper-time method. The finite one-loop leading low-energy effective action is calculated and analyzed. Also in the abelian case we discuss the gauge dependence of the quantum corrections and present its precise form for the one-loop divergent part of the effective action.


title: $\beta$-VAEs can retain label information even at high compression

abstract: In this paper, we investigate the degree to which the encoding of a $\beta$-VAE captures label information across multiple architectures on Binary Static MNIST and Omniglot. Even though they are trained in a completely unsupervised manner, we demonstrate that a $\beta$-VAE can retain a large amount of label information, even when asked to learn a highly compressed representation.


title: Optimized heat transfer at exceptional points in quantum circuits

abstract: Superconducting quantum circuits are potential candidates to realize a large-scale quantum computer. The envisioned large density of integrated components, however, requires a proper thermal management and control of dissipation. To this end, it is advantageous to utilize tunable dissipation channels and to exploit the optimized heat flow at exceptional points (EPs). Here, we experimentally realize an EP in a superconducting microwave circuit consisting of two resonators. The EP is a singularity point of the Hamiltonian, and corresponds to the most efficient heat transfer between the resonators without oscillation of energy. We observe a crossover from underdamped to overdamped coupling via the EP by utilizing photon-assisted tunneling as an \emph{in situ} tunable dissipative element in one of the resonators. The methods studied here can be applied to different circuits to obtain fast dissipation, for example, for initializing qubits to their ground states. In addition, these results pave the way towards thorough investigation of parity--time ($\mathcal{PT}$) symmetric systems and the spontaneous symmetry breaking in superconducting microwave circuits operating at the level of single energy quanta.


title: Range-separated tensor representation of the discretized multidimensional Dirac delta and elliptic operator inverse

abstract: In this paper, we introduce the operator dependent range-separated tensor approximation of the discretized Dirac delta in $\mathbb{R}^d$. It is constructed by application of the discrete elliptic operator to the range-separated decomposition of the associated Green kernel discretized on the Cartesian grid in $\mathbb{R}^d$. The presented operator dependent local-global splitting of the Dirac delta can be applied for solving the potential equations in non-homogeneous media when the density in the right-hand side is given by the large sum of pointwise singular charges. We show how the idea of the operator dependent RS splitting of the Dirac delta can be extended to the closely related problem on the range separated tensor representation of the elliptic resolvent. The numerical tests confirm the expected localization properties of the obtained operator dependent approximation of the Dirac delta represented on a tensor grid. As an example of application, we consider the regularization scheme for solving the Poisson-Boltzmann equation for modeling the electrostatics in bio-molecules.


title: Graph similarity drives zeolite diffusionless transformations and intergrowth

abstract: Predicting and directing polymorphic transformations is a critical challenge in zeolite synthesis. Although interzeolite transformations enable selective crystallization, their design lacks predictions to connect framework similarity and experimental observations. Here, computational and theoretical tools are combined to data-mine, analyze and explain interzeolite relations. It is observed that building units are weak predictors of topology interconversion and insufficient to explain intergrowth. By introducing a supercell-invariant metric that compares crystal structures using graph theory, we show that topotactic and reconstructive (diffusionless) transformations occur only between graph-similar pairs. Furthermore, all known instances of intergrowth occur between either structurally-similar or graph-similar frameworks. Backed with exhaustive literature results, we identify promising pairs for realizing novel diffusionless transformations and intergrowth. Hundreds of low-distance pairs are identified among known zeolites, and thousands of hypothetical frameworks are connected to known zeolites counterparts. The theory opens a venue to understand and control zeolite polymorphism.


title: Innermost stable circular orbits of neutron stars in dilatonic-Einstein-Gauss-Bonnet theory

abstract: The innermost stable circular orbits (ISCOs) around rapidly rotating neutron stars are studied in dilatonic Einstein-Gauss-Bonnet theory. Universal relations for properly scaled ISCO properties are extended from General Relativity to dilatonic Einstein-Gauss-Bonnet theory and additional relations are obtained.


title: Adaptive multicenter designs for continuous response clinical trials in the presence of an unknown sensitive subgroup

abstract: The partial effectiveness of drugs is of importance to the pharmaceutical industry. Randomized controlled trials (RCTs) assuming the existence of a subgroup sensitive to the treatment are already used. These designs, however, are available only if there is a known marker for identifying subjects in the subgroup. In this paper we investigate a model in which the response in the treatment group $Z^T$ has a two-component mixture density $(1-p)\mathcal N(\mu^C, \sigma^2)+p\mathcal N(\mu^T, \sigma^2)$ representing the treatment responses of \emph{placebo responders} and \emph{drug responders}. The treatment-specific effect is $\mu = \frac{\mu^T-\mu^C}{\sigma}$ and $p$ is the prevalence of the drug responders in the population. Other patients in the treatment group react as if they had received a placebo. We develop one- and two-stage RCT designs that are able to detect a sensitive subgroup based solely on the responses. We also extend them to a multicenter RCTs using Hochberg's step-up procedure. We avoid extensive simulations and use simple and quick numerical optimization methods.


title: Existence of stationary fronts in a system of two coupled wave equations with spatial inhomogeneity

abstract: We investigate the existence of stationary fronts in a coupled system of two sine-Gordon equations with a smooth, "hat-like" spatial inhomogeneity. The spatial inhomogeneity corresponds to a spatially dependent scaling of the sine-Gordon potential term. The uncoupled inhomogeneous sine-Gordon equation has stable stationary front solutions that persist in the coupled system. Carrying out a numerical investigation it is found that these inhomogeneous sine-Gordon fronts loose stability, provided the coupling between the two inhomogeneous sine-Gordon equations is strong enough, with new stable fronts bifurcating. In order to analytically study the bifurcating fronts, we first approximate the smooth spatial inhomogeneity by a piecewise constant function. With this approximation, we prove analytically the existence of a pitchfork bifurcation. To complete the argument, we prove that transverse fronts for a piecewise constant inhomogeneity persist for the smooth "hat-like" spatial inhomogeneity by introducing a fast-slow structure and using geometric singular perturbation theory.


title: Existence, uniqueness and coalescence of directed planar geodesics: proof via the increment-stationary growth process

abstract: We present a proof of the almost sure existence, uniqueness and coalescence of directed semi-infinite geodesics in planar growth models that is based on properties of an increment-stationary version of the growth process. The argument is developed in the context of the exponential corner growth model. It uses coupling, planar monotonicity, and properties of the stationary growth process to derive the existence of Busemann functions, which in turn control geodesics. This soft approach is in some situations an alternative to the much-applied 20-year-old arguments of C. Newman and co-authors. Along the way we derive some related results such as the distributional equality of the directed geodesic tree and its dual, originally due to L. Pimentel.


title: Provably Efficient Maximum Entropy Exploration

abstract: Suppose an agent is in a (possibly unknown) Markov Decision Process in the absence of a reward signal, what might we hope that an agent can efficiently learn to do? This work studies a broad class of objectives that are defined solely as functions of the state-visitation frequencies that are induced by how the agent behaves. For example, one natural, intrinsically defined, objective problem is for the agent to learn a policy which induces a distribution over state space that is as uniform as possible, which can be measured in an entropic sense. We provide an efficient algorithm to optimize such such intrinsically defined objectives, when given access to a black box planning oracle (which is robust to function approximation). Furthermore, when restricted to the tabular setting where we have sample based access to the MDP, our proposed algorithm is provably efficient, both in terms of its sample and computational complexities. Key to our algorithmic methodology is utilizing the conditional gradient method (a.k.a. the Frank-Wolfe algorithm) which utilizes an approximate MDP solver.


title: Baryon-antiBaryon Charge Asymmetry at LHC: String Junction Transfer and SJ Torus as DM Candidate

abstract: The asymmetry of baryon/antibaryon production has been measured in many proton-proton, pion-proton and electron-proton experiments. In the framework of Quark-Gluon String Model (QGSM) the energy dependence of asymmetry tells us about the value of $\alpha_{SJ}$(0), the intercept of String Junction Regge trajectory. In previous QGSM study, the value of intercept has been estimated as 0.5 < $\alpha_{SJ}$(0) < 0.9. Here, SJ behaviors are accumulated in the model based on topological expansion in order to build a neutral object with zero baryon charge. By the way, QCD mass falling under the event horizon of Black Hole (BH) should be symmetric, or in other words, have no charge information. The baryon junctions are easily combinable with antibaryon ones in hexagons. Topologically, hexagon net can coherently cover only the torus surface. The net on the torus has discrete number of baryon/antibaryon junctions. This is only parameter that marks the mass/energy level of this object. It looks like DM particle, is not it? In high energy collisions at LHC, such pomeron loops are to be produced approximately in 1.2 percent of inelastic events. Furthermore, the torus configurations of matter have been revealed in many bright events in space. As an example, Chandra experiment has detected such dense "doughnut" near the event horizon of Super Massive Black Hole (SMBH), which X-ray radiation is screened on 40 percent's. This topological symmetry model of DM seems rather realistic and can help us to deal with an "arm wrestling" between the stiffness of toroid structure of QCD matter and the pressure of gravitational singularity at extremely heavy masses. On the other hand, the instabilities in structure of matter in SMBHs can cause the bursts of giant relativistic hadron jets with the masses of order the own BH mass.


title: Termination of Solar Cycles and Correlated Tropospheric Variability

abstract: The Sun provides the energy required to sustain life on Earth and drive our planet's atmospheric circulation. However, establishing a solid physical connection between solar and tropospheric variability has posed a considerable challenge across the spectrum of Earth-system science. The canon of solar variability, the solar fiducial clock, lies almost exclusively with the 400 years of human telescopic observations that demonstrates the waxing and waning number of sunspots, over an 11(ish) year period. Recent research has demonstrated the critical importance of the underlying 22-year magnetic polarity cycle in establishing the shorter sunspot cycle. Integral to the manifestation of the latter is the spatio-temporal overlapping and migration of oppositely polarized magnetic bands. The points when these bands emerge at high solar latitudes and cancel at the equator are separated by almost 20 years. Here we demonstrate the impact of these "termination" points on the Sun's radiative output and particulate shielding of our atmosphere through the dramatically rapid reconfiguration of solar magnetism. These events reset the Sun's fiducial clock and present a new portal to explore the Sun-Earth connection. Using direct observation and proxies of solar activity going back six decades we can, with high statistical significance, demonstrate an apparent correlation between the solar cycle terminations and the largest swings of Earth's oceanic indices---a previously overlooked correspondence. Forecasting the Sun's global behavior places the next solar termination in early 2020; should a major oceanic swing follow, our challenge becomes: when does correlation become causation and how does the process work?


title: Quantifying error and leakage in an encoded Si/SiGe triple-dot qubit

abstract: Quantum computation requires qubits that satisfy often-conflicting criteria, including scalable control and long-lasting coherence. One approach to creating a suitable qubit is to operate in an encoded subspace of several physical qubits. Though such encoded qubits may be particularly susceptible to leakage out of their computational subspace, they can be insensitive to certain noise processes and can also allow logical control with a single type of entangling interaction while maintaining favorable features of the underlying physical system. Here we demonstrate a qubit encoded in a subsystem of three coupled electron spins confined in gated, isotopically enhanced silicon quantum dots. Using a modified "blind" randomized benchmarking protocol that determines both computational and leakage errors, we show that unitary operations have an average total error of 0.35%, with 0.17% of that coming from leakage driven by interactions with substrate nuclear spins. This demonstration utilizes only the voltage-controlled exchange interaction for qubit manipulation and highlights the operational benefits of encoded subsystems, heralding the realization of high-quality encoded multi-qubit operations.


title: Characters and invariant random subgroups of the finitary symmetric group

abstract: We will describe the relationship between the indecomposable characters of the finitary symmetric group and its ergodic invariant random subgroups; and we will interpret each Thoma character as an asymptotic limit of a naturally associated sequence of characters induced from linear characters of Young subgroups of finite symmetric groups.


title: Robust Modeling of Acoustic Phonon Transmission in Nanomechanical Structures

abstract: The transmission of acoustic phonons is an important element in the design and performance of nano-mechanical devices operating in the mesoscopic limit. Analytic expressions for the power transmission coefficient, T, exist only in the low-frequency (quasi-static) limit described by thin-plate elastic theory, and for well-defined elastic wave-guiding geometries. We compare two numerical techniques based on finite-element computations to determine the frequency dependence of T for arbitrary phonon scattering structures. Both methods take into account acoustic mode conversion to acoustic and optical modes. In one case, phase and amplitude of complex-valued reflected waves are determined and related to transmission through a Fresnel equation, while in the other the magnitude of the transmitted mechanical power is directly calculated. The numerical robustness of these methods is demonstrated by considering the transmission across an abrupt junction in a rectangular elastic beam, a well-known problem of considerable importance in mesoscopic device physics. The simulations presented extend the standard results for acoustic phonon transmission at an abrupt junction, and are in good agreement with analytic predictions from thin-plate elastic theory in the long-wavelength limit. More generally, the numerical methods developed provide an effective tool for calculating acoustic mode energy loss in nano-mechanical resonators through mode conversion and heat transfer in arbitrary mesoscopic structures.


title: Differentially Private Fair Learning

abstract: Motivated by settings in which predictive models may be required to be non-discriminatory with respect to certain attributes (such as race), but even collecting the sensitive attribute may be forbidden or restricted, we initiate the study of fair learning under the constraint of differential privacy. We design two learning algorithms that simultaneously promise differential privacy and equalized odds, a 'fairness' condition that corresponds to equalizing false positive and negative rates across protected groups. Our first algorithm is a private implementation of the equalized odds post-processing approach of [Hardt et al., 2016]. This algorithm is appealingly simple, but must be able to use protected group membership explicitly at test time, which can be viewed as a form of 'disparate treatment'. Our second algorithm is a differentially private version of the oracle-efficient in-processing approach of [Agarwal et al., 2018] that can be used to find the optimal fair classifier, given access to a subroutine that can solve the original (not necessarily fair) learning problem. This algorithm is more complex but need not have access to protected group membership at test time. We identify new tradeoffs between fairness, accuracy, and privacy that emerge only when requiring all three properties, and show that these tradeoffs can be milder if group membership may be used at test time. We conclude with a brief experimental evaluation.


title: On rotational-vibrational spectrum of diatomic beryllium molecule

abstract: The eigenvalue problem for second-order ordinary differential equation (SOODE) in a finite interval with the boundary conditions of the first, second and third kind is formulated. A computational scheme of the finite element method (FEM) is presented that allows the solution of the eigenvalue problem for a SOODE with the known potential function using the programs ODPEVP and KANTBP 4M that implement FEM in the Fortran and Maple, respectively. Numerical analysis of the solution using the KANTBP 4M program is performed for the SOODE exactly solvable eigenvalue problem. The discrete energy eigenvalues and eigenfunctions are analyzed for vibrational-rotational states of the diatomic beryllium molecule solving the eigenvalue problem for the SOODE numerically with the table-valued potential function approximated by interpolation Lagrange and Hermite polynomials and its asymptotic expansion for large values of the independent variable specified as Fortran function. The efficacy of the programs is demonstrated by the calculations of twelve eigenenergies of vibrational bound states with the required accuracy, in comparison with those known from literature, and the vibrational-rotational spectrum of the diatomic beryllium molecule.


title: A Numerical Routine for the Crossed Vertex Diagram with a Massive-Particle Loop

abstract: We present an evaluation of the two master integrals for the crossed vertex diagram with a closed loop of top quarks that allows for an easy numerical implementation. The differential equations obeyed by the master integrals are used to generate power series expansions centered around all the singular points. The different series are then matched numerically with high accuracy in intermediate points. The expansions allow a fast and precise numerical calculation of the two master integrals in all the regions of the phase space. A numerical routine that implements these expansions is presented.


title: Online Model Distillation for Efficient Video Inference

abstract: High-quality computer vision models typically address the problem of understanding the general distribution of real-world images. However, most cameras observe only a very small fraction of this distribution. This offers the possibility of achieving more efficient inference by specializing compact, low-cost models to the specific distribution of frames observed by a single camera. In this paper, we employ the technique of model distillation (supervising a low-cost student model using the output of a high-cost teacher) to specialize accurate, low-cost semantic segmentation models to a target video stream. Rather than learn a specialized student model on offline data from the video stream, we train the student in an online fashion on the live video, intermittently running the teacher to provide a target for learning. Online model distillation yields semantic segmentation models that closely approximate their Mask R-CNN teacher with 7 to 17$\times$ lower inference runtime cost (11 to 26$\times$ in FLOPs), even when the target video's distribution is non-stationary. Our method requires no offline pretraining on the target video stream, achieves higher accuracy and lower cost than solutions based on flow or video object segmentation, and can exhibit better temporal stability than the original teacher. We also provide a new video dataset for evaluating the efficiency of inference over long running video streams.


title: H\"ormander spaces on manifolds, and their application to elliptic boundary-value problems

abstract: We introduce an extended Sobolev scale on a smooth compact manifold with boundary. The scale is formed by inner-product H\"ormander spaces for which an RO-varying radial function serves as a regularity index. These spaces do not depend on a choice of local charts on the manifold. The scale consists of all Hilbert spaces that are interpolation ones for pairs of inner-product Sobolev spaces, is obtained by the interpolation with a function parameter of these pairs, and is closed with respect to this interpolation. As an application of the scale introduced, we give a theorem on the Fredholm property of a general elliptic boundary-value problem on appropriate H\"ormander spaces and find sufficient conditions under which its generalized solutions belong to the space of $p\geq0$ times continuously differential functions.


title: Kinetic theory of spin diffusion and superdiffusion in XXZ spin chains

abstract: We address the nature of spin transport in the integrable XXZ spin chain, focusing on the isotropic Heisenberg limit. We calculate the diffusion constant using a kinetic picture based on generalized hydrodynamics combined with Gaussian fluctuations: we find that it diverges, and show that a self-consistent treatment of this divergence gives superdiffusion, with an effective time-dependent diffusion constant that scales as $D(t) \sim t^{1/3}$. This exponent had previously been observed in large-scale numerical simulations, but had not been theoretically explained. We briefly discuss XXZ models with easy-axis anisotropy $\Delta > 1$. Our method gives closed-form expressions for the diffusion constant $D$ in the infinite-temperature limit for all $\Delta > 1$. We find that $D$ saturates at large anisotropy, and diverges as the Heisenberg limit is approached, as $D \sim (\Delta - 1)^{-1/2}$.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el quant-ph


title: Selective measurement of charge dynamics in an ensemble of nitrogen-vacancy centers in nano- and bulk diamonds

abstract: Nitrogen-vacancy (NV) centers in diamond have attracted considerable interest in sensing of weak magnetic fields, such as those created by biological systems. Detecting such feeble signals requires near-surface NV centers, to reduce the distance between NVs and sources. Moreover, dense ensembles of NVs are highly desirable to reduce measurement time. However, robust charge state switching is often observed in these systems, resulting in a complex interplay between charge and spin dynamics that can reduce the attainable level of spin polarization, and consequently, sensitivity. Understanding the mechanisms behind charge state switching is, therefore, crucial to developing NV based sensors. Here, we demonstrate a novel method to selectively measure charge dynamics in an ensemble of NVs by quenching the spin polarization using an off-axis magnetic field. Utilizing this technique, we show that, in nanodiamonds, charge state instability increases with increasing NV density. In the case of bulk single crystal diamond, we show that NV centers located near the surface are more stable in the neutral (NV$^{0}$) charge state, while the negatively charged (NV$^{-}$) form is more stable in bulk.


title: Higher-order Stein kernels for Gaussian approximation

abstract: We introduce higher-order Stein kernels relative to the standard Gaussian measure, which generalize the usual Stein kernels by involving higher-order derivatives of test functions. We relate the associated discrepancies to various metrics on the space of probability measures and prove new functional inequalities involving them. As an application, we obtain new explicit improved rates of convergence in the classical multidimensional CLT under higher moment and regularity assumptions.


title: A Radial Velocity Error Budget for Single-mode Fiber Doppler Spectrographs

abstract: Single-mode fiber (SMF) spectrographs fed with adaptive optics (AO) offer a unique path for achieving extremely precise radial velocity (EPRV) measurements. We present a radial velocity (RV) error budget based on end-to-end numerical simulations of an instrument named iLocater that is being developed for the Large Binocular Telescope(LBT). Representing the first AO-fed, SMF spectrograph, iLocater's design is used to quantify and assess the relative advantages and drawbacks of precise Doppler time series measurements made at the diffraction limit. This framework can be applied for trade-study work to investigate the impact of instrument design decisions on systematic uncertainties encountered in the regime of sub-meter-per-second precision. We find that working at the diffraction-limit through the use of AO and SMF's allows for high spectral resolution and improved instrument stability at the expense of limiting magnitude and longer integration times. Large telescopes equipped with AO alleviates the primary challenges of SMF spectrographs.


title: Frequency Tracking: LMS and RLS Applied to Speech Formant Estimation (2000)

abstract: Introduction Several speech processing algorithms assume the signal is stationary during short intervals (approximately 20 to 30 ms). This assumption is valid for several applications, but it is too restrictive in some contexts. This work investigates the application of adaptive signal processing to the problem of estimating the formant frequencies of speech. Two algorithms were implemented and tested. The first one is the conventional Least-Mean-Square (LMS) algorithm, and the second is the conventional Recursive Least-Squares (RLS) algorithm. The formant frequencies are the resonant frequencies of the vocal tract. The speech is the result of the convolution between the excitation and the vocal tract impulse response [Rabiner, 78], thus a kind of "deconvolution" is required to recover the formants. This is not an easy problem because one does not have the excitation signal available. There are several algorithms for formant estimation [Rabiner, 78], [Snell, 93], [Laprie, 94


title: Coarse-Graining Auto-Encoders for Molecular Dynamics

abstract: Molecular dynamics simulations provide theoretical insight into the microscopic behavior of materials in condensed phase and, as a predictive tool, enable computational design of new compounds. However, because of the large temporal and spatial scales involved in thermodynamic and kinetic phenomena in materials, atomistic simulations are often computationally unfeasible. Coarse-graining methods allow simulating larger systems, by reducing the dimensionality of the simulation, and propagating longer timesteps, by averaging out fast motions. Coarse-graining involves two coupled learning problems; defining the mapping from an all-atom to a reduced representation, and the parametrization of a Hamiltonian over coarse-grained coordinates. Multiple statistical mechanics approaches have addressed the latter, but the former is generally a hand-tuned process based on chemical intuition. Here we present Autograin, an optimization framework based on auto-encoders to learn both tasks simultaneously. Autograin is trained to learn the optimal mapping between all-atom and reduced representation, using the reconstruction loss to facilitate the learning of coarse-grained variables. In addition, a force-matching method is applied to variationally determine the coarse-grained potential energy function. This procedure is tested on a number of model systems including single-molecule and bulk-phase periodic simulations.


title: Video Action Transformer Network

abstract: We introduce the Action Transformer model for recognizing and localizing human actions in video clips. We repurpose a Transformer-style architecture to aggregate features from the spatiotemporal context around the person whose actions we are trying to classify. We show that by using high-resolution, person-specific, class-agnostic queries, the model spontaneously learns to track individual people and to pick up on semantic context from the actions of others. Additionally its attention mechanism learns to emphasize hands and faces, which are often crucial to discriminate an action - all without explicit supervision other than boxes and class labels. We train and test our Action Transformer network on the Atomic Visual Actions (AVA) dataset, outperforming the state-of-the-art by a significant margin using only raw RGB frames as input.


title: Ternary mixed-anion semiconductors with tunable band gaps from machine-learning and crystal structure prediction

abstract: We report the computational investigation of a series of ternary X$_4$Y$_2$Z and X$_5$Y$_2$Z$_2$ compounds with X={Mg, Ca, Sr, Ba}, Y={P, As, Sb, Bi}, and Z={S, Se, Te}. The compositions for these materials were predicted through a search guided by machine learning, while the structures were resolved using the minima hopping crystal structure prediction method. Based on $\textit{ab initio}$ calculations, we predict that many of these compounds are thermodynamically stable. In particular, 21 of the X$_4$Y$_2$Z compounds crystallize in a tetragonal structure with $\textit{I-42d}$ symmetry, and exhibit band gaps in the range of 0.3 and 1.8 eV, well suited for various energy applications. We show that several candidate compounds (in particular X$_4$Y$_2$Te and X$_4$Sb$_2$Se) exhibit good photo absorption in the visible range, while others (e.g., Ba$_4$Sb$_2$Se) show excellent thermoelectric performance due to a high power factor and extremely low lattice thermal conductivities.


title: On stochastic gradient Langevin dynamics with dependent data streams in the logconcave case

abstract: We study the problem of sampling from a probability distribution $\pi$ on $\rset^d$ which has a density \wrt\ the Lebesgue measure known up to a normalization factor $x \mapsto \rme^{-U(x)} / \int_{\rset^d} \rme^{-U(y)} \rmd y$. We analyze a sampling method based on the Euler discretization of the Langevin stochastic differential equations under the assumptions that the potential $U$ is continuously differentiable, $\nabla U$ is Lipschitz, and $U$ is strongly concave. We focus on the case where the gradient of the log-density cannot be directly computed but unbiased estimates of the gradient from possibly dependent observations are available. This setting can be seen as a combination of a stochastic approximation (here stochastic gradient) type algorithms with discretized Langevin dynamics. We obtain an upper bound of the Wasserstein-2 distance between the law of the iterates of this algorithm and the target distribution $\pi$ with constants depending explicitly on the Lipschitz and strong convexity constants of the potential and the dimension of the space. Finally, under weaker assumptions on $U$ and its gradient but in the presence of independent observations, we obtain analogous results in Wasserstein-2 distance.


title: On Critical Infrastructures, Their Security and Resilience - Trends and Vision

abstract: This short paper is presented in observance and promotion of November, the National Month of Critical Infrastructure Security and Resilience (CISR), established by the United States Department of Homeland Security in 2013. The CISR term focuses on essential assets (critical infrastructures) and two ultimate goals of making them secure and resilient. These assets and goals were put together in 2013 in the now well-known Presidential Policy Directive on CISR (PPD-21). This paper presents easy-to-ready material laying down the building blocks of CISR - what it means to you as a regular citizen, professional, or government worker. This paper presents concepts behind security and resilience pertinent to various types of activities - from every day to field-specific activities. This paper also presents basic elements to the field: 1. high-level introduction to the organizational units dealing with CISR in the United States; 2. explanation of basic terms and a list of further reading material; and 3. several discussion topics on the vision and future of CISR in critical infrastructure cyber-physical systems.


title: Formal Synthesis of Analytic Controllers for Sampled-Data Systems via Genetic Programming

abstract: This paper presents an automatic formal controller synthesis method for nonlinear sampled-data systems with safety and reachability specifications. Fundamentally, the presented method is not restricted to polynomial systems and controllers. We consider periodically switched controllers based on a Control Lyapunov Barrier-like functions. The proposed method utilizes genetic programming to synthesize these functions as well as the controller modes. Correctness of the controller are subsequently verified by means of a Satisfiability Modulo Theories solver. Effectiveness of the proposed methodology is demonstrated on multiple systems.


title: Analyzing the hate and counter speech accounts on Twitter

abstract: The online hate speech is proliferating with several organization and countries implementing laws to ban such harmful speech. While these restrictions might reduce the amount of such hateful content, it does so by restricting freedom of speech. Thus, an promising alternative supported by several organizations is to counter such hate speech with more speech. In this paper, We analyze hate speech and the corresponding counters (aka counterspeech) on Twitter. We perform several lexical, linguistic and psycholinguistic analysis on these user accounts and obverse that counter speakers employ several strategies depending on the target community. The hateful accounts express more negative sentiments and are more profane. We also find that the hate tweets by verified accounts have much more virality as compared to a tweet by a non-verified account. While the hate users seem to use words more about envy, hate, negative emotion, swearing terms, ugliness, the counter users use more words related to government, law, leader. We also build a supervised model for classifying the hateful and counterspeech accounts on Twitter and obtain an F-score of 0.77. We also make our dataset public to help advance the research on hate speech.


title: PartNet: A Large-scale Benchmark for Fine-grained and Hierarchical Part-level 3D Object Understanding

abstract: We present PartNet: a consistent, large-scale dataset of 3D objects annotated with fine-grained, instance-level, and hierarchical 3D part information. Our dataset consists of 573,585 part instances over 26,671 3D models covering 24 object categories. This dataset enables and serves as a catalyst for many tasks such as shape analysis, dynamic 3D scene modeling and simulation, affordance analysis, and others. Using our dataset, we establish three benchmarking tasks for evaluating 3D part recognition: fine-grained semantic segmentation, hierarchical semantic segmentation, and instance segmentation. We benchmark four state-of-the-art 3D deep learning algorithms for fine-grained semantic segmentation and three baseline methods for hierarchical semantic segmentation. We also propose a novel method for part instance segmentation and demonstrate its superior performance over existing methods.


title: Defining the Proton Radius: a Unified Treatment

abstract: Background: There is significant current interest in knowing the value of the proton radius and also its proper definition. Purpose: Combine the disparate literatures of hydrogen spectroscopy and diverse modern parton distributions to show that the quantity $r_p^2\equiv -6 G_E'(0)$ is the relativistically proper definition that originates from the separate bodies of work. Methods: Use perturbation theory, light-front dynamics and elementary techniques to find relativistically correct definitions of the proton radius and charge density. Results: It is found that the very same proton radius is accessed by measurements of hydrogen spectroscopy and elastic lepton scattering. The derivation of the mean-square radius as a moment of a spherically symmetric three-dimensional density is shown to be incorrect. A relativistically-correct, two-dimensional charge density is related to the diverse modern literature of various parton distributions. Relativistically invariant moments thereof are derived in a new moment expansion, the RME.


title: An Improved Cost Function for Hierarchical Cluster Trees

abstract: Hierarchical clustering has been a popular method in various data analysis applications. It partitions a data set into a hierarchical collection of clusters, and can provide a global view of (cluster) structure behind data across different granularity levels. A hierarchical clustering (HC) of a data set can be naturally represented by a tree, called a HC-tree, where leaves correspond to input data and subtrees rooted at internal nodes correspond to clusters. Many hierarchical clustering algorithms used in practice are developed in a procedure manner. Dasgupta proposed to study the hierarchical clustering problem from an optimization point of view, and introduced an intuitive cost function for similarity-based hierarchical clustering with nice properties as well as natural approximation algorithms. We observe that while Dasgupta's cost function is effective at differentiating a good HC-tree from a bad one for a fixed graph, the value of this cost function does not reflect how well an input similarity graph is consistent to a hierarchical structure. In this paper, we present a new cost function, which is developed based on Dasgupta's cost function, to address this issue. The optimal tree under the new cost function remains the same as the one under Dasgupta's cost function. However, the value of our cost function is more meaningful. The new way of formulating the cost function also leads to a polynomial time algorithm to compute the optimal cluster tree when the input graph has a perfect HC-structure, or an approximation algorithm when the input graph 'almost' has a perfect HC-structure. Finally, we provide further understanding of the new cost function by studying its behavior for random graphs sampled from an edge probability matrix.


title: Cross-Domain 3D Equivariant Image Embeddings

abstract: Spherical convolutional networks have been introduced recently as tools to learn powerful feature representations of 3D shapes. Spherical CNNs are equivariant to 3D rotations making them ideally suited to applications where 3D data may be observed in arbitrary orientations. In this paper we learn 2D image embeddings with a similar equivariant structure: embedding the image of a 3D object should commute with rotations of the object. We introduce a cross-domain embedding from 2D images into a spherical CNN latent space. This embedding encodes images with 3D shape properties and is equivariant to 3D rotations of the observed object. The model is supervised only by target embeddings obtained from a spherical CNN pretrained for 3D shape classification. We show that learning a rich embedding for images with appropriate geometric structure is sufficient for tackling varied applications, such as relative pose estimation and novel view synthesis, without requiring additional task-specific supervision.


title: PROSA PDFs and astrophysical applications

abstract: The PROSA parton distribution function fit was the first one appeared in the literature incorporating data on open charm and open bottom production at LHCb, in order to reduce the uncertainties on gluons and sea quarks at low $x$'s ($x$ $<$ $10^{-4}$). We discuss aspects of the PROSA PDFs of particular relevance for their usage in the field of Neutrino Astronomy, and their application in the computation of prompt neutrino fluxes.


title: Chandra Early-Type Galaxy Atlas

abstract: The hot ISM in early type galaxies (ETGs) plays a crucial role in understanding their formation and evolution. The structural features of the hot gas identified by Chandra observations point to key evolutionary mechanisms, (e.g., AGN and stellar feedback, merging history). In our Chandra Galaxy Atlas (CGA) project, taking full advantage of the Chandra capabilities, we systematically analyzed the archival Chandra data of 70 ETGs and produced uniform data products for the hot gas properties. The primary data products are spatially resolved 2D spectral maps of the hot gas from individual galaxies. We emphasize that new features can be identified in the spectral maps which are not readily visible in the surface brightness maps. The high-level images can be viewed at the dedicated CGA website, and the CGA data products can be downloaded to compare with data at other wavelengths and to perform further analyses. Utilizing our data products, we address a few focused science topics.


title: Light hidden photon production in high energy collisions

abstract: The visible and dark sectors of particle physics can be connected via kinetic mixing between ordinary and hidden photons. If the latter is light its production in high energy collisions of ordinary particles proceeds via oscillations with ordinary photons similarly to the neutrino processes. Generically, the experiments are insensitive to mass of the hidden vector, if it is lighter than 1\,MeV, and it does not decay into $e^+e^-$-pair. Still, one can use the missing energy and scattering off the detector material as signatures to search for the light vectors. Presence of media suppresses production of the light vectors making the experiments insensitive to the entire model. We present analytic formulas for the light hidden photon production, propagation and detection valid for searches at colliders and beam-target experiments and apply them to estimate the impact on the sensitivities of a set of experiments --- NA64, FASER, MATHUSLA, SHiP, T2K, DUNE, NA62 --- in a zero-background case.


title: Deconfined criticality in the QED$_{3}$ Gross-Neveu-Yukawa model: The $1/N$ expansion revisited

abstract: The critical properties of the QED$_{3}$ Gross-Neveu-Yukawa (GNY) model in 2+1 dimensions with $N$ flavors of two-component Dirac fermions are computed to first order in the $1/N$ expansion. For the specific case of $N=2$, the critical point is conjectured to be dual to the N\'eel-to-valence-bond-solid (VBS) deconfined critical point of quantum antiferromagnets on the square lattice. It is found that Aslamazov-Larkin diagrams, missed by previous $\epsilon$- and $1/N$-expansion studies with four-component fermions, give important contributions to the scaling dimensions of various operators. With the inclusion of these diagrams, the resummed scaling dimensions of the adjoint fermion bilinear and scalar field at the QED$_{3}$ GNY critical point are in reasonable agreement with numerical studies of the N\'eel-to-VBS transition, in support of the duality conjecture.


title: Classifying local fractal subsystem symmetry protected topological phases

abstract: We study symmetry-protected topological (SPT) phases of matter in 2D protected by symmetries acting on fractal subsystems of a certain type. Despite the total symmetry group of such systems being subextensively large, we show that only a small number of phases are actually realizable by local Hamiltonians. Which phases are possible depends crucially on the spatial structure of the symmetries, and we show that in many cases no non-trivial SPT phases are possible at all. In cases where non-trivial SPT phases do exist, we give an exhaustive enumeration of them in terms of their locality.


title: Project Rosetta: A Childhood Social, Emotional, and Behavioral Developmental Ontology

abstract: There is a wide array of existing instruments used to assess childhood behavior and development for the evaluation of social, emotional and behavioral disorders. Many of these instruments either focus on one diagnostic category or encompass a broad set of childhood behaviors. We built an extensive ontology of the questions associated with key features that have diagnostic relevance for child behavioral conditions, such as Autism Spectrum Disorder (ASD), attention-deficit/hyperactivity disorder (ADHD), and anxiety, by incorporating a subset of existing child behavioral instruments and categorizing each question into clinical domains. Each existing question and set of question responses were then mapped to a new unique Rosetta question and set of answer codes encompassing the semantic meaning and identified concept(s) of as many existing questions as possible. This resulted in 1274 existing instrument questions mapping to 209 Rosetta questions creating a minimal set of questions that are comprehensive of each topic and subtopic. This resulting ontology can be used to create more concise instruments across various ages and conditions, as well as create more robust overlapping datasets for both clinical and research use.


title: Possible Bilepton Resonances in Like-Sign Pairs

abstract: We consider pair production of bileptons $Y^{++}Y^{--}$ at the LHC for the presently accumulated integrated luminosity of $150/fm$. It is shown that the entire mass range 800 GeV $\leq$ M(Y) $\leq$ 2000 GeV can be successfully searched. A bilepton resonance will have an exceptionally large ratio of signal to background because the standard model prediction is so infinitesimal. A $5\sigma$ discovery is quite feasible.


title: Structural Damage Detection Using Ensemble Empirical Mode Decomposition, Hilbert Transform and Artificial Neural Networks

abstract: Civil structures are on the verge of changing which leads energy dissipation capacity to decline. Structural Health Monitoring (SHM) as a process in order to implement a damage detection strategy and assess the condition of structure plays a key role in structural reliability. Earthquake is a recognized factor in variation of structures condition, inasmuch as inelastic behavior of a building subjected to design level earthquakes is plausible. In this study Hilbert Huang Transformation (HHT) is superseded by Ensemble Empirical Mode decomposition (EEMD) and Hilbert Transform (HT) together. Albeit analogous, EEMD brings more appropriate Intrinsic Mode Functions (IMFs) than Empirical Mode Decomposition (EMD). IMFs are employed to assess first mode frequency and mode shape. Afterward, Artificial Neural Networks (ANN) is applied to predict story acceleration based on acceleration of structure during previous moments. ANN functions precisely. Therefore, any congruency between predicted and measured acceleration provides onset of damage. Then another ANN method is applied to estimate stiffness matrix. Though first mode shape and frequency is calculated in advance, it essentially requires an inverse problem to be solved in order to find stiffness matrix. This task is done by ANN. In other words, these two ANN methods are exercised to forecast location and measure severity of damage respectively. This algorithm is implemented on one nonlinear moment-resisting steel frame and the results are acceptable.


title: Visual Object Networks: Image Generation with Disentangled 3D Representation

abstract: Recent progress in deep generative models has led to tremendous breakthroughs in image generation. However, while existing models can synthesize photorealistic images, they lack an understanding of our underlying 3D world. We present a new generative model, Visual Object Networks (VON), synthesizing natural images of objects with a disentangled 3D representation. Inspired by classic graphics rendering pipelines, we unravel our image formation process into three conditionally independent factors---shape, viewpoint, and texture---and present an end-to-end adversarial learning framework that jointly models 3D shapes and 2D images. Our model first learns to synthesize 3D shapes that are indistinguishable from real shapes. It then renders the object's 2.5D sketches (i.e., silhouette and depth map) from its shape under a sampled viewpoint. Finally, it learns to add realistic texture to these 2.5D sketches to generate natural images. The VON not only generates images that are more realistic than state-of-the-art 2D image synthesis methods, but also enables many 3D operations such as changing the viewpoint of a generated image, editing of shape and texture, linear interpolation in texture and shape space, and transferring appearance across different objects and viewpoints.


title: Simulation of Stylized Facts in Agent-Based Computational Economic Market Models

abstract: We study the qualitative and quantitative appearance of stylized facts in several agent-based computational economic market (ABCEM) models. We perform our simulations with the SABCEMM (Simulator for Agent-Based Computational Economic Market Models) tool recently introduced by the authors (Trimborn et al. 2019). Furthermore, we present novel ABCEM models created by recombining existing models and study them with respect to stylized facts as well. This can be efficiently performed by the SABCEMM tool thanks to its object-oriented software design. The code is available on GitHub (Trimborn et al. 2018), such that all results can be reproduced by the reader.


title: Partial compositeness and baryon matrix elements on the lattice

abstract: Partial compositeness is a mechanism for the generation of fermion masses which replaces a direct Higgs coupling to the fermions by a linear mixing with heavy composite partners. We present the first calculation of the relevant matrix element in a lattice model which is very close to a candidate theory containing a composite Higgs boson and a partially composite top quark. Specifically, our model is an SU(4) gauge theory coupled to dynamical fermions in the fundamental and two-index antisymmetric (sextet) representations. The matrix element we obtain is small and hence our result disfavors the scenario of obtaining a realistic top mass in this model.


title: On Decoupling the Integrals of Cosmological Perturbation Theory

abstract: Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher-order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high-dimensional and naively require increasing computational resources the higher the order of the corrections. Here we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1-D integrals. Our approach is related to a commonly used method for calculating multi-loop Feynman integrals in Quantum Field Theory, the Gegenbauer Polynomial $x$-Space Technique (GPxT). We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9-D integrations, to sums over successive 1-D radial integrals. These 1-D integrals can further be performed as convolutions, rendering the scaling of this method $N_{\rm g} \log N_{\rm g}$ with $N_{\rm g}$ the number of grid points used for each Fast Fourier Transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fit values.


title: Geometric variational principles for computational homogenization

abstract: The homogenization of periodic elastic composites is addressed through the reformulation of the local equations of the mechanical problem in a geometric functional setting. This relies on the definition of Hilbert spaces of kinematically and statically admissible tensor fields, whose orthogonality and duality properties are recalled. These are endowed with specific energetic scalar products that make use of a reference and uniform elasticity tensor. The corresponding strain and stress Green's operators are introduced and interpreted as orthogonal projection operators in the admissibility spaces. In this context and as an alternative to classical minimum energy principles, two geometric variational principles are investigated with the introduction of functionals that aim at measuring the discrepancy of arbitrary test fields to the kinematic, static or material admissibility conditions of the problem. By relaxing the corresponding local equations, this study aims in particular at laying the groundwork for the homogenization of composites whose constitutive properties are only partially known or uncertain. The local fields in the composite and their macroscopic responses are computed through the minimization of the proposed geometric functionals. To do so, their gradients are computed using the Green's operators and gradient-based optimization schemes are discussed. A FFT-based implementation of these schemes is proposed and they are assessed numerically on a canonical example for which analytical solutions are available.


title: Joint cosmological inference of standard sirens and gravitational wave weak lensing

abstract: We present the first joint inference of standard sirens and gravitational wave weak lensing by filtering of the same dataset. We imagine a post-LISA scenario emerging around the late 2030s when LISA will have accumulated a number of detections at high redshift; LIGO-VIRGO will have finished observing at low redshift, and Einstein Telescope will have started making new observations out to redshifts possibly overlapping with LISA. Euclid and other cosmological probes will have provided constraints at percent level by then, but mostly exhausted their ability to improve any further. We derive forecasts assuming $\sim1\,\text{deg}^{-2}$ detected sources, in conjunction with a spectroscopic follow-up (e.g. by Euclid, DESI, or ATHENA). Thanks to the statistical power of standard sirens as a geometry probe -- lifting key degeneracies in the gravitational wave weak lensing -- and no external priors assumed, the constraints on dark matter and its clustering, namely $\Omega_m$ and $\sigma_8$, could be achieved to $2\%$ and $3\%$. The Hubble constant could be constrained to better than $1\%$ in all cases; the dark energy density, $\Omega_\Lambda$, to $2\%$, and curvature, $\Omega_K$, to $0.02$; the amplitude and spectral tilt of the scalar fluctuations, $\ln(10^{10}A_s)$ and $n_s$, to $2\%$ and $7\%$. As a completely independent cosmological probe, with less calibration requirements, the joint inference of standard sirens and gravitational wave weak lensing might help solve the tensions currently observed between other cosmological probes, such as CMB, galaxy lensing, and Type Ia SNs, and distinguish between residual systematics and new physics.


title: A new scale in the bias expansion

abstract: The fact that the spatial nonlocality of galaxy formation is controlled by some short length scale like the Lagrangian radius is the cornerstone of the bias expansion for large-scale-structure tracers. However, the first sources of ionizing radiation between $z\approx 15$ and $z\approx 6$ are expected to have significant effects on the formation of galaxies we observe at lower redshift, at least on low-mass galaxies. These radiative-transfer effects introduce a new scale in the clustering of galaxies, i.e. the finite distance which ionizing radiation travels until it reaches a given galaxy. This mean free path can be very large, of order $100\,h^{-1}\,{\rm Mpc}$. Consequently, higher-derivative terms in the bias expansion could turn out to be non-negligible even on these scales: treating them perturbatively would lead to a massive loss in predictivity and, for example, could spoil the determination of the BAO feature or constraints on the neutrino mass. Here, we investigate under what assumptions an explicit non-perturbative model of radiative-transfer effects can maintain the robustness of large-scale galaxy clustering as a cosmological probe.


title: The simultaneous formation of cored, tangentially biased, and kinematically decoupled centers in massive early-type galaxies

abstract: We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of massive ETGs, namely flat stellar density distributions (cores), tangentially biased velocity distributions and kinematically decoupled (counter-)rotation can naturally result from a single process $-$ the scouring by SMBHs. Major mergers with mass ratios of $q>1/3$ produce flat, tangentially biased cores with kinematically distinct components. These features originate from spin reversals during the decay of the SMBH orbits caused by gravitational torques from tidally expelled nuclear material. Minor mergers ($q\lesssim1/3$) on the other hand, form non-rotating cores and the tidal reversal becomes less important. Low-density stellar cores scoured in (multiple) minor mergers are less tangentially biased. This implies that the nuclear stellar properties of massive ETGs can be solely explained by stellar dynamical processes during their final assembly without any need for `feedback' from accreting black holes. We predict a strong correlation between decoupled cores, central anisotropy and merger history: decoupled cores form in binary mergers and we predict them to occur in elliptical galaxies with the strongest central anisotropy. Measurements of the central orbital structure are the key to understanding the number of mergers a given galaxy has experienced.


title: Spin transfer torque and exchange coupling in Josephson junctions with ferromagnetic superconductor reservoirs

abstract: In this paper, the spin transfer torque (STT) and the exchange coupling of the Josephson junctions containing interesting cases of diffusive/ballistic-triplet/singlet ferromagnetic superconductor materials together with diffusive Josephson junction of the form S$_{1}$/F$_{1}$/I$_1$/N/I$_2$/F$_2$ with I being insulating barrier are investigated. Using the Nazarov quantum circuit theory, it is found that for the diffusive FS$_1$/N/FS$_2$ structure the only emerged torque in normal direction to the plane of the exchange fields of F$_1$ and F$_2$ results in antiparallel/parallel or vice versa parallel/antiparallel transition of favorable exchange coupling depending on the considered parameters of the system, including the nonmagnetic spacer thickness, the superconducting phase difference, the length and the exchange field of the ferromagnets. Furthermore, the analyze of the width of the transitions, the phase difference interval in which an interlayer length-induced antiparallel/parallel transition can be occurred, is performed. For instance, as the exchange field or the temperature increases, the interval of phase difference gets larger. On the other hand, ballistic Josephson junction containing the triplet ferromagnetic superconductor reservoirs solving the 16$\times$16 Bogoliubov-de-Gennes equation is studied. It is found that although the exchange fields of the FS are laid in the z and y direction, the STT interestingly exists in all three directions of x, y and z. This exciting finding suggests that the favorable equilibrium configuration concerning the least exchange coupling occurs in the relative exchange field direction different from 0 or $\pi$. To the best of our knowledge it is for the first time that the occurrence of the in-plane STT is reported. Moreover, the occurrence of the beat like behavior with two oscillation period for the out-of-plane STT is interestingly acquired.


title: Learning to Design Circuits

abstract: Analog IC design relies on human experts to search for parameters that satisfy circuit specifications with their experience and intuitions, which is highly labor intensive, time consuming and suboptimal. Machine learning is a promising tool to automate this process. However, supervised learning is difficult for this task due to the low availability of training data: 1) Circuit simulation is slow, thus generating large-scale dataset is time-consuming; 2) Most circuit designs are propitiatory IPs within individual IC companies, making it expensive to collect large-scale datasets. We propose Learning to Design Circuits (L2DC) to leverage reinforcement learning that learns to efficiently generate new circuits data and to optimize circuits. We fix the schematic, and optimize the parameters of the transistors automatically by training an RL agent with no prior knowledge about optimizing circuits. After iteratively getting observations, generating a new set of transistor parameters, getting a reward, and adjusting the model, L2DC is able to optimize circuits. We evaluate L2DC on two transimpedance amplifiers. Trained for a day, our RL agent can achieve comparable or better performance than human experts trained for a quarter. It first learns to meet hard-constraints (eg. gain, bandwidth), and then learns to optimize good-to-have targets (eg. area, power). Compared with grid search-aided human design, L2DC can achieve $\mathbf{250}\boldsymbol{\times}$ higher sample efficiency with comparable performance. Under the same runtime constraint, the performance of L2DC is also better than Bayesian Optimization.


title: Discriminants of stable rank two sheaves on some general type surfaces

abstract: We prove sharp bounds on the discriminants of stable rank two sheaves on surfaces in three-dimensional projective space. The key technical ingredient is to study them as torsion sheaves in projective space via tilt stability in the derived category. We then proceed to describe the surface itself as a moduli space of rank two vector bundles on it. Lastly, we give a proof of the Bogomolov inequality for semistable rank two sheaves on integral surfaces in three-dimensional projective space in all characteristics.


title: A Technical Survey on Statistical Modelling and Design Methods for Crowdsourcing Quality Control

abstract: Online crowdsourcing provides a scalable and inexpensive means to collect knowledge (e.g. labels) about various types of data items (e.g. text, audio, video). However, it is also known to result in large variance in the quality of recorded responses which often cannot be directly used for training machine learning systems. To resolve this issue, a lot of work has been conducted to control the response quality such that low-quality responses cannot adversely affect the performance of the machine learning systems. Such work is referred to as the quality control for crowdsourcing. Past quality control research can be divided into two major branches: quality control mechanism design and statistical models. The first branch focuses on designing measures, thresholds, interfaces and workflows for payment, gamification, question assignment and other mechanisms that influence workers' behaviour. The second branch focuses on developing statistical models to perform effective aggregation of responses to infer correct responses. The two branches are connected as statistical models (i) provide parameter estimates to support the measure and threshold calculation, and (ii) encode modelling assumptions used to derive (theoretical) performance guarantees for the mechanisms. There are surveys regarding each branch but they lack technical details about the other branch. Our survey is the first to bridge the two branches by providing technical details on how they work together under frameworks that systematically unify crowdsourcing aspects modelled by both of them to determine the response quality. We are also the first to provide taxonomies of quality control papers based on the proposed frameworks. Finally, we specify the current limitations and the corresponding future directions for the quality control research.


title: On Configurable Defense against Adversarial Example Attacks

abstract: Machine learning systems based on deep neural networks (DNNs) have gained mainstream adoption in many applications. Recently, however, DNNs are shown to be vulnerable to adversarial example attacks with slight perturbations on the inputs. Existing defense mechanisms against such attacks try to improve the overall robustness of the system, but they do not differentiate different targeted attacks even though the corresponding impacts may vary significantly. To tackle this problem, we propose a novel configurable defense mechanism in this work, wherein we are able to flexibly tune the robustness of the system against different targeted attacks to satisfy application requirements. This is achieved by refining the DNN loss function with an attack sensitive matrix to represent the impacts of different targeted attacks. Experimental results on CIFAR-10 and GTSRB data sets demonstrate the efficacy of the proposed solution.


title: On the Modelling of Soft-robots as Quasi-Continuum Lagrangian Dynamical Systems with Well-posed Input Matrix

abstract: In this paper, considering a braided continuum soft-robot, whose radial deformation is constrained but elongation is assumed, a quasi-Lagrangian model is proposed that meets the Lagrangian models properties, including a well-posed input matrix. Actuation is considered throughout three inner pressure cambers, and torsional effects are neglected. The closed-form analytical model is obtained using a scalar varying mass density field, previously neglected in the literature, which produces on one hand a varying center of mass, which generally does not lay in the backbone curve, and one the other hand a coordinate-dependent inertial tensor. The Lagrangian approach enforces the basic skew symmetric property, thus exhibiting passivity. The advantage of dealing with all these effects together display the following distinct features: \textit{i)} the Lagrangian soft-robot dynamic model is similar to the Lagrangian rigid-robot case; \textit{ii)} the non-linear system is affine in the control input; \textit{iii)} the continuum deformable body stands for a segment of constant curvature, when interconnected with other segments of different constant curvature each, would leads to a quasi-continuum $n$-segments variable curvature soft-robot, yet preserving the aforementioned previous features of one segment. \redc{Representative simulations and videos are shown, in open- and closed-loop.


title: Circular orbits of a ball on a rotating conical turntable

abstract: This paper is concerned with the study of the rolling without slipping of a dynamically symmetric (in particular, homogeneous) heavy ball on a cone which rotates uniformly about its symmetry axis. The equations of motion of the system are obtained, partial periodic solutions are found and their stability is analyzed.


title: Gauge invariance for a class of tree diagrams in the standard model

abstract: For gauge theory, the matrix element for any physical process is independent of the gauge used. Since this is a formal statement and examples are known where gauge invariance is violated, for any specific process this gauge invariance needs to be checked by explicit calculation. In this paper, gauge invariance is found to hold for a large non-trivial class of processes described by tree diagrams in the standard model -- tree diagrams with two external $W$ bosons and any number of external Higgs bosons. This verification of gauge invariance is quite complicated, and is based on a direct study of the difference between different gauges through induction on the number of external Higgs bosons.


title: Consistent dust and gas models for protoplanetary disks III. Models for selected objects from the FP7 DIANA project

abstract: The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).


title: Gamma positivity of the Excedance based Eulerian polynomial in positive elements of Classical Weyl Groups

abstract: The classical Eulerian polynomials $A_n(t)$ are known to be gamma positive. Define the positive Eulerian polynomial $\mathsf{AExc^{+}}_n(t)$ as the polynomial obtained when we sum excedances over the alternating group. We show that $\mathsf{AExc^{+}}_n(t)$ is gamma positive iff $n \geq 5$ and $n \equiv 1$ (mod 2). When $n \geq 4$, and $n \equiv 0$ (mod 2) we show that $\mathsf{AExc^{+}}_n(t)$ can be written as a sum of two gamma positive polynomials. Similar results are shown when we consider the positive type-D and type-D Eulerian polynomials. Finally, we show gamma positivity results when we sum excedances over derangements with positive and negative sign. Our main resuls is that the polynomial obtained by summing excedance over a conjugacy class indexed by $\lambda$ is gamma positive.


title: Optimization on fractal sets

abstract: We outline necessary and sufficient condition to the existence of extrmas of a function on a self-similar set, and we describe discrete gradient algorithm to find the extrema.


title: Effective-one-body multipolar waveform for tidally interacting binary neutron stars up to merger

abstract: Gravitational-wave astronomy with coalescing binary neutron star sources requires the availability of gravitational waveforms with tidal effects accurate up to merger. This article presents an improved version of \TEOBResum, a nonspinning effective-one-body (EOB) waveform model with enhanced analytical information in its tidal sector. The tidal potential governing the conservative dynamics employs resummed expressions based on post-Newtonian (PN) and gravitational self-force (GSF) information. In particular, we compute a GSF-resummed expression for the leading-order octupolar gravitoelectric term and incorporate the leading-order gravitomagnetic term (either in PN-expanded or GSF-resummed form). The multipolar waveform and fluxes are augmented with gravitoelectric and magnetic terms recently obtained in PN. The new analytical information enhances tidal effects toward merger accelerating the coalescence. We quantify the impact on the gravitational-wave phasing of each physical effect. The most important contribution is given by the resummed gravitoelectric octupolar term entering the EOB interaction potential, that can yield up to 1~rad of dephasing (depending on the NS model) with respect to its nonresummed version. The model's energetics and the gravitational wave phasing are validated with eccentricity-reduced and multi-resolution numerical relativity simulations with different equations of state and mass ratios. We also present EOB-NR waveform comparisons for higher multipolar modes beyond the dominant quadrupole one.


title: SDSS-IV MaStar -- A Large and Comprehensive Empirical Stellar Spectral Library: First Release

abstract: We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.


title: Bang-bang control as a design principle for classical and quantum optimization algorithms

abstract: Physically motivated classical heuristic optimization algorithms such as simulated annealing (SA) treat the objective function as an energy landscape, and allow walkers to escape local minima. It has been argued that quantum properties such as tunneling may give quantum algorithms advantage in finding ground states of vast, rugged cost landscapes. Indeed, the Quantum Adiabatic Algorithm (QAO) and the recent Quantum Approximate Optimization Algorithm (QAOA) have shown promising results on various problem instances that are considered classically hard. Here, we argue that the type of control strategy used by the optimization algorithm may be crucial to its success. Along with SA, QAO and QAOA, we define a new, bang-bang version of simulated annealing, BBSA, and study the performance of these algorithms on two well-studied problem instances from the literature. Both classically and quantumly, the successful control strategy is found to be bang-bang, exponentially outperforming the quasistatic analogues on the same instances. Lastly, we construct O(1)-depth QAOA protocols for a class of symmetric cost functions, and provide an accompanying physical picture.


title: Gravitational Waves from a Pati-Salam Phase Transition

abstract: We analyse the gravitational wave and low energy signatures of a Pati-Salam phase transition. For a Pati-Salam scale of $M_{PS} \sim 10^5$ GeV, we find a stochastic power spectrum within reach of the next generation of ground-based interferometer experiments such as the Einstein Telescope, in parts of the parameter space. We study the lifetime of the proton in this model, as well as complementarity with low energy constraints including electroweak precision data, neutrino mass measurements, lepton flavour violation, and collider constraints.


title: High-Resolution Transit Spectroscopy of Warm Saturns

abstract: We present high-resolution optical transmission spectroscopy of two sub-Saturn mass transiting exoplanets, HAT-P-12b and WASP-69b. With relatively low densities and high atmospheric scale heights, these planets are particularly well-suited to characterization through transit spectroscopy, and serve as ideal candidates for extending previously-tested methods to lower planetary masses. Using a single transit for each planet, we take advantage of the Doppler cross-correlation technique to search for sodium, potassium, and water absorption features. Our analysis reveals a likely ($3.2\sigma$) detection of sodium absorption features in the atmosphere of HAT-P-12b, and enables us to place constraints on the presence of alkaline and molecular species in the atmospheres of both planets. With our results, we highlight the efficacy of ground-based campaigns for characterizing exoplanetary atmospheres, and pave the way for future analyses of low-mass planets.


title: Be it therefore resolved: Cosmological Simulations of Dwarf Galaxies with Extreme Resolution

abstract: We study a suite of extremely high-resolution cosmological FIRE simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}$$M_{\odot}$), run to $z=0$ with $30 M_{\odot}$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6} M_{\odot}$ is populated by a resolved {\em stellar} galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast}<10^{5}\,M_{\odot}$) have their star formation truncated early ($z\gtrsim2$), likely by reionization, while classical dwarfs ($M_{\ast}>10^{5} M_{\odot}$) continue forming stars to $z<0.5$. The systems have bursty star formation (SF) histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with $M_{\ast}/M_{\rm halo} > 10^{-4}$ to form a dark matter core $>200$pc, while lower-mass UFDs exhibit cusps down to $\lesssim100$pc, as expected from energetic arguments. Our dwarfs with $M_{\ast}>10^{4}\,M_{\odot}$ have half-mass radii ($R_{\rm 1/2}$) in agreement with Local Group (LG) dwarfs; dynamical mass vs. $R_{1/2}$ and the degree of rotational support also resemble observations. The lowest-mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).