title: More on Five Dimensional EVH Black Rings

abstract: In this paper we continue our analysis of arXiv:1308.1478[hep-th] and study in detail the parameter space of three families of doubly spinning black ring solutions: balanced black ring, unbalanced ring and dipole-charged balanced black rings. In all these three families the Extremal Vanishing Horizon (EVH) ring appears in the vanishing limit of the dimensionful parameter of the solution which measures the ring size. We study the near horizon limit of the EVH black rings and for all three cases we find a (pinching orbifold) AdS$_3$ throat with the AdS$_3$ radius $\ell^2=8 G_5 M/(3\pi)$ where $M$ is the ring mass and $G_5$ is the 5d Newton constant. We also discuss the near horizon limit of near-EVH black rings and show that the AdS$_3$ factor is replaced with a generic BTZ black hole. We use these results to extend the EVH/CFT correspondence for black rings, a 2d CFT dual to near-EVH black rings.


title: The Scaling, Splitting and Squaring Method for the Exponential of Perturbed Matrices

abstract: We propose splitting methods for the computation of the exponential of perturbed matrices which can be written as the sum $A=D+\varepsilon B$ of a sparse and efficiently exponentiable matrix $D$ with sparse exponential $e^D$ and a dense matrix $\varepsilon B$ which is of small norm in comparison with $D$. The predominant algorithm is based on scaling the large matrix $A$ by a small number $2^{-s}$, which is then exponentiated by efficient Pad\'e or Taylor methods and finally squared in order to obtain an approximation for the full exponential. In this setting, the main portion of the computational cost arises from dense-matrix multiplications and we present a modified squaring which takes advantage of the smallness of the perturbed matrix $B$ in order to reduce the number of squarings necessary. Theoretical results on local error and error propagation for splitting methods are complemented with numerical experiments and show a clear improvement over existing methods when medium precision is sought.


title: Parallel Worldline Numerics: Implementation and Error Analysis

abstract: We give an overview of the worldline numerics technique, and discuss the parallel CUDA implementation of a worldline numerics algorithm. In the worldline numerics technique, we wish to generate an ensemble of representative closed-loop particle trajectories, and use these to compute an approximate average value for Wilson loops. We show how this can be done with a specific emphasis on cylindrically symmetric magnetic fields. The fine-grained, massive parallelism provided by the GPU architecture results in considerable speedup in computing Wilson loop averages. Furthermore, we give a brief overview of uncertainty analysis in the worldline numerics method. There are uncertainties from discretizing each loop, and from using a statistical ensemble of representative loops. The former can be minimized so that the latter dominates. However, determining the statistical uncertainties is complicated by two subtleties. Firstly, the distributions generated by the worldline ensembles are highly non-Gaussian, and so the standard error in the mean is not a good measure of the statistical uncertainty. Secondly, because the same ensemble of worldlines is used to compute the Wilson loops at different values of $T$ and $x_\mathrm{ cm}$, the uncertainties associated with each computed value of the integrand are strongly correlated. We recommend a form of jackknife analysis which deals with both of these problems.


title: Hopf fibrations for turbulent pipe flows

abstract: We propose a generalization of Hopf fibrations to quotient the streamwise translation symmetry of turbulent pipe flows viewed as dynamical systems. In particular, we exploit the geometric structure of the associate high dimensional state space, which is that of a principal fiber bundle. The relation between the comoving frame velocity $U_{d}$ associated with the dynamical phase of an orbit in the bundle and the Taylor's hypothesis is investigated. As an application, Laser-Induced-Fluorescence techniques are exploited to capture planar fluorescent dye concentration fields tracing a turbulent pipe flow at the bulk Reynolds number $\mathfrak{\mathsf{Re}}=3200$. The symmetry reduction analysis of the experimental data reveals that the speed $u$ of dye concentration bursts is associated with the dynamical and geometric phases of the corresponding orbits in the fiber bundle. In particular, in the symmetry-reduced frame we unveil a pattern-changing dynamics of the passive scalar structures, which explains the observed speed $u\approx U_{d}+U_{g}$ of intense bursting events in terms of the geometric phase velocity $U_{g}\approx0.43U_{d}$ associated with the orbits in the bundle.


title: Statistical Mechanics of Hamiltonian Adaptive Resolution Simulations

abstract: The Adaptive Resolution Scheme (AdResS) is a hybrid scheme that allows one to treat a molecular system with different levels of resolution depending on the location of the molecules. The construction of a Hamiltonian based on the this idea (H-AdResS) allows one to formulate the usual tools of ensembles and statistical mechanics. We present a number of exact and approximate results that provide a statistical mechanics foundation for this simulation method. We also present simulation results that illustrate the theory.


title: Bent crystal spectrometer for both frequency and wavenumber resolved x-ray scattering at a seeded free-electron laser

abstract: We present a cylindrically curved GaAs x-ray spectrometer with energy resolution $\Delta E/E = 1.1\cdot 10^{-4}$ and wave-number resolution of $\Delta k/k = 3\cdot 10^{-3}$, allowing plasmon scattering at the resolution limits of the Linac Coherent Light Source (LCLS) x-ray free-electron laser. It spans scattering wavenumbers of 3.6 to $5.2/$\AA\ in 100 separate bins, with only 0.34\% wavenumber blurring. The dispersion of 0.418~eV/$13.5\,\mu$m agrees with predictions within 1.3\%. The reflection homogeneity over the entire wavenumber range was measured and used to normalize the amplitude of scattering spectra. The proposed spectrometer is superior to a mosaic HAPG spectrometer when the energy resolution needs to be comparable to the LCLS seeded bandwidth of 1~eV and a significant range of wavenumbers must be covered in one exposure.


title: Casimir Interactions between Magnetic Flux Tubes in a Dense Lattice

abstract: We use the worldline numerics technique to study a cylindrically symmetric model of magnetic flux tubes in a dense lattice and the non-local Casimir forces acting between regions of magnetic flux. Within a superconductor the magnetic field is constrained within magnetic flux tubes and if the background magnetic field is on the order the quantum critical field strength, $B_k = \frac{m^2}{e} = 4.4 \times 10^{13}$ Gauss, the magnetic field is likely to vary rapidly on the scales where \acs{QED} effects are important. In this paper, we construct a cylindrically symmetric toy model of a flux tube lattice in which the non-local influence of \acs{QED} on neighbouring flux tubes is taken into account. We compute the effective action densities using the worldline numerics technique. The numerics predict a greater effective energy density in the region of the flux tube, but a smaller energy density in the regions between the flux tubes compared to a locally-constant-field approximation. We also compute the interaction energy between a flux tube and its neighbours as the lattice spacing is reduced from infinity. Because our flux tubes exhibit compact support, this energy is entirely non-local and predicted to be zero in local approximations such as the derivative expansion. This Casimir-Polder energy can take positive or negative values depending on the distance between the flux tubes, and it may cause the flux tubes in neutron stars to form bunches. In addition to the above results we also discuss two important subtleties of determining the statistical uncertainties within the worldline numerics technique and recommend a form of jackknife analysis.


title: Minimum Time Optimal Synthesis for a Control System on SU(2)

abstract: For the time optimal control on an invariant system on SU(2), with two independent controls and a bound on the norm of the control, the extremals of the maximum principle are explicit functions of time and the resulting differential equations can be explicitly integrated. We use this fact here to perform the optimal synthesis for these systems, i.e., find all optimal trajectories. As a consequence, we describe a simple method to find the minimum time control for every desired final condition. Although the Lie group SU(2) is three dimensional, optimal trajectories can be described in the unit disk of the complex plane. We find that a circular trajectory separates optimal trajectories that reach the boundary of the unit disk from the others. Inside this separatrix circle another trajectory (the critical trajectory) plays an important role in that all optimal trajectories end at an intersection with this curve. Our results are of interest to find the minimum time needed to achieve a given evolution of a two level quantum system.


title: Synchrotron and inverse-Compton emissions from pairs formed in GRB afterglows (analytical treatment)

abstract: We calculate the synchrotron and inverse-Compton emissions from pairs formed in GRB afterglows from high-energy photons (above 100 MeV), assuming a power-law photon spectrum C_nu ~ nu^{-2} and considering only the pairs generated from primary high-energy photons. The essential properties of these pairs (number, minimal energy, cooling energy, distribution with energy) and of their emission (peak flux, spectral breaks, spectral slope) are set by the observables GeV fluence Phi (t) = Ft and spectrum, and by the Lorentz factor Gamma and magnetic field B of the source of high-energy photons, at observer-time t. Optical and X-ray pseudo--light-curves F_nu (Gamma) are calculated for given B; proper synchrotron self-Compton light-curves are calculated by setting the dynamics Gamma(t) of the high-energy photons source to be that of a decelerating, relativistic shock. It is found that the emission from pairs can accommodate the flux and decays of the optical flashes measured during the prompt (GRB) phase and of the faster-decaying X-ray plateaus observed during the delayed (afterglow) phase. The brightest pair optical emission is obtained for 100 < Gamma < 500, and depends mostly on the GeV fluence, being independent of the source redshift. Emission from pairs formed during the GRB phase offers an alternate explanation to reverse-shock optical flashes. These two models may be distinguished based on their corresponding flux decay index--spectral slope relations, different correlations with the LAT fluence, or through modeling of the afterglow multiwavelength data.


title: Biphasic, Lyotropic, Active Nematics

abstract: We perform dynamical simulations of a two-dimensional active nematic fluid in coexistence with an isotropic fluid. Drops of active nematic become elongated, and an effective anchoring develops at the nematic-isotropic interface. The activity also causes an undulatory instability of the interface. This results in defects of positive topological charge being ejected into the nematic, leaving the interface with a diffuse negative charge. Quenching the active lyotropic fluid results in a steady state in which phase-separating domains are elongated and then torn apart by active stirring.


title: Search for new particles in events with one lepton and missing transverse momentum in $pp$ collisions at $\sqrt{s}$ = 8 TeV with the ATLAS detector

abstract: This paper presents a search for new particles in events with one lepton (electron or muon) and missing transverse momentum using 20.3 $fb^{-1}$ of proton-proton collision data at $\sqrt{s}=8$ TeV recorded by the ATLAS experiment at the Large Hadron Collider. No significant excess beyond Standard Model expectations is observed. A $W'$ with Sequential Standard Model couplings is excluded at the 95% confidence level for masses up to 3.24 TeV. Excited chiral bosons ($W^{*}$) with equivalent coupling strengths are excluded for masses up to 3.21 TeV. In the framework of an effective field theory limits are also set on the dark matter-nucleon scattering cross-section as well as the mass scale $M_{*}$ of the unknown mediating interaction for dark matter pair production in association with a leptonically decaying $W$.


title: When is a space Menger at infinity?

abstract: We try to characterize those Tychonoff spaces $X$ such that $\beta X \setminus X$ has the Menger property.


title: A correlation between the amount of dark matter in elliptical galaxies and their shape

abstract: We discuss the correlation between the dark matter content of elliptical galaxies and their ellipticities. We then explore a mechanism for which the correlation would emerge naturally. Such mechanism leads to identifying the dark matter particles to gravitons. A similar mechanism is known in Quantum Chromodynamics (QCD) and is essential to our understanding of the mass and structure of baryonic matter.


title: Multiple positive solutions of parabolic systems with nonlinear, nonlocal initial conditions

abstract: In this paper we study the existence, localization and multiplicity of positive solutions for parabolic systems with nonlocal initial conditions. In order to do this, we extend an abstract theory that was recently developed by the authors jointly with Radu Precup, related to the existence of fixed points of nonlinear operators satisfying some upper and lower bounds. Our main tool is the Granas fixed point index theory. We also provide a non-existence result and an example to illustrate our theory.


title: Directed Multicut with linearly ordered terminals

abstract: Motivated by an application in network security, we investigate the following "linear" case of Directed Mutlicut. Let $G$ be a directed graph which includes some distinguished vertices $t_1, \ldots, t_k$. What is the size of the smallest edge cut which eliminates all paths from $t_i$ to $t_j$ for all $i < j$? We show that this problem is fixed-parameter tractable when parametrized in the cutset size $p$ via an algorithm running in $O(4^p p n^4)$ time.


title: CPT, T, and CP Violation in B Mixing

abstract: Recent measurements of observables sensitive to CPT, T and CP violation in neutral B meson mixing by the B-factory experiments Belle and BaBar are presented.


title: Rigidity and bifurcation results for CMC hypersurfaces in warped product spaces

abstract: In this paper, we deduce some rigidity results in warped product spaces under normal variations of CMC hypersurfaces. In particular, we prove the existence of one-parameter families locally rigid on the spatial fiber of Anti-de Sitter Schwarzschild spacetime and one-parameter families with bifurcation points on the spatial fiber of de Sitter Schwarzschild spacetime.


title: Interaction corrections to the minimal conductivity of graphene via dimensional regularization

abstract: We compute the two-loop interaction correction to the minimal conductivity of disorder-free intrinsic graphene with the help of dimensional regularization. The calculation is done in two different ways: via density-density and via current-current correlation functions. Upon properly renormalizing the perturbation theory, in both cases, we find that: $\sigma = \sigma_0\,( 1 + \al\,(19-6\pi)/12) \approx \sigma_0 \,(1 + 0.01\, \al)$, where $\al = e^2 / (4 \pi \hbar v)$ is the renormalized fine structure constant and $\sigma_0 = e^2 / (4 \hbar)$. Our results are consistent with experimental uncertainties and resolve a theoretical dispute.


title: Understanding Random Forests: From Theory to Practice

abstract: Data analysis and machine learning have become an integrative part of the modern scientific methodology, offering automated procedures for the prediction of a phenomenon based on past observations, unraveling underlying patterns in data and providing insights about the problem. Yet, caution should avoid using machine learning as a black-box tool, but rather consider it as a methodology, with a rational thought process that is entirely dependent on the problem under study. In particular, the use of algorithms should ideally require a reasonable understanding of their mechanisms, properties and limitations, in order to better apprehend and interpret their results. Accordingly, the goal of this thesis is to provide an in-depth analysis of random forests, consistently calling into question each and every part of the algorithm, in order to shed new light on its learning capabilities, inner workings and interpretability. The first part of this work studies the induction of decision trees and the construction of ensembles of randomized trees, motivating their design and purpose whenever possible. Our contributions follow with an original complexity analysis of random forests, showing their good computational performance and scalability, along with an in-depth discussion of their implementation details, as contributed within Scikit-Learn. In the second part of this work, we analyse and discuss the interpretability of random forests in the eyes of variable importance measures. The core of our contributions rests in the theoretical characterization of the Mean Decrease of Impurity variable importance measure, from which we prove and derive some of its properties in the case of multiway totally randomized trees and in asymptotic conditions. In consequence of this work, our analysis demonstrates that variable importances [...].


title: $\Lambda_{\bar{\textrm{MS}}}^{(n_f=2)}$ from a momentum space analysis of the quark-antiquark static potential

abstract: We determine $\Lambda_{\bar{\textrm{MS}}}^{(n_f=2)}$ by fitting perturbative expressions for the quark-antiquark static potential to lattice results for QCD with $n_f=2$ dynamical quark flavors. To this end we use the perturbative static potential at the presently best known accuracy, i.e. up to ${\cal O}(\alpha_s^4)$, in momentum space. The lattice potential is computed on a fine lattice with $a \approx 0.042 \, \textrm{fm}$ in position space. To allow for a comparison and matching of both results, the lattice potential is transformed into momentum space by means of a discrete Fourier transform. The value of $\Lambda_{\bar{\textrm{MS}}}^{(n_f=2)}$ is extracted in momentum space. All sources of statistical and systematic errors are discussed. The uncertainty in the value of $\Lambda_{\bar{\textrm{MS}}}^{(n_f=2)}$ is found to be smaller than that obtained in a recent position space analysis of the static potential based on the same lattice data.


title: A Fast Hierarchical Method for Multi-script and Arbitrary Oriented Scene Text Extraction

abstract: Typography and layout lead to the hierarchical organisation of text in words, text lines, paragraphs. This inherent structure is a key property of text in any script and language, which has nonetheless been minimally leveraged by existing text detection methods. This paper addresses the problem of text segmentation in natural scenes from a hierarchical perspective. Contrary to existing methods, we make explicit use of text structure, aiming directly to the detection of region groupings corresponding to text within a hierarchy produced by an agglomerative similarity clustering process over individual regions. We propose an optimal way to construct such an hierarchy introducing a feature space designed to produce text group hypotheses with high recall and a novel stopping rule combining a discriminative classifier and a probabilistic measure of group meaningfulness based in perceptual organization. Results obtained over four standard datasets, covering text in variable orientations and different languages, demonstrate that our algorithm, while being trained in a single mixed dataset, outperforms state of the art methods in unconstrained scenarios.


title: Nonlinear gravitational waves as dark energy in warped spacetimes

abstract: On a warped five-dimensional Friedmann-Lema\^{\i}tre-Robertson-Walker(FLRW) spacetime, dark energy can be induced by a U(1) scalar-gauge field on the brane. We consider a zero effective cosmological constant, i.e., the Randall-Sundrum(RS) fine-tuning and no bulk matter fields. The standard model fields interact via the bulk Weyl tensor and cause brane fluctuations. Due to the warp factor, disturbances don't fade away during the expansion of the universe. The late-time behavior could be significant deviate from the standard evolution of the universe. The effect is triggered by the time-dependent part of the warp factor. The self-gravitating cosmic string builds up a huge mass per unit length in the bulk and can induce massive KK-modes felt on the brane. From a nonlinear perturbation analysis it is found that the effective Einstein equations contain a "back-reaction" term on the righthand side caused by the projected 5D Weyl tensor and can act as a dark energy term. The propagation equations to first order for the metric components and scalar-gauge fields show explicit $\varphi$-dependency.


title: A Projected Preconditioned Conjugate Gradient Algorithm for Computing Many Extreme Eigenpairs of a Hermitian Matrix

abstract: We present an iterative algorithm for computing an invariant subspace associated with the algebraically smallest eigenvalues of a large sparse or structured Hermitian matrix A. We are interested in the case in which the dimension of the invariant subspace is large (e.g., over several hundreds or thousands) even though it may still be small relative to the dimension of A. These problems arise from, for example, density functional theory based electronic structure calculations for complex materials. The key feature of our algorithm is that it performs fewer Rayleigh--Ritz calculations compared to existing algorithms such as the locally optimal precondition conjugate gradient or the Davidson algorithm. It is a block algorithm, hence can take advantage of efficient BLAS3 operations and be implemented with multiple levels of concurrency. We discuss a number of practical issues that must be addressed in order to implement the algorithm efficiently on a high performance computer.


title: SB-Labelings, Distributivity, and Bruhat Order on Sortable Elements

abstract: In this article, we investigate the set of $\gamma$-sortable elements, associated with a Coxeter group $W$ and a Coxeter element $\gamma\in W$, under Bruhat order, and we denote this poset by $\mathcal{B}_{\gamma}$. We show that this poset belongs to the class of SB-lattices recently introduced by Hersh and M\'esz\'aros, by proving a more general statement, namely that all join-distributive lattices are SB-lattices. The observation that $\mathcal{B}_{\gamma}$ is join-distributive is due to Armstrong. Subsequently, we investigate for which finite Coxeter groups $W$ and which Coxeter elements $\gamma\in W$ the lattice $\mathcal{B}_{\gamma}$ is in fact distributive. It turns out that this is the case for the "coincidental" Coxeter groups, namely the groups $A_{n},B_{n},H_{3}$ and $I_{2}(k)$. We conclude this article with a conjectural characteriziation of the Coxeter elements $\gamma$ of said groups for which $\mathcal{B}_{\gamma}$ is distributive in terms of forbidden orientations of the Coxeter diagram.


title: Efficient Regularized Regression for Variable Selection with L0 Penalty

abstract: Variable (feature, gene, model, which we use interchangeably) selections for regression with high-dimensional BIGDATA have found many applications in bioinformatics, computational biology, image processing, and engineering. One appealing approach is the L0 regularized regression which penalizes the number of nonzero features in the model directly. L0 is known as the most essential sparsity measure and has nice theoretical properties, while the popular L1 regularization is only a best convex relaxation of L0. Therefore, it is natural to expect that L0 regularized regression performs better than LASSO. However, it is well-known that L0 optimization is NP-hard and computationally challenging. Instead of solving the L0 problems directly, most publications so far have tried to solve an approximation problem that closely resembles L0 regularization. In this paper, we propose an efficient EM algorithm (L0EM) that directly solves the L0 optimization problem. $L_0$EM is efficient with high dimensional data. It also provides a natural solution to all Lp p in [0,2] problems. The regularized parameter can be either determined through cross-validation or AIC and BIC. Theoretical properties of the L0-regularized estimator are given under mild conditions that permit the number of variables to be much larger than the sample size. We demonstrate our methods through simulation and high-dimensional genomic data. The results indicate that L0 has better performance than LASSO and L0 with AIC or BIC has similar performance as computationally intensive cross-validation. The proposed algorithms are efficient in identifying the non-zero variables with less-bias and selecting biologically important genes and pathways with high dimensional BIGDATA.


title: Matrix nil-clean factorizations over abelian rings

abstract: A ring $R$ is nil-clean if every element in $R$ is the sum of an idempotent and a nilpotent. A ring $R$ is abelian if every idempotent is central. We prove that if $R$ is abelian then $M_n(R)$ is nil-clean if and only if $R/J(R)$ is Boolean and $M_n(J(R))$ is nil. This extend the main results of Breaz et al. ~\cite{BGDT} and that of Ko\c{s}an et al.~\cite{KLZ}.


title: Photon-photon gate via the interaction between two collective Rydberg excitations

abstract: We propose a scheme for a deterministic controlled-phase gate between two photons based on the strong interaction between two stationary collective Rydberg excitations in an atomic ensemble. The distance-dependent character of the interaction causes both a momentum displacement of the collective excitations and unwanted entanglement between them. We show that these effects can be overcome by swapping the collective excitations in space and by optimizing the geometry, resulting in a photon-photon gate with high fidelity and efficiency.


title: Little String Amplitudes (and the Unreasonable Effectiveness of 6D SYM)

abstract: We study tree level scattering amplitudes of four massless states in the double scaled little string theory, and compare them to perturbative loop amplitudes in six-dimensional super-Yang-Mills theory. The little string amplitudes are computed from correlators in the cigar coset CFT and in N=2 minimal models. The results are expressed in terms of integrals of conformal blocks and evaluated numerically in the alpha' expansion. We find striking agreements with up to 2-loop scattering amplitudes of massless gluons in 6D SU(k) SYM at a Z_k invariant point on the Coulomb branch. We comment on the issue of UV divergence at higher loop orders in the gauge theory and discuss the implication of our results.


title: Production of rovibronic ground state RbCs molecules via two-photon cascade decay

abstract: We report the production of ultracold RbCs molecules in the rovibronic ground state, i.e., $X^1\Sigma^+ (v=0, J=0)$, by short-range photoassociation to the $2^3\Pi_{0}$ state followed by spontaneous emission. We use narrowband depletion spectroscopy to probe the distribution of rotational levels formed in the $X ^1\Sigma^+(v=0)$ state. We conclude, based on selection rules, that the primary decay route to $X ^1\Sigma^+(v=0)$ is a two-step cascade decay that leads to as much as $33\%$ branching into the $J=0$ rotational level. The experimental simplicity of our scheme opens up the possibility of easier access to the study and manipulation of ultracold heteronuclear molecules in the rovibronic ground state.


title: Incidence Bounds for Block Designs

abstract: We prove three theorems giving extremal bounds on the incidence structures determined by subsets of the points and blocks of a balanced incomplete block design (BIBD). These results generalize and strengthen known bounds on the number of incidences between points and $m$-flats in affine geometries over finite fields. First, we show an upper bound on the number of incidences between sufficiently large subsets of the points and blocks of a BIBD. Second, we show that a sufficiently large subset of the points of a BIBD determines many t-rich blocks. Third, we show that a sufficiently large subset of the blocks of a BIBD determines many t-rich points. These last two results are new even in the special case of incidences between points and $m$-flats in an affine geometry over a finite field. As a corollary we obtain a tight bound on the number of t-rich points determined by a set of points in a plane over a finite field, and use it to sharpen a result of Iosevich, Rudnev, and Zhai on the number of triangles with distinct areas determined by a set of points in a plane over a finite field.


title: Simulating flexible fiber suspensions using a scalable immersed boundary algorithm

abstract: We present an approach for numerically simulating the dynamics of flexible fibers in a three-dimensional shear flow using a scalable immersed boundary (IB) algorithm based on Guermond and Minev's pseudo-compressible fluid solver. The fibers are treated as one-dimensional Kirchhoff rods that resist stretching, bending, and twisting, within the generalized IB framework. We perform a careful numerical comparison against experiments on single fibers performed by S. G. Mason and co-workers, who categorized the fiber dynamics into several distinct orbit classes. We show that the orbit class may be determined using a single dimensionless parameter for low Reynolds flows. Lastly, we simulate dilute suspensions containing up to hundreds of fibers using a distributed- memory computer cluster. These simulations serve as a stepping stone for studying more complex suspension dynamics including non-dilute suspensions and aggregation of fibers (also known as flocculation).


title: Proximity induced superconductivity in Weyl semi-metals

abstract: We introduce superconducting proximity effects in Weyl semi-metals (WSM) with broken time reversal symmetry by tunnel coupling one of its surfaces to an s-wave superconductor using the Green's function approach. We find that the band structure develops coherence peaks, but despite the presence of metallic states in the bulk, the coherence peaks do not extend far into the bulk and remain confined to a few layers close to the interface, similar to the proximity effect in the topological insulators (TI) which are gapped in the bulk. The Weyl nodes remain unaffected, and in that sense, no true gap develops. We also study the various induced p and s-wave pairing amplitudes classified by their symmetries, as a function of the various relevant parameters of the theory and note the exponential decay of the induced pairings in the bulk both in the TI and the WSM, even at finite chemical potential.


title: Asteroseismic inference on the spin-orbit misalignment and stellar parameters of HAT-P-7

abstract: The measurement of obliquities in star-planet systems is of great importance for the understanding of planet system formation and evolution. The bright and well studied HAT-P-7 system is intriguing as several Rossiter-McLaughlin (RM) measurements found a large projected obliquity in this system, but it was so far not possible to determine if the orbit is polar and/or retrograde. The goal of this study is to measure the stellar inclination and hereby the full 3D obliquity of the HAT-P-7 system instead of only the 2D projection as measured by the RM effect. In addition we provide an updated set of stellar parameters for the star. We use the full set of available observations from Kepler spanning Q0-Q17 to produce the power spectrum of HAT-P-7. We extract oscillation mode frequencies via an MCMC peak-bagging routine, and use the results from this to estimate the stellar inclination angle. Combining this with the projected obliquity from RM and the inclination of the orbital plane allows us to determine the stellar obliquity. We use asteroseismology to model the star from the extracted frequencies using two different approaches to the modelling where either the MESA or the GARSTEC stellar evolution codes are adopted. Using our updated asteroseismic modelling we find, i.a., the following stellar parameters for HAT-P-7: M=1.51{+0.04}{-0.05}Msun, $R=2.00{+0.01}{-0.02}Rsun, and age = 2.07{+0.28}{-0.23} Gyr. Our asteroseismic modelling offers a high precision on the stellar parameters, for instance is the uncertainty on age of the order ~11%. For the stellar inclination we estimate i_star<36.5 deg., which translates to an obliquity between 83 and 111 deg. We find that the planet HAT-P-7b is likely retrograde in its orbit, and that the orbit is close to being polar. The new parameters for the star gives an updated planetary density of 0.65+-0.03 g cm^{-3}, which is lower than previous estimates.


title: Security bound of cheat sensitive quantum bit commitment

abstract: Cheat sensitive quantum bit commitment (CSQBC) loosens the security requirement of quantum bit commitment (QBC), so that the existing impossibility proofs of unconditionally secure QBC can be evaded. But here we analyze the common features in all existing CSQBC protocols, and show that in any CSQBC having these features, the receiver can always learn a non-trivial amount of information on the sender's committed bit before it is unveiled, while his cheating can pass the security check with a probability not less than 50%. The sender's cheating is also studied. The optimal CSQBC protocols that can minimize the sum of the cheating probabilities of both parties are found to be trivial, as they are practically useless. We also discuss the possibility of building a fair protocol in which both parties can cheat with equal probabilities.


title: An estimate of the average number of recessive lethal mutations carried by humans

abstract: The effects of inbreeding on human health depend critically on the number and severity of recessive, deleterious mutations carried by individuals. In humans, existing estimates of these quantities are based on comparisons between consanguineous and non-consanguineous couples, an approach that confounds socioeconomic and genetic effects of inbreeding. To circumvent this limitation, we focused on a founder population with almost complete Mendelian disease ascertainment and a known pedigree. By considering all recessive lethal diseases reported in the pedigree and simulating allele transmissions, we estimated that each haploid set of human autosomes carries on average 0.29 (95% credible interval [0.10, 0.83]) autosomal, recessive alleles that lead to complete sterility or severe disorders at birth or before reproductive age when homozygous. Comparison to existing estimates of the deleterious effects of all recessive alleles suggests that a substantial fraction of the burden of autosomal, recessive variants is due to single mutations that lead to death between birth and reproductive age. In turn, the comparison to estimates from other eukaryotes points to a surprising constancy of the average number of recessive lethal mutations across organisms with markedly different genome sizes.


title: Majorana fermions in ferromagnetic chains on the surface of bulk spin-orbit coupled $s$-wave superconductors

abstract: Majorana fermion (MF) excitations in solid state system have non-Abelian statistics which is essential for topological quantum computation. Previous proposals to realize MF, however, generally requires fine-tuning of parameters. Here we explore a platform which avoids the fine-tuning problem, namely a ferromagnetic chain deposited on the surface of a spin-orbit coupled $s$-wave superconductor. We show that it generically supports zero-energy topological MF excitations near the two ends of the chain with minimal fine-tuning. Depending on the strength of the ferromagnetic moment in the chain, the number of MFs at each end, $n$, can be either one or two, and should be revealed by a robust zero-bias peak (ZBP) of height $2ne^2/h$ in scanning tunneling microscopy (STM) measurements which would show strong (weak) signals at the ends (middle) of the chain. The role of an approximate chiral symmetry which gives an integer topological invariant to the system is discussed.


title: Analysis of Sunyaev-Zel'dovich Effect Mass-Observable Relations using South Pole Telescope Observations of an X-ray Selected Sample of Low Mass Galaxy Clusters and Groups

abstract: (Abridged) We use 95, 150, and 220GHz observations from the SPT to examine the SZE signatures of a sample of 46 X-ray selected groups and clusters drawn from ~6 deg^2 of the XMM-BCS. These systems extend to redshift z=1.02, have characteristic masses ~3x lower than clusters detected directly in the SPT data and probe the SZE signal to the lowest X-ray luminosities (>10^42 erg s^-1) yet. We develop an analysis tool that combines the SZE information for the full ensemble of X-ray-selected clusters. Using X-ray luminosity as a mass proxy, we extract selection-bias corrected constraints on the SZE significance- and Y_500-mass relations. The SZE significance- mass relation is in good agreement with an extrapolation of the relation obtained from high mass clusters. However, the fit to the Y_500-mass relation at low masses, while in good agreement with the extrapolation from high mass SPT clusters, is in tension at 2.8 sigma with the constraints from the Planck sample. We examine the tension with the Planck relation, discussing sample differences and biases that could contribute. We also present an analysis of the radio galaxy point source population in this ensemble of X-ray selected systems. We find 18 of our systems have 843 MHz SUMSS sources within 2 arcmin of the X-ray centre, and three of these are also detected at significance >4 by SPT. Of these three, two are associated with the group brightest cluster galaxies, and the third is likely an unassociated quasar candidate. We examine the impact of these point sources on our SZE scaling relation analyses and find no evidence of biases. We also examine the impact of dusty galaxies using constraints from the 220 GHz data. The stacked sample provides 2.8$\sigma$ significant evidence of dusty galaxy flux, which would correspond to an average underestimate of the SPT Y_500 signal that is (17+-9) per cent in this sample of low mass systems.


title: Congruences for Fishburn numbers modulo prime powers

abstract: The Fishburn numbers $\xi (n)$ are defined by the formal power series \[ \sum_{n \geq 0} \xi (n) q^n = \sum_{n \geq 0} \prod_{j = 1}^n (1 - (1 - q)^j). \] Recently, G. Andrews and J. Sellers discovered congruences of the form $\xi (p m + j) \equiv 0$ modulo $p$, valid for all $m \geq 0$. These congruences have then been complemented and generalized to the case of $r$-Fishburn numbers by F. Garvan. In this note, we answer a question of Andrews and Sellers regarding an extension of these congruences to the case of prime powers. We show that, under a certain condition, all these congruences indeed extend to hold modulo prime powers.


title: Effect of strain on electronic and thermoelectric properties of few layers to bulk MoS$_{2}$

abstract: The sensitive dependence of electronic and thermoelectric properties of MoS$_2$ on the applied strain opens up a variety of applications in the emerging area of straintronics. Using first principles based density functional theory calculations, we show that the band gap of few layers of MoS$_2$ can be tuned by applying i) normal compressive (NC), ii) biaxial compressive (BC), and iii) biaxial tensile (BT) strain. A reversible semiconductor to metal transition (S-M transition) is observed under all three types of strain. In the case of NC strain, the threshold strain at which S-M transition occurs increases with increasing number of layers and becomes maximum for the bulk. On the other hand, the threshold strain for S-M transition in both BC and BT strain decreases with the increase in number of layers. The difference in the mechanisms for the S-M transition is explained for different types of applied strain. Furthermore, the effect of strain type and number of layers on the transport properties are also studied using Botzmann transport theory. We optimize the transport properties as a function of number of layers and applied strain. 3L- and 2L-MoS$_2$ emerge as the most efficient thermoelectric material under NC and BT strain, respectively. The calculated thermopower is large and comparable to some of the best thermoelectric materials. A comparison between the feasibility of these three types of strain is also discussed.


title: Dynamics of swollen fractal networks

abstract: The dynamics of swollen fractal networks (Rouse model) has been studied through computer simulations. The fluctuation-relaxation theorem was used instead of the usual Langevin approach to Brownian dynamics. We measured the equivalent of the mean square displacement $\langle \vec r^{\,2} \rangle$ and the coefficient of self-diffusion $D$ of two-and three-dimensional Sierpinski networks and of the two-dimensional percolation network. The results showed an anomalous diffusion, i. e., a power law for $D$, decreasing with time with an exponent proportional to the spectral dimension of the network.


title: Coherent neutrino-nucleus scattering detection with a CsI[Na] scintillator at the SNS spallation source

abstract: We study the possibility of using CsI[Na] scintillators as an advantageous target for the detection of coherent elastic neutrino-nucleus scattering (CENNS), using the neutrino emissions from the SNS spallation source at Oak Ridge National Laboratory. The response of this material to low-energy nuclear recoils like those expected from this process is characterized. Backgrounds are studied using a 2 kg low-background prototype crystal in a dedicated radiation shield. The conclusion is that a planned 14 kg detector should measure approximately 550 CENNS events per year above a demonstrated $\sim7$ keVnr low-energy threshold, with a signal-to-background ratio sufficient for a first measurement of the CENNS cross-section. The cross-section for the $^{208}$Pb($\nu_{e},e^{-}$)$^{208}$Bi reaction, of interest for future supernova neutrino detection, can be simultaneously obtained.

categories: physics.ins-det astro-ph.IM hep-ex hep-ph nucl-ex


title: No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~2 in CANDELS / 3D-HST

abstract: We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3<z<2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z~2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z~2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z~1.85 - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.


title: Holographic p-wave Superconductor with Disorder

abstract: We implement the effects of disorder on a holographic p-wave superconductor by introducing a random chemical potential which defines the local energy of the charge carriers. Since there are various possibilities for the orientation of the vector order parameter, we explore the behavior of the condensate in the parallel and perpendicular directions to the introduced disorder. We clarify the nature of various branches representing competing solutions and construct the disordered phase diagram. We find that moderate disorder enhances superconductivity as determined by the value of the condensate. Though we mostly focus on uncorrelated noise, we also consider a disorder characterized by its spectral properties and study in detail its influence on the spectral properties of the condensate and charge density. We find fairly universal responses of the resulting power spectra characterized by linear functions of the disorder power spectrum.

categories: hep-th cond-mat.dis-nn cond-mat.str-el cond-mat.supr-con


title: MUSE sneaks a peek at extreme ram-pressure stripping events. I. A kinematic study of the archetypal galaxy ESO137-001

abstract: We present MUSE observations of ESO137-001, a spiral galaxy infalling towards the center of the massive Norma cluster at z~0.0162. During the high-velocity encounter of ESO137-001 with the intracluster medium, a dramatic ram-pressure stripping event gives rise to an extended gaseous tail, traced by our MUSE observations to >30 kpc from the galaxy center. By studying the H-alpha surface brightness and kinematics in tandem with the stellar velocity field, we conclude that ram pressure has completely removed the interstellar medium from the outer disk, while the primary tail is still fed by gas from the inner regions. Gravitational interactions do not appear to be a primary mechanism for gas removal. The stripped gas retains the imprint of the disk rotational velocity to ~20 kpc downstream, without a significant gradient along the tail, which suggests that ESO137-001 is fast moving along a radial orbit in the plane of the sky. Conversely, beyond ~20 kpc, a greater degree of turbulence is seen, with velocity dispersion up to >100 km/s. For a model-dependent infall velocity of ~3000 km/s, we conclude that the transition from laminar to turbulent flow in the tail occurs on timescales >6.5 Myr. Our work demonstrates the terrific potential of MUSE for detailed studies of how ram-pressure stripping operates on small scales, providing a deep understanding of how galaxies interact with the dense plasma of the cluster environment.


title: The Calibration of the WISE W1 and W2 Tully-Fisher Relation

abstract: In order to explore local large-scale structures and velocity fields, accurate galaxy distance measures are needed. We now extend the well-tested recipe for calibrating the correlation between galaxy rotation rates and luminosities -- capable of providing such distance measures -- to the all-sky, space-based imaging data from the Wide-field Infrared Survey Explorer (WISE) W1 ($3.4\mu$m) and W2 ($4.6\mu$m) filters. We find a linewidth to absolute magnitude correlation (known as the Tully-Fisher Relation, TFR) of $\mathcal{M}^{b,i,k,a}_{W1} = -20.35 - 9.56 (\log W^i_{mx} - 2.5)$ (0.54 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.76 - 9.74 (\log W^i_{mx} - 2.5)$ (0.56 magnitudes rms) from 310 galaxies in 13 clusters. We update the I-band TFR using a sample 9% larger than in Tully & Courtois (2012). We derive $\mathcal{M}^{b,i,k}_I = -21.34 - 8.95 (\log W^i_{mx} - 2.5)$ (0.46 magnitudes rms). The WISE TFRs show evidence of curvature. Quadratic fits give $\mathcal{M}^{b,i,k,a}_{W1} = -20.48 - 8.36 (\log W^i_{mx} - 2.5) + 3.60 (\log W^i_{mx} - 2.5)^2$ (0.52 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.91 - 8.40 (\log W^i_{mx} - 2.5) + 4.32 (\log W^i_{mx} - 2.5)^2$ (0.55 magnitudes rms). We apply an I-band -- WISE color correction to lower the scatter and derive $\mathcal{M}_{C_{W1}} = -20.22 - 9.12 (\log W^i_{mx} - 2.5)$ and $\mathcal{M}_{C_{W2}} = -19.63 - 9.11 (\log W^i_{mx} - 2.5)$ (both 0.46 magnitudes rms). Using our three independent TFRs (W1 curved, W2 curved and I-band), we calibrate the UNION2 supernova Type Ia sample distance scale and derive $H_0 = 74.4 \pm 1.4$(stat) $\pm\ 2.4$(sys) kms$^{-1}$ Mpc$^{-1}$ with 4% total error.


title: Revisiting monotop production at the LHC

abstract: Scenarios of new physics where a single top quark can be produced in association with large missing energy (monotop) have been recently studied both from the theoretical point of view and by experimental collaborations. We revisit the originally proposed monotop setup by embedding the effective couplings of the top quark in an SU(2)L invariant formalism. We show that minimality selects one model for each of the possible production mechanisms: a scalar field coupling to a right-handed top quark and an invisible fermion when the monotop system is resonantly produced, and a vector field mediating the interactions of a dark sector to right-handed quarks for the non-resonant production mode. We study in detail constraints on the second class of scenarios, originating from contributions to standard single top processes when the mediator is lighter than the top quark and from the dark matter relic abundance when the mediator is heavier than the top quark.


title: Finding $\eta$ Car Analogs in Nearby Galaxies Using Spitzer: II. Identification of An Emerging Class of Extragalactic Self-Obscured Stars

abstract: Understanding the late-stage evolution of the most massive stars such as $\eta$ Carinae is challenging because no true analogs of $\eta$ Car have been clearly identified in the Milky Way or other galaxies. In Khan et. al. (2013), we utilized Spitzer IRAC images of $7$ nearby ($\lesssim4$ Mpc) galaxies to search for such analogs, and found $34$ candidates with flat or red mid-IR spectral energy distributions. Here, in Paper II, we present our characterization of these candidates using multi-wavelength data from the optical through the far-IR. Our search detected no true analogs of $\eta$ Car, which implies an eruption rate that is a fraction $0.01\lesssim F \lesssim 0.19$ of the ccSN rate. This is roughly consistent with each $M_{ZAMS} \gtrsim 70M_\odot$ star undergoing $1$ or $2$ outbursts in its lifetime. However, we do identify a significant population of $18$ lower luminosity $\left(\log(L/L_\odot)\simeq5.5-6.0\right)$ dusty stars. Stars enter this phase at a rate that is fraction $0.09 \lesssim F \lesssim 0.55$ of the ccSN rate, and this is consistent with all $25 < M_{ZAMS} < 60M_\odot$ stars undergoing an obscured phase at most lasting a few thousand years once or twice. These phases constitute a negligible fraction of post-main sequence lifetimes of massive stars, which implies that these events are likely to be associated with special periods in the evolution of the stars. The mass of the obscuring material is of order $\sim M_\odot$, and we simply do not find enough heavily obscured stars for theses phases to represent more than a modest fraction ($\sim 10\%$ not $\sim 50\%$) of the total mass lost by these stars. In the long term, the sources that we identified will be prime candidates for detailed physical analysis with JWST.


title: Chaotic cold accretion on to black holes in rotating atmospheres

abstract: Chaotic cold accretion (CCA) profoundly differs from classic black hole accretion models. Using 3D high-resolution simulations, we probe the impact of rotation on the hot and cold accretion flow in a typical massive galaxy. In the hot mode, with or without turbulence, the pressure-dominated flow forms a geometrically thick rotational barrier, suppressing the accretion rate to ~1/3 of the Bondi rate. When radiative cooling is dominant, the gas loses pressure support and quickly circularizes in a cold thin disk. In the more common state of a turbulent and heated atmosphere, CCA drives the dynamics if the gas velocity dispersion exceeds the rotational velocity, i.e., turbulent Taylor number < 1. Extended multiphase filaments condense out of the hot phase via thermal instability and rain toward the black hole, boosting the accretion rate up to 100 times the Bondi rate. Initially, turbulence broadens the angular momentum distribution of the hot gas, allowing the cold phase to condense with prograde or retrograde motion. Subsequent chaotic collisions between the cold filaments, clouds, and a clumpy variable torus promote the cancellation of angular momentum, leading to high accretion rates. The simulated sub-Eddington accretion rates cover the range inferred from AGN cavity observations. CCA predicts inner flat X-ray temperature and $r^{-1}$ density profiles, as recently discovered in M 87 and NGC 3115. The synthetic H{\alpha} images reproduce the main features of cold gas observations in massive ellipticals, as the line fluxes and the filaments versus disk morphology. Such dichotomy is key for the long-term AGN feedback cycle. As gas cools, filamentary CCA develops and boosts AGN heating; the cold mode is thus reduced and the rotating disk remains the sole cold structure. Its consumption leaves the atmosphere in hot mode with suppressed accretion and feedback, reloading the cycle.


title: A Variable PV Broad Absorption Line and Quasar Outflow Energetics

abstract: Broad absorption lines (BALs) in quasar spectra identify high velocity outflows that might exist in all quasars and could play a major role in feedback to galaxy evolution. The viability of BAL outflows as a feedback mechanism depends on their kinetic energies, as derived from the outflow velocities, column densities, and distances from the central quasar. We estimate these quantities for the quasar, Q1413+1143 (redshift $z_e = 2.56$), aided by the first detection of PV $\lambda\lambda$1118,1128 BAL variability in a quasar. In particular, PV absorption at velocities where the CIV trough does not reach zero intensity implies that the CIV BAL is saturated and the absorber only partially covers the background continuum source (with characteristic size <0.01 pc). With the assumption of solar abundances, we estimate that the total column density in the BAL outflow is log N_H > 22.3 (cm^-2). Variability in the PV and saturated CIV BALs strongly disfavors changes in the ionization as the cause of the BAL variability, but supports models with high-column density BAL clouds moving across our lines of sight. The observed variability time of 1.6 yr in the quasar rest frame indicates crossing speeds >750 km/s and a radial distance from the central black hole of <3.5 pc, if the crossing speeds are Keplerian. The total outflow mass is ~4100 M_solar, the kinetic energy ~4x10^54 erg, and the ratio of the outflow kinetic energy luminosity to the quasar bolometric luminosity is ~0.02 (at the minimum column density and maximum distance), which might be sufficient for important feedback to the quasar's host galaxy.


title: Fractionalized gapless quantum vortex liquids

abstract: The standard theoretical approach to gapless spin liquid phases of two-dimensional frustrated quantum antiferromagnets invokes the concept of fermionic slave particles into which the spin fractionalizes. As an alternate we explore new kinds of gapless spin liquid phases in frustrated quantum magnets with XY anisotropy where the vortex of the spin fractionalizes into gapless itinerant fermions. The resulting gapless fractionalized vortex liquid phases are studied within a slave particle framework that is dual to the usual one. We demonstrate the stability of some such phases and describe their properties. We give an explicit construction in an XY-spin-1 system on triangular lattice, and interpret it as a critical phase in the vicinity of spin-nematic states.


title: Is the CNMSSM more credible than the CMSSM?

abstract: With Bayesian statistics, we investigate the full parameter space of the constrained "next-to-minimal" supersymmetric Standard Model (CNMSSM) with naturalness priors, which were derived in a previous work. In the past, most Bayesian analyses of the CNMSSM ignored naturalness of the electroweak (EW) scale by making prejudicial assumptions for parameters defined at the EW scale. We test the CNMSSM against the CMSSM with Bayesian evidence, which, with naturalness priors, incorporates a penalty for fine-tuning of the EW scale. With the evidence, we measure credibility with respect to all measurements, including the EW scale and LHC direct searches. We find that the evidence in favor of the CNMSSM versus the CMSSM is "positive" to "strong" but that if one ignores the $\mu$-problem, the evidence is "barely worth mentioning" to "positive." The $\mu$-problem significantly influences our findings. Unless one considers the $\mu$-problem, the evidence in favor of the CNMSSM versus the CMSSM is at best "positive," which is two grades below "very strong." We, furthermore, identify the most probable regions of the CMSSM and CNMSSM parameter spaces and examine prospects for future discovery at hadron colliders.


title: Absence of diffusion in an interacting system of spinless fermions on a one-dimensional disordered lattice

abstract: We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems.


title: Describing the Observed Cosmic Neutrinos by Interactions of Nuclei with Matter

abstract: IceCube have observed neutrinos which are presumably of extra-galactic origin. Since specific sources have not yet been identified, we discuss what could be learned from the conceptual point of view. We use a simple model for neutrino production from the interactions between nuclei and matter, and we focus on the description of the spectral shape and flavor composition observed by IceCube. Our main parameters are spectral index, maximal energy, magnetic field, and composition of the accelerated nuclei. We show that a cutoff at PeV energies can be achieved by soft enough spectra, a cutoff of the primary energy, or strong enough magnetic fields. These options, however, are difficult to reconcile with the hypothesis that these neutrinos originate from the same sources as the ultra-high energy cosmic rays. We demonstrate that heavier nuclei accelerated in the sources may be a possible way out if the maximal energy scales appropriately with the mass number of the nuclei. In this scenario, neutrino observations can actually be used to test the UHECR acceleration mechanism. We also emphasize the need for a volume upgrade of the IceCube detector for future precision physics, for which the flavor information becomes a statistical meaningful model discriminator as qualitatively new ingredient.


title: Mass production of volume phase holographic gratings for the VIRUS spectrograph array

abstract: The Visible Integral-field Replicable Unit Spectrograph (VIRUS) is a baseline array of 150 copies of a simple, fiber-fed integral field spectrograph that will be deployed on the Hobby-Eberly Telescope (HET). VIRUS is the first optical astronomical instrument to be replicated on an industrial scale, and represents a relatively inexpensive solution for carrying out large-area spectroscopic surveys, such as the HET Dark Energy Experiment (HETDEX). Each spectrograph contains a volume phase holographic (VPH) grating with a 138 mm diameter clear aperture as its dispersing element. The instrument utilizes the grating in first-order for 350-550 nm. Including witness samples, a suite of 170 VPH gratings has been mass produced for VIRUS. Here, we present the design of the VIRUS VPH gratings and a discussion of their mass production. We additionally present the design and functionality of a custom apparatus that has been used to rapidly test the first-order diffraction efficiency of the gratings for various discrete wavelengths within the VIRUS spectral range. This device has been used to perform both in-situ tests to monitor the effects of adjustments to the production prescription as well as to carry out the final acceptance tests of the gratings' diffraction efficiency. Finally, we present the as-built performance results for the entire suite of VPH gratings.


title: HIghMass -- High HI Mass, HI-rich Galaxies at z~0: Sample Definition, Optical and Halpha Imaging, and Star Formation Properties

abstract: We present first results of the study of a set of exceptional HI sources identified in the 40% ALFALFA extragalactic HI survey catalog alpha.40 as being both HI massive (M_HI > 10^10 Msun) and having high gas fractions for their stellar masses: the HIghMass galaxy sample. We analyze UV- and optical-broadband and Halpha images to understand the nature of their relatively underluminous disks in optical and to test whether their high gas fractions can be tracked to higher dark matter halo spin parameters or late gas accretion. Estimates of their star formation rates (SFRs) based on SED-fitting agree within uncertainties with the Halpha luminosity inferred SFRs. The HII region luminosity functions have standard slopes at the luminous end. The global SFRs demonstrate that the HIghMass galaxies exhibit active ongoing star formation (SF) with moderate SF efficiency, but relative to normal spirals, a lower integrated SFR in the past. Because the SF activity in these systems is spread throughout their extended disks, they have overall lower SFR surface densities and lower surface brightness in the optical bands. Relative to normal disk galaxies, the majority of HIghMass galaxies have higher Halpha equivalent widths and are bluer in their outer disks, implying an inside-out disk growth scenario. Downbending double exponential disks are more frequent than upbending disks among the gas-rich galaxies, suggesting that SF thresholds exist in the downbending disks, probably as a result of concentrated gas distribution.


title: The origin of the 'local' 1/4 keV X-ray flux in both charge exchange and a hot bubble

abstract: The Solar neighborhood is the closest and most easily studied sample of the Galactic interstellar medium, an understanding of which is essential for models of star formation and galaxy evolution. Observations of an unexpectedly intense diffuse flux of easily-absorbed 1/4 keV X rays, coupled with the discovery that interstellar space within ~100 pc of the Sun is almost completely devoid of cool absorbing gas led to a picture of a "local cavity" filled with X-ray emitting hot gas dubbed the local hot bubble. This model was recently upset by suggestions that the emission could instead be produced readily within the solar system by heavy solar wind ions charge exchanging with neutral H and He in interplanetary space, potentially removing the major piece of evidence for the existence of million-degree gas within the Galactic disk. Here we report results showing that the total solar wind charge exchange contribution is 40% +/- 5% (stat) +/- 5% (sys) of the 1/4 keV flux in the Galactic plane. The fact that the measured flux is not dominated by charge exchange supports the notion of a million-degree hot bubble of order 100 pc extent surrounding the Sun.


title: The Euler current and relativistic parity odd transport

abstract: For a spacetime of odd dimensions endowed with a unit vector field, we introduce a new topological current that is identically conserved and whose charge is equal to the Euler character of the even dimensional spacelike foliations. The existence of this current allows us to introduce new Chern-Simons-type terms in the effective field theories describing relativistic quantum Hall states and (2+1) dimensional superfluids. Using effective field theory, we calculate various correlation functions and identify transport coefficients. In the quantum Hall case, this current provides the natural relativistic generalization of the Wen-Zee term, required to characterize the shift and Hall viscosity in quantum Hall systems. For the superfluid case this term is required to have nonzero Hall viscosity and to describe superfluids with non s-wave pairing.


title: An apparent metal insulator transition in high mobility 2D InAs heterostructures

abstract: We report on the first experimental observation of an apparent metal insulator transition in a 2D electron gas confined in an InAs quantum well. At high densities we find that the carrier mobility is limited by background charged impurities and the temperature dependence of the resistivity shows a metallic behavior with resistivity increasing with increasing temperature. At low densities we find an insulating behavior below a critical density of $n_{c} = 5 \times 10^{10}$ cm$^{-2}$ with the resistivity decreasing with increasing temperature. We analyze this transition using a percolation model arising from the failure of screening in random background charged impurities. We also examine the percolation transition experimentally by introducing remote ionized impurities at the surface. Using a bias during cool-down, we modify the screening charge at the surface which strongly affects the critical density. Our study shows that transition from a metallic to an insulating phase in our system is due to percolation transition.


title: Generalised geometry for string corrections

abstract: We present a general formalism for incorporating the string corrections in generalised geometry, which necessitates the extension of the generalised tangent bundle. Not only are such extensions obstructed, string symmetries and the existence of a well-defined effective action require a precise choice of the (generalised) connection. The action takes a universal form given by a generalised Lichnerowitz--Bismut theorem. As examples of this construction we discuss the corrections linear in $\alpha'$ in heterotic strings and the absence of such corrections for type II theories.


title: Deep Chandra Observations of HCG 16 - I. Active Nuclei, Star formation and Galactic Winds

abstract: We present new, deep Chandra X-ray and Giant Metrewave Radio Telescope 610~MHz observations of the spiral-galaxy-rich compact group HCG 16, which we use to examine nuclear activity, star formation and the high luminosity X-ray binary populations in the major galaxies. We confirm the presence of obscured active nuclei in NGC 833 and NGC 835, and identify a previously unrecognized nuclear source in NGC 838. All three nuclei are variable on timescales of months to years, and for NGC 833 and NGC 835 this is most likely caused by changes in accretion rate. The deep Chandra observations allow us to detect for the first time an Fe-K$\alpha$ emission line in the spectrum of the Seyfert 2 nucleus of NGC 835. We find that NGC 838 and NGC 839 are both starburst-dominated systems, with only weak nuclear activity, in agreement with previous optical studies. We estimate the star formation rates in the two galaxies from their X-ray and radio emission, and compare these results with estimates from the infra-red and ultra-violet bands to confirm that star formation in both galaxies is probably declining after galaxy-wide starbursts were triggered ~400-500 Myr ago. We examine the physical properties of their galactic superwinds, and find that both have temperatures of ~0.8 keV. We also examine the X-ray and radio properties of NGC 848, the fifth largest galaxy in the group, and show that it is dominated by emission from its starburst.


title: Re-Examining Astrophysical Constraints on the Dark Matter Model

abstract: Recent high-resolution simulations that include Cold Dark Matter (CDM) and baryons have shown that baryonic physics can dramatically alter the dark matter structure of galaxies. These results modify our predictions for observed galaxy evolution and structure. Given these updated expectations, it is timely to re-examine observational constraints on the dark matter model. A few observations are reviewed that may indirectly trace dark matter, and may help confirm or deny possible dark matter models. Warm Dark Matter (WDM) and Self-Interacting Dark Matter (SIDM) are currently the favorite alternative models to CDM. Constraints on the WDM particle mass require it to be so heavy that WDM is nearly indistinguishable from CDM. The best observational test of SIDM is likely to be in the dark matter distribution of faint dwarf galaxies, but there is a lack of theoretical predictions for galaxy structure in SIDM that account for the role of baryons.


title: Cosmological Impact of Population III Binaries

abstract: We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $\sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.

categories: astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.SR


title: Deep Chandra Observations of HCG 16 - II. The Development of the Intra-group Medium in a Spiral-Rich Group

abstract: We use a combination of deep Chandra X-ray observations and radio continuum imaging to investigate the origin and current state of the intra-group medium in the spiral-rich compact group HCG 16. We confirm the presence of a faint ($L_{X,{\rm bolo}}$=1.87$^{+1.03}_{-0.66}$$\times$10$^{41}$ erg/s), low temperature (0.30$^{+0.07}_{-0.05}$ keV) intra-group medium (IGM) extending throughout the ACIS-S3 field of view, with a ridge linking the four original group members and extending to the southeast, as suggested by previous Rosat and XMM-Newton observations. This ridge contains 6.6$^{+3.9}_{-3.3}$$\times$10$^9$ solar masses of hot gas and is at least partly coincident with a large-scale HI tidal filament, indicating that the IGM in the inner part of the group is highly multi-phase. We present evidence that the group is not yet virialised, and show that gas has probably been transported from the starburst winds of NGC 838 and NGC 839 into the surrounding IGM. Considering the possible origin of the IGM, we argue that material ejected by galactic winds may have played a significant role, contributing 20-40% of the observed hot gas in the system.


title: Probing the Terrestrial Regions of Planetary Systems: Warm Debris Disks with Emission Features

abstract: Observations of debris disks allow for the study of planetary systems, even where planets have not been detected. However, debris disks are often only characterized by unresolved infrared excesses that resemble featureless blackbodies, and the location of the emitting dust is uncertain due to a degeneracy with the dust grain properties. Here we characterize the Spitzer IRS spectra of 22 debris disks exhibiting 10 micron silicate emission features. Such features arise from small warm dust grains, and their presence can significantly constrain the orbital location of the emitting debris. We find that these features can be explained by the presence of an additional dust component in the terrestrial zones of the planetary systems, i.e. an exozodiacal belt. Aside from possessing exozodiacal dust, these debris disks are not particularly unique; their minimum grain sizes are consistent with the blowout sizes of their systems, and their brightnesses are comparable to those of featureless warm debris disks. These disks are in systems with a range of ages, although the older systems with features are found only around A-type stars. The features in young systems may be signatures of terrestrial planet formation. Analyzing the spectra of unresolved debris disks with emission features may be one of the simplest and most accessible ways to study the terrestrial regions of planetary systems.


title: Raytracing simulations of coupled dark energy models

abstract: Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $\Lambda$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $\Lambda$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $\Lambda$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $\Lambda$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.


title: Bayesian parameter estimation of core collapse supernovae using gravitational wave simulations

abstract: Using the latest numerical simulations of rotating stellar core collapse, we present a Bayesian framework to extract the physical information encoded in noisy gravitational wave signals. We fit Bayesian principal component regression models with known and unknown signal arrival times to reconstruct gravitational wave signals, and subsequently fit known astrophysical parameters on the posterior means of the principal component coefficients using a linear model. We predict the ratio of rotational kinetic energy to gravitational energy of the inner core at bounce by sampling from the posterior predictive distribution, and find that these predictions are generally very close to the true parameter values, with $90\%$ credible intervals $\sim 0.04$ and $\sim 0.06$ wide for the known and unknown arrival time models respectively. Two supervised machine learning methods are implemented to classify precollapse differential rotation, and we find that these methods discriminate rapidly rotating progenitors particularly well. We also introduce a constrained optimization approach to model selection to find an optimal number of principal components in the signal reconstruction step. Using this approach, we select 14 principal components as the most parsimonious model.


title: Pair-Instability Supernovae of Non-Zero Metallicity Stars

abstract: Observational evidence suggests that some very massive stars in the local Universe may die as pair-instability supernovae. We present 2D simulations of the pair-instability supernova of a non-zero metallicity star. We find that very little mixing occurs in this explosion because metals in the stellar envelope drive strong winds that strip the hydrogen envelope from the star prior to death. Consequently, a reverse shock cannot form and trigger fluid instabilities during the supernova. Only weak mixing driven by nuclear burning occurs in the earliest stages of the supernova, and it is too weak to affect the observational signatures of the explosion.


title: Free function theory through matrix invariants

abstract: In this article we introduce powerful tools and techniques from invariant theory to free analysis. This enables us to study free maps with involution. These maps are free noncommutative analogs of real analytic functions of several variables. With examples we demonstrate that they do not exhibit strong rigidity properties of their involution-free free counterparts. We present a characterization of polynomial free maps via properties of their finite-dimensional slices. This is used to establish power series expansions for analytic free maps about scalar and non-scalar points; the latter are given by series of generalized polynomials for which we obtain an invariant-theoretic characterization. Finally, we give an inverse and implicit function theorem for free maps with involution.


title: Reengineering observatory operations for the time domain

abstract: Observatories are complex scientific and technical institutions serving diverse users and purposes. Their telescopes, instruments, software, and human resources engage in interwoven workflows over a broad range of timescales. These workflows have been tuned to be responsive to concepts of observatory operations that were applicable when various assets were commissioned, years or decades in the past. The astronomical community is entering an era of rapid change increasingly characterized by large time domain surveys, robotic telescopes and automated infrastructures, and - most significantly - of operating modes and scientific consortia that span our individual facilities, joining them into complex network entities. Observatories must adapt and numerous initiatives are in progress that focus on redesigning individual components out of the astronomical toolkit. New instrumentation is both more capable and more complex than ever, and even simple instruments may have powerful observation scripting capabilities. Remote and queue observing modes are now widespread. Data archives are becoming ubiquitous. Virtual observatory standards and protocols and astroinformatics data-mining techniques layered on these are areas of active development. Indeed, new large-aperture ground-based telescopes may be as expensive as space missions and have similarly formal project management processes and large data management requirements. This piecewise approach is not enough. Whatever challenges of funding or politics facing the national and international astronomical communities it will be more efficient - scientifically as well as in the usual figures of merit of cost, schedule, performance, and risks - to explicitly address the systems engineering of the astronomical community as a whole.


title: Experimental excitation of multiple surface-plasmon-polariton waves with 2D gratings

abstract: The excitation of multiple SPP waves as Floquet harmonics was demonstrated in structures fabricated as one-dimensional photonic crystals (PCs) on top of two-dimensional gold gratings. Each period of the PC comprised nine layers of silicon oxynitrides of different compositions, and each PC had either two or three periods. Absorptances for obliquely incident $p$- and $s$-polarized light were measured in the 500--1000-nm wavelength regime and the sharp bands in the absorptance spectra were compared with the solutions of the underlying canonical boundary-value problem. The excitation of multiple surface-plasmon-polariton (SPP) waves as Floquet harmonics was confirmed. The structures demonstrated broadband absorption with overall weak dependences on the incidence angle and the polarization state of the incident light, and has potential application for harvesting solar energy.


title: Substrate Sensitive Mid-Infrared Photoresponse in Graphene

abstract: We report mid-infrared photocurrent spectra of graphene nanoribbon arrays on SiO2 dielectrics showing dual signatures of the substrate interaction. First, hybrid polaritonic modes of graphene plasmons and dielectric surface polar phonons produce a thermal photocurrent in graphene with spectral features that are tunable by gate voltage, nanoribbon width, and light polarization. Secondly, phonon-polaritons associated with the substrate are excited, which indirectly heat up the graphene leading to a graphene photocurrent with fixed spectral features. Models for other commonly used substrates show that the responsivity of graphene infrared photodetectors can be tailored to specific mid-IR frequency bands by the choice of the substrate.


title: Galaxy And Mass Assembly (GAMA): Stellar mass functions by Hubble type

abstract: We present an estimate of the galaxy stellar mass function and its division by morphological type in the local (0.025 < z < 0.06) Universe. Adopting robust morphological classifications as previously presented (Kelvin et al.) for a sample of 3,727 galaxies taken from the Galaxy And Mass Assembly survey, we define a local volume and stellar mass limited sub-sample of 2,711 galaxies to a lower stellar mass limit of M = 10^9.0 M_sun. We confirm that the galaxy stellar mass function is well described by a double Schechter function given by M* = 10^10.64 M_sun, {\alpha}1 = -0.43, {\phi}*1 = 4.18 dex^-1 Mpc^-3, {\alpha}2 = -1.50 and {\phi}*2 = 0.74 dex^-1 Mpc^-3. The constituent morphological-type stellar mass functions are well sampled above our lower stellar mass limit, excepting the faint little blue spheroid population of galaxies. We find approximately 71+3-4% of the stellar mass in the local Universe is found within spheroid dominated galaxies; ellipticals and S0-Sas. The remaining 29+4-3% falls predominantly within late type disk dominated systems, Sab-Scds and Sd-Irrs. Adopting reasonable bulge-to-total ratios implies that approximately half the stellar mass today resides in spheroidal structures, and half in disk structures. Within this local sample, we find approximate stellar mass proportions for E : S0-Sa : Sab-Scd : Sd-Irr of 34 : 37 : 24 : 5.


title: Entropic one-class classifiers

abstract: The one-class classification problem is a well-known research endeavor in pattern recognition. The problem is also known under different names, such as outlier and novelty/anomaly detection. The core of the problem consists in modeling and recognizing patterns belonging only to a so-called target class. All other patterns are termed non-target, and therefore they should be recognized as such. In this paper, we propose a novel one-class classification system that is based on an interplay of different techniques. Primarily, we follow a dissimilarity representation based approach; we embed the input data into the dissimilarity space by means of an appropriate parametric dissimilarity measure. This step allows us to process virtually any type of data. The dissimilarity vectors are then represented through a weighted Euclidean graphs, which we use to (i) determine the entropy of the data distribution in the dissimilarity space, and at the same time (ii) derive effective decision regions that are modeled as clusters of vertices. Since the dissimilarity measure for the input data is parametric, we optimize its parameters by means of a global optimization scheme, which considers both mesoscopic and structural characteristics of the data represented through the graphs. The proposed one-class classifier is designed to provide both hard (Boolean) and soft decisions about the recognition of test patterns, allowing an accurate description of the classification process. We evaluate the performance of the system on different benchmarking datasets, containing either feature-based or structured patterns. Experimental results demonstrate the effectiveness of the proposed technique.


title: Collective effects in light-heavy ion collisions

abstract: We present results for the azimuthal anisotropy of charged hadron distributions in A+A, p+A, d+A, and He3+A collisions within the IP-Glasma+MUSIC model. Obtained anisotropies are due to the fluid dynamic response of the system to the fluctuating initial geometry of the interaction region. While the elliptic and triangular anisotropies in peripheral Pb+Pb collisions at root-s=2.76 TeV are well described by the model, the same quantities in root-s=5.02 TeV p+Pb collisions underestimate the experimental data. This disagreement can be due to neglected initial state correlations or the lack of a detailed description of the fluctuating spatial structure of the proton, or both. We further present predictions for azimuthal anisotropies in p+Au, d+Au, and He3+Au collisions at root-s=200 GeV. For d+Au and 3He+Au collisions we expect the detailed substructure of the nucleon to become less important.


title: Unexpected Universality in the Viscosity of Metallic Liquids

abstract: The range of the magnitude of the liquid viscosity as a function of the temperature (T) is one of the most impressive of any physical property, changing by approximately 17 orders of magnitude from its extrapolated value at infinite temperature to that at the glass transition. We present experimental measurements of containerlessly processed metallic liquids that reveal that the ratio of the viscosity to its extrapolated infinite temperature value follows a universal function of Tcoop/T. The temperature Tcoop corresponds to the onset of cooperative motion and is strongly correlated with the glass transition temperature. On average the extrapolated infinite temperature viscosity is found to be nh, where h is Planck's constant and n is the particle number density. A surprising universality in the viscosity of metallic liquids and its relation to the glass transition is demonstrated.


title: Toward a multilevel representation of protein molecules: comparative approaches to the aggregation/folding propensity problem

abstract: This paper builds upon the fundamental work of Niwa et al. [34], which provides the unique possibility to analyze the relative aggregation/folding propensity of the elements of the entire Escherichia coli (E. coli) proteome in a cell-free standardized microenvironment. The hardness of the problem comes from the superposition between the driving forces of intra- and inter-molecule interactions and it is mirrored by the evidences of shift from folding to aggregation phenotypes by single-point mutations [10]. Here we apply several state-of-the-art classification methods coming from the field of structural pattern recognition, with the aim to compare different representations of the same proteins gathered from the Niwa et al. data base; such representations include sequences and labeled (contact) graphs enriched with chemico-physical attributes. By this comparison, we are able to identify also some interesting general properties of proteins. Notably, (i) we suggest a threshold around 250 residues discriminating "easily foldable" from "hardly foldable" molecules consistent with other independent experiments, and (ii) we highlight the relevance of contact graph spectra for folding behavior discrimination and characterization of the E. coli solubility data. The soundness of the experimental results presented in this paper is proved by the statistically relevant relationships discovered among the chemico-physical description of proteins and the developed cost matrix of substitution used in the various discrimination systems.


title: Towards Automatic Migration of ROS Components from Software to Hardware

abstract: The use of the ROS middleware is a growing trend in robotics in general, ROS and hard real-time embedded systems have however not been easily uniteable while retaining the same overall communication and processing methodology at all levels. In this paper we present an approach aimed at tackling the schism between high-level, flexible software and low-level, real-time software. The key idea of our approach is to enable software components written for a high-level publish-subscribe software architecture to be automatically migrated to a dedicated hardware architecture implemented using programmable logic. Our approach is based on the Unity framework, a unified software/hardware framework based on FPGAs for quickly interfacing high-level software to low-level robotics hardware.


title: Supernovae at the cosmic dawn

abstract: Modern cosmological simulations predict that the first generation of stars formed with a mass scale around 100 solar masses about 300-400 million years after the Big Bang. When the first stars reached the end of their lives, many of them might have died as energetic supernovae that could have significantly affected the early Universe via injecting large amounts of energy and metals into the primordial intergalactic medium. In this paper, we review the current models of the first supernovae by discussing on the relevant background physics, computational methods, and the latest results.

categories: astro-ph.HE astro-ph.CO astro-ph.SR physics.comp-ph


title: Volume Integral of the Pressure Gradient and Archimedes' Principle

abstract: The theorems of vector analysis (divergence theorem, etc.) are typically first applied in the undergraduate physics curriculum in the context of the electromagnetic field and the differential forms of Maxwell's equations. However, these tools are analyzed in depth several courses later in the junior-senior level. I discuss here a "bridge" problem, using the language of vector calculus in a mechanics setting to understand Archimedes' principle as a consequence of hydrostatic equilibrium and the superposition of the external forces. It is my hope that this treatment will help students better integrate and understand understand these and similar vector analysis results in contexts beyond electromagnetism.


title: The Time Domain Lippmann-Schwinger Equation and Convolution Quadrature

abstract: We consider time domain acoustic scattering from a penetrable medium with a variable sound speed. This problem can be reduced to solving a time domain volume Lippmann-Schwinger integral equation. Using convolution quadrature in time and trigonometric collocation in space we can compute an approximate solution. We prove that the time domain Lippmann-Schwinger equation has a unique solution and prove conditional convergence and error estimates for the fully discrete solution for smooth sound speeds. Preliminary numerical results show that the method behaves well even for discontinuous sound speeds.


title: Continuity of the Perron Root

abstract: That the Perron root of a square nonnegative matrix A varies continuously with the entries in A is a corollary of theorems regarding continuity of eigenvalues or roots of polynomial equations, the proofs of which necessarily involve complex numbers. But since continuity of the Perron root is a question that is entirely in the field of real numbers, it seems reasonable that there should exist a development involving only real analysis. This article presents a simple and completely self-contained development that depends only on real numbers and first principles.


title: A restriction on proper actions on homogeneous spaces of reductive type

abstract: Let L be a reductive subgroup of a reductive Lie group G. Let G/H be a homogeneous space of reductive type. We provide a necessary condition for the properness of the action of L on G/H. As an application we give examples of spaces that do not admit standard compact Clifford-Klein forms.


title: Dependence versus Conditional Dependence in Local Causal Discovery from Gene Expression Data

abstract: Motivation: Algorithms that discover variables which are causally related to a target may inform the design of experiments. With observational gene expression data, many methods discover causal variables by measuring each variable's degree of statistical dependence with the target using dependence measures (DMs). However, other methods measure each variable's ability to explain the statistical dependence between the target and the remaining variables in the data using conditional dependence measures (CDMs), since this strategy is guaranteed to find the target's direct causes, direct effects, and direct causes of the direct effects in the infinite sample limit. In this paper, we design a new algorithm in order to systematically compare the relative abilities of DMs and CDMs in discovering causal variables from gene expression data. Results: The proposed algorithm using a CDM is sample efficient, since it consistently outperforms other state-of-the-art local causal discovery algorithms when samples sizes are small. However, the proposed algorithm using a CDM outperforms the proposed algorithm using a DM only when sample sizes are above several hundred. These results suggest that accurate causal discovery from gene expression data using current CDM-based algorithms requires datasets with at least several hundred samples. Availability: The proposed algorithm is freely available at https://github.com/ericstrobl/DvCD.


title: Dynamical Casimir Effect in Quantum Information Processing

abstract: We demonstrate, in the regime of ultrastrong matter-field coupling, the strong connection between the dynamical Casimir effect (DCE) and the performance of quantum information protocols. Our results are illustrated by means of a realistic quantum communication channel and show that the DCE is a fundamental limit for quantum computation and communication and that novel schemes are required to implement ultrafast and reliable quantum gates. Strategies to partially counteract the DCE are also discussed.


title: Transitive factorizations of permutations and geometry

abstract: We give an account of our work on transitive factorizations of permutations. The work has had impact upon other areas of mathematics such as the enumeration of graph embeddings, random matrices, branched covers, and the moduli spaces of curves. Aspects of these seemingly unrelated areas are seen to be related in a unifying view from the perspective of algebraic combinatorics. At several points this work has intertwined with Richard Stanley's in significant ways.


title: Constructing Krall-Hahn orthogonal polynomials

abstract: Given a sequence of polynomials $(p_n)_n$, an algebra of operators $\mathcal A$ acting in the linear space of polynomials and an operator $D_p\in \mathcal A$ with $D_p(p_n)=\theta_np_n$, where $\theta_n$ is any arbitrary eigenvalue, we construct a new sequence of polynomials $(q_n)_n$ by considering a linear combination of $m+1$ consecutive $p_n$: $q_n=p_n+\sum_{j=1}^m\beta_{n,j}p_{n-j}$. Using the concept of $\mathcal{D}$-operator, we determine the structure of the sequences $\beta_{n,j}, j=1,\ldots,m,$ in order that the polynomials $(q_n)_n$ are eigenfunctions of an operator in the algebra $\mathcal A$. As an application, from the classical discrete family of Hahn polynomials we construct orthogonal polynomials $(q_n)_n$ which are also eigenfunctions of higher-order difference operators.


title: A complexity approach to the soliton resolution conjecture

abstract: The soliton resolution conjecture is one of the most interesting open problems in the theory of nonlinear dispersive equations. Roughly speaking it asserts that a solution with generic initial condition converges to a finite number of solitons plus a radiative term. In this paper we use the complexity of a finite object, a notion introduced in Algorithmic Information Theory, to show that the soliton resolution conjecture is equivalent to the analogous of the second law of thermodynamics for the complexity of a solution of a dispersive equation.


title: Energy current and energy fluctuations in driven quantum wires

abstract: We discuss the energy current and the energy fluctuations in an isolated quantum wire driven far from equilibrium. The system consists of interacting spinless fermions and is driven by a time--dependent magnetic flux. The energy current is defined by the continuity equation for the energy density which is derived both for homogeneous as well as for inhomogeneous systems. Since the total energy is not conserved in the driven system, the continuity equation includes the source terms which are shown to represent the Joule heating effects. For short times and weak drivings the energy current agrees with the linear response theory. For stronger fields or longer times of driving the system enters the quasiequilibrium regime when the energy current gradually diminishes due to the heating effects. Finally, for even stronger driving the energy current is shown to undergo a damped Bloch oscillations. The energy spread also increases upon driving. However, the time--dependence of this quantity in the low field regime is quite unexpected since it is determined mostly by the time of driving being quite independent of the instantaneous energy of the system.