title: On spectral asymptotic of quasi-exactly solvable quartic

abstract: Motivated by the earlier results, we study theoretically and numerically the asymptotics and the monodromy of the quasi-exactly solvable part of the spectrum of the quasi-exactly solvable quartic introduced by C.~M.~Bender and S.~Boettcher. In particular, we formulate a conjecture on the coincidence of the asymptotic shape of the configuration of the branching points of the latter quartic with the asymptotic shape of zeros of the Yablonskii-Vorob'ev polynomials recently described and present its (conjectural) alternative description. Further we present a numerical study of the spectral monodromy for the os- cillator in question.


title: Gate-Tunable Resonant Tunneling in Double Bilayer Graphene Heterostructures

abstract: We demonstrate gate-tunable resonant tunneling and negative differential resistance in the interlayer current-voltage characteristics of rotationally aligned double bilayer graphene heterostructures separated by hexagonal boron-nitride (hBN) dielectric. An analysis of the heterostructure band alignment using individual layer densities, along with experimentally determined layer chemical potentials indicates that the resonance occurs when the energy bands of the two bilayer graphene are aligned. We discuss the tunneling resistance dependence on the interlayer hBN thickness, as well as the resonance width dependence on mobility and rotational alignment.


title: Atomic limits in the search for galactic dark matter

abstract: Direct searches for low mass dark matter particles via scattering off target nuclei require detection of recoiling atoms with energies of ~1 keV or less. The amount of electronic excitation produced by such atoms is quenched relative to a recoiling electron of the same energy. The Lindhard model of this quenching, as originally formulated, remains widely used after more than 50 years. The present work shows that for very small energies, a simplifying approximation of that model must be removed. Implications for the sensitivity of direct detection experiments are discussed.


title: k=0 magnetic structure and absence of ferroelectricity in SmFeO3

abstract: SmFeO3 has attracted considerable attention very recently due to the reported multiferroic properties above room-temperature. We have performed powder and single crystal neutron diffraction as well as complementary polarization dependent soft X-ray absorption spectroscopy measurements on floating-zone grown SmFeO3 single crystals in order to determine its magnetic structure. We found a k=0 G-type collinear antiferromagnetic structure that is not compatible with inverse Dzyaloshinskii-Moriya interaction driven ferroelectricity. While the structural data reveals a clear sign for magneto-elastic coupling at the N\'eel-temperature of ~675 K, the dielectric measurements remain silent as far as ferroelectricity is concerned.


title: Singular behavior of fluctuations in a relaxation process

abstract: Carrying out explicitly the computation in a paradigmatic model of non-interacting systems, the Gaussian Model, we show the existence of a singular point in the probability distribution $P(M)$ of an extensive variable $M$. Interpreting $P(M)$ as a thermodynamic potential of a dual system obtained from the original one by applying a constraint, we discuss how the non-analytical point of $P(M)$ is the counterpart of a phase-transition in the companion system. We show the generality of such mechanism by considering both the system in equilibrium or in the non-equilibrium state following a temperature quench.


title: Electromagnetically induced transparency of a single-photon in dipole-coupled one-dimensional atomic clouds

abstract: We investigate the propagation of a single photon under conditions of electromagnetically induced transparency in two parallel one-dimensional atomic clouds which are coupled via Rydberg dipole-dipole interaction. Initially the system is prepared with a single delocalized Rydberg excitation shared between the two ensembles and the photon enters both of them in an arbitrary path-superposition state. By properly aligning the transition dipoles of the atoms of each cloud we show that the photon can be partially transferred from one cloud to the other via the dipole-dipole interaction. This coupling leads to the formation of dark and bright superpositions of the light which experience different absorption and dispersion. We show that this feature can be exploited to filter the incident photon in such a way that only a desired path-superposition state is transmitted transparently. Finally, we generalize the analysis to the case of N coupled one-dimensional clouds. We show that within some approximations the dynamics of the propagating photon can be mapped on that of a free particle with complex mass.


title: Truncated Barsotti-Tate Groups and Displays

abstract: We define truncated displays over rings in which a prime p is nilpotent, we associate crystals to truncated displays, and we define functors from truncated displays to truncated Barsotti-Tate groups.


title: A complete characterization of connected Lie groups with the Approximation Property

abstract: We give a complete characterization of connected Lie groups with the Approximation Property for groups (AP). To this end, we introduce a strengthening of property (T), that we call property (T*), which is a natural obstruction to the AP. In order to define property (T*), we first prove that for every locally compact group G, there exists a unique left invariant mean on the space of completely bounded Fourier multipliers of G. A locally compact group G is said to have property (T*) if this mean is a weak* continuous functional. After proving that the groups SL(3,R), Sp(2,R), and the universal covering of Sp(2,R) have property (T*), we address the question which connected Lie groups have the AP. A technical problem that arises when considering this question from the point of view of the AP is that the semisimple part of the global Levi decomposition of a connected Lie group need not be closed. Because of an important permanence property of property (T*), this problem vanishes. It follows that a connected Lie group has the AP if and only if all simple factors in the semisimple part of its Levi decomposition have real rank 0 or 1. Finally, we are able to establish property (T*) for all connected simple higher rank Lie groups with finite center.


title: A Framework for Transforming Departmental Culture to Support Educational Innovation

abstract: This paper provides a research-based framework for promoting institutional change in higher education. To date, most educational change efforts have focused on relatively narrow subsets of the university system (e.g., faculty teaching practices or administrative policies) and have been largely driven by implicit change logics; both of these features have limited the success of such efforts at achieving sustained, systemic change. Drawing from the literature on organizational and cultural change, our framework encourages change agents to coordinate their activities across three key levels of the university and to ground their activities in the various change perspectives that emerge from that literature. We use examples from a change project that we have been carrying out at a large research university to illustrate how our framework can be used as a basis for planning and implementing holistic change.


title: Realizability of Tropical Curves in a Plane in the Non-Constant Coefficient Case

abstract: Let X be a plane in a torus over an algebraically closed field K, with tropicalization the matroidal fan Sigma. In this paper we present an algorithm which completely solves the question whether a given one-dimensional balanced polyhedral complex in Sigma is relatively realizable, i.e. whether it is the tropicalization of an algebraic curve, over the field of Puiseux series over K, in X. The algorithm implies that the space of all such relatively realizable curves of fixed degree is an abstract polyhedral set. In the case when X is a general plane in 3-space, we use the idea of this algorithm to prove some necessary and some sufficient conditions for relative realizability. For 1-dimensional polyhedral complexes in Sigma that have exactly one bounded edge, passing through the origin, these necessary and sufficient conditions coincide, so that they give a complete non-algorithmic solution of the relative realizability problem.


title: Plasma probe characteristics in low density hydrogen pulsed plasmas

abstract: Probe theories are only applicable in the regime where the probe's perturbation of the plasma can be neglected. However, it is not always possible to know, a priori, that a particular probe theory can be successfully applied, especially in low density plasmas. This is especially difficult in the case of transient, low density plasmas. Here, we applied probe diagnostics in combination with a 2D particle-in-cell model, to an experiment with a pulsed low density hydrogen plasma. The calculations took into account the full chamber geometry, including the plasma probe as an electrode in the chamber. It was found that the simulations reproduce the time evolution of the probe IV characteristics with good accuracy. The disagreement between the simulated and probe measured plasma density is attributed to the limited applicability of probe theory to measurements of low density pulsed plasmas. Indeed, in the case studied here, probe measurements would lead to a large overestimate of the plasma density. In contrast, the simulations of the plasma evolution and the probe characteristics do not suffer from such strict applicability limits. These studies show that probe theory cannot be justified through probe measurements.


title: Towards a Model of Systemic Change in University STEM Education

abstract: Despite numerous calls for the transformation of undergraduate STEM education, there is still a lack of successful models for creating large-scale, systemic cultural changes in STEM departments. To date, change efforts have generally focused on one of three areas: developing reflective teachers, disseminating curricula and pedagogy, or enacting institutional policy. These efforts illustrate many of the challenges of departmental change; in particular, they highlight the need for a holistic approach that integrates across all three of these levels: individual faculty, whole departments, and university policymakers. To address these challenges, as part of our campus-wide AAU-sponsored effort in STEM education transformation, we import and integrate models of change from multiple perspectives. We draw from models in organizational change, from departmental and disciplinary change in STEM education, and from efforts to support individual efforts such as the development and dissemination model. As a result, our departmental cultural change efforts are an attempt at holistic reform. We will discuss our theoretical underpinnings and ground this theory in a sample of approaches in two departments.


title: Model-based Path Integral Stochastic Control: A Bayesian Nonparametric Approach

abstract: Over the last few years, sampling-based stochastic optimal control (SOC) frameworks have shown impressive performances in reinforcement learning (RL) with applications in robotics. However, such approaches require a large amount of samples from many interactions with the physical systems. To improve learning efficiency, we present a novel model-based and data-driven SOC framework based on path integral formulation and Gaussian processes (GPs). The proposed approach learns explicit and time-varying optimal controls autonomously from limited sampled data. Based on this framework, we propose an iterative control scheme with improved applicability in higher-dimensional and more complex control tasks. We demonstrate the effectiveness and efficiency of the proposed framework using two nontrivial examples. Compared to state-of-the-art RL methods, the proposed framework features superior control learning efficiency.


title: Multilevel Monte Carlo for stochastic differential equations with small noise

abstract: We consider the problem of numerically estimating expectations of solutions to stochastic differential equations driven by Brownian motions in the commonly occurring small noise regime. We consider (i) standard Monte Carlo methods combined with numerical discretization algorithms tailored to the small noise setting, and (ii) a multilevel Monte Carlo method combined with a standard Euler-Maruyama implementation. Under the assumptions we make on the underlying model, the multilevel method combined with Euler-Maruyama is often found to be the most efficient option. Moreover, under a wide range of scalings the multilevel method is found to give the same asymptotic complexity that would arise in the idealized case where we have access to exact samples of the required distribution at a cost of $O(1)$ per sample. A key step in our analysis is to analyze the variance between two coupled paths directly, as opposed to their $L^2$ distance. Careful simulations are provided to illustrate the asymptotic results.


title: Rapid Sampling for Visualizations with Ordering Guarantees

abstract: Visualizations are frequently used as a means to understand trends and gather insights from datasets, but often take a long time to generate. In this paper, we focus on the problem of rapidly generating approximate visualizations while preserving crucial visual proper- ties of interest to analysts. Our primary focus will be on sampling algorithms that preserve the visual property of ordering; our techniques will also apply to some other visual properties. For instance, our algorithms can be used to generate an approximate visualization of a bar chart very rapidly, where the comparisons between any two bars are correct. We formally show that our sampling algorithms are generally applicable and provably optimal in theory, in that they do not take more samples than necessary to generate the visualizations with ordering guarantees. They also work well in practice, correctly ordering output groups while taking orders of magnitude fewer samples and much less time than conventional sampling schemes.


title: Connection between dark matter abundance and primordial tensor perturbations

abstract: Primordial inflation and Dark Matter (DM) could both belong to the hidden sector. It is therefore plausible that the inflaton, which drives inflation, could couple to the DM either directly or indirectly, thus providing a common origin for both luminous and non-luminous matter. We explore this interesting possibility and show that, in certain scenarios, the DM mass can be correlated with the tensor-to-scalar ratio. This correlation might provide us with a window of opportunity for unravelling the properties of DM beyond the standard freeze-out paradigm.


title: Nonequilibrium plasmons with gain in graphene

abstract: Graphene supports strongly confined transverse-magnetic sheet plasmons whose spectral characteristics depend on the energetic distribution of Dirac particles. The question arises whether plasmons can become amplified when graphene is pumped into a state of inversion. In establishing a theory for the dynamic non-equilibrium polarizability, we are able to determine the exact complex-frequency plasmon dispersion of photo-inverted graphene and study the impact of doping, collision loss, and temperature on the plasmon gain. We calculate the spontaneous emission spectra and carrier recombination rates self-consistently and compare the results with approximations based on Fermi's golden rule. Our results show that amplification of plasmons is possible under realistic conditions but inevitably competes with ultrafast spontaneous emission, which for intrinsic graphene, is a factor 5 faster than previously estimated. This work casts new light on the nature of non-equilibrium plasmons and may aid the experimental realization of active plasmonic devices based on graphene.


title: Single-color two-photon spectroscopy of Rydberg states in electric fields

abstract: Rydberg states of atomic helium with principal quantum numbers ranging from n=20 to n=100 have been prepared by non-resonance-enhanced single-color two-photon excitation from the metastable 2 {^3}S{_1} state. Photoexcitation was carried out using linearly and circularly polarized pulsed laser radiation. In the case of excitation with circularly polarized radiation, Rydberg states with azimuthal quantum number |m_{\ell}|=2 were prepared in zero electric field, and in homogeneous electric fields oriented parallel to the propagation axis of the laser radiation. In sufficiently strong electric fields, individual Rydberg-Stark states were resolved spectroscopically, highlighting the suitability of non-resonance-enhanced multiphoton excitation schemes for the preparation of long-lived high-|m_{\ell}| hydrogenic Rydberg states for deceleration and trapping experiments. Applications of similar schemes for Doppler-free excitation of positronium atoms to Rydberg states are also discussed.


title: Plastic yielding in nanocrystalline Pd-Au alloys mimics universal behavior of metallic glasses

abstract: We studied solid solution effects on the mechanical properties of nanocrystalline (NC) $\mathrm{Pd}_{\mathrm{100-x}} \mathrm{Au}_{\mathrm{x}}$ alloys ($0 \leq \mathrm{x} < 50 \mathrm{at.\%}$) at the low end of the nanoscale. Concentration has been used as control parameter to tune material properties (elastic moduli, Burgers vector, stacking fault energies) at basically unaltered microstructure (grain size $D\approx 10\mathrm{ nm}$). In stark contrast to coarse grained fcc alloys, we observe solid solution softening for increasing Au-content. The available predictions from models and theories taking explicitly into account the effect of the nanoscale microstructure on the concentration-dependent shear strength have been disproved without exception. As a consequence, it is implied that dislocation activity contributes only marginally to strength. In fact, we find a linear correlation between shear strength and shear modulus which quantitatively agrees with the universal behavior of metallic glasses discovered by Johnson and Samwer [W.L. Johnson and K. Samwer, PRL 95, 195501 (2005)].


title: Soft X-ray Pulsations in Solar Flares

abstract: The soft X-ray emissions of solar flares come mainly from the bright coronal loops at the highest temperatures normally achieved in the flare process. Their ubiquity has led to their use as a standard measure of flare occurrence and energy, although the overwhelming bulk of the total flare energy goes elsewhere. Recently Dolla et al. (2012) noted quasi-periodic pulsations (QPP) in the soft X-ray signature of the X-class flare SOL2011-02-15, as observed by the standard photometric data from the GOES (Geostationary Operational Environmental Satellite) spacecraft. In this paper we analyze the suitability of the GOES data for this kind of analysis and find them to be generally valuable after September, 2010 (GOES-15). We then extend the Dolla et al. result to a complete list of X-class flares from Cycle 24, and show that most of them (80%) display QPPs in the impulsive phase. The pulsations show up cleanly in both channels of the GOES data, making use of time-series of irradiance differences (the digital time derivative on the 2-s sampling). We deploy different techniques to characterize the periodicity of GOES pulsations, considering the red-noise properties of the flare signals, finding a range of chracteristic time scales of the QPPs for each event, but usually with no strong signature of a single period dominating in the power spectrum. The QPP may also appear on somewhat longer time scales during the later gradual phase, possibly with a greater tendency towards coherence, but the sampling noise in GOES difference data for large irradiance values (X-class flares) makes these more uncertain. We show that there is minimal phase difference between the differenced GOES energy channels, or between them and the hard X-ray variations on short time scales. During the impulsive phase the footpoints of the newly-forming flare loops may also contribute to the observed soft X-ray variations.


title: Provable Tensor Methods for Learning Mixtures of Generalized Linear Models

abstract: We consider the problem of learning mixtures of generalized linear models (GLM) which arise in classification and regression problems. Typical learning approaches such as expectation maximization (EM) or variational Bayes can get stuck in spurious local optima. In contrast, we present a tensor decomposition method which is guaranteed to correctly recover the parameters. The key insight is to employ certain feature transformations of the input, which depend on the input generative model. Specifically, we employ score function tensors of the input and compute their cross-correlation with the response variable. We establish that the decomposition of this tensor consistently recovers the parameters, under mild non-degeneracy conditions. We demonstrate that the computational and sample complexity of our method is a low order polynomial of the input and the latent dimensions.


title: Nonadiabatic nuclear dynamics of the ammonia cation studied by surface hopping classical trajectory calculations

abstract: The Landau--Zener (LZ) type classical-trajectory surface-hopping algorithm is applied to the nonadiabatic nuclear dynamics of the ammonia cation after photoionization of the ground-state neutral molecule to the excited states of the cation. The algorithm employs the recently proposed formula for nonadiabatic LZ transition probabilities derived from the adiabatic potential energy surfaces. The evolution of the populations of the ground state and the two lowest excited adiabatic states is calculated up to 200 fs. The results agree well with quantum simulations available for the first 100 fs based on the same potential energy surfaces. Four different time scales are detected for the nuclear dynamics: Ultrafast Jahn--Teller dynamics between the excited states on a 5 fs time scale; fast transitions between the excited state and the ground state within a time scale of 20 fs; relatively slow partial conversion of a first-excited-state population to the ground state within a time scale of 100 fs; and nearly constant populations after roughly 120 fs due to a dynamical equilibrium between all three states. The latter provides a possible explanation of the experimental evidence that ammonia cation is nonfluorescent.


title: Howson's property for semidirect products of semilattices by groups

abstract: An inverse semigroup $S$ is a Howson inverse semigroup if the intersection of finitely generated inverse subsemigroups of $S$ is finitely generated. Given a locally finite action $\theta$ of a group $G$ on a semilattice $E$, it is proved that $E \ast_{\theta} G$ is a Howson inverse semigroup if and only if $G$ is a Howson group. It is also shown that this equivalence fails for arbitrary actions.


title: Bi-Lipschitz homeomorphic subanalytic sets have bi-Lipschitz homeomorphic tangent cones

abstract: We prove that if there exists a bi-Lipschitz homeomorphism (not necessarily subanalytic) between two subanalytic sets, then their tangent cones are bi-Lipschitz homeomorphic. As a consequence of this result, we show that any Lipschitz regular complex analytic set, i.e any complex analytic set which is locally bi-lipschitz homeomorphic to an Euclidean ball must be smooth. Finally, we give an alternative proof of S. Koike and L. Paunescu's result about the bi-Lipschitz invariance of directional dimensions of subanalytic sets.


title: A Bayesian model selection approach for identifying differentially expressed transcripts from RNA-Seq data

abstract: Recent advances in molecular biology allow the quantification of the transcriptome and scoring transcripts as differentially or equally expressed between two biological conditions. Although these two tasks are closely linked, the available inference methods treat them separately: a primary model is used to estimate expression and its output is post-processed using a differential expression model. In this paper, both issues are simultaneously addressed by proposing the joint estimation of expression levels and differential expression: the unknown relative abundance of each transcript can either be equal or not between two conditions. A hierarchical Bayesian model builds upon the BitSeq framework and the posterior distribution of transcript expression and differential expression is inferred using Markov Chain Monte Carlo (MCMC). It is shown that the proposed model enjoys conjugacy for fixed dimension variables, thus the full conditional distributions are analytically derived. Two samplers are constructed, a reversible jump MCMC sampler and a collapsed Gibbs sampler, and the latter is found to perform best. A cluster representation of the aligned reads to the transcriptome is introduced, allowing parallel estimation of the marginal posterior distribution of subsets of transcripts under reasonable computing time. The proposed algorithm is benchmarked against alternative methods using synthetic datasets and applied to real RNA-sequencing data. Source code is available online (https://github.com/mqbssppe/cjBitSeq).


title: POPE: Post Optimization Posterior Evaluation of Likelihood Free Models

abstract: In many domains, scientists build complex simulators of natural phenomena that encode their hypotheses about the underlying processes. These simulators can be deterministic or stochastic, fast or slow, constrained or unconstrained, and so on. Optimizing the simulators with respect to a set of parameter values is common practice, resulting in a single parameter setting that minimizes an objective subject to constraints. We propose a post optimization posterior analysis that computes and visualizes all the models that can generate equally good or better simulation results, subject to constraints. These optimization posteriors are desirable for a number of reasons among which easy interpretability, automatic parameter sensitivity and correlation analysis and posterior predictive analysis. We develop a new sampling framework based on approximate Bayesian computation (ABC) with one-sided kernels. In collaboration with two groups of scientists we applied POPE to two important biological simulators: a fast and stochastic simulator of stem-cell cycling and a slow and deterministic simulator of tumor growth patterns.


title: Effects of a tilted magnetic field in a Dirac double layer

abstract: We calculate the energy spectrum of a Dirac double layer, where each layer has the Dirac electronic dispersion, in the presence of a tilted magnetic field and small interlayer tunneling. We show that the energy splitting between the Landau levels has an oscillatory dependence on the in-plane magnetic field and vanishes at a series of special tilt angles of the magnetic field. Using a semiclassical analysis, we show that these special tilt angles are determined by the Berry phase of the Dirac Hamiltonian. The interlayer tunneling conductance also exhibits an oscillatory dependence on the magnetic field tilt angle, known as the angular magnetoresistance oscillations (AMRO). Our results are applicable to graphene double layers and thin films of topological insulators.


title: Bulk Property on Cayley Tree with Smooth Boundary Condition

abstract: We study a nearest-neighbor hopping model on the Cayley tree under the smooth boundary condition with the modulation function $f_s=\sin^2[\pi s/(2M+1)]$, where $s$ is a distance from the central site, and $M$ is the number of shells on the tree. As a result of this smoothing, the particle density in the ground state becomes nearly uniform in the bulk region even when $M$ is relatively small. We compare the calculated particle density at the center with exact result on the Bethe lattice, and they show a good agreement. The calculated bond energy at the center also agrees with that on the Bethe lattice.


title: Unitary Cayley Graphs of Dedekind Domain Quotients

abstract: If $X$ is a commutative ring with unity, then the unitary Cayley graph of $X$, denoted $G_X$, is defined to be the graph whose vertex set is $X$ and whose edge set is $\{\{a,b\}\colon a-b\in X^\times\}$. When $R$ is a Dedekind domain and $I$ is an ideal of $R$ such that $R/I$ is finite and nontrivial, we refer to $G_{R/I}$ as a \emph{generalized totient graph}. We study generalized totient graphs as generalizations of the graphs $G_{\mathbb{Z}/(n)}$, which have appeared recently in the literature, sometimes under the name \emph{Euler totient Cayley graphs}. We begin by generalizing to Dedekind domains the arithmetic functions known as Schemmel totient functions, and we use one of these generalizations to provide a simple formula, for any positive integer $m$, for the number of cliques of order $m$ in a generalized totient graph. In particular, we prove that the number of cliques of order $m$ in $G_{\mathbb Z/(n)}$ is \[\prod_{k=1}^m\frac{S_{k-1}(n)}{k},\] where $S_r$ is the $r^{\text{th}}$ Schemmel totient function. We then proceed to determine many properties of generalized totient graphs such as their clique numbers, chromatic numbers, chromatic indices, clique domination numbers, and (in many, but not all cases) girths. We also determine the diameter of each component of a generalized totient graph. We correct one erroneous claim about the clique domination numbers of Euler totient Cayley graphs that has appeared in the literature and provide a counterexample to a second claim about the strong domination numbers of these graphs.


title: Transport properties of a two-dimensional electron gas dressed by light

abstract: We show theoretically that the strong interaction of a two-dimensional electron gas (2DEG) with a dressing electromagnetic field drastically changes its transport properties. Particularly, the dressing field leads to the giant increase of conductivity (which can reach thousands of percents), results in nontrivial oscillating dependence of conductivity on the field intensity, and suppresses the weak localization of 2DEG. As a consequence, the developed theory opens an unexplored way to control transport properties of 2DEG by a strong high-frequency electromagnetic field. From experimental viewpoint, this theory is applicable directly to quantum wells exposed to a laser-generated electromagnetic wave.


title: Phishing Detection in IMs using Domain Ontology and CBA - An innovative Rule Generation Approach

abstract: User ignorance towards the use of communication services like Instant Messengers, emails, websites, social networks etc. is becoming the biggest advantage for phishers. It is required to create technical awareness in users by educating them to create a phishing detection application which would generate phishing alerts for the user so that phishing messages are not ignored. The lack of basic security features to detect and prevent phishing has had a profound effect on the IM clients, as they lose their faith in e-banking and e-commerce transactions, which will have a disastrous impact on the corporate and banking sectors and businesses which rely heavily on the internet. Very little research contributions were available in for phishing detection in Instant messengers. A context based, dynamic and intelligent phishing detection methodology in IMs is proposed, to analyze and detect phishing in Instant Messages with relevance to domain ontology (OBIE) and utilizes the Classification based on Association (CBA) for generating phishing rules and alerting the victims. A PDS Monitoring system algorithm is used to identify the phishing activity during exchange of messages in IMs, with high ratio of precision and recall. The results have shown improvement by the increased percentage of precision and recall when compared to the existing methods.


title: Adaptive finite difference methods for nonlinear elliptic and parabolic partial differential equations with free boundaries

abstract: Monotone finite difference methods provide stable convergent discretizations of a class of degenerate elliptic and parabolic Partial Differential Equations (PDEs). These methods are best suited to regular rectangular grids, which leads to low accuracy near curved boundaries or singularities of solutions. In this article we combine monotone finite difference methods with an adaptive grid refinement technique to produce a PDE discretization and solver which is applied to a broad class of equations, in curved or unbounded domains which include free boundaries. The grid refinement is flexible and adaptive. The discretization is combined with a fast solution method, which incorporates asynchronous time stepping adapted to the spatial scale. The framework is validated on linear problems in curved and unbounded domains. Key applications include the obstacle problem and the one-phase Stefan free boundary problem.


title: On the formation of shocks for quasilinear wave equations

abstract: The paper is devoted to the study of shock formation of the 3-dimensional quasilinear wave equation \begin{equation}\label{Main Equation} - \big(1+3G^{\prime\prime}(0) (\partial_t\phi)^2\big)\partial^2_t \phi +\Delta\phi=0,\tag{\textbf{$\star$}} \end{equation} where $G^{\prime\prime}(0)$ is a non-zero constant. We will exhibit a family of smooth initial data and show that the foliation of the incoming characteristic hypersurfaces collapses. Similar to 1-dimensional conservational laws, we refer this specific type breakdown of smooth solutions as shock formation. Since $(\star)$ satisfies the classical null condition, it admits global smooth solutions for small data. Therefore, we will work with large data (in energy norm). Moreover, no symmetry condition is imposed on the initial datum. We emphasize the geometric perspectives of shock formations in the proof. More specifically, the key idea is to study the interplay between the following two objects: (1) the energy estimates of the linearized equations of $(\star)$; (2) the differential geometry of the Lorentzian metric $g=-\dfrac{1}{\left(1+3G^{\prime\prime}(0) (\partial_t\phi)^2\right)} d t^2+dx_1^2+dx_2^2+dx_3^2$. Indeed, the study of the characteristic hypersurfaces (implies shock formation) is the study of the null hypersurfaces of $g$. The techniques in the proof are inspired by the work \cite{Ch-Shocks} in which the formation of shocks for $3$-dimensional relativistic compressible Euler equations with small initial data is established. We also use the short pulse method which is introduced in the study of formation of black holes in general relativity in \cite{Ch-BlackHoles} and generalized in \cite{K-R-09}.


title: The role of homology in fluid vortices I: non-relativistic flow

abstract: The methods of singular and de Rham homology and cohomology are reviewed to the extent that they are applicable to the structure and motion of vortices. In particular, they are first applied to the concept of integral invariants. After a brief review of the elements of fluid mechanics, when expressed in the language of exterior differential forms and homology theory, the basic laws of vortex theory are shown to be statements that are rooted in the homology theory of integral invariants.


title: Holographic Double Diffraction of Higgs and the AdS Graviton/Pomeron

abstract: The holographic approach to double diffractive Higgs production is presented in terms of exchanging the AdS graviton/Pomeron. The goal is to provide a simple framework for central exclusive production from a dual strong coupling perspective.


title: On the physical inconsistency of a new statistical scaling symmetry in incompressible Navier-Stokes turbulence

abstract: A detailed theoretical investigation is given which demonstrates that a recently proposed statistical scaling symmetry is physically void. Although this scaling is mathematically admitted as a unique symmetry transformation by the underlying statistical equations for incompressible Navier-Stokes turbulence on the level of the functional Hopf equation, by closer inspection, however, it leads to physical inconsistencies and erroneous conclusions in the theory of turbulence. The new statistical symmetry is thus misleading in so far as it forms within an unmodelled theory an analytical result which at the same time lacks physical consistency. Our investigation will expose this inconsistency on different levels of statistical description, where on each level we will gain new insights for its non-physical transformation behavior. With a view to generate invariant turbulent scaling laws, the consequences will be finally discussed when trying to analytically exploit such a symmetry. In fact, a mismatch between theory and numerical experiment is conclusively quantified. We ultimately propose a general strategy on how to not only track unphysical statistical symmetries, but also on how to avoid generating such misleading invariance results from the outset. All the more so as this specific study on a physically inconsistent scaling symmetry only serves as a representative example within the broader context of statistical invariance analysis. In this sense our investigation is applicable to all areas of statistical physics in which symmetries get determined in order to either characterize complex dynamical systems, or in order to extract physically useful and meaningful information from the underlying dynamical process itself.

categories: physics.flu-dyn cond-mat.stat-mech math-ph math.MP


title: The Burgess inequality and the least k-th power non-residue

abstract: The Burgess inequality is the best upper bound we have for the character sum $S_{\chi}(M,N) = \sum_{M<n\le M+N} \chi(n).$ Until recently, no explicit estimates had been given for the inequality. In 2006, Booker gave an explicit estimate for quadratic characters which he used to calculate the class number of a 32-digit discriminant. McGown used an explicit estimate to show that there are no norm-Euclidean Galois cubic fields with discriminant greater than $10^{140}$. Both of their explicit estimates are on restricted ranges. In this paper we prove an explicit estimate that works for any $M$ and $N$. We also improve McGown's estimates in a slightly narrower range, getting explicit estimates for characters of any order. We apply the estimates to the question of how large must a prime $p$ be to ensure that there is a $k$-th power non-residue less than $p^{1/6}$.


title: Modified chalcogenide glasses for optical device applications

abstract: This thesis focuses on two different, but complementary, aspects of the modification of gallium lanthanum sulphide (GLS) glasses. Firstly the addition of transition metal ions as dopants is examined and their potential for use as active optical materials is explored. It is also argued that the spectroscopic analysis of transition metal ions is a useful tool for evaluating the local environment of their host. Secondly femtosecond (fs) laser modification of GLS is investigated as a method for waveguide formation.


title: Homomorphisms of connectome graphs

abstract: We propose to study homomorphisms of connectome graphs. Homomorphisms can be studied as sequences of elementary homomorphisms - folds, which identify pairs of vertices. Several fold types are defined. Initial computation results for some connectome graphs are described.


title: Optimization of a Non-arsenic Iron-based Superconductor for Wire Fabrication

abstract: We report on the optimization of synthesis of iron-selenide (non-arsenic) superconducting powders that are based on '122' composition, with optimal Tc = 38 K and Jc = 10^5 A/cm2 (4 K). We also report on the wire proof-of concept for these materials, by producing ~ 40 ft of wire that produce Ic. The 122 selenides are more difficult to synthesize and have more complex crystal structures compared to '11' selenides (FeSe and FeSe1-xTex), but they do offer higher Tc and might provoke a natural extension for 11 work.


title: An Anti-Ramsey Problem Concerning Complete Bipartite Graphs

abstract: We consider quadruples of positive integers $(a,b,m,n)$ with $a\leq b$ and $m\leq n$ such that any proper edge-coloring of the complete bipartite graph $K_{m,n}$ contains a rainbow $K_{a,b}$ subgraph. We show that any such quadruple with $a\leq m$ and $n>(a^2-a+1)(b-1)$ satisfies this property. We also show that the quadruple $(2,3,3,6)$ satisfies this property. We end with a conjecture.


title: Covers in Partitioned Intersecting Hypergraphs

abstract: Given an integer $r$ and a vector $\vec{a}=(a_1, \ldots ,a_p)$ of positive numbers with $\sum_{i \le p} a_i=r$, an $r$-uniform hypergraph $H$ is said to be $\vec{a}$-partitioned if $V(H)=\bigcup_{i \le p}V_i$, where the sets $V_i$ are disjoint, and $|e \cap V_i|=a_i$ for all $e \in H,~~i \le p$. A $\vec{1}$-partitioned hypergraph is said to be $r$-partite. Let $t(\vec{a})$ be the maximum, over all intersecting $\vec{a}$-partitioned hypergraphs $H$, of the minimal size of a cover of $H$. A famous conjecture of Ryser is that $t(\vec{1})\le r-1$. Tuza conjectured that if $r>2$ then $t(\vec{a})=r$ for every two components vector $\vec{a}=(a,b)$. We prove this conjecture whenever $a\neq b$, and also for $\vec{a}=(2,2)$ and $\vec{a}=(4,4)$.


title: Decomposition theorems for a generalization of the holonomy Lie algebra of an arrangement

abstract: In [7, Papadima and Suciu, When does the associated graded Lie algebra of an arrangement group decompose? Comment. Math. Helv. {\bf 81:4} (2006), 859--875] it is proved that the holonomy Lie algebra of an arrangement of hyperplanes through origo decomposes as a direct product of Lie algebras in degree at least two if and only if a certain (computable) condition is fulfilled. We prove similar results for a class of Lie algebras which is a generalization of the holonomy Lie algebras. The proof methods are the same as in [7].


title: Is the log-law a first principle result from Lie-group invariance analysis?

abstract: The invariance method of Lie-groups in the theory of turbulence carries the high expectation of being a first principle method for generating statistical scaling laws. The purpose of this comment is to show that this expectation has not been met so far. In particular for wall-bounded turbulent flows, the prospects for success are not promising in view of the facts we will present herein. Although the invariance method of Lie-groups is able to generate statistical scaling laws for wall-bounded turbulent flows, like the log-law for example, these invariant results yet not only fail to fulfil the basic requirements for a first principle result, but also are strongly misleading. The reason is that not the functional structure of the log-law itself is misleading, but that its invariant Lie-group based derivation yielding this function is what is misleading. By revisiting the study of Oberlack (2001) we will demonstrate that all Lie-group generated scaling laws derived therein do not convince as first principle solutions. Instead, a rigorous derivation reveals complete arbitrariness rather than uniqueness in the construction of invariant turbulent scaling laws. Important to note here is that the key results obtained in Oberlack (2001) are based on several technical errors, which all will be revealed, discussed and corrected. The reason and motivation why we put our focus solely on Oberlack (2001) is that it still marks the core study and central reference point when applying the method of Lie-groups to turbulence theory. Hence it is necessary to shed the correct light onto that study. Nevertheless, even if the method of Lie-groups in its full extent is applied and interpreted correctly, strong natural limits of this method within the theory of turbulence exist, which, as will be finally discussed, constitute insurmountable obstacles in the progress of achieving a significant breakthrough.


title: Superconductivity in Semimetallic Bi3O2S3

abstract: Here we report further investigation on the thermodynamic, and transport properties, and the first assessment of theoretical calculations for the BiS2-layered Bi3O2S3 superconductor. The polycrystalline sample is synthesized with the superconducting transition temperature of Tc_onset = 5.75 K, that drops to 3.3 K by applying hydrostatic pressure of 6 kbar. Density-of-states calculations give substantial hybridization between Bi, O and S, with the largest Bi component at DOS, which supports the idea that BiS2 layer is relevant for producing electron-phonon coupling. The electron charge carrier concentration of only 1.5 x 10^19 cm-3 for Bi3O2S3 is additionally suggestive of strong electron-phonon interaction in Bi-O-S system. The analysis of Seebeck coefficient results strongly suggests that Bi3O2S3 is a semimetal. In fact, the semimetallic or narrow band gap behavior may be general in the BiS2-layered class of materials, including that of Bi4O4S3.


title: Generalizing $p$-goodness to ordered graphs

abstract: It is known that the connected graphs that are $p$-good for all $p$ are the trees. In this paper, we ask which ordered trees are order-$p$-good for all $p$.


title: On 2-powerfully Perfect Numbers in Three Quadratic Rings

abstract: Using an extension of the abundancy index to imaginary quadratic rings with unique factorization, we define what we call $n$-powerfully perfect numbers in these rings. This definition serves to extend the concept of perfect numbers that have been defined and studied in the integers. We investigate the properties of $2$-powerfully perfect numbers in the rings $\mathcal O_{\mathbb{Q}(\sqrt{-1})}$, $\mathcal O_{\mathbb{Q}(\sqrt{-2})}$, and $\mathcal O_{\mathbb{Q}(\sqrt{-7})}$, the three imaginary quadratic rings with unique factorization in which $2$ is not a prime.


title: Local integrals of motion and the logarithmic lightcone in many-body localized systems

abstract: We propose to define full many-body localization in terms of the recently introduced integrals of motion[Chandran et al., arXiv:1407.8480], which characterize the time-averaged response of the system to a local perturbation. The quasi-locality of such integrals of motion implies an effective lightcone that grows logarithmically in time. This subsequently implies that (i) the average entanglement entropy can grow at most logarithmically in time for a global quench from a product state, and (ii) with high probability, the time evolution of a local operator for a time interval $|t|$ can be classically simulated with a resource that scales polynomially in $|t|$.


title: Dynamical Generation of a Repulsive Vector Contribution to the Quark Pressure

abstract: Lattice QCD results for the coefficient $c_2$ appearing in the Taylor expansion of the pressure show that this quantity raises with the temperature towards the Stefan-Boltzmann limit. On the other hand, model approximations predict that when a vector repulsion, parametrized by $G_V$, is present this coefficient reaches a maximum just after $T_c$ and then deviates from the lattice predictions. Recently, this discrepancy has been used as a guide to constrain the (presently unknown) value of $G_V$ within the framework of effective models at large-$N_c$ (LN). In the present investigation we show that, due to finite $N_c$ effects, $c_2$ may also develop a maximum even when $G_V=0$ since a vector repulsive term can be dynamically generated by exchange type of radiative corrections. Here we apply the the Optimized Perturbation Theory (OPT) method to the two flavor Polyakov--Nambu--Jona-Lasinio model (at $G_V=0$) and compare the results with those furnished by lattice simulations an by the LN approximation at $G_V=0$ and also at $G_V \ne 0$. The OPT numerical results for $c_2$ are impressively accurate for $T \lesssim 1.2\, T_c$ but, as expected, predict that this quantity develops a maximum at high-$T$. After identifying the mathematical origin of this extremum we argue that such a discrepant behavior may naturally arise within these effective quark models (at $G_V=0$) whenever the first $1/N_c$ corrections are taken into account. We then interpret this hypothesis as an indication that beyond the large-$N_c$ limit the correct high temperature (perturbative) behavior of $c_2$ will be faithfully described by effective models only if they also mimic the asymptotic freedom phenomenon.


title: Spectrophotometric support of spectral observations with the telescope BTA

abstract: We report on the development of the medium-resolution spectrophotometer to accompany (in real time) spectral observations performed with the high resolution spectrographs on the 6-m telescope BTA.


title: The Computational Complexity of Structure-Based Causality

abstract: Halpern and Pearl introduced a definition of actual causality; Eiter and Lukasiewicz showed that computing whether X=x is a cause of Y=y is NP-complete in binary models (where all variables can take on only two values) and\ Sigma_2^P-complete in general models. In the final version of their paper, Halpern and Pearl slightly modified the definition of actual cause, in order to deal with problems pointed by Hopkins and Pearl. As we show, this modification has a nontrivial impact on the complexity of computing actual cause. To characterize the complexity, a new family D_k^P, k= 1, 2, 3, ..., of complexity classes is introduced, which generalizes the class DP introduced by Papadimitriou and Yannakakis (DP is just D_1^P). %joe2 %We show that the complexity of computing causality is $\D_2$-complete %under the new definition. Chockler and Halpern \citeyear{CH04} extended the We show that the complexity of computing causality under the updated definition is $D_2^P$-complete. Chockler and Halpern extended the definition of causality by introducing notions of responsibility and blame. The complexity of determining the degree of responsibility and blame using the original definition of causality was completely characterized. Again, we show that changing the definition of causality affects the complexity, and completely characterize it using the updated definition.


title: Witnessing gas mixing in the metal distribution of the Hickson Compact Group HCG 31

abstract: We present, for the first time, direct evidence that, in a merger of disk galaxies, the preexisting central metallicities will mix as a result of gas being transported in the merger interface region, along the line that joins the two coalescing nuclei. This is shown using detailed 2D kinematics as well as metallicity measurements for the nearby ongoing merger in the center of the compact group HCG 31. We focus on the emission line gas, which is extensive in the system. The two coalescing cores display similar oxygen abundances while in between the two nuclei the metallicity changes smoothly from one nucleus to the other, indicating a mix of metals in this region, which is confirmed by the high resolution H$\alpha$ kinematics (R=45900). This nearby system is especially important because it involves the merging of two fairly low mass and clumpy galaxies (LMC-like galaxies), making it an important system for comparison with high redshift galaxies.


title: Hierarchical Mixture-of-Experts Model for Large-Scale Gaussian Process Regression

abstract: We propose a practical and scalable Gaussian process model for large-scale nonlinear probabilistic regression. Our mixture-of-experts model is conceptually simple and hierarchically recombines computations for an overall approximation of a full Gaussian process. Closed-form and distributed computations allow for efficient and massive parallelisation while keeping the memory consumption small. Given sufficient computing resources, our model can handle arbitrarily large data sets, without explicit sparse approximations. We provide strong experimental evidence that our model can be applied to large data sets of sizes far beyond millions. Hence, our model has the potential to lay the foundation for general large-scale Gaussian process research.


title: Computoser - rule-based, probability-driven algorithmic music composition

abstract: This paper presents the Computoser hybrid probability/rule based algorithm for music composition (http://computoser.com) and provides a reference implementation. It addresses the issues of unpleasantness and lack of variation exhibited by many existing approaches by combining the two methods (basing the parameters of the rules on data obtained from preliminary analysis). A sample of 500+ musical pieces was analyzed to derive probabilities for musical characteristics and events (e.g. scale, tempo, intervals). The algorithm was constructed to produce musical pieces using the derived probabilities combined with a large set of composition rules, which were obtained and structured after studying established composition practices. Generated pieces were published on the Computoser website where evaluation was performed by listeners. The feedback was positive (58.4% approval), asserting the merits of the undertaken approach. The paper compares this hybrid approach to other approaches to algorithmic composition and presents a survey of the pleasantness of the resulting music.


title: On Sparsely Schemmel Totient Numbers

abstract: For each positive integer $r$, let $S_r$ denote the $r^{th}$ Schemmel totient function, a multiplicative arithmetic function defined by \[S_r(p^{\alpha})=\begin{cases} 0, & \mbox{if } p\leq r; \\ p^{\alpha-1}(p-r), & \mbox{if } p>r \end{cases}\] for all primes $p$ and positive integers $\alpha$. The function $S_1$ is simply Euler's totient function $\phi$. Masser and Shiu have established several fascinating results concerning sparsely totient numbers, positive integers $n$ satisfying $\phi(n)<\phi(m)$ for all integers $m>n$. We define a sparsely Schemmel totient number of order $r$ to be a positive integer $n$ such that $S_r(n)>0$ and $S_r(n)<S_r(m)$ for all $m>n$ with $S_r(m)>0$. We then generalize some of the results of Masser and Shiu.


title: Quantum Monte Carlo methods for nuclear physics

abstract: Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

categories: nucl-th astro-ph.SR cond-mat.quant-gas nucl-ex quant-ph


title: The principle of stationary nonconservative action for classical mechanics and field theories

abstract: We further develop a recently introduced variational principle of stationary action for problems in nonconservative classical mechanics and extend it to classical field theories. The variational calculus used is consistent with an initial value formulation of physical problems and allows for time-irreversible processes, such as dissipation, to be included at the level of the action. In this formalism, the equations of motion are generated by extremizing a nonconservative action $\mathcal{S}$, which is a functional of a doubled set of degrees of freedom. The corresponding nonconservative Lagrangian contains a potential $K$ which generates nonconservative forces and interactions. Such a nonconservative potential can arise in several ways, including from an open system interacting with inaccessible degrees of freedom or from integrating out or coarse-graining a subset of variables in closed systems. We generalize Noether's theorem to show how Noether currents are modified and no longer conserved when $K$ is non-vanishing. Consequently, the nonconservative aspects of a physical system are derived solely from $K$. We show how to use the formalism with examples of nonconservative actions for discrete systems including forced damped harmonic oscillators, radiation reaction on an accelerated charge, and RLC circuits. We present examples for nonconservative classical field theories. Our approach naturally allows for irreversible thermodynamic processes to be included in an unconstrained variational principle. We present the nonconservative action for a Navier-Stokes fluid including the effects of viscous dissipation and heat diffusion, as well as an action that generates the Maxwell model for viscoelastic materials, which can be easily generalized to more realistic rheological models. We show that the nonconservative action can be derived as the classical limit of a more complete quantum theory.

categories: math-ph astro-ph.HE cond-mat.stat-mech math.MP physics.flu-dyn physics.plasm-ph


title: Resonances in ultracold dipolar atomic and molecular gases

abstract: A previously developed approach for the numerical treatment of two particles that are confined in a finite optical-lattice potential and interact via an arbitrary isotropic interaction potential has been extended to incorporate an additional anisotropic dipole-dipole interaction. The interplay of a model but realistic short-range Born-Oppenheimer potential and the dipole-dipole interaction for two confined particles is investigated. A variation of the strength of the dipole-dipole interaction leads to diverse resonance phenomena. In a harmonic confinement potential some resonances show similarities to $s$-wave scattering resonances while in an anharmonic trapping potential like the one of an optical lattice inelastic confinement-induced dipolar resonances occur. The latter are due to a coupling of the relative and center-of-mass motion caused by the anharmonicity of the external confinement.


title: Clique-Relaxed Competitive Graph Coloring

abstract: We investigate a variation of the graph coloring game, as studied in [2]. In the original coloring game, two players, Alice and Bob, alternate coloring vertices on a graph with legal colors from a fixed color set, where a color {\alpha} is legal for a vertex if said vertex has no neighbors colored {\alpha}. Other variations of the game change this definition of a legal color. For a fixed color set, Alice wins the game if all vertices are colored when the game ends, while Bob wins if there is a point in the game in which a vertex cannot be assigned a legal color. The least number of colors needed for Alice to have a winning strategy on a graph G is called the game chromatic number of G, and is denoted \c{hi}g(G). A well studied variation is the d-relaxed coloring game [5] in which a legal coloring of a graph G is defined as any assignment of colors to V (G) such that the subgraph of G induced by any color class has maximum degree d. We focus on the k-clique-relaxed n-coloring game. A k-clique-relaxed n-coloring of a graph G is an n-coloring in which the subgraph of G induced by any color class has maximum clique size k or less. In other words, a k-clique-relaxed n-coloring of G is an assignment of n colors to V (G) in which there are no monochromatic (k + 1)-cliques.


title: Matrix models, Toeplitz determinants and recurrence times for powers of random unitary matrices

abstract: The purpose of this article is to study the eigenvalues $u_1^{\, t}=e^{it\theta_1},\dots,u_N^{\,t}=e^{it\theta_N}$ of $U^t$ where $U$ is a large $N\times N$ random unitary matrix and $t>0$. In particular we are interested in the typical times $t$ for which all the eigenvalues are simultaneously close to $1$ in different ways thus corresponding to recurrence times in the issue of quantum measurements. Our strategy consists in rewriting the problem as a random matrix integral and use loop equations techniques to compute the first orders of the large $N$ asymptotic. We also connect the problem to the computation of a large Toeplitz determinant whose symbol is the characteristic function of several arc segments of the unit circle. In particular in the case of a single arc segment we recover Widom's formula. Eventually we explain why the first return time is expected to converge towards an exponential distribution when $N$ is large. Numeric simulations are provided along the paper to illustrate the results.


title: Letter of Intent to Construct a nuPRISM Detector in the J-PARC Neutrino Beamline

abstract: As long-baseline neutrino experiments enter the precision era, the difficulties associated with understanding neutrino interaction cross sections on atomic nuclei are expected to limit experimental sensitivities to oscillation parameters. In particular, the ability to relate experimental observables to neutrino energy in previous experiments has relied solely on theoretical models of neutrino-nucleus interactions, which currently suffer from very large theoretical uncertainties. By observing charged current $\nu_\mu$ interactions over a continuous range of off-axis angles from 1 to 4 degrees, the nuPRISM water Cherenkov detector can provide a direct measurement of the far detector lepton kinematics for any given set of oscillation parameters, which largely removes neutrino interaction modeling uncertainties from T2K oscillation measurements. This naturally provides a direct constraint on the relationship between lepton kinematics and neutrino energy. In addition, nuPRISM is a sensitive probe of sterile neutrino oscillations with multiple energy spectra, which provides unique constraints on possible background-related explanations of the MiniBooNE anomaly. Finally, high-precision measurements of neutrino cross sections on water are possible, including $\nu_e$ measurements and the first ever measurements of neutral current interactions as a function of neutrino energy. The nuPRISM detector also benefits the proposed Hyper-Kamiokande project. A demonstration that neutrino interaction uncertainties can be controlled will be important to understanding the physics reach of Hyper-K. In addition, nuPRISM will provide an easily accessible prototype detector for many of the new hardware components currently under consideration for Hyper-K. The following document presents the configuration, physics impact, and preliminary cost estimates for a nuPRISM detector in the J-PARC neutrino beamline.


title: Critical behavior and out-of-equilibrium dynamics of a two-dimensional Ising model with dynamic couplings

abstract: We study the critical behavior and the out-of-equilibrium dynamics of a two-dimensional Ising model with non-static interactions. In our model, bonds are dynamically changing according to a majority rule depending on the set of closest neighbors of each spin pair, which prevents the system from ordering in a full ferromagnetic or antiferromagnetic state. Using a parallel-tempering Monte Carlo algorithm, we find that the model undergoes a continuous phase transition at finite temperature, which belongs to the Ising universality class. The properties of the bond structure and the ground-state entropy are also studied. Finally, we analyze the out-of-equilibrium dynamics which displays typical glassy characteristics at a temperature well below the critical one.


title: Filters and Matrix Factorization

abstract: We give a number of explicit matrix-algorithms for analysis/synthesis in multi-phase filtering; i.e., the operation on discrete-time signals which allow a separation into frequency-band components, one for each of the ranges of bands, say $N$, starting with low-pass, and then corresponding filtering in the other band-ranges. If there are $N$ bands, the individual filters will be combined into a single matrix action; so a representation of the combined operation on all $N$ bands by an $N \times N$ matrix, where the corresponding matrix-entries are periodic functions; or their extensions to functions of a complex variable. Hence our setting entails a fixed $N \times N$ matrix over a prescribed algebra of functions of a complex variable. In the case of polynomial filters, the factorizations will always be finite. A novelty here is that we allow for a wide family of non-polynomial filter-banks. Working modulo $N$ in the time domain, our approach also allows for a natural matrix-representation of both down-sampling and up-sampling. The implementation encompasses the combined operation on input, filtering, down-sampling, transmission, up-sampling, an action by dual filters, and synthesis, merges into a single matrix operation. Hence our matrix-factorizations break down the global filtering-process into elementary steps. To accomplish this, we offer a number of adapted matrix factorization-algorithms, such that each factor in our product representation implements in a succession of steps the filtering across pairs of frequency-bands; and so it is of practical significance in implementing signal processing, including filtering of digitized images. Our matrix-factorizations are especially useful in the case of the processing a fixed, but large, number of bands.


title: On Schemmel Nontotient Numbers

abstract: For each positive integer $r$, let $S_r$ denote the $r^{th}$ Schemmel totient function, a multiplicative arithmetic function defined by \[S_r(p^{\alpha})=\begin{cases} 0, & \mbox{if } p\leq r; \\ p^{\alpha-1}(p-r), & \mbox{if } p>r \end{cases}\] for all primes $p$ and positive integers $\alpha$. The function $S_1$ is simply Euler's totient function $\phi$. We define a Schemmel nontotient number of order $r$ to be a positive integer that is not in the range of the function $S_r$. In this paper, we modify several proofs due to Zhang in order to illustrate how many of the results currently known about nontotient numbers generalize to results concerning Schemmel nontotient numbers. We also invoke Zsigmondy's Theorem in order to generalize a result due to Mendelsohn.


title: Strong micro-macro entanglement from a weak cross-Kerr nonlinearity

abstract: We study the entanglement generated by a weak cross-Kerr nonlinearity between two initial coherent states, one of which has an amplitude close to the single-photon level, while the other one is macroscopic. We show that strong micro-macro entanglement is possible for weak phase shifts by choosing the amplitude of the macroscopic beam sufficiently large. We analyze the effects of loss and discuss possible experimental demonstrations of the micro-macro entanglement based on homodyne tomography and on a new entanglement witness.


title: The (black hole)-bulge mass scaling relation at low masses

abstract: Several recent papers have reported on the occurrence of active galactic nuclei (AGN) containing under-massive black holes relative to a linear scaling relation between black hole mass (M_bh) and host spheroid stellar mass (M_sph,*). Dramatic revisions to the M_bh-M_sph,* and M_bh-L_sph relations, based on samples containing predominantly inactive galaxies, have however recently identified a new steeper relation at M_bh < (2-10)x10^8 M_Sun, roughly corresponding to M_sph,* < (0.3-1)x10^{11} M_Sun. We show that this steeper, quadratic-like M_bh-M_sph,* relation defined by the Sersic galaxies, i.e. galaxies without partially depleted cores, roughly tracks the apparent offset of the AGN having 10^5 < M_bh/M_Sun < 0.5x10^8. That is, these AGN are not randomly offset with low black hole masses, but also follow a steeper (non-linear) relation. As noted by Busch et al., confirmation or rejection of a possible AGN offset from the steeper M_bh-M_sph,* relation defined by the Sersic galaxies will benefit from improved stellar mass-to-light ratios for the spheroids hosting these AGN. Several implications for formation theories are noted. Furthermore, reasons for possible under- and over-massive black holes, the potential existence of intermediate mass black holes (<10^5 M_Sun), and the new steep (black hole)--(nuclear star cluster) relation, M_bh ~ (M_nc)^{2.7+/-0.7}, are also discussed.


title: Eigenvalue Density of the Doubly Correlated Wishart Model: Exact Results

abstract: Data sets collected at different times and different observing points can possess correlations at different times $and$ at different positions. The doubly correlated Wishart model takes both into account. We calculate the eigenvalue density of the Wishart correlation matrices using supersymmetry. In the complex case we obtain a new closed form expression which we compare to previous results in the literature. In the more relevant and much more complicated real case we derive an expression for the density in terms of a fourfold integral. Finally, we calculate the density in the limit of large correlation matrices.


title: Two Modes of Partially Screened Gap

abstract: The analysis of X-ray observations suggest an ultrastrong ($B\gtrsim 10^{14} \,{\rm G}$) surface magnetic field at the polar cap of pulsars (Szary, 2013). On the other hand, the temperature of the polar caps is about a few millions Kelvin. Based on these two facts we use the Partially Screened Gap (PSG) model to describe the Inner Acceleration Region (IAR). The PSG model assumes that the temperature of the actual polar cap is equal to the so-called critical value, i.e. the temperature at which the outflow of thermal ions from the surface screens the gap completely. We have found that, depending on the conditions above the polar cap, the generation of high energetic photons in IAR can be caused either by Curvature Radiation (CR) or by Inverse Compton Scattering (ICS). Completely different properties of both processes result in two different scenarios of breaking the acceleration gap: the so-called PSG-off mode for the gap dominated by CR and the PSG-on mode for the gap dominated by ICS. The existence of two different mechanisms of gap breakdown naturally explains the mode-changing phenomenon. Different characteristics of plasma generated in the acceleration region for both modes also explain the pulse nulling phenomenon.


title: Distinguishing Di-jet Resonances at the LHC

abstract: Anticipating that a di-jet resonance could be discovered at the 14 TeV LHC, we present two different strategies to reveal the nature of such a particle; in particular to discern whether it is a quark-antiquark (qqbar), quark-gluon (qg), or gluon-gluon (gg) resonance. The first method relies on the color discriminant variable, which can be calculated at the LHC from the measurements of the di-jet signal cross section, the resonance mass and the resonance width. Including estimated statistical uncertainties and experimental resolution, we show that a qg excited quark resonance can be efficiently distinguished from either a qqbar coloron or a gg color-octet scalar using the color discriminant variable at LHC-14. The second strategy is based on the study of the energy profiles of the two leading jets in the di-jet channel. Including statistical uncertainties in the signal and the QCD backgrounds, we show that one can distinguish, in a model-independent way, between gg, qg, and qqbar resonances; an evaluation of systematic uncertainties in the measurement of the jet energy profile will require a detailed detector study once sufficient 14 TeV di-jet data is in hand.


title: Scheduling with two non-unit task lengths is NP-complete

abstract: The non-preemptive job scheduling problem with release times and deadlines on a single machine is fundamental to many scheduling problems. We parameterize this problem by the set of job lengths the jobs can have. The case where all job lengths are identical is known to be solvable in polynomial time. We prove that the problem with two job lengths is NP-complete, except for the case in which the short jobs have unit job length, which was already known to be efficiently solvable. The proof uses a reduction from satisfiability to an auxiliary scheduling problem that includes a set of paired jobs that each have both an early and a late deadline, and of which at least one should be scheduled before the early deadline. This reduction is enabled by not only these pairwise dependencies between jobs, but also by dependencies introduced by specifically constructed sets of jobs which have deadlines close to each other. The auxiliary scheduling problem in its turn can be reduced to the scheduling problem with two job lengths by representing each pair of jobs with two deadlines by four different jobs.


title: N-valid trees in wavelet theory on Vilenkin groups

abstract: We consider a class of $(N,M)$-elementary step functions on the $p$-adic Vilenkin group. We prove that $(N,M)$-elementary step function generates a MRA on $p$-adic Vilenkin group iff it is generated by a special $N$-valid rooted tree on the set of vertices $\{0,1,\dots p-1\}$ with the vector $(0,...,0)\in \mathbb Z^N$ as a root. Bibliography: 15 titles.


title: Properties of Holographic Mesons on Dense Medium

abstract: We study the energy dispersions of holographic light mesons and their decay constants on dense nuclear medium. As the spatial momenta of mesons along the boundary direction increase, both observables of the mesons not only increase but also split according to the isospin charges. The decay constant of the negative meson is more large than that of the positive meson of the same type due to the chemical potentials of the background nucleons.


title: On the number of active links in random wireless networks

abstract: This paper presents results on the typical number of simultaneous point-to-point transmissions above a minimum rate that can be sustained in a network with $n$ transmitter-receiver node pairs when all transmitting nodes can potentially interfere with all receivers. In particular we obtain a scaling law when the fading gains are independent Rayleigh distributed random variables and the transmitters over different realizations are located at the points of a stationary Poisson field in the plane. We show that asymptotically with probability approaching 1, the number of simultaneous transmissions (links that can transmit at greater than a minimum rate) is of the order of $O(n^{\frac{1}{4}})$. These asymptotic results are confirmed from simulations.


title: Phonon thermal conductivity by non-local non-equilibrium molecular dynamics

abstract: Non-equilibrium (NE) molecular dynamics (MD), or NEMD, gives a "direct" simulation of thermal conductivity kappa. Heat H(x) is added and subtracted in equal amounts at different places x. After steady state is achieved, the temperature T(x) is found by averaging over finite sections. Usually the aim is to extract a value of dT/dx from a place distant from sources and sinks of heat. This yields an effective kappa(L) for the thermal conductivity, L being the system size. The result is then studied as a function of L, to extract the bulk limit kappa. Here instead, our heat is H(x)~sin(qx), where q=2pi/L. This causes a steady-state temperature T_0 + Delta T sin(2pi x/L). A thermal conductivity kappa(q) is extracted, which is well converged at the chosen q (or L). Bulk conductivity kappa requires taking the q to 0 limit. The method is tested for liquid and crystalline argon. One advantage is reduced computational noise at a given total MD run time. Another advantage is that kappa(q) has a more physical meaning than kappa(L). It can be easily studied using Peierls-Boltzmann transport theory. New formulas for kappa(q) in simplified Debye-type models give new insight about extrapolation to q to 0 or 1/L to 0. In particular, it is shown that kappa(L$ is unlikely to behave as kappa -C/L, and much more likely to behave as kappa-C'/sqrt(L). Convergence problems encountered in computational cells with very large aspect ratios L(parallel)/L(perp) are also analyzed. Some details are contained in the "Supplemental Material" file.


abstract: We derive a family of linear inference algorithms that generalize existing graph-based label propagation algorithms by allowing them to propagate generalized assumptions about "attraction" or "compatibility" between classes of neighboring nodes (in particular those that involve heterophily between nodes where "opposites attract"). We thus call this formulation Semi-Supervised Learning with Heterophily (SSLH) and show how it generalizes and improves upon a recently proposed approach called Linearized Belief Propagation (LinBP). Importantly, our framework allows us to reduce the problem of estimating the relative compatibility between nodes from partially labeled graph to a simple optimization problem. The result is a very fast algorithm that -- despite its simplicity -- is surprisingly effective: we can classify unlabeled nodes within the same graph in the same time as LinBP but with a superior accuracy and despite our algorithm not knowing the compatibilities.


title: Fermi Orbital Derivatives in Self-Interaction Corrected Density Functional Theory: Applications to Closed Shell Atoms

abstract: A recent modification of the Perdew-Zunger self-interaction-correction (SIC) to the density-functional formalism (Pederson, Ruzsinszky, Perdew) has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Lowdin orthonormalized Fermi-orbitals (Luken et al) which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested here on atoms. Total energies and ionization energies in closed-shell atoms, where correlation is less important, using the PW92 LDA functional are in very good to excellent agreement with experiment and non-relativistic Quantum-Monte-Carlo (QMC) results.


title: Analysis of Average Travel Time for Stateless Opportunistic Routing Techniques

abstract: Wireless network applications, such as, searching, routing, self stabilization and query processing can be modeled as random walks on graphs. Stateless Opportunistic routing technique is a robust distributed routing technique based on random walk approach, where nodes transfer the packets to one of their direct neighbors uniformly, until the packets reach their destinations. Simplicity in execution, fault tolerance, low overhead and robustness to topology changes made it more suitable to wireless sensor networks scenarios. But the main disadvantage of stateless opportunistic routing is estimating and studying the effect of network parameters on the packet latency. In this work, we derived the analytical expressions for mean latency or average packet travel time for r-nearest neighbor cycle, r-nearest neighbor torus networks. Further, we derived the generalized expression for mean latency for m-dimensional r- nearest neighbor torus networks and studied the effect of number of nodes, nearest neighbors and network dimension on average packet travel time.


title: Sequential Hypothesis Tests for Adaptive Locality Sensitive Hashing

abstract: All pairs similarity search is a problem where a set of data objects is given and the task is to find all pairs of objects that have similarity above a certain threshold for a given similarity measure-of-interest. When the number of points or dimensionality is high, standard solutions fail to scale gracefully. Approximate solutions such as Locality Sensitive Hashing (LSH) and its Bayesian variants (BayesLSH and BayesLSHLite) alleviate the problem to some extent and provides substantial speedup over traditional index based approaches. BayesLSH is used for pruning the candidate space and computation of approximate similarity, whereas BayesLSHLite can only prune the candidates, but similarity needs to be computed exactly on the original data. Thus where ever the explicit data representation is available and exact similarity computation is not too expensive, BayesLSHLite can be used to aggressively prune candidates and provide substantial speedup without losing too much on quality. However, the loss in quality is higher in the BayesLSH variant, where explicit data representation is not available, rather only a hash sketch is available and similarity has to be estimated approximately. In this work we revisit the LSH problem from a Frequentist setting and formulate sequential tests for composite hypothesis (similarity greater than or less than threshold) that can be leveraged by such LSH algorithms for adaptively pruning candidates aggressively. We propose a vanilla sequential probability ration test (SPRT) approach based on this idea and two novel variants. We extend these variants to the case where approximate similarity needs to be computed using fixed-width sequential confidence interval generation technique.


title: Raman microscopy as a defect microprobe for hydrogen bonding characterization in materials used in fusion applications

abstract: We present the Raman microscopy ability to detect and characterize the way hydrogen is bonded with elements that will be used for ITER's plasma facing components. For this purpose we first use hydrogenated amorphous carbon samples, formed subsequently to plasma-wall interactions (hydrogen implantation, erosion, deposition...) occurring inside tokamaks, to demonstrate how this technique can be used to retrieve useful information. We pay attention in identifying which spectroscopic parameters are sensitive to the local structure (sp 3 /sp 2) and which gives information on the hydrogen content using isothermal and linear temperature ramp studies on reference samples produced by plasma enhanced chemical vapor deposition. We then focus on the possibility to use this fast, non-destructive and non-contact technique to characterize the influence of hydrogen isotope implantation in few nanometers of graphite and beryllium as C is still used in the JT-60 tokamak and Be is used in JET and will be used as plasma-facing component in the future reactor ITER. We also pay attention on implantation in tungsten oxide which may be formed accidently in the machine.


title: Unitary Multiperfect Numbers in Certain Quadratic Rings

abstract: A unitary divisor $c$ of a positive integer $n$ is a positive divisor of $n$ that is relatively prime to $\displaystyle{\frac{n}{c}}$. For any integer $k$, the function $\sigma_k^*$ is a multiplicative arithmetic function defined so that $\sigma_k^*(n)$ is the sum of the $k^{th}$ powers of the unitary divisors of $n$. We provide analogues of the functions $\sigma_k^*$ in imaginary quadratic rings that are unique factorization domains. We then explore properties of what we call $n$-powerfully unitarily $t$-perfect numbers, analogues of the unitary multiperfect numbers that have been defined and studied in the integers. We end with a list of several opportunities for further research.


title: The X-ray Flux Distribution of Sagittarius A* as Seen by Chandra

abstract: We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-\xi},$ $\xi=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $\sigma=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism.


title: Precision Natural SUSY at CEPC, FCC-ee, and ILC

abstract: Testing the idea of naturalness is and will continue to be one of the most important goals of high energy physics experiments. It will play a central role in the physics program of future colliders. In this paper, we present projections of the reach of natural SUSY at future lepton colliders: CEPC, FCC-ee and ILC. We focus on the observables which give the strongest reach, the electroweak precision observables (for left-handed stops), and Higgs to gluon and photon decay rates (for both left- and right-handed stops). There is a "blind spot" when the stop mixing parameter Xt is approximately equal to the average stop mass. We argue that in natural scenarios, bounds on the heavy Higgs bosons from tree-level mixing effects that modify the bottom Yukawa coupling together with bounds from b to s gamma play a complementary role in probing the blind spot region. For specific natural SUSY scenarios such as folded SUSY in which the top partners do not carry Standard Model color charges, electroweak precision observables could be the most sensitive probe. In all the scenarios discussed in this paper, the combined set of precision measurements will probe down to a few percent in fine-tuning.


title: Soft Theorems from Effective Field Theory

abstract: The singular limits of massless gauge theory amplitudes are described by an effective theory, called soft-collinear effective theory (SCET), which has been applied most successfully to make all-orders predictions for observables in collider physics and weak decays. At tree-level, the emission of a soft gauge boson at subleading order in its energy is given by the Low-Burnett-Kroll theorem, with the angular momentum operator acting on a lower-point amplitude. For well separated particles at tree-level, we prove the Low-Burnett-Kroll theorem using matrix elements of subleading SCET Lagrangian and operator insertions which are individually gauge invariant. These contributions are uniquely determined by gauge invariance and the reparametrization invariance (RPI) symmetry of SCET. RPI in SCET is connected to the infinite-dimensional asymptotic symmetries of the S-matrix. The Low-Burnett-Kroll theorem is generically spoiled by on-shell corrections, including collinear loops and collinear emissions. We demonstrate this explicitly both at tree-level and at one-loop. The effective theory correctly describes these configurations, and we generalize the Low-Burnett-Kroll theorem into a new one-loop subleading soft theorem for amplitudes. Our analysis is presented in a manner that illustrates the wider utility of using effective theory techniques to understand the perturbative S-matrix.


title: Numerical Simulations of Internal Wave Generation by Convection in Water

abstract: Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. In order to isolate the physical mechanism exciting internal waves, we use the novel spectral code Dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations over excite high frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal wave generation by convection.


title: Vast planes of satellites in a high resolution simulation of the Local Group: comparison to Andromeda

abstract: We search for vast planes of satellites (VPoS) in a high resolution simulation of the Local Group performed by the CLUES project, which improves significantly the resolution of former similar studies. We use a simple method for detecting planar configurations of satellites, and validate it on the known plane of M31. We implement a range of prescriptions for modelling the satellite populations, roughly reproducing the variety of recipes used in the literature, and investigate the occurence and properties of planar structures in these populations. The structure of the simulated satellite systems is strongly non-random and contains planes of satellites, predominantly co-rotating, with, in some cases, sizes comparable to the plane observed in M31 by Ibata et al.. However the latter is slightly richer in satellites, slightly thinner and has stronger co-rotation, which makes it stand out as overall more exceptional than the simulated planes, when compared to a random population. Although the simulated planes we find are generally dominated by one real structure, forming its backbone, they are also partly fortuitous and are thus not kinematically coherent structures as a whole. Provided that the simulated and observed planes of satellites are indeed of the same nature, our results suggest that the VPoS of M31 is not a coherent disc and that one third to one half of its satellites must have large proper motions perpendicular to the plane.


title: The variation in molecular gas depletion time among nearby galaxies: II the impact of galaxy internal structures

abstract: We use a data set of nearby galaxies drawn from the HERACLES, ATLAS3D, and COLD GASS surveys to study variations in molecular gas depletion time (Tdep) in galaxy structures such as bulges, grand-design spiral arms, bars and rings. Molecular gas is traced by CO line emission and star formation rate (SFR) is derived using the combination of far-ultraviolet and mid-infrared (MIR) data. The contribution of old stars to MIR emission for the ATLAS3D sample is corrected using 2MASS K-band images. We apply a two-dimensional image decomposition algorithm to decompose galaxies into bulges and discs. Spiral arms, bars and rings are identified in the residual maps, and molecular gas depletion times are derived on a square grid of 1 kpc^2 size. In previous work, we showed that Tdep correlates strongly with specific star formation rate (sSFR). We now find that at a given sSFR, the bulge has shorter Tdep than the disc. The shift to shorter depletion times is most pronounced in the inner bulge (R < 0.1Re). Grids from galaxies with bars and rings are similar to those from galactic bulges in that they have reduced Tdep at a given sSFR. In contrast, the Tdep versus sSFR relation in the discs of galaxies with spiral arms is displaced to longer Tdep at fixed sSFR. We then show that the differences in the Tdep-sSFR relation for bulges, discs, arms, bars and rings can be linked to variations in "stellar", rather than gas surface density between different structures. Our best current predictor for Tdep, both globally and for 1 kpc grids, is given by Tdep= -0.36log(Sigma_SFR)-0.5log(Sigma_*)+5.87.


title: An extreme gravitationally redshifted iron line at 4.8 keV in Mrk 876

abstract: X-ray spectral lines at unforeseen energies are important because they can shed light on the extreme physical conditions of the environment around the supermassive black holes of active galactic nuclei (AGN). Mrk 876 displays such a line at 4.8 keV rest-frame energy. A possible interpretation of its origin can be found in the hotspot scenario. In this scenario the primary radiation from a flare in the hot corona of an AGN illuminates a limited portion of the accretion disk that emits by fluorescence. In this context the line can represent an extreme gravitationally redshifted Fe line originating on the accretion disk below 6 gravitational radii from a rotating supermassive black hole. The correct estimate of the line significance requires a dedicated approach. Based on an existing rigorous approach, we have performed extensive Monte Carlo simulations. We determine that the line is a real feature at ~99 confidence level.


title: Reducing cosmological small scale structure via a large dark matter-neutrino interaction: constraints and consequences

abstract: Cold dark matter explains a wide range of data on cosmological scales. However, there has been a steady accumulation of evidence for discrepancies between simulations and observations at scales smaller than galaxy clusters. Solutions to these small scale structure problems may indicate that simulations need to improve how they include feedback from baryonic matter, or may imply that dark matter properties differ from the standard cold, noninteracting scenario. One promising way to affect structure formation on small scales is a relatively strong coupling of dark matter to neutrinos. We construct an experimentally viable, simple, renormalizable, model with new interactions between neutrinos and dark matter. We show that addressing the small scale structure problems requires dark matter with a mass that is tens of MeV, and a present-day density determined by an initial particle-antiparticle asymmetry in the dark sector. Generating a sufficiently large dark matter-neutrino coupling requires a new heavy neutrino with a mass around 100 MeV. The heavy neutrino is mostly sterile but has a substantial $\tau$ neutrino component, while the three nearly massless neutrinos are partly sterile. We provide the first discussion of how such dark matter-neutrino interactions affect neutrino (especially $\tau$ neutrino) phenomenology. This model can be tested by future astrophysical, particle physics, and neutrino oscillation data. A feature in the neutrino energy spectrum and flavor content from a future nearby supernova would provide strong evidence of neutrino-dark matter interactions. Promising signatures include anomalous matter effects in neutrino oscillations due to nonstandard interactions and a component of the $\tau$ neutrino with mass around 100 MeV.


title: Chemically Distinct Nuclei and Outflowing Shocked Molecular Gas in Arp 220

abstract: We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.


title: Secular dynamics of hierarchical quadruple systems: the case of a triple system orbited by a fourth body

abstract: We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.