cssi/SciSciGPT-SciSciNet huggingface dataset:
An Optimal Lower Bound for Monotonicity Testing over Hypergrids.

For positive integers $n, d$, consider the hypergrid $[n]^d$ with the coordinate-wise product partial ordering denoted by $\prec$. A function $f: [n]^d \mapsto \mathbb{N}$ is monotone if $\forall x \prec y$, $f(x) \leq f(y)$. A function $f$ is $\eps$-far from monotone if at least an $\eps$-fraction of values must be changed to make $f$ monotone. Given a parameter $\eps$, a \emph{monotonicity tester} must distinguish with high probability a monotone function from one that is $\eps$-far. We prove that any (adaptive, two-sided) monotonicity tester for functions $f:[n]^d \mapsto \mathbb{N}$ must make $\Omega(\eps^{-1}d\log n - \eps^{-1}\log \eps^{-1})$ queries. Recent upper bounds show the existence of $O(\eps^{-1}d \log n)$ query monotonicity testers for hypergrids. This closes the question of monotonicity testing for hypergrids over arbitrary ranges. The previous best lower bound for general hypergrids was a non-adaptive bound of $\Omega(d \log n)$.
Self-Supervised Feature Learning via Exploiting Multi-Modal Data for Retinal Disease Diagnosis

The automatic diagnosis of various retinal diseases from fundus images is important to support clinical decision-making. However, developing such automatic solutions is challenging due to the requirement of a large amount of human-annotated data. Recently, unsupervised/self-supervised feature learning techniques receive a lot of attention, as they do not need massive annotations. Most of the current self-supervised methods are analyzed with single imaging modality and there is no method currently utilize multi-modal images for better results. Considering that the diagnostics of various vitreoretinal diseases can greatly benefit from another imaging modality, e.g ., FFA, this paper presents a novel self-supervised feature learning method by effectively exploiting multi-modal data for retinal disease diagnosis. To achieve this, we first synthesize the corresponding FFA modality and then formulate a patient feature-based softmax embedding objective. Our objective learns both modality-invariant features and patient-similarity features. Through this mechanism, the neural network captures the semantically shared information across different modalities and the apparent visual similarity between patients. We evaluate our method on two public benchmark datasets for retinal disease diagnosis. The experimental results demonstrate that our method clearly outperforms other self-supervised feature learning methods and is comparable to the supervised baseline. Our code is available at GitHub.
Backward Reachability Using Integral Quadratic Constraints for Uncertain Nonlinear Systems

A method is proposed to compute robust inner-approximations to the backward reachable set of nonlinear systems and to generate a robust control law that drives trajectories starting in these inner-approximations to a target set. The method merges dissipation inequalities and integral quadratic constraints (IQCs) with both hard and soft IQC factorizations, allowing for a variety of perturbations including parametric uncertainty, unmodeled dynamics, and uncertain time delays. Computational algorithms are developed using the generalized S-procedure and sum-of-squares techniques, and illustrated on a 6-state quadrotor with actuator uncertainties.
Near-Optimal Methods for Minimizing Star-Convex Functions and Beyond

In this paper, we provide near-optimal accelerated first-order methods for minimizing a broad class of smooth nonconvex functions that are strictly unimodal on all lines through a minimizer. This function class, which we call the class of smooth quasar-convex functions, is parameterized by a constant $\gamma \in (0,1]$, where $\gamma = 1$ encompasses the classes of smooth convex and star-convex functions, and smaller values of $\gamma$ indicate that the function can be "more nonconvex." We develop a variant of accelerated gradient descent that computes an $\epsilon$-approximate minimizer of a smooth $\gamma$-quasar-convex function with at most $O(\gamma^{-1} \epsilon^{-1/2} \log(\gamma^{-1} \epsilon^{-1}))$ total function and gradient evaluations. We also derive a lower bound of $\Omega(\gamma^{-1} \epsilon^{-1/2})$ on the number of gradient evaluations required by any deterministic first-order method in the worst case, showing that, up to a logarithmic factor, no deterministic first-order algorithm can improve upon ours.
Facial Expression and Peripheral Physiology Fusion to Decode Individualized Affective Experience

In this paper, we present a multimodal approach to simultaneously analyze facial movements and several peripheral physiological signals to decode individualized affective experiences under positive and negative emotional contexts, while considering their personalized resting dynamics. We propose a person-specific recurrence network to quantify the dynamics present in the person's facial movements and physiological data. Facial movement is represented using a robust head vs. 3D face landmark localization and tracking approach, and physiological data are processed by extracting known attributes related to the underlying affective experience. The dynamical coupling between different input modalities is then assessed through the extraction of several complex recurrent network metrics. Inference models are then trained using these metrics as features to predict individual's affective experience in a given context, after their resting dynamics are excluded from their response. We validated our approach using a multimodal dataset consists of (i) facial videos and (ii) several peripheral physiological signals, synchronously recorded from 12 participants while watching 4 emotion-eliciting video-based stimuli. The affective experience prediction results signified that our multimodal fusion method improves the prediction accuracy up to 19% when compared to the prediction using only one or a subset of the input modalities. Furthermore, we gained prediction improvement for affective experience by considering the effect of individualized resting dynamics.
Separating Effect From Significance in Markov Chain Tests

We give qualitative and quantitative improvements to theorems which enable significance testing in Markov chains, with a particular eye toward the goal of enabling strong, interpretable, and statis...
Hepatic Pseudoaneurysm Incidence After Liver Trauma.

Abstract Background Posttraumatic hepatic artery pseudoaneurysm is a potentially devastating complication after complex liver injury. Increasing computed tomography (CT) use may lead to more frequent identification of posttraumatic hepatic complications. This study was designed to determine the rate of hepatic pseudoaneurysm after traumatic liver injury. Methods We conducted a retrospective review of patients at an urban level 1 trauma center over 5 y (2012-2016). Injury characteristics, patient management, and complications were extracted from trauma registry data and chart review. Results Six hundred thirty-four hepatic injuries (11 no grade/no CT, 159 grade I, 154 grade II, 165 grade III, 93 grade IV, and 52 grade V) were identified from our trauma registry. No patient with a grade I or II injury had a subsequent bleeding complication. Eighteen patients had a documented hepatic pseudoaneurysm: grade III n = 3 (1.8%), grade IV n = 6 (6.5%), grade V n = 9 (17.3%). The median time to pseudoaneurysm identification was 6.5 d. Seven pseudoaneurysms were found on asymptomatic surveillance CT-angiography on average 5 d after injury. Eleven patients were symptomatic at the time of CT-angiography performed at a median of 9 d after admission. Of the 11 symptomatic patients, four were in hemorrhagic shock, and two died from hepatic-related complications. Conclusions The incidence of hepatic artery pseudoaneurysm increases with higher grade liver injury. Aggressive surveillance for hepatic pseudoaneurysm with interval CT-angiography 5-7 d postinjury may be warranted, especially for grade IV and V injuries.
Does Telemedicine Reduce Emergency Room Congestion? Evidence from New York State

Overcrowding in emergency rooms (ERs) is a common yet nagging problem. It not only is costly for hospitals but also compromises care quality and patient experience. Our paper provides solid evidenc...
A Dual-Port 8-T CAM-Based Network Intrusion Detection Engine for IoT

This letter presents an energy- and memory-efficient pattern-matching engine for a network intrusion detection system (NIDS) in the Internet of Things. Tightly coupled architecture and circuit co-designs are proposed to fully exploit the statistical behaviors of NIDS pattern matching. The proposed engine performs pattern matching in three phases, where the phase-1 prefix matching employs reconfigurable pipelined automata processing to minimize memory footprint without loss of throughput and efficiency. The processing elements utilize 8-T content-addressable memory (CAM) cells for dual-port search by leveraging proposed fixed-1s encoding. A 65-nm prototype demonstrates best-in-class 1.54-fJ energy per search per pattern byte and 0.9-byte memory usage per pattern byte.
Discovery of Drug-like Ligands for the Mac1 Domain of SARS-CoV-2 Nsp3

Small molecules that bind the SARS-CoV-2 nonstructural protein 3 Mac1 domain in place of ADP-ribose could be useful as molecular probes or scaffolds for COVID-19 antiviral drug discovery because Mac1 has been linked to the ability of coronaviruses to evade cellular detection. A high-throughput assay based on differential scanning fluorimetry (DSF) was therefore optimized and used to identify possible Mac1 ligands in small libraries of drugs and drug-like compounds. Numerous promising compounds included nucleotides, steroids, β-lactams, and benzimidazoles. The main drawback to this approach was that a high percentage of compounds in some libraries were found to influence the observed Mac1 melting temperature. To prioritize DSF screening hits, the shapes of the observed melting curves and initial assay fluorescence were examined, and the results were compared with virtual screens performed using AutoDock Vina. The molecular basis for alternate ligand binding was also examined by determining a structure of one of the hits, cyclic adenosine monophosphate, with atomic resolution.
Defending Hardware-Based Malware Detectors Against Adversarial Attacks

In the era of Internet of Things (IoT), Malware has been proliferating exponentially over the past decade. Traditional anti-virus software are ineffective against modern complex Malware. In order to address this challenge, researchers have proposed hardware-assisted Malware detection (HMD) using hardware performance counters (HPCs). The HPCs are used to train a set of machine learning (ML) classifiers, which in turn, are used to distinguish benign programs from Malware. Recently, adversarial attacks have been designed by introducing perturbations in the HPC traces using an adversarial sample predictor to misclassify a program for specific HPCs. These attacks are designed with the basic assumption that the attacker is aware of the HPCs being used to detect Malware. Since modern processors consist of hundreds of HPCs, restricting to only a few of them for Malware detection aids the attacker. In this article, we propose a moving target defense (MTD) for this adversarial attack by designing multiple ML classifiers trained on different sets of HPCs. The MTD randomly selects a classifier; thus, confusing the attacker about the HPCs or the number of classifiers applied. We have developed an analytical model which proves that the probability of an attacker to guess the perfect HPC-classifier combination for MTD is extremely low (in the range of $10^{-1864}$ for a system with 20 HPCs). Our experimental results prove that the proposed defense is able to improve the classification accuracy of HPC traces that have been modified through an adversarial sample generator by up to 31.5%, for a near perfect (99.4%) restoration of the original accuracy.
Q* Approximation Schemes for Batch Reinforcement Learning: A Theoretical Comparison

We prove performance guarantees of two algorithms for approximating $Q^\star$ in batch reinforcement learning. Compared to classical iterative methods such as Fitted Q-Iteration---whose performance loss incurs quadratic dependence on horizon---these methods estimate (some forms of) the Bellman error and enjoy linear-in-horizon error propagation, a property established for the first time for algorithms that rely solely on batch data and output stationary policies. One of the algorithms uses a novel and explicit importance-weighting correction to overcome the infamous "double sampling" difficulty in Bellman error estimation, and does not use any squared losses. Our analyses reveal its distinct characteristics and potential advantages compared to classical algorithms.
Identifying causal effects in maximally oriented partially directed acyclic graphs.

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

The Wide-field Imager for Solar Probe (WISPR) onboard the Parker Solar Probe (PSP), observing in white light, has a fixed angular field of view, extending from 13.5∘ to 108∘ from the Sun and approximately 50∘ in the transverse direction. Because of the highly elliptical orbit of PSP, the physical extent of the imaged coronal region varies directly as the distance from the Sun, requiring new techniques for analysis of the motions of observed density features. Here, we present a technique for determining the 3D trajectory of CMEs and other coronal ejecta moving radially at a constant velocity by first tracking the motion in a sequence of images and then applying a curve-fitting procedure to determine the trajectory parameters (distance vs. time, velocity, longitude, and latitude). To validate the technique, we have determined the trajectory of two CMEs observed by WISPR that were also observed by another white-light imager, either the Solar and Heliophysics Observatory (SOHO) / Large Angle and Spectrometric COronagraph (LASCO)-C3 or the Solar Terrestrial Relations Observatory (STEREO)-A/HI1. The second viewpoint was used to verify the trajectory results from this new technique and help determine its uncertainty.
On the Nonuniform Berry–Esseen Bound

Abstract Due to the effort of a number of authors, the value c u of the absolute constant factor in the uniform Berry–Esseen (BE) bound for sums of independent random variables has been gradually reduced to 0.4748 in the iid case and 0.5600 in the general case; both these values were recently obtained by Shevtsova. On the other hand, Esseen had shown that c u cannot be less than 0.4097. Thus, the gap factor between the best known upper and lower bounds on (the least possible value of) c u is now rather close to 1. The situation is quite different for the absolute constant factor c nu in the corresponding nonuniform BE bound. Namely, the best correctly established upper bound on c nu in the iid case is over 17 times the corresponding best known lower bound, and this gap factor is greater than 21 in the general case. In the present paper, improvements to the prevailing method (going back to S. Nagaev) of obtaining nonuniform BE bounds are suggested. Moreover, two new methods are presented, of a rather purely Fourier kind, based on two families of smoothing inequalities, which work better in the tail zones. As illustrations, two corresponding quick proofs of Nagaev’s nonuniform BE bound are given. Some further refinements in the application of the methods are shown as well.
What are Ophthalmology Patients Asking Online? An Analysis of the Eye Triage Subreddit

Importance Ophthalmology patients are seeking medical advice on social media websites like Reddit, where users are able to post comments and discuss issues pertaining to different topics that are organized in 'subreddits'. Understanding which issues are most pertinent will guide ophthalmic providers in delivering more effective patient education. Methods This cross-sectional study assessed a systematic sample of the first 22 posts and their top 3 comments from each month since January 27th, 2019, the subreddit's creation. Information was gathered from reddit.com/r/eyetriage in October 2019 and analyzed in November 2019. Main Outcomes The posts were characterized by date and time, inclusion of an image, type, content, emotional tone, and number of upvotes and comments. The comments were categorized based on content, emotional tone, time of comment, and user background. Post and comment content codes were categorized in an iterative manner with differences resolved by author consensus. Categorical statistics were compiled. Results Two hundred posts and 456 comments were analyzed since the creation of r/eyetriage, a forum created exclusively for patients to seek advice from health-care professionals. Twenty-six (13%) of the total posts included an image. On average, comments received 1.76 ± 2.17 upvotes along with 4.50 ± 4.47 replies. The most common content codes among the posts were 42 (21.0%) seeking diagnoses, 23 (11.5%) surgical complications, and 13 (6.50%) alternative medication options. Eighty-two (41%) posts conveyed a clear emotional tone, most notably 11 (13.4%) with anxiety and 10 (12.2%) with worry. The top comments came from 165 (36.2%) self-identified patients, 151 (33.1%) optometrists, and 49 (10.8%) ophthalmologists. The top comment codes for replies included 158 (34.7%) with treatment advice, 70 (15.4%) with advice deferred to follow-up appointment with other health-care specialists, and 60 (13.2%) with sharing information. Conclusions Patients are asking ophthalmology-related questions on the Eye Triage subreddit, and they are more likely to receive information from other patients or optometrists than from self-identified ophthalmologists. When emotions were revealed, patients were often anxious and worried. Opportunities exist for ophthalmologists to take a more active role on this subreddit and help educate patients.
Random band matrices in the delocalized phase, III: averaging fluctuations

We consider a general class of symmetric or Hermitian random band matrices $$H=(h_{xy})_{x,y \in \llbracket 1,N\rrbracket ^d}$$ in any dimension $$d\ge 1$$ , where the entries are independent, centered random variables with variances $$s_{xy}=\mathbb {E}|h_{xy}|^2$$ . We assume that $$s_{xy}$$ vanishes if $$|x-y|$$ exceeds the band width W, and we are interested in the mesoscopic scale with $$1\ll W\ll N$$ . Define the generalized resolvent of H as $$G(H,Z):=(H - Z)^{-1}$$ , where Z is a deterministic diagonal matrix with entries $$Z_{xx}\in \mathbb {C}_+$$ for all x. Then we establish a precise high-probability bound on certain averages of polynomials of the resolvent entries. As an application of this fluctuation averaging result, we give a self-contained proof for the delocalization of random band matrices in dimensions $$d\ge 2$$ . More precisely, for any fixed $$d\ge 2$$ , we prove that the bulk eigenvectors of H are delocalized in certain averaged sense if $$N\le W^{1+\frac{d}{2}}$$ . This improves the corresponding results in He and Marcozzi (Diffusion profile for random band matrices: a short proof, 2018. arXiv:1804.09446 ) that imposed the assumption $$N\ll W^{1+\frac{d}{d+1}}$$ , and the results in Erdős and Knowles (Ann Henri Poincare12(7):1227–1319, 2011; Commun Math Phys 303(2): 509–554, 2011) that imposed the assumption $$N\ll W^{1+\frac{d}{6}}$$ . For 1D random band matrices, our fluctuation averaging result was used in Bourgade et al. (J Stat Phys 174:1189–1221, 2019; Random band matrices in the delocalized phase, I: quantum unique ergodicity and universality, 2018. arXiv:1807.01559 ) to prove the delocalization conjecture and bulk universality for random band matrices with $$N\ll W^{4/3}$$ .
ScratchR: sharing user-generated programmable media

In this paper, I describe a platform for sharing programmable media on the web called ScratchR. As the backbone of an on-line community of creative learners, ScratchR will give members access to an audience and inspirational ideas from each other. ScratchR seeks to support different states of participation: from passive consumption to active creation. This platform is being evaluated with a group of middle-school students and a larger community of beta testers.
The Dynamical Inverse Problem for Axisymmetric Stellar Systems

The standard method of modelling axisymmetric stellar systems begins from the assumption that mass follows light. The gravitational potential is then derived from the luminosity distribution, and the unique two-integral distribution function f(E,Lz) that generates the stellar density in this potential is found. It is shown that the gravitational potential can instead be generated directly from the velocity data in a two-integral galaxy, thus allowing one to drop the assumption that mass follows light. The two-dimensional rotational velocity field can also be recovered in a model-independent way. Regularized algorithms for carrying out the inversions are presented and tested by application to pseudo-data from a family of oblate models.
THE GBT 67–93.6 GHz SPECTRAL LINE SURVEY OF ORION-KL

We present a 67–93.6 GHz spectral line survey of Orion-KL with the new 4 mm Receiver on the Green Bank Telescope (GBT). The survey reaches unprecedented depths and covers the low-frequency end of the 3 mm atmospheric window which has been relatively unexplored previously. The entire spectral-line survey is published electronically for general use by the astronomical community. The calibration and performance of the 4 mm Receiver on the GBT is also summarized.
Invariant barriers to reactive front propagation in fluid flows

We present theory and experiments on the dynamics of reaction fronts in two-dimensional, vortex-dominated flows, for both time-independent and periodically driven cases. We find that the front propagation process is controlled by one-sided barriers that are either fixed in the laboratory frame (time-independent flows) or oscillate periodically (periodically driven flows). We call these barriers burning invariant manifolds (BIMs), since their role in front propagation is analogous to that of invariant manifolds in the transport and mixing of passive impurities under advection. Theoretically, the BIMs emerge from a dynamical systems approach when the advection-reaction-diffusion dynamics is recast as an ODE for front element dynamics. Experimentally, we measure the location of BIMs for several laboratory flows and confirm their role as barriers to front propagation.
The 2012 July 23 Backside Eruption: An Extreme Energetic Particle Event?

The backside coronal mass ejection (CME) of 2012 July 23 had a short Sun-to-Earth shock transit time (18.5 hr).The associated solar energetic particle (SEP) event had a greater than 10 MeV proton flux peaking at approximately 5000 pfu, and the energetic storm particle event was an order of magnitude larger, making it the most intense event in the space era at these energies. By a detailed analysis of the CME, shock, and SEP characteristics, we find that the July 23 event is consistent with a high-energy SEP event (accelerating particles to giga-electron volt energies). The times of maximum and fluence spectra in the range 10100 MeV were very hard, similar to those of ground-level enhancement (GLE) events. We found a hierarchical relationship between the CME initial speeds and the fluence spectral indices: CMEs with low initial speeds had SEP events with the softest spectra, while those with the highest initial speeds had SEP events with the hardest spectra. CMEs attaining intermediate speeds result in moderately hard spectra. The July 23 event was in the group of hard-spectrum events. During the July 23 event, the shock speed greater than (2000 km s(exp -1), the initial acceleration (approximately 1.70 km s(exp -2), and the shock-formation height (approximately 1.5 solar radii)were all typical of GLE events. The associated type II burst had emission components from meter to kilometer wavelengths, suggesting a strong shock. These observations confirm that the 2012 July 23 event is likely to be an extreme event in terms of the energetic particles it accelerated.
The life of a vortex knot

The idea that the knottedness (hydrodynamic Helicity) of a fluid flow is conserved has a long history in fluid mechanics. The quintessential example of a knotted flow is a knotted vortex filament, however, owing to experimental difficulties, it has not been possible until recently to directly generate knotted vortices in real fluids. Using 3D printed hydrofoils and high-speed laser scanning tomography, we generate vortex knots and links and measure their subsequent evolution. In both cases, we find that the vortices deform and stretch until a series of vortex reconnections occurs, eventually resulting several disjoint vortex rings. This article accompanies a fluid dynamics video entered into the Gallery of Fluid Motion at the 66th Annual Meeting of the APS Division of Fluid Dynamics.
EFFECTS OF MERIDIONAL FLOW VARIATIONS ON SOLAR CYCLES 23 AND 24

The faster meridional flow that preceded the solar cycle 23/24 minimum is thought to have led to weaker polar field strengths, producing the extended solar minimum and the unusually weak cycle 24. To determine the impact of meridional flow variations on the sunspot cycle, we have simulated the Sun's surface magnetic field evolution with our newly developed surface flux transport model. We investigate three different cases: a constant average meridional flow, the observed time-varying meridional flow, and a time-varying meridional flow in which the observed variations from the average have been doubled. Comparison of these simulations shows that the variations in the meridional flow over cycle 23 have a significant impact (~20%) on the polar fields. However, the variations produced polar fields that were stronger than they would have been otherwise. We propose that the primary cause of the extended cycle 23/24 minimum and weak cycle 24 was the weakness of cycle 23 itself—with fewer sunspots, there was insufficient flux to build a big cycle. We also find that any polar counter-cells in the meridional flow (equatorward flow at high latitudes) produce flux concentrations at mid-to-high latitudes that are not consistent with observations.
Additive quantile regression for clustered data with an application to children's physical activity

Additive models are flexible regression tools that handle linear as well as non-linear terms. The latter are typically modelled via smoothing splines. Additive mixed models extend additive models to include random terms when the data are sampled according to cluster designs (e.g. longitudinal).These models find applications in the study of phenomena like growth, certain disease mechanisms and energy expenditure in humans, when repeated measurements are available. We propose a novel additive mixed model for quantile regression. Our methods are motivated by an application to physical activity based on a data set with more than half a million accelerometer measurements in children of the UK Millennium Cohort Study. In a simulation study, we assess the proposed methods against existing alternatives.
Hamiltonian simulation with nearly optimal dependence on spectral norm

We present a quantum algorithm for approximating the real time evolution e−iHt of an arbitrary d-sparse Hamiltonian to error є, given black-box access to the positions and b-bit values of its non-zero matrix entries. The query complexity of our algorithm is O((t√d||H||1 → 2)1+o(1)/єo(1)) with respect to the largest Euclidean row norm ||H||1 → 2, which is shown to be optimal up to subpolynomial factors through a matching lower bound, and it uses a factor O(b) more gates. This provides a polynomial speedup in sparsity for the common case where the spectral norm is known, and generalizes previous approaches which achieve optimal scaling, but with respect to more restrictive parameters. By exploiting knowledge of the spectral norm, our algorithm solves the black-box unitary implementation problem – O(d1/2+o(1)) queries suffice to approximate any d-sparse unitary in the black-box setting, which matches the quantum search lower bound of Ω(√d) queries and improves upon prior art [Berry and Childs, QIP 2010] of O(d2/3) queries. Combined with known techniques, we also solve systems of sparse linear equations with condition number κ using O((κ √d)1+o(1)/єo(1)) queries, which is a quadratic improvement in sparsity.
Blazar Sheath Illumination of the Outer Molecular Torus: A Resolution of the Seed Photon Problem for the Far-GeV Blazar Flares

Recent multi-wavelength work led by the Boston University blazar group (e.g., Marscher et al. 2010) strongly suggests that a fraction of the blazar flares seen by the Fermi Large Area Telescope (LAT) take place a few to several pc away from the central engine. However, at such distances from the central engine, there is no adequate external photon field to provide the seed photons required for producing the observed GeV emission under leptonic inverse Compton (IC) models. A possible solution is a spine-sheath geometry for the emitting region (MacDonald et al. 2015, but see Nalewajko et al. 2014). Here we use the current view of the molecular torus (e.g., Elitzur, 2012; Netzer 2015) in which the torus extends a few pc beyond the dust sublimation radius with dust clouds distributed with a declining density for decreasing polar angle. We show that for a spine-sheath blazar jet embedded in the torus, the wide beaming pattern of the synchrotron radiation of the relatively slow sheath will heat molecular clouds whose subsequent IR radiation will be seen highly boosted in the spine comoving frame, and that under reasonable conditions this photon field can dominate over the sheath photons directly entering the spine. If the sheath is sufficiently luminous it will sublimate the dust, and if the sheath synchrotron radiation extends to optical-UV energies (as may happen during flares), this will illuminate the sublimated dust clouds to produce emission lines that will vary in unison with the optical-UV continuum, as has been very recently reported for blazar CTA 102 (Jorstad et al. 2017).
SEEDING THE FORMATION OF COLD GASEOUS CLOUDS IN MILKY WAY-SIZE HALOS

We use one of the highest resolution cosmological smoothed particle hydrodynamic simulations to date to demonstrate that cold gaseous clouds form around Milky Way-size galaxies. We further explore mechanisms responsible for their formation and show that a large fraction of clouds originate as a consequence of late-time filamentary "cold mode" accretion. Here, filaments that are still colder and denser than the surrounding halo gas are not able to connect directly to galaxies, as they do at high redshift, but are instead susceptible to the combined action of cooling and Rayleigh-Taylor instabilities at intermediate radii within the halo leading to the production of cold, dense pressure-confined clouds, without an associated dark matter component. This process is aided through the compression of the incoming filaments by the hot halo gas and expanding shocks during the halo buildup. Our mechanism directly seeds clouds from gas with substantial local overdensity, unlike in previous models, and provides a channel for the origin of cloud complexes. These clouds can later "rain" onto galaxies, delivering fuel for star formation. Owing to the relatively large cross-section of filaments and the net angular momentum carried by the gas, the clouds will be distributed in a modestly flattened region around a galaxy.
Wave propagation in periodic networks of thin fibers

We will discuss a one-dimensional approximation for the problem of wave propagation in networks of thin fibers. The main objective here is to describe the boundary (gluing) conditions at branching points of the limiting one-dimensional graph. The results will be applied to Mach–Zehnder interferometers on chips and to periodic chains of the interferometers. The latter allows us to find parameters which guarantee the transparency and slowing down of wave packets.
Applications of differential algebra to algebraic independence of arithmetic functions

We generalize and unify the proofs of several results on algebraic in- dependence of arithmetic functions and Dirichlet series by a theorem of Ax on differential Schanuel conjecture. Along the way, we find counter-examples to some results in the literature.
Phase transition in the rich-get-richer mechanism due to finite-size effects

The rich-get-richer mechanism (agents increase their 'wealth' randomly at a rate proportional to their holdings) is often invoked to explain the Pareto power-law distribution observed in many physical situations, such as the degree distribution of growing scale-free nets. We use two different analytical approaches, as well as numerical simulations, to study the case where the number of agents is fixed and finite (but large), and the rich-get-richer mechanism is invoked a fraction r of the time (the remainder of the time wealth is disbursed by a homogeneous process). At short times, we recover the Pareto law observed for an unbounded number of agents. In later times, the (moving) distribution can be scaled to reveal a phase transition with a Gaussian asymptotic form for , and a Pareto-like tail (on the positive side) and a novel stretched exponential decay (on the negative side) for .
Supersymmetry without Grassmann variables

Supersymmetry transformations may be represented by unitary operators in a formulation of supersymmetry without numbers that anti-commute. The physical relevance of this formulation hinges on whether or not one may add states of even and odd fermion number, a question which soon may be settled by experiment.
Graded cofinite rings of differential operators

In this paper we study subalgebras A of the algebra D(X) of differential operators on a smooth variety X which are big in the following sense: using the order of a differential operator, the ring D(X) is equipped with a filtration. Its associated graded algebra D(X) is commutative and can be regarded as the set of regular functions on the cotangent bundle ofX . The subalgebra A inherits a filtration from D(X) and its associated graded algebraA is a subalgebra of D(X). We call A graded cofinite in D(X) if D(X) is a finitely generated A-module. Our guiding example of a graded cofinite subalgebra is the algebra of invariants D(X) where W is a finite group acting on X . Other examples can be constructed as follows. Let φ : X → Y be a finite dominant morphism onto a normal variety Y . Then we put
The 2-Channel Kondo Model: II. CFT Calculation of Non-equilibrium Conductance through a Nanoconstriction Containing 2-Channel Kondo Impurities

Abstract Recent experiments by Ralph and Buhrman on zero-bias anomalies in quenched Cu nanoconstrictions (reviewed in a preceding paper) are in accord with the assumption that the interaction between electrons and nearly degenerate two-level systems in the constriction can be described, for sufficiently small voltages and temperatures (V, T
Effects of the Dissolution of Low-Concentration Globular Clusters on the Evolution of Globular Cluster Systems

We investigate the role of dissolution of low-concentration clusters due to mass loss through stellar evolution on the evolution of the properties of globular cluster systems (GCSs) in elliptical galaxies. Our simulations show that, for an initial mass-concentration relationship based on that inferred from Galactic globular clusters, dissolution of low-concentration clusters leads to the disruption of a large number of clusters. A power-law initial globular cluster system mass function (GCMF) similar to that observed in young cluster systems in merging galaxies is transformed by this dissolution into a bell-shaped GCMF with a mean mass similar to that of old GCSs for all the galaxies investigated. Two-body relaxation and dynamical friction, which are also included in our simulations, subsequently lead to an additional significant evolution and disruption of the population of clusters. As shown previously, when these processes act on a bell-shaped GCMF with a mean mass similar to that of old GCSs, they do not significantly alter the value of the mean mass. The final GCMFs are bell-shaped with similar peaks at different radii within galaxies and between different galaxies, in agreement with current observations.
On Stochastic Sea of the Standard Map

Consider a generic one-parameter unfolding of a homoclinic tangency of an area preserving surface diffeomorphism. We show that for many parameters (residual subset in an open set approaching the critical value) the corresponding diffeomorphism has a transitive invariant set Ω of full Hausdorff dimension. The set Ω is a topological limit of hyperbolic sets and is accumulated by elliptic islands.
Comprehensive Solution to the Cosmological Constant, Zero-Point Energy, and Quantum Gravity Problems

We present a solution to the cosmological constant, the zero-point energy, and the quantum gravity problems within a single comprehensive framework. We show that in quantum theories of gravity in which the zero-point energy density of the gravitational field is well-defined, the cosmological constant and zero-point energy problems solve each other by mutual cancellation between the cosmological constant and the matter and gravitational field zero-point energy densities. Because of this cancellation, regulation of the matter field zero-point energy density is not needed, and thus does not cause any trace anomaly to arise. We exhibit our results in two theories of gravity that are well-defined quantum-mechanically. Both of these theories are locally conformal invariant, quantum Einstein gravity in two dimensions and Weyl-tensor-based quantum conformal gravity in four dimensions (a fourth-order derivative quantum theory of the type that Bender and Mannheim have recently shown to be ghost-free and unitary). Central to our approach is the requirement that any and all departures of the geometry from Minkowski are to be brought about by quantum mechanics alone. Consequently, there have to be no fundamental classical fields, and all mass scales have to be generated by dynamical condensates. In such a situation the trace of the matter field energy-momentum tensor is zero, a constraint that obliges its cosmological constant and zero-point contributions to cancel each other identically, no matter how large they might be. In our approach quantization of the gravitational field is caused by its coupling to quantized matter fields, with the gravitational field not needing any independent quantization of its own. With there being no a priori classical curvature, one does not have to make it compatible with quantization.
The Exceptional Zero Conjecture for Symmetric Powers of CM Modular Forms: the Ordinary Case

We prove the exceptional zero conjecture for the symmetric powers of CM cuspidal eigenforms at ordinary primes. In other words, we determine the trivial zeroes of the associated p-adic L-functions, compute the L-invariants, and show that they agree with Greenberg's L-invariants.
New Description of Evolution of Magnetic Phases in Artificial Honeycomb Lattice.

Artificial magnetic honeycomb lattice provides a two-dimensional archetypal system to explore novel phenomena of geometrically frustrated magnets. According to theoretical reports, an artificial magnetic honeycomb lattice is expected to exhibit several phase transitions to unique magnetic states as a function of reducing temperature. Experimental investigations of permalloy artificial honeycomb lattice of connected ultra-small elements, $${\boldsymbol{\simeq }}$$ ≃ 12 nm, reveal a more complicated behavior. First, upon cooling the sample to intermediate temperature, $${\bf{T}}{\boldsymbol{\simeq }}$$ T ≃ 175 K, the system manifests a non-unique state where the long range order co-exists with short-range magnetic charge order and weak spin ice state. Second, at much lower temperature, $${\bf{T}}{\boldsymbol{\simeq }}$$ T ≃ 6 K, the long-range spin solid state exhibits a re-entrant behavior. Both observations are in direct contrast to the present understanding of this system. New theoretical approaches are needed to develop a comprehensive formulation of this two dimensional magnet.
Shorter Stabilizer Circuits via Bruhat Decomposition and Quantum Circuit Transformations

In this paper, we improve the layered implementation of arbitrary stabilizer circuits introduced by Aaronson and Gottesman in Phys. Rev. A 70 (052328), 2004: to implement a general stabilizer circuit, we reduce their 11-stage computation -H-C-P-C-P-C-H-P-C-P-C- over the gate set consisting of Hadamard, controlled-NOT, and phase gates, into a 7-stage computation of the form -C-CZ-P-H-P-CZ-C-. We show arguments in support of using -CZ- stages over the -C- stages: not only the use of -CZ- stages allows a shorter layered expression, but -CZ- stages are simpler and appear to be easier to implement compared to the -C- stages. Based on this decomposition, we develop a two-qubit gate depth- $(14n{-}4)$ implementation of stabilizer circuits over the gate library $\{ \mathrm{H}, \mathrm{P}, \mathrm{CNOT}\}$ , executable in the Linear Nearest Neighbor (LNN) architecture, improving best previously known depth- $25n$ circuit, also executable in the LNN architecture. Our constructions rely on Bruhat decomposition of the symplectic group and on folding arbitrarily long sequences of the form (-P-C-) $^{m}$ into a three-stage computation -P-CZ-C-. Our results include the reduction of the 11-stage decomposition -H-C-P-C-P-C-H-P-C-P-C- into a 9-stage decomposition of the form -C-P-C-P-H-C-P-C-P-. This reduction is based on the Bruhat decomposition of the symplectic group. This result also implies a new normal form for stabilizer circuits. We show that a circuit in this normal form is optimal in the number of Hadamard gates used. We also show that the normal form has an asymptotically optimal number of parameters.
Future Directions in the Microwave Cavity Search for Dark Matter Axions

The axion is a light pseudoscalar particle which suppresses CP-violating effects in strong interactions and also happens to be an excellent dark matter candidate. Axions constituting the dark matter halo of our galaxy may be detected by their resonant conversion to photons in a microwave cavity permeated by a magnetic field. The current generation of the microwave cavity experiment has demonstrated sensitivity to plausible axion models, and upgrades in progress should achieve the sensitivity required for a definitive search, at least for low mass axions. However, a comprehensive strategy for scanning the entire mass range, from 1–1000 μeV, will require significant technological advances to maintain the needed sensitivity at higher frequencies. Such advances could include sub-quantum-limited amplifiers based on squeezed vacuum states, bolometers, and/or superconducting microwave cavities. The Axion Dark Matter eXperiment at High Frequencies (ADMX-HF) represents both a pathfinder for first data in the 20–100 μeV range (~5–25 GHz), and an innovation test-bed for these concepts.
Next-to-leading corrections to the Higgs boson transverse momentum spectrum in gluon fusion

We present a fully analytic calculation of the Higgs boson transverse momentum and rapidity distributions, for nonzero Higgs p⊥, at next-to-leading order in the infinite-top-mass approximation. We separate the cross section into a part that contains the dominant soft, virtual, collinear, and small-p⊥-enhanced contributions, and the remainder, which is organized by the contributions due to different parton helicities. We use this cross section to investigate analytically the small-p⊥ limit and compare with the expectation from the resummation of large logarithms of the type ln mH/p⊥. We also compute numerically the cross section at moderate p⊥ where a fixed-order calculation is reliable. We find a K-factor that varies from ≈ 1.6−1.8, and a reduction in the scale dependence, as compared to leading order. Our analysis suggests that the contribution of current parton distributions to the total uncertainty on this cross section at the LHC is probably less than that due to uncalculated higher orders.
The Stellar Initial Mass Function in Early-Type Galaxies from Absorption Line Spectroscopy. III. Radial Gradients

There is good evidence that the centers of massive early-type galaxies have a bottom-heavy stellar initial mass function (IMF) compared to the IMF of the Milky Way. Here we study the radial variation of the IMF within such galaxies, using a combination of high quality Keck spectroscopy and a new suite of stellar population synthesis models that cover a wide range in metallicity. As in the previous studies in this series, the models are fitted directly to the spectra and treat all elemental abundance ratios as free parameters. Using newly obtained spectroscopy for six galaxies, including deep data extending to ~1Re for the galaxies NGC1407, NGC1600, and NGC2695, we find that the IMF varies strongly with galactocentric radius. For all six galaxies the IMF is bottom-heavy in the central regions, with average mass-to-light ratio "mismatch" parameter a~2.5 at R=0. The IMF rapidly becomes more bottom-light with increasing radius, flattening off near the Milky Way value (a~1.1) at R>0.4Re. A consequence is that the luminosity-weighted average IMF depends on the measurement aperture: within R=Re we find =1.3-1.5, consistent with recent lensing and dynamical results from SLACS and ATLAS-3D. Our results are also consistent with several earlier studies that were based on analyses of radial gradients of line indices. The observed IMF gradients support galaxy formation models in which the central regions of massive galaxies had a different formation history than their outer parts. Finally, we make use of the high signal-to-noise central spectra of NGC1407 and NGC2695 to demonstrate how we can disentangle IMF effects and abundance effects.
CHEAP SPACE-BASED MICROLENS PARALLAXES FOR HIGH-MAGNIFICATION EVENTS

We show that for high-magnification (A{sub max} {approx}> 100) microlensing events, accurate microlens parallaxes can be obtained from three or fewer photometric measurements from a small telescope on a satellite in solar orbit at O(AU) from Earth. This is 1-2 orders of magnitude less observing resources than are required for standard space-based parallaxes. Such microlens parallax measurements would yield accurate mass and distance measurements to the lens for all cases in which finite-source effects were observed from the ground over peak. This would include virtually all high-magnification events with detected planets and a substantial fraction of those without. Hence, it would permit accurate estimates of the Galactic distribution of planets.