title: SAME but Different: Fast and High-Quality Gibbs Parameter Estimation

abstract: Gibbs sampling is a workhorse for Bayesian inference but has several limitations when used for parameter estimation, and is often much slower than non-sampling inference methods. SAME (State Augmentation for Marginal Estimation) \cite{Doucet99,Doucet02} is an approach to MAP parameter estimation which gives improved parameter estimates over direct Gibbs sampling. SAME can be viewed as cooling the posterior parameter distribution and allows annealed search for the MAP parameters, often yielding very high quality (lower loss) estimates. But it does so at the expense of additional samples per iteration and generally slower performance. On the other hand, SAME dramatically increases the parallelism in the sampling schedule, and is an excellent match for modern (SIMD) hardware. In this paper we explore the application of SAME to graphical model inference on modern hardware. We show that combining SAME with factored sample representation (or approximation) gives throughput competitive with the fastest symbolic methods, but with potentially better quality. We describe experiments on Latent Dirichlet Allocation, achieving speeds similar to the fastest reported methods (online Variational Bayes) and lower cross-validated loss than other LDA implementations. The method is simple to implement and should be applicable to many other models.


title: Deformable Part Models are Convolutional Neural Networks

abstract: Deformable part models (DPMs) and convolutional neural networks (CNNs) are two widely used tools for visual recognition. They are typically viewed as distinct approaches: DPMs are graphical models (Markov random fields), while CNNs are "black-box" non-linear classifiers. In this paper, we show that a DPM can be formulated as a CNN, thus providing a novel synthesis of the two ideas. Our construction involves unrolling the DPM inference algorithm and mapping each step to an equivalent (and at times novel) CNN layer. From this perspective, it becomes natural to replace the standard image features used in DPM with a learned feature extractor. We call the resulting model DeepPyramid DPM and experimentally validate it on PASCAL VOC. DeepPyramid DPM significantly outperforms DPMs based on histograms of oriented gradients features (HOG) and slightly outperforms a comparable version of the recently introduced R-CNN detection system, while running an order of magnitude faster.


title: Motivic classes of some classifying stacks

abstract: We prove that the class of the classifying stack $B PGL_n$ is the multiplicative inverse of the class of the projective linear group $PGL_n$ in the Grothendieck ring of stacks for $n = 2$ and $n = 3$ under mild conditions on the base field $k$. In particular, although it is known that the multiplicativity relation $\{T\} = \{S\}\cdot\{PGL_n\}$ does not hold for all $PGL_n$-torsors $T \to S$, it holds for the universal $PGL_n$-torsors for said $n$.


title: Lattice Structures for Attractors II

abstract: The algebraic structure of the attractors in a dynamical system determine much of its global dynamics. The collection of all attractors has a natural lattice structure, and this structure can be detected through attracting neighborhoods, which can in principle be computed. Indeed, there has been much recent work on developing and implementing general computational algorithms for global dynamics, which are capable of computing attracting neighborhoods efficiently. Here we address the question of whether all of the algebraic structure of attractors can be captured by these methods.


title: The HI Chronicles of LITTLE THINGS BCDs II: The Origin of IC 10's HI Structure

abstract: In this paper we analyze Very Large Array (VLA) telescope and Green Bank Telescope (GBT) atomic hydrogen (HI) data for the LITTLE THINGS(1) blue compact dwarf galaxy IC 10. The VLA data allow us to study the detailed HI kinematics and morphology of IC 10 at high resolution while the GBT data allow us to search the surrounding area at high sensitivity for tenuous HI. IC 10's HI appears highly disturbed in both the VLA and GBT HI maps with a kinematically distinct northern HI extension, a kinematically distinct southern plume, and several spurs in the VLA data that do not follow the general kinematics of the main disk. We discuss three possible origins of its HI structure and kinematics in detail: a current interaction with a nearby companion, an advanced merger, and accretion of intergalactic medium. We find that IC 10 is most likely an advanced merger or a galaxy undergoing accretion. 1:Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey; https://science.nrao.edu/science/surveys/littlethings


title: Exceeding the leading spike intensity and fluence limits in backward Raman amplifiers

abstract: The leading amplified spike in backward Raman amplifiers can reach nearly relativistic intensities before the saturation by the relativistic electron nonlinearity, which sets an upper limit to the largest achievable leading spike intensity. It is shown here that this limit can be substantially exceeded by the initially sub-dominant spikes, which surprisingly outgrow the leading spike after its nonlinear saturation. Furthermore, an initially negligible group velocity of the amplified pulse in strongly under-critical plasma appears to be capable of delaying the filamentation instability in the nonlinear saturation regime. This enables further amplification of the pulse to even larger output fluences.


title: Different approaches to TMD Evolution with scale

abstract: Many apparently contradictory approaches to TMD factorization and its non-perturbative content exist. This talk evaluated the different methods and proposed tools for resolving the contradictions and experimentally adjudicating the results.


title: Fractal analysis of the galaxy distribution in the redshift range 0.45 < z < 5.0

abstract: Evidence is presented that the galaxy distribution can be described as a fractal system in the redshift range of the FDF galaxy survey. The fractal dimension $D$ was derived using the FDF galaxy volume number densities in the spatially homogeneous standard cosmological model with $\Omega_{m_0}=0.3$, $\Omega_{\Lambda_0}=0.7$ and $H_0=70 \; \mbox{km} \; {\mbox{s}}^{-1} \; {\mbox{Mpc}}^{-1}$. The ratio between the differential and integral number densities $\gamma$ and $\gamma^\ast$ obtained from the red and blue FDF galaxies provides a direct method to estimate $D$, implying that $\gamma$ and $\gamma^\ast$ vary as power-laws with the cosmological distances. The luminosity distance $d_{\scriptscriptstyle L}$, galaxy area distance $d_{\scriptscriptstyle G}$ and redshift distance $d_z$ were plotted against their respective number densities to calculate $D$ by linear fitting. It was found that the FDF galaxy distribution is characterized by two single fractal dimensions at successive distance ranges. Two straight lines were fitted to the data, whose slopes change at $z \approx 1.3$ or $z \approx 1.9$ depending on the chosen cosmological distance. The average fractal dimension calculated using $\gamma^\ast$ changes from $\langle D \rangle=1.4^{\scriptscriptstyle +0.7}_{\scriptscriptstyle -0.6}$ to $\langle D \rangle=0.5^{\scriptscriptstyle +1.2}_{\scriptscriptstyle -0.4}$ for all galaxies, and $D$ decreases as $z$ increases. Small values of $D$ at high $z$ mean that in the past galaxies were distributed much more sparsely and the large-scale galaxy structure was then possibly dominated by voids. Results of Iribarrem et al. (2014, arXiv:1401.6572) indicating similar fractal features with $\langle D \rangle =0.6 \pm 0.1$ in the far-infrared sources of the Herschel/PACS evolutionary probe (PEP) at $1.5 \lesssim z \lesssim 3.2$ are also mentioned.


title: The Binomial Theorem and motivic classes of universal quasi-split tori

abstract: Taking symmetric powers of varieties can be seen as a functor from the category of varieties to the category of varieties with an action by the symmetric group. We study a corresponding map between the Grothendieck groups of these categories. In particular, we derive a binomial formula and use it to give explicit expressions for the classes of universal quasi-split tori in the equivariant Grothendieck group of varieties.


title: Normalizations of Eisenstein integrals for reductive symmetric spaces

abstract: We construct minimal Eisenstein integrals for a reductive symmetric space G/H as matrix coefficients of the minimal principal series of G. The Eisenstein integrals thus obtained include those from the \sigma-minimal principal series. In addition, we obtain related Eisenstein integrals, but with different normalizations. Specialized to the case of the group, this wider class includes Harish-Chandra's minimal Eisenstein integrals.


title: Non-Hermitian ${\cal PT}$-Symmetric Dirac-Pauli Hamiltonians with Real Energy Eigenvalues in the Magnetic Field

abstract: The modified Dirac-Pauli equations, which is entered by means of ${\gamma_5}$-mass extension of Hamiltonian operators, are considered. We also take into account the interaction of fermions with the intensive homogenous magnetic field focusing attention on (g-2) gyromagnetic factor of particles with spin $1/2$. Without the use of perturbation theory in the external field the exact energy spectra are deduced with regard to spin effects of fermions. We discuss the possible proposals of experimental measuring of properties of new particles which arising in this model.


title: An in situ measurement of the radio-frequency attenuation in ice at Summit Station, Greenland

abstract: We report an in situ measurement of the electric field attenuation length at radio frequencies for the bulk ice at Summit Station, Greenland, made by broadcasting radio-frequency signals vertically through the ice and measuring the relative power in the return ground bounce signal. We find the depth-averaged field attenuation length to be 947 +92/-85 meters at 75 MHz. While this measurement has clear radioglaciological applications, the radio clarity of the ice also has implications for the detection of ultra-high energy (UHE) astrophysical particles via their radio emission in dielectric media such as ice. Assuming a reliable extrapolation to higher frequencies, the measured attenuation length at Summit Station is comparable to previously measured radio-frequency attenuation lengths at candidate particle detector sites around the world, and strengthens the case for Summit Station as a promising northern site for UHE neutrino detection.


title: Differentially Private Linear Algebra in the Streaming Model

abstract: Numerical linear algebra plays an important role in computer science. In this paper, we initiate the study of performing linear algebraic tasks while preserving privacy when the data is streamed online. Our main focus is the space requirement of the privacy-preserving data-structures. We give the first {\em sketch-based} algorithm for differential privacy. We give optimal, up to logarithmic factor, space data-structures that can compute low rank approximation, linear regression, and matrix multiplication, while preserving differential privacy with better additive error bounds compared to the known results. Notably, we match the best known space bound in the non-private setting by Kane and Nelson (J. ACM, 61(1):4). Our mechanism for differentially private low-rank approximation {\em reuses} the random Gaussian matrix in a specific way to provide a single-pass mechanism. We prove that the resulting distribution also preserve differential privacy. This can be of independent interest. We do not make any assumptions, like singular value separation or normalized row assumption, as made in the earlier works. The mechanisms for matrix multiplication and linear regression can be seen as the private analogues of the known non-private algorithms. All our mechanisms, in the form presented, can also be computed in the distributed setting.


title: Ground state energy of large polaron systems

abstract: The last unsolved problem about the many-polaron system, in the Pekar-Tomasevich approximation, is the case of bosons with the electron-electron Coulomb repulsion of strength exactly 1 (the 'neutral case'). We prove that the ground state energy, for large $N$, goes exactly as $-N^{7/5}$, and we give upper and lower bounds on the asymptotic coefficient that agree to within a factor of $2^{2/5}$.


title: Crowded-field image simulator for WSO-UV/ISSIS: first functional version developed by the Glendama team

abstract: We are developing a web-based interactive software to simulate crowded-field imaging with the Imaging and Slitless Spectrograph Instrument for Surveys (ISSIS) on board the future World Space Observatory - Ultraviolet (WSO-UV). This new tool is aimed to prepare WSO-UV/ISSIS proposals to observe multicomponent targets and dense fields. For a given combination of UV channel, filters and exposure time, the user creates a set of point-like and extended sources (source model). This source model produces a final image, which takes into account a pixelated field of view, a realistic conversion between physical flux and counts per second, the convolution with the expected point spread function, a sky background and noise fluctuations. The current version of the simulator is available at the Gravitational LENses and DArk MAtter (Glendama) website, and it allows users to specify all relevant parameters of each point-like or extended source, drag-and-drop sources by using a mouse or a fingertip/stylus on a touchscreen, change the frame size or the brightness scale, etc.


title: Visualizing operators of coupled spin systems

abstract: The state of quantum systems, their energetics, and their time evolution is modeled by abstract operators. How can one visualize such operators for coupled spin systems? A general approach is presented which consists of several shapes representing linear combinations of spherical harmonics. It is applicable to an arbitrary number of spins and can be interpreted as a generalization of Wigner functions. The corresponding visualization transforms naturally under non-selective spin rotations as well as spin permutations. Examples and applications are illustrated for the case of three spins 1/2.


title: Decoherence imaging of spin ensembles using a scanning single-electron spin in diamond

abstract: The nitrogen-vacancy (NV) defect center in diamond has demonstrated great capability for nanoscale magnetic sensing and imaging for both static and periodically modulated target fields. However, it remains a challenge to detect and image randomly fluctuating magnetic fields. Recent theoretical and numerical works have outlined detection schemes that exploit changes in decoherence of the detector spin as a sensitive measure for fluctuating fields. Here we experimentally monitor the decoherence of a scanning NV center in order to image the fluctuating magnetic fields from paramagnetic impurities on an underlying diamond surface. We detect a signal corresponding to roughly 800 ${\mu}_B$ in 2 s of integration time, without any control on the target spins, and obtain magnetic-field spectral information using dynamical decoupling techniques. The extracted spatial and temporal properties of the surface paramagnetic impurities provide insight to prolonging the coherence of near-surface qubits for quantum information and metrology applications.


title: Footprints and footprint analysis for atmospheric dispersion problems

abstract: Footprint analysis, also known as the study of Influence areas, is a first order method for solving inverse atmospheric dispersion problems. We revisit the concept of footprints giving a rigorous definition of the concept (denoted posterior footprints and posterior zero footprints) in terms of spatio-temporal domains. The notion of footprints is then augmented the to the forward dispersion problem by defining prior footprints and prior zero footprints. We then study how posterior footprints and posterior zero footprints can be combined to reveal more information about the source, and how prior footprints and prior footprints can be combined to yield more information about the measurements.


title: High Power Test Of A 3.9 GHz 5-cell Deflecting-mode Cavity In A Cryogenic Operation

abstract: A 3.9 GHz deflecting mode (pi, TM110) cavity has been long used for six-dimensional phase-space beam manipulation tests at the A0 Photo-Injector Lab (16 MeV) in Fermilab and their extended applications with vacuum cryomodules are currently planned at the Advanced Superconducting Test Accelerator (ASTA) user facility (> 50 MeV). Despite the successful test results, the cavity, however, demonstrated limited RF performance during liquid nitrogen (LN2) ambient operation that was inferior to theoretical prediction. We have been performing full analysis of the designed cavity by analytic calculation and comprehensive system simulation analysis to solve complex thermodynamics and mechanical stresses. The re-assembled cryomodule is currently under the test with a 50 kW klystron at the Fermilab A0 beamline, which will benchmark the modeling analysis. The test result will be used to design vacuum cryomodules for the 3.9 GHz deflecting mode cavity that will be employed at the ASTA facility for beam diagnostics and phase-space control.


title: Operator analysis of physical states on magnetized $T^{2}/Z_{N}$ orbifolds

abstract: We discuss an effective way for analyzing the system on the magnetized twisted orbifolds in operator formalism, especially in the complicated cases $T^{2}/Z_{3}$, $T^{2}/Z_{4}$ and $T^{2}/Z_{6}$. We can obtain the exact and analytical results which can be applicable for any larger values of the quantized magnetic flux M, and show that the (non-diagonalized) kinetic terms are generated via our formalism and the number of the surviving physical states are calculable in a rigorous manner by simply following usual procedures in linear algebra in any case. Our approach is very powerful when we try to examine properties of the physical states on (complicated) magnetized orbifolds $T^{2}/Z_{3}$, $T^{2}/Z_{4}$, $T^{2}/Z_{6}$ (and would be in other cases on higher-dimensional torus) and could be an essential tool for actual realistic model construction based on these geometries.


title: Evidence for Kinetic Limitations as a Controlling Factor of Ge Pyramid Formation: a Study of Structural Features of Ge/Si(001) Wetting Layer Formed by Ge Deposition at Room Temperature Followed by Annealing at 600 {\deg}C

abstract: The article presents an experimental study of an issue of whether the formation of arrays of Ge quantum dots on the Si(001) surface is an equilibrium process or it is kinetically controlled. We deposited Ge on Si(001) at the room temperature and explored crystallization of the disordered Ge film as a result of annealing at 600 {\deg}C. The experiment has demonstrated that the Ge/Si(001) film formed in the conditions of an isolated system consists of the standard patched wetting layer and large droplike clusters of Ge rather than of huts or domes which appear when a film is grown in a flux of Ge atoms arriving on its surface. We conclude that the growth of the pyramids appearing at temperatures greater than 600 {\deg}C is controlled by kinetics rather than thermodynamic equilibrium whereas the wetting layer is an equilibrium structure.


title: A trivariate interpolation algorithm using a cube-partition searching procedure

abstract: In this paper we propose a fast algorithm for trivariate interpolation, which is based on the partition of unity method for constructing a global interpolant by blending local radial basis function interpolants and using locally supported weight functions. The partition of unity algorithm is efficiently implemented and optimized by connecting the method with an effective cube-partition searching procedure. More precisely, we construct a cube structure, which partitions the domain and strictly depends on the size of its subdomains, so that the new searching procedure and, accordingly, the resulting algorithm enable us to efficiently deal with a large number of nodes. Complexity analysis and numerical experiments show high efficiency and accuracy of the proposed interpolation algorithm.


title: Using SOS for Analysis of Zeno Stability in Hybrid systems with Nonlinearity and Uncertainy

abstract: Hybrid systems exhibit phenomena which do not occur in systems with continuous vector fields. One such phenomenon - Zeno executions - is characterized by an infinite number of discrete events or transitions occurring over a finite interval of time. This phenomenon is not necessarily undesirable and may indeed be used to capture physical phenomena. In this paper, we examine the problem of proving the existence and stability of zero executions. Our approach is to develop a polynomial-time algorithm - based on the sum-of-squares methodology - for verifying the stability of a Zeno execution. We begin by stating Lyapunov-like theorems for local Zeno stability based on existing results. Then, for hybrid systems with polynomial vector fields, we use polynomial Lyapunov functions and semialgebraic geometry (Positivstellensatz results) to reduce the local Lyapunov-like conditions to a convex feasibility problem in polynomial variables. The feasibility problem is then tested using an algorithm for sum-of-squares programming - SOSTOOLS. We also extend these results to hybrid system with parametric uncertainty, where the uncertain parameters lie in a semialgebriac set. We also provide several examples illustrating the use of our technique.


title: Sharp entropy decay for hypocoercive and non-symmetric Fokker-Planck equations with linear drift

abstract: We investigate the existence of steady states and exponential decay for hypocoercive Fokker--Planck equations on the whole space with drift terms that are linear in the position variable. For this class of equations, we first establish that hypoellipticity of its generator and confinement of the system is equivalent to the existence of a unique normalised steady state. These two conditions also imply hypocoercivity, i.e. exponential convergence of the solution to equilibrium. Since the standard entropy method does not apply to degenerate parabolic equations, we develop a new modified entropy method (based on a modified, non-degenerate entropy dissipation-like functional) to prove this exponential decay in relative entropy (logarithmic till quadratic) - with a sharp rate. Furthermore, we compute the spectrum and eigenspaces of the generator as well as flow-invariant manifolds of Gaussian functions. Next, we extend our method to kinetic Fokker-Planck equations with a class of non-quadratic potentials. And, finally, we apply this new method to non-symmetric, uniformly parabolic Fokker-Planck equations with linear drift. At least in 2D this always yields the sharp exponential envelopes for the entropy function. In this case, we obtain even a sharp multiplicative constant in the decay estimate for the non-symmetric semigroup.


title: Issues and R&D Required for the Intensity Frontier Accelerators

abstract: Operation, upgrade and development of accelerators for Intensity Frontier face formidable challenges in order to satisfy both the near-term and long-term Particle Physics program. Here we discuss key issues and R&D required for the Intensity Frontier accelerators.


title: Modular matrices from universal wave function overlaps in Gutzwiller-projected parton wave functions

abstract: We implement the universal wave function overlap (UWFO) method to extract modular $S$ and $T$ matrices for topological orders in Gutzwiller-projected parton wave functions (GPWFs). The modular $S$ and $T$ matrices generate a projective representation of $SL(2,\mathbb{Z})$ on the degenerate-ground-state Hilbert space on a torus and may fully characterize the 2+1D topological orders, i.e. the quasi-particle statistics and chiral central charge (up to $E_8$ bosonic quantum Hall states). We used the variational Monte Carlo method to computed the $S$ and $T$ matrices of the chiral spin liquid (CSL) constructed by the GPWF on the square lattice, and confirm that the CSL carries the same topological order as the $\nu=\frac{1}{2}$ bosonic Laughlin state. We find that the non-universal exponents in UWFO can be small and direct numerical computation is able to be applied on relatively large systems. We also discuss the UWFO method for GPWFs on other Bravais lattices in two and three dimensions by using the Monte Carlo method. UWFO may be a powerful method to calculate the topological order in GPWFs.


title: Connecting the interstellar magnetic field at the heliosphere to the Loop I superbubble

abstract: The local interstellar magnetic field affects both the heliosphere and the surrounding cluster of interstellar clouds (CLIC). Measurements of linearly polarized starlight provide the only test of the magnetic field threading the CLIC. Polarization measurements of the CLIC magnetic field show multiple local magnetic structures, one of which is aligned with the magnetic field traced by the center of the "ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX). Comparisons between the bulk motion of the CLIC through the local standard of rest, the magnetic field direction, the geometric center of Loop I, and the polarized dust bridge extending from the heliosphere toward the North Polar Spur direction all suggest that the CLIC is part of the rim region of the Loop I superbubble.


title: Antifragile Electronic Warfare

abstract: This letter introduces the concept of antifragile electronic warfare (EW), which we define as the ability to allow a communications link to improve performance due to the presence of a jammer. This concept should not be confused with jamming countermeasures (a.k.a. anti-jamming or electronic protection). Rather, antifragile EW can be thought of as the next step beyond simply avoiding or mitigating jamming. After introducing the concept we narrow down the subset of jammers this concept can be applied to, and provide a brief example of an antifragile EW strategy.


title: Density-wave instabilities of fractionalized Fermi liquids

abstract: Recent experiments in the underdoped regime of the hole-doped cuprates have found evidence for an incommensurate charge density wave state. We present an analysis of the charge ordering instabilities in a metal with antiferromagnetic correlations, where the electronic excitations are coupled to the fractionalized excitations of a quantum fluctuating antiferromagnet on the square lattice. The resulting charge density wave state emerging out of such a fractionalized Fermi-liquid (FL*) has wavevectors of the form $(\pm Q_0,0), (0,\pm Q_0)$, with a predominantly $d$-form factor, in agreement with experiments on a number of different families of the cuprates. In contrast, as previously shown, the charge density wave instability of a nearly antiferromagnetic metal with a large Fermi surface, interacting via short-range interactions, has wavevectors of the type $(\pm Q_0,\pm Q_0)$. Our results show that the observed charge density wave appears as a low-energy instability of a fractionalized metallic state linked to the proximity to an antiferromagnetic insulator, and the pseudogap regime can be described by such a metal at least over intermediate length and energy scales.


title: Three-Dimensional Simulations of Core-Collapse Supernovae: From Shock Revival to Shock Breakout

abstract: We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.


title: Formation of Dark Matter Torii Around Supermassive Black Holes Via The Eccentric Kozai-Lidov Mechanism

abstract: We explore the effects of long term secular perturbations on the distribution of dark matter particles around Supermassive Black Hole (BH) binaries. We show that in the hierarchical (in separation) three-body problem, one of the BHs and a dark matter particle form an inner binary. Gravitational perturbations from the BH companion, on a much wider orbit, can cause the dark matter particle to reach extremely high eccentricities and even get accreted onto the BH, by what is known as the Eccentric Kozai-Lidov (EKL) mechanism. We show that this may produce a torus-like configuration for the dark matter distribution around the less massive member of the BH binary. We first consider an Intermediate BH (IMBH) in the vicinity of our Galactic Center, which may be a relic of a past minor merger. We show that if the IMBH is close enough (i.e., near the stellar disk) the EKL mechanism is very efficient in exciting the eccentricity of dark matter particles in near-polar configurations to extremely high values where they are accreted by the IMBH. We show that this mechanism is even more effective if the central BH grows in mass, where we have assumed adiabatic growth. Since near-polar configurations are disrupted, a torus-like shape is formed. We also show that this behavior is also likely to be relevant for Supermassive BH binaries. We suggest that if the BHs are spinning, the accreted dark matter particles may linger in the ergosphere and thereby may generate self-annihilations and produce an indirect signature of potential interest.


title: The Herschel view of the dominant mode of galaxy growth from z=4 to the present day

abstract: We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.


title: Galaxy Zoo: Are Bars Responsible for the Feeding of Active Galactic Nuclei at 0.2 < z < 1.0?

abstract: We present a new study investigating whether active galactic nuclei (AGN) beyond the local universe are preferentially fed via large-scale bars. Our investigation combines data from Chandra and Galaxy Zoo: Hubble (GZH) in the AEGIS, COSMOS, and GOODS-S surveys to create samples of face-on, disc galaxies at 0.2 < z < 1.0. We use a novel method to robustly compare a sample of 120 AGN host galaxies, defined to have 10^42 erg/s < L_X < 10^44 erg/s, with inactive control galaxies matched in stellar mass, rest-frame colour, size, Sersic index, and redshift. Using the GZH bar classifications of each sample, we demonstrate that AGN hosts show no statistically significant enhancement in bar fraction or average bar likelihood compared to closely-matched inactive galaxies. In detail, we find that the AGN bar fraction cannot be enhanced above the control bar fraction by more than a factor of two, at 99.7% confidence. We similarly find no significant difference in the AGN fraction among barred and non-barred galaxies. Thus we find no compelling evidence that large-scale bars directly fuel AGN at 0.2<z<1.0. This result, coupled with previous results at z=0, implies that moderate-luminosity AGN have not been preferentially fed by large-scale bars since z=1. Furthermore, given the low bar fractions at z>1, our findings suggest that large-scale bars have likely never directly been a dominant fueling mechanism for supermassive black hole growth.


title: CoMaLit III. Literature Catalogs of weak Lensing Clusters of galaxies (LC^2)

abstract: The measurement of the mass of clusters of galaxies is crucial for their use in cosmology and astrophysics. Masses can be efficiently determined with weak lensing (WL) analyses. I compiled Literature Catalogs of weak Lensing Clusters (LC$^2$). Cluster identifiers, coordinates, and redshifts have been standardised. WL masses were reported to over-densities of 2500, 500, 200, and to the virial one in the reference $\Lambda$CDM model. Duplicate entries were carefully handled. I produced three catalogs: LC$^2$-single, with 485 unique groups and clusters analysed with the single-halo model; LC$^2$-substructure, listing substructures in complex systems; LC$^2$-all, listing all the 822 WL masses found in literature. The catalogs and future updates are publicly available at http://pico.bo.astro.it/~sereno/CoMaLit/LC2/.


title: Classifying symmetry-protected topological phases through the anomalous action of the symmetry on the edge

abstract: It is well known that (1+1)-D bosonic symmetry-protected topological (SPT) phases with symmetry group $G$ can be identified by the projective representation of the symmetry at the edge. Here, we generalize this result to higher dimensions. We assume that the representation of the symmetry on the spatial edge of a ($d+1$)-D SPT is /local/ but not necessarily /on-site/, such that there is an obstruction to its implementation on a region with boundary. We show that such obstructions are classified by the cohomology group $H^{d+1}(G, U(1))$, in agreement with the classification of bosonic SPT phases proposed in [Chen et al, Science 338, 1604 (2012)]. Our analysis allows for a straightforward calculation of the element of $H^{d+1}(G, U(1))$ corresponding to physically meaningful models such as non-linear sigma models with a theta term in the action. SPT phases outside the classification of Chen et al are those in which the symmetry cannot be represented locally on the edge. With some modifications, our framework can also be applied to fermionic systems in (2+1)-D.


title: LOFAR low-band antenna observations of the 3C295 and Bootes fields: source counts and ultra-steep spectrum sources

abstract: We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.


title: Simultaneous B and L Violation: New Signatures from RPV-SUSY

abstract: Studies of R-parity violating (RPV) supersymmetry typically assume that nucleon stability is protected by approximate baryon number (B) or lepton number (L) conservation. We present a new class of RPV models that violate B and L simultaneously (BLRPV), without inducing rapid nucleon decay. These models feature an approximate $Z_2^e \times Z_2^\mu \times Z_2^\tau$ flavor symmetry, which forbids 2-body nucleon decay and ensures that flavor antisymmetric $L L E^c$ couplings are the only non-negligible L-violating operators. Nucleons are predicted to decay through $N \rightarrow K e \mu \nu$ and $n \rightarrow e \mu \nu$; the resulting bounds on RPV couplings are rather mild. Novel collider phenomenology arises because the superpartners can decay through both L-violating and B-violating couplings. This can lead to, for example, final states with high jet multiplicity and multiple leptons of different flavor, or a spectrum in which depending on the superpartner, either B or L violating decays dominate. BLRPV can also provide a natural setting for displaced $\tilde{\nu} \rightarrow \mu e$ decays, which evade many existing collider searches for RPV supersymmetry.


title: Updated fit to three neutrino mixing: status of leptonic CP violation

abstract: We present a global analysis of solar, atmospheric, reactor and accelerator neutrino data in the framework of three-neutrino oscillations based on data available in summer 2014. We provide the allowed ranges of the six oscillation parameters and show that their determination is stable with respect to uncertainties related to reactor neutrino and solar neutrino flux predictions. We find that the maximal possible value of the Jarlskog invariant in the lepton sector is $0.0329 \pm 0.0009$ ($\pm 0.0027$) at the $1\sigma$ ($3\sigma$) level and we use leptonic unitarity triangles to illustrate the ability of global oscillation data to obtain information on CP violation. We discuss "tendencies and tensions" of the global fit related to the octant of $\theta_{23}$ as well as the CP violating phase $\delta_\mathrm{CP}$. The favored values of $\delta_\mathrm{CP}$ are around $3\pi/2$ while values around $\pi/2$ are disfavored at about $\Delta\chi^2 \simeq 6$. We comment on the non-trivial task to assign a confidence level to this $\Delta\chi^2$ value by performing a Monte Carlo study of T2K data.


title: Non-Gaussian Shape Discrimination with Spectroscopic Galaxy Surveys

abstract: [Abridged] We consider how galaxy clustering data, from Mpc to Gpc scales, from upcoming large scale structure surveys, such as Euclid and DESI, can provide discriminating information about the bispectrum shape arising from a variety of inflationary scenarios. Through exploring in detail the weighting of shape properties in the calculation of the halo bias and halo mass function we show how they probe a broad range of configurations, beyond those in the squeezed limit, that can help distinguish between shapes with similar large scale bias behaviors. We assess the impact, on constraints for a diverse set of non-Gaussian shapes, of galaxy clustering information in the mildly non-linear regime, and surveys that span multiple redshifts and employ different galactic tracers of the dark matter distribution. Fisher forecasts are presented for a Euclid-like spectroscopic survey of H$\alpha$-selected emission line galaxies (ELGs) using recent revisions of the expected H$\alpha$ luminosity function, and a DESI-like survey, of luminous red galaxies (LRGs) and [O-II] doublet-selected ELGs, in combination with Planck-like CMB temperature and polarization data. While ELG samples provide better probes of shapes that are divergent in the squeezed limit, LRG constraints, centered below $z<1$, yield stronger constraints on shapes with scale-independent large-scale halo biases, such as the equilateral template. The ELG and LRG samples provide complementary degeneracy directions for distinguishing between different shapes. If the Gaussian galaxy bias is constrained to better than a percent level, such as can be determined from the galaxy bispectrum or weak lensing, then the LSS and CMB data could provide complementary constraints that will enable differentiation of bispectra with distinct theoretical origins but with similar large scale, squeezed-limit properties.


title: Polar Kerr effect from chiral-nematic charge order

abstract: We analyze the polar Kerr effect in an itinerant electron system on a square lattice in the presence of a composite charge order proposed for the pseudogap state in underdoped cuprates. This composite charge order preserves discrete translational symmetries, and is "chiral-nematic" in the sense that it breaks time-reversal symmetry, mirror symmetries in $x$ and $y$ directions, and $C_4$ lattice rotation symmetry. The Kerr angle $\theta_K$ in $C_4$-symmetric system is proportional to the antisymmetric component of the anomalous Hall conductivity $\sigma_{xy}-\sigma_{yx}$. We show that this result holds when $C_4$ symmetry is broken. We show that in order for $\sigma_{xy}$ and $\sigma_{yx}$ to be non-zero the mirror symmetries in $x$ and $y$ directions have to be broken, and that for $\sigma_{xy}-\sigma_{yx}$ to be non-zero time-reversal symmetry has to be broken. The chiral-nematic charge order satisfies all these conditions, such that a non-zero signal in a polar Kerr effect experiment is symmetry allowed. We further show that to get a non-zero $\theta_K$ in a one-band spin-fluctuation scenario, in the absence of disorder, one has to extend the spin-mediated interaction to momenta away from $(\pi,\pi)$ and has to include particle-hole asymmetry. Alternatively, in the presence of disorder one can get a non-zero $\theta_K$ from impurity scattering: either due to skew scattering (with non-Gaussian disorder) or due to particle-hole asymmetry in case of Gaussian disorder. The impurity analysis in our case is similar to that in earlier works on Kerr effect in $p_x+ip_y$ superconductor, however in our case the magnitude of $\theta_K$ is enhanced by the flattening of the Fermi surface in the "hot" regions which mostly contribute to charge order.


title: About the Absence of Exotics and the Coulomb Branch Formula

abstract: The absence of exotics is a conjectural property of the spectrum of BPS states of four--dimensional $\mathcal{N}=2$ supersymmetric QFT's. In this letter we revisit the precise statement of this conjecture, and develop a general strategy that, if applicable, entails the absence of exotic BPS states. Our method is based on the Coulomb branch formula and on quiver mutations. In particular, we obtain the absence of exotic BPS states for all pure SYM theories with simple simply--laced gauge group $G$, and, as a corollary, of infinitely many other lagrangian $\mathcal{N}=2$ theories.


title: Exploratory Analysis of a Terabyte Scale Web Corpus

abstract: In this paper we present a preliminary analysis over the largest publicly accessible web dataset: the Common Crawl Corpus. We measure nine web characteristics from two levels of granularity using MapReduce and we comment on the initial observations over a fraction of it. To the best of our knowledge two of the characteristics, the language distribution and the HTML version of pages have not been analyzed in previous work, while the specific dataset has been only analyzed on page level.


title: Darboux Transformations for (2+1)-Dimensional Extensions of the KP Hierarchy

abstract: New extensions of the KP and modified KP hierarchies with self-consistent sources are proposed. The latter provide new generalizations of $(2+1)$-dimensional integrable equations, including the DS-III equation and the $N$-wave problem. Furthermore, we recover a system that contains two types of the KP equation with self-consistent sources as special cases. Darboux and binary Darboux transformations are applied to generate solutions of the proposed hierarchies.


title: ASTA at Fermilab: Accelerator Physics and Accelerator Education Programs at the Modern Accelerator R&D Users Facility for HEP and Accelerator Applications

abstract: We present the current and planned beam physics research program and accelerator education program at Advanced Superconducting Test Accelerator (ASTA) at Fermilab.


title: Quantifying (dis)agreement between direct detection experiments in a halo-independent way

abstract: We propose an improved method to study recent and near-future dark matter direct detection experiments with small numbers of observed events. Our method determines in a quantitative and halo-independent way whether the experiments point towards a consistent dark matter signal and identifies the best-fit dark matter parameters. To achieve true halo independence, we apply a recently developed method based on finding the velocity distribution that best describes a given set of data. For a quantitative global analysis we construct a likelihood function suitable for small numbers of events, which allows us to determine the best-fit particle physics properties of dark matter considering all experiments simultaneously. Based on this likelihood function we propose a new test statistic that quantifies how well the proposed model fits the data and how large the tension between different direct detection experiments is. We perform Monte Carlo simulations in order to determine the probability distribution function of this test statistic and to calculate the p-value for both the dark matter hypothesis and the background-only hypothesis.


title: Event-Horizon-Telescope Evidence for Alignment of the Black Hole in the Center of the Milky Way with the Inner Stellar Disk

abstract: Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.


title: CIV and CIII] reverberation mapping of the luminous quasar PG 1247+267

abstract: So far the masses of about 50 active galactic nuclei have been measured through the reverberation mapping technique (RM). Most measurements have been performed for objects of moderate luminosity and redshift, based on H$\beta$, which is also used to calibrate the scaling relation which allows single-epoch (SE) mass determination based on AGN luminosity and the width of different emission lines. The SE mass obtained from CIV$(1549 {\rm\AA})$ line shows a large spread around mean values, due to complex structure and gas dynamics of the relevant emission region. Direct RM measures of CIV exist for only 6 AGNs of low luminosity and redshift, and only one luminous quasar (Kaspi et al 2007). We have collected since 2003 photometric and spectroscopic observations of PG1247+267, the most luminous quasar ever analyzed for RM. We provide light curves for the continuum and for CIV$(1549 {\rm\AA})$ and CIII]$(1909{\rm\AA})$, and measures of the reverberation time lags based on the SPEAR method (Zu et al. 2011). The sizes of the line emission regions are in a ratio $R_{CIII]}/R_{CIV}\sim 2$, similar to the case of Seyfert galaxies, indicating for the first time a similar ionization stratification in a luminous quasar and low luminosity nuclei. Due to relatively small broad line region size and relatively narrow line widths, we estimate a small mass and an anomalously high Eddington ratio. We discuss the possibility that either the shape of the emission region or an amplification of the luminosity caused by gravitational lensing may be in part responsible of the result.


title: Jet substructure and probes of CP violation in Vh production

abstract: We analyse the hVV (V = W, Z) vertex in a model independent way using Vh production. To that end, we consider possible corrections to the Standard Model Higgs Lagrangian, in the form of higher dimensional operators which parametrise the effects of new physics. In our analysis, we pay special attention to linear observables that can be used to probe CP violation in the same. By considering the associated production of a Higgs boson with a vector boson (W or Z), we use jet substructure methods to define angular observables which are sensitive to new physics effects, including an asymmetry which is linearly sensitive to the presence of CP odd effects. We demonstrate how to use these observables to place bounds on the presence of higher dimensional operators, and quantify these statements using a log likelihood analysis. Our approach allows one to probe separately the hZZ and hWW vertices, involving arbitrary combinations of BSM operators, at the Large Hadron Collider.


title: Improving Reliability of Subject-Level Resting-State fMRI Parcellation with Shrinkage Estimators

abstract: A recent interest in resting state functional magnetic resonance imaging (rsfMRI) lies in subdividing the human brain into anatomically and functionally distinct regions of interest. For example, brain parcellation is often used for defining the network nodes in connectivity studies. While inference has traditionally been performed on group-level data, there is a growing interest in parcellating single subject data. However, this is difficult due to the low signal-to-noise ratio of rsfMRI data, combined with typically short scan lengths. A large number of brain parcellation approaches employ clustering, which begins with a measure of similarity or distance between voxels. The goal of this work is to improve the reproducibility of single-subject parcellation using shrinkage estimators of such measures, allowing the noisy subject-specific estimator to "borrow strength" in a principled manner from a larger population of subjects. We present several empirical Bayes shrinkage estimators and outline methods for shrinkage when multiple scans are not available for each subject. We perform shrinkage on raw intervoxel correlation estimates and use both raw and shrinkage estimates to produce parcellations by performing clustering on the voxels. Our proposed method is agnostic to the choice of clustering method and can be used as a pre-processing step for any clustering algorithm. Using two datasets---a simulated dataset where the true parcellation is known and is subject-specific and a test-retest dataset consisting of two 7-minute rsfMRI scans from 20 subjects---we show that parcellations produced from shrinkage correlation estimates have higher reliability and validity than those produced from raw estimates. Application to test-retest data shows that using shrinkage estimators increases the reproducibility of subject-specific parcellations of the motor cortex by up to 30%.


title: On Hilbert bases of cuts

abstract: A Hilbert basis is a set of vectors X such that the integer cone (semigroup) generated by X is the intersection of the lattice generated by X with the cone generated by X. Define a graph to be (cut) Hilbert if its set of cuts forms a Hilbert basis. We show that the Hilbert property is not closed under edge deletions, subdivisions, nor 2-sums. Furthermore, no graph having K_6-e as a minor is Hilbert. This corrects an error in [M. Laurent. Hilbert bases of cuts. Discrete Math., 150(1-3):257-279 (1996)]. For positive results, we give conditions under which the 2-sum of two graphs produces a Hilbert graph. Using these conditions we show that all H-minor-free graphs are Hilbert , where H is the unique 3-connected graph obtained by uncontracting an edge of K_5. We also establish a relationship between edge deletion and subdivision. Namely, if G' is obtained from a Hilbert graph G by subdividing an edge e two or more times, then G-e is Hilbert if and only if G' is Hilbert.


title: Windows into Geometric Events: Data Structures for Time-Windowed Querying of Temporal Point Sets

abstract: We study geometric data structures for sets of point-based temporal events, answering time-windowed queries, i.e., given a contiguous time interval we answer common geometric queries about the point events with time stamps in this interval. The geometric queries we consider include queries based on the skyline, convex hull, and proximity relations of the point set. We provide space efficient data structures which answer queries in polylogarithmic time.


title: Simultaneous submicrometric 3D imaging of the micro-vascular network and the neuronal system in a mouse spinal cord

abstract: Defaults in vascular (VN) and neuronal networks of spinal cord are responsible for serious neurodegenerative pathologies. Because of inadequate investigation tools, the lacking knowledge of the complete fine structure of VN and neuronal systems is a crucial problem. Conventional 2D imaging yields incomplete spatial coverage leading to possible data misinterpretation, whereas standard 3D computed tomography imaging achieves insufficient resolution and contrast. We show that X-ray high-resolution phase-contrast tomography allows the simultaneous visualization of three-dimensional VN and neuronal systems of mouse spinal cord at scales spanning from millimeters to hundreds of nanometers, with neither contrast agent nor a destructive sample-preparation. We image both the 3D distribution of micro-capillary network and the micrometric nerve fibers, axon-bundles and neuron soma. Our approach is a crucial tool for pre-clinical investigation of neurodegenerative pathologies and spinal-cord-injuries. In particular, it should be an optimal tool to resolve the entangled relationship between VN and neuronal system.


title: Simulating (p)reheating after inflation via the DCE?

abstract: We note some close parallels between preheating/perturbative reheating, (p)reheating, models in post-inflationary cosmology and the dynamical Casimir effect (DCE) in quantum optics. For the plasma-mirror model we show how the effective plasma mass (arising from conduction electrons) behaves like an oscillating inflaton field, while created photons behave like a scalar field coupled quadratically to the inflaton. Furthermore, the effect of spacetime expansion can also be incorporated by varying the dielectric function. We propose an experiment that could mimic (p)reheating for both narrow and broad parametric resonance, by employing technology already being used in attempts to detect DCE photons via plasma-mirrors.


title: Permutation Entropy and Statistical Complexity Analysis of Turbulence in Laboratory Plasmas and the Solar Wind

abstract: The Bandt-Pompe permutation entropy and the Jensen-Shannon statistical complexity are used to analyze fluctuating time series of three different plasmas: the magnetohydrodynamic (MHD) turbulence in the plasma wind tunnel of the Swarthmore Spheromak Experiment (SSX), drift-wave turbulence of ion saturation current fluctuations in the edge of the Large Plasma Device (LAPD) and fully-developed turbulent magnetic fluctuations of the solar wind taken from the WIND spacecraft. The entropy and complexity values are presented as coordinates on the CH plane for comparison among the different plasma environments and other fluctuation models. The solar wind is found to have the highest permutation entropy and lowest statistical complexity of the three data sets analyzed. Both laboratory data sets have larger values of statistical complexity, suggesting these systems have fewer degrees of freedom in their fluctuations, with SSX magnetic fluctuations having slightly less complexity than the LAPD edge fluctuations. The CH plane coordinates are compared to the shape and distribution of a spectral decomposition of the waveforms. These results suggest that fully developed turbulence (solar wind) occupies the lower-right region of the CH plane, and that other plasma systems considered to be turbulent have less permutation entropy and more statistical complexity. This paper presents the first use of this statistical analysis tool on solar wind plasma, as well as on an MHD turbulent experimental plasma.


title: Gemini Planet Imager Spectroscopy of the HR 8799 planets c and d

abstract: During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.


title: Collective Transport Properties of Driven Skyrmions with Random Disorder

abstract: We use particle-based simulations to examine the static and driven collective phases of skyrmions interacting with random quenched disorder. We show that non-dissipative effects due to the Magnus term reduce the depinning threshold and strongly affect the skyrmion motion and the nature of the dynamic phases. The quenched disorder causes the Hall angle to become drive-dependent in the moving skyrmion phase, while different flow regimes produce distinct signatures in the transport curves. For weak disorder, the skyrmions form a pinned crystal and depin elastically, while for strong disorder the system forms a pinned amorphous state that depins plastically. At high drives the skyrmions can dynamically reorder into a moving crystal, with the onset of reordering determined by the strength of the Magnus term.


title: Test of Optical Stochastic Cooling in the IOTA Ring

abstract: A new 150 MeV electron storage ring is being built at Fermilab. The construction of a new machine pursues two goals: A test of highly non-linear integrable optics and a test of optical stochastic cooling. This paper discusses details of OSC arrangements, choice of major parameters of the cooling scheme and incoming experimental tests of the optical amplifier prototype which uses highly doped Ti-sapphire crystal as amplification medium.


title: The molecular clock of neutral evolution can be accelerated or slowed by asymmetric spatial structure

abstract: Over time, a population acquires neutral genetic substitutions as a consequence of random drift. A famous result in population genetics asserts that the rate, $K$, at which these substitutions accumulate in the population coincides with the mutation rate, $u$, at which they arise in individuals: $K=u$. This identity enables genetic sequence data to be used as a "molecular clock" to estimate the timing of evolutionary events. While the molecular clock is known to be perturbed by selection, it is thought that $K=u$ holds very generally for neutral evolution. Here we show that asymmetric spatial population structure can alter the molecular clock rate for neutral mutations, leading to either $K<u$ or $K>u$. Deviations from $K=u$ occur because mutations arise unequally at different sites and have different probabilities of fixation depending on where they arise. If birth rates are uniform across sites, then $K \leq u$. In general, $K$ can take any value between 0 and $Nu$. Our model can be applied to a variety of population structures. In one example, we investigate the accumulation of genetic mutations in the small intestine. In another application, we analyze over 900 Twitter networks to study the effect of network topology on the fixation of neutral innovations in social evolution.


title: Spin-transfer torque induced spin waves in antiferromagnetic insulators

abstract: We explore the possibility of exciting spin waves in insulating antiferromagnetic films by injecting spin current at the surface. We analyze both magnetically compensated and uncompensated interfaces. We find that the spin current induced spin-transfer torque can excite spin waves in insulating antiferromagnetic materials and that the chirality of the excited spin wave is determined by the polarization of the injected spin current. Furthermore, the presence of magnetic surface anisotropy can greatly increase the accessibility of these excitations.


title: Storage Ring And Interaction Region Magnets For A {\mu}+{\mu}- Higgs Factory

abstract: A low-energy Muon Collider (MC) offers unique opportunities to study the recently found Higgs boson. However, due to a relatively large beam emittance with moderate cooling in this machine, large-aperture high- field superconducting (SC) magnets are required. The magnets need also an adequate margin to operate at a large radiation load from the muon decay showers. General specifications of the SC dipoles and quadrupoles for the 125 GeV c.o.m. Higgs Factory with an average luminosity of ~2x10**31 cm-2s-1 are formulated. Magnet conceptual designs and parameters are reported. The impact of the magnet fringe fields on the beam dynamics as well as the IR and lattice magnet protection from radiation are also reported and discussed.


title: Fast and deterministic computation of the determinant of a polynomial matrix

abstract: Given a square, nonsingular matrix of univariate polynomials $\mathbf{F}\in\mathbb{K}[x]^{n\times n}$ over a field $\mathbb{K}$, we give a deterministic algorithm for finding the determinant of $\mathbf{F}$. The complexity of the algorithm is $\bigO \left(n^{\omega}s\right)$ field operations where $s$ is the average column degree or the average row degree of $\mathbf{F}$. Here $\bigO$ notation is Big-$O$ with log factors omitted and $\omega$ is the exponent of matrix multiplication.


title: Some spectral properties for generalized derivations

abstract: Given Banach spaces $\mathcal{X}$ and $\mathcal{Y}$ and Banach space operators $A\in L(\mathcal{X})$ and $B\in L(\mathcal{Y}).$ The generalized derivation $\delta_{A,B} \in L(L(\mathcal{Y},\mathcal{X}))$ is defined by $\delta_{A,B}(X)=(L_{A}-R_{B})(X)=AX-XB$. This paper is concerned with the problem of the transferring the left polaroid property, from operators $A$ and $B^{*}$ to the generalized derivation $\delta_{A,B}$. As a consequence, we give necessary and sufficient conditions for $\delta_{A,B}$ to satisfy generalized a-Browder's theorem and generalized a-Weyl's theorem. As application, we extend some recent results concerning Weyl type theorems.


title: High-energy Particle Colliders: Past 20 Years, Next 20 Years, And Beyond

abstract: Particle colliders for high-energy physics have been in the forefront of scientific discoveries for more than half a century. The accelerator technology of the colliders has progressed immensely, while the beam energy, luminosity, facility size, and cost have grown by several orders of magnitude. The method of colliding beams has not fully exhausted its potential but has slowed down considerably in its progress. This paper briefly reviews the colliding beam method and the history of colliders, discusses the development of the method over the last two decades in detail, and examines near-term collider projects that are currently under development. The paper concludes with an attempt to look beyond the current horizon and to find what paradigm changes are necessary for breakthroughs in the field.


title: Harish-Chandra bimodules over rational Cherednik algebras

abstract: We study Harish-Chandra bimodules over the rational Cherednik algebra $H_{c}(W)$ associated to a complex reflection group $W$ with parameter $c$. Our results allow us to partially reduce the study of these bimodules to smaller algebras. We classify those pairs of parameters $(c,c')$ for which there exist fully supported Harish-Chandra bimodules, and give a description of the category of all Harish-Chandra bimodules modulo those without full support. When $W$ is a symmetric group we are able to classify all irreducible Harish-Chandra bimodules. Our proofs are based on localization techniques, the action of the Namikawa-Weyl group on the set of parameters, and the study of partial KZ functors.


title: Higher-Order Triangular-Distance Delaunay Graphs: Graph-Theoretical Properties

abstract: We consider an extension of the triangular-distance Delaunay graphs (TD-Delaunay) on a set $P$ of points in the plane. In TD-Delaunay, the convex distance is defined by a fixed-oriented equilateral triangle $\triangledown$, and there is an edge between two points in $P$ if and only if there is an empty homothet of $\triangledown$ having the two points on its boundary. We consider higher-order triangular-distance Delaunay graphs, namely $k$-TD, which contains an edge between two points if the interior of the homothet of $\triangledown$ having the two points on its boundary contains at most $k$ points of $P$. We consider the connectivity, Hamiltonicity and perfect-matching admissibility of $k$-TD. Finally we consider the problem of blocking the edges of $k$-TD.


title: Monopole correlations in holographically flavored liquids

abstract: Many-body systems with a conserved U(1) current in (2+1) dimensions may be probed by weakly gauging this current and studying correlation functions of magnetic monopole operators in the resulting dynamical gauge theory. We study such monopole correlations in holographic liquids with fundamental flavor, where the monopole operator is dual to a magnetically charged particle in the bulk. In charge-gapped phases the monopole operator is expected to condense. We show that this condensation is holographically dual to the capping off of the bulk flavor brane and compute the monopole condensate. We argue that from the lower-dimensional point of view this may be understood as a simple example of confinement of a gauge field in the bulk. In a compressible finite-density phase we present a novel calculation of the monopole correlation in space and time: the correlation is power law in time but is Gaussian in space due to interaction with the background charge density.


title: The Star Formation Relation for Regions in the Galactic Plane: The Effect of Spatial Resolution

abstract: We examined the relations between molecular gas surface density and star formation rate surface density in a 11 square degree region of the Galactic Plane. Dust continuum at 1.1 mm from the Bolocam Galactic Plane Survey and 22 micron emission from the WISE All-sky survey were used as tracers of molecular gas and star formation rate, respectively, across Galactic longitude of 31.5 > l > 20.5 and Galactic latitude of 0.5 > b > -0.5. The relation was studied over a range of resolutions from 33 arcsecond to 20 arcminute by convolving images to larger scales. The pixel-by-pixel correlation between 1.1 mm and 22 micron increases rapidly at small scales and levels off at the scale of 5 - 8 arcminute. We studied the star formation relation based on pixel-by-pixel analysis and 1.1 mm and 22 micron peaks analysis. The star formation relation was found to be nearly linear with no significant changes in the form of the relation across all spatial scales and lie above the extragalactic relation from Kennicutt (1998). The average gas depletion time is approximately 200 Myr and does not change significantly at different scales, but the scatter in the depletion time decreases as the scale increases.


title: The nature of the glass and gel transitions in sticky spheres

abstract: Glasses and gels are the two dynamically arrested, disordered states of matter. Despite their importance, their similarities and differences remain elusive, especially at high density. We identify dynamical and structural signatures which distinguish the gel and glass transitions in a colloidal model system of hard and "sticky" spheres. Gelation is induced by crossing the gas-liquid phase-separation line and the resulting rapid densification of the colloid-rich phase leads to a sharp change in dynamics and local structure. Thus, we find that gelation is first-order-like and can occur at much higher densities than previously thought: far from being low-density networks, gels have a clear "thermodynamic" definition which nevertheless leads to a non-equilibrium state with a distinct local structure characteristic of a rapidly quenched glass. In contrast, approaching the glass transition, the dynamics slow continuously accompanied by the emergence of local five-fold symmetric structure. Our findings provide a general thermodynamic, kinetic, and structural basis upon which to distinguish gelation from vitrification.


title: The Project-x Injector Experiment: A Novel High Performance Front-end For A Future High Power Proton Facility At Fermilab

abstract: A multi-MW proton facility, Project X, has been proposed and is currently under development at Fermilab. We are carrying out a program of research and development aimed at integrated systems testing of critical components comprising the front end of Project X. This program, known as the Project X Injector Experiment (PXIE), is being undertaken as a key component of the larger Project X R&D program. The successful completion of this program will validate the concept for the Project X front end, thereby minimizing a primary technical risk element within Project X. PXIE is currently under construction at Fermilab and will be completed over the period FY12-17. PXIE will include an H- ion source, a CW 2.1-MeV RFQ and two superconductive RF (SRF) cryomodules providing up to 25 MeV energy gain at an average beam current of 1 mA (upgradable to 2 mA). Successful systems testing will also demonstrate the viability of novel front end technologies that are expected find applications beyond Project X.


title: Leveraging a Relationship with Biology to Expand a Relationship with Physics

abstract: This work examines how experiences in one disciplinary domain (biology) can impact the relationship a student builds with another domain (physics). We present a model for disciplinary relationships using the constructs of identity, affect, and epistemology. With these constructs we examine an ethnographic case study of a student who experienced a significant shift in her relationship with physics. We describe how this shift demonstrates (1) a stronger identification with physics, (2) a more mixed affective stance towards physics, and (3) more expert-like ways of knowing in physics. We argue that recruiting the students relationship with biology into experiences of learning physics impacted her relationship with physics as well as her sense of how physics and biology are linked.


title: Type-IV Pilus Deformation Can Explain Retraction Behavior

abstract: Polymeric filament like type IV Pilus (TFP) can transfer forces in excess of 100pN during their retraction before stalling, powering surface translocation(twitching). Single TFP level experiments have shown remarkable nonlinearity in the retraction behavior influenced by the external load as well as levels of PilT molecular motor protein. This includes reversal of motion near stall forces when the concentration of the PilT protein is lowered significantly. In order to explain this behavior, we analyze the coupling of TFP elasticity and interfacial behavior with PilT kinetics. We model retraction as reaction controlled and elongation as transport controlled process. The reaction rates vary with TFP deformation which is modeled as a compound elastic body consisting of multiple helical strands under axial load. Elongation is controlled by monomer transport which suffer entrapment due to excess PilT in the cell periplasm. Our analysis shows excellent agreement with a host of experimental observations and we present a possible biophysical relevance of model parameters through a mechano-chemical stall force map


title: Diffraction, forward physics and soft QCD results from CMS

abstract: We report on the recent CMS measurements of soft hadron-hadron production, including the measurements of inclusive single- and double-diffractive cross sections, as well as the measurement of pseudorapidity distributions, and of the leading charged particle and leading jet cross sections. We present the cross sections for low-pT forward jets production, and production of dijets with large rapidity separation (central-forward, and forward-backward jets). Studies of double parton scattering using four-jet and W+2-jet events, and the cross sections for Drell-Yan with associated jet production are discussed as well. The results, corresponding to the proton-proton center-of-mass energy of 7 or 8 TeV, are compared to predictions of various Monte Carlo models.


title: Computational Gravitational Dynamics with Modern Numerical Accelerators

abstract: We review the recent optimizations of gravitational $N$-body kernels for running them on graphics processing units (GPUs), on single hosts and massive parallel platforms. For each of the two main $N$-body techniques, direct summation and tree-codes, we discuss the optimization strategy, which is different for each algorithm. Because both the accuracy as well as the performance characteristics differ, hybridizing the two algorithms is essential when simulating a large $N$-body system with high-density structures containing few particles, and with low-density structures containing many particles. We demonstrate how this can be realized by splitting the underlying Hamiltonian, and we subsequently demonstrate the efficiency and accuracy of the hybrid code by simulating a group of 11 merging galaxies with massive black holes in the nuclei.


title: A Lattice Path Interpretation of the Diamond Product

abstract: The diamond product is the poset operation that when applied to the face lattices of two polytopes results in the face lattice of the Cartesian product of the polytopes. Application of the diamond product to two Eulerian posets is a bilinear operation on the cd-indices of the two posets, yielding a product on cd-polynomials. A lattice path interpretation is provided for this product of two cd-monomials.


title: Maximizing graph probability under conditionally exponential models

abstract: Designing reliable networks consists in finding topological structures, which are able to successfully carry out desired processes and operations. When this set of activities performed within a network are unknown and the only available information is a probabilistic model reflecting topological network features, highly probable networks are regarded as "reliable", in the sense of being consistent with those probabilistic model. In this paper we are studying the reliability maximization, based on the Exponential Random Graph Model (ERGM), whose statistical properties has been widely used to capture complex topological feature of real-world networks. Under such models the probability of a network is maximized when specified structural properties appear in the network. However, the search of maximally reliable (highly probable) networks might result in difficult combinatorial optimization problems and an important goal of this work is to translate them into solvable systems of linear constraints. Analytical and numerical results are provided, using exact optimization techniques and efficient computer implementation.


title: SN2014J gamma-rays from the 56Ni decay chain

abstract: The measurement of gamma-ray lines from the decay chain of 56Ni provides unique information about the explosion in supernovae. The 56Ni freshly-produced in the supernova powers the optical light curve, as it emits gamma-rays upon its radioactive decay to 56Co and then 56Fe. Gamma-ray lines from 56Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the 56Ni ejecta. With the SPI spectrometer on INTEGRAL and using an improved instrumental background method, we detect the two main lines from 56Co decay at 847 and 1238 keV from SN2014J at 3.3 Mpc, significantly Doppler-broadened, and at intensities (3.65+/-1.21) 10^-4 and (2.27+/-0.69) 10^-4 ph cm^-2s^-1, respectively, at brightness maximum. We measure their rise towards a maximum after about 60-100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from 56Co decay (0.62+/-0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of 56Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of 56Ni and of the unburnt material which occults the 56Co emission. Comparing light curves from different plausible models, the resulting 56Ni mass is determined as 0.49+/-0.09 M_sol.


title: The ample cone of moduli spaces of sheaves on the plane

abstract: Fix a Chern character of a stable sheaf on the plane. Assume either the rank is at most 6 or the rank and first Chern class are coprime and the discriminant is sufficiently large. We use recent results of Bayer and Macri on Bridgeland stability to compute the ample cone of the moduli space of Gieseker semistable sheaves with the given Chern character. We recover earlier results, such as those by Stromme and Yoshioka, as special cases.


title: Ionization Cooling for Muon Experiments

abstract: Possible application for muon experiments such as mu2e is discussed of the initial part of the ionization cooling channel originally developed for muon collider. It is shown that with the FNAL Booster as the proton driver the mu2e sensitivity can be increased by two orders of magnitude compared to the presently considered experiment.


title: Pre-main-sequence isochrones -- III. The Cluster Collaboration isochrone server

abstract: We present an isochrone server for semi-empirical pre-main-sequence model isochrones in the following systems: Johnson-Cousins, Sloan Digital Sky Survey, Two-Micron All-Sky Survey, Isaac Newton Telescope (INT) Wide-Field Camera, and INT Photometric H$\alpha$ Survey (IPHAS)/UV-Excess Survey (UVEX). The server can be accessed via the Cluster Collaboration webpage {http://www.astro.ex.ac.uk/people/timn/isochrones/}. To achieve this we have used the observed colours of member stars in young clusters with well-established age, distance and reddening to create fiducial loci in the colour-magnitude diagram. These empirical sequences have been used to quantify the discrepancy between the models and data arising from uncertainties in both the interior and atmospheric models, resulting in tables of semi-empirical bolometric corrections (BCs) in the various photometric systems. The model isochrones made available through the server are based on existing stellar interior models coupled with our newly derived semi-empirical BCs. As part of this analysis we also present new cluster parameters for both the Pleiades and Praesepe, yielding ages of $135^{+20}_{-11}$ and $665^{+14}_{-7}\,\rm{Myr}$ as well as distances of $132 \pm 2$ and $184 \pm 2\,\rm{pc}$ respectively (statistical uncertainty only).


title: Iterated socles and integral dependence in regular rings

abstract: Let $R$ be a formal power series ring over a field, with maximal ideal $\mathfrak m$, and let $I$ be an ideal of $R$ such that $R/I$ is Artinian. We study the iterated socles of $I$, that is the ideals which are defined as the largest ideal $J$ with $J\mathfrak m^s\subset I$ for a fixed positive integer $s$. We are interested in these ideals in connection with the notion of integral dependence of ideals. In this article we show that the iterated socles are integral over $I$, with reduction number one, provided $s \leq \text{o}(I_1(\varphi_d))-1$, where $\text{o}(I_1(\varphi_d))$ is the order of the ideal of entries of the last map in a minimal free $R$-resolution of $R/I$. In characteristic zero, we also provide formulas for the generators of iterated socles whenever $s\leq \text{o}(I_1(\varphi_d))$. This result generalizes previous work of Herzog, who gave formulas for the socle generators of any ${\mathfrak m}$-primary homogeneous ideal $I$ in terms of Jacobian determinants of the entries of the matrices in a minimal homogeneous free $R$-resolution of $R/I$. Applications are given to iterated socles of determinantal ideals with generic height. In particular, we give surprisingly simple formulas for iterated socles of height two ideals in a power series ring in two variables. These generators are suitable determinants obtained from the Hilbert-Burch matrix.


title: Strong Far-IR Cooling Lines, Peculiar CO Kinematics and Possible Star Formation Suppression in Hickson Compact Group 57

abstract: We present [C II] and [O I] observations from Herschel and CO(1-0) maps from the Combined Array for{\dag} Research in Millimeter Astronomy (CARMA) of the Hickson Compact Group HCG 57, focusing on the galaxies HCG 57a and HCG 57d. HCG 57a has been previously shown to contain enhanced quantities of warm molecular hydrogen consistent with shock and/or turbulent heating. Our observations show that HCG 57d has strong [C II] emission compared to L$_{\rm FIR}$ and weak CO(1-0), while in HCG 57a, both the [C II] and CO(1-0) are strong. HCG 57a lies at the upper end of the normal distribution of [C II]/CO and [C II]/FIR ratios, and its far-IR cooling supports a low density warm diffuse gas that falls close to the boundary of acceptable PDR models. However, the power radiated in the [C II] and warm H$_2$ emission have similar magnitudes, as seen in other shock-dominated systems and predicted by recent models. We suggest that shock-heating of the [C II] is a viable alternative to photoelectric heating in violently disturbed diffuse gas. The existence of shocks is also consistent with peculiar CO kinematics in the galaxy, indicating highly non-circular motions are present. These kinematically disturbed CO regions also show evidence of suppressed star formation, falling a factor of 10-30 below normal galaxies on the Kennicutt-Schmidt relation. We suggest that the peculiar properties of both galaxies are consistent with a highly dissipative off-center collisional encounter between HCG 57d and 57a, creating ring-like morphologies in both systems. Highly dissipative gas-on-gas collisions may be more common in dense groups because of the likelihood of repeated multiple encounters. The possibility of shock-induced SF suppression may explain why a subset of these HCG galaxies have been found previously to fall in the mid-infrared green valley.


title: Computing Eigen-Emittances from Tracking Data

abstract: In a strongly nonlinear system the particle distribution in the phase space may develop long tails which contribution to the covariance (sigma) matrix should be suppressed for a correct estimate of the beam emittance. A method is offered based on Gaussian approximation of the original particle distribution in the phase space (Klimontovich distribution) which leads to an equation for the sigma matrix which provides efficient suppression of the tails and cannot be obtained by introducing weights. This equation is easily solved by iterations in the multi-dimensional case. It is also shown how the eigen-emittances and coupled optics functions can be retrieved from the sigma matrix in a strongly coupled system. Finally, the developed algorithm is applied to 6D ionization cooling of muons in HFOFO channel.


title: Ground-state Bethe root densities and quantum phase transitions

abstract: Exactly solvable models provide a unique method, via qualitative changes in the distribution of the ground-state roots of the Bethe Ansatz equations, to identify quantum phase transitions. Here we expand on this approach, in a quantitative manner, for two models of Bose--Einstein condensates. The first model deals with the interconversion of bosonic atoms and molecules. The second is the two-site Bose--Hubbard model, widely used to describe tunneling phenomena in Bose--Einstein condensates. For these systems we calculate the ground-state root density. This facilitates the determination of analytic forms for the ground-state energy, and associated correlation functions through the Hellmann--Feynman theorem. These calculations provide a clear identification of the quantum phase transition in each model. For the first model we obtain an expression for the molecular fraction expectation value. For the two-site Bose--Hubbard model we find that there is a simple characterisation of condensate fragmentation.

categories: nlin.SI cond-mat.quant-gas math-ph math.MP


title: Laser-modified Angular Distribution of Muon Decay

abstract: We show theoretically that the angular distribution of decay rate of muon can be changed dramatically by embedding the decaying muon in a strong linearly polarized laser field. Evaluating the S-matrix elements taking all electronic multiphoton processes into account. The results suggest the muon may have internal structures instead structureless as in the standard model.


title: A Learning Based Approach to Control Synthesis of Markov Decision Processes for Linear Temporal Logic Specifications

abstract: We propose to synthesize a control policy for a Markov decision process (MDP) such that the resulting traces of the MDP satisfy a linear temporal logic (LTL) property. We construct a product MDP that incorporates a deterministic Rabin automaton generated from the desired LTL property. The reward function of the product MDP is defined from the acceptance condition of the Rabin automaton. This construction allows us to apply techniques from learning theory to the problem of synthesis for LTL specifications even when the transition probabilities are not known a priori. We prove that our method is guaranteed to find a controller that satisfies the LTL property with probability one if such a policy exists, and we suggest empirically with a case study in traffic control that our method produces reasonable control strategies even when the LTL property cannot be satisfied with probability one.


title: Efficient construction of homological Seifert surfaces

abstract: Let $\Omega$ be a bounded domain of $\mathbb{R}^3$ whose closure $\overline{\Omega}$ is polyhedral, and let $\mathcal{T}$ be a triangulation of $\overline{\Omega}$. Assuming that the boundary of $\Omega$ is sufficiently regular, we provide an explicit formula for the computation of homological Seifert surfaces of any $1$-boundary $\gamma$ of $\mathcal{T}$; namely, $2$-chains of $\mathcal{T}$ whose boundary is $\gamma$. It is based on the existence of special spanning trees of the complete dual graph of $\mathcal{T}$, and on the computation of certain linking numbers associated with those spanning trees. If the triangulation $\mathcal{T}$ is fine, the explicit formula is too expensive to be used directly. For this reason, making also use of a simple elimination procedure, we devise a fast algorithm for the computation of homological Seifert surfaces. Some numerical experiments illustrate the efficiency of this algorithm.


title: Connecting the Hadron Mass Scale to the Fundamental Mass Scale of Quantum Chromodynamics

abstract: Establishing an explicit connection between the long distance physics of confinement and the dynamical interactions of quarks and gluons at short distances has been a long-sought goal of quantum chromodynamics. Using holographic QCD, we derive a direct analytic relation between the scale $\kappa$ which determines the masses of hadrons and the scale $\Lambda_{s}$ which controls the predictions of perturbative QCD at very short distances. The resulting prediction $\Lambda_{s}=0.341\pm0.032$ GeV in the $\overline{MS}$ scheme agrees well with the experimental average $0.339\pm0.016$ GeV. We also derive a relation between $\Lambda_{s}$ and the QCD string tension $\sigma$. This connection between the fundamental hadronic scale underlying the physics of quark confinement and the perturbative QCD scale controlling hard collisions can be carried out in any renormalization scheme.


title: Sound Focusing by Acoustic Luneburg Lens

abstract: Luneburg lens is a gradient index lens that focuses the incoming wave on the opposite side of the lens without aberration. We developed a two-dimensional acoustic Luneburg lens by changing the refractive index of the medium in the lens. It has a cylindrical shape and is composed of 701 aluminum columns with various radii of less than 1cm. It focuses the incoming sound wave on the edge of the opposite side of the lens without aberration as well in the frequency range of 1,000 - 3,200Hz. It increases the amplitude of the incoming sound wave by 3 - 4 times, and the amplification corresponds to the sound level difference of 10 - 15dB. The ability of the acoustic Luneburg lens as sonar was checked by VU meters in the air. Acoustic Luneburg lens has a simple structure and easy to build. It could be a strong candidate of a next generation of sonar.


title: Tidal Waves in $^{102}$Pd: A Phenomenological Analysis

abstract: Rotational and electromagnetic properties of the yrast band in $^{102}$Pd are analyzed in terms of a phenomenological phonon model that includes anharmonic terms. Both the moment of inertia and $B(E2)$'s are well reproduced by the model, providing an independent confirmation of the multi-phonon picture recently proposed. The (empirical) dependence of the phonon-phonon interaction on the phonon frequency, in Ru, Pd, and Ru isotopes, follows the expectations from particle-vibration coupling.


title: Images encryption using AES and variable permutations

abstract: This work proposes a different procedure to encrypt images of 256 grey levels and colour, using the symmetric system Advanced Encryption Standard with a variable permutation in the first round, after the x-or operation. Variable permutation means using a different one for each input block of 128 bits. In this vein, an algorithm is constructed that defines a Bijective function between sets Nm = {n in N, 0 <= n < fac(m)} with n >= 2 and Pm = {pi, pi is a permutation of 0, 1, ..., m-1}. This algorithm calculates permutations on 128 positions with 127 known constants. The transcendental numbers are used to select the 127 constants in a pseudo-random way. The proposed encryption quality is evaluated by the following criteria: Correlation; horizontal, vertical and diagonal, Entropy and Discrete Fourier Transform. The latter uses the NIST standard 800-22. Also, a sensitivity analysis was performed in encrypted figures. Furthermore, an additional test is proposed which considers the distribution of 256 shades of the three colours; red, green and blue for colour images. On the other hand, it is important to mention that the images are encrypted without loss of information because many banking companies and some safety area countries do not allow the figures to go through a compression process with information loss. i.e., it is forbidden to use formats such as JPEG.


title: Relaxation Principle in (n+1) Dimensions

abstract: Why do we live in a (3+1) dimensional universe? In this paper, I review the "relaxation principle" first introduced by Karch and Randall (arXiv:hep-th/0506053v2) and generalize its ideas to an arbitrary (n+1) dimensions. This is done by referring to the Friedman equation and the scaling solution to derive the energy densities of the non-interacting and interacting branes, and then formulating their dimensionalities. I also demonstrate that the largest interacting d-brane in n spatial dimensions is always the (n-2)-brane, and that such dimensionality constraint "relaxed" our universe from nine dimensions to three dimensions.


title: K I 7699 {\AA} and related shell lines during the recent eclipse of $\epsilon$ Aurigae

abstract: We report high-resolution (R = 30\ 000, 45\ 000 and 75\ 000) echelle and medium-resolution (R = 22\ 000 and 10\ 000) spectroscopic observations of the long-period, eclipsing binary $\epsilon$ Aurigae during the 2009 - 2011 eclipse. Low-excitation shell lines, viz, the K {\sevensize \rm I} line at 7699 \AA\ (with 346 data points), Cr {\sevensize \rm I} lines at 5345.807 \AA\ and 5348.326 \AA\ and Fe {\sevensize \rm I} line at 5110.435 \AA\ which originated from the disk shaped secondary, H$\alpha$ and the shell components of the Na D$_{1}$ and D$_{2}$ lines show significant variation in their shapes and radial velocities during the eclipse. The equivalent width curve shown by the K {\sevensize \rm I} line around the ingress and egress phases indicates that the gas density in the trailing edge is about a factor of two higher than the density in the leading edge. Using a geometrical model, in which a homogeneous, cylindrical Keplarian disk eclipses the F0Ia primary star and the shell absorption lines originate from the gaseous atmosphere around an opaque disk, we fit the equivalent width and the radial velocity curves of the K {\sevensize \rm I} line covering the full eclipse. A reasonably good fit can be achieved by a low-mass binary model where the mass of the central star of the disk is 5.4 M$_{\odot}$ and the mass of the primary is 2.5 M$_{\odot}$ and a disk size of 8.9 AU. The low-mass of the primary, with enhanced $s$-process elements found by Sadakane et al. (2010), supports that it is a post-AGB F supergiant. For the high-mass binary model, the modelled radial velocity curve deviates significantly from the observations.


title: Calorimetry of a harmonically trapped Bose gas

abstract: We experimentally study the energy-temperature relationship of a harmonically trapped Bose-Einstein condensate by transferring a known quantity of energy to the condensate and measuring the resulting temperature change. We consider two methods of heat transfer, the first using a free expansion under gravity and the second using an optical standing wave to diffract the atoms in the potential. We investigate the effect of interactions on the thermodynamics and compare our results to various finite temperature theories.


title: Efficient Feature Group Sequencing for Anytime Linear Prediction

abstract: We consider \textit{anytime} linear prediction in the common machine learning setting, where features are in groups that have costs. We achieve anytime (or interruptible) predictions by sequencing the computation of feature groups and reporting results using the computed features at interruption. We extend Orthogonal Matching Pursuit (OMP) and Forward Regression (FR) to learn the sequencing greedily under this group setting with costs. We theoretically guarantee that our algorithms achieve near-optimal linear predictions at each budget when a feature group is chosen. With a novel analysis of OMP, we improve its theoretical bound to the same strength as that of FR. In addition, we develop a novel algorithm that consumes cost $4B$ to approximate the optimal performance of \textit{any} cost $B$, and prove that with cost less than $4B$, such an approximation is impossible. To our knowledge, these are the first anytime bounds at \textit{all} budgets. We test our algorithms on two real-world data-sets and evaluate them in terms of anytime linear prediction performance against cost-weighted Group Lasso and alternative greedy algorithms.