cssi/SciSciGPT-SciSciNet huggingface dataset:
MAPPING THE ASYMMETRIC THICK DISK. II. DISTANCE, SIZE, AND MASS OF THE HERCULES THICK DISK CLOUD

The Hercules Thick Disk Cloud was initially discovered as an excess in the number of faint blue stars between Quadrants 1 and 4 of the Galaxy. The origin of the Cloud could be an interaction with the disk bar, a triaxial Thick Disk, or a merger remnant or stream. To better map the spatial extent of the Cloud along the line of sight, we have obtained multi-color UBVR photometry for 1.2 million stars in 63 fields each of approximately 1 deg{sup 2}. Our analysis of the fields beyond the apparent boundaries of the excess has already ruled out a triaxial Thick Disk as a likely explanation. In this paper, we present our results for the star counts over all of our fields, determine the spatial extent of the overdensity across and along the line of sight, and estimate the size and mass of the Cloud. Using photometric parallaxes, the stars responsible for the excess are between 1 and 6 kpc from the Sun, 0.5-4 kpc above the Galactic plane, and extend approximately 3-4 kpc across our line of sight. The Cloud is thus a major substructure in the Galaxy. The distribution of the excess along our sight lines corresponds withmore » the density contours of the bar in the Disk, and its most distant stars are directly over the bar. We also see through the Cloud to its far side. Over the entire 500 deg{sup 2} of the sky containing the Cloud, we estimate more than 5.6 million stars and 1.9 million solar masses of material. If the overdensity is associated with the bar, it would exceed 1.4 billion stars and more than 50 million solar masses. Finally, we argue that the Hercules-Aquila Cloud is actually the Hercules Thick Disk Cloud.« less
Viewing A Program Transformation System At Work

How to decrease labor and improve reliability in the development of efficient implementations of nonnumerical algorithms and labor intensive software is an increasingly important problem as the demand for computer technology shifts from easier applications to more complex algorithmic ones; e.g., optimizing compilers for supercomputers, intricate data structures to implement efficient solutions to operations research problems, search and analysis algorithms in genetic engineering, complex software tools for workstations, design automation, etc. It is also a difficult problem that is not solved by current CASE tools and software management disciplines, which are oriented towards data processing and other applications, where the implementation and a prediction of its resource utilization follow more directly from the specification.
A large N duality via a geometric transition

Abstract We propose a large N dual of 4d, N =1 supersymmetric, SU(N) Yang–Mills with adjoint field Φ and arbitrary superpotential W(Φ). The field theory is geometrically engineered via D-branes partially wrapped over certain cycles of a non-trivial Calabi–Yau geometry. The large N, or low-energy, dual arises from a geometric transition of the Calabi–Yau, where the branes have disappeared and have been replaced by suitable fluxes. This duality yields highly non-trivial exact results for the gauge theory. The predictions indeed agree with expected results in cases where it is possible to use standard techniques for analyzing the strongly coupled, supersymmetric gauge theories. Moreover, the proposed large N dual provides a simpler and more unified approach for obtaining exact results for this class of supersymmetric gauge theories.
Using astronomical photographs to investigate misconceptions about galaxies and spectra: Question development for clicker use

Development of multiple-choice questions to identify student misconceptions about galaxies and spectra.
Gaussian Random Field: Physical Origin of Sersic Profiles

While the Sersic profile family provides adequate fits for the surface brightness profiles of observed galaxies, its physical origin is unknown. We show that if the cosmological density field is seeded by random Gaussian fluctuations, as in the standard cold dark matter model, galaxies with steep central profiles have simultaneously extended envelopes of shallow profiles in the outskirts, whereas galaxies with shallow central profiles are accompanied by steep density profiles in the outskirts. These properties are in accord with those of the Sersic profile family. Moreover, galaxies with steep central profiles form their central regions in smaller denser subunits that possibly merge subsequently, which naturally leads to the formation of bulges. In contrast, galaxies with shallow central profiles form their central regions in a coherent fashion without significant substructure, a necessary condition for disk galaxy formation. Thus, the scenario is self-consistent with respect to the correlation between observed galaxy morphology and the Sersic index. We further predict that clusters of galaxies should display a similar trend, which should be verifiable observationally.
Asymptotics of quantum channels: conserved quantities, an adiabatic limit, and matrix product states

This work derives an analytical formula for the asymptotic state---the quantum state resulting from an infinite number of applications of a general quantum channel on some initial state. For channels admitting multiple fixed or rotating points, conserved quantities---the left fixed/rotating points of the channel---determine the dependence of the asymptotic state on the initial state. The formula stems from a Noether-like theorem stating that, for any channel admitting a full-rank fixed point, conserved quantities commute with that channel’s Kraus operators up to a phase. The formula is applied to adiabatic transport of the fixed-point space of channels, revealing cases where the dissipative/spectral gap can close during any segment of the adiabatic path. The formula is also applied to calculate expectation values of noninjective matrix product states (MPS) in the thermodynamic limit, revealing that those expectation values can also be calculated using an MPS with reduced bond dimension and a modified boundary.
Tachyon condensation using the disc partition function

It has been recently proposed that the background independent open superstring field theory action is given by the disc partition function with all possible open string operators inserted at the boundary of the disc. We use this proposal to study tachyon condensation in the D0-D2 system. We evaluate the disc partition function for the D0-D2 system in presence of a large Neveu-Schwarz B-field using perturbation theory. This perturbative expansion of the disc partition function makes sense as the boundary tachyon operator for the large Neveu-Schwarz B-field is almost marginal. We find that the mass defect for the formation of the D0-D2 bound state agrees exactly with the expected result in the large B-field limit.
Ghera: A Repository of Android App Vulnerability Benchmarks

Security of mobile apps affects the security of their users. This has fueled the development of techniques to automatically detect vulnerabilities in mobile apps and help developers secure their apps; specifically, in the context of Android platform due to openness and ubiquitousness of the platform. Despite a slew of research efforts in this space, there is no comprehensive repository of up-to-date and lean benchmarks that contain most of the known Android app vulnerabilities and, consequently, can be used to rigorously evaluate both existing and new vulnerability detection techniques and help developers learn about Android app vulnerabilities. In this paper, we describe Ghera, an open source repository of benchmarks that capture 25 known vulnerabilities in Android apps (as pairs of exploited/benign and exploiting/malicious apps). We also present desirable characteristics of vulnerability benchmarks and repositories that we uncovered while creating Ghera.
Domain wall fermion study of scaling in non-perturbative renormalization of quark composite operators

Abstract We compute non-perturbatively the renormalization coefficients of scalar and pseudoscalar operators, local vector and axial currents, conserved vector and axial current, and O Δ S = 2 LL over a wide range of energy scales using a scaling technique that connects the results of simulations at different values of coupling β. We use the domain wall fermion formulation in the quenched approximation at a series of three values of β = 6.0, 6.45, 7.05 corresponding to lattice spacing scaling by factors of two.
Hierarchical Predictive Control Algorithms for Optimal Design and Operation of Microgrids

In recent years, microgrids, i.e., disconnected distribution systems, have received increasing interest from power system utilities to support the economic and resiliency posture of their systems. The economics of long distance transmission lines prevent many remote communities from connecting to bulk transmission systems and these communities rely on off-grid microgrid technology. Furthermore, communities that are connected to the bulk transmission system are investigating microgrid technologies that will support their ability to disconnect and operate independently during extreme events. In each of these cases, it is important to develop methodologies that support the capability to design and operate microgrids in the absence of transmission over long periods of time. Unfortunately, such planning problems tend to be computationally difficult to solve and those that are straightforward to solve often lack the modeling fidelity that inspires confidence in the results. To address these issues, we first develop a high fidelity model for design and operations of a microgrid that include component efficiencies, component operating limits, battery modeling, unit commitment, capacity expansion, and power flow physics; the resulting model is a mixed-integer quadratically-constrained quadratic program (MIQCQP). We then develop an iterative algorithm, referred to as the Model Predictive Control (MPC) algorithm, that allows us to solve the resulting MIQCQP. We show, through extensive computational experiments, that the MPC-based method can scale to problems that have a very long planning horizon and provide high quality solutions that lie within 5 % of optimal.
Asteroseismic signatures of magnetic activity variations in solar-type stars

Observations of magnetic activity cycles in other stars provide a broader context for our understanding of the 11-year sunspot cycle. The discovery of short activity cycles in a few stars, and the recognition of analogous variability in the Sun, suggest that there may be two distinct dynamos operating in different regions of the interior. Consequently, there is a natural link between studies of magnetic activity and asteroseismology, which can characterize some of the internal properties that are relevant to dynamos. I provide a brief historical overview of the connection between these two fields (including prescient work by Wojtek Dziembowski in 2007), and I highlight some exciting results that are beginning to emerge from the Kepler mission.
Walking dynamics are symmetric (enough)

Many biological phenomena such as locomotion, circadian cycles and breathing are rhythmic in nature and can be modelled as rhythmic dynamical systems. Dynamical systems modelling often involves neglecting certain characteristics of a physical system as a modelling convenience. For example, human locomotion is frequently treated as symmetric about the sagittal plane. In this work, we test this assumption by examining human walking dynamics around the steady state (limit-cycle). Here, we adapt statistical cross-validation in order to examine whether there are statistically significant asymmetries and, even if so, test the consequences of assuming bilateral symmetry anyway. Indeed, we identify significant asymmetries in the dynamics of human walking, but nevertheless show that ignoring these asymmetries results in a more consistent and predictive model. In general, neglecting evident characteristics of a system can be more than a modelling convenience—it can produce a better model.
The rank of a hypergeometric system

The holonomic rank of the A -hypergeometric system M A ( β ) is the degree of the toric ideal I A for generic parameters; in general, this is only a lower bound. To the semigroup ring of A we attach the ranking arrangement and use this algebraic invariant and the exceptional arrangement of non-generic parameters to construct a combinatorial formula for the rank jump of M A ( β ). As consequences, we obtain a refinement of the stratification of the exceptional arrangement by the rank of M A ( β ) and show that the Zariski closure of each of its strata is a union of translates of linear subspaces of the parameter space. These results hold for generalized A -hypergeometric systems as well, where the semigroup ring of A is replaced by a non-trivial weakly toric module M ⊆ℂ[ℤ A ] . We also provide a direct proof of the main result in [M. Saito, Isomorphism classes of A-hypergeometric systems , Compositio Math. 128 (2001), 323–338] regarding the isomorphism classes of M A ( β ) .
Polyhedral approximation in mixed-integer convex optimization

Generalizing both mixed-integer linear optimization and convex optimization, mixed-integer convex optimization possesses broad modeling power but has seen relatively few advances in general-purpose solvers in recent years. In this paper, we intend to provide a broadly accessible introduction to our recent work in developing algorithms and software for this problem class. Our approach is based on constructing polyhedral outer approximations of the convex constraints, resulting in a global solution by solving a finite number of mixed-integer linear and continuous convex subproblems. The key advance we present is to strengthen the polyhedral approximations by constructing them in a higher-dimensional space. In order to automate this extended formulation we rely on the algebraic modeling technique of disciplined convex programming (DCP), and for generality and ease of implementation we use conic representations of the convex constraints. Although our framework requires a manual translation of existing models into DCP form, after performing this transformation on the MINLPLIB2 benchmark library we were able to solve a number of unsolved instances and on many other instances achieve superior performance compared with state-of-the-art solvers like Bonmin, SCIP, and Artelys Knitro.
Reionization history and CMB parameter estimation

We study how uncertainty in the reionization history of the universe affects estimates of other cosmological parameters from the Cosmic Microwave Background. We analyze WMAP7 data and synthetic Planck-quality data generated using a realistic scenario for the reionization history of the universe obtained from high-resolution numerical simulation. We perform parameter estimation using a simple sudden reionization approximation, and using the Principal Component Analysis (PCA) technique proposed by Mortonson and Hu. We reach two main conclusions: (1) Adopting a simple sudden reionization model does not introduce measurable bias into values for other parameters, indicating that detailed modeling of reionization is not necessary for the purpose of parameter estimation from future CMB data sets such as Planck. (2) PCA analysis does not allow accurate reconstruction of the actual reionization history of the universe in a realistic case.
Learning Colour Representations of Search Queries

Image search engines rely on appropriately designed ranking features that capture various aspects of the content semantics as well as the historic popularity. In this work, we consider the role of colour in this relevance matching process. Our work is motivated by the observation that a significant fraction of user queries have an inherent colour associated with them. While some queries contain explicit colour mentions (such as 'black car' and 'yellow daisies'), other queries have implicit notions of colour (such as 'sky' and 'grass'). Furthermore, grounding queries in colour is not a mapping to a single colour, but a distribution in colour space. For instance, a search for 'trees' tends to have a bimodal distribution around the colours green and brown. We leverage historical clickthrough data to produce a colour representation for search queries and propose a recurrent neural network architecture to encode unseen queries into colour space. We also show how this embedding can be learnt alongside a cross-modal relevance ranker from impression logs where a subset of the result images were clicked. We demonstrate that the use of a query-image colour distance feature leads to an improvement in the ranker performance as measured by users' preferences of clicked versus skipped images.
Fano collective resonance as complex mode in a two-dimensional planar metasurface of plasmonic nanoparticles

Fano resonances are features in transmissivity/reflectivity/absorption that owe their origin to the interaction between a broad bright resonance and a dark (i.e., sub-radiant) narrower one. They may emerge in the optical properties of planar two-dimensional (2D) periodic arrays (metasurfaces) of plasmonic nanoparticles. In this letter, we provide a thorough assessment of their nature for the general case of normal and oblique plane wave incidence, highlighting when a Fano resonance is affected by the mutual coupling in an array and its capability to support free modal solutions. We analyze the representative case of a metasurface of plasmonic nanoshells at ultraviolet frequencies and compute its absorption under TE- and TM-polarized, oblique plane-wave incidence. In particular, we find that plasmonic metasurfaces display two distinct types of resonances observable as absorption peaks: one is related to the Mie electric dipolar resonance of each nanoparticle and the other is due to the forced excitation of free modes with small attenuation constant, usually found at oblique incidence. The latter is thus an array-induced collective Fano resonance. This realization opens up to manifold flexible designs at optical frequencies mixing individual and collective resonances. We explain the physical origin of such Fano resonances using the modal analysis through which we calculate the free modes with complex wavenumber supported by the metasurface. We define equivalent array dipolar polarizabilities that are directly related to the absorption physics at oblique incidence and show a direct dependence between array modal phase and attenuation constants and Fano resonances. We thus provide a more complete picture of Fano resonances that may lead to the design of filters, energy-harvesting devices, photodetectors, and sensors at ultraviolet frequencies. Similar resonances may be also extended to the visible range with an appropriate choice of geometries and materials.
Channel activation of CHSH nonlocality

Quantum channels that break CHSH nonlocality on all input states are known as CHSH-breaking channels. In quantum networks, such channels are useless for distributing correlations that can violate the CHSH Inequality. Motivated by previous work on activation of nonlocality in quantum states, here we demonstrate an analogous activation of CHSH-breaking channels. That is, we show that certain pairs of CHSH-breaking channels are no longer CHSH-breaking when used in combination. We find that this type of activation can emerge in both uni-directional and bi-directional communication scenarios.
Constrained $H^1$-Regularization Schemes for Diffeomorphic Image Registration

We propose regularization schemes for deformable registration and efficient algorithms for their numerical approximation. We treat image registration as a variational optimal control problem. The deformation map is parametrized by its velocity. Tikhonov regularization ensures well-posedness. Our scheme augments standard smoothness regularization operators based on $H^1$- and $H^2$-seminorms with a constraint on the divergence of the velocity field, which resembles variational formulations for Stokes incompressible flows. In our formulation, we invert for a stationary velocity field and a mass source map. This allows us to explicitly control the compressibility of the deformation map and by that the determinant of the deformation gradient. We also introduce a new regularization scheme that allows us to control shear. We use a globalized, preconditioned, matrix-free, reduced space (Gauss--)Newton--Krylov scheme for numerical optimization. We exploit variable elimination techniques to reduce the number of un...
Spanning trees of graphs on surfaces and the intensity of loop-erased random walk on planar graphs

We show how to compute the probabilities of various connection topologies for uniformly random spanning trees on graphs embedded in surfaces. As an application, we show how to compute the "intensity" of the loop-erased random walk in ${\mathbb Z}^2$, that is, the probability that the walk from (0,0) to infinity passes through a given vertex or edge. For example, the probability that it passes through (1,0) is 5/16; this confirms a conjecture from 1994 about the stationary sandpile density on ${\mathbb Z}^2$. We do the analogous computation for the triangular lattice, honeycomb lattice and ${\mathbb Z} \times {\mathbb R}$, for which the probabilities are 5/18, 13/36, and $1/4-1/\pi^2$ respectively.
Polynomial reduction and evaluation of tree- and loop-level CHY amplitudes

We develop a polynomial reduction procedure that transforms any gauge fixed CHY amplitude integrand for n scattering particles into a σ-moduli multivariate polynomial of what we call the standard form. We show that a standard form polynomial must have a specific ladder type monomial structure, which has finite size at any n, with highest multivariate degree given by (n − 3)(n − 4)/2. This set of monomials spans a complete basis for polynomials with rational coefficients in kinematic data on the support of scattering equations. Subsequently, at tree and one-loop level, we employ the global residue theorem to derive a prescription that evaluates any CHY amplitude by means of collecting simple residues at infinity only. The prescription is then applied explicitly to some tree and one-loop amplitude examples.
Online learning with dynamics: A minimax perspective

We study the problem of online learning with dynamics, where a learner interacts with a stateful environment over multiple rounds. In each round of the interaction, the learner selects a policy to deploy and incurs a cost that depends on both the chosen policy and current state of the world. The state-evolution dynamics and the costs are allowed to be time-varying, in a possibly adversarial way. In this setting, we study the problem of minimizing policy regret and provide non-constructive upper bounds on the minimax rate for the problem. Our main results provide sufficient conditions for online learnability for this setup with corresponding rates. The rates are characterized by 1) a complexity term capturing the expressiveness of the underlying policy class under the dynamics of state change, and 2) a dynamics stability term measuring the deviation of the instantaneous loss from a certain counterfactual loss. Further, we provide matching lower bounds which show that both the complexity terms are indeed necessary. Our approach provides a unifying analysis that recovers regret bounds for several well studied problems including online learning with memory, online control of linear quadratic regulators, online Markov decision processes, and tracking adversarial targets. In addition, we show how our tools help obtain tight regret bounds for a new problems (with non-linear dynamics and non-convex losses) for which such bounds were not known prior to our work.
Roles for the RNA polymerase III regulator MAFR-1 in regulating sperm quality in Caenorhabditis elegans

The negative regulator of RNA polymerase (pol) III mafr-1 has been shown to affect RNA pol III transcript abundance, lipid biosynthesis and storage, progeny output, and lifespan. We deleted mafr-1 from the Caenorhabditis elegans genome and found that animals lacking mafr-1 replicated many phenotypes from previous RNAi-based studies and discovered a new sperm-specific role. Utilizing a yeast two-hybrid assay, we discovered several novel interactors of MAFR-1 that are expressed in a sperm- and germline-enriched manner. In support of a role for MAFR-1 in the male germline, we found mafr-1 null males have smaller spermatids that are less capable in competition for fertilization; a phenotype that was dependent on RNA pol III activity. Restoration of MAFR-1 expression specifically in the germline rescued the spermatid-related phenotypes, suggesting a cell autonomous role for MAFR-1 in nematode male fertility. Based on the high degree of conservation of Maf1 activity across species, our study may inform similar roles for Maf1 and RNA pol III in mammalian male fertility.
ELECTROMAGNETIC EXTRACTION OF ENERGY FROM BLACK-HOLE–NEUTRON-STAR BINARIES

The coalescence of black-hole-neutron-star binaries is expected to be a principal source of gravitational waves for the next generation of detectors, Advanced LIGO and Advanced Virgo. For black hole masses not much larger than the neutron star mass, the tidal disruption of the neutron star by the black hole provides one avenue for generating an electromagnetic counterpart. However, in this work, we demonstrate that, for all black-hole-neutron-star binaries observable by Advanced LIGO/Virgo, the interaction of the black hole with the magnetic field of the neutron star will generate copious luminosity, comparable to supernovae and active galactic nuclei. This novel effect may have already been observed as a new class of very short gamma-ray bursts by the Swift Gamma-Ray Burst Telescope. These events may be observable to cosmological distances, so that any black-hole-neutron-star coalescence detectable with gravitational waves by Advanced LIGO/Virgo could also be detectable electromagnetically.
Small generating sets for the Torelli group

Proving a conjecture of Dennis Johnson, we show that the Torelli subgroup Ig of the genus g mapping class group has a finite generating set whose size grows cubically with respect to g . Our main tool is a new space called the handle graph on which Ig acts cocompactly. 20F05; 20F38, 57M07, 57N05 Let Ug;n be a compact connected oriented genus g surface with n boundary components. The mapping class group of Ug;n , denoted Modg;n , is the group of orientationpreserving homeomorphisms of Ug;n that fix the boundary pointwise modulo isotopies that fix the boundary pointwise. We will often omit the n if it vanishes. For n 1, the Torelli group, denoted Ig;n , is the kernel of the action of Modg;n on H1.Ug;nIZ/. The Torelli group has been the object of intensive study ever since the seminal work of Dennis Johnson in the early ’80’s. See [10] for a survey of Johnson’s work.
Higher cohomology of parabolic actions on certain homogeneous spaces

We show that for a parabolic Rd-action on PSL(2,R)d/Γ, the cohomologies in degrees 1 through d − 1 trivialize, and we give the obstructions to solving the degree-d coboundary equation, along with bounds on Sobolev norms of primitives. In previous papers, we have established these results for certain Anosov systems. This work extends the methods of those papers to systems that are not Anosov. The main new idea is defining special elements of representation spaces that allow us to modify the arguments from the previous papers. We discuss how to generalize this strategy to Rd-systems coming from a product of Lie groups, as in the systems analyzed here.
IRON OPACITY BUMP CHANGES THE STABILITY AND STRUCTURE OF ACCRETION DISKS IN ACTIVE GALACTIC NUCLEI

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be larger than the electron scattering opacity due to the large number of bound–bound transitions in iron. We study the effects of this iron opacity “bump” on the thermal stability and vertical structure of radiation-pressure-dominated accretion disks, utilizing three-dimensional radiation magnetohydrodynamic (MHD) simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a 5 × 10(8) solar mass black hole with ∼3% of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than 10 thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free–free opacity are included as in the standard thin disk model, the disk collapses on the thermal timescale. The difference is caused by a combination of (1) an anti-correlation between the total optical depth and the midplane pressure, and (2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of active galactic nucleus (AGN) accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.
A broadband fiber-optic nonlinear interferometer

We describe an all-fiber nonlinear interferometer based on four-wave mixing in highly nonlinear fiber. Our configuration realizes phase-sensitive interference with 97% peak visibility and >90% visibility over a broad 554 GHz optical band. By comparing the output noise power to the shot-noise level, we confirm noise cancellation at dark interference fringes, as required for quantum-enhanced sensitivity. Our device extends nonlinear interferometry to the important platform of highly nonlinear optical fiber and could find applications in a variety of fiber-based sensors.
Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

Author(s): Medford, MS; Lu, JR; Dawson, WA; Lam, CY; Golovich, NR; Schlafly, EF; Nugent, P | Abstract: © 2020. The American Astronomical Society. All rights reserved. Microlensing surveys have discovered thousands of events, with almost all events discovered within the Galactic bulge or toward the Magellanic Clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky, including the Galactic plane, every few nights. The ZTF observes ∼109 stars in the g and r bands and can significantly contribute to the observed microlensing population. We predict that the ZTF will observe ∼1100 microlensing events in 3 yr of observing within 10° latitude of the Galactic plane, with ∼500 events in the outer Galaxy (l ≥ 10°). This yield increases to ∼1400 (∼800) events by combining every three ZTF exposures, ∼1800 (∼900) events if the ZTF observes for a total of 5 yr, and ∼2400 (∼1300) events for a 5 yr survey with postprocessing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. The ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.
Photosynthesis on a Planet Orbiting an M Dwarf: Enhanced Effectiveness during Flares

On planets near M dwarfs, oxygenic photosynthesis (PS) will occur with an effectiveness that depends on the supply of visible photons with wavelengths between 400 and 700 nm. In this paper, we quantify the effectiveness of PS in two contexts that are relevant for M dwarfs. First, using photons from an M dwarf in its quiescent (nonflaring) state, we find that PS on an M dwarf planet in the habitable zone (HZ) of its parent star is less effective than on Earth by a factor of 10 for a flare star with mid-M spectral type. For a flare star with late-M spectral type, PS effectiveness (PSE) is smaller than on Earth by a factor of 100 or more. Second, using photons that are incident on the HZ planet during flares, we find that PSE can increase by factors of 5–20 above the quiescent values. In the case of a flare star with mid-M spectral type, we find that the PSE during a flare can increase up to as much as 50%–60% of the values on Earth. However, for a late-M flare star, even during flares, the PSE remains almost one order of magnitude smaller than on Earth. We suggest that for biological processes on M dwarf planets, the stellar activity cycle may replace the orbital period as the "year."
AFFINE GEOMETRY OF STRATA OF DIFFERENTIALS

Affine varieties among all algebraic varieties have simple structures. For example, an affine variety does not contain any complete algebraic curve. In this paper, we study affine-related properties of strata of , then the corresponding stratum does not contain any complete curve. Moreover, we explore the amusing question whether affine invariant manifolds arising from Teichmuller dynamics are affine varieties, and confirm the answer for Teichmuller curves, Hurwitz spaces of torus coverings, hyperelliptic strata as well as some low genus strata.
Heavy quark limit in the light front quark model

An explicit relativistic light-front quark model is presented which gives the momentum transfer dependent form factors of weak hadronic currents among heavy pseudoscalar and vector mesons in the whole accessible kinematical region 0 ≤ q2 ≤ q2max. It is shown that in the limit of infinite masses of active quarks these form factors can be expressed in terms of universal Isgur- Wise function. The explisit expression for this function is obtained. It is shown that neglect of pair creation from the vacuum in calculations of form factors does not violate Luke's theorem.
The Acceleration of Thermal Protons at Parallel Collisionless Shocks: Three-dimensional Hybrid Simulations

We present three-dimensional hybrid simulations of collisionless shocks that propagate parallel to the background magnetic field to study the acceleration of protons that forms a high-energy tail on the distribution. We focus on the initial acceleration of thermal protons and compare it with results from one-dimensional simulations. We find that for both one- and three-dimensional simulations, particles that end up in the high-energy tail of the distribution later in the simulation gained their initial energy right at the shock. This confirms previous results but is the first to demonstrate this using fully three-dimensional fields. The result is not consistent with the "thermal leakage" model. We also show that the gyrocenters of protons in the three-dimensional simulation can drift away from the magnetic field lines on which they started due to the removal of ignorable coordinates that exist in one- and two-dimensional simulations. Our study clarifies the injection problem for diffusive shock acceleration.
Measuring galaxy clustering and the evolution of [C II] mean intensity with far-IR line intensity mapping during 0.5 < z < 1.5

Infrared fine-structure emission lines from trace metals are powerful diagnostics of the interstellar medium in galaxies. We explore the possibility of studying the redshifted far-IR fine-structure line emission using the three-dimensional (3D) power spectra obtained with an imaging spectrometer. The intensity mapping approach measures the spatio-spectral fluctuations due to line emission from all galaxies, including those below the individual detection threshold. The technique provides 3D measurements of galaxy clustering and moments of the galaxy luminosity function. Furthermore, the linear portion of the power spectrum can be used to measure the total line emission intensity including all sources through cosmic time with redshift information naturally encoded. Total line emission, when compared to the total star formation activity and/or other line intensities, reveals evolution of the interstellar conditions of galaxies in aggregate. As a case study, we consider measurement of [C II] autocorrelation in the 0.5 < z < 1.5 epoch, where interloper lines are minimized, using far-IR/submillimeter balloon-borne and future space-borne instruments with moderate and high sensitivity, respectively. In this context, we compare the intensity mapping approach to blind galaxy surveys based on individual detections. We find that intensity mapping is nearly always the best way to obtain the total linemore » emission because blind, wide-field galaxy surveys lack sufficient depth and deep pencil beams do not observe enough galaxies in the requisite luminosity and redshift bins. Also, intensity mapping is often the most efficient way to measure the power spectrum shape, depending on the details of the luminosity function and the telescope aperture.« less
What Can (and Can't) we Do with Sparse Polynomials?

Simply put, a sparse polynomial is one whose zero coefficients are not explicitly stored. Such objects are ubiquitous in exact computing, and so naturally we would like to have efficient algorithms to handle them. However, with this compact storage comes new algorithmic challenges, as fast algorithms for dense polynomials may no longer be efficient. In this tutorial we examine the state of the art for sparse polynomial algorithms in three areas: arithmetic, interpolation, and factorization. The aim is to highlight recent progress both in theory and in practice, as well as opportunities for future work.
Many-Body Quantum Spin Dynamics with Monte Carlo Trajectories on a Discrete Phase Space

Following the time evolution of a large and high-dimensional quantum system is often computationally intractable. A new semiclassical method of modeling quantum dynamics yields surprisingly good results for nonequilibrium dynamics.
Coupling a Surface Acoustic Wave to an Electron Spin in Diamond via a Dark State

Coupling artificial atoms and acoustic waves may be key in future quantum information processing efforts. An experimental breakthrough in coupling nitrogen vacancy centers strongly to acoustic waves in a way that still preserves their spin coherence is reported.
Quadrupole topological photonic crystals

Quadrupole topological phases, exhibiting protected boundary states that are themselves topological insulators of lower dimensions, have recently been of great interest. Extensions of these ideas from current tight binding models to continuum theories for realistic materials require the identification of quantized invariants describing the bulk quadrupole order. Here we identify the analog of quadrupole order in Maxwell’s equations for a gyromagnetic photonic crystal (PhC) through a double-band-inversion process. The quadrupole moment is quantized by the simultaneous presence of crystalline symmetry and broken time-reversal symmetry, which is confirmed using three independent methods: analysis of symmetry eigenvalues, numerical calculations of the nested Wannier bands and the expectation value of the quadrupole operator. Furthermore, we reveal the boundary manifestations of quadrupole phases as quantized edge polarizations and fractional corner charges. The latter are the consequence of a filling anomaly of energy bands as first predicted in electronic systems. Most higher order topological phases are realized by emulations of tight binding models. Extending these concepts to continuum theories requires the identification of invariants describing the bulk multipole order. Here the authors realize the analog of quadrupole order for a gyromagnetic photonic crystal.
CONFORMAL MODELS OF TWO-DIMENSIONAL TURBULENCE

Polyakov recently showed how to use conformal field theory to describe two-dimensional turbulence. Here we construct an infinite hierarchy of solutions, both for the constant enstrophy flux cascade, and the constant energy flux cascade. We conclude with some speculations concerning the stability and physical meaning of these solutions.
Maximum Entropy Method for Solving the Turbulent Channel Flow Problem

There are two components in this work that allow for solutions of the turbulent channel flow problem: One is the Galilean-transformed Navier-Stokes equation which gives a theoretical expression for the Reynolds stress (u'v'); and the second the maximum entropy principle which provides the spatial distribution of turbulent kinetic energy. The first concept transforms the momentum balance for a control volume moving at the local mean velocity, breaking the momentum exchange down to its basic components, u'v', u'2, pressure and viscous forces. The Reynolds stress gradient budget confirms this alternative interpretation of the turbulence momentum balance, as validated with DNS data. The second concept of maximum entropy principle states that turbulent kinetic energy in fully-developed flows will distribute itself until the maximum entropy is attained while conforming to the physical constraints. By equating the maximum entropy state with maximum allowable (viscous) dissipation at a given Reynolds number, along with other constraints, we arrive at function forms (inner and outer) for the turbulent kinetic energy. This allows us to compute the Reynolds stress, then integrate it to obtain the velocity profiles in channel flows. The results agree well with direct numerical simulation (DNS) data at Reτ = 400 and 1000.
Improving Low Compute Language Modeling with In-Domain Embedding Initialisation

Many NLP applications, such as biomedical data and technical support, have 10-100 million tokens of in-domain data and limited computational resources for learning from it. How should we train a language model in this scenario? Most language modeling research considers either a small dataset with a closed vocabulary (like the standard 1 million token Penn Treebank), or the whole web with byte-pair encoding. We show that for our target setting in English, initialising and freezing input embeddings using in-domain data can improve language model performance by providing a useful representation of rare words, and this pattern holds across several different domains. In the process, we show that the standard convention of tying input and output embeddings does not improve perplexity when initializing with embeddings trained on in-domain data.
Safety Control Synthesis with Input Limits: a Hybrid Approach

We introduce a hybrid (discrete–continuous) safety controller which enforces strict state and input constraints on a system–but only acts when necessary, preserving transparent operation of the original system within some safe region of the state space. We define this space using a Min-Quadratic Barrier function, which we construct along the equilibrium manifold using the Lyapunov functions which result from linear matrix inequality controller synthesis for locally valid uncertain linearizations. We also introduce the concept of a barrier pair, which makes it easy to extend the approach to include trajectory-based augmentations to the safe region, in the style of LQR-Trees. We demonstrate our controller and barrier pair synthesis method in simulation-based examples.
Controlled Hallucinations: Learning to Generate Faithfully from Noisy Data

Neural text generation (data- or text-to-text) demonstrates remarkable performance when training data is abundant which for many applications is not the case. To collect a large corpus of parallel data, heuristic rules are often used but they inevitably let noise into the data, such as phrases in the output which cannot be explained by the input. Consequently, models pick up on the noise and may hallucinate–generate fluent but unsupported text. Our contribution is a simple but powerful technique to treat such hallucinations as a controllable aspect of the generated text, without dismissing any input and without modifying the model architecture. On the WikiBio corpus (Lebret et al., 2016), a particularly noisy dataset, we demonstrate the efficacy of the technique both in an automatic and in a human evaluation.
Stellar Populations and the Star Formation Histories of LSB Galaxies: IV Spitzer Surface Photometry of LSB Galaxies

Surface photometry at 3.6 μm is presented for 61 low surface brightness (LSB) galaxies (μ o <19 3.6 μm mag arcsecs −2 ). The sample covers a range of luminosity from −11 to −22 in M 3.6 and size from 1 to 25 kpc. The morphologies in the mid-IR are comparable to those in the optical with 3.6 μm imaging reaches similar surface brightness depth as ground-based optical imaging. A majority of the resulting surface brightness profiles are single exponential in shape with very few displaying upward or downward breaks. The mean V − 3.6 colour of LSB is 2.3 with a standard deviation of 0.5. Colour-magnitude and two-colour diagrams are well matched to models of constant star formation, where the spread in colour is due to small changes in the star formation rate (SFR) over the last 0.5 Gyrs as also suggested by the specific SFR measured by Hα.
New gravitational memories

The conventional gravitational memory effect is a relative displacement in the position of two detectors induced by radiative energy flux. We find a new type of gravitational ‘spin memory’ in which beams on clockwise and counterclockwise orbits acquire a relative delay induced by radiative angular momentum flux. It has recently been shown that the displacement memory formula is a Fourier transform in time of Weinberg’s soft graviton theorem. Here we see that the spin memory formula is a Fourier transform in time of the recently-discovered subleading soft graviton theorem.
Detection of KS -band Thermal Emission from WASP-3b

We report the detection of thermal emission from the hot Jupiter WASP-3b in the K{sub S} band, using a newly developed guiding scheme for the WIRC instrument at the Palomar Hale 200 inch telescope. Our new guiding scheme has improved the telescope guiding precision by a factor of {approx}5-7, significantly reducing the correlated systematics in the measured light curves. This results in the detection of a secondary eclipse with depth of 0.181% {+-} 0.020% (9{sigma})-a significant improvement in WIRC's photometric precision and a demonstration of the capability of Palomar/WIRC to produce high-quality measurements of exoplanetary atmospheres. Our measured eclipse depth cannot be explained by model atmospheres with heat redistribution but favors a pure radiative equilibrium case with no redistribution across the surface of the planet. Our measurement also gives an eclipse phase center of 0.5045 {+-} 0.0020, corresponding to an ecos {omega} of 0.0070 {+-} 0.0032. This result is consistent with a circular orbit, although it also suggests that the planet's orbit might be slightly eccentric. The possible non-zero eccentricity provides insight into the tidal circularization process of the star-planet system, but might also have been caused by a second low-mass planet in the system, as suggested by a previousmore » transit timing variation study. More secondary eclipse observations, especially at multiple wavelengths, are necessary to determine the temperature-pressure profile of the planet's atmosphere and shed light on its orbital eccentricity.« less
Anticyclotomic p-adic L-function of Central Critical Rankin–Selberg L-value

Let M be an imaginary quadratic field, f a Hecke eigenform on GL2(Q) and \pi the unitary base-change to M of the automorphic representation associated to f. Take a unitary arithmetic Hecke character \chi of M inducing the inverse of the central character of f. The celebrated formula of Waldspurger relates the square of an integral L(f) of f and \chi over the idele class group of M to the central critical value L(1/2,{\pi}\otimes{\chi}). In this paper, we present a construction of a new p-adic L-function that interpolates L(f) over arithmetic characters for a cusp form f in the spirit of the landmark result of Katz where he did that for Eisenstein series.