cssi/SciSciGPT-SciSciNet huggingface dataset:
Classifying antimicrobial and multifunctional peptides with Bayesian network models

Bayesian network models are finding success in characterizing enzyme-catalyzed reactions, slow conformational changes, predicting enzyme inhibition, and genomics. In this work, we apply them to statistical modeling of peptides by simultaneously identifying amino acid sequence motifs and using a motif-based model to clarify the role motifs may play in antimicrobial activity. We construct models of increasing sophistication, demonstrating how chemical knowledge of a peptide system may be embedded without requiring new derivation of model fitting equations after changing model structure. These models are used to construct classifiers with good performance (94% accuracy, Matthews correlation coefficient of 0.87) at predicting antimicrobial activity in peptides, while at the same time being built of interpretable parameters. We demonstrate use of these models to identify peptides that are potentially both antimicrobial and antifouling, and show that the background distribution of amino acids could play a greater role in activity than sequence motifs do. This provides an advancement in the type of peptide activity modeling that can be done and the ease in which models can be constructed.
Dispersive estimates for the wave equation on Riemannian manifolds of bounded curvature

We prove space-time dispersive estimates for solutions to the wave equation on compact Riemannian manifolds with bounded curvature tensor, where we assume that the metric tensor is of W 1 , p regularity for some p > d , which ensures that the curvature tensor is well-defined in the weak sense. The estimates are established for the same range of Lebesgue and Sobolev exponents that hold in the case of smooth metrics. Our results are for bounded time intervals, so by finite propagation velocity they hold also on noncompact manifolds under appropriate uniform geometry conditions.
Brunnian links are determined by their complements.

If L 1 and L 2 are two Brunnian links with all pairwise linking numbers 0, then we show that L 1 and L 2 are equivalent if and only if they have homeomorphic complements. In particular, this holds for all Brunnian links with at least three components. If L 1 is a Brunnian link with all pairwise linking numbers 0, and the complement of L 2 is homeomorphic to the complement of L 1 , then we show that L 2 may be obtained from L 1 by a sequence of twists around unknotted components. Finally, we show that for any positive integer n, an algorithm for detecting an n-component unlink leads immediately to an algorithm for detecting an unlink of any number of components. This algorithmic generalization is conceptually simple, but probably computationally impractical.
SUZAKU OBSERVATIONS OF THE X-RAY BRIGHTEST FOSSIL GROUP ESO 3060170

'Fossil' galaxy groups, each dominated by a relatively isolated giant elliptical galaxy, have many properties intermediate between groups and clusters of galaxies. We used the Suzaku X-ray observatory to observe the X-ray brightest fossil group, ESO 3060170, out to R{sub 200}, in order to better elucidate the relation between fossil groups, normal groups, and clusters. We determined the intragroup gas temperature, density, and metal abundance distributions and derived the entropy, pressure, and mass profiles for this group. The entropy and pressure profiles in the outer regions are flatter than in simulated clusters, similar to what is seen in observations of massive clusters. This may indicate that the gas is clumpy and/or the gas has been redistributed. Assuming hydrostatic equilibrium, the total mass is estimated to be ∼1.7 × 10{sup 14} M{sub ☉} within a radius R{sub 200} of ∼1.15 Mpc, with an enclosed baryon mass fraction of 0.13. The integrated iron mass-to-light ratio of this fossil group is larger than in most groups and comparable to those of clusters, indicating that this fossil group has retained the bulk of its metals. A galaxy luminosity density map on a scale of 25 Mpc shows that this fossil group resides in a more » relatively isolated environment, unlike the filamentary structures in which typical groups and clusters are embedded. « less
HOW SPIRALS AND GAPS DRIVEN BY COMPANIONS IN PROTOPLANETARY DISKS APPEAR IN SCATTERED LIGHT AT ARBITRARY VIEWING ANGLES

Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.
AN ENHANCED RATE OF TIDAL DISRUPTIONS IN THE CENTRALLY OVERDENSE E+A GALAXY NGC 3156

Time domain optical surveys have discovered roughly a dozen candidate stellar tidal disruption flares in the last five years, and future surveys like the {\it Large Synoptic Survey Telescope} will likely find hundreds to thousands more. These tidal disruption events (TDEs) present an interesting puzzle: a majority of the current TDE sample is hosted by rare post-starburst galaxies, and tens of percent are hosted in even rarer E+A galaxies, which make up $\sim 0.1\%$ of all galaxies in the local universe. E+As are therefore overrepresented among TDE hosts by 1-2 orders of magnitude, a discrepancy unlikely to be accounted for by selection effects. We analyze {\it Hubble Space Telescope} photometry of one of the nearest E+A galaxies, NGC~3156, to estimate the rate of stellar tidal disruption produced as two-body relaxation diffuses stars onto orbits in the loss cone of the central supermassive black hole. The rate of TDEs produced by two-body relaxation in NGC~3156 is large when compared to other galaxies with similar black hole mass: $\dot{N}_{\rm TDE}\sim 1\times 10^{-3}~{\rm yr}^{-1}$. This suggests that the preference of TDEs for E+A hosts may be due to central stellar overdensities produced in recent starbursts.
Holographic entanglement entropy and Fermi surfaces

The entanglement entropy in theories with a Fermi surface is known to pro- duce a logarithmic violation of the usual area law behavior. We explore the possibility of producing this logarithmic violation holographically by analyzing the IR regions of the bulk geometries dual to such theories. The geometry of Ogawa, Takayanagi, and Ugajin is explored and shown to have a null curvature singularity for all values of parameters, except for dynamical critical exponent z = 3=2 (the hyperscaling violation parameter = 1 for this geometry) in four dimensions. The results are extended to general hyperscaling viola- tion exponent. We explore strings propagating through the singularity and show that they become innitely excited, suggesting the singularity is not resolved by stringy eects and
Shilnikov Homoclinic Bifurcation of Mixed-Mode Oscillations

The Koper model is a three-dimensional vector field that was developed to study complex electrochemical oscillations arising in a diffusion process. Koper and Gaspard described paradoxical dynamics in the model: they discovered complicated, chaotic behavior consistent with a homoclinic orbit of Shilnikov type but were unable to locate the orbit itself. The Koper model has since served as a prototype to study the emergence of mixed-mode oscillations (MMOs) in slow-fast systems, but only in this paper is the existence of these elusive homoclinic orbits established. They are found first in a larger family that has been used to study singular Hopf bifurcation in multiple time scale systems with two slow variables and one fast variable. A curve of parameters with homoclinic orbits in this larger family is obtained by continuation and shown to cross the submanifold of the Koper system. The strategy used to compute the homoclinic orbits is based upon systematic investigation of intersections of invariant manifol...
Lessons learned from (and since) the Voyager 2 flybys of Uranus and Neptune.

More than 30 years have passed since the Voyager 2 flybys of Uranus and Neptune. This paper outlines a range of lessons learned from Voyager, broadly grouped into process, planning and people. In t...
Star-shaped oscillations of Leidenfrost drops

An experimental investigation of oscillations of Leidenfrost drops finds that the oscillations are driven by capillary waves beneath the drop.
Less is more: Sampling chemical space with active learning

The development of accurate and transferable machine learning (ML) potentials for predicting molecular energetics is a challenging task. The process of data generation to train such ML potentials is a task neither well understood nor researched in detail. In this work, we present a fully automated approach for the generation of datasets with the intent of training universal ML potentials. It is based on the concept of active learning (AL) via Query by Committee (QBC), which uses the disagreement between an ensemble of ML potentials to infer the reliability of the ensemble's prediction. QBC allows the presented AL algorithm to automatically sample regions of chemical space where the ML potential fails to accurately predict the potential energy. AL improves the overall fitness of ANAKIN-ME (ANI) deep learning potentials in rigorous test cases by mitigating human biases in deciding what new training data to use. AL also reduces the training set size to a fraction of the data required when using naive random sampling techniques. To provide validation of our AL approach, we develop the COmprehensive Machine-learning Potential (COMP6) benchmark (publicly available on GitHub) which contains a diverse set of organic molecules. Active learning-based ANI potentials outperform the original random sampled ANI-1 potential with only 10% of the data, while the final active learning-based model vastly outperforms ANI-1 on the COMP6 benchmark after training to only 25% of the data. Finally, we show that our proposed AL technique develops a universal ANI potential (ANI-1x) that provides accurate energy and force predictions on the entire COMP6 benchmark. This universal ML potential achieves a level of accuracy on par with the best ML potentials for single molecules or materials, while remaining applicable to the general class of organic molecules composed of the elements CHNO.
Model Realization and Numerical Studies of a Three-Dimensional Bosonic Topological Insulator and Symmetry-Enriched Topological Phases

We study a topological phase of interacting bosons in (3 + 1) dimensions that is protected by charge conservation and time-reversal symmetry. We present an explicit lattice model that realizes this phase and that can be studied in sign-free Monte Carlo simulations. The idea behind our model is to bind bosons to topological defects called hedgehogs. We determine the phase diagram of the model and identify a phase where such bound states are proliferated. In this phase, we observe a Witten effect in the bulk whereby an external monopole binds half of the elementary boson charge, which confirms that it is a bosonic topological insulator. We also study the boundary between the topological insulator and a trivial insulator. We find a surface phase diagram that includes exotic superfluids, a topologically ordered phase, and a phase with a Hall effect quantized to one-half of the value possible in a purely two-dimensional system. We also present models that realize symmetry-enriched topologically ordered phases by binding multiple hedgehogs to each boson; these phases show charge fractionalization and intrinsic topological order as well as a fractional Witten effect.
OBSCURATION BY GAS AND DUST IN LUMINOUS QUASARS

We explore the connection between absorption by neutral gas and extinction by dust in mid-infrared (IR) selected luminous quasars. We use a sample of 33 quasars at redshifts 0.7 150 counts) in Chandra observations. We divide the quasars into dust-obscured and unobscured samples based on their optical to mid-IR color, and measure the neutral hydrogen column density N {sub H} through fitting of the X-ray spectra. We find that all subsets of quasars have consistent power law photon indices Γ ≈ 1.9 that are uncorrelated with N {sub H}. We classify the quasars as gas-absorbed or gas-unabsorbed if N {sub H} > 10{sup 22} cm{sup –2} or N {sub H} < 10{sup 22} cm{sup –2}, respectively. Of 24 dust-unobscured quasars in the sample, only one shows clear evidence for significant intrinsic N {sub H}, while 22 have column densities consistent with N {sub H} < 10{sup 22} cm{sup –2}. In contrast, of the nine dust-obscured quasars, six show evidence for intrinsic gas absorption, and three are consistent with Nmore » {sub H} < 10{sup 22} cm{sup –2}. We conclude that dust extinction in IR-selected quasars is strongly correlated with significant gas absorption as determined through X-ray spectral fitting. These results suggest that obscuring gas and dust in quasars are generally co-spatial, and confirm the reliability of simple mid-IR and optical photometric techniques for separating quasars based on obscuration.« less
Oort groups and lifting problems

Let k be an algebraically closed field of positive characteristic p. We consider which finite groups G have the property that every faithful action of G on a connected smooth projective curve over k lifts to characteristic zero. Oort conjectured that cyclic groups have this property. We show that if a cyclic-by-p group G has this property, then G must be either cyclic or dihedral, with the exception of A4 in characteristic two. This proves one direction of a strong form of the Oort conjecture.
Progenitor Mass Distribution for Core-collapse Supernova Remnants in M31 and M33

Using the star formation histories (SFHs) near 94 supernova remnants (SNRs), we infer the progenitor mass distribution for core-collapse supernovae. We use Bayesian inference and model each SFH with multiple bursts of star formation (SF), one of which is assumed to be associated with the SNR. Assuming single-star evolution, the minimum mass of CCSNe is $7.33^{+0.02}_{-0.16}$ $\text{M}_\odot$, the slope of the progenitor mass distribution is $\alpha = -2.96^{+0.45}_{-0.25}$, and the maximum mass is greater than $\text{M}_\textrm{max} > 59$ $\text{M}_\odot$ with a 68% confidence. While these results are consistent with previous inferences, they also provide tighter constraints. The progenitor distribution is somewhat steeper than a Salpeter initial mass function ($\alpha$ = -2.35). This suggests that either SNR catalogs are biased against the youngest SF regions, or the most massive stars do not explode as easily as lower mass stars. If SNR catalogs are biased, it will most likely affect the slope but not the minimum mass. The uncertainties are dominated by three primary sources of uncertainty, the SFH resolution, the number of SF bursts, and the uncertainty on SF rate in each age bin. We address the first two of these uncertainties, with an emphasis on multiple bursts. The third will be addressed in future work.
Cohomology for Drinfeld doubles of some infinitesimal group schemes

Consider a field k of characteristic p > 0 , the r -th Frobenius kernel G ( r ) of a smooth algebraic group G , the Drinfeld double D G ( r ) of G ( r ) , and a finite dimensional D G ( r ) -module M . We prove that the cohomology algebra H ∗ ( D G ( r ) , k ) is finitely generated and that H ∗ ( D G ( r ) , M ) is a finitely generated module over this cohomology algebra. We exhibit a finite map of algebras θ r : H ∗ ( G ( r ) , k ) ⊗ S ( g ) → H ∗ ( D G ( r ) , k ) , which offers an approach to support varieties for D G ( r ) -modules. For many examples of interest, θ r is injective and induces an isomorphism of associated reduced schemes. For M an irreducible D G ( r ) -module, θ r enables us to identify the support variety of M in terms of the support variety of M viewed as a G ( r ) -module.
Optimal experimental design under irreducible uncertainty for linear inverse problems governed by PDEs

We present a method for computing A-optimal sensor placements for infinite-dimensional Bayesian linear inverse problems governed by PDEs with irreducible model uncertainties. Here, irreducible uncertainties refers to uncertainties in the model that exist in addition to the parameters in the inverse problem, and that cannot be reduced through observations. Specifically, given a statistical distribution for the model uncertainties, we compute the optimal design that minimizes the expected value of the posterior covariance trace. The expected value is discretized using Monte Carlo leading to an objective function consisting of a sum of trace operators and a binary-inducing penalty. Minimization of this objective requires a large number of PDE solves in each step. To make this problem computationally tractable, we construct a composite low-rank basis using a randomized range finder algorithm to eliminate forward and adjoint PDE solves. We also present a novel formulation of the A-optimal design objective that requires the trace of an operator in the observation rather than the parameter space. The binary structure is enforced using a weighted regularized $\ell_0$-sparsification approach. We present numerical results for inference of the initial condition in a subsurface flow problem with inherent uncertainty in the flow fields and in the initial times.
Fields and fluids on curved non-relativistic spacetimes

We consider non-relativistic curved geometries and argue that the background structure should be generalized from that considered in previous works. In this approach the derivative operator is defined by a Galilean spin connection valued in the Lie algebra of the Galilean group. This includes the usual spin connection plus an additional “boost connection” which parameterizes the freedom in the derivative operator not fixed by torsion or metric compatibility. As an example we write down the most general theory of dissipative fluids consistent with the second law in curved non-relativistic geometries and find significant differences in the allowed transport coefficients from those found previously. Kubo formulas for all response coefficients are presented. Our approach also immediately generalizes to systems with independent mass and charge currents as would arise in multicomponent fluids. Along the way we also discuss how to write general locally Galilean invariant non-relativistic actions for multiple particle species at any order in derivatives. A detailed review of the geometry and its relation to non-relativistic limits may be found in a companion paper.
Homology operations and cosimplicial iterated loop spaces

If X is a cosimplical $E_{n+1}$ space then Tot(X) is an $E_{n+1}$ space and its mod 2 homology $H_*(Tot(X))$ has Dyer-Lashof and Browder operations. It's natural to ask if the spectral sequence converging to $H_*(Tot(X))$ admits compatible operations. In this paper I give a positive answer to this question.
A multiscale continuous Galerkin method for stochastic simulation and robust design of photonic crystals

Abstract We present a multiscale continuous Galerkin (MSCG) method for the fast and accurate stochastic simulation and optimization of time-harmonic wave propagation through photonic crystals. The MSCG method exploits repeated patterns in the geometry to drastically decrease computational cost and incorporates the following ingredients: (1) a reference domain formulation that allows us to treat geometric variability resulting from manufacturing uncertainties; (2) a reduced basis approximation to solve the parametrized local subproblems; (3) a gradient computation of the objective function; and (4) a model and variance reduction technique that enables the accelerated computation of statistical outputs by exploiting the statistical correlation between the MSCG solution and the reduced basis approximation. The proposed method is thus well suited for both deterministic and stochastic simulations, as well as robust design of photonic crystals. We provide convergence and cost analysis of the MSCG method, as well as a simulation results for a waveguide T-splitter and a Z-bend to illustrate its advantages for stochastic simulation and robust design.
(Small) Resonant non-Gaussianities: Signatures of a Discrete Shift Symmetry in the Effective Field Theory of Inflation

We apply the Effective Field Theory of Inflation to study the case where the continuous shift symmetry of the Goldstone boson π is softly broken to a discrete subgroup. This case includes and generalizes recently proposed String Theory inspired models of Inflation based on Axion Monodromy. The models we study have the property that the 2-point function oscillates as a function of the wavenumber, leading to oscillations in the CMB power spectrum. The non-linear realization of time diffeomorphisms induces some self-interactions for the Goldstone boson that lead to a peculiar non-Gaussianity whose shape oscillates as a function of the wavenumber. We find that in the regime of validity of the effective theory, the oscillatory signal contained in the n−point correlation functions, with n > 2, is smaller than the one contained in the 2-point function, implying that the signature of oscillations, if ever detected, will be easier to find first in the 2-point function, and only then in the higher order correlation functions. Still the signal contained in higher-order correlation functions, that we study here in generality, could be detected at a subleading level, providing a very compelling consistency check for an approximate discrete shift symmetry being realized during inflation.
Echoes from a Warped Dimension

The Randall-Sundrum (RS) model, based on a slice of warped 5D spacetime, was originally introduced to explain the apparent hierarchy between the scales of weak and gravitational interactions. Over the past decade, this model has been extended to provide a predictive theory of flavor, as well as to address various constraints from precision data. In this talk, we will present a brief review of the RS model and some of its extensions. A survey of some key signatures of realistic warped models and the experimental challenges they pose will be presented. We will also consider truncated warped scenarios that address smaller hierarchies and discuss why their discovery can have improved prospects, at the LHC.
Extinction Map of Baade’s Window

Recently Woźniak & Stanek proposed a new method to investigate interstellar extinction, based on two-band photometry, which uses red clump stars as a means to construct the reddening curve. I apply this method to the color-magnitude diagrams obtained by the Optical Gravitational Lensing Experiment to construct an extinction map of a 40' × 40' region of Baade's window, with resolution of ~30''. Such a map should be useful for studies of this frequently observed region of the Galactic bulge. The map and software useful for its applications are available via anonymous ftp. The total extinction AV varies from 1.26 to 2.79 mag within the 40' × 40' field of view centered on (α2000, δ2000) = (18:03:20.9, -30:02:06), i.e., (l, b) = (1.001, -3.885). The ratio AV/E(V - I) = 2.49 ± 0.02 is determined with this new method.
Over-limiting Current and Control of Dendritic Growth by Surface Conduction in Nanopores

Understanding over-limiting current (faster than diffusion) is a long-standing challenge in electrochemistry with applications in desalination and energy storage. Known mechanisms involve either chemical or hydrodynamic instabilities in unconfined electrolytes. Here, it is shown that over-limiting current can be sustained by surface conduction in nanopores, without any such instabilities, and used to control dendritic growth during electrodeposition. Copper electrodeposits are grown in anodized aluminum oxide membranes with polyelectrolyte coatings to modify the surface charge. At low currents, uniform electroplating occurs, unaffected by surface modification due to thin electric double layers, but the morphology changes dramatically above the limiting current. With negative surface charge, growth is enhanced along the nanopore surfaces, forming surface dendrites and nanotubes behind a deionization shock. With positive surface charge, dendrites avoid the surfaces and are either guided along the nanopore centers or blocked from penetrating the membrane.
Screw-Dislocation-Driven Growth of Two-Dimensional Few-Layer and Pyramid-like WSe2 by Sulfur-Assisted Chemical Vapor Deposition

Two-dimensional (2D) layered tungsten diselenides (WSe2) material has recently drawn a lot of attention due to its unique optoelectronic properties and ambipolar transport behavior. However, direct chemical vapor deposition (CVD) synthesis of 2D WSe2 is not as straightforward as other 2D materials due to the low reactivity between reactants in WSe2 synthesis. In addition, the growth mechanism of WSe2 in such CVD process remains unclear. Here we report the observation of a screw-dislocation-driven (SDD) spiral growth of 2D WSe2 flakes and pyramid-like structures using a sulfur-assisted CVD method. Few-layer and pyramid-like WSe2 flakes instead of monolayer were synthesized by introducing a small amount of sulfur as a reducer to help the selenization of WO3, which is the precursor of tungsten. Clear observations of steps, helical fringes, and herringbone contours under atomic force microscope characterization reveal the existence of screw dislocations in the as-grown WSe2. The generation and propagation mec...
Reconstruction of a compactly supported sound profile in the presence of a random background medium

In this paper, we present algorithms for reconstructing an unknown compact scatterer embedded in a random noisy background medium, given measurements of the scattered field and information about the background medium and the sound profile. We present six different methods for the solution of this inverse problem using different amounts of scattered data and prior information about the random background medium and the scatterer. The different inversion algorithms are defined by a combination of stochastic programming methods and Bayesian formulation. Our basic results show that if we have data for just one instance of the random background medium the best strategy is to invert for both random medium and unknown scatterer with appropriate regularization. However, if we have data for multiple instances of the medium it may be worth solving a coupled set of multiple inverse problems. We present several numerical results for inverting for various scatterer geometries under different inversion scenarios. The main take-away of our study is that one should invert for both unknown scatterer and random medium, with appropriate, prior-information based regularization. Furthermore, if data from multiple realizations of the background medium is available, then it may be beneficial to combine results from multiple inversions.
Regge gravity from spinfoams

We consider spinfoam quantum gravity in the flipped limit, which is the double scaling limit γ → 0, j → ∞ with γj = const., where γ is the Immirzi parameter, j is the spin and γj gives the physical area in Planck units. In this regime the amplitude for a 2-complex becomes effectively an integral over Regge-like metrics and seems to enforce Einstein equations in the semiclassical regime. The Immirzi parameter must be considered as dynamical in the sense that it runs to zero when the fine structure of the foam is averaged. In addition to quantum corrections which vanish for ℏ → 0, we find new corrections due to the discreteness of geometric spectra.
The chiral SYK model

We study the generalization of the Sachdev-Ye-Kitaev (SYK) model to a 1 + 1 dimensional chiral SYK model of N flavors of right-moving chiral Majorana fermions with all-to-all random 4-fermion interactions. The interactions in this model are exactly marginal, leading to an exact scaling symmetry. We show the Schwinger-Dyson equation of this model in the large N limit is exactly solvable. In addition, we show this model is integrable for small N ≤ 6 by bosonization. Surprisingly, the two point function in the large N limit has exactly the same form as that for N = 4, although the four point functions of the two cases are quite different. The ground state entropy in the large N limit is the same as that of N free chiral Majorana fermions, leading to a zero ground state entropy density. The OTOC of the model in the large N limit exhibits a non-trivial spacetime structure reminscent of that found by Gu and Kitaev [1] for generic SYK-like models. Specifically we find a Lyapunov regime inside an asymmetric butterfly cone, which are signatures of quantum chaos, and that the maximal velocity dependent Lyapunov exponent approaches the chaos bound 2π/β as the interaction strength approaches its physical upper bound. Finally, the model is integrable for (at least) N ≤ 6 but chaotic in the large N limit, leading us to conjecture that there is a transition from integrability to chaos as N increases past a critical value.
An Extended Magnetohydrodynamics Model for Relativistic Weakly Collisional Plasmas

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.
Quark and Glue Components of the Proton Spin from Lattice Calculation

The status of lattice calculations of the quark spin, the quark orbital angular momentum, the glue angular momentum and glue spin in the nucleon is summarized. The quark spin calculation is recently carried out from the anomalous Ward identity with chiral fermions and is found to be small mainly due to the large negative anomaly term which is believed to be the source of the ‘proton spin crisis’. We also present the first calculation of the glue spin at finite nucleon momenta.
A continuation of solutions to convolution equations with the loss of smoothness

In the present paper the smoothness loss of a continuation of solutions to convolution equations is studied. Also examples for some kinds of convolvers are given.
K -theory Schubert calculus of the affine Grassmannian

We construct the Schubert basis of the torus-equivariant K-homology of the ane Grassmannian of a simple algebraic group G, using the K-theoretic NilHecke ring of Kostant and Kumar. This is the K-theoretic analogue of a construction of Peterson in equivariant homology. For the case where G = SLn, the K-homology of the ane Grassmannian is identied with a sub-Hopf algebra of the ring of symmetric functions. The Schubert basis is represented by inhomogeneous symmetric functions, called K-k-Schur functions, whose highest-degree term is a k-Schur function. The dual basis in K-cohomology is given by the ane stable Grothendieck polynomials, verifying a conjecture of Lam. In addition, we give a Pieri rule in K-homology. Many of our constructions have geometric interpretations by means of Kashiwara’s thick ane ag manifold.
RENORMALIZATION GROUP TRAJECTORIES BETWEEN TWO FIXED POINTS

We report on our recent rigorous construction of complete renormalization group trajectories between two fixed points for the three-dimensional phi-four model with modified propagator considered by Brydges, Mitter and Scoppola (BMS). These are discrete critical trajectories which connect the ultraviolet Gaussian fixed point to the nontrivial BMS infrared fixed point which is an analogue of the Wilson–Fisher fixed point. The renormalization group map is defined rigorously and nonperturbatively, without using the hierarchical approximation. The trajectories are constructed by a fixed point argument in a suitable Banach space of sequences, where one perturbs a nonlinear one-dimensional iteration.
Stability of Nitrogen in Planetary Atmospheres in Contact with Liquid Water

Molecular nitrogen is the most commonly assumed background gas that supports habitability on rocky planets. Despite its chemical inertness, nitrogen molecule is broken by lightning, hot volcanic vents, and bolide impacts, and can be converted into soluble nitrogen compounds and then sequestered in the ocean. The very stability of nitrogen, and that of nitrogen-based habitability, is thus called into question. Here we determine the lifetime of molecular nitrogen vis-a-vis aqueous sequestration, by developing a novel model that couples atmospheric photochemistry and oceanic chemistry. We find that HNO, the dominant nitrogen compounds produced in anoxic atmospheres, is converted to N2O in the ocean, rather than oxidized to nitrites or nitrates as previously assumed. This N2O is then released back into the atmosphere and quickly converted to N2. We also find that the deposition rate of NO is severely limited by the kinetics of the aqueous-phase reaction that converts NO to nitrites in the ocean. Putting these insights together, we conclude that the atmosphere must produce nitrogen species at least as oxidized as NO2 and HNO2 to enable aqueous sequestration. The lifetime of molecular nitrogen in anoxic atmospheres is determined to be >1 billion years on temperate planets of both Sun-like and M dwarf stars. This result upholds the validity of molecular nitrogen as a universal background gas on rocky planets.
Synchrotron radio emission in radio-quiet AGNs

The basic mechanism responsible for radio emission in radio-loud active galactic nuclei (AGNs) is assumed to be synchrotron radiation. We suggest here that radio emission in radio-quiet objects is also due to synchrotron radiation of particles accelerated in shocks. We consider generic shocks and study the resulting synchrotron properties. We estimate the synchrotron radio luminosity and compare it with the X-ray component produced by inverse Compton emission, and find that the radio-to-X-ray luminosity ratio is much less than unity, with values typical of radio-quiet sources. The predicted trends in source parameters, black hole mass, and accretion rate may account for the anticorrelation between radio loudness and Eddington ratio observed in different AGN samples.
A low-rank coordinate-descent algorithm for semidefinite programming relaxations of optimal power flow

The alternating-current optimal power flow (ACOPF) is one of the best known non-convex nonlinear optimization problems. We present a novel re-formulation of ACOPF, which is based on lifting the rectangular power-voltage rank-constrained formulation, and makes it possible to derive alternative semidefinite programming relaxations. For those, we develop a first-order method based on the parallel coordinate descent with a novel closed-form step based on roots of cubic polynomials.
EXACT UPPER AND LOWER BOUNDS ON THE DIFFERENCE BETWEEN THE ARITHMETIC AND GEOMETRIC MEANS

Exact upper and lower bounds on the difference between the arithmetic and geometric means are obtained. The inequalities providing these bounds may be viewed, respectively, as a reverse Jensen inequality and an improvement of the direct Jensen inequality, in the case when the convex function is the exponential.
Magnetophoretic separation of blood cells at the microscale

We present a method and model for the direct and continuous separation of red and white blood cells in plasma. The method is implemented at the microscale using a microfluidic system that consists of an array of integrated soft-magnetic elements embedded adjacent to a microfluidic channel. The microsystem is passive and is activated via application of a bias field that magnetizes the elements. Once magnetized, the elements produce a nonuniform magnetic field distribution in the microchannel, which gives rise to a force on blood cells as they pass through the microsystem. In whole blood, white blood cells behave as diamagnetic microparticles while red blood cells exhibit diamagnetic or paramagnetic behaviour depending on the oxygenation of their haemoglobin. We develop a mathematical model for predicting the motion of blood cells in the microsystem that takes into account the dominant magnetic, fluidic and buoyant forces on the cells. We use the model to study red/white blood cell transport, and our analysis indicates that the microsystem is capable of rapid and efficient red/white blood cell separation.
Capacity of 1-to-K broadcast packet erasure channels with channel output feedback

This paper focuses on the 1-to-K broadcast packet erasure channel (PEC), which is a generalization of the broadcast binary erasure channel from the binary symbol to that of arbitrary finite fields GF(q) with sufficiently large q. We consider the setting in which the source node has instant feedback of the channel outputs of the K receivers after each transmission. Such a setting directly models network coded packet transmission in the downlink direction with integrated feedback mechanisms (such as Automatic Repeat reQuest (ARQ)). The main results of this paper are: (i) The capacity region for general 1-to-3 broadcast PECs, and (ii) The capacity region for two classes of 1-to-K broadcast PECs: the symmetric PECs, and the spatially independent PECs with one-sided fairness constraints. This paper also develops (iii) A pair of outer and inner bounds of the capacity region for arbitrary 1-to-K broadcast PECs, which can be evaluated by any linear programming solver. For most practical scenarios, the outer and inner bounds meet and thus jointly characterize the capacity.
Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to ${t}^{1/4}$ for surface diffusion and ${t}^{1/3}$ for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.
Gender and collaboration patterns in a temporal scientific authorship network

One can point to a variety of historical milestones for gender equality in STEM (science, technology, engineering, and mathematics), however, practical effects are incremental and ongoing. It is important to quantify gender differences in subdomains of scientific work in order to detect potential biases and monitor progress. In this work, we study the relevance of gender in scientific collaboration patterns in the Institute for Operations Research and the Management Sciences (INFORMS), a professional society with sixteen peer-reviewed journals. Using their publication data from 1952 to 2016, we constructed a large temporal bipartite network between authors and publications, and augmented the author nodes with gender labels. We characterized differences in several basic statistics of this network over time, highlighting how they have changed with respect to relevant historical events. We find a steady increase in participation by women (e.g., fraction of authorships by women and of new women authors) starting around 1980. However, women still comprise less than 25% of the INFORMS society and an even smaller fraction of authors with many publications. Moreover, we describe a methodology for quantifying the structural role of an authorship with respect to the overall connectivity of the network, using it to measure subtle differences between authorships by women and by men. Specifically, as measures of structural importance of an authorship, we use effective resistance and contraction importance, two measures related to diffusion throughout a network. As a null model, we propose a degree-preserving temporal and geometric network model with emergent communities. Our results suggest the presence of systematic differences between the collaboration patterns of men and women that cannot be explained by only local statistics.
Sequential multiple testing with generalized error control: An asymptotic optimality theory

The sequential multiple testing problem is considered under two generalized error metrics. Under the first one, the probability of at least $k$ mistakes, of any kind, is controlled. Under the second, the probabilities of at least $k_{1}$ false positives and at least $k_{2}$ false negatives are simultaneously controlled. For each formulation, the optimal expected sample size is characterized, to a first-order asymptotic approximation as the error probabilities go to 0, and a novel multiple testing procedure is proposed and shown to be asymptotically efficient under every signal configuration. These results are established when the data streams for the various hypotheses are independent and each local log-likelihood ratio statistic satisfies a certain strong law of large numbers. In the special case of i.i.d. observations in each stream, the gains of the proposed sequential procedures over fixed-sample size schemes are quantified.
Effects of spatial randomness on locating a point source with distributed sensors

Most studies that consider the problem of estimating the location of a point source in wireless sensor networks assume that the source location is estimated by a set of spatially distributed sensors, whose locations are fixed. Motivated by the fact that the observation quality and performance of the localization algorithm depend on the location of the sensors, which could be randomly distributed, this paper investigates the performance of a recently proposed energy-based source-localization algorithm under the assumption that the sensors are positioned according to a uniform clustering process. Practical considerations such as the existence and size of the exclusion zones around each sensor and the source will be studied. By introducing a novel performance measure called the estimation outage, it will be shown how parameters related to the network geometry such as the distance between the source and the closest sensor to it as well as the number of sensors within a region surrounding the source affect the localization performance.
Shells in the Magellanic System

The Magellanic System harbors >800 expanding shells of neutral hydrogen, providing a unique opportunity for statistical investigations. Most of these shells are surprisingly young, 2–10 Myr old, and correlate poorly with young stellar populations. I summarize what we have learned about shell properties and particularly focus on the puzzling correlation between the shell radius and expansion velocity. In the framework of the standard, adiabatic model for shell evolution this tight correlation suggests a coherent burst of star formation across the whole Magellanic System. However, more than one mechanism for shell formation may be taking place.
Cramped subgroups and generalized Harish-Chandra modules

Let G be a reductive complex Lie group with Lie algebra g. We call a subgroup H C G cramped if there is an integer b(G, H) such that each finite-dimensional representation of G has a non-trivial invariant subspace of dimension less than b(G, H). We show that a subgroup is cramped if and only if the moment map T*(K/L) - t* is surjective, where K and L are compact forms of G and H. We will use this in conjunction with sufficient conditions for crampedness given by Willenbring and Zuckerman (2004) to prove a geometric proposition on the intersections between adjoint orbits and Killing orthogonals to subgroups. We will also discuss applications of the techniques of symplectic geometry to the generalized Harish-Chandra modules introduced by Penkov and Zuckerman (2004), of which our results on crampedness are special cases.
Eventual regularization for the slightly supercritical quasi-geostrophic equation

Abstract We prove that weak solutions of the slightly supercritical quasi-geostrophic equation become smooth for large time. The proof uses ideas from a recent article of Caffarelli and Vasseur and is based on an argument in the style of De Giorgi.
Estimates for the number of visible galaxy-spanning civilizations and the cosmological expansion of life

If advanced civilizations appear in the universe with an ability and desire to expand, the entire universe can become saturated with life on a short timescale, even if such expanders appear rarely. Our presence in an apparently untouched Milky Way thus constrains the appearance rate of galaxy-spanning Kardashev type III (K3) civilizations, if it is assumed that some fraction of K3 civilizations will continue their expansion at intergalactic distances. We use this constraint to estimate the appearance rate of K3 civilizations for 81 cosmological scenarios by specifying the extent to which humanity is a statistical outlier. We find that in nearly all plausible scenarios, the distance to the nearest visible K3 is cosmological. In searches for K3 galaxies where the observable range is limited, we also find that the most likely detections tend to be expanding civilizations who have entered the observable range from farther away. An observation of K3 clusters is thus more likely than isolated K3 galaxies.