cssi/SciSciGPT-SciSciNet huggingface dataset:
Construction of a Calibrated Probabilistic Classification Catalog: Application to 50k Variable Sources in the All-Sky Automated Survey

With growing data volumes from synoptic surveys, astronomers necessarily must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of classification purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All-Sky Automated Survey (ASAS), and release the Machine-learned ASAS Classification Catalog (MACC), a 28 class probabilistic classification catalog of 50,124 ASAS sources in the ASAS Catalog of Variable Stars. We estimate that MACC achieves a sub-20% classification error rate and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes.
Beam induced electron cloud resonances in dipole magnetic fields

The buildup of low energy electrons in an accelerator, known as electron cloud, can be severely detrimental to machine performance. Under certain beam conditions, the beam can become resonant with the cloud dynamics, accelerating the buildup of electrons. This paper will examine two such effects: multipacting resonances, in which the cloud development time is resonant with the bunch spacing, and cyclotron resonances, in which the cyclotron period of electrons in a magnetic field is a multiple of bunch spacing. Both resonances have been studied directly in dipole fields using retarding field analyzers installed in the Cornell Electron Storage Ring. These measurements are supported by both analytical models and computer simulations.
Gravitational production of super-Hubble-mass particles: an analytic approach

Through a mechanism similar to perturbative particle scattering, particles of mass mχ larger than the Hubble expansion rate Hinf during inflation can be gravitationally produced at the end of inflation without the exponential suppression powers of exp(−mχ/Hinf ). Here we develop an analytic formalism for computing particle production for such massive particles. We apply our formalism to specific models that have been previously been studied only numerically, and we find that our analytical approximations reproduce those numerical estimates well.
We don't need another hero?: the impact of "heroes" on software development

A software project has "Hero Developers" when 80% of contributions are delivered by 20% of the developers. Are such heroes a good idea? Are too many heroes bad for software quality? Is it better to have more/less heroes for different kinds of projects? To answer these questions, we studied 661 open source projects from Public open source software (OSS) Github and 171 projects from an Enterprise Github.We find that hero projects are very common. In fact, as projects grow in size, nearly all projects become hero projects. These findings motivated us to look more closely at the effects of heroes on software development. Analysis shows that the frequency to close issues and bugs are not significantly affected by the presence of heroes or project type (Public or Enterprise). Similarly, the time needed to resolve an issue/bug/enhancement is not affected by heroes or project type. This is a surprising result since, before looking at the data, we expected that increasing heroes on a project will slow down how fast that project reacts to change. However, we do find a statistically significant association between heroes, project types, and enhancement resolution rates. Heroes do not affect enhancement resolution rates in Public projects. However, in Enterprise projects, heroes increase the rate at which projects complete enhancements.In summary, our empirical results call for a revision of a long-held truism in software engineering. Software heroes are far more common and valuable than suggested by the literature, particularly for medium to large Enterprise developments. Organizations should reflect on better ways to find and retain more of these heroes.
Thermal convection in a liquid metal battery

Generation of thermal convection flow in the liquid metal battery, a device recently proposed as a promising solution for the problem of the short-term energy storage, is analyzed using a numerical model. It is found that convection caused by Joule heating of electrolyte during charging or discharging is virtually unavoidable. It exists in laboratory prototypes larger than a few centimeters in size and should become much stronger in larger-scale batteries. The phenomenon needs further investigation in view of its positive (enhanced mixing of reactants) and negative (loss of efficiency and possible disruption of operation due to the flow-induced deformation of the electrolyte layer) effects.
Upper bounds for the number of orbital topological types of planar polynomial vector fields modulo limit cycles

The purpose of this paper is to find an upper bound for the number of orbital topological types of nth-degree polynomial fields in the plane. An obstacle to obtaining such a bound is related to the unsolved second part of the Hilbert 16th problem. This obstacle is avoided by introducing the notion of equivalence modulo limit cycles. Earlier, the author obtained a lower bound of the form \(2^{cn^2 } \). In the present paper, an upper bound of the same form but with a different constant is found. Moreover, for each planar polynomial vector field with finitely many singular points, a marked planar graph is constructed that represents a complete orbital topological invariant of this field.
A Geometric Construction of Solutions to 11D Supergravity

Necessary and sufficient conditions are provided for a class of warped product manifolds with non-vanishing flux to be supersymmetric solutions of 11D supergravity. Many non-compact, but complete solutions can be obtained in this manner, including the multi-membrane solution initially found by Duff and Stelle. In a different direction, an explicit 5-parameter moduli space of solutions to 11D supergravity is also constructed which can be viewed as non-supersymmetric deformations of the Duff–Stelle solution.
Whitehead moves for G-trees

We generalize the familiar notion of a Whitehead move from Culler and Vogtmann's Outer space to the setting of deformation spaces of G-trees. Specifically, we show that there are two moves, each of which transforms a reduced G-tree into another reduced G-tree, that suffice to relate any two reduced trees in the same deformation space. These two moves further factor into three moves between reduced trees that have simple descriptions in terms of graphs of groups. This result has several applications.
Robustness and modular structure in networks

Complex networks have recently attracted much interest due to their prevalence in nature and our daily lives (Vespignani, 2009; Newman, 2010). A critical property of a network is its resilience to random breakdown and failure (Albert et al., 2000; Cohen et al., 2000; Callaway et al., 2000; Cohen et al., 2001), typically studied as a percolation problem (Stauffer & Aharony, 1994; Achlioptas et al., 2009; Chen & D'Souza, 2011) or by modeling cascading failures (Motter, 2004; Buldyrev et al., 2010; Brummitt, et al. 2012). Many complex systems, from power grids and the Internet to the brain and society (Colizza et al., 2007; Vespignani, 2011; Balcan & Vespignani, 2011), can be modeled using modular networks comprised of small, densely connected groups of nodes (Girvan & Newman, 2002). These modules often overlap, with network elements belonging to multiple modules (Palla et al. 2005; Ahn et al. 2010). Yet existing work on robustness has not considered the role of overlapping, modular structure. Here we study the robustness of these systems to the failure of elements. We show analytically and empirically that it is possible for the modules themselves to become uncoupled or non-overlapping well before the network disintegrates. If overlapping modular organization plays a role in overall functionality, networks may be far more vulnerable than predicted by conventional percolation theory.
Proposed high-power beta cells from MgAlB14-type icosahedral-boron semiconductors

Beta cells generate electric power as carrier-producing beta irradiation from incorporated radioisotopes bombard a series of p-n-junctions. However, radiation damage to the semiconductors commonly used in solar cells limits beta cells to extremely weak irradiations that generate concomitantly miniscule electric powers, e.g. micro-Watts. By contrast, beta cells that generate many orders-of-magnitude larger powers are possible with icosahedral boron-rich semiconductors since their bombardment-induced atomic displacements spontaneously self-heal. Furthermore, substitutions for Mg and Al atoms of icosahedral-boron-rich semiconductors based on the MgAlB14 structure can produce p-n junctions as electron transfers from doping-induced interstitial extra-icosahedral atoms convert some normally p-type materials to n-type. Moreover, electron-phonon interactions of the resulting readily displaceable interstitial cations with charge carriers foster their forming large polarons. Oppositely charged polarons repel one another at short range. These repulsions suppress the recombination of n-type with p-type polarons thereby increasing the beta-cell efficiency. All told, use of these icosahedral boron-rich semiconductors could enable beta cells with electric powers that are many orders of magnitude larger than those of existing beta cells. This development opens a new avenue for generating electricity from nuclear decays.
MINIMALITY OF THE DATA IN WAVELET FILTERS

Abstract Orthogonal wavelets, or wavelet frames, for L2( R ) are associated with quadrature mirror filters (QMF), a set of complex numbers which relate the dyadic scaling of functions on R to the Z -translates. In this paper, we show that generically, the data in the QMF-systems of wavelets are minimal, in the sense that the data cannot be nontrivially reduced. The minimality property is given a geometric formulation in the Hilbert space l2( Z ), and it is then shown that minimality corresponds to irreducibility of a wavelet representation of the algebra O 2; and so our result is that this family of representations of O 2 on the Hilbert space l2( Z ) is irreducible for a generic set of values of the parameters which label the wavelet representations.
Generating synthetic data for task-oriented semantic parsing with hierarchical representations

Modern conversational AI systems support natural language understanding for a wide variety of capabilities. While a majority of these tasks can be accomplished using a simple and flat representation of intents and slots, more sophisticated capabilities require complex hierarchical representations supported by semantic parsing. State-of-the-art semantic parsers are trained using supervised learning with data labeled according to a hierarchical schema which might be costly to obtain or not readily available for a new domain. In this work, we explore the possibility of generating synthetic data for neural semantic parsing using a pretrained denoising sequence-to-sequence model (i.e., BART). Specifically, we first extract masked templates from the existing labeled utterances, and then fine-tune BART to generate synthetic utterances conditioning on the extracted templates. Finally, we use an auxiliary parser (AP) to filter the generated utterances. The AP guarantees the quality of the generated data. We show the potential of our approach when evaluating on the Facebook TOP dataset for navigation domain.
Geometric k-nearest neighbor estimation of entropy and mutual information

Nonparametric estimation of mutual information is used in a wide range of scientific problems to quantify dependence between variables. The k-nearest neighbor (knn) methods are consistent, and therefore expected to work well for a large sample size. These methods use geometrically regular local volume elements. This practice allows maximum localization of the volume elements, but can also induce a bias due to a poor description of the local geometry of the underlying probability measure. We introduce a new class of knn estimators that we call geometric knn estimators (g-knn), which use more complex local volume elements to better model the local geometry of the probability measures. As an example of this class of estimators, we develop a g-knn estimator of entropy and mutual information based on elliptical volume elements, capturing the local stretching and compression common to a wide range of dynamical system attractors. A series of numerical examples in which the thickness of the underlying distribution and the sample sizes are varied suggest that local geometry is a source of problems for knn methods such as the Kraskov-Stogbauer-Grassberger estimator when local geometric effects cannot be removed by global preprocessing of the data. The g-knn method performs well despite the manipulation of the local geometry. In addition, the examples suggest that the g-knn estimators can be of particular relevance to applications in which the system is large, but the data size is limited.
Searching for dark photon dark matter with cosmic ray antideuterons

Low energy antideuteron detection presents a unique channel for indirect detection, targeting dark matter that annihilates into hadrons in a relatively background-free way. Since the idea was first proposed, many WIMP-type models have already been disfavored by direct detection experiments, and current constraints indicate that any thermal relic candidates likely annihilate through some hidden sector process. In this paper, we show that cosmic ray antideuteron detection experiments represent one of the best ways to search for hidden sector thermal relic dark matter, and in particular investigate a vector portal dark matter that annihilates via a massive dark photon. We find that the parameter space with thermal relic annihilation and $$ {m}_{\chi }>{m}_{A\prime}\underset{\sim }{>} $$ 20 GeV is largely unconstrained, and near future antideuteron experiment GAPS will be able to probe models in this space with mχ ≈ mA′ up to masses of O (100 GeV). Specifically the dark matter models favored by the Fermi Galactic center excess is expected to be detected or constrained at the 5(3) −σ level assuming a optimistic (conservative) propagation model.
High Energy Physics and Cosmology at the Unification Frontier: Opportunities and Challenges in the coming years

We give here an overview of recent developments in high energy physics and cosmology and their interconnections that relate to unification, and discuss prospects for the future. Thus there are curr...
Dynamically distinguishing polynomials

A polynomial with integer coefficients yields a family of dynamical systems indexed by primes as follows: For any prime p, reduce its coefficients mod p and consider its action on the field $${\mathbb {F}}_p$$ . We say a subset of $${\mathbb {Z}}[x]$$ is dynamically distinguishable mod p if the associated mod p dynamical systems are pairwise non-isomorphic. For any $$k,M\in {\mathbb {Z}}_{>1}$$ , we prove that there are infinitely many sets of integers $${\mathcal {M}}$$ of size M such that $$\left\{ x^k+m\mid m\in {\mathcal {M}}\right\} $$ is dynamically distinguishable mod p for most p (in the sense of natural density). Our proof uses the Galois theory of dynatomic polynomials largely developed by Morton, who proved that the Galois groups of these polynomials are often isomorphic to a particular family of wreath products. In the course of proving our result, we generalize Morton’s work and compute statistics of these wreath products.
Recursive representation of Wronskians in confluent supersymmetric quantum mechanics

A recursive form of arbitrary-order Wronskian associated with transformation functions in the confluent algorithm of supersymmetric quantum mechanics (SUSY) is constructed. With this recursive form regularity conditions for the generated potentials can be analyzed. Moreover, as byproducts we obtain new representations of solutions to Schrodinger equations that underwent a confluent SUSY-transformation.
Directivity Modes of Earthquake Populations with Unsupervised Learning

We present a novel approach for resolving modes of rupture directivity in large populations of earthquakes. A seismic spectral decomposition technique is used to first produce relative measurements of radiated energy for earthquakes in a spatially-compact cluster. The azimuthal distribution of energy for each earthquake is then assumed to result from one of several distinct modes of rupture propagation. Rather than fitting a kinematic rupture model to determine the most likely mode of rupture propagation, we instead treat the modes as latent variables and learn them with a Gaussian mixture model. The mixture model simultaneously determines the number of events that best identify with each mode. The technique is demonstrated on four datasets in California with several thousand earthquakes. We show that the datasets naturally decompose into distinct rupture propagation modes that correspond to different rupture directions, and the fault plane is unambiguously identified for all cases. We find that these small earthquakes exhibit unilateral ruptures 53-74% of the time on average. The results provide important observational constraints on the physics of earthquakes and faults.
Orientation, flow, and clogging in a two-dimensional hopper: Ellipses vs. disks

Two-dimensional (2D) hopper flow of disks has been extensively studied. Here, we investigate hopper flow of ellipses with aspect ratio , and we contrast that behavior to the flow of disks. We use a quasi-2D hopper containing photoelastic particles to obtain stress/force information. We simultaneously measure the particle motion and stress. We determine several properties, including discharge rates, clogging probabilities, and the number of particles in clogging arches. For both particle types, the size of the opening, D, relative to the size of particles, l, is an important dimensionless measure. The orientation of the ellipses plays an important role in flow rheology and clogging. The alignment of contacting ellipses enhances the probability of forming stable arches. This study offers insights into applications involving the flow of granular materials consisting of ellipsoidal shapes, and possibly other non-spherical shapes.
What broad emission lines tell us about how active galactic nuclei work

Abstract I review progress made in understanding the nature of the broad-line region (BLR) of active galactic nuclei (AGNs) and the role BLRs play in the AGN phenomenon. The high equivalent widths of the lines imply a high BLR covering factor, and the absence of clear evidence for absorption by the BLR means that the BLR has a flattened distribution and that we always view it near pole-on. The BLR gas is strongly self-shielding near the equatorial plane. Velocity-resolved reverberation mapping has long strongly excluded significant outflow of the BLR and shows instead that the predominant motions are Keplerian with large turbulence and a significant net inflow. The rotation and turbulence are consistent with the inferred geometry. The blueshifting of high-ionization lines is a consequence of scattering off inflowing material rather than the result of an outflowing wind. The rate of inflow of the BLR is sufficient to provide the accretion rate needed to power the AGN. Because the motions of the BLR are gravitationally dominated, and the BLR structure is very similar in most AGNs, consistent black hole masses can be determined. The good correlation between these estimates and masses predicted from the bulge luminosities of host galaxies provides strong support for the similarity of AGN continuum shapes and the correctness of the BLR picture presented. It is concluded that although many mysteries remain about the details of how AGNs work, a general overall picture of the torus and BLR is becoming clear.
Exponential Decay for the Asymptotic Geometry of the Hitchin Metric

We consider Hitchin’s hyperkahler metric $$g_{L^2}$$ on the SU(n)-Hitchin moduli space over a compact Riemann surface. We prove that the difference between the metric $$g_{L^2}$$ and a simpler “semiflat” hyperkahler metric $$g_{\mathrm {sf}}$$ is exponentially-decaying along generic rays in the Hitchin moduli space, as conjectured by Gaiotto–Moore–Neitzke.
Secure Relaying in Non-Orthogonal Multiple Access: Trusted and Untrusted Scenarios

A downlink single-input single-output non-orthogonal multiple access setting is considered, in which a base station (BS) is communicating with two legitimate users in two possible scenarios of unsecure environments: existence of an external eavesdropper and communicating through an untrusted relay. For the first scenario, a number of trusted cooperative half-duplex relays is employed to assist with the BS’s transmission and secure its signals from the external eavesdropper. Various relaying schemes are proposed and analyzed for that matter: cooperative jamming , decode-and-forward , and amplify-and-forward . For each scheme, secure beamforming signals are devised at the relays to maximize the achievable secrecy rate regions. For the second scenario, with the untrusted relay, achievable secrecy rate regions are derived for two different relaying schemes, compress-and-forward and amplify-and-forward , under two different modes of operation. In the first mode, coined passive user mode, the users receive signals from both the BS and the untrusted relay and combine them to decode their messages. In the second mode, termed the active user mode, the users transmit a cooperative jamming signal simultaneously with the BS’s transmission to further confuse the relay. Focusing on half-duplex nodes, the users cannot receive the BS’s signal while jamming the relay, i.e., while being active, and rely only on the signals forwarded to them by the relay. It is shown that the best relaying scheme highly depends on the system parameters, in particular the distances between the nodes, and also on the part of the secrecy rate region at which the system is to operate.
Computational study of flow dynamics from a dc arc plasma jet

Plasma jets produced by direct-current (dc) non-transferred arc plasma torches, at the core of technologies ranging from spray coating to pyrolysis, present intricate dynamics due to the coupled interaction of fluid flow, thermal and electromagnetic phenomena. The flow dynamics from an arc discharge plasma jet are investigated using time-dependent three-dimensional simulations encompassing the dynamics of the arc inside the torch, the evolution of the jet through the discharge environment, and the subsequent impingement of the jet over a flat substrate. The plasma is described by a chemical equilibrium and thermodynamic nonequilibrium (two-temperature) model. The numerical formulation of the physical model is based on a monolithic and fully coupled treatment of the fluid and electromagnetic equations using a variational multiscale finite element method. Simulation results uncover distinct aspects of the flow dynamics, including the jet forcing due to the movement of the electric arc, the prevalence of deviations between heavy-species and electron temperatures in the plasma fringes, the development of shear-flow instabilities around the jet, the occurrence of localized regions with high electric fields far from the arc, and the formation and evolution of coherent flow structures.
Mechanics-Based Solution Verification for Porous Media Models

This paper presents a new approach to verify the accuracy of computational simulations. We develop mathematical theorems which can serve as robust a posteriori error estimation techniques to identify numerical pollution, check the performance of adaptive meshes, and verify numerical solutions. We demonstrate performance of this methodology on problems from flow thorough porous media. However, one can extend it to other models. We construct mathematical properties such that the solutions to Darcy and Darcy-Brinkman equations satisfy them. The mathematical properties include the total minimum mechanical power, minimum dissipation theorem, reciprocal relation, and maximum principle for the vorticity. All the developed theorems have firm mechanical bases and are independent of numerical methods. So, these can be utilized for solution verification of finite element, finite volume, finite difference, lattice Boltzmann methods and so forth. In particular, we show that, for a given set of boundary conditions, Darcy velocity has the minimum total mechanical power of all the kinematically admissible vector fields. We also show that a similar result holds for Darcy-Brinkman velocity. We then show for a conservative body force, the Darcy and Darcy-Brinkman velocities have the minimum total dissipation among their respective kinematically admissible vector fields. Using numerical examples, we show that the minimum dissipation and total mechanical power theorems can be utilized to identify pollution errors in numerical solutions. The solutions to Darcy and Darcy-Brinkman equations are shown to satisfy a reciprocal relation, which has the potential to identify errors in the numerical implementation of boundary conditions. It is also shown that the vorticity under both steady and transient Darcy-Brinkman equations satisfy maximum principles if the body force is conservative and the permeability is homogeneous and isotropic. A discussion on the nature of vorticity under steady and transient Darcy equations is also presented. Using several numerical examples, we will demonstrate the predictive capabilities of the proposed a posteriori techniques in assessing the accuracy of numerical solutions for a general class of problems, which could involve complex domains and general computational grids.
Finite temperature resolution of the Klebanov–Tseytlin singularity

Naked singularities in the gravitational backgrounds dual to gauge theories can be hidden behind the black hole horizon. We present an exact black hole solution in the Klebanov–Tseytlin geometry [I.R. Klebanov, A.A. Tseytlin, hep-th/0002159]. Our solution realizes Maldacena dual of the finite temperature N=1 duality cascade of [I.R. Klebanov, M.J. Strassler, hep-th/0007191] above the temperature of the chiral symmetry breaking. We compare the Bekenstein–Hawking entropy of the black hole with the entropy of the SU(N+M)×SU(N) gauge theory.
Results from EDGES High-band. I. Constraints on Phenomenological Models for the Global 21 cm Signal

We report constraints on the global $21$ cm signal due to neutral hydrogen at redshifts $14.8 \geq z \geq 6.5$. We derive our constraints from low foreground observations of the average sky brightness spectrum conducted with the EDGES High-Band instrument between September $7$ and October $26$, $2015$. Observations were calibrated by accounting for the effects of antenna beam chromaticity, antenna and ground losses, signal reflections, and receiver parameters. We evaluate the consistency between the spectrum and phenomenological models for the global $21$ cm signal. For tanh-based representations of the ionization history during the epoch of reionization, we rule out, at $\geq2\sigma$ significance, models with duration of up to $\Delta z = 1$ at $z\approx8.5$ and higher than $\Delta z = 0.4$ across most of the observed redshift range under the usual assumption that the $21$ cm spin temperature is much larger than the temperature of the cosmic microwave background (CMB) during reionization. We also investigate a `cold' IGM scenario that assumes perfect Ly$\alpha$ coupling of the $21$ cm spin temperature to the temperature of the intergalactic medium (IGM), but that the IGM is not heated by early stars or stellar remants. Under this assumption, we reject tanh-based reionization models of duration $\Delta z \lesssim 2$ over most of the observed redshift range. Finally, we explore and reject a broad range of Gaussian models for the $21$ cm absorption feature expected in the First Light era. As an example, we reject $100$ mK Gaussians with duration (full width at half maximum) $\Delta z \leq 4$ over the range $14.2\geq z\geq 6.5$ at $\geq2\sigma$ significance.
Prodirect manipulation: bidirectional programming for the masses

Software interfaces today generally fall at either end of a spectrum. On one end are programmable systems, which allow expert users (i.e. programmers) to write software artifacts that describe complex abstractions, but programs are disconnected from their eventual output. On the other end are domain-specific graphical user interfaces (GUIs), which allow end users (i.e. non-programmers) to easily create varied content but present insurmountable walls when a desired feature is not built-in. Both programmatic and direct manipulation have distinct strengths, but users must typically choose one over the other or use some ad-hoc combination of systems. Our goal, put simply, is to bridge this divide. We envision novel software systems that tightly couple programmatic and direct manipulation --- a combination we dub prodirect manipulation --- for a variety of use cases. This will require advances in a broad range of software engineering disciplines, from program analysis and program synthesis technology to user interface design and evaluation. In this extended abstract, we propose two general strategies --- real-time program synthesis and domain-specific synthesis of general-purpose programs --- that may prove fruitful for overcoming the technical challenges. We also discuss metrics that will be important in evaluating the usability and utility of prodirect manipulation systems.
Quantifying the Annular Mode Dynamics in an Idealized Atmosphere

AbstractThe linear response function (LRF) of an idealized GCM, the dry dynamical core with Held–Suarez physics, is used to accurately compute how eddy momentum and heat fluxes change in response t...
Night Vision for Small Telescopes

We explore the feasibility of using current generation, off-the-shelf, indium gallium arsenide (InGaAs) near-infrared (NIR) detectors for astronomical observations. Light-weight InGaAs cameras, developed for the night vision industry and operated at or near room temperature, enable cost-effective new paths for observing the NIR sky, particularly when paired with small telescopes. We have tested an InGaAs camera in the laboratory and on the sky using 12 and 18-inch telescopes. The camera is a small-format, 320x240 pixels of 40$\mu$m pitch, Short Wave Infra-Red (SWIR) device from Sensors Unlimited. Although the device exhibits a room-temperature dark current of $5.7 \times 10^4$ $e^-s^{-1}$ per pixel, we find observations of bright sources and low-positional-resolution observations of faint sources remain feasible. We can record unsaturated images of bright ($J=3.9$) sources due to the large pixel well-depth and resulting high dynamic range. When mounted on an 18-inch telescope, the sensor is capable of achieving milli-magnitude precision for sources brighter than $J=8$. Faint sources can be sky-background-limited with modest thermoelectric cooling. We can detect faint sources ($J=16.4$ at $10\sigma$) in a one-minute exposure when mounted to an 18-inch telescope. From laboratory testing, we characterize the noise properties, sensitivity, and stability of the camera in a variety of different operational modes and at different operating temperatures. Through sky testing, we show that the (unfiltered) camera can enable precise and accurate photometry, operating like a filtered $J$-band detector, with small color corrections. In the course of our sky testing, we successfully measured sub-percent flux variations in an exoplanet transit. We have demonstrated an ability to detect transient sources in dense fields using image subtraction of existing reference catalogs.
Resonant Capture by Inward-migrating Planets

We investigate resonant capture of small bodies by planets that migrate inward, using analytic arguments and three-body integrations. If the orbits of the planet and the small body are initially circular and coplanar, the small body is captured when it crosses the 2 : 1 resonance with the planet. As the orbit shrinks, it becomes more eccentric until, by the time its semimajor axis has shrunk by a factor of 4, its eccentricity reaches nearly unity (1 - e 10-4). In typical planetary systems, bodies in this high-eccentricity phase are likely to be consumed by the central star. If they can avoid this fate, as migration continues the inclination flips from 0° to i = 180°; thereafter the eccentricity declines until the semimajor axis is a factor of 9 smaller than at capture, at which point the small body is released from the 2 : 1 resonance on a nearly circular retrograde orbit. Small bodies captured into resonance from initially inclined or eccentric orbits can also be ejected from the system, or released from the resonance on highly eccentric polar orbits (i 90°) that are stabilized by a secular resonance. We conclude that migration could drive much of the inner planetesimal disk into the star, and that postmigration multiplanet systems may not be coplanar.
ON THE STAR FORMATION LAW FOR SPIRAL AND IRREGULAR GALAXIES

A dynamical model for star formation on a galactic scale is proposed in which the interstellar medium is constantly condensing to star-forming clouds on the dynamical time of the average midplane density, and the clouds are constantly being disrupted on the dynamical timescale appropriate for their higher density. In this model, the areal star formation rate scales with the 1.5 power of the total gas column density throughout the main regions of spiral galaxies, and with a steeper power, 2, in the far outer regions and in dwarf irregular galaxies because of the flaring disks. At the same time, there is a molecular star formation law that is linear in the main and outer parts of disks and in dIrrs because the duration of individual structures in the molecular phase is also the dynamical timescale, canceling the additional 0.5 power of surface density. The total gas consumption time scales directly with the midplane dynamical time, quenching star formation in the inner regions if there is no accretion, and sustaining star formation for ~100 Gyr or more in the outer regions with no qualitative change in gas stability or molecular cloud properties. The ULIRG track follows from high densities in galaxy collisions.
Boltzmann-Gibbs entropy is sufficient but not necessary for the likelihood factorization required by Einstein

In 1910 Einstein published a work on a crucial aspect of his understanding of the Boltzmann entropy. He essentially argued that the likelihood function of any system composed by two probabilistically independent subsystems ought to be factorizable into the likelihood functions of each of the subsystems. Consistently he was satisfied by the fact that the Boltzmann (additive) entropy fulfills this epistemologically fundamental requirement. We show here that entropies (e.g., the q-entropy on which nonextensive statistical mechanics is based) which generalize the BG one through violation of its well-known additivity can also fulfill the same requirement. This important fact sheds light on the very foundations of the connection between the micro- and macro-scopic worlds, and consistently supports that the classical thermodynamical Legendre structure is more powerful than the role to it reserved by the Boltzmann-Gibbs statistical mechanics.
Photometric Orbits of Extrasolar Planets

We define and analyze the photometric orbit (PhO) of an extrasolar planet observed in reflected light. In our definition, the PhO is a Keplerian entity with six parameters: semimajor axis, eccentricity, mean anomaly at some particular time, argument of periastron, inclination angle, and effective radius, which is the square root of the geometric albedo times the planetary radius. Preliminarily, we assume a Lambertian phase function. We study in detail the case of short-period giant planets (SPGPs) and observational parameters relevant to the Kepler mission: 20 ppm photometry with normal errors, 6.5 hr cadence, and three-year duration. We define a relevant "planetary population of interest" in terms of probability distributions of the PhO parameters. We perform Monte Carlo experiments to estimate the ability to detect planets and to recover PhO parameters from light curves. We calibrate the completeness of a periodogram search technique, and find structure caused by degeneracy. We recover full orbital solutions from synthetic Kepler data sets and estimate the median errors in recovered PhO parameters. We treat in depth a case of a Jupiter body-double. For the stated assumptions, we find that Kepler should obtain orbital solutions for many of the 100-760 SPGP that Jenkins & Doyle estimate Kepler will discover. Because most or all of these discoveries will be followed up by ground-based radial velocity observations, the estimates of inclination angle from the PhO may enable the calculation of true companion masses: Kepler photometry may break the "msin i" degeneracy. PhO observations may be difficult. There is uncertainty about how low the albedos of SPGPs actually are, about their phase functions, and about a possible noise floor due to systematic errors from instrumental and stellar sources. Nevertheless, simple detection of SPGPs in reflected light should be robust in the regime of Kepler photometry, and estimates of all six orbital parameters may be feasible in at least a subset of cases.
Vortex guidance and transport in channeled pinning arrays

We numerically examine vortices interacting with pinning arrays where a portion of the pinning sites have been removed in order to create coexisting regions of strong and weak pinning. The region w...
Letter Features as Predictors of Letter-Name Acquisition in Four Languages with Three Scripts

To expand our understanding of script-general and script-specific principles in the learning of letter names, we examined how three characteristics of alphabet letters – their frequency in printed ...
Low-temperature series for renormalized operators: The ferromagnetic square-lattice Ising model

A method for computing low-temperature series for renormalized operators in the two-dimensional Ising model is proposed. These series are applied to the study of the properties of the truncated renormalized Hamiltonians when we start at very low temperature and zero field. The truncated Hamiltonians for majority rule, Kadanoff transformation, and decimation for 2×2 blocks depend on the how we approach the first-order phase-transition line. The renormalization group transformations are multivalued and discontinuous at this first-order transition line when restricted to some finite-dimensional interaction space.
On the martingale problem for degenerate-parabolic partial differential operators with unbounded coefficients and a mimicking theorem for Itô processes

Using results from our companion article [26] on a Schauder approach to existence of solutions to a degenerate-parabolic partial differential equation, we solve three intertwined problems, motivated by probability theory and mathematical finance, concerning degenerate diffusion processes. We show that the martingale problem associated with a degenerate-elliptic differential operator with unbounded, locally Holder continuous coefficients on a half-space is wellposed in the sense of Stroock and Varadhan. Second, we prove existence, uniqueness, and the strong Markov property for weak solutions to a stochastic differential equation with degenerate diffusion and unbounded coefficients with suitable Holder continuity properties. Third, for an Ito process with degenerate diffusion and unbounded but appropriately regular coefficients, we prove existence of a strong Markov process, unique in the sense of probability law, whose onedimensional marginal probability distributions match those of the given Ito process.
Exponential decay in the mapping class group

We show that the probability that a finitely supported random walk on a non-elementary subgroup of the the mapping class group gives a non-pseudo-Anosov element decays exponentially in the length of the random walk. More generally, we show that if R is a set of mapping class group elements with an upper bound on their translation lengths on the complex of curves, then the probability that a random walk lies in R decays exponentially in the length of the random walk.
A Constrained Variable Projection Reconstruction Method for Photoacoustic Computed Tomography Without Accurate Knowledge of Transducer Responses

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the absorbed optical energy density within tissue. When the imaging system employs conventional piezoelectric ultrasonic transducers, the ideal photoacoustic (PA) signals are degraded by the transducers' acousto-electric impulse responses (EIRs) during the measurement process. If unaccounted for, this can degrade the accuracy of the reconstructed image. In principle, the effect of the EIRs on the measured PA signals can be ameliorated via deconvolution; images can be reconstructed subsequently by application of a reconstruction method that assumes an idealized EIR. Alternatively, the effect of the EIR can be incorporated into an imaging model and implicitly compensated for during reconstruction. In either case, the efficacy of the correction can be limited by errors in the assumed EIRs. In this work, a joint optimization approach to PACT image reconstruction is proposed for mitigating errors in reconstructed images that are caused by use of an inaccurate EIR. The method exploits the bi-linear nature of the imaging model and seeks to refine the measured EIR during the process of reconstructing the sought-after absorbed optical energy density. Computer-simulation and experimental studies are conducted to investigate the numerical properties of the method and demonstrate its value for mitigating image distortions and enhancing the visibility of fine structures.
Diagnosing added value of convection-permitting regional models using precipitation event identification and tracking

Dynamical downscaling with high-resolution regional climate models may offer the possibility of realistically reproducing precipitation and weather events in climate simulations. As resolutions fall to order kilometers, the use of explicit rather than parametrized convection may offer even greater fidelity. However, these increased resolutions both allow and require increasingly complex diagnostics for evaluating model fidelity. In this study we focus on precipitation evaluation and analyze five 2-month-long dynamically downscaled model runs over the continental United States that employ different convective and microphysics parameterizations, including one high-resolution convection-permitting simulation. All model runs use the Weather Research and Forecasting Model driven by National Center for Environmental Prediction reanalysis data. We show that employing a novel rainstorm identification and tracking algorithm that allocates essentially all rainfall to individual precipitation events (Chang et al. in J Clim 29(23):8355–8376, 2016) allows new insights into model biases. Results include that, at least in these runs, model wet bias is driven by excessive areal extent of individual precipitating events, and that the effect is time-dependent, producing excessive diurnal cycle amplitude. This amplified cycle is driven not by new production of events but by excessive daytime enlargement of long-lived precipitation events. We further show that in the domain average, precipitation biases appear best represented as additive offsets. Of all model configurations evaluated, convection-permitting simulations most consistently reduced biases in precipitation event characteristics.
The generic quantum superintegrable system on the sphere and Racah operators

We consider the generic quantum superintegrable system on the d-sphere with potential $$V(y)=\sum _{k=1}^{d+1}\frac{b_k}{y_k^2}$$ , where $$b_k$$ are parameters. Appropriately normalized, the symmetry operators for the Hamiltonian define a representation of the Kohno–Drinfeld Lie algebra on the space of polynomials orthogonal with respect to the Dirichlet distribution. The Gaudin subalgebras generated by Jucys–Murphy elements are diagonalized by families of Jacobi polynomials in d variables on the simplex. We define a set of generators for the symmetry algebra, and we prove that their action on the Jacobi polynomials is represented by the multivariable Racah operators introduced in Geronimo and Iliev (Constr Approx 31(3):417–457, 2010). The constructions also yield a new Lie-theoretic interpretation of the bispectral property for Tratnik’s multivariable Racah polynomials.
Web-Based VR Experiments Powered by the Crowd

We build on the increasing availability of Virtual Reality (VR) devices and Web technologies to conduct behavioral experiments in VR using crowdsourcing techniques. A new recruiting and validation method allows us to create a panel of eligible experiment participants recruited from Amazon Mechanical Turk. Using this panel, we ran three different crowdsourced VR experiments, each reproducing one of three VR illusions: place illusion, embodiment illusion, and plausibility illusion. Our experience and worker feedback on these experiments show that conducting Web-based VR experiments using crowdsourcing is already feasible, though some challenges---including scale---remain. Such crowdsourced VR experiments on the Web have the potential to finally support replicable VR experiments with diverse populations at a low cost.
FIRST SIMULTANEOUS DETECTION OF MOVING MAGNETIC FEATURES IN PHOTOSPHERIC INTENSITY AND MAGNETIC FIELD DATA

The formation and the temporal evolution of a bipolar moving magnetic feature (MMF) was studied with high-spatial and temporal resolution. The photometric properties were observed with the New Solar Telescope at Big Bear Solar Observatory using a broadband TiO filter (705.7 nm), while the magnetic field was analyzed using the spectropolarimetric data obtained by Hinode. For the first time, we observed a bipolar MMF simultaneously in intensity images and magnetic field data, and studied the details of its structure. The vector magnetic field and the Doppler velocity of the MMF were also studied. A bipolar MMF with its positive polarity closer to the negative penumbra formed, accompanied by a bright, filamentary structure in the TiO data connecting the MMF and a dark penumbral filament. A fast downflow (≤2 km s–1) was detected at the positive polarity. The vector magnetic field obtained from the full Stokes inversion revealed that a bipolar MMF has a U-shaped magnetic field configuration. Our observations provide a clear intensity counterpart of the observed MMF in the photosphere, and strong evidence of the connection between the MMF and the penumbral filament as a serpentine field.